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.展开更多
A new approach is proposed in this study for accountable capability improvement based on interpretable capability evaluation using the belief rule base(BRB).Firstly,a capability evaluation model is constructed and opt...A new approach is proposed in this study for accountable capability improvement based on interpretable capability evaluation using the belief rule base(BRB).Firstly,a capability evaluation model is constructed and optimized.Then,the key sub-capabilities are identified by quantitatively calculating the contributions made by each sub-capability to the overall capability.Finally,the overall capability is improved by optimizing the identified key sub-capabilities.The theoretical contributions of the proposed approach are as follows.(i)An interpretable capability evaluation model is constructed by employing BRB which can provide complete access to decision-makers.(ii)Key sub-capabilities are identified according to the quantitative contribution analysis results.(iii)Accountable capability improvement is carried out by only optimizing the identified key sub-capabilities.Case study results show that“Surveillance”,“Positioning”,and“Identification”are identified as key sub-capabilities with a summed contribution of 75.55%in an analytical and deducible fashion based on the interpretable capability evaluation model.As a result,the overall capability is improved by optimizing only the identified key sub-capabilities.The overall capability can be greatly improved from 59.20%to 81.80%with a minimum cost of 397.Furthermore,this paper also investigates how optimizing the BRB with more collected data would affect the evaluation results:only optimizing“Surveillance”and“Positioning”can also improve the overall capability to 81.34%with a cost of 370,which thus validates the efficiency of the proposed approach.展开更多
We systemically investigate optical trapping capability of a kind of tornado waves on Rayleigh particles.Such tornado waves are named as tornado circular Pearcey beams(TCPBs)and produced by combining two circular Pear...We systemically investigate optical trapping capability of a kind of tornado waves on Rayleigh particles.Such tornado waves are named as tornado circular Pearcey beams(TCPBs)and produced by combining two circular Pearcey beams with different radii.Our theoretical exploration delves into various aspects,including the propagation dynamics,energy flux,orbital angular momentum,trapping force,and torque characteristics of TCPBs.The results reveal that the orbital angular momentum,trapping force,and torque of these beams can be finely tuned through the judicious manipulation of their topological charges(l_(1)and l_(2)).Notably,we observe a precise control mechanism wherein the force diminishes with|l_(1)+l_(2)|and|l_(1)-l_(2)|,while the torque exhibits enhancement by decreasing solely with|l_(1)+l_(2)|or increasing with|l_(1)-l_(2)|.These results not only provide quantitative insights into the optical trapping performance of TCPBs but also serve as a valuable reference for the ongoing development of innovative photonic tools.展开更多
Concomitant with the advancement of contemporary medical technology,the significance of perioperative nursing has been increasingly accentuated,necessitating elevated standards for the pedagogy of perioperative nursin...Concomitant with the advancement of contemporary medical technology,the significance of perioperative nursing has been increasingly accentuated,necessitating elevated standards for the pedagogy of perioperative nursing.Presently,the PBL(problem-based learning)pedagogical approach,when integrated with CBL(case-based learning),has garnered considerable interest.An extensive literature review has been conducted to analyze the application of the PBL-CBL fusion in the education of perioperative nursing.Findings indicate that this integrative teaching methodology not only enhances students’theoretical knowledge,practical competencies,and collaborative skills but also contributes to the elevation of teaching quality.In conclusion,the PBL-CBL teaching approach holds immense potential for broader application in perioperative nursing education.Nevertheless,it is imperative to continually refine this combined pedagogical strategy to further enhance the caliber of perioperative nursing instruction and to cultivate a greater number of exceptional nursing professionals in the operating room setting.展开更多
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 emergence of Y6-type nonfullerene acceptors has greatly enhanced the power conversion efficiency(PCE)of organic solar cells(OSCs).However,which structural feature is responsible for the excellent photovoltaic perf...The emergence of Y6-type nonfullerene acceptors has greatly enhanced the power conversion efficiency(PCE)of organic solar cells(OSCs).However,which structural feature is responsible for the excellent photovoltaic performance is still under debate.In this study,two Y6-like acceptors BDOTP-1 and BDOTP-2 were designed.Different from previous Y6-type acceptors featuring an A–D–Aʹ–D–A structure,BDOTP-1,and BDOTP-2 have no electron-deficient Aʹfragment in the core unit.Instead,there is an electron-rich dibenzodioxine fragment in the core.Although this modification leads to a marked change in the molecular dipole moment,electrostatic potential,frontier orbitals,and energy levels,BDOTP acceptors retain similar three-dimensional packing capability as Y6-type acceptors due to the similar banana-shaped molecular configuration.BDOTP acceptors show good performance in OSCs.High PCEs of up to 18.51%(certified 17.9%)are achieved.This study suggests that the banana-shaped configuration instead of the A–D–Aʹ–D–A structure is likely to be the determining factor in realizing high photovoltaic performance.展开更多
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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry...The effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry rock.At present,how temporary plugging agents plug artificial fractures under high temperature remains unclear.In this paper,by establishing an improved experimental system for the evaluation of temporary plugging performance at high temperature,we clarified the effects of high temperature,injection rate,and fracture width on the pressure response and plugging efficiency of the fracture.The results revealed that the temporary plugging process of artificial fractures in hot dry rock can be divided into four main stages:the initial stage of temporary plugging,the bridging stage of the particles,the plugging formation stage,and the high-pressure dense plugging stage.As the temperature increases,the distribution distance of the temporary plugging agent,the number of pressure fluctuations,and the time required for crack plugging increases.Particularly,when the temperature increases by 100℃,the complete plugging time increases by 90.7%.展开更多
Antimony(Sb) is an attractive cathode for liquid metal batteries(LMBs) because of its high theoretical voltage and low cost.The main obstacles associated with the Sb-based cathodes are unsatisfactory energy density an...Antimony(Sb) is an attractive cathode for liquid metal batteries(LMBs) because of its high theoretical voltage and low cost.The main obstacles associated with the Sb-based cathodes are unsatisfactory energy density and poor rate-capability.Herein,we propose a novel Sb_(64)Cu_(36)cathode that effectively tackles these issues.The Sb_(64)Cu_(36)(melting point:525℃) cathode presents a novel lithiation mechanism involving sequentially the generation of Li_(2)CuSb,the formation of Li_(3)Sb,and the conversion reaction of Li_(2)CuSb to Li_(3)Sb and Cu.The generated intermetallic compounds show a unique microstructure of the upper floated Li_(2)CuSb layer and the below cross-linked structure with interpenetrated Li_(2)CuSb and Li_(3)Sb phases.Compared with Li_(3)Sb,the lower Li migration energy barrier(0.188 eV) of Li_(2)CuSb significantly facilitates the lithium diffusion across the intermediate compounds and accelerates the reaction kinetics.Consequently,the Li‖Sb_(64)Cu_(36)cell delivers a more excellent electrochemical performance(energy density:353 W h kg^(-1)at 0.4 A cm^(-2);rate capability:0.59 V at 2.0 A cm^(-2)),and a much lower energy storage cost of only 38.45 $ kW h^(-1)than other previously reported Sb-based LMBs.This work provides a novel cathode design concept for the development of high-performance LMBs in applications for large-scale energy storage.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(72471067,72431011,72471238,72231011,62303474,72301286)the Fundamental Research Funds for the Provincial Universities of Zhejiang(GK239909299001-010).
文摘A new approach is proposed in this study for accountable capability improvement based on interpretable capability evaluation using the belief rule base(BRB).Firstly,a capability evaluation model is constructed and optimized.Then,the key sub-capabilities are identified by quantitatively calculating the contributions made by each sub-capability to the overall capability.Finally,the overall capability is improved by optimizing the identified key sub-capabilities.The theoretical contributions of the proposed approach are as follows.(i)An interpretable capability evaluation model is constructed by employing BRB which can provide complete access to decision-makers.(ii)Key sub-capabilities are identified according to the quantitative contribution analysis results.(iii)Accountable capability improvement is carried out by only optimizing the identified key sub-capabilities.Case study results show that“Surveillance”,“Positioning”,and“Identification”are identified as key sub-capabilities with a summed contribution of 75.55%in an analytical and deducible fashion based on the interpretable capability evaluation model.As a result,the overall capability is improved by optimizing only the identified key sub-capabilities.The overall capability can be greatly improved from 59.20%to 81.80%with a minimum cost of 397.Furthermore,this paper also investigates how optimizing the BRB with more collected data would affect the evaluation results:only optimizing“Surveillance”and“Positioning”can also improve the overall capability to 81.34%with a cost of 370,which thus validates the efficiency of the proposed approach.
基金Project supported by the National Natural Science Foundation of China(Grant No.11604058)the Guangxi Natural Science Foundation(Grant Nos.2020GXNSFAA297041 and 2023JJA110112)the Innovation Project of Guangxi Graduate Education(Grant No.YCSW2023083)。
文摘We systemically investigate optical trapping capability of a kind of tornado waves on Rayleigh particles.Such tornado waves are named as tornado circular Pearcey beams(TCPBs)and produced by combining two circular Pearcey beams with different radii.Our theoretical exploration delves into various aspects,including the propagation dynamics,energy flux,orbital angular momentum,trapping force,and torque characteristics of TCPBs.The results reveal that the orbital angular momentum,trapping force,and torque of these beams can be finely tuned through the judicious manipulation of their topological charges(l_(1)and l_(2)).Notably,we observe a precise control mechanism wherein the force diminishes with|l_(1)+l_(2)|and|l_(1)-l_(2)|,while the torque exhibits enhancement by decreasing solely with|l_(1)+l_(2)|or increasing with|l_(1)-l_(2)|.These results not only provide quantitative insights into the optical trapping performance of TCPBs but also serve as a valuable reference for the ongoing development of innovative photonic tools.
文摘Concomitant with the advancement of contemporary medical technology,the significance of perioperative nursing has been increasingly accentuated,necessitating elevated standards for the pedagogy of perioperative nursing.Presently,the PBL(problem-based learning)pedagogical approach,when integrated with CBL(case-based learning),has garnered considerable interest.An extensive literature review has been conducted to analyze the application of the PBL-CBL fusion in the education of perioperative nursing.Findings indicate that this integrative teaching methodology not only enhances students’theoretical knowledge,practical competencies,and collaborative skills but also contributes to the elevation of teaching quality.In conclusion,the PBL-CBL teaching approach holds immense potential for broader application in perioperative nursing education.Nevertheless,it is imperative to continually refine this combined pedagogical strategy to further enhance the caliber of perioperative nursing instruction and to cultivate a greater number of exceptional nursing professionals in the operating room setting.
基金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.
基金the open research fund of the Songshan Lake Materials Laboratory(2021SLABFK02)the National Key Research and Development Program of China(2017YFA0206600)the National Natural Science Foundation of China(51922032 and 21961160720).
文摘The emergence of Y6-type nonfullerene acceptors has greatly enhanced the power conversion efficiency(PCE)of organic solar cells(OSCs).However,which structural feature is responsible for the excellent photovoltaic performance is still under debate.In this study,two Y6-like acceptors BDOTP-1 and BDOTP-2 were designed.Different from previous Y6-type acceptors featuring an A–D–Aʹ–D–A structure,BDOTP-1,and BDOTP-2 have no electron-deficient Aʹfragment in the core unit.Instead,there is an electron-rich dibenzodioxine fragment in the core.Although this modification leads to a marked change in the molecular dipole moment,electrostatic potential,frontier orbitals,and energy levels,BDOTP acceptors retain similar three-dimensional packing capability as Y6-type acceptors due to the similar banana-shaped molecular configuration.BDOTP acceptors show good performance in OSCs.High PCEs of up to 18.51%(certified 17.9%)are achieved.This study suggests that the banana-shaped configuration instead of the A–D–Aʹ–D–A structure is likely to be the determining factor in realizing high photovoltaic performance.
基金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.
基金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 (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 (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.
基金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 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.
基金supported financially by the Beijing Natural Science Foundation Project(No.3222030)the National Natural Science Foundation of China(No.51936001,No.52274002 and No.52192622)+1 种基金the PetroChina Science and Technology Innovation Foundation Project(2021DQ02–0201)Award Cultivation Foundation from Beijing Institute of Petrochemical Technology(No.BIPTACF-002).
文摘The effective plugging of artificial fractures is key to the success of temporary plugging and diverting fracturing technology,which is one of the most promising ways to improve the heat recovery efficiency of hot dry rock.At present,how temporary plugging agents plug artificial fractures under high temperature remains unclear.In this paper,by establishing an improved experimental system for the evaluation of temporary plugging performance at high temperature,we clarified the effects of high temperature,injection rate,and fracture width on the pressure response and plugging efficiency of the fracture.The results revealed that the temporary plugging process of artificial fractures in hot dry rock can be divided into four main stages:the initial stage of temporary plugging,the bridging stage of the particles,the plugging formation stage,and the high-pressure dense plugging stage.As the temperature increases,the distribution distance of the temporary plugging agent,the number of pressure fluctuations,and the time required for crack plugging increases.Particularly,when the temperature increases by 100℃,the complete plugging time increases by 90.7%.
基金financially supported by the National Natural Science Foundation of China(52074023)the Beijing Natural Science Foundation(2222062)+1 种基金the National Key R&D Program of China(2018YFB0905600)the Interdisciplinary Research Project for Young Teachers of USTB(Fundamental Research Funds for the Central Universities)(FRF-IDRY-21-023)。
文摘Antimony(Sb) is an attractive cathode for liquid metal batteries(LMBs) because of its high theoretical voltage and low cost.The main obstacles associated with the Sb-based cathodes are unsatisfactory energy density and poor rate-capability.Herein,we propose a novel Sb_(64)Cu_(36)cathode that effectively tackles these issues.The Sb_(64)Cu_(36)(melting point:525℃) cathode presents a novel lithiation mechanism involving sequentially the generation of Li_(2)CuSb,the formation of Li_(3)Sb,and the conversion reaction of Li_(2)CuSb to Li_(3)Sb and Cu.The generated intermetallic compounds show a unique microstructure of the upper floated Li_(2)CuSb layer and the below cross-linked structure with interpenetrated Li_(2)CuSb and Li_(3)Sb phases.Compared with Li_(3)Sb,the lower Li migration energy barrier(0.188 eV) of Li_(2)CuSb significantly facilitates the lithium diffusion across the intermediate compounds and accelerates the reaction kinetics.Consequently,the Li‖Sb_(64)Cu_(36)cell delivers a more excellent electrochemical performance(energy density:353 W h kg^(-1)at 0.4 A cm^(-2);rate capability:0.59 V at 2.0 A cm^(-2)),and a much lower energy storage cost of only 38.45 $ kW h^(-1)than other previously reported Sb-based LMBs.This work provides a novel cathode design concept for the development of high-performance LMBs in applications for large-scale energy storage.