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 metho...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.展开更多
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. Meanwh...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.展开更多
This paper presents, an efficient approach for solving Euler-Lagrange Equation which arises from calculus of variations. Homotopy analysis method to find an approximate solution of variational problems is proposed. An...This paper presents, an efficient approach for solving Euler-Lagrange Equation which arises from calculus of variations. Homotopy analysis method to find an approximate solution of variational problems is proposed. An optimal value of the convergence control parameter is given through the square residual error. By minimizing the the square residual error, the optimal convergence-control parameters can be obtained. It is showed that the homotopy analysis method was valid and feasible to the study of variational problems.展开更多
In this paper the fractional Euler-Lagrange equation is considered.The fractional equation with the left and right Caputo derivatives of order a∈(0,1]is transformed into its corresponding integral form.Next,we presen...In this paper the fractional Euler-Lagrange equation is considered.The fractional equation with the left and right Caputo derivatives of order a∈(0,1]is transformed into its corresponding integral form.Next,we present a numerical solution of the integral form of the considered equation.On the basis of numerical results,the convergence of the proposed method is determined.Examples of numerical solutions of this equation are shown in the final part of this paper.展开更多
Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based...Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.展开更多
By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by si...By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by simulating Ai-and Bi-soliton solutions of the cylindrical Korteweg-de Vries(cKdV)equation.展开更多
In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied tho...In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.展开更多
In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitatio...In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.展开更多
In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation e...In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.展开更多
On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear ...On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].展开更多
This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and M...This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and Marchand(Control Cybernet.2011,40:971-988),we establish an observability inequality of the conservative problem,and then discuss the equivalence between the exponential stabilization of a dissipative system and the internal observational inequality of the corresponding conservative system.展开更多
Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained thro...Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained through the Hirota method and symbolic computation. Breather-like and bound-state solitons are constructed in which the signs of the imaginary parts of the complex wave numbers and the initial separations of the two parallel solitons are important factors for the interaction patterns. The periodic structures and position-induced phase shift of some solutions are introduced.展开更多
The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane ...The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane significantly influences the value of the geodetic coordinates.This research analyzes this influence and how it can contribute to their solutions.The study investigates the mathematical relation between them and presents an exact formula relating to the geodetic height and the ellipse equation.In addition,a heuristic formula for the relation between the geodetic height and the ellipse equation is proposed,which is independent of the geodetic latitude and has a relative accuracy better than 99.9 %.The calculation is stable,and the cost is low.展开更多
The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macrosc...The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.展开更多
Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of ...Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.展开更多
The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid...The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.展开更多
We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the k...We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the kinetic TCS equation for a particle ensemble and the Stokes equations for a fluid via a drag force.In this paper,we present a complete analysis of the existence of global-in-time strong solutions to the coupled model without any smallness restrictions on the initial data.展开更多
A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploratio...A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploration of the time fractional Schrodinger equation within the context of a non-Markovian environment.By leveraging a two-level atom as an illustrative case,we find that the choice to raise i to the order of the time derivative is inappropriate.In contrast to the conventional approach used to depict the dynamic evolution of quantum states in a non-Markovian environment,the time fractional Schrodinger equation,when devoid of fractional-order operations on the imaginary unit i,emerges as a more intuitively comprehensible framework in physics and offers greater simplicity in computational aspects.Meanwhile,we also prove that it is meaningless to study the memory of time fractional Schrodinger equation with time derivative 1<α≤2.It should be noted that we have not yet constructed an open system that can be fully described by the time fractional Schrodinger equation.This will be the focus of future research.Our study might provide a new perspective on the role of time fractional Schrodinger equation.展开更多
基金Project supported by the National Natural Science Foundation of China (No. 19902006).
文摘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.
文摘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.
文摘This paper presents, an efficient approach for solving Euler-Lagrange Equation which arises from calculus of variations. Homotopy analysis method to find an approximate solution of variational problems is proposed. An optimal value of the convergence control parameter is given through the square residual error. By minimizing the the square residual error, the optimal convergence-control parameters can be obtained. It is showed that the homotopy analysis method was valid and feasible to the study of variational problems.
文摘In this paper the fractional Euler-Lagrange equation is considered.The fractional equation with the left and right Caputo derivatives of order a∈(0,1]is transformed into its corresponding integral form.Next,we present a numerical solution of the integral form of the considered equation.On the basis of numerical results,the convergence of the proposed method is determined.Examples of numerical solutions of this equation are shown in the final part of this paper.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11931017 and 12071447)。
文摘Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.
基金supported by the National Natural Science Foundation of China(Grant Nos.12175111 and 12235007)the K.C.Wong Magna Fund in Ningbo University。
文摘By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by simulating Ai-and Bi-soliton solutions of the cylindrical Korteweg-de Vries(cKdV)equation.
文摘In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.
文摘In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.
基金supported by the NSFC(12101012)the PhD Scientific Research Start-up Foundation of Anhui Normal University.Zeng’s research was supported by the NSFC(11961160716,11871054,12131017).
文摘In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.
基金Supported by the National Natural Science Foundation of China(12261023,11861023)the Foundation of Science and Technology project of Guizhou Province of China([2018]5769-05)。
文摘On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].
基金Supported by the National Natural Science Foundation of China(11771216)the Key Research and Development Program of Jiangsu Province(Social Development)(BE2019725)the Qing Lan Project of Jiangsu Province。
文摘This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and Marchand(Control Cybernet.2011,40:971-988),we establish an observability inequality of the conservative problem,and then discuss the equivalence between the exponential stabilization of a dissipative system and the internal observational inequality of the corresponding conservative system.
基金Project supported by the National Natural Science Foundation of China (Grant No. 12161061)the Fundamental Research Funds for the Inner Mongolia University of Finance and Economics (Grant No. NCYWT23036)+2 种基金the Young Innovative and Entrepreneurial Talents of the Inner Mongolia Grassland Talents Project in 2022,Autonomous Region “Five Major Tasks” Research Special Project for the Inner Mongolia University of Finance and Economics in 2024 (Grant No. NCXWD2422)High Quality Research Achievement Cultivation Fund for the Inner Mongolia University of Finance and Economics in 2024 (Grant No. GZCG2426)the Talent Development Fund of Inner Mongolia Autonomous Region, China。
文摘Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained through the Hirota method and symbolic computation. Breather-like and bound-state solitons are constructed in which the signs of the imaginary parts of the complex wave numbers and the initial separations of the two parallel solitons are important factors for the interaction patterns. The periodic structures and position-induced phase shift of some solutions are introduced.
文摘The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane significantly influences the value of the geodetic coordinates.This research analyzes this influence and how it can contribute to their solutions.The study investigates the mathematical relation between them and presents an exact formula relating to the geodetic height and the ellipse equation.In addition,a heuristic formula for the relation between the geodetic height and the ellipse equation is proposed,which is independent of the geodetic latitude and has a relative accuracy better than 99.9 %.The calculation is stable,and the cost is low.
基金Project supported by the National Natural Science Foundation of China (Grant No.12172226)。
文摘The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.
文摘Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.
文摘The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.
基金supported by the National Natural Science Foundation of China (12001033)。
文摘We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the kinetic TCS equation for a particle ensemble and the Stokes equations for a fluid via a drag force.In this paper,we present a complete analysis of the existence of global-in-time strong solutions to the coupled model without any smallness restrictions on the initial data.
基金Project supported by the National Natural Science Foun dation of China(Grant No.11274398).
文摘A significant obstacle impeding the advancement of the time fractional Schrodinger equation lies in the challenge of determining its precise mathematical formulation.In order to address this,we undertake an exploration of the time fractional Schrodinger equation within the context of a non-Markovian environment.By leveraging a two-level atom as an illustrative case,we find that the choice to raise i to the order of the time derivative is inappropriate.In contrast to the conventional approach used to depict the dynamic evolution of quantum states in a non-Markovian environment,the time fractional Schrodinger equation,when devoid of fractional-order operations on the imaginary unit i,emerges as a more intuitively comprehensible framework in physics and offers greater simplicity in computational aspects.Meanwhile,we also prove that it is meaningless to study the memory of time fractional Schrodinger equation with time derivative 1<α≤2.It should be noted that we have not yet constructed an open system that can be fully described by the time fractional Schrodinger equation.This will be the focus of future research.Our study might provide a new perspective on the role of time fractional Schrodinger equation.