The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbit...The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.展开更多
The structure-preserving property, in both the time domain and the frequency domain, is an important index for evaluating validity of a numerical method. Even in the known structure-preserving methods such as the symp...The structure-preserving property, in both the time domain and the frequency domain, is an important index for evaluating validity of a numerical method. Even in the known structure-preserving methods such as the symplectic method, the inherent conser- vation law in the frequency domain is hardly conserved. By considering a mathematical pendulum model, a Stormer-Verlet scheme is first constructed in a Hamiltonian frame- work. The conservation law of the StSrmer-Verlet scheme is derived, including the total energy expressed in the time domain and periodicity in the frequency domain. To track the structure-preserving properties of the Stormer-Verlet scheme associated with the con- servation law, the motion of the mathematical pendulum is simulated with different time step lengths. The numerical results illustrate that the StSrmer-Verlet scheme can preserve the total energy of the model but cannot preserve periodicity at all. A phase correction is performed for the StSrmer-Verlet scheme. The results imply that the phase correction can improve the conservative property of periodicity of the Stormer-Verlet scheme.展开更多
The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background,in which the longitudinal wave dissipation determines some important performances of the slender struc...The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background,in which the longitudinal wave dissipation determines some important performances of the slender structure.To reproduce the longitudinal wave dissipation effects on an elastic rod with a variable cross-section,a structure-preserving approach is developed based on the dynamic symmetry breaking theory.For the dynamic model controlling the longitudinal wave propagating in the elastic rod with the variable cross-section,the approximate multi-symplectic form is deduced based on the multi-symplectic method,and the expression of the local energy dissipation for the longitudinal wave propagating in the rod is presented,referring to the dynamic symmetry breaking theory.A structure-preserving method focusing on the residual of the multi-symplectic structure and the local energy dissipation of the dynamic model is constructed by using the midpoint difference discrete method.The longitudinal wave propagating in an elastic rod fixed at one end is simulated,and the local/total energy dissipations of the longitudinal wave are investigated by the constructed structure-preserving scheme in two typical cases in detail.展开更多
A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the ...A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the symplectic structure of a finite dimension Hamiltonian system, and a multi-symplectic method can preserve the multi-symplectic structure of an infinite dimension Hamiltonian system. In this paper, the structure-preserving properties of three differential schemes for an oscillator system are investigated in detail. Both the theoretical results and the numerical results show that the results obtained by the standard forward Euler scheme lost all the three geometric properties of the oscillator system, i.e., periodicity, boundedness, and total energy, the symplectic scheme can preserve the first two geometric properties of the oscillator system, and the St？rmer-Verlet scheme can preserve the three geometric properties of the oscillator system well. In addition, the relative errors for the Hamiltonian function of the symplectic scheme increase with the increase in the step length, suggesting that the symplectic scheme possesses good structure-preserving properties only if the step length is small enough.展开更多
In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, b...In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradient- Hamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms （SPAs） for the Duffing equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffing equation. However, considering the explicit Runge Kutta methods, this paper finds that there is an error term of order p＋l for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method （a second-order Runge-Kutta method） are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when e is small or equal to zero.展开更多
The current structure-preserving theory, including the symplectic method and the multisymplectic method, pays most attention on the conservative properties of the continuous systems because that the conservative prope...The current structure-preserving theory, including the symplectic method and the multisymplectic method, pays most attention on the conservative properties of the continuous systems because that the conservative properties of the conservative systems can be formulated in the mathematical form. But, the nonconservative characteristics are the nature of the systems existing in engineering. In this letter, the structure-preserving approach for the infinite dimensional nonconservative systems is proposed based on the generalized multi-symplectic method to broaden the application fields of the current structure-preserving idea. In the numerical examples,two nonconservative factors, including the strong excitation on the string and the impact on the cantilever, are considered respectively. The vibrations of the string and the cantilever are investigated by the structure-preserving approach and the good long-time numerical behaviors as well as the high numerical precision of which are illustrated by the numerical results presented.展开更多
ABSTRACT The main idea of the structure-preserving method is to preserve the intrinsic geo- metric properties of the continuous system as much as possible in numerical algorithm design. The geometric constraint in the...ABSTRACT The main idea of the structure-preserving method is to preserve the intrinsic geo- metric properties of the continuous system as much as possible in numerical algorithm design. The geometric constraint in the multi-body systems, one of the difficulties in the numerical methods that are proposed for the multi-body systems, can also be regarded as a geometric property of the multi-body systems. Based on this idea, the symplectic precise integration method is applied in this paper to analyze the kinematics problem of folding and unfolding process of nose under- carriage. The Lagrange governing equation is established for the folding and unfolding process of nose undercarriage with the generalized defined displacements firstly. And then, the constrained Hamiltonian canonical form is derived from the Lagrange governing equation based on the Hamil- tonian variational principle. Finally, the symplectic precise integration scheme is used to simulate the kinematics process of nose undercarriage during folding and unfolding described by the con- strained Hamiltonian canonical formulation. From the numerical results, it can be concluded that the geometric constraint of the undercarriage system can be preserved well during the numerical simulation on the folding and unfolding process of undercarriage using the symplectic precise integration method.展开更多
Recent development of structure-preserving geometric particle-in-cell （PIC） algorithms for Vlasov-Maxwell systems is summarized. With the arrival of 100 petaflop and exaflop computing power, it is now possible to ca...Recent development of structure-preserving geometric particle-in-cell （PIC） algorithms for Vlasov-Maxwell systems is summarized. With the arrival of 100 petaflop and exaflop computing power, it is now possible to carry out direct simulations of multi-scale plasma dynamics based on first-principles. However, standard algorithms currently adopted by the plasma physics community do not possess the long-term accuracy and fidelity required for these large-scale simulations. This is because conventional simulation algorithms are based on numerically solving the underpinning differential （or integro-differential） equations, and the algorithms used in general do not preserve the geometric and physical structures of the systems, such as the local energy-momentum conservation law, the symplectic structure, and the gauge symmetry. As a consequence, numerical errors accumulate coherently with time and long-term simulation results are not reliable. To overcome this difficulty and to harness the power of exascale computers, a new generation of structure-preserving geometric PIC algorithms have been developed. This new generation of algorithms utilizes modem mathematical techniques, such as discrete manifolds, interpolating differential forms, and non-canonical symplectic integrators, to ensure gauge symmetry, space-time symmetry and the conservation of charge, energy-momentum, and the symplectic structure. These highly desired properties are difficult to achieve using the conventional PIC algorithms. In addition to summarizing the recent development and demonstrating practical implementations, several new results are also presented, including a structure-preserving geometric relativistic PIC algorithm, the proof of the correspondence between discrete gauge symmetry and discrete charge conservation law, and a reformulation of the explicit non-canonical symplectic algorithm for the discrete Poisson bracket using the variational approach. Numerical examples are given to verify the advantages of the structure- preserving geometric PIC algorithms in comparison with the conventional PIC methods.展开更多
The variational calculus of time-scale non-shifted systems includes both the traditional continuous and traditional significant discrete variational calculus.Not only can the combination ofand∇derivatives be beneficia...The variational calculus of time-scale non-shifted systems includes both the traditional continuous and traditional significant discrete variational calculus.Not only can the combination ofand∇derivatives be beneficial to obtaining higher convergence order in numerical analysis,but also it prompts the timescale numerical computational scheme to have good properties,for instance,structure-preserving.In this letter,a structure-preserving algorithm for time-scale non-shifted Hamiltonian systems is proposed.By using the time-scale discrete variational method and calculus theory,and taking a discrete time scale in the variational principle of non-shifted Hamiltonian systems,the corresponding discrete Hamiltonian principle can be obtained.Furthermore,the time-scale discrete Hamilton difference equations,Noether theorem,and the symplectic scheme of discrete Hamiltonian systems are obtained.Finally,taking the Kepler problem and damped oscillator for time-scale non-shifted Hamiltonian systems as examples,they show that the time-scale discrete variational method is a structure-preserving algorithm.The new algorithm not only provides a numerical method for solving time-scale non-shifted dynamic equations but can be calculated with variable step sizes to improve the computational speed.展开更多
The Pfaff-Birkhoff variational principle is discretized, and based on the discrete variational principle the discrete Birkhoffian equations are obtained. Taking the discrete equations as an algorithm, the correspondin...The Pfaff-Birkhoff variational principle is discretized, and based on the discrete variational principle the discrete Birkhoffian equations are obtained. Taking the discrete equations as an algorithm, the corresponding discrete flow is proved to be symplectic. That means the algorithm preserves the symplectic structure of Birkhofflan systems. Finally, simulation results of the given example indicate that structure-preserving algorithms have great advantage in stability and energy conserving.展开更多
In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct effi...In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct efficient and accurate time discretization scheme for a large class of Hamiltonian partial differential equations(PDEs).The proposed scheme is a linear system,and can be solved more efficient than the original energy-preserving ex-ponential integrator scheme which usually needs nonlinear iterations.Various experiments are performed to verify the conservation,efficiency and good performance at relatively large time step in long time computations.展开更多
In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries （KdV） equation, including eight multi-symplec...In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries （KdV） equation, including eight multi-symplectic algorithms, eight local energy-conserving algo- rithms and eight local momentum-conserving algorithms. Among these algorithms, some have been discussed and widely used while the most are new. The outstanding advantage of these proposed algorithms is that they conserve the local structures in any time-space re- gion exactly. Therefore, the local structure-preserving algorithms overcome the restriction of global structure-preserving algorithms on the boundary conditions. Numerical experiments are conducted to show the performance of the proposed methods. Moreover, the unified framework can be easily applied to many other equations.展开更多
This paper gives several structure-preserving schemes for the Degasperis-Procesi equation which has bi-Hamiltonian structures consisted of both complex and non-local Hamiltonian differential operators. For this sake, ...This paper gives several structure-preserving schemes for the Degasperis-Procesi equation which has bi-Hamiltonian structures consisted of both complex and non-local Hamiltonian differential operators. For this sake, few structure-preserving schemes have been proposed so far. In our work, based on one of the bi-Hamiltonian structures, a symplectic scheme and two new energy-preserving schemes are constructed. The symplecticity comes straightly from the application of the implicit midpoint method on the semi-discrete system which is proved to remain Hamiltonian, while the energy conservation is derived by the combination of the averaged vector field method of second and fourth order, respectively. Some numerical results are presented to show that the three schemes do have the advantages in numerical stability, accuracy in long time computing and ability to preserve the invariants of the DP equation.展开更多
For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a s...For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a structure-preserving method based on a discrete variational derivative method.Furthermore,to prevent the vertex concentration that may lead to numerical instability,we discretely introduce Deckelnick’s tangential velocity.Here,a modification term is introduced in the process of adding tangential velocity.This modified term enables the method to reproduce the equations’properties while preventing vertex concentration.Numerical experiments demonstrate that the proposed approach captures the equations’properties with high accuracy and avoids the concentration of vertices.展开更多
This paper introduces two novel conformal structure-preserving algorithms for solving the coupled damped nonlinear Schr¨odinger(CDNLS)system,which are based on the conformal multi-symplectic Hamiltonian formulati...This paper introduces two novel conformal structure-preserving algorithms for solving the coupled damped nonlinear Schr¨odinger(CDNLS)system,which are based on the conformal multi-symplectic Hamiltonian formulation and its conformal conservation laws.The proposed algorithms can preserve corresponding conformal multi-symplectic conservation lawand conformalmomentum conservation lawin any local time-space region,respectively.Moreover,it is further shown that the algorithms admit the conformal charge conservation law,and exactly preserve the dissipation rate of charge under appropriate boundary conditions.Numerical experiments are presented to demonstrate the conformal properties and effectiveness of the proposed algorithms during long-time numerical simulations and validate the analysis.展开更多
Presents a study which examined the structure-preserving algorithms to phase space volume for linear dynamical systems. Preservation of phase space volume for source-free dynamical systems; Description of a volume-pre...Presents a study which examined the structure-preserving algorithms to phase space volume for linear dynamical systems. Preservation of phase space volume for source-free dynamical systems; Description of a volume-preserving scheme for linear system with canonical form; Information on structure-preserving schemes for linear dynamical systems.展开更多
The rheological properties and limited flow velocities of solvent-free nanofluids are crucial for their technologically significant applications.In particular,the flow in a solvent-free nanofluid system is steady only...The rheological properties and limited flow velocities of solvent-free nanofluids are crucial for their technologically significant applications.In particular,the flow in a solvent-free nanofluid system is steady only when the flow velocity is lower than a critical value.In this paper,we establish a rigid-flexible dynamic model to investigate the existence of the upper bound on the steady flow velocities for three solvent-free nanofluid systems.Then,the effects of the structural parameters on the upper bound on the steady flow velocities are examined with the proposed structure-preserving method.It is found that each of these solvent-free nanofluid systems has an upper bound on the steady flow velocity,which exhibits distinct dependence on their structural parameters,such as the graft density of branch chains and the size of the cores.In addition,among the three types of solvent-free nanofluids,the magnetic solvent-free nanofluid poses the largest upper bound on the steady flow velocity,demonstrating that it is a better choice when a large flow velocity is required in real applications.展开更多
We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations ...We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.展开更多
The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on t...The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on the dynamic behaviors of the ultra-large spatial structure.Taking the spatial beam as an example,the minimum control energy defined by the difference between the initial total energy and the final total energy in the assumed stable attitude state of the beam is investigated by the structure-preserving method proposed in our previous studies in two cases:the spatial beam considering the flexibility as well as the damping effect,and the spatial beam ignoring both the flexibility and the damping effect.In the numerical experiments,the assumed simulation interval of three months is evaluated on whether or not it is long enough for the spatial flexible damping beam to arrive at the assumed stable attitude state.And then,taking the initial attitude angle and the initial attitude angle velocity as the independent variables,respectively,the minimum control energies of the mentioned two cases are investigated in detail.From the numerical results,the following conclusions can be obtained.With the fixed initial attitude angle velocity,the minimum control energy of the spatial flexible damping beam is higher than that of the spatial rigid beam when the initial attitude angle is close to or far away from the stable attitude state.With the fixed initial attitude angle,ignoring the flexibility and the damping effect will underestimate the minimum control energy of the spatial beam.展开更多
connecting wires are the main manifestations of the coupling dynamic effects on the orbit evolution,the attitude adjusting and the flexible vibration of the tethered satellite system.To investigate attitude evolution ...connecting wires are the main manifestations of the coupling dynamic effects on the orbit evolution,the attitude adjusting and the flexible vibration of the tethered satellite system.To investigate attitude evolution of the tethered system and the mechanical energy transfer/loss characteristics between the bus system and the solar sail via the connecting wires,a structure-preserving method is developed in this paper.Simplifying the tethered satellite system as a composite structure consisting of a particle and a flexible thin panel connected by four special springs,the dynamic model is deduced via the Hamiltonian variational principle firstly.Then,a structure-preserving approach that connects the symplectic Runge-Kutta method and the multi-symplectic method is developed.The excellent structure-preserving property of the numerical scheme constructed is presented to illustrate the credibility of the numerical results obtained by the constructed structure-preserving approach.From the numerical results on the mechanical energy transfer/loss in the composite structure,it can be found that the mechanical energy transfer tendency in the tethered system is dependent of the initial attitude angle of the system while the total mechanical energy loss of the system is almost independent of the initial attitude angle.In addition,the special stiffness range of the spring is found in the attitude angle evolution of the system,which provides a structural parameter design window for the connecting wires,that is,the duration needed to arrive the stable attitude is short when the stiffness of the wire is designed in this special range.展开更多
基金supported by National Natural Science Foundation of China (Nos. 11975068 and 11925501)the National Key R&D Program of China (No. 2022YFE03090000)the Fundamental Research Funds for the Central Universities (No. DUT22ZD215)。
文摘The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.
基金the National Natural Science Foundation of China(Nos.11672241,11372253,and 11432010)the Astronautics Supporting Technology Foundation of China(No.2015-HT-XGD)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(Nos.GZ1312 and GZ1605)
文摘The structure-preserving property, in both the time domain and the frequency domain, is an important index for evaluating validity of a numerical method. Even in the known structure-preserving methods such as the symplectic method, the inherent conser- vation law in the frequency domain is hardly conserved. By considering a mathematical pendulum model, a Stormer-Verlet scheme is first constructed in a Hamiltonian frame- work. The conservation law of the StSrmer-Verlet scheme is derived, including the total energy expressed in the time domain and periodicity in the frequency domain. To track the structure-preserving properties of the Stormer-Verlet scheme associated with the con- servation law, the motion of the mathematical pendulum is simulated with different time step lengths. The numerical results illustrate that the StSrmer-Verlet scheme can preserve the total energy of the model but cannot preserve periodicity at all. A phase correction is performed for the StSrmer-Verlet scheme. The results imply that the phase correction can improve the conservative property of periodicity of the Stormer-Verlet scheme.
基金Projected supported by the National Natural Science Foundation of China(Nos.11872303,12172281,11972284)the Fund for Distinguished Young Scholars of Shaanxi Province of China(No.2019JC-29)+2 种基金the Foundation Strengthening Programme Technical Area Fund(No.2021-JCJQ-JJ-0565)the Fund of the Youth Innovation Team of Shaanxi Universitiesthe Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(No.GZ19103)。
文摘The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background,in which the longitudinal wave dissipation determines some important performances of the slender structure.To reproduce the longitudinal wave dissipation effects on an elastic rod with a variable cross-section,a structure-preserving approach is developed based on the dynamic symmetry breaking theory.For the dynamic model controlling the longitudinal wave propagating in the elastic rod with the variable cross-section,the approximate multi-symplectic form is deduced based on the multi-symplectic method,and the expression of the local energy dissipation for the longitudinal wave propagating in the rod is presented,referring to the dynamic symmetry breaking theory.A structure-preserving method focusing on the residual of the multi-symplectic structure and the local energy dissipation of the dynamic model is constructed by using the midpoint difference discrete method.The longitudinal wave propagating in an elastic rod fixed at one end is simulated,and the local/total energy dissipations of the longitudinal wave are investigated by the constructed structure-preserving scheme in two typical cases in detail.
基金supported by the National Natural Science Foundation of China(Nos.1117223911002115+4 种基金and 11372253)Doctoral Program Foundation of Education Ministry of China(No.20126102110023)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(Nos.GZ0802 and GZ1312)the Special Fund for Basic Scientific Researchof Central CollegesChang’an University(No.CHD2011JC040)
文摘A numerical method for the Hamiltonian system is required to preserve some structure-preserving properties. The current structure-preserving method satisfies the requirements that a symplectic method can preserve the symplectic structure of a finite dimension Hamiltonian system, and a multi-symplectic method can preserve the multi-symplectic structure of an infinite dimension Hamiltonian system. In this paper, the structure-preserving properties of three differential schemes for an oscillator system are investigated in detail. Both the theoretical results and the numerical results show that the results obtained by the standard forward Euler scheme lost all the three geometric properties of the oscillator system, i.e., periodicity, boundedness, and total energy, the symplectic scheme can preserve the first two geometric properties of the oscillator system, and the St？rmer-Verlet scheme can preserve the three geometric properties of the oscillator system well. In addition, the relative errors for the Hamiltonian function of the symplectic scheme increase with the increase in the step length, suggesting that the symplectic scheme possesses good structure-preserving properties only if the step length is small enough.
基金Project supported by the National Natural Science Foundation of China (Grant No 10572021)the Doctoral Programme Foundation of Institute of Higher Education of China (Grant No 20040007022)
文摘In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradient- Hamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms （SPAs） for the Duffing equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffing equation. However, considering the explicit Runge Kutta methods, this paper finds that there is an error term of order p＋l for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method （a second-order Runge-Kutta method） are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when e is small or equal to zero.
基金supported by the National Natural Science Foundation of China (Grant 11672241)the Seed Foundation of Qian Xuesen Laboratory of Space Technologythe Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment (Grant GZ1605)
文摘The current structure-preserving theory, including the symplectic method and the multisymplectic method, pays most attention on the conservative properties of the continuous systems because that the conservative properties of the conservative systems can be formulated in the mathematical form. But, the nonconservative characteristics are the nature of the systems existing in engineering. In this letter, the structure-preserving approach for the infinite dimensional nonconservative systems is proposed based on the generalized multi-symplectic method to broaden the application fields of the current structure-preserving idea. In the numerical examples,two nonconservative factors, including the strong excitation on the string and the impact on the cantilever, are considered respectively. The vibrations of the string and the cantilever are investigated by the structure-preserving approach and the good long-time numerical behaviors as well as the high numerical precision of which are illustrated by the numerical results presented.
基金supported by the National Natural Science Foundation of China(Nos.11672241,11372253 and 11432010)the Shanxi National Science Foundation(No.2015JM1026)+3 种基金the Astronautics Supporting Technology Foundation of China(No.2015-HT-XGD)111 Project(No.B07050) to the Northwestern Polytechnical Universitythe fund of the State Key Laboratory of Solidification Processing in NWPU(No.SKLSP201643)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(No.GZ1605)
文摘ABSTRACT The main idea of the structure-preserving method is to preserve the intrinsic geo- metric properties of the continuous system as much as possible in numerical algorithm design. The geometric constraint in the multi-body systems, one of the difficulties in the numerical methods that are proposed for the multi-body systems, can also be regarded as a geometric property of the multi-body systems. Based on this idea, the symplectic precise integration method is applied in this paper to analyze the kinematics problem of folding and unfolding process of nose under- carriage. The Lagrange governing equation is established for the folding and unfolding process of nose undercarriage with the generalized defined displacements firstly. And then, the constrained Hamiltonian canonical form is derived from the Lagrange governing equation based on the Hamil- tonian variational principle. Finally, the symplectic precise integration scheme is used to simulate the kinematics process of nose undercarriage during folding and unfolding described by the con- strained Hamiltonian canonical formulation. From the numerical results, it can be concluded that the geometric constraint of the undercarriage system can be preserved well during the numerical simulation on the folding and unfolding process of undercarriage using the symplectic precise integration method.
基金supported by National Natural Science Foundation of China (NSFC-11775219, 11775222, 11505186, 11575185 and 11575186)the National Key Research and Development Program (2016YFA0400600, 2016YFA0400601 and 2016YFA0400602)+3 种基金the ITER-China Program (2015GB111003, 2014GB124005)Chinese Scholar Council (201506340103)China Postdoctoral Science Foundation (2017LH002)the GeoA lgorithmic Plasma Simulator (GAPS) Project
文摘Recent development of structure-preserving geometric particle-in-cell （PIC） algorithms for Vlasov-Maxwell systems is summarized. With the arrival of 100 petaflop and exaflop computing power, it is now possible to carry out direct simulations of multi-scale plasma dynamics based on first-principles. However, standard algorithms currently adopted by the plasma physics community do not possess the long-term accuracy and fidelity required for these large-scale simulations. This is because conventional simulation algorithms are based on numerically solving the underpinning differential （or integro-differential） equations, and the algorithms used in general do not preserve the geometric and physical structures of the systems, such as the local energy-momentum conservation law, the symplectic structure, and the gauge symmetry. As a consequence, numerical errors accumulate coherently with time and long-term simulation results are not reliable. To overcome this difficulty and to harness the power of exascale computers, a new generation of structure-preserving geometric PIC algorithms have been developed. This new generation of algorithms utilizes modem mathematical techniques, such as discrete manifolds, interpolating differential forms, and non-canonical symplectic integrators, to ensure gauge symmetry, space-time symmetry and the conservation of charge, energy-momentum, and the symplectic structure. These highly desired properties are difficult to achieve using the conventional PIC algorithms. In addition to summarizing the recent development and demonstrating practical implementations, several new results are also presented, including a structure-preserving geometric relativistic PIC algorithm, the proof of the correspondence between discrete gauge symmetry and discrete charge conservation law, and a reformulation of the explicit non-canonical symplectic algorithm for the discrete Poisson bracket using the variational approach. Numerical examples are given to verify the advantages of the structure- preserving geometric PIC algorithms in comparison with the conventional PIC methods.
基金This work was supported by the National Natural Science Foundation of China(Nos.11972241,11572212)the Natural Science Foundation of Jiangsu Province(No.BK20191454)the Postgraduate Research&Practice Innovation Program of Jiangsu Province(No.KYCX20_0251).
文摘The variational calculus of time-scale non-shifted systems includes both the traditional continuous and traditional significant discrete variational calculus.Not only can the combination ofand∇derivatives be beneficial to obtaining higher convergence order in numerical analysis,but also it prompts the timescale numerical computational scheme to have good properties,for instance,structure-preserving.In this letter,a structure-preserving algorithm for time-scale non-shifted Hamiltonian systems is proposed.By using the time-scale discrete variational method and calculus theory,and taking a discrete time scale in the variational principle of non-shifted Hamiltonian systems,the corresponding discrete Hamiltonian principle can be obtained.Furthermore,the time-scale discrete Hamilton difference equations,Noether theorem,and the symplectic scheme of discrete Hamiltonian systems are obtained.Finally,taking the Kepler problem and damped oscillator for time-scale non-shifted Hamiltonian systems as examples,they show that the time-scale discrete variational method is a structure-preserving algorithm.The new algorithm not only provides a numerical method for solving time-scale non-shifted dynamic equations but can be calculated with variable step sizes to improve the computational speed.
基金Supported by the National Natural Science Foundation of China (10932002,10972031)
文摘The Pfaff-Birkhoff variational principle is discretized, and based on the discrete variational principle the discrete Birkhoffian equations are obtained. Taking the discrete equations as an algorithm, the corresponding discrete flow is proved to be symplectic. That means the algorithm preserves the symplectic structure of Birkhofflan systems. Finally, simulation results of the given example indicate that structure-preserving algorithms have great advantage in stability and energy conserving.
基金supported by the National Natural Science Foundation of China(Grant Nos.12171245,11971416,11971242,12301508)by the Natural Science Foundation of Henan Province(Grant No.222300420280)+1 种基金by the Natural Science Foundation of Hunan Province(Grant No.2023JJ40656)by the Scientific Research Fund of Xuchang University(Grant No.2024ZD010).
文摘In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct efficient and accurate time discretization scheme for a large class of Hamiltonian partial differential equations(PDEs).The proposed scheme is a linear system,and can be solved more efficient than the original energy-preserving ex-ponential integrator scheme which usually needs nonlinear iterations.Various experiments are performed to verify the conservation,efficiency and good performance at relatively large time step in long time computations.
文摘In this paper, based on the concatenating method, we present a unified framework to construct a series of local structure-preserving algorithms for the Korteweg-de Vries （KdV） equation, including eight multi-symplectic algorithms, eight local energy-conserving algo- rithms and eight local momentum-conserving algorithms. Among these algorithms, some have been discussed and widely used while the most are new. The outstanding advantage of these proposed algorithms is that they conserve the local structures in any time-space re- gion exactly. Therefore, the local structure-preserving algorithms overcome the restriction of global structure-preserving algorithms on the boundary conditions. Numerical experiments are conducted to show the performance of the proposed methods. Moreover, the unified framework can be easily applied to many other equations.
基金the National Natural Science Foundation of China (Grant No. 11771213)the National Key Research and Development Project of China (Grant No. 2016YFC0600310)the Major Projects of Natural Sciences of University in Jiangsu Province of China (Grant No. 15KJA110002).
文摘This paper gives several structure-preserving schemes for the Degasperis-Procesi equation which has bi-Hamiltonian structures consisted of both complex and non-local Hamiltonian differential operators. For this sake, few structure-preserving schemes have been proposed so far. In our work, based on one of the bi-Hamiltonian structures, a symplectic scheme and two new energy-preserving schemes are constructed. The symplecticity comes straightly from the application of the implicit midpoint method on the semi-discrete system which is proved to remain Hamiltonian, while the energy conservation is derived by the combination of the averaged vector field method of second and fourth order, respectively. Some numerical results are presented to show that the three schemes do have the advantages in numerical stability, accuracy in long time computing and ability to preserve the invariants of the DP equation.
基金This work was supported by JSPS KAKENHI Grant Nos.19K14590,21K18301,Japan.
文摘For fourth-order geometric evolution equations for planar curves with the dissipation of the bending energy,including the Willmore and the Helfrich flows,we consider a numerical approach.In this study,we construct a structure-preserving method based on a discrete variational derivative method.Furthermore,to prevent the vertex concentration that may lead to numerical instability,we discretely introduce Deckelnick’s tangential velocity.Here,a modification term is introduced in the process of adding tangential velocity.This modified term enables the method to reproduce the equations’properties while preventing vertex concentration.Numerical experiments demonstrate that the proposed approach captures the equations’properties with high accuracy and avoids the concentration of vertices.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11501570,91530106 and 11571366)Research Fund ofNUDT(Grant No.JC15-02-02)the fund from HPCL.
文摘This paper introduces two novel conformal structure-preserving algorithms for solving the coupled damped nonlinear Schr¨odinger(CDNLS)system,which are based on the conformal multi-symplectic Hamiltonian formulation and its conformal conservation laws.The proposed algorithms can preserve corresponding conformal multi-symplectic conservation lawand conformalmomentum conservation lawin any local time-space region,respectively.Moreover,it is further shown that the algorithms admit the conformal charge conservation law,and exactly preserve the dissipation rate of charge under appropriate boundary conditions.Numerical experiments are presented to demonstrate the conformal properties and effectiveness of the proposed algorithms during long-time numerical simulations and validate the analysis.
文摘Presents a study which examined the structure-preserving algorithms to phase space volume for linear dynamical systems. Preservation of phase space volume for source-free dynamical systems; Description of a volume-preserving scheme for linear system with canonical form; Information on structure-preserving schemes for linear dynamical systems.
基金Project supported by the National Natural Science Foundation of China(No.12172281)the Fund of the Science and Technology Innovation Team of Shaanxi Province of China(No.2022TD-61)+2 种基金the Open Project of State Key Laboratory of Performance Monitoring and Protecting of Rail Transit InfrastructureEast China Jiaotong University(No.HJGZ2023102)the Shaanxi Province Key Research and Development Project(No.2024SFYBXM-531)。
文摘The rheological properties and limited flow velocities of solvent-free nanofluids are crucial for their technologically significant applications.In particular,the flow in a solvent-free nanofluid system is steady only when the flow velocity is lower than a critical value.In this paper,we establish a rigid-flexible dynamic model to investigate the existence of the upper bound on the steady flow velocities for three solvent-free nanofluid systems.Then,the effects of the structural parameters on the upper bound on the steady flow velocities are examined with the proposed structure-preserving method.It is found that each of these solvent-free nanofluid systems has an upper bound on the steady flow velocity,which exhibits distinct dependence on their structural parameters,such as the graft density of branch chains and the size of the cores.In addition,among the three types of solvent-free nanofluids,the magnetic solvent-free nanofluid poses the largest upper bound on the steady flow velocity,demonstrating that it is a better choice when a large flow velocity is required in real applications.
文摘We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.
基金The research is supported by the National Natural Science Foundation of China(11672241,11972284,11432010)Fund for Distinguished Young Scholars of Shaanxi Province(2019JC-29)Fund of the Youth Innovation Team of Shaanxi Universities,the Seed Foundation of Qian Xuesen Laboratory of Space Technology,and the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(GZ1605).
文摘The dynamic analysis on the ultra-large spatial structure can be simplified drastically by ignoring the flexibility and damping of the structure.However,these simplifications will result in the erroneous estimate on the dynamic behaviors of the ultra-large spatial structure.Taking the spatial beam as an example,the minimum control energy defined by the difference between the initial total energy and the final total energy in the assumed stable attitude state of the beam is investigated by the structure-preserving method proposed in our previous studies in two cases:the spatial beam considering the flexibility as well as the damping effect,and the spatial beam ignoring both the flexibility and the damping effect.In the numerical experiments,the assumed simulation interval of three months is evaluated on whether or not it is long enough for the spatial flexible damping beam to arrive at the assumed stable attitude state.And then,taking the initial attitude angle and the initial attitude angle velocity as the independent variables,respectively,the minimum control energies of the mentioned two cases are investigated in detail.From the numerical results,the following conclusions can be obtained.With the fixed initial attitude angle velocity,the minimum control energy of the spatial flexible damping beam is higher than that of the spatial rigid beam when the initial attitude angle is close to or far away from the stable attitude state.With the fixed initial attitude angle,ignoring the flexibility and the damping effect will underestimate the minimum control energy of the spatial beam.
基金was supported by the National Natural Science Foundation of China(Grants 11972284,11872303)the Fund for Distinguished Young Scholars of Shaanxi Province(2019JC-29)the Fund of the Youth Innovation Team of Shaanxi Universities,and the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment(Grant GZ19103).
文摘connecting wires are the main manifestations of the coupling dynamic effects on the orbit evolution,the attitude adjusting and the flexible vibration of the tethered satellite system.To investigate attitude evolution of the tethered system and the mechanical energy transfer/loss characteristics between the bus system and the solar sail via the connecting wires,a structure-preserving method is developed in this paper.Simplifying the tethered satellite system as a composite structure consisting of a particle and a flexible thin panel connected by four special springs,the dynamic model is deduced via the Hamiltonian variational principle firstly.Then,a structure-preserving approach that connects the symplectic Runge-Kutta method and the multi-symplectic method is developed.The excellent structure-preserving property of the numerical scheme constructed is presented to illustrate the credibility of the numerical results obtained by the constructed structure-preserving approach.From the numerical results on the mechanical energy transfer/loss in the composite structure,it can be found that the mechanical energy transfer tendency in the tethered system is dependent of the initial attitude angle of the system while the total mechanical energy loss of the system is almost independent of the initial attitude angle.In addition,the special stiffness range of the spring is found in the attitude angle evolution of the system,which provides a structural parameter design window for the connecting wires,that is,the duration needed to arrive the stable attitude is short when the stiffness of the wire is designed in this special range.
