Integrating factors and conservation laws for nonconservative dynamical systems

A general approach to the construction of conservation laws for classical nonconservative dynamical systems is presented. The conservation laws are constructed by finding corresponding integrating factors for the equations of motion. Necessary conditions for existence of the conservation laws are studied in detail. A connection between an a priori known conservation law and the corresponding integrating factors is established. The theory is applied to two particular problems.     

Space conservation law in finite volume calculations of fluid flow

In the numerical solutions of fluid flow problems in moving co-ordinates, an additional conservation equation, namely the space conservation law, has to be solved simultaneously with the mass, momentum and energy conservation equations. In this paper a method of incorporating the space conservation law into a finite volume procedure is proposed and applied to a number of test cases. The results show that the method is efficient and produces accurate results for all grid velocities and time steps for which temporal accuracy suffices. It is also demonstrated, by analysis and test calculations, that not satisfying the space conservation law in a numerical solution procedure introduces errors in the form of artificial mass sources. These errors can be made negligible only by choosing a sufficiently small time step, which sometimes may be smaller than required by the temporal discretization accuracy.     

Conservation laws and solutions of a quantum drift-diffusion model for semiconductors

A non-linear system of partial differential equations describing a quantum drift-diffusion model for semiconductor devices is investigated by methods of group analysis. An infinite number of conservation laws associated with symmetries of the model are found. These conservation laws are used for representing the system of equations under consideration in the conservation form. Exact solutions provided by the method of conservation laws are discussed. These solutions are different from invariant solutions.     

Method of <Emphasis Type="Italic">A</Emphasis>-operators and conservation laws for the equations of gas dynamics

An algorithm is proposed which allows all conservation laws for a system of differential equations to be to obtained from its one zero-order conservation law for which the general rank of the Jacobi matrix is equal to the number of independent variables of the system. The efficiency of the algorithm is shown by examples of the equations of gas dynamics, for which new conservation laws are derived. For the equations considered, additional symmetry properties related to these conservation laws are established. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 53–60, March–April, 2009.     

Variational principles and conservation laws to the Burridge–Knopoff equation

We generate conservation laws for the Burridge–Knopoff equation which model nonlinear dynamics of earthquake faults by a new conservation theorem proposed recently by Ibragimov. One can employ this new general theorem for every differential equation (or systems) and derive new local and nonlocal conservation laws. Nonlocal conservation laws comprise nonlocal variables defined by the adjoint equations to the Burridge–Knopoff equation.     

On Time Reversible Description of the Process of Coagulation and Fragmentation

The process of coagulation is associated with scalar conservation laws, where the adhesion particle dynamics results from shock waves. Conversely, the fragmentation of a massive particle into a number of smaller ones, or into a continuous (dust) distribution, is associated with rarefaction waves. It is generally agreed that a reversible solution of a conservation law can include neither shock waves nor the spontaneous emergence of rarefaction waves. The present paper is an attempt to demonstrate that both coagulation and fragmentation may coexist for a reversible solution, under a natural generalization of the system of conservation law. This is done by introducing an action principle which includes, in addition to the inertial (kinetic energy) term, also an appropriately defined internal energy. The above generalization of the system of conservation law appears as the Euler–Lagrange equations for this action.     

Conservation laws in linear elasticity based upon divergence transformations

The existence of conservation laws in linear elasticity based upon divergence transformations of the Lagrangian density function is investigated. It is found that there exist a set of conservation laws which correspond to infinitesimal homogeneous perturbations of the strain and velocity fields. These conservation laws have a unique feature not shared by other conservation laws in linear elasticity in that they contain an arbitrary free parameter.     

一般坐标系下的二维Euler方程时—空守恒格式

将文「１－３」中的一维时－空守恒格式推广到了二维情形,得到了一般坐标系下的二维Ｅｕｌｅｒ方程时－空守恒格式,并用几个典型算例进行了检验计算,结果表明：本文得到的二维时－空守恒格式保留了一维格式所有的优点,格式简单,通用性强,而且对激波等间断具有很高的分辨率。     

