基于Lagrange方法封闭的两相湍流场方程模型

两相湍流场方程模型采用基于Euler方法的一阶矩方程,而二阶矩方程由Lagrange方法得到,新模型的封闭不需要附加其它假设.首先基于概率密度函数给出颗粒运动的连续方程和动量方程,其次由基于平均Langevin方程的Lagrange模型推导得到颗粒二阶矩方程,最终获得封闭的二阶矩模型.将新模型用于气固两相壁面射流的数值模拟,结果表明新模型合理有效.     

关于Birkhoff表示的Lagrange像的研究

研究运动微分方程Birkhoff表示的Lagrange像.得出二阶Lagrange函数应满足的条件,在此条件下广义Lagrange方程为二阶微分方程组;提出新的求解Lagrange力学逆问题路线;指出在此问题研究中曾发生过的失误.举例说明所得结果的应用.     

Variational iteration method for solving time-fractional diffusion equations in porous the medium

<正>The variational iteration method is successfully extended to the case of solving fractional differential equations, and the Lagrange multiplier of the method is identified in a more accurate way.Some diffusion models with fractional derivatives are investigated analytically,and the results show the efficiency of the new Lagrange multiplier for fractional differential equations of arbitrary order.     

完整力学系统的高阶运动微分方程

从质点系的牛顿动力学方程出发,引入系统的高阶速度能量,导出完整力学系统的高阶Lagrange方程、高阶Nielsen方程以及高阶Appell方程,并证明了完整系统三种形式的高阶运动微分方程是等价的.结果表明,完整系统高阶运动微分方程揭示了系统运动状态的改变与力的各阶变化率之间的联系,这是牛顿动力学方程以及传统分析力学方程不能直接反映的.因此,完整系统高阶运动微分方程是对牛顿动力学方程及传统Lagrange方程、Nielsen方程、Appell方程等二阶运动微分方程的进一步补充. 关键词： 高阶速度能量 高阶Lagrange方程 高阶 Nielsen方程 高阶Appell方程     

Integration method for the Dirac equation

An integration method for the Dirac equation is proposed. The method, based on diagonalization, reduces the problem to one of integration of independent second-order differential equations.     

LEPAGE EQUIVALENTS OF SECOND-ORDER EULER–LAGRANGE FORMS AND THE INVERSE PROBLEM OF THE CALCULUS OF VARIATIONS

In the calculus of variations, Lepage (n + 1)-forms are closed differential forms, representing Euler–Lagrange equations. They are fundamental for investigation of variational equations by means of exterior differential systems methods, with important applications in Hamilton and Hamilton–Jacobi theory and theory of integration of variational equations. In this paper, Lepage equivalents of second-order Euler–Lagrange quasi-linear PDE's are characterised explicitly. A closed (n + 1)-form uniquely determined by the Euler–Lagrange form is constructed, and used to find a geometric solution of the inverse problem of the calculus of variations.     

Global Lagrange stability for neutral-type inertial neural networks with discrete and distributed time delays

This paper focus on the problem of global Lagrange stability for neutral-type inertial neural networks with discrete and distributed time delays. By choosing a proper variable substitution, an inertial neural network consisting of second-order differential equations can be converted into a first-order differential model. The sufficient conditions of the inertial neural network with neutral delay are derived by constructing suitable Lyapunov-Krasovskii functional candidates, introducing new free weighting matrices, utilizing inequality techniques and analytical method. Through the LMI condition, we analyze the global exponential stability of the delayed inertial neural networks in Lagrange sense. Meanwhile, the global exponential attractive set is also given. Finally, some example is given to illustrate our theoretical results.     

First integrals and stability of second-order differential equations

The stability of second-order differential equations is studied by using their integrals. A system of second-order differential equations can be considered as a mechanical system with holonomic constraints. A conserved quantity of the mechanical system or a part of the system is obtained by using the Noether theory. It is possible that the conserved quantity becomes a Liapunov function and the stability of the system is proved by the Liapunov theorem.     

