One and two-parameter bifurcations to divergence and flutter in the three-dimensional motions of a fluid conveying viscoelastic tube with D4-symmetry

The loss of stability of the trivial downhanging equilibrium position of a slender circular tube conveying incompressible fluid flow is studied. The tube is clamped at its upper end and free at its lower end. Inbetween the three-dimensional transversal motion is constrained by an elastic support which is considered to beD ₄-symmetric, that is, has the symmetry of the square (Figure 1). Kirchhoff's rod theory and the Kelvin-Voigt viscoelastic law are used to derive the tube equations under the assumption of large displacement but small strain.The stability analysis is performed making use of the methods of equivariant bifurcation theory, that is, making use of the symmetry properties of the original system in deriving the amplitude equations of the critical modes. All cases of loss of stability which are possible for generic one-parameter bifurcations and the coincident case of a zero root and a purely imaginary pair of roots are investigated.Dedicated to Professor P. R. Sethna on the Occasion of His 70th Birthday     

Hopf Bifurcation in the presence of symmetry

Using group theoretic techniques, we obtain a generalization of the Hopf Bifurcation Theorem to differential equations with symmetry, analogous to a static bifurcation theorem of Cicogna. We discuss the stability of the bifurcating branches, and show how group theory can often simplify stability calculations. The general theory is illustrated by three detailed examples: O(2) acting on R ², O(n) on R ⁿ , and O(3) in any irreducible representation on spherical harmonics.The work of second author was also supported by a visiting position in the Department of Mathematics, University of Houston     

Effects of nanoparticles mean diameter on the particle migration and thermo-hydraulic behavior of laminar mixed convection of a nanofluid in an inclined tube

Laminar mixed convection of a nanofluid consisting of Al₂O₃ and water through an inclined tube has been investigated numerically. As mathematical model two-phase mixture model has been adopted, thus three dimensional elliptical governing equations have been solved to understand the flow behavior at different Re–Gr combinations. Control volume technique is used for discretization of the governing equations. For the convective and diffusive terms the second order upwind method was used while the SIMPLEC procedure was adopted for the velocity–pressure coupling. For different nanoparticle mean diameters and tube inclinations thermo-fluid parameters such as secondary flow, axial velocity profiles, nanoparticles distribution at the tube cross section, axial evolution of peripheral average convective heat transfer coefficient and pressure drop along the tube, have been presented and discussed. Maximum enhancement on the heat transfer coefficient is seen at tube inclination of 45°.     

Mechanical quadrature methods and extrapolation for solving nonlinear boundary Helmholtz integral equations

This paper presents mechanical quadrature methods (MQMs) for solving nonlinear boundary Helmholtz integral equations. The methods have high accuracy of order O(h ³) and low computation complexity. Moreover, the mechanical quadrature methods are simple without computing any singular integration. A nonlinear system is constructed by discretizing the nonlinear boundary integral equations. The stability and convergence of the system are proved based on an asymptotical compact theory and the Stepleman theorem. Using the h ³-Richardson extrapolation algorithms (EAs), the accuracy to the order of O(h ⁵) is improved. To slove the nonlinear system, the Newton iteration is discussed extensively by using the Ostrowski fixed point theorem. The efficiency of the algorithms is illustrated by numerical examples.     

Some problems in the Z-C-X space

The new concepts of the Z-C-X space and excellent cone are introduced. Some problems of random semiclosed 1-set-contractive operator are investigated in the Z-C-X space. At first, an important inequality is proved. Secondly, several new conclusions are proved by means of random fixed point index in the theory of random topological degree. A random solution of a class of random operator equations under conditions of imitating the parallelogram law is obtained, famous Altman’s theorem is generalized in partially ordered Z-C-X space. Therefore, some new results are obtained.     

THEORY OF GRAVITATIONAL RADIATION OF THE(Ω,Aμr)-FIELD

Further exploration of the Ω-field theory as first proposed by Yu (1989) is here presented to cover the equation of motion of a test particle which induces gravitational radiation. The same theory is shown to contain an exact gravitational radiation equation derived as a logical consequence of field equations without extra postulates. In this general dynamic context the theory is renamed "The (Ω,A_μr) field Theory".     