Theory of multiphase mixtures—A thermomechanical formulation

The paper presents a theory of mixtures with the nonzero interfacial area between the constituents of the mixture. The conservation laws are physically motivated by utilizing a volume averaging procedure and by the definition of a mapping transformation. It is shown that the theory constructed in this manner is consistent with the theory of mixtures with a vanishingly small interfacial area and that a second law of thermodynamics can be assigned for each phase of the mixture. The conservation laws are examined for invariance properties with the principle of the material frame indifference, and a particular constitutive assumption is discussed. Also presented in the paper are the conservation laws in the integral form and the jump conditions for the singular surfaces in the multiphase mixture.     

二维时-空守恒格式及其在激波传播计算中的应用

将改进后的一维时－空守恒格式推广到了二维情形,得到了一个新的一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式,并用该格式对几个具有复杂波系的流场进行了数值模拟。结果表明,该格式保留了一维格式通用性好、结构简单的优点,其计算结果精度高,对激波等间断具有很强的分辨率。     

A unified theory for turbulent wake flows described by eddy viscosity

In this paper, we consider the conservation laws for the far downstream wake equations described by eddy viscosity. A basis of conserved vectors is constructed. The well-known conserved quantities for the turbulent classical wake and the turbulent wake of a self-propelled body are obtained by integrating the corresponding conservation law across the wake and imposing the boundary conditions. For the wake of a self-propelled body the additional condition that the drag on the body is zero and is required to obtain the conserved quantity. A third conservation law, which possibly belongs to another type of wake, is discovered. The Lie point symmetry associated with the conserved vector is used to obtain the invariant solution and a typical velocity profile for this wake is provided. This wake appears to have common properties with the other two well-known wakes. We then analyse the invariant solutions to all three wake problems and prove that a simple mathematical relationship exists between them thus unifying the theory for turbulent wake flows.     

On the uniqueness of symmetric algebraic consequences implied by a system of conservation laws consistent with the additional conservation equation

In mathematical physics, one often encounters systems of conservation laws which are consistent with an additional conservation equation. Such systems are of particular interest from the point of view of phenomenological thermodynamics where the additional conservation equation is often interpreted as the entropy law. The systems of conservation laws which imply the additional conservation law are strongly related to symmetric systems. These relations are exploited in thermodynamical theories where the system of field equations consistent with the balance of entropy is often assumed to be symmetric.In this paper we use an invariant definition of symmetric system in order to show that the system of balance laws implies the additional balance law if and only if it implies a symmetric system of a certain kind (see Section 2) and that such a symmetric system is uniquely defined.This property is interesting in the context of a more general question; what conditions for a given system of conservation laws are necessary and/or sufficient to ensure the existence of the additional conservation law.     

CONSERVATION LAW AND APPLICATION OF J-INTEGRAL IN MULTI-MATERIALS

ＩｎｔｒｏｄｕｃｔｉｏｎＦｏｒｍａｔｅｒｉａｌｓｃｉｅｎｃｅａｎｄｅｎｇｉｎｅｅｒｉｎｇｓｔｒｕｃｔｕｒｅａｎａｌｙｓｉｓ,ＳｔｒｅｓｓＩｎｔｅｎｓｉｔｙＦａｃｔｏｒ (ＳＩＦ)ａｎｄＥｎｅｒｇｙＲｅｌｅａｓｅＲａｔｉｏ (ＥＲＲ)ａｒｅｔｗｏｉｍｐｏｒｔａｎｔｉｎｄｅｘｅｓ.Ｉｆｍｅｄｉｕｍｉｓｌｉｎｅａｒｅｌａｓｔｉｃ,ｔｈｅＥＲＲｉｓｅｑｕａｌｔｏＪ＿ｉｎｔｅｇｒａｌ.Ａｎｄ ,ｔｈｉｓｅｎｅｒｇｙｉｎｔｅｇｒａｌｉｓｎｏｔｏｎｌｙａｐｐｌｉｃａｂｌｅｔｏｌｉｎｅａｒｅｌａｓｔｉｃｍｅｄｉａ ,ｂｕｔａ…     

