Blowup Criterion for Viscous Baratropic Flows with Vacuum States

We prove that the maximum norm of the deformation tensor of velocity gradients controls the possible breakdown of smooth(strong) solutions for the 3-dimensional (3D) barotropic compressible Navier-Stokes equations. More precisely, if a solution of the 3D barotropic compressible Navier-Stokes equations is initially regular and loses its regularity at some later time, then the loss of regularity implies the growth without bound of the deformation tensor as the critical time approaches. Our result is the same as Ponce’s criterion for 3-dimensional incompressible Euler equations (Ponce in Commun Math Phys 98:349–353, 1985). In addition, initial vacuum states are allowed in our cases.     

Remarks on the breakdown of smooth solutions for the 3-D Euler equations

The authors prove that the maximum norm of the vorticity controls the breakdown of smooth solutions of the 3-D Euler equations. In other words, if a solution of the Euler equations is initially smooth and loses its regularity at some later time, then the maximum vorticity necessarily grows without bound as the critical time approaches; equivalently, if the vorticity remains bounded, a smooth solution persists.Partially supported by O.N.R. Contract No. N00014-76-C-0316 and N.S.F. Grant No. MCS-81-01639Partially supported by N.S.F. Grant No. MCS-82-00171Partially supported by N.S.F. Grant No. MCS-81-02360     

Vanishing Viscosity Solutions of the Compressible Euler Equations with Spherical Symmetry and Large Initial Data

We are concerned with spherically symmetric solutions of the Euler equations for multidimensional compressible fluids, which are motivated by many important physical situations. Various evidences indicate that spherically symmetric solutions of the compressible Euler equations may blow up near the origin at a certain time under some circumstance. The central feature is the strengthening of waves as they move radially inward. A longstanding open, fundamental problem is whether concentration could be formed at the origin. In this paper, we develop a method of vanishing viscosity and related estimate techniques for viscosity approximate solutions, and establish the convergence of the approximate solutions to a global finite-energy entropy solution of the isentropic Euler equations with spherical symmetry and large initial data. This indicates that concentration is not formed in the vanishing viscosity limit, even though the density may blow up at a certain time. To achieve this, we first construct global smooth solutions of appropriate initial-boundary value problems for the Euler equations with designed viscosity terms, approximate pressure function, and boundary conditions, and then we establish the strong convergence of the viscosity approximate solutions to a finite-energy entropy solution of the Euler equations.     

On the Convergence Rate of the Euler-α, an Inviscid Second-Grade Complex Fluid, Model to the Euler Equations

We study the convergence rate of the solutions of the incompressible Euler-α, an inviscid second-grade complex fluid, equations to the corresponding solutions of the Euler equations, as the regularization parameter α approaches zero. First we show the convergence in H ^s , s>n/2+1, in the whole space, and that the smooth Euler-α solutions exist at least as long as the corresponding solution of the Euler equations. Next we estimate the convergence rate for two-dimensional vortex patch with smooth boundaries.     

Post-Newtonian perturbed theory of elastic astronomical bodies in rotating spherical coordinates

In this paper the dynamical equations for an elastic deformable body in the first post-Newtonian approximation of Einstein theory of gravity are derived in rotating spherical coordinates. The unperturbed rotating body (the relaxed ground state) is described as uniformly rotating, stationary and axisymmetric configuration in an asymptotically flat space-time manifold. Deviations from the equilibrium configuration are described by means of a displacement field. By making use of the schemes developed by Damour, Soffel and Xu, and by Carter and Quintana we calculate the post-Newtonian Lagrangian strain tensor and symmetric trace-free shear tensor. Considering the Euler variations of Einstein's energy-momentum conservation law, we derive the post- Newtonian energy equation and Euler equations of elastic deformable bodies in rotating spherical coordinates.     

Classification of Equilibrium Solutions of the Cometary Flow Equation and Explicit Solutions of the Euler Equations for Monatomic Ideal Gases

The set of smooth equilibrium solutions of a kinetic model for cometary flows is split into equivalence classes according to similarity transformations. For each equivalence class in the two- and three-dimensional cases a normal form is computed. Each such equilibrium solution gives rise to an explicit solution of the compressible Euler equations for monatomic gases. The set of these solutions is discussed with special emphasis on solutions containing vacuum regions.     

欧拉运动学方程的另一种推导方法

欧拉运动学方程的推导,关键是确定本体坐标系与空间坐标系之间的转换关系.通过引进方位矢量可以较轻易地得到反映两坐标系转换关系的转换张量,从而能清晰直观地推导出欧拉运动学方程,而这是有别于传统的推导方法的.     

Blowup of Smooth Solutions for Relativistic Euler Equations

We study the singularity formation of smooth solutions of the relativistic Euler equations in (3 + 1)-dimensional spacetime for both finite initial energy and infinite initial energy. For the finite initial energy case, we prove that any smooth solution, with compactly supported non-trivial initial data, blows up in finite time. For the case of infinite initial energy, we first prove the existence, uniqueness and stability of a smooth solution if the initial data is in the subluminal region away from the vacuum. By further assuming the initial data is a smooth compactly supported perturbation around a non-vacuum constant background, we prove the property of finite propagation speed of such a perturbation. The smooth solution is shown to blow up in finite time provided that the radial component of the initial ``generalized' momentum is sufficiently large.     