流体饱和多孔介质动力问题的显式时域解法

为了简化分析,Zienkiewicz等基于Biot理论,在忽略流体相对于土骨架运动的加速度条件下,建立了以土骨架位移u和孔隙流体压力p为基本变量的u-p格式饱和两相介质动力方程。针对该u-p方程,在空间上,采用伽辽金法有限元离散,并结合对角化形式的质量矩阵和流体压缩矩阵,忽略相邻结点间的惯性和流体压缩量间的耦合作用。在时域内,基于杜修力等提出的显式算法和Euler预估-校正法,建立了一种具有二阶精度的全显式时域积分法。采用一维饱和土模型,对比提出算法的数值解与Simon方法的解析解,发现两者吻合良好,验证了本文方法的正确性。并分析了饱和土二维动力问题,以及渗透系数和排水条件对饱和土动力响应的影响。     

Stability for manifolds of equilibrium states of fractional generalized Hamiltonian systems

For a fractional generalized Hamiltonian system, in terms of Riesz derivatives, stability theory for the manifolds of equilibrium states is presented. The gradient representation and second order gradient representation of a fractional generalized Hamiltonian system are studied, and the conditions under which the system can be considered as a gradient system and a second order gradient system are given, respectively. Then, equilibrium equations, disturbance equations, and first approximate equations of a fractional generalized Hamiltonian system are obtained. A theorem for the stability of the manifolds of equilibrium states of the general autonomous system is used to a fractional generalized Hamiltonian system, and three propositions on the stability of the manifolds of equilibrium states of the system are investigated. As the special cases of this article, the conditions which a fractional generalized Hamiltonian system can be reduced to a generalized Hamiltonian system, a fractional Hamiltonian system and a Hamiltonian system are given, respectively, and the stability theory for the manifolds of equilibrium states of these systems are obtained. Further, a fractional dynamical system and a fractional Volterra model of the three species groups are given to illustrate the method and results of the application. Finally, by using the method in this paper, we construct a new kind of fractional dynamical model, i.e. the fractional Hénon–Heiles model, and we study its stability of the manifolds of equilibrium states.     

The numerical solution of the flow in a general bifurcating channel at moderately high Reynolds number using boundary-fitted co-ordinates,primitive variables and Newton iteration

A numerical code has been implemented for the numerical solution of the steady, incompressible Navier–Stokes equations using primitive variables in a bifurcating channel. A boundary-fitted, numerically generated grid is placed onto the domain of the channel which is transformed into either a rectilinear C- or T-shaped region. The differenced equations are solved using Newton's iteration which makes upwinding at high Reynolds number unnecessary. Practical implications of inverting the huge Jacobian matrix of Newton's method are discussed. The results have relative error of 2–3 × 10^?3 at Reynolds number 100, with T-geometry being marginally but significantly more accurate than C-geometry. Results have been obtained for Reynolds numbers up to 1000 for three bifurcations one of which models the carotid arterial bifurcation in the human head. For this latter bifurcation the wall shear stress is calculated in connection with the onset of atherosclerosis. Finally, the results of flows having different daughter tube end pressures are presented.     

RENEWAL OF BASIC LAWS AND PRINCIPLES FOR POLAR CON-TINUUM THEORIES (Ⅱ)—MICROMORPHIC CONTINUUM THEORY AND COUPLE STRESS THEORY

The purpose is to reestablish the balance laws of momentum, angular momentum and energy and to derive the corresponding local and nonlocal balance equations for micromorphic continuum mechanics and couple stress theory. The desired results for micromorphic continuum mechanics and couple stress theory are naturally obtained via direct transitions and reductions from the coupled conservation law of energy for micropolar continuum theory, respectively. The basic balance laws and equations for micromorphic continuum mechanics and couple stress theory are constituted by combining these results derived here and the traditional conservation laws and equations of mass and microinertia and the entropy inequality. The incomplete degrees of the former related continuum theories are clarified. Finally, some special cases are conveniently derived. Foundation items: the National Natural Science Foundation of China (10072024); the Research Foundation of Liaoning Education Committee (990111001) Biography: DAI Tian-min (1931≈)     

