Short-Time Structural Stability of Compressible Vortex Sheets with Surface Tension

A thermodynamically consistent two-phase Stefan problem with temperature-dependent surface tension and with or without kinetic undercooling is studied. It is shown that these problems generate local semiflows in well-defined state manifolds. If a solution does not exhibit singularities, it is proved that it exists globally in time and converges towards an equilibrium of the problem. In addition, stability and instability of equilibria is studied. In particular, it is shown that multiple spheres of the same radius are unstable if surface heat capacity is small; however, if kinetic undercooling is absent, they are stable if surface heat capacity is sufficiently large.     

How to build master equations for complex systems

Typical complex systems, e. g., complex chemical reactions, reaction-diffusion systems, and turbulent fluids are described on a macroscopic level, that is, neglecting fluctuations, with the help of deterministic equations for corresponding variables. In this article it is shown on a phenomenological level, that these systems can be described in terms of integer- or real-valued Markov processes as well, which are governed by master equations. The latter are constructed such that the macroscopic law and the fluctuations around it are reproduced correctly. Stochastic processes defined through master equations can easily be simulated. The efficiency, the stability and the parallelization of the algorithms for stochastic simulations are discussed for some examples. In the last part of the paper it is shown that the same phenomenological approach can be successfully applied to open quantum systems. The wave function is assumed to be a complex valued stochastic process in Hilbert space and the quantum master equation for the statistical operator is regarded as the equation of motion for the two-point correlation function.     

Constitutive equations for micropolar hyper-elastic materials

In this paper, the concept of hyper-elasticity in the micropolar continuum theory is investigated. The restrictions on the fourth-order elasticity tensors are investigated. Using the representation theorems, a general form of constitutive equations for micropolar hyper-elastic isotropic materials is presented. As some special cases, generalizations of the neo-Hookean and Mooney-Rivlin type materials to the micropolar continuum theory are presented. The generalized constitutive equations reduce to those of the microplar linear elasticity theory when the deformations are infinitesimal. Also, Updated Lagrangian finite element formulations for the micropolar hyper-elastic materials are presented. Considering two planar examples, it is shown that an increase in the micropolar parameter results in the reduction of the deformation of the bodies. Also, it is shown that for a specimen with very small dimensions, e.g. in the micron level, the micropolar effects are more sensible. Furthermore, it is shown that the influence of the micropolar parameters is dependent not only on the size of the body, but also to its geometry and loading conditions. For the problems in which the deformation is very close to a homogeneous state, the micropolar effects are negligible.     

Necessary and sufficient conditions for complex balancing in chemical kinetics

Summary In a recent publication (Horn & Jackson [1]) it was shown that complex balancing together with mass action type rate laws ensures certain stability properties of a kinetic system, thereby precluding sustained oscillations, bistability and other types of irregular dynamics. In this paper a necessary condition for complex balancing in general kinetics and necessary and sufficient conditions for complex balancing in mass action systems are derived. A theorem is stated which excludes the occurence of equilibria in certain composition regions of general kinetic systems. For mass action systems it is shown that it is sometimes true that the algebraic structure of the reactions suffices to ensure complex balancing, while in other cases complex balancing occurs only if certain relations between the rate constants are satisfied. The number of these relations, called the deficiency of the mass action system is determined by the algebraic structure of the set of reactions underlying that system.     

Magnetoelastic interactions in a soft ferromagnetic body with a nonlinear law of magnetization: Some applications

Linearized equations and boundary conditions of a magnetoelastic ferromagnetic body are obtained with the nonlinear law of magnetization. Magnetoelastic interactions in a multi-domain ferromagnetic materials are considered for magneto soft materials, i.e. the case when the magnetic field intensity vector and magnetization vector are parallel. As a special case, the following two problems are considered: (1) the magnetoelastic stability of a ferromagnetic plate-strip in a homogeneous transverse magnetic field; (2) the stress–strain state of a ferromagnetic plane with a moving crack in a transverse magnetic field. It is shown that the modeling of magnetoelastic equations with a nonlinear law of magnetization provides qualitative and quantitative predictions on physical quantities including critical loads and stresses. In particular, it is shown that the critical magnetic field in plate stability problems found with the nonlinear law of magnetization is in better agreement with the experimental finding than the one found with a linear law. Furthermore, it is also shown that the stress concentration factor around a crack predicted with the nonlinear law of magnetization is more accurate than the one obtained with a linear counterpart. Numerical results are presented for above mentioned two problems and for various forms of nonlinear laws of magnetization.     

Large Eddy Simulation of the sound field of a round turbulent jet

In this paper we will use Large Eddy Simulation (LES) to obtain the flow field of a turbulent round jet at a Reynolds number based on the jet orifice velocity of 11000. In the simulations it is assumed that the flow field is incompressible. The acoustic field of the jet is calculated with help of the Lighthill acoustic analogy. The coupling between the flow solver and the acoustic solver is discussed in detail. The Mach number used in the acoustic calculation was equal to 0.6. It is shown that the decay of the jet centerline velocity and centerline rms are in good agreement with experimental data of [12]. Furthermore, it is shown that the influence of the LES modeling on the acoustic field is very small, if the dynamic subgrid model is used.     

