Modulation equations and Reynolds averaging for finite-amplitude non-linear waves in an incompressible fluid

A formal perturbation scheme is developed to determine originalmodulation equations for laminar finite-amplitude non-linearwaves in an incompressible fluid. Three idealized problems areanalysed. The modulation equations comprise conservation ofwaves, averaged conditions for conservation of mass, momentum,kinetic energy and angular momentum and the averaged projectionof the Navier–Stokes equations onto the vorticity vector.The last of these modulation equations, which is related tovortex stretching, only appears in 3D problems. The techniqueof Reynolds averaging is also employed to obtain equations forthe mean velocities and pressure. The Reynolds-averaged Navier–Stokesequations correspond to the modulation equations for conservationof mass and momentum. However, the Reynolds stress transportequations are shown to be inconsistent with the other necessarymodulation equations. In two further idealized problems, exactsolutions of the Navier–Stokes equations are obtainedby employing the modulation equations.     

THE IMPROVED FOURIER SPLITTING METHOD AND DECAY ESTIMATES OF THE GLOBAL SOLUTIONS OF THE CAUCHY PROBLEMS FOR NONLINEAR SYSTEMS OF FLUID DYNAMICS EQUATIONS

Consider the Cauchy problems for an n-dimensional nonlinear system of fluid dynamics equations. The main purpose of this paper is to improve the Fourier splitting method to accomplish the decay estimates with sharp rates of the global weak solutions of the Cauchy problems. We will couple together the elementary uniform energy estimates of the global weak solutions and a well known Gronwall''s inequality to improve the Fourier splitting method. This method was initiated by Maria Schonbek in the 1980''s to study the optimal long time asymptotic behaviours of the global weak solutions of the nonlinear system of fluid dynamics equations. As applications, the decay estimates with sharp rates of the global weak solutions of the Cauchy problems for $n$-dimensional incompressible Navier-Stokes equations, for the $n$-dimensional magnetohydrodynamics equations and for many other very interesting nonlinear evolution equations with dissipations can be established.     

On the T 3-Gowdy Symmetric Einstein-Maxwell Equations

Recently, progress has been made in the analysis of the expanding direction of Gowdy spacetimes. The purpose of the present paper is to point out that some of the techniques used in the analysis can be applied to other problems. The essential equations in the case of the Gowdy spacetimes can be considered as a special case of a wider class of variational problems. Here we are interested in the asymptotic behaviour of solutions to this class of equations. Two particular members arise when considering the T³-Gowdy symmetric Einstein-Maxwell equations and when considering T³-Gowdy symmetric IIB superstring cosmology. The main result concerns the rate of decay of a naturally defined energy. A subclass of the variational problems can be interpreted as wave map equations, and in that case one gets the following picture. The non-linear wave equations one ends up with have as a domain the positive real line in Cartesian product with the circle. For each point in time, the wave map can thus be seen as a loop in some Riemannian manifold. As a consequence of the decay of the energy mentioned above, the length of the loop converges to zero at a specific rate. Communicated by Sergiu Klainerman submitted 14/02/05, accepted 21/04/05     

THE GENERAL SCHEME OF THE GENERALIZED WRONSKIAN TECHNIQUE

The general scheme of the generalized Wronskian technique ms given and is applied tofour energy-dependent potential's eigenvalue problems to generate four classes of nonlinearevolution equations solvable by the inverse spectral transform.     

对带有微结构的弹性固体理论的再研究

对现有的带有微结构的弹性固体理论进行了再研究,并指出由于引进许多记号而使诸基本方程的推导过程复杂化的原因。针对该理论中存在的问题,提出更为普遍的功率能率原理,并借助此原理和广义Piola定理即可自然地推导出局部和非局部理论的运动方程、能率和能量的均衡方程、本构方程和边界条件。这些基本方程和边界条件连同初始条件一起即可用来解决带有微结构的弹性固体动力学理论的混合问题。     

大气环流闭合方程组的总能量守恒格式及全球预报误差的严格估计

本文讨论大气环流闭合方程组,由于同时考虑了热传导效应,内摩擦效应及表达动能向内能转化的耗散项,因此符合总能量守恒律,文中对这一方程组建立了加权平均守恒型差分格式,并证明当选择最优参数时,它满足离散形式的总能量守恒律,通常的二次守恒格式是其次优的情况,文中还综合应用了Jessen不等式,Hardy不等式等等,从而严格证明了在一定条件下,存在t₀>0,当t     

发展方程的计算稳定性问题

