Global Well–Posedness of the 3<Emphasis Type="Italic">D</Emphasis> Primitive Equations with Partial Vertical Turbulence Mixing Heat Diffusion

The three–dimensional incompressible viscous Boussinesq equations, under the assumption of hydrostatic balance, govern the large scale dynamics of atmospheric and oceanic motion, and are commonly called the primitive equations. To overcome the turbulence mixing a partial vertical diffusion is usually added to the temperature advection (or density stratification) equation. In this paper we prove the global regularity of strong solutions to this model in a three-dimensional infinite horizontal channel, subject to periodic boundary conditions in the horizontal directions, and with no-penetration and stress-free boundary conditions on the solid, top and bottom boundaries. Specifically, we show that short time strong solutions to the above problem exist globally in time, and that they depend continuously on the initial data.     

3D Stochastic Primitive Equations of the Large-Scale Ocean: Global Well-Posedness and Attractors

In this paper, we consider the global well-posedness and long-time dynamics for the three-dimensional viscous primitive equations describing the large-scale oceanic motion under a random forcing, which is an additive white in time noise. We firstly prove the existence and uniqueness of global strong solutions to the initial boundary value problem for the stochastic primitive equations. Subsequently, by studying the asymptotic behavior of strong solutions, we obtain the existence of random attractors for the corresponding random dynamical system.     

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.     

Incompressible Limit of the Compressible Magnetohydrodynamic Equations with Periodic Boundary Conditions

This paper is concerned with the incompressible limit of the compressible magnetohydrodynamic equations with periodic boundary conditions. It is rigorously shown that the weak solutions of the compressible magnetohydrodynamic equations converge to the strong solution of the viscous or inviscid incompressible magnetohydrodynamic equations as long as the latter exists both for the well-prepared initial data and general initial data. Furthermore, the convergence rates are also obtained in the case of the well-prepared initial data.     

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.     

On Blowup of Classical Solutions to the Compressible Navier-Stokes Equations

In this paper, we study the finite time blow up of smooth solutions to the Compressible Navier-Stokes system when the initial data contain vacuums. We prove that any classical solutions of viscous compressible fluids without heat conduction will blow up in finite time, as long as the initial data has an isolated mass group (see Definition 2.2). The results hold regardless of either the size of the initial data or the far fields being vacuum or not. This improves the blowup results of Xin (Comm Pure Appl Math 51:229–240, 1998) by removing the crucial assumptions that the initial density has compact support and the smooth solution has finite total energy. Furthermore, the analysis here also yields that any classical solutions of viscous compressible fluids without heat conduction in bounded domains or periodic domains will blow up in finite time, if the initial data have an isolated mass group satisfying some suitable conditions.     

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.     

Global Well-Posedness Issues for the Inviscid Boussinesq System with Yudovich’s Type Data

The present paper is dedicated to the study of the global existence for the inviscid two-dimensional Boussinesq system. We focus on finite energy data with bounded vorticity and we find out that, under quite a natural additional assumption on the initial temperature, there exists a global unique solution. No smallness conditions are imposed on the data. The global existence issues for infinite energy initial velocity, and for the Bénard system are also discussed.     

Finite-Wavelength Stability¶of Capillary-Gravity Solitary Waves

We consider the Euler equations describing nonlinear waves on the free surface of a two-dimensional inviscid, irrotational fluid layer of finite depth. For large surface tension, Bond number larger than 1/3, and Froude number close to 1, the system possesses a one-parameter family of small-amplitude, traveling solitary wave solutions. We show that these solitary waves are spectrally stable with respect to perturbations of finite wave-number. In particular, we exclude possible unstable eigenvalues of the linearization at the soliton in the long-wavelength regime, corresponding to small frequency, and unstable eigenvalues with finite but bounded frequency, arising from non-adiabatic interaction of the infinite-wavelength soliton with finite-wavelength perturbations. Received: 7 February 2001 / Accepted: 6 October 2001     

An explicit construction of a class of suspensions and autonomous differential equations for diffeomorphisms in the plane

