On the Semi-geostrophic System in Physical Space with General Initial Data

Archive for Rational Mechanics and Analysis - In order to accommodate general initial data, an appropriately relaxed notion of renormalized Lagrangian solutions for the Semi-Geostrophic system in...     

A Generalized Theory of Stress and Strain Measures in the Classical Continuum Mechanics

A generalized theory of stress and strain tensor measures in the classical continuum mechanics is discussed: the main axioms of the theory are proposed, the general formulas for new tensor measures are derived, arid an energy conjugate theorem is formulated to distinguish the complete Lagrangian class of measures. As a subclass, a simple Lagrangian class of energy conjugate measures of stresses and finite strains is constructed in which the families of holonomic and corotational measures are distinguished. The characteristics of holonomic and corotational measures are studied by comparing the tensor measures of the simple Lagrangian class with one another and with logarithmic measures. For the simple Lagrangian class and its families, their completeness and closure are shown with respect to the choice of a generating pair of energetically conjugate measures. The applications of the new tensor measures in modeling the properties of plasticity, viscoelasticity, and shape memory are mentioned.     

The suspended elastic cable under the action of concentrated vertical loads

This investigation treats the static response of a single elastic cable which is suspended between two points that are not necessarily at the same level. The cable is loaded by its self-weight and any number of concentrated vertical loads which may be arbitrarily placed along its length. The analysis presented uses a Lagrangian approach. For the strained cable profile, the tension and displacements are given as functions of a single Lagrangian co-ordinate. A specific application of the general analysis is made and compared with a simple experiment.     

拉格朗日松弛算法在平台罗经多故障诊断技术中的应用研究

最优多故障诊断问题是一个NP-hard问题。针对平台罗经这一复杂系统，采用有向图描述元件与测试点间的因果依赖关系，并建立系统的多信号模型。在考虑元件发生故障的先验概率的前提下，提出一种基于拉格朗日算法（LRA）和子梯度优化算法（SOA）近最优多故障诊断算法，并在某型平台罗经的方位稳定系统多故障诊断中得到应用。结果表明：该方法相对于传统的单故障诊断方法诊断速度快，适用于平台罗经这类大型复杂系统的多故障诊断。     

A lagrangian density for the dynamics of elastic dielectrics

The equations governing the behaviour of an elastic dielectric in the quasi-static approximation are analysed in connection with the inverse problem of the calculus of variations. First, through compatibility with the potential conditions, the most general constitutive equations are derived which allow the model to admit a variational formulation. Then the explicit form of the Lagrangian density is determined. Consistent with the adoption of Lagrangian coordinates, the Lagrangian density ascribes to the dielectric the structure of a particle-like system.     

The effects of control‐volume distributed multi‐point flux approximation (CVD‐MPFA) on upscaling—A study using the CVD‐MPFA schemes

A family of flux‐continuous, locally conservative, finite‐volume schemes has been developed for solving the general geometry‐permeability tensor (petroleum reservoir‐simulation) pressure equation on structured and unstructured grids and are control‐volume distributed (textit Comput. Geo. 1998; 2 :259–290; Comput. Geo. 2002; 6 :433–452). The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir‐simulation schemes (two‐point flux approximation) when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization (Int. J. Numer. Meth. Fluids 2006; 51 :1177–1203). Improved convergence (for two‐ and three‐dimensional formulation) using the quadrature parameterization has been observed for the family of flux‐continuous control‐volume distributed multi‐point flux approximation (CVD‐MPFA) schemes (Ph.D. Thesis, University of Wales, Swansea, U.K., 2007). In this paper family of flux‐continuous (CVD‐MPFA) schemes are used as a part of numerical upscaling procedure for upscaling the fine‐scale grid information (permeability) onto a coarse grid scale. A series of data‐sets (SPE, 2001) are tested where the upscaled permeability tensor is computed on a sequence of grid levels using the same fixed range of quadrature points in each case. The refinement studies presented involve:

• (i) Refinement comparison study: In this study, permeability distribution for cells at each grid level is obtained by upscaling directly from the fine‐scale permeability field as in standard simulation practice.
• (ii) Refinement study with renormalized permeability: In this refinement comparison, the local permeability is upscaled to the next grid level hierarchically, so that permeability values are renormalized to each coarser level. Hence, showing only the effect of increased grid resolution on upscaled permeability, compared with that obtained directly from the fine‐scale solution.
• (iii) Refinement study with invariant permeability distribution: In this study, a classical mathematical convergence test is performed. The same coarse‐scale underlying permeability map is preserved on all grid levels including the fine‐scale reference solution.