Fractional differential equations of motion in terms of combined Riemann—Liouville derivatives

In this paper,we focus on studying the fractional variational principle and the differential equations of motion for a fractional mechanical system.A combined Riemann-Liouville fractional derivative operator is defined,and a fractional Hamilton principle under this definition is established.The fractional Lagrange equations and the fractional Hamilton canonical equations are derived from the fractional Hamilton principle.A number of special cases are given,showing the universality of our conclusions.At the end of the paper,an example is given to illustrate the application of the results.     

非惯性系下考虑剪切变形的柔性梁的动力学建模

研究非惯性坐标系下考虑剪切变形的柔性梁的动力学建模. 首先借鉴Euler-Bernoulli梁的几何非线性变形模式,考虑了Timoshenko梁弯曲以及剪切变形产生的几何非线性效应对纵向、横向变形位移的影响,在考虑两个方向的变形耦合项后,利用有限元法对柔性梁进行了离散,采用Lagrange方程建立了柔性梁的动力学模型,首次建立了包含变形二次耦合量的Timoshenko梁的动力学方程. 关键词： 非惯性坐标系 剪切变形 柔性梁 动力学建模     

Majorana Transformation for Differential Equations

We present a method for reducing the order of ordinary differential equations satisfying a given scaling relation (Majorana scale-invariant equations). We also develop a variant of this method, aimed to reduce the degree of nonlinearity of the lower order equation. Some applications of these methods are carried out and, in particular, we show that second-order Emden–Fowler equations can be transformed into first-order Abel equations. The work presented here is a generalization of a method used by Majorana in order to solve the Thomas–Fermi equation.     

精确Cosserat弹性杆动力学的分析力学方法

以脱氧核糖核酸和工程中的细长结构为背景, 大变形大范围运动的弹性杆动力学受到关注. 将分析力学方法运用到精确Cosserat弹性杆动力学, 旨在为前者拓展新的应用领域, 为后者提供新的研究方法. 基于平面截面假定, 在弯扭基础上再计及拉压和剪切变形形成精确Cosserat弹性杆模型. 用刚体运动的概念描述弹性杆的变形, 导出弹性杆变形和运动的几何关系; 在定义截面虚位移及其变分法则的基础上, 建立用矢量表达的d’Alembert-Lagrange原理, 在线性本构关系下化作分析力学形式, 并导出Lagrange方程和Nielsen方程, 定义正则变量后化作Hamilton正则方程; 对于只在端部受力的弹性杆静力学, 导出了将守恒量预先嵌入的Lagrange方程, 并讨论了其首次积分. 从弹性杆的d’Alembert-Lagrange原理导出积分变分原理, 在线性本构关系下化作Hamilton原理. 形成的分析力学方法使弹性杆的全部动力学方程具有统一的形式, 为弹性杆动力学的对称性和守恒量的研究及其数值计算铺平道路. 关键词： 精确Cosserat弹性杆 分析动力学方法 变分原理 Lagrange方程     

弹性介质的Lagrange动力学与地震波方程

本文将地球介质看作是弹性介质, 从弹性体的Navier方程出发, 建立均匀弹性介质和非均匀弹性介质的分析动力学方程——Lagrange方程, 利用弹性介质的Lagrange方程导出匀弹性介质和非均匀弹性介质的地震波方程, 为用Lagrange分析动力学研究地球介质中地震波传播规律和解决地震勘探中的有关问题提供基础. 关键词： 地震勘探 弹性介质 Lagrange 方程 地震波方程     

基于广义交替数值通量的局部间断Galerkin方法求解二维波动方程

波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播.     

AN OLD METHOD OF JACOBI TO FIND LAGRANGIANS