Ray integrals of integrated photoelasticity at weak optical anisotropy

As initial data, optical polarization tomography of tensor and vector fields uses ray integrals the values of which are determined from polarization measurements. Most algorithms of integrated photoelasticity are based on the linear approximation of the solution of the birefringence equations. In this work, the approximation of smooth rotation of quasi-principal axes of the polarization tensor is proposed for use in approximating values of ray integrals.     

Note on loss of regularity for solutions of the 3—D incompressible euler and related equations

One of the central problems in the mathematical theory of turbulence is that of breakdown of smooth (indefinitely differentiable) solutions to the equations of motion. In 1934 J. Leray advanced the idea that turbulence may be related to the spontaneous appearance of singularities in solutions of the 3—D incompressible Navier-Stokes equations. The problem is still open. We show in this report that breakdown of smooth solutions to the 3—D incompressible slightly viscous (i.e. corresponding to high Reynolds numbers, or highly turbulent) Navier-Stokes equations cannot occur without breakdown in the corresponding solution of the incompressible Euler (ideal fluid) equation. We prove then that solutions of distorted Euler equations, which are equations closely related to the Euler equations for short term intervals, do breakdown.Sponsored by the United States Army under Contract No. DAAG29-80-C-0041, and partially supported by the National Science Foundation under Grant No. MCS-82-01599     

On the Spectral Dynamics of the Deformation Tensor and New A Priori Estimates for the 3D Euler Equations

In this paper we study the dynamics of eigenvalues of the deformation tensor for solutions of the 3D incompressible Euler equations. Using the evolution equation of the L² norm of spectra, we deduce new a priori estimates of the L² norm of vorticity. As an immediate corollary of the estimate we obtain a new sufficient condition of L² norm control of vorticity. We also obtain decay in time estimates of the ratios of the eigenvalues. In the remarks we discuss what these estimates suggest in the study of searching initial data leading to possible finite time singularities. We find that the dynamical behaviors of L² norm of vorticity are controlled completely by the second largest eigenvalue of the deformation tensor. Part of this work was done while the author was visiting CSCAMM, University of Maryland, USA. The author would like to thank to Professor E. Tadmor for his hospitality     

Zero Viscosity Limit for Analytic Solutions of the Navier-Stokes Equation on a Half-Space.¶ II. Construction of the Navier-Stokes Solution

This is the second of two papers on the zero-viscosity limit for the incompressible Navier-Stokes equations in a half-space in either 2D or 3D. Under the assumption of analytic initial data, we construct solutions of Navier-Stokes for a short time which is independent of the viscosity. The Navier-Stokes solution is constructed through a composite asymptotic expansion involving the solutions of the Euler and Prandtl equations, which were constructed in the first paper, plus an error term. This shows that the Navier-Stokes solution goes to an Euler solution outside a boundary layer and to a solution of the Prandtl equations within the boundary layer. The error term is written as a sum of first order Euler and Prandtl corrections plus a further error term. The equation for the error term is weakly nonlinear; its linear part is the time dependent Stokes equation. This error equation is solved by inversion of the Stokes equation, through expressing the solution as a regular (Euler-like) part plus a boundary layer (Prandtl-like) part. The main technical tool in this analysis is the Abstract Cauchy-Kowalewski Theorem. Received: 5 September 1996 / Accepted: 14 July 1997     

Study of nonlinear system stability using eigenvalue analysis: Gyroscopic motion

General computational multibody system (MBS) algorithms allow for the linearization of the highly nonlinear equations of motion at different points in time in order to obtain the eigenvalue solution. This eigenvalue solution of the linearized equations is often used to shed light on the system stability at different configurations that correspond to different time points. Different MBS algorithms, however, employ different sets of orientation coordinates, such as Euler angles and Euler parameters, which lead to different forms of the dynamic equations of motion. As a consequence, the forms of the linearized equations and the eigenvalue solution obtained strongly depend on the set of orientation coordinates used. This paper addresses this fundamental issue by examining the effect of the use of different orientation parameters on the linearized equations of a gyroscope. The nonlinear equations of motion of the gyroscope are formulated using two different sets of orientation parameters: Euler angles and Euler parameters. In order to obtain a set of linearized equations that can be used to define the eigenvalue solution, the algebraic equations that describe the MBS constraints are systematically eliminated leading to a nonlinear form of the equations of motion expressed in terms of the system degrees of freedom. Because in MBS applications the generalized forces can be highly nonlinear and can depend on the velocities, a state space formulation is used to solve the eigenvalue problem. It is shown in this paper that the independent state equations formulated using Euler angles and Euler parameters lead to different eigenvalue solutions. This solution is also different from the solution obtained using a form of the Newton-Euler matrix equation expressed in terms of the angular accelerations and angular velocities. A time-domain solution of the linearized equations is also presented in order to compare between the solutions obtained using two different sets of orientation parameters and also to shed light on the important issue of using the eigenvalue analysis in the study of MBS stability. The validity of using the eigenvalue analysis based on the linearization of the nonlinear equations of motion in the study of the stability of railroad vehicle systems, which have known critical speeds, is examined. It is shown that such an eigenvalue analysis can lead to wrong conclusions regarding the stability of nonlinear systems.     

