Gauge principle and variational formulation for flows of an ideal fluid

A gauge principle is applied to mass flows of an ideal compressible fluid subject to Galilei transformation. A free-field Lagrangian defined at the outset is invariant with respeet to global SO(3) gauge transformations as well as Galilei transformations. The action principle leads to the equation of potential flows under constraint of a continuity equation. However, the irrotational flow is not invariant with respect to local SO(3) gauge transformations. According to the gauge principle, a gauge-covariant derivative is defined by introducing a new gauge field. Galilei invariance of the derivative requires the gauge field to coincide with the vorticity, i.e. the curl of the velocity field. A full gauge-covariant variational formulation is proposed on the basis of the Hamilton‘‘s principle and an assoicated Lagrangian. By means of an isentropic material variation taking into account individual particle motion, the Euler‘‘s equation of motion is derived for isentropic flows by using the covariant derivative. Noether‘‘s law associated with global SO(3) gauge invariance leads to the conservation of total angular momentum. In addition, the Lagrangian has a local symmetry of particle permutation which results in local conservation law equivalent to the vorticity equation.     

VARIATIONAL PRINCIPLE AND GENERALIZED VARIATIONAL PRINCIPLE IN HYDRODYNAMICS OF A CLASS OF NON－NEWTONIAN FLUID

ＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＡＮＤＧＥＮＥＲＡＬＩＺＥＤＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＩＮＨＹＤＲＯＤＹＮＡＭＩＣＳＯＦＡＣＬＡＳＳＯＦＮＯＮ－ＮＥＷＴＯＮＩＡＮＦＬＵＩＤＳｈｅｎＭｉｎ（沈敏）,ＳｕｎＱｉ－ｒｅｎ（孙其仁）（Ｓ...     

Variational formulation of the static Levinson beam theory

In this communication, we provide a consistent variational formulation for the static Levinson beam theory. First, the beam equations according to the vectorial formulation by Levinson are reviewed briefly. By applying the Clapeyron's theorem, it is found that the stresses on the lateral end surfaces of the beam are an integral part of the theory. The variational formulation is carried out by employing the principle of virtual displacements. As a novel contribution, the formulation includes the external virtual work done by the stresses on the end surfaces of the beam. This external virtual work contributes to the boundary conditions in such a way that artificial end effects do not appear in the theory. The obtained beam equations are the same as the vectorially derived Levinson equations. Finally, the exact Levinson beam finite element is developed.     

Averaged equations of ideal fluid turbulence

A system of equations for averaged velocity of turbulent motion of ideal fluid is derived from the assumption that this motion is ergodic. Several fluid flows are considered. A remarkable feature of the averaged equations is that in some cases they have a form of an eigenvalue problem similar to that for Schrődinger’s equation. The equations derived suggest that there are no universal equations for averaged velocity of turbulent flows.      

Variational formulation of a rod under torsional vibration for crack identification

A variational formulation is developed for the torsional vibration of a cylindrical shaft with a circumferential crack. The work is compared with existing methods. The Hu–Washizu–Barr variational formulation was used to develop the differential equation and the boundary conditions of the cracked rod. The general variational principle and the independent assumptions about the displacement, the momentum, the strain and the stress fields of the cracked rod, and the equations of motion for a uniform rod in torsional vibration, are derived. The crack was modelled as a continuous flexibility using the displacement field in the vicinity of the crack, found with fracture mechanics methods. Rayleigh quotient was used to approximate the natural frequencies of the cracked rod. Independent evaluations of crack identification methods in rotating shafts are reported and compared with methods using the continuous crack flexibility theory.     

Variational principles in hydrodynamics of a non-Newtonian fluid

I.IntroductionTileresearchworkonvariationalprincipleinhydrodynamicsstartedasearlyas1940's.Referencesl:~6]aretheoneswhichmakemoreimportantcolltributions.l-lowevcr,mostoftheseworksstudied11on-viscousflow,putemphasisontheproblemsofexternalllcld,andworkedontilebasisofBernoulli'sequation.Prot'essorChienWeizang11I.slesttlblishcdvariationalprincipleforhydrodynalllicpl-oblemofviscousfluid,thatistheprillciplcofmaxilnLllllpowerlosses,alldgeneralizedvariational(stationary)pl.inciplesoiltilebiLsisofNav…     

Variational formulation of the equivalent eigenstrain method with an application to a problem with radial eigenstrains

In this paper a variational formulation of the equivalent eigenstrain method is established. A functional of the Hashin–Shtrikman type is proposed such that the solution of the equivalent eigenstrain equation is a unique minimizer of the functional. Moreover, it is also shown that the equivalent eigenstrain equation is the Euler–Lagrange equation of the potential energy of the inclusions. An approximate solution of the equivalent eigenstrain equation is then found as a minimizer of the functional on a finite dimensional span of basic eigenstrains. Special attention is paid to possible symmetries of the problem. The variational formulation is illustrated by determination of effective linear elastic properties. In particular, material with a simple cubic microstructure is considered in detail. A solution for the polynomial radial basic eigenstrains approximation is found. In particular, for the homogeneous eigenstrain approximation, the effective moduli are derived in an exact closed form.     

