LaNets:Hybrid Lagrange Neural Networks for Solving Partial Differential Equations

We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

A Modified Lagrange Method for Solving Convex Quadratic Optimization Problems

In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

On Form Invariance,Lie Symmetry and Three Kinds of Conserved Quantities of Generalized Lagrange's Equations

In this paper,the form invariance and the Lie symmetry of Lagrange's equations for nonconservativesystem in generalized classical mechanics under the infinitesimal transformations of group are studied,and the Noether'sconserved quantity,the new form conserved quantity,and the Hojman's conserved quantity of system are derived fromthem.Finally,an example is given to illustrate the application of the result.

Lagrange equations of nonholonomic systems with fractional derivatives

This paper obtains Lagrange equations of nonholonomic systems with fractional derivatives. First, the exchanging relationships between the isochronous variation and the fractional derivatives are derived. Secondly, based on these exchanging relationships, the Hamilton＇s principle is presented for non-conservative systems with fractional derivatives. Thirdly, Lagrange equations of the systems are obtained. Furthermore, the d＇Alembert-Lagrange principle with fractional derivatives is presented, and the Lagrange equations of nonholonomic systems with fractional derivatives are studied. An example is designed to illustrate these results.

Discrete Fractional Lagrange Equations of Nonconservative Systems

In order to study discrete nonconservative system,Hamilton's principle within fractional difference operators of Riemann-Liouville type is given. Discrete Lagrange equations of the nonconservative system as well as the nonconservative system with dynamic constraint are established within fractional difference operators of Riemann-Liouville type from the view of time scales. Firstly,time scale calculus and fractional calculus are reviewed.Secondly,with the help of the properties of time scale calculus,discrete Lagrange equation of the nonconservative system within fractional difference operators of Riemann-Liouville type is presented. Thirdly,using the Lagrange multipliers,discrete Lagrange equation of the nonconservative system with dynamic constraint is also established.Then two special cases are discussed. Finally,two examples are devoted to illustrate the results.

Conformal invariance and conserved quantity of third-order Lagrange equations for non-conserved mechanical systems

This paper studies conformal invariance and conserved quantity of third-order Lagrange equations for non- conserved mechanical systems. Third-order Lagrange equations, the definition and a determining equation of conformal invariance of the system are presented. The conformal factor expression is deduced from conformal invariance and Lie symmetry. The necessary and sufficient condition that conformal invaxiance of the system would have Lie symmetry under single-parameter infinitesimal transformations is obtained. The corresponding conserved quantity of conformal invariance is derived with the aid of a structure equation. Lastly, an example is given to illustrate the application of the results.

Noether conserved quantities and Lie point symmetries of difference Lagrange-Maxwell equations and lattices

This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.

能量泛函及Euler-Lagrange方程在图像降噪中的应用研究

研究了自适应分数阶偏微分方程修正模型的能量泛函及Euler-Lagrange方程。首先,定义了自适应分数阶偏微分方程修正模型的能量泛函,其中包含未知函数和拉格朗日乘子的集合。然后,通过求解能量泛函的极值方程,推导出了Euler-Lagrange方程。最后,讨论了Euler-Lagrange方程在自适应分数阶偏微分方程修正模型中的应用。

Lagrange-Noether method for solving second-order differential equations

The purpose of this paper is to provide a new method called the Lagrange-Noether method for solving second-order differential equations. The method is, firstly, to write the second-order differential equations completely or partially in the form of Lagrange equations, and secondly, to obtain the integrals of the equations by using the Noether theory of the Lagrange system. An example is given to illustrate the application of the result.

Comparison of Numerical Approximations of One-Dimensional Space Fractional Diffusion Equation Using Different Types of Collocation Points in Spectral Method Based on Lagrange’s Basis Polynomials

Recently many research works have been conducted and published regarding fractional order differential equations. There are several approaches available for numerical approximations of the solution of fractional order diffusion equations. Spectral collocation method based on Lagrange’s basis polynomials to approximate numerical solutions of one-dimensional (1D) space fractional diffusion equations are introduced in this research paper. The proposed form of approximate solution satisfies non-zero Dirichlet’s boundary conditions on both boundaries. Collocation scheme produce a system of first order Ordinary Differential Equations (ODE) from the fractional diffusion equation. We applied this method with four different sets of collocation points to compare their performance.

CONSTRAINT VIOLATION STABILIZATION OF EULER-LAGRANGE EQUATIONS WITH NON-HOLONOMIC CONSTRAINTS

Two constraint violation stabilization methods are presented to solve the Euler Lagrange equations of motion of a multibody system with nonholonomic constraints. Compared to the previous works, the newly devised methods can deal with more complicated problems such as those with nonholonomic constraints or redundant constraints, and save the computation time. Finally a numerical simulation of a multibody system is conducted by using the methods given in this paper.

LAGRANGE EQUATION OF ANOTHER CLASS OF NONHOLONOMIC SYSTEMS

Using the method of [1], the present paper derives the Lagrange equation without multipliers for another class of first-order nonholonomic dynamical systems by means of variational principle. This kind of equations is also new.

The Sufficient and Necessary Condition of Lagrange Stability of Quasi-periodic Pendulum Type Equations

The quasi-periodic pendulum type equations are considered. A sufficient and necessary condition of Lagrange stability for this kind of equations is obtained. The result obtained answers a problem proposed by Moser under the quasi-periodic case.

The Hamiltonian Canonical Form for Euler-Lagrange Equations

Based on the theory of calculus of variation, some suffcient conditions are given for some Euler-Lagrangcequations to be equivalently represented by finite or even infinite many Hamiltonian canonical equations. Meanwhile,some further applications for equations such as the KdV equation, MKdV equation, the general linear Euler Lagrangeequation and the cylindric shell equations are given.

Numerical Solution of Two Dimensional Fredholm Integral Equations of the Second Kind by the Barycentric Lagrange Function

This paper solves the two dimensional linear Fredholm integral equations of the second kind by combining the meshless barycentric Lagrange interpolation functions and the Gauss-Legendre quadrature formula. Inspired by this thought, we convert the equations into the associated algebraic equations. The results of the numerical examples are given to illustrate that the approximated method is feasible and efficient.