Singular rigid/plastic solutions in anisotropic plasticity under plane strain conditions

Assuming a rigid plastic material model with arbitrary smooth yield criterion, it is shown that the plane strain solutions are singular in the vicinity of maximum friction surfaces. In particular, some components of the strain rate tensor and thus the equivalent strain rate approach infinity. It is also shown that the exact asymptotic representation of the solution near maximum friction surfaces depends on the shape of the yield contour in the Mohr stress plane.     

Simple Torsion of Isotropic, Hyperelastic, Incompressible Materials with Limiting Chain Extensibility

The purpose of this research is to investigate the simple torsion problem for a solid circular cylinder composed of isotropic hyperelastic incompressible materials with limiting chain extensibility. Three popular models that account for hardening at large deformations are examined. These models involve a strain-energy density which depends only on the first invariant of the Cauchy–Green tensor. In the limit as a polymeric chain extensibility tends to infinity, all of these models reduce to the classical neo-Hookean form. The main mechanical quantities of interest in the torsion problem are obtained in closed form. In this way, it is shown that the torsional response of all three materials is similar. While the predictions of the models agree qualitatively with experimental data, the quantitative agreement is poor as is the case for the neo-Hookean material. In fact, by using a global universal relation, it is shown that the experimental data cannot be predicted quantitatively by any strain-energy density which depends solely on the first invariant. It is shown that a modification of the strain energies to include a term linear in the second invariant can be used to remedy this defect. Whether the modified strain-energies, which reflect material hardening, are a feasible alternative to the classic Mooney–Rivlin model remains an open question which can be resolved only by large strain experiments. This revised version was published online in August 2006 with corrections to the Cover Date.     

Chaos and Entropy for Interval Maps

In this paper, various chaotic properties and their relationships for interval maps are discussed. It is shown that the proximal relation is an equivalence relation for any zero entropy interval map. The structure of the set of f-nonseparable pairs is well demonstrated and so is its relationship to Li-Yorke chaos. For a zero entropy interval map, it is shown that a pair is a sequence entropy pair if and only if it is f-nonseparable. Moreover, some equivalent conditions of positive entropy which relate to the number “3” are obtained. It is shown that for an interval map if it is topological null, then the pattern entropy of every open cover is of polynomial order, answering a question by Huang and Ye when the space is the closed unit interval.     

An analysis of overlapping Schwarz method for a weakly coupled system of singularly perturbed convection-diffusion equations

In this article, we have developed an overlapping Schwarz method for a weakly coupled system of convection-diffusion equations. The method splits the original domain into two overlapping subdomains. A hybrid difference scheme is proposed in which on the boundary layer region, we use the central finite difference scheme on a uniform mesh, whereas on the nonlayer region, we use the mid-point difference scheme on a uniform mesh. It is shown that the numerical approximations converge in the maximum norm to the exact solution. We have proved that, when appropriate subdomains are used, the method produces almost second-order convergence. Furthermore, it is shown that two iterations are sufficient to achieve the expected accuracy. Numerical examples are presented to support the theoretical results. The main advantage of this method used with the proposed scheme is that it reduces iteration counts very much and easily identifies in which iteration the Schwarz iterate terminates.     

Simulating the fluid forces and fluid-elastic instabilities of a clamped–clamped cylinder in turbulent axial flow

In this paper, the fluid forces and the dynamics of a flexible clamped–clamped cylinder in turbulent axial flow are computed numerically. In the presented numerical model, there is no need to tune parameters for each specific case or to obtain coefficients from experiments. The results are compared with the dynamics measured in experiments available in the literature. The specific case studied here consists of a silicone cylinder mounted in axial water flow. Computationally it is found that the cylinder loses stability first by buckling. The threshold for buckling is in quantitative agreement with experimental results and weakly nonlinear theory. At higher flow speed a fluttering motion is predicted, in agreement with experimental results. It is also shown that even a small misalignment between the flow and the structure can have a significant impact on the dynamical behavior. To provide insight in the results of these fluid–structure interaction simulations, forces are computed on rigid inclined and curved cylinders, showing the existence of two different flow regimes. Furthermore it is shown that the inlet turbulence state has a non-negligible effect on these forces and thus on the dynamics of the cylinder.     

Time and spatially resolved measurements of interfacial waves in vertical annular flow

Film thickness measurements have been performed in a vertical air/water annular flow in a pipe of 0.05 m diameter. A sensor has been built which allows to measure the film thickness evolution in time at 320 positions, such that the interface of the vertical annular flow can be reconstructed. The large-scale structures moving on the interface are described statistically, with a special attention to the disturbance waves. Probability density functions and mean statistics are given for the height, length, velocity, frequency and spatial distribution of the disturbance waves. In particular, it is shown that the disturbance waves are three-dimensional structures with large height fluctuations in the circumferential and axial direction, giving a meandering path between the maximum height around the circumference. It is also shown that the disturbance waves can flow with a slight inclination with respect to the axial direction. Finally, the disturbance waves are shown to be located randomly in space, within a Gamma distribution whose order only depends on the liquid superficial velocity. Due to the nature of the Gamma distribution, it could indicate that the spatial distribution of the disturbance waves results from a cascade of coalescence processes between the original disturbance waves on the film.     