Nonlinear statics and dynamics of a simply supported nonuniform tube conveying an incompressible inviscid fluid

Nonlinear static and dynamic behaviour of a simply supported fluid-conveying tube, which has a constant inner diameter and a variable thickness is analysed analytically and numerically. Nonlinear static bending is considered in two loading cases: (i) a tube subjected to supercritical axial compressive forces acting at its edges or (ii) a tube loaded by concentrated bending moments, which provide a symmetrical (with respect to the mid-span) shape of a tube. The nonlinear governing equations of motions are derived by using Hamilton's principle. The elementary plug flow theory of an incompressible inviscid fluid is adopted for modelling a fluid–structure interaction. The flow velocity is taken as the sum of a principal constant ‘mean’ velocity component and a fairly small pulsating component. Firstly, eigenfrequencies and eigenmodes of a deformed tube are found from linearised equations of motions. Then resonant nonlinear oscillations of a tube about its deformed static equilibrium position in a plane of static bending are considered. A multiple scales method is used and a weak resonant excitation by the flow pulsation is considered in a single-mode regime and in a bi-modal regime (in the case of an internal parametric resonance) and the stability of each of them is examined. The brief parametric study of these regimes of motions is carried out.     

Mechanical Quadrature Methods and Extrapolation Algorithms for Boundary Integral Equations with Linear Boundary Conditions in Elasticity

By potential theory, elastic problems with linear boundary conditions are converted into boundary integral equations (BIEs) with logarithmic and Cauchy singularity. In this paper, a mechanical quadrature method (MQMs) is presented to deal with the logarithmic and the Cauchy singularity simultaneously for solving the boundary integral equations. The convergence and stability are proved based on Anselone??s collective compact and asymptotical compact theory. Furthermore, an asymptotic expansion with odd powers of errors is presented, which possesses high accuracy order O(h ³). Using h ³?Richardson extrapolation algorithms (EAs), the accuracy order of the approximation can be greatly improved to O(h ⁵), and an a posteriori error estimate can be obtained for constructing a self-adaptive algorithm. The efficiency of the algorithm is illustrated by examples.     

The axisymmetric mixed boundary-value problem of the vertical vibration of a rigid foundation on saturated layered soil subgrade

Based on the theory of elastic wave propagation in saturated soil subgrade established by the author of this paper, the axisymmetric vertical vibration of a rigid circular foundation resting on partially saturated soil subgrade which is composed of a dry elastic layer and a saturated substratum is studied. The analysis relied on the use of integral transform techniques and a pair of dual integral equations governing the vertical vibration of the rigid foundation is listed under the consideration of mixed boundary-value condition. The results are reduced to the case for saturated half-space. The set of dual integral equations are reduced to a Fredholm integral equation of the second kind and solved by numerical procedures. Numerical examples are given at the end of the paper and plots of the dynamic compliance coefficient C_v versus the dimensionless frequency a₀ are presented. Foundation item: the National Natural Science Foundation of China (59579018) Biography: CHEN Sheng-li (1974-)     

Optimization study of spray detonation initiation by electric discharges

Development of air-breathing pulse detonation engines is faced with a challenging problem of detonation initiation in fuel sprays at distances feasible for propulsion applications. Extensive experimental study on initiation of a confined n-hexane spray detonation in air by electric discharges is reported. It is found that for direct initiation of spray detonation with minimal energy requirements (1) it is worth to use one discharger located near the closed end of a detonation tube and at least one additional discharger downstream from it to be triggered in-phase with primary shock wave arrival; (2) the discharge area should be properly insulated to avoid electric loss to metal tube walls; (3) discharge duration should be minimized to at least 50 μs; (4) discharge channel should preferably occupy a large portion of a tube cross-section; (5) test tube should be preferably of a diameter close to the limiting tube diameter; (6) gradual transition between the volume with electric discharger and the tube should be used; and (7) a powerful electric discharger utilized for generating a primary shock wave can be replaced by a primary shock wave generator comprising a relatively low-energy electric discharger, Shchelkin spiral, and tube coil. With all these principles implemented, the rated electric energy of about 100 J was required to initiate n-hexane spray–air detonation in a 28-mm tube at a distance of about 1 m from the atomizer. PACS 02.60.Cb; 05.10.Ln; 47.11.+j; 47.15.Cb; 47.40.Nm This paper was based on work that was presented at the 19th Inter-national Colloquium on the Dynamics of Explosions and Reactive Systems, Hakone, Japan, July 27 - August 1, 2003     