一、演变过程方程及差分格式 在数值天气预报中以及求解非定常流体运动时,必须设计计算稳定的格式,所以关于计算稳定性问题的理论研究是很有意义的.在这一类问题中,所要求解的问题大都可以     

Gradient stability of numerical algorithms in local nonequilibrium problems of critical dynamics

The critical dynamics of a spatially inhomogeneous system are analyzed with allowance for local nonequilibrium, which leads to a singular perturbation in the equations due to the appearance of a second time derivative. An extension is derived for the Eyre theorem, which holds for classical critical dynamics described by first-order equations in time and based on the local equilibrium hypothesis. It is shown that gradient-stable numerical algorithms can also be constructed for second-order equations in time by applying the decomposition of the free energy into expansive and contractive parts, which was suggested by Eyre for classical equations. These gradient-stable algorithms yield a monotonically nondecreasing free energy in simulations with an arbitrary time step. It is shown that the gradient stability conditions for the modified and classical equations of critical dynamics coincide in the case of a certain time approximation of the inertial dynamics relations introduced for describing local nonequilibrium. Model problems illustrating the extended Eyre theorem for critical dynamics problems are considered.     

Dynamical problems for geometrically exact theories of nonlinearly viscoelastic rods

Summary This paper surveys recent results and open problems for the equations of motion for geometrically exact theories of nonlinearly viscoelastic and elastic rods. These rods can deform in space by undergoing not only flexure and torsion, but also extension and shear. The paper begins with a derivation of the governing equations, which for viscoelastic rods form a quasilinear system of hyperbolic-parabolic partial differential equations of high order. It then derives the energy equation and discusses difficulties that can arise in getting useful energy estimates. The paper next treats constitutive assumptions precluding total compression. The paper then discusses the curious asymptotic problems that arise when the inertia of the rod is small relative to that of a rigid body attached to its end. The paper concludes with discussions of traveling waves and shock structure, Hopf bifurcation problems, and problems of control. This paper is dedicated to the memory of Juan C. Simo This paper was solicited by the editors to be part of a volume dedicated to the memory of Juan Simo.     

A priori bounds for degenerate and singular evolutionary partial integro-differential equations

We study quasilinear evolutionary partial integro-differential equations of second order which include time fractional p-Laplace equations of time order less than one. By means of suitable energy estimates and De Giorgi’s iteration technique we establish results asserting the global boundedness of appropriately defined weak solutions of these problems. We also show that a maximum principle is valid for such equations.     

Explicit-implicit difference scheme for the joint solution of the radiative transfer and energy equations by the splitting method

High-order accurate explicit and implicit conservative predictor-corrector schemes are presented for the radiative transfer and energy equations in the multigroup kinetic approximation solved together by applying the splitting method with respect to physical processes and spatial variables. The original system of integrodifferential equations is split into two subsystems: one of partial differential equations without sources and one of ordinary differential equations (ODE) with sources. The general solution of the ODE system and the energy equation is written in quadratures based on total energy conservation in a cell. A feature of the schemes is that a new approximation is used for the numerical fluxes through the cell interfaces. The fluxes are found along characteristics with the interaction between radiation and matter taken into account. For smooth solutions, the schemes approximating the transfer equations on spatially uniform grids are second-order accurate in time and space. As an example, numerical results for Fleck’s test problems are presented that confirm the increased accuracy and efficiency of the method.     

Refined Orthotropic Plate Motion Equations for Acoustoelasticity Problem Statement

The formulation of the acoustoelasticity problem is given on the basis of refined motion equations of orthotropic plates. These equations are constructed in the first approximation by reducing the three-dimensional equations of the theory of elasticity to the two-dimensional equations of the theory of plates, where the approximation of the transverse tangential stresses and the transverse reduction stress is made with the help of trigonometric basis functions in the thickness direction. Wherein at the points of the boundary (front) surfaces, the static boundary conditions of the problem for tangential stresses are satisfied exactly and for transverse normal stress — approximately. Accounting for internal energy dissipation in the plate material is based on the Thompson—Kelvin—Voigt hysteresis model. In case of formulating problems on dynamic processes of plate deformation in vacuum, the equations are divided into two separate systems of equations. The first of these systems describes non-classical shear-free, longitudinal-transverse forms of movement, accompanied by a distortion of the flat form of cross sections, and the second system describes transverse bending-shear forms of movement. The latter are practically equivalent in quality and content to the analogous equations of the well-known variants of refined theories, but, unlike them, with a decrease in the relative thickness parameter, they lead to solutions according to the classical theory of plates. The motion of the surrounding the plate acoustic media is described by the generalized Helmholtz wave equations, constructed with account of energy dissipation by introducing into consideration the complex sound velocity according to Skudrzyk.     