Research on the hysteresis properties of unsteady aerodynamics about the oscillating wings

ＩｎｔｒｏｄｕｃｔｉｏｎＡｓｔｈｅｒｅｑｕｉｒｅｍｅｎｔｏｆｈｉｇｈｐｅｒｆｏｒｍａｎｃｅａｎｄｍａｎｅｕｖｅｒａｂｉｌｉｔｙ,ｔｈｅｎｅｘｔｇｅｎｅｒａｔｉｏｎｏｆｔｈｅｆｉｇｈｔｅｒａｉｒｃｒａｆｔｉｓｂｅｉｎｇｄｅｓｉｇｎｅｄｔｏｆｌｙａｎｄｂ...     

On the behavior of a three-dimensional fractional viscoelastic constitutive model

In this paper a three-dimensional isotropic fractional viscoelastic model is examined. It is shown that if different time scales for the volumetric and deviatoric components are assumed, the Poisson ratio is time varying function; in particular viscoelastic Poisson ratio may be obtained both increasing and decreasing with time. Moreover, it is shown that, from a theoretical point of view, one-dimensional fractional constitutive laws for normal stress and strain components are not correct to fit uniaxial experimental test, unless the time scale of deviatoric and volumetric are equal. Finally, the model is proved to satisfy correspondence principles also for the viscoelastic Poisson’s ratio and some issues about thermodynamic consistency of the model are addressed.     

On the instability of an axially moving elastic plate

Problems of stability of an axially moving elastic band travelling at constant velocity between two supports and experiencing small transverse vibrations are considered in a 2D formulation. The model of a thin elastic plate subjected to bending and tension is used to describe the bending moment and the distribution of membrane forces. The stability of the plate is investigated with the help of an analytical approach. In the frame of a general dynamic analysis, it is shown that the onset of instability takes place in the form of divergence (buckling). Then the static forms of instability are investigated, and critical regimes are studied as functions of geometric and mechanical problem parameters. It is shown that in the limit of a narrow strip, the 2D formulation reduces to the classical 1D model. In the limit of a wide band, there is a small but finite discrepancy between the results given by the 1D model and the full 2D formulation, where the discrepancy depends on the Poisson ratio of the material. Finally, the results are illustrated via numerical examples, and it is observed that the transverse displacement becomes localised in the vicinity of free boundaries.     

基于单元的静力区间有限元法

在许多工程问题中 ,结构参数和荷载具有某种程度的误差或不确定性。若不将它们定量化或模型化加以考虑 ,就不能作出合理的分析和设计。考虑到有限元法在科学界和工程界的广泛应用 ,本文以连续梁结构为例 ,建立了基于单元的静力区间有限元法。为了说明本方法的有效性 ,本文给出了一个数值例子 ,并把所得结果与文献 [1 2 ]进行了比较。     

Buckling behaviour of imperfect elastic and linearly viscoelastic structures

In a previous paper, the author has shown that the behaviour of imperfect elastic structures subjected to buckling forces could be predicted on the basis of the eigenvalues and eigenmodes. After a brief recall of these properties, it is first shown—in Appendix—that they extend to cases of spatial buckling like the buckling of flexure and torsion. Then, it is shown that the correspondence principle valid in first-order behaviour of linearly viscoelastic structures can be extended in full generality—by the use of Laplace transformation—to buckling problems, whatever be the constitutive equations of the material.Finally, several examples of eulerian buckling or flexural-torsional buckling of a bar and of buckling of a plate, are treated in detail.     

The visco-elastic Ekman layer

Solutions are obtained for the flow in an Ekman layer in a visco-elastic fluid using various constitutive relations. For some models (those which exhibit no shear thinning in simple shear flow) it is found that the solution is identical with the Newtonian solution for all Weissenberg numbers. For other models it is shown that the solution differs from its Newtonian counterpart and may cease to exist when a critical Weissenberg number is exceeded. In these cases, when the solution exists it is found that the visco-elastic Ekman layer is thinner than the Newtonian Ekman layer.     

Version of the linear elasticity theory with a structural parameter

It is shown that, in the planar case, the system of constitutive equations of the linear elasticity theory should contain five independent equations. In the classical theory, only three equations are formulated, while the other two are contained in implicit form in the postulate of diffeomorphism, which is the assumption of smoothness of the displacement field. A closed elasticity model is constructed without the assumption of diffeomorphism, and it contains a structural parameter having a dimension of length. It is shown that, in a static version, macrodeformations depend on stresses and second derivatives of stresses with respect to the coordinates, while there is dispersion of longitudinal and transverse waves in the dynamics.