In a recent paper by Ibragimov a method was presented in order to find Lagrangians of certain second-order ordinary differential equations admitting a two-dimensional Lie symmetry algebra. We present a method devised by Jacobi which enables one to derive (many) Lagrangians of any second-order differential equation. The method is based on the search of the Jacobi Last Multipliers for the equations. We exemplify the simplicity and elegance of Jacobi's method by applying it to the same two equations as Ibragimov did. We show that the Lagrangians obtained by Ibragimov are particular cases of some of the many Lagrangians that can be obtained by Jacobi's method.     

A spectral fictitious domain method with internal forcing for solving elliptic PDEs

A fictitious domain method is presented for solving elliptic partial differential equations using Galerkin spectral approximation. The fictitious domain approach consists in immersing the original domain into a larger and geometrically simpler one in order to avoid the use of boundary fitted or unstructured meshes. In the present study, boundary constraints are enforced using Lagrange multipliers and the novel aspect is that the Lagrange multipliers are associated with smooth forcing functions, compactly supported inside the fictitious domain. This allows the accuracy of the spectral method to be preserved, unlike the classical discrete Lagrange multipliers method, in which the forcing is defined on the boundaries. In order to have a robust and efficient method, equations for the Lagrange multipliers are solved directly with an influence matrix technique. Using a Fourier–Chebyshev approximation, the high-order accuracy of the method is demonstrated on one- and two-dimensional elliptic problems of second- and fourth-order. The principle of the method is general and can be applied to solve elliptic problems using any high order variational approximation.     

基于二阶离散变分方法的非线性映射参数估计

提出一种估计非线性映射未知参数的二阶离散变分方法.首先针对非线性离散混沌系统, 利用变分方法导出了伴随方程和目标泛函梯度, 以此为基础利用二阶离散变分方法给出了二阶伴随方程和精确计算Hessian矩阵-向量乘积的显式表达式; 其次设计了估计非线性映射未知参数的新算法, 并以此对Hyperhenón映射和二维抛物映射中的未知参数进行了精确的估计. 数值仿真结果表明了该方法的有效性和优点. 关键词： 非线性映射 参数估计 二阶离散变分方法 伴随方程     

An alternative approach for the analysis of nonlinear vibrations of pipes conveying fluid

Based on a novel extended version of the Lagrange equations for systems containing non-material volumes, the nonlinear equations of motion for cantilever pipe systems conveying fluid are deduced. An alternative to existing methods utilizing Newtonian balance equations or Hamilton's principle is thus provided. The application of the extended Lagrange equations in combination with a Ritz method directly results in a set of nonlinear ordinary differential equations of motion, as opposed to the methods of derivation previously published, which result in partial differential equations. The pipe is modeled as a Euler elastica, where large deflections are considered without order-of-magnitude assumptions. For the equations of motion, a dimensional reduction with arbitrary order of approximation is introduced afterwards and compared with existing lower-order formulations from the literature. The effects of nonlinearities in the equations of motion are studied numerically. The numerical solutions of the extended Lagrange equations of the cantilever pipe system are compared with a second approach based on discrete masses and modeled in the framework of the multibody software HOTINT/MBS. Instability phenomena for an increasing number of discrete masses are presented and convergence towards the solution for pipes conveying fluid is shown.     

Noether＇s and Poisson＇s methods for solving differential equation xs^（m） = Fs（t,xk^（m-2） , xk（m-1）

This paper studies integration of a higher-order differential equation which can be reduced to a second-order ordinary differential equation. The solution of the second-order equation can be obtained by the Noether method and the Poisson method. Then the solution of the higher-order equation can be obtained by integrating the solution of the second-order equation.     

Perturbation theory in a framework of iteration methods

In a previous paper [J. Phys. A 36 (2003) 11807], we introduced the ‘asymptotic iteration method’ for solving second-order homogeneous linear differential equations. In this Letter, we study perturbed problems in quantum mechanics and we use the method to find the coefficients in the perturbation series for the eigenvalues and eigenfunctions directly, without first solving the unperturbed problem.