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

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

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 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.     

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.     

Unified one‐fluid formulation for incompressible flexible solids and multiphase flows: Application to hydrodynamics using the immersed structural potential method (ISPM)

In this paper, we present a two‐dimensional computational framework for the simulation of fluid‐structure interaction problems involving incompressible flexible solids and multiphase flows, further extending the application range of classical immersed computational approaches to the context of hydrodynamics. The proposed method aims to overcome shortcomings such as the restriction of having to deal with similar density ratios among different phases or the restriction to solve single‐phase flows. First, a variation of classical immersed techniques, pioneered with the immersed boundary method (IBM), is presented by rearranging the governing equations, which define the behaviour of the multiple physics involved. The formulation is compatible with the “one‐fluid” formulation for two‐phase flows and can deal with large density ratios with the help of an anisotropic Poisson solver. Second, immersed deformable structures and fluid phases are modelled in an identical manner except for the computation of the deviatoric stresses. The numerical technique followed in this paper builds upon the immersed structural potential method developed by the authors, by adding a level set–based method for the capturing of the fluid‐fluid interfaces and an interface Lagrangian‐based meshless technique for the tracking of the fluid‐structure interface. The spatial discretisation is based on the standard marker‐and‐cell method used in conjunction with a fractional step approach for the pressure/velocity decoupling, a second‐order time integrator, and a fixed‐point iterative scheme. The paper presents a wide d range of two‐dimensional applications involving multiphase flows interacting with immersed deformable solids, including benchmarking against both experimental and alternative numerical schemes.     

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.     

真实虚拟流方法在多介质可压缩流动模拟中的应用

为了克服原始虚拟流方法(ghost fluid method,GFM)在处理激波与大密度比流体-流体（气-水）界面相互作用时遇到的困难,采用真实虚拟流法(real ghost fluid method,RGFM)处理流体界面附近的虚拟点,结合HLLC(Harten-Lax-Van Leer with contact discontinuities)格式求解Euler方程,采用五阶WENO(weighted essentially nonoscillatory)格式求解level set输运方程。通过一维和二维算例的物质界面捕捉研究,证明RGFM在处理小密度比界面问题时优于GFM,同时RGFM还可用于求解激波与大密度比物质界面相互作用问题。计算表明,将RGFM引入到本文算法中,可精确捕捉到激波与界面（气-气、气-水界面）相互作用的变化细节,包括大密度比界面的剧烈变形和破碎,并具有较高的计算分辨率。     

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.