Variational formulation of three-dimensional viscous free-surface flows: Natural imposition of surface tension boundary conditions

We present a new surface-intrinsic linear form for the treatment of normal and tangential surface tension boundary conditions in C⁰-geometry variational discretizations of viscous incompressible free-surface flows in three space dimensions. The new approach is illustrated by a finite (spectral) element unsteady Navier-Stokes analysis of the stability of a falling liquid film.     

A finite element formulation for laminar flow of a fluid with microstructure

A finite element formulation for the steady laminar flow of an incompressible fluid with microstructure has been developed. The particular fluids considered are commonly known as micropolar fluids, in which case suspended particulate microstructures are modelled by an ‘extended’ continuum formulation. The particle microspin is a new kinematic variable which is independent of the classical vorticity vector and thereby allows relative rotation between particles and the surrounding fluid. This formulation also gives rise to couple stresses in addition to classical force or traction stresses. The finite element formulation utilizes a variational approach and imposes conservation of mass through a penalty function. A general boundary condition for microspin has been incorporated whereby microspin at a solid boundary is constrained to be proportional to the fluid vorticity. The proportionality constant in this case can vary from zero to unity. Sample solutions are presented for fully developed flow through a straight tube and compared with an analytical solution. Results are also generated for flow through a constricted tube and compared with a Newtonian fluid solution.     

Nonstationary ideal incompressible fluid flows: Conditions of existence and uniqueness of solutions

The problem of formulating minimal conditions on input data that can guarantee the existence and uniqueness of solutions of the boundary value problems describing non-one-dimensional ideal incompressible fluid flow is considered using as an example the initial boundary value problem in a space-time cylinder constructed on a bounded flow domain with the nonpenetration condition on its boundary (which corresponds to fluid flow in a closed vessel). The existence problems are considered only for plane flows, and the uniqueness issues for three-dimensional flows as well. The required conditions are obtained in the form of conditions specifying that the vorticity belongs to definite functional Orlicz spaces. The results are compared with well-known results. Examples are given of admissible types of singularities for which the obtained results are valid, which is a physical interpretation of these results. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 130–145, July–August, 2008.     

Dynamic behaviour of an axially moving plate undergoing small cylindrical deformation submerged in axially flowing ideal fluid

The out-of-plane dynamic response of a moving plate, travelling between two rollers at a constant velocity, is studied, taking into account the mutual interaction between the vibrating plate and the surrounding, axially flowing ideal fluid. Transverse displacement of the plate (assumed cylindrical) is described by an integro-differential equation that includes a local inertia term, Coriolis and centrifugal forces, the aerodynamic reaction of the external medium, the vertical projection of membrane tension, the bending resistance, and external perturbation forces. In the two-dimensional model thus set up, the aerodynamic reaction is found analytically as a functional of the cylindrical displacement, using the techniques of complex analysis. The resulting integro-differential problem is discretized in space with the Fourier-Galerkin method, and integrated in time with the diagonalization method. Examples are computed with physical parameters corresponding to air and some paper materials. The effects of the surrounding fluid on the critical velocity and first natural frequency are investigated, for stationary air, for an air mass moving with the plate, and for some arbitrary axial fluid velocities. The obtained results are applicable for both an ideal membrane and a plate with nonzero bending rigidity.     

On an attempt to simplify the Quartapelle–Napolitano approach to computation of hydrodynamic forces in open flows

In this paper we are interested in the Quartapelle–Napolitano approach to calculation of forces in viscous incompressible flows in exterior domains. We study the possibility of deriving a simpler formulation of this approach which might lead to a more convenient expression for the hydrodynamic force, but conclude that such a simplification is, within the family of approaches considered, impossible. This shows that the original Quartapelle–Napolitano formula is in fact “optimal” within this class of approaches.     

Development of a compressive surface capturing formulation for modelling free‐surface flow by using the volume‐of‐fluid approach

With the aim of accurately modelling free‐surface flow of two immiscible fluids, this study presents the development of a new volume‐of‐fluid free‐surface capturing formulation. By building on existing volume‐of‐fluid approaches, the new formulation combines a blended higher resolution scheme with the addition of an artificial compressive term to the volume‐of‐fluid equation. This reduces the numerical smearing of the interface associated with explicit higher resolution schemes while limiting the contribution of the artificial compressive term to ensure the integrity of the interface shape is maintained. Furthermore, the computational efficiency of the the higher resolution scheme is improved through the reformulation of the normalised variable approach and the implementation of a new higher resolution blending function. The volume‐of‐fluid equation is discretised via an unstructured vertex‐centred finite volume method and solved via a Jacobian‐type dual time‐stepping approach. Copyright © 2012 John Wiley & Sons, Ltd.     

Onset of self-oscillations upon the loss of stability of spatially periodic two-dimensional viscous fluid flows relative to long-wave perturbations

The long-wave asymptotics of the secondary flow that arises after a steady, spatially periodic flow loses stability are studied when one of the periods tends to infinity and the rate of base flow along the longer period is equal to zero. It is shown that if certain non-degeneracy conditions are satisfied, then from the base solution a self-oscillatory regime branches off and both hard and soft stability loss is possible. For the leading terms of the asymptotics explicit formulas are obtained. Examples of self-oscillations calculated for specific flows are presented and the behavior of the fluid particle trajectories in the self-oscillatory regime branching off from the base flow is investigated.     