The dynamic response of an incompressible non-linearly elastic membrane tube subjected to a dynamic extension

The dynamic response of an isotropic hyperelastic membrane tube, subjected to a dynamic extension at its one end, is studied. In the first part of the paper, an asymptotic expansion technique is used to derive a non-linear membrane theory for finite axially symmetric dynamic deformations of incompressible non-linearly elastic circular cylindrical tubes by starting from the three-dimensional elasticity theory. The equations governing dynamic axially symmetric deformations of the membrane tube are obtained for an arbitrary form of the strain-energy function. In the second part of the paper, finite amplitude wave propagation in an incompressible hyperelastic membrane tube is considered when one end is fixed and the other is subjected to a suddenly applied dynamic extension. A Godunov-type finite volume method is used to solve numerically the corresponding problem. Numerical results are given for the Mooney-Rivlin incompressible material. The question how the present numerical results are related to those obtained in the literature is discussed.     

On the Ill-Posedness of the Prandtl Equations in Three-Dimensional Space

In this paper, we give an instability criterion for the Prandtl equations in three-dimensional space, which shows that the monotonicity condition on tangential velocity fields is not sufficient for the well-posedness of the three-dimensional Prandtl equations, in contrast to the classical well-posedness theory of the two-dimensional Prandtl equations under the Oleinik monotonicity assumption. Both linear stability and nonlinear stability are considered. This criterion shows that the monotonic shear flow is linearly stable for the three-dimensional Prandtl equations if and only if the tangential velocity field direction is invariant with respect to the normal variable, and this result is an exact complement to our recent work (A well-posedness theory for the Prandtl equations in three space variables. arXiv:1405.5308, 2014) on the well-posedness theory for the three-dimensional Prandtl equations with a special structure.     

Experimental study of ethylene glycol-based Al2O3 nanofluid turbulent heat transfer enhancement in the corrugated tube with twisted tapes

In this study, fluid flow of the Al₂O₃/ethylene glycol (EG) nanofluid in a corrugated tube fitted with twisted tapes were experimentally studied under turbulent flow conditions. The experiments with different twists ratio and different nanofluid concentration were performed under similar operation condition. The investigated ranges are (1) three different Al₂O₃ concentrations: 0.5, 1 and 1.5 % by volume (2) three different twist ratios of twisted tape: y/w = 2, 3.6 and 5 and (3) Reynolds number from 6000 to 30,000. Regarding the experimental data, utilization of twists together with nanofluids tends to increase heat transfer and friction factor as compared with the base fluid. In addition, heat transfer performances were weakened by using for high nanoparticle concentration. The thermal performances of the heat exchanger with nanofluid and twisted tapes were evaluated for the assessment of overall improvement in thermal behavior. Over the range studied, the maximum thermal performance factor 4.2 is found with the use of Al₂O₃/EG nanofluid at concentration of 0.5 % by volume in corrugated tube together with twisted tape at twist ratio of 2.     

The viscous air flow pattern in the Stefan diffusion tube