The study is carried out for the discretization of the scheme in physical space. The benefit of using specific quadrature points is demonstrated for upscaling in this study and superconvergence is observed. Copyright © 2010 John Wiley & Sons, Ltd.     

Certain Aspects of an “Incompressible Fluid“ Model in Gas Dynamics

An “incompressible fluid” model in gas dynamics is developed in the linear approximation. Using the dissipative relaxation time as a characteristic scale, we arrive at another form of the dimensionless Boltzmann equation. In the limiting case of small Knudsen numbers an approximate solution is obtained in the form of a Hilbert multiple-scale asymptotic expansion. It is revealed that for slow, weakly nonisothermal processes the asymptotic expansion for the linearized Boltzmann equation leads in a first stage to equations for the velocity, pressure and temperature that do not contain the density (quasi-incompressible approximation). The density depends on the temperature and can, if necessary, be found from the equation of state. The next-approximation equations contain the Burnett effects, the velocity calculation being reduced to the general problem of finding a vector field from a given divergence and rotation. With reference to a simple case of the heating of a stationary gas in a half-space it is shown that the temperature establishment process is accompanied by gas flow from the wall.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, 2005, pp. 170–178.Original Russian Text Copyright © 2005 by Chekmarev.     

Forms of null Lagrangians in field theories of continuum mechanics

The divergence representation of a null Lagrangian that is regular in a star-shaped domain is used to obtain its general expression containing field gradients of order ≤ 1 in the case of spacetime of arbitrary dimension. It is shown that for a static three-component field in the three-dimensional space, a null Lagrangian can contain up to 15 independent elements in total. The general form of a null Lagrangian in the four-dimensional Minkowski spacetime is obtained (the number of physical field variables is assumed arbitrary). A complete theory of the null Lagrangian for the n-dimensional spacetime manifold (including the four-dimensional Minkowski spacetime as a special case) is given. Null Lagrangians are then used as a basis for solving an important variational problem of an integrating factor. This problem involves searching for factors that depend on the spacetime variables, field variables, and their gradients and, for a given system of partial differential equations, ensure the equality between the scalar product of a vector multiplier by the system vector and some divergence expression for arbitrary field variables and, hence, allow one to formulate a divergence conservation law on solutions to the system.     

Higher-Order Terms in the Capillary Limit Approximation of Viscous-Capillary Flow

We propose a new approximation to the solution of the steady state equations describing two-phase immiscible flow in a porous medium. It is demonstrated that the general procedure contains the capillary equilibrium approximation as a special case. The solution is approximated by a perturbation series in a parameter related to the capillary number. The expansion of the solution results in a sequence of decoupled linear elliptic boundary value problems. This sequence is solved numerically by a Finite Element method, and the accuracy of the approximations is evaluated.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

Augmented Lagrangian technique applied to a spatially periodic,harmonic Navier-Stokes problem

An application and an extension (to complex variables) of the classical augmented Lagrangian method is performed. Finite element computations are realized in the two-dimensional case of an harmonic Navier-Stokes problem with periodic boundary conditions. A formulation (extended from the traditional Stokes problem) involving a simple Lagrangian, solved by the Uzawa algorithim, was previously used.¹ This treatment proved unsatisfactory for large frequencies. The efficient and well-known augmented Lagrangian technique solved by the Uzawa algorithm is used to overcome these shortcomings. Other, better techniques could be used. Nevertheless the simple method used here is efficient. Moreover the numerical implementation needs little memory storage, which is an important factor in this particular case. The well-known conditioning technique employed is shown to be well-adapted in this case, a fact which emerges from the study of the non-symmetric problem involved. Finally, many tests, computations and experimental data are presented.     

Finite element analysis of incompressible laminar boundary layer flows

A numerical procedure was developed to solve the two-dimensional and axisymmetric incompressible laminar boundary layer equations using the semi-discrete Galerkin finite element method. Linear Lagrangian, quadratic Lagrangian, and cubic Hermite interpolating polynomials were used for the finite element discretization; the first-order, the second-order backward difference approximation, and the Crank-Nicolson method were used for the system of non-linear ordinary differential equations; the Picard iteration and the Newton-Raphson technique were used to solve the resulting non-linear algebraic system of equations. Conservation of mass is treated as a constraint condition in the procedure; hence, it is integrated numerically along the solution line while marching along the time-like co-ordinate. Among the numerical schemes tested, the Picard iteration technique used with the quadratic Lagrangian polynomials and the second-order backward difference approximation case turned out to be the most efficient to achieve the same accuracy. The advantages of the method developed lie in its coarse grid accuracy, global computational efficiency, and wide applicability to most situations that may arise in incompressible laminar boundary layer flows.     