Numerical simulations of inviscid and viscous free‐surface flows with application to nonlinear water waves

The purpose of the present study is to establish a numerical model appropriate for solving inviscid/viscous free‐surface flows related to nonlinear water wave propagation. The viscous model presented herein is based on the Navier–Stokes equations, and the free‐surface is calculated through an arbitrary Lagrangian–Eulerian streamfunction‐vorticity formulation. The streamfunction field is governed by the Poisson equation, and the vorticity is obtained on the basis of the vorticity transport equation. For computing the inviscid flow the Laplace streamfunction equation is used. These equations together with the respective (appropriate) fully nonlinear free‐surface boundary conditions are solved using a finite difference method. To demonstrate the model feasibility, in the present study we first simulate collision processes of two solitary waves of different amplitudes, and compute the phenomenon of overtaking of such solitary waves. The developed model is subsequently applied to calculate (both inviscid and the viscous) flow field, as induced by passing of a solitary wave over submerged rectangular structures and rigid ripple beds. Our study provides a reasonably good understanding of the behavior of (inviscid/viscous) free‐surface flows, within the framework of streamfunction‐vorticity formulation. The successful simulation of the above‐mentioned test cases seems to suggest that the arbitrary Lagrangian–Eulerian/streamfunction‐vorticity formulation is a potentially powerful approach, capable of effectively solving the fully nonlinear inviscid/viscous free‐surface flow interactions. Copyright © 2009 John Wiley & Sons, Ltd.     

On a difficulty in the formulation of initial and boundary conditions for eigenfunction expansion solutions for the start-up of fluid flow

Most mathematics and engineering textbooks describe the process of “subtracting off” the steady state of a linear parabolic partial differential equation as a technique for obtaining a boundary-value problem with homogeneous boundary conditions that can be solved by separation of variables (i.e., eigenfunction expansions). While this method produces the correct solution for the start-up of the flow of, e.g., a Newtonian fluid between parallel plates, it can lead to erroneous solutions to the corresponding problem for a class of non-Newtonian fluids. We show that the reason for this is the non-rigorous enforcement of the start-up condition in the textbook approach, which leads to a violation of the principle of causality. Nevertheless, these boundary-value problems can be solved correctly using eigenfunction expansions, and we present the formulation that makes this possible (in essence, an application of Duhamel's principle). The solutions obtained by this new approach are shown to agree identically with those obtained by using the Laplace transform in time only, a technique that enforces the proper start-up condition implicitly (hence, the same error cannot be committed).     

MHD flows of a second grade fluid between two side walls perpendicular to a plate through a porous medium

Exact analytical solutions for magnetohydrodynamic (MHD) flows of an incompressible second grade fluid in a porous medium are developed. The modified Darcy's law for second grade fluid has been used in the flow modelling. The Hall effect is taken into account. The exact solutions for the unsteady flow induced by the time-dependent motion of a plane wall between two side walls perpendicular to the plane has been constructed by means of Fourier sine transforms. The similar solutions for a Newtonian fluid, performing the same motion, appear as limiting cases of the solutions obtained here. The influence of various parameters of interest on the velocity and shear stress at the bottom wall has been shown and discussed through several graphs. A comparison between a Newtonian and a second grade fluids is also made.     

Viscometric flows of polymeric melts treated according to the rational theory of Eckart's continuum. I. General theory

The purely rational theory of Eckart continua (i.e. elastic bodies with a variable relaxed state) is applied to viscometric flows of polymeric melts. The main assumptions are thermodynamic non-interaction of inelastic behaviour and of non-elastic stress, as well as elastic isotropy. After establishing the time-dependent differential equations of viscometric flow, these equations are simplified to a set of algebraic equations describing steady-state flow. From this we deduce two general equations connecting the three elastic steady-state viscometric functions which do not depend upon the elastic behaviour. The law of rubber elasticity used in this paper is described in the Appendix.     

Numerical modelling of fluid and solid thermomechanics in additive manufacturing by powder-bed fusion: Continuum and level set formulation applied to track- and part-scale simulations

The thermomechanical analysis of powder-bed fusion using a laser beam is simulated in both meso- and macroscales within a framework combining continuum assumption and level-set formulation. The mesoscale simulation focuses on laser interaction with the powder bed, and on subsequent melting and solidification. Modelling is conducted at the scale of material deposition, in which powder-bed fusion, hydrodynamics in the melt pool, and thermal stress are simulated. The macroscale model focuses on part construction and post-deposition. During construction, by contrast with the mesoscale approach, the fluid flow in the fusion zone is ignored and material addition is simplified by modelling it at the scale of entire layers, or fractions of layers. The modelling of the energy input is adapted accordingly. This thermomechanical model addresses heat exchange, residual stress, and distortion at the part's scale. In both approaches, adaptive remeshing is applied, providing a good compromise between the needs to provide accurate prediction and maintaining sustainable computation times.