From a large class of diffeomorphisms in the plane, which are known to produce chaotic dynamics, we explicitly construct their continuous suspension on a three dimensional cylinder. This suspension is smooth (C ¹) and can be characterized by the choice of two smooth functions on the unit interval, which have to fulfill certain boundary conditions. For the case of entire Cremona transformations, we are able to construct the corresponding autonomous differential equations of the flow explicitly. Thus it is possible to relate properties of discrete maps to those of ordinary differential equations in a quantitative manner. Furthermore, our construction makes it possible to study the exact solutions of chaotic differential-equations directly.     

Initial-boundary value problems for a class of nonlinear thermoelastic plate equations

This paper studies initial-boundary value problems for a class of nonlinear thermoelastic plate equations. Under some certain initial data and boundary conditions, it obtains an existence and uniqueness theorem of global weak solutions of the nonlinear thermoelstic plate equations, by means of the Galerkin method. Moreover, it also proves the existence of strong and classical solutions.     

The Complex Geometry of Weak Piecewise Smooth Solutions of Integrable Nonlinear PDE's¶of Shallow Water and Dym Type

An extension of the algebraic-geometric method for nonlinear integrable PDE's is shown to lead to new piecewise smooth weak solutions of a class of N-component systems of nonlinear evolution equations. This class includes, among others, equations from the Dym and shallow water equation hierarchies. The main goal of the paper is to give explicit theta-functional expressions for piecewise smooth weak solutions of these nonlinear PDE's, which are associated to nonlinear subvarieties of hyperelliptic Jacobians. The main results of the present paper are twofold. First, we exhibit some of the special features of integrable PDE's that admit piecewise smooth weak solutions, which make them different from equations whose solutions are globally meromorphic, such as the KdV equation. Second, we blend the techniques of algebraic geometry and weak solutions of PDE's to gain further insight into, and explicit formulas for, piecewise-smooth finite-gap solutions. The basic technique used to achieve these aims is rather different from earlier papers dealing with peaked solutions. First, profiles of the finite-gap piecewise smooth solutions are linked to certain finite dimensional billiard dynamical systems and ellipsoidal billiards. Second, after reducing the solution of certain finite dimensional Hamiltonian systems on Riemann surfaces to the solution of a nonstandard Jacobi inversion problem, this is resolved by introducing new parametrizations. Amongst other natural consequences of the algebraic-geometric approach, we find finite dimensional integrable Hamiltonian dynamical systems describing the motion of peaks in the finite-gap as well as the limiting (soliton) cases, and solve them exactly. The dynamics of the peaks is also obtained by using Jacobi inversion problems. Finally, we relate our method to the shock wave approach for weak solutions of wave equations by determining jump conditions at the peak location. Received: 16 February 1999 / Accepted: 10 April 2001     

Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

Convergence of the Vanishing Viscosity Approximation for Superpositions of Confined Eddies

A confined eddy is a circularly symmetric flow with vorticity of compact support and zero net circulation. Confined eddies with disjoint supports can be superimposed to generate stationary weak solutions of the two-dimensional incompressible inviscid Euler equations. In this work, we consider the unique weak solution of the two-dimensional incompressible Navier-Stokes equations having a disjoint superposition of very singular confined eddies as the initial datum. We prove the convergence of these weak solutions back to the initial configuration, as the Reynolds number goes to infinity. This implies that the stationary superposition of confined eddies with disjoint supports is the unique physically correct weak solution of the two-dimensional incompressible Euler equations.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

Wave propagation according to higher order field equations I