Multi-symplectic methods for membrane free vibration equation

In this paper,the multi-symplectic formulations of the membrane free vi- bration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws—a multi-symplectic conservation law(CLS),a local energy conservation law(ECL)as well as a local momentum conservation law(MCL)—is constructed to discrete the PDEs that are derived from the membrane free vibra- tion equation.The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior.     

A 3D staggered Lagrangian scheme for ideal magnetohydrodynamics on unstructured meshes

In this paper, we propose a 3D staggered Lagrangian scheme for the ideal magnetohydrodynamics (MHD) on unstructured meshes. All the thermal variables and the magnetic induction are defined in the cell centers while the fluid velocity is located at the nodes. The meshes are compatibly discretized to ensure the geometric conservation laws in Lagrangian computation by the classical subcell method, then the momentum equation is discretized using the subcell forces and the specific internal energy equation is obtained by the total energy conservation. Invoking the Galilean invariance, magnetic flux conservation, and the thermodynamic consistency, the expressions of subcell force as well as the cell-centered velocity are derived. Besides, the magnetic divergence-free constraint is fulfilled by a projection method after each time step. Various numerical tests are presented to assert the robustness and accuracy of our scheme.     

一般形式的二维时—空守恒格式

本文将经作者改进后的一维时－空守恒格式推广到了二维情形，得到了一个一般形式的二维Ｅｕｌｅｒ方程时－空守恒格式，该格式对各种不规则几何区域内的流动问题具有很强的适应性，同时它还保留了一维格式的优点。几个典型算例的计算结果表明，本文格式不仅精度高，通用性好，而且对激波等间断具有很高的分辨率。     

WKBJ近似保辛吗?

WKBJ短波近似是最常用的有效求解方法之一。保守体系的微分方程可用Hamilton体系的方法描述,其特点是保辛。保辛给出保守体系结构最重要的特性。但WKBJ短波近似却未曾考虑保辛的问题。本文给出验证近似解保辛的条件,并指出WKBJ近似难于保辛。然后给出正则变换的摄动保辛方法。数值例题展示了提出的保辛算法的有效性。     

固体炸药爆轰与惰性介质相互作用的一种扩散界面模型

提出一种保持热力学一致性的扩散界面模型，用来数值模拟固体炸药爆轰与惰性介质的相互作用问题。基于混合网格内各组分物质间可以达到力学平衡状态而不能达到热学平衡状态的假设，由混合网格能量守恒以及压力相等条件，推导出每种组分物质的体积分数演化方程。由此获得的扩散界面模型包括组分物质的质量守恒方程、混合物质的动量及总能量守恒方程，同时包括组分物质的体积分数演化方程和混合物质的压力演化方程。该扩散界面模型的主要特点是考虑了化学反应以及热学非平衡的影响。提出的扩散界面模型在物质界面附近不会出现物理量的非物理振荡现象、适用于任意表达形式的物质状态方程以及任意数目的惰性介质。     

Stress intensity factor calculations based on a conservation integral

It is observed that one of the integral conservation laws of elastostatics, the so-called M-integral conservation law, has certain special features which make it possible to apply this conservation law for a class of plane elastic crack problems in order to calculate the elastic stress intensity factor in each case without solving the corresponding boundary value problem. The main characteristics which a problem must have in order for the approach to be useful are (1) for points very near to the origin of coordinates, the known elastic stresses are 0(r^?r) where r is the radial coordinate and γ ? 1, (2) for points very far from the origin, the known elastic stresses are 0(r^?r) where γ ? 1, and (3) the boundary of the body is made up of radial lines on which certain traction and/or displacement conditions are satisfied. The approach is demonstrated by determining the stress intensity factors for four familiar elastic crack problems directly from the conservation law, and then four similar additional applications of the M-integral conservation law are discussed.