Computation of aeolian tone from a circular cylinder using source models

The generation of aeolian tones from a two-dimensional circular cylinder situated in a uniform cross-flow is investigated. The major emphasis here is placed on identifying the important noise generation mechanisms. Acoustic-viscous splitting techniques are utilized to compute modelled acoustic source terms and their corresponding acoustic fields. The incompressible Reynolds averaged Navier-Stokes equation is used to compute the near-field viscous flow solution, from which modelled acoustic source terms are extracted based on an approximation to the Lighthill’s stress tensor. Acoustic fields are then computed with an acoustic solver to solve the linearized Euler equations forced by the modelled source terms. Computations of the acoustic field based on the approximated Lighthill’s stress tensor are shown to be in good agreement with those computed from the surface dipole sources obtained using Curle’s solution to the acoustic analogy. It is shown in this paper that the stress tensor source term in the streamwise direction makes a comparable, but slightly larger contribution to the overall radiated field, compared with that due to the stress tensor in the direction normal to the mean flow. In addition, it is shown that shear sources, which arise due to the interaction between the fluctuating velocity and the background steady mean velocity, make the greatest contribution to the acoustic field, while the self-noise sources, which represents the interaction between the fluctuating velocities, is shown to be comparably negligible.     

An alternative f(R, T) gravity theory and the dark energy problem

Recently, a generalized gravity theory was proposed by Harko et al. where the Lagrangian density is an arbitrary function of the Ricci scalar R and the trace of the stress-energy tensor T, known as F(R,T) gravity. In their derivation of the field equations, they have not considered conservation of the stress-energy tensor. In the present work, we have shown that a part of the arbitrary function f(R,T) can be determined if we take into account of the conservation of stress-energy tensor, although the form of the field equations remain similar. For homogeneous and isotropic model of the universe the field equations are solved and corresponding cosmological aspects has been discussed. Finally, we have studied the energy conditions in this modified gravity theory both generally and a particular case of perfect fluid with constant equation of state.     

Remarks on the Helicity of the 3-D Incompressible Euler Equations

In this note we obtain a sufficient condition on the regularity of the weak solutions to guarantee the conservation of helicity for the 3-D incompressible Euler equations. As a corollary we obtain a lower bound of the vorticities for a weak solution of the Euler equations.     

Fluid dynamic limits of kinetic equations. I. Formal derivations

The connection between kinetic theory and the macroscopic equations of fluid dynamics is described. In particular, our results concerning the incompressible Navier-Stokes equations are based on a formal derivation in which limiting moments are carefully balanced rather than on a classical expansion such as those of Hilbert or Chapman-Enskog. The moment formalism shows that the limit leading to the incompressible Navier-Stokes equations, like that leading to the compressible Euler equations, is a natural one in kinetic theory and is contrasted with the systematics leading to the compressible Navier-Stokes equations. Some indications of the validity of these limits are given. More specifically, the connection between the DiPerna-Lions renormalized solution of the classical Boltzmann equation and the Leray solution of the Navier-Stokes equations is discussed.This paper is dedicated to Joel Lebowitz on his 60th-birthday.     

On the Partial Regularity of a 3D Model of the Navier-Stokes Equations

We study the partial regularity of a 3D model of the incompressible Navier-Stokes equations which was recently introduced by the authors in [11]. This model is derived for axisymmetric flows with swirl using a set of new variables. It preserves almost all the properties of the full 3D Euler or Navier-Stokes equations except for the convection term which is neglected in the model. If we add the convection term back to our model, we would recover the full Navier-Stokes equations. In [11], we presented numerical evidence which seems to support that the 3D model develops finite time singularities while the corresponding solution of the 3D Navier-Stokes equations remains smooth. This suggests that the convection term play an essential role in stabilizing the nonlinear vortex stretching term. In this paper, we prove that for any suitable weak solution of the 3D model in an open set in space-time, the one-dimensional Hausdorff measure of the associated singular set is zero. The partial regularity result of this paper is an analogue of the Caffarelli-Kohn-Nirenberg theory for the 3D Navier-Stokes equations.     

Mach reflection for the two-dimensional Burgers equation

We study shock reflection for the two 2D Burgers equation. This model equation is an asymptotic limit of the Euler equations, and retains many of the features of the full equations. A von Neumann type analysis shows that the 2D Burgers equation has detachment, sonic, and Crocco points in complete analogy with gas dynamics. Numerical solutions support the detachment/sonic criterion for transition from regular to Mach reflection. There is also strong numerical evidence that the reflected shock in the 2D Burgers Mach reflection forms a smooth wave near the Mach stem, as proposed by Colella and Henderson in their study of the Euler equations.