A detailed study is made of wave propagation according to a sixth-order partial differential equation with complex masses proposed by Swieca and Marques, which presents a kind of generalized Klein-Gordon equation. The choice of definite Green's functions in the corresponding Yang-Feldman integral equation corresponds to a certain choice of boundary conditions for the allowed solutions of the corresponding partial differential equation. The advanced and retarded Green's functions used possess the anomalous feature of having non-zero values in the neighbourhoods of those, past or future parts of the light cone, for which traditional advanced and retarded Green's functions are zero. However, it is shown that a suitable averaging procedure provides the possibility of defining sets of functions, such that solutions of the Yang-Feldman equations belonging to this set possess the property that the future behaviour of the solution is determined by its asymptotic initial conditions. Certain features of the wave propagation, according to the equations considered, can be usefully compared with the properties of the solutions of the ordinary differential equation - and corresponding integral equation - which represents the equation of motion of a charged particle including the force for radiation reaction. The particle then has a certain “size”. Analogously the “non-local field equations” have solutions characterized by a certain “fundamental length” indicating the space-time distances for which averaging occurs. The admitted solutions of the field equations seem to represent a relativistic field with a “finite a number of degrees of freedom” within a finite volume.     

Singularities of Euler Flow? Not Out of the Blue!

Does three-dimensional incompressible Euler flow with smooth initial conditions develop a singularity with infinite vorticity after a finite time? This blowup problem is still open. After briefly reviewing what is known and pointing out some of the difficulties, we propose to tackle this issue for the class of flows having analytic initial data for which hypothetical real singularities are preceded by singularities at complex locations. We present some results concerning the nature of complex space singularities in two dimensions and propose a new strategy for the numerical investigation of blowup.     

A pseudospectral method for solution of the three-dimensional incompressible Navier-Stokes equations

A new Chebyshev pseudospectral technique (based on the projection method that was previously applied by the authors to the solution of two-dimensional incompressible Navier-Stokes equations in primitive variables for nonperiodic boundary conditions) is extended to solve the three-dimensional Navier-Stokes equations. The crucial point of the method is the requirement that the continuity equation be satisfied everywhere in the domain, on the boundaries as well as in the interior. The key feature of the work presented in this paper is that the computer storage requirements of the full matrix inversion resulting from direct solution of the pressure Poisson equation in three dimensions is greatly reduced by considering an eigenfunction decomposition. The method was tested on a two-dimensional driven cavity flow and the results were compared with those of the most accurate finite-difference calculation. The three-dimensional driven cavity flow was then calculated at the same Reynolds numbers as the two-dimensional cases, i.e., Re = 100, 400, and 1000. In the calculated results, three-dimensional boundary effects were observed in all cases and became more apparent with increasing Reynolds number.     

Solving second-order nonlinear evolution partial differential equations using deep learning

Solving nonlinear evolution partial differential equations has been a longstanding computational challenge. In this paper, we present a universal paradigm of learning the system and extracting patterns from data generated from experiments. Specifically, this framework approximates the latent solution with a deep neural network, which is trained with the constraint of underlying physical laws usually expressed by some equations. In particular, we test the effectiveness of the approach for the Burgers' equation used as an example of second-order nonlinear evolution equations under different initial and boundary conditions. The results also indicate that for soliton solutions, the model training costs significantly less time than other initial conditions.     

FREE VIBRATIONS OF MULTILAYERED PLATES BASED ON A MIXED VARIATIONAL APPROACH IN CONJUNCTION WITH GLOBAL PIECEWISE-SMOOTH FUNCTIONS

This paper presents a two-dimensional model for the analysis of freely vibrating laminated plates. The governing differential equations, associated boundary conditions and constitutive equations are derived from Reissner's mixed variational theorem. Both the governing differential equations and the related boundary conditions are presented in terms of resultant stresses and displacements. The model is able to provide the results for the corresponding three-dimensional theory. Such a performance is guaranteed from an appropriate expansion of relevant kinetic and stress quantities through the thickness of the multilayered plate. The expansion is realized by using a novel selection of global piecewise-smooth functions (GPSFs). The number of GPSFs can be arbitrarily increased to achieve a two-dimensional plate theory which is, at least, as accurate as that of a full layerwise theory. It is also shown that GPSFs permit to deal a multilayered plate as if it was virtually made of a single layer. Indeed, the theory need not explicitly introduce continuity conditions for both displacements and relevant stresses. The performance of the present two-dimensional model in conjunction with the global piecewise-smooth functions is tested and discussed by comparing its resulting eigen-parameters, for a class of simply supported plates, with those of other two-dimensional models and with those existing of the exact three-dimensional theory.