Some non-local problems for the parabolic-hyperbolic type equation with non-characteristic line of changing type

In this work two non-local problems for the parabolic-hyperbolic type equation with non-characteristic line of changing type are considered. Unique solvability of these problems is proven. The uniqueness of the solution is proven by the method of energy integrals and the existence is proven by the method of integral equations.     

Some Regularity Results In Ferromagnetism

The magnetism of a rigid ferromagnet occupying a spatial region Ωis described by a unit vectorfield m on Ω. The total energy of m involves several terms: the anisotropy energy imposed by the lattice structure of the material, the exchange energy discouraging very rapid local changes in m , the applied energy due to external magnetic sources, and the induced magnetic field energy. Here, while incorporating all energy terms, we show that minimizers have at most isolated singularities, usually in the interior of Ω, and that there is nice asymptotic behavior at such singularities. In contrast to related harmonic map problems, the field energy is a nonlocal term, involving a solution of Maxwell's equations with coefficients depending on m     

Some non‐local problems for the parabolic‐hyperbolic type equation with complex spectral parameter

In the present paper the unique solvability of two non‐local problems for the mixed parabolic‐hyperbolic type equation with complex spectral parameter is proved. Sectors for values of the spectral parameter where these problems have unique solutions are shown. Uniqueness of the solution is proved by the method of energy integral and existence is proved by the method of integral equations. In particular cases, eigenvalues and corresponding eigenfunctions of the studied problems are found. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimization problems on general classes of rearrangements

This paper is concerned with maximization and minimization problems of the energy integral associated to p-Laplace equations depending on functions that belong to a class of rearrangements. We prove existence and uniqueness results, and present some features of optimal solutions. The radial case is discussed in detail. We also prove a result of uniqueness for a class of p-Laplace equations under non-standard assumptions.     

Energetic BEM-FEM coupling for the numerical solution of the damped wave equation

Time-dependent problems modeled by hyperbolic partial differential equations can be reformulated in terms of boundary integral equations and solved via the boundary element method. In this context, the analysis of damping phenomena that occur in many physics and engineering problems is a novelty. Starting from a recently developed energetic space-time weak formulation for the coupling of boundary integral equations and hyperbolic partial differential equations related to wave propagation problems, we consider here an extension for the damped wave equation in layered media. A coupling algorithm is presented, which allows a flexible use of finite element method and boundary element method as local discretization techniques. Stability and convergence, proved by energy arguments, are crucial in guaranteeing accurate solutions for simulations on large time intervals. Several numerical benchmarks, whose numerical results confirm theoretical ones, are illustrated and discussed.     

Difference spline scheme for modeling of convective flows using full Navier — Stokes equations

A difference scheme is constructed, in which enhanced stability is achieved by simultaneous solutions of the equations of motion, energy, and continuity. Spline approximations of spatial derivatives (with the original equations written in divergence form) substantially improve the accuracy of the scheme compared with the standard difference scheme using symmetric differences. The efficiency of the scheme is demonstrated for some problems of convective flow of compressible gas with lateral and bottom heating.Translated from Matematicheskoe Modelirovanie i Reshenie Obratnykh Zadach. Matematicheskoi Fiziki, pp. 38–45, 1993.     

Semi-implicit method for thermodynamically linked equations in phase change problems (SIMTLE)

Thermodynamic coupling of temperature and composition fields in phase-change problems has been a challenge for decades. A compromise has been always desired between numerical efficiency and detailed physical consideration, toward a general scheme. In the present work, a macro–micro numerical method is proposed to link the conservation equations of energy and species with the thermodynamics of the solidification problems. Firstly, the basic structure of the method, simplified with a local equilibrium assumption, is presented. The method is then extended to a multi-phase model, demonstrating a three-phase approach to the solidification of a eutectic binary alloy. Relaxing the limitations imposed by the equilibrium assumption, non-equilibrium and microscale considerations was also included subsequently by a suggested modification to the macroscopic mathematical model. Advantages gained through the general algorithm proposed are concerned with two features of the method; (a) consistency with the energy and species equations. (b) No need of a predefined solidification path; that allows for the usage of raw phase diagram curves and offers simplicity and generality for extension through complex problems (i.e. microscopic, multi-phase or non-equilibrium). A benchmark problem was employed to test the performance of the proposed method in two cases of local equilibrium and Scheil-like solidification. The obtained results were validated in comparison with available semi-analytical solution.