In this paper, we consider the problem of the binary viscous diffusion of vapour through a Stefan tube, which is the model of an elementary capillary. While some preceding results in particular cases supposed parabolic velocity profiles and showed air recirculation, we treat here the general problem of a tube of finite length, submitted to a double viscous diffusion of vapour and air from a liquid surface. The movement of gas is expressed with conservation equations and ideal gas equations. The following added restrictions: constant temperature, no buyoancy effect, no inertial forces, are compatible with a capillary. A numerical solution based on the control volume method is obtained at every point in the tube. The results give the vapour and air flux, describe the circulation pattern and show that the vapour profile of concentration is level. In the lower part of the cylindrical tube space, over a distance of the length of a radius an important radial movement occurs, due to the recirculation of air which changes direction once it reaches the liquid surface. The velocity profile of the gas flow then becomes parabolic in the upper part of the tube.In order to easily obtain a numerical solution, the system of dimensionless equations is expanded to a series and transformed into a set of sub-systems. The little parameter used for this expansion is tied to the vapour concentration on the liquid surface. The solution of the sub-system of order zero, which is easier to compute, represents a good approximation of the complete solution. These solutions are situated in comparison with the Stefan diffusion and show that the influence of the viscous effect on the vapour flux is limited to a few percents.In order to apply the results to porous media where the pores are not so regular, we consider at last the diffusion in a tube including a contracted section in the middle of the tube. Since the diffusion paths are longer, the vapour flux is reduced, while the viscosity effect becomes more considerable. The reduction of the air flux is more significant than that of the vapour. This part of the study provides a better understanding of the diffusion through the pits at the wall fiber, and gives data for the air flux which permeates into the oak wood and produces tannin oxidation and thus discolouration.Nomenclature m _v/ vapour concentration - m _l vapour concentration in equilibrium with its liquid phase - D coefficient of molecular diffusion of a vapour in air (m²/s) - J _a vector density of mass flux of dry air (kg/m²s) - J _v vector density of mass flux of vapour (kg/m²s) - L Capillary length (m) - M _a dry air molar mass (kg/mole) - M _v Vapour molar mass (kg/mole) - P _atm atmospheric pressure - P gas mixture total pressure (Pascal) - R Ideal gas constant (J/mole K) - r _a= R/M _a,r _v=R/M _v - r ₀ tube or capillary radius (m) - T Temperature (K) - u axial component ofV - V gas mixture velocity vector (m/s) - v radial component ofV Greek Letters density of gas (kg/m³) - gas mixture dynamic viscosity (kg/ms) Numerical Values of Parameters D 3×10^–5m²/s (water vapour in air) - 2×10^–5kg/ms - M _a 29×10^–3kg/mole - M _v 18×10^–3kg/mole - T 323K - R 8.31 J/moleK - P _atm 10⁵Pa     

On the adjacent equilibrium method in the stability theory of conservative elastic continua

The relationship of the adjacent equilibrium method, the regular perturbation method and the energy method for neutral equilibrium is studied. It is shown that unlike the adjacent equilibrium method, the regular perturbation method yields, for the problems under consideration, non-homogeneous perturbation equations and that adjacent states of equilibrium do not exist at the bifurcation point. These results are then compared with the result of the energy criterion for neutral equilibrium V₂[u] = 0. It is found that although the physical arguments are different in the three methods, the resulting stability equations are identical; thus showing why the adjacent equilibrium argument, even for cases when it is incorrect, yields correct critical loads. This is followed by a discussion of an incorrect derivation of a stability condition and a notion about a load type introduced in the stability literature, which are consequences of the assumption of the general existence of adjacent equilibrium states at bifurcation points.     

Intermittency as a codimension-three phenomenon

We analyze the interaction of three Hopf modes and show that locally a bifurcation gives rise to intermittency between three periodic solutions. This phenomenon can occur naturally in three-parameter families. Consider a vector fieldf with an equilibrium and suppose that the linearization off about this equilibrium has three rationally independent complex conjugate pairs of eigenvalues on the imaginary axis. As the parameters are varied, generically three branches of periodic solutions bifurcate from the steady-state solution. Using Birkhoff normal form, we can approximatef close to the bifurcation point by a vector field commuting with the symmetry group of the three-torus. The resulting system decouples into phase amplitude equations. The main part of the analysis concentrates on the amplitude equations in R³ that commute with an action ofZ ₂+Z ₂+Z ₂. Under certain conditions, there exists an asymptotically stable heteroclinic cycle. A similar example of such a phenomenon can be found in recent work by Guckenheimer and Holmes. The heteroclinic cycle connects three fixed points in the amplitude equations that correspond to three periodic orbits of the vector field in Birkhoff normal form. We can considerf as being an arbitrarily small perturbation of such a vector field. For this perturbation, the heteroclinic cycle disappears, but an invariant region where it was is still stable. Thus, we show that nearby solutions will still cycle around among the three periodic orbits.