Fourth-order quasi-harmonic equations of state for solids of cubic and tetragonal symmetry

General expressions are given for the dependence of the pressure and the effective elastic moduli on deformation and temperature in the form of a Taylor series expansion with respect to elastic and thermal strains. The temperature dependence of these expressions is derived within the quasi-harmonic approximation of lattice dynamics. The expressions are developed in terms of the Lagrangian strain and an alternative strain measure identical with the Eulerian strain for a pure deformation. They are then used to obtain the third- and fourth-order equations of state for crystals of cubic and tetragonal symmetry and to relate the parameters entering these equations to quantities which are commonly (or may be potentially) measured experimentally. It is shown that available ultrasonic data are not completely sufficient to evaluate the parameters of fourth-order equations of state. For tetragonal symmetry, this problem is still in abeyance; while in the cubic case, it is possible to estimate the fourth-order parameters from shock-wave data and so to give illustrative numerical applications of our equations. Finally, the third- and fourth-order Hugoniots and isotherms of Cu and Ag are calculated in terms of both the Lagrangian and Eulerian strain measures.     

Développement en polynômes de chaos d'un modèle lagrangien d'écoulement autour d'un profilPolynomial chaos expansion of a Lagrangian model for the flow around an airfoil

A polynomial chaos (PC) expansion a the Lagrangian model for the stochastic incompressible inviscid flow around an airfoil is presented. The flow field is modeled using a distribution of lumped vortices on the airfoil surface while the wake is modeled with Lagrangian point vortices. An original technique is proposed for the computation of the PC coefficients of the velocities induced by the vortices. Two computational examples for random airfoil motions are provided to illustrate the capability of the method to deal with complex situations. To cite this article: O. Le Maître, C. R. Mecanique 334 (2006).     

Uniformly valid approximations for non-linear oscillations with small damping

A uniformly valid zeroth-order approximation is obtained for the general equation y + εH(y)y + M(y)y = 0, where ε is a small parameter. The notion of multiple scaling is utilized to set up a systematic approximation scheme. Examples are given for simple polynomials for H(y) and M(y), which lead to results involving elliptic integrals. Further restrictions allow progress to be made in terms of gamma functions.     

Incremental elastic motions superimposed on a finite deformation in the presence of an electromagnetic field

The equations describing the interaction of an electromagnetic sensitive elastic solid with electric and magnetic fields under finite deformations are summarized, both for time-independent deformations and, in the non-relativistic approximation, time-dependent motions. The equations are given in both Eulerian and Lagrangian form, and the latter are then used to derive the equations governing incremental motions and electromagnetic fields superimposed on a configuration with a known static finite deformation and time-independent electromagnetic field. As a first application the equations are specialized to the quasimagnetostatic approximation and in this context the general equations governing time-harmonic plane-wave disturbances of an initial static configuration are derived. For a prototype model of an incompressible isotropic magnetoelastic solid a specific formula for the acoustic shear wave speed is obtained, which allows results for different relative orientations of the underlying magnetic field and the direction of wave propagation to be compared. The general equations are then used to examine two-dimensional motions, and further expressions for the wave speed are obtained for a general incompressible isotropic magnetoelastic solid.     

Pseudolocalized three-dimensional solitary waves as quasi-particles

A higher-order dispersive equation is introduced as a candidate for the governing equation of a field theory. A new class of solutions of the three-dimensional field equation are considered, which are not localized functions in the sense of the integrability of the square of the profile over an infinite domain. For this new class of solutions, the gradient and/or the Hessian/Laplacian are square integrable. In the linear limiting case, an analytical expression for the pseudolocalized solution is found and the method of variational approximation is applied to find the dynamics of the centers of the quasi-particles (QPs) corresponding to these solutions. A discrete Lagrangian can be derived due to the localization of the gradient and the Laplacian of the profile. The equations of motion of the QPs are derived from the discrete Lagrangian. The pseudomass (“wave mass”) of a QP is defined as well as the potential of interaction. The most important trait of the new QPs is that, at large distances, the force of attraction is proportional to the inverse square of the distance between the QPs. This can be considered analogous to the gravitational force in classical mechanics.