The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varyin...The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varying dynamics,making the exact solutions challenging to obtain.As a result,numerical integration methods are typically employed.However,conventional methods often suffer from low computational efficiency.To address this,this paper explores the application of the parameter freezing precise exponential integrator to vehicle-road coupling models.The model accounts for road roughness irregularities,incorporating all terms unrelated to the linear part into the algorithm's inhomogeneous vector.The general construction process of the algorithm is detailed.The validity of numerical results is verified through approximate analytical solutions(AASs),and the advantages of this method over traditional numerical integration methods are demonstrated.Multiple parameter freezing precise exponential integrator schemes are constructed based on the Runge-Kutta framework,with the fourth-order four-stage scheme identified as the optimal one.The study indicates that this method can quickly and accurately capture the dynamic system's vibration response,offering a new,efficient approach for numerical studies of high-dimensional vehicle-road coupling systems.展开更多
In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential ...In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential integrator.The scheme is comprised of two explicit stages.One is the exponential integrator type stage,and the second is the Runge-Kutta type stage.The spatial-dependent terms are discretized using the compact technique.The compact scheme can achieve fourth or sixth-order spatial accuracy,while the proposed scheme attains second-order temporal accuracy.Also,a mathematical model for the electro-osmotic flow of Carreau fluid over the stationary sheet is presented with heat and mass transfer effects.The governing equations are transformed into dimensionless partial differential equations and solved by the proposed scheme.Simulation results reveal that increasing the Helmholtz-Smoluchowski velocityUHS by 400%leads to a 60%-75%rise in peak flowvelocity,while the electro-osmotic parameter me enhances near-wall acceleration.Conversely,velocity decreases significantly with higher Weissenberg numbers,indicating the Carreau fluid’s elastic resistance and increased magnetic field strength due to improved Lorentz forces.Temperature rises with the thermal conductivity parameter 2,while higher reaction ratesγdiminish concentration and local Sherwood number values.The simulation findings show the scheme’s correctness and efficacy in capturing the complicated interactions in non-Newtonian electro-osmotic transport by revealing the notable impact of electrokinetic factors on flowbehaviour.Theproposedmodel is particularly relevant for BiologicalMicro-Electro-Mechanical Systems(BioMEMS)applications,where precise control of electro-thermal transport in non-Newtonian fluids is critical for lab-on-a-chip diagnostics,drug delivery,and micro-scale thermal management.展开更多
In this paper, by Laplace transform version of the Trotter-Kato approximation theorem and the integrated C-semigroup introduced by Myadera, the authors obtained some Trotter-Kato approximation theorems on exponentiall...In this paper, by Laplace transform version of the Trotter-Kato approximation theorem and the integrated C-semigroup introduced by Myadera, the authors obtained some Trotter-Kato approximation theorems on exponentially bounded C-semigroups, where the range of C (and so the domain of the generator) may not be dense. The authors deduced the corresponding results on exponentially bounded integrated semigroups with nondensely generators. The results of this paper extended and perfected the results given by Lizama, Park and Zheng.展开更多
An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
Exponential integral for real arguments is evaluated by employing a fast-converging power series originally developed for the resolution of Grandi’s paradox. Laguerre’s historic solution is first recapitulated and t...Exponential integral for real arguments is evaluated by employing a fast-converging power series originally developed for the resolution of Grandi’s paradox. Laguerre’s historic solution is first recapitulated and then the new solution method is described in detail. Numerical results obtained from the present series solution are compared with the tabulated values correct to nine decimal places. Finally, comments are made for the further use of the present approach for integrals involving definite functions in denominator.展开更多
In this work,we propose a low-regularity Fourier integrator with almost mass conservation to solve the Davey-StewartsonⅡsystem(hyperbolic-elliptic case).Arbitrary order mass convergence could be achieved by the suita...In this work,we propose a low-regularity Fourier integrator with almost mass conservation to solve the Davey-StewartsonⅡsystem(hyperbolic-elliptic case).Arbitrary order mass convergence could be achieved by the suitable addition of correction terms,while keeping the first order accuracy in H~γ×H^(γ+1)for initial data in H^(γ+1)×H^(γ+1)withγ>1.The main theorem is that,up to some fixed time T,there exist constantsτ_(0)and C depending only on T and‖u‖_(L^(∞)((0,T);H^(γ+1)))such that,for any 0<τ≤τ_(0),we have that‖u(t_(n),·)-u^(n)‖H_γ≤C_(τ),‖v(t_(n),·)-v^(n)‖_(Hγ+1)≤C_(τ),where u^(n)and v^(n)denote the numerical solutions at t_(n)=nτ.Moreover,the mass of the numerical solution M(u^(n))satisfies that|M(u^(n))-M(u_0)|≤Cτ~5.展开更多
Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonl...Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonlinearity controlled by a small valueε∈(0,1],an exponential wave integrator Fourier pseudo-spectral(EWIFP)discretization has been developed(Guo et al.,2021)and proved to be uniformly accurate aboutεup to the time atΟ(1/ε^(2))However,the EWIFP method is not time symmetric and can not preserve the discrete energy.As we know,the time symmetry and energy-preservation are the important structural features of the true solution and we hope that this structure can be inherited along the numerical solution.In this work,we propose a time symmetric and energy-preserving exponential wave integrator Fourier pseudo-spectral(SEPEWIFP)method for the NLSW with periodic boundary conditions.Through rigorous error analysis,we establish uniform error bounds of the numerical solution atΟ(h^(mo)+ε^(2-βτ2))up to the time atΟ(1/ε^(β))forβ∈[0,2]where h andτare the mesh size and time step,respectively,and m0 depends on the regularity conditions.The tools for error analysis mainly include cut-off technique and the standard energy method.We also extend the results on error bounds,energy-preservation and time symmetry to the oscillatory NLSW with wavelength atΟ(1/ε^(2))in time which is equivalent to the NLSW with weak nonlinearity.Numerical experiments confirm that the theoretical results in this paper are correct.Our method is novel because that to the best of our knowledge there has not been any energy-preserving exponential wave integrator method for the NLSW.展开更多
In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materi...In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.展开更多
This study develops a high-order computational scheme for analyzing unsteady tangent hyperbolic fluid flow with variable thermal conductivity,thermal radiation,and coupled heat andmass transfer effects.Amodified twost...This study develops a high-order computational scheme for analyzing unsteady tangent hyperbolic fluid flow with variable thermal conductivity,thermal radiation,and coupled heat andmass transfer effects.Amodified twostage Exponential Time Integrator is introduced for temporal discretization,providing second-order accuracy in time.A compact finite difference method is employed for spatial discretization,yielding sixth-order accuracy at most grid points.The proposed framework ensures numerical stability and convergence when solving stiff,nonlinear parabolic systems arising in fluid flow and heat transfer problems.The novelty of the work lies in combining exponential integrator schemes with compact high-order spatial discretization,enabling accurate and efficient simulations of tangent hyperbolic fluids under complex boundary conditions,such as oscillatory plates and varying thermal conductivity.This approach addresses limitations of classical Euler,Runge–Kutta,and spectral methods by significantly reducing numerical errors(up to 45%)and computational cost.Comprehensive parametric studies demonstrate how viscous dissipation,chemical reactions,the Weissenberg number,and the Hartmann number influence flow behaviour,heat transfer,and mass transfer.Notably,heat transfer rates increase by 18.6%with stronger viscous dissipation,while mass transfer rates rise by 21.3%with more intense chemical reactions.The real-world relevance of the study is underscored by its direct applications in polymer processing,heat exchanger design,radiative thermal management in aerospace,and biofluid transport in biomedical systems.The proposed scheme thus provides a robust numerical framework that not only advances the mathematical modelling of non-Newtonian fluid flows but also offers practical insights for engineering systems involving tangent hyperbolic fluids.展开更多
We study a class of stochastic semilinear damped wave equations driven by additive Wiener noise.Owing to the damping term,under appropriate conditions on the nonlinearity,the solution admits a unique invariant distrib...We study a class of stochastic semilinear damped wave equations driven by additive Wiener noise.Owing to the damping term,under appropriate conditions on the nonlinearity,the solution admits a unique invariant distribution.We apply semi-discrete and fully-discrete methods in order to approximate this invariant distribution,using a spectral Galerkin method and an exponential Euler integrator for spatial and temporal discretization respectively.We prove that the considered numerical schemes also admit unique invariant distributions,and we prove error estimates between the approximate and exact invariant distributions,with identification of the orders of convergence.To the best of our knowledge this is the first result in the literature concerning numerical approximation of invariant distributions for stochastic damped wave equations.展开更多
In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation method...In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.展开更多
We study an explicit exponential scheme for the time discretisation of stochastic SchrS- dinger Equations Driven by additive or Multiplicative It6 Noise. The numerical scheme is shown to converge with strong order 1 i...We study an explicit exponential scheme for the time discretisation of stochastic SchrS- dinger Equations Driven by additive or Multiplicative It6 Noise. The numerical scheme is shown to converge with strong order 1 if the noise is additive and with strong order 1/2 for multiplicative noise. In addition, if the noise is additive, we show that the exact solutions of the linear stochastic Sehr6dinger equations satisfy trace formulas for the expected mass, energy, and momentum （i. e., linear drifts in these quantities）. Furthermore, we inspect the behaviour of the numerical solutions with respect to these trace formulas. Several numerical simulations are presented and confirm our theoretical results.展开更多
An efficient and accurate exponential wave integrator Fourier pseudospectral （EWI-FP） method is proposed and analyzed for solving the symmetric regularized-long-wave （SRLW） equation, which is used for modeling the...An efficient and accurate exponential wave integrator Fourier pseudospectral （EWI-FP） method is proposed and analyzed for solving the symmetric regularized-long-wave （SRLW） equation, which is used for modeling the weakly nonlinear ion acoustic and space-charge waves. The numerical method here is based on a Gautschi-type exponential wave integrator for temporal approximation and the Fourier pseudospectral method for spatial discretization. The scheme is fully explicit and efficient due to the fast Fourier transform. Numerical analysis of the proposed EWI-FP method is carried out and rigorous error estimates are established without CFL-type condition by means of the mathematical induction. The error bound shows that EWI-FP has second order accuracy in time and spectral accuracy in space. Numerical results are reported to confirm the theoretical studies and indicate that the error bound here is optimal.展开更多
In this paper,we consider the derivatives of intersection local time for two independent d-dimensional symmetricα-stable processes X^(α) and X^(α)with respective indices α and α.We first study the sufficient cond...In this paper,we consider the derivatives of intersection local time for two independent d-dimensional symmetricα-stable processes X^(α) and X^(α)with respective indices α and α.We first study the sufficient condition for the existence of the derivatives,which makes us obtain the exponential integrability and H?lder continuity.Then we show that this condition is also necessary for the existence of derivatives of intersection local time at the origin.Moreover,we also study the power variation of the derivatives.展开更多
We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0...We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O（ ε^2） and O（1） in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference （CNFD） method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain ＇correct＇ numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability： τ- = O（ε3） and h = O（√ε）. Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O（ε2） and h = O（1） when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.展开更多
In this paper,we consider the numerics of the dispersion-managed Kortewegde Vries(DM-KdV)equation for describingwave propagations in inhomogeneous media.The DM-KdV equation contains a variable dispersion map with disc...In this paper,we consider the numerics of the dispersion-managed Kortewegde Vries(DM-KdV)equation for describingwave propagations in inhomogeneous media.The DM-KdV equation contains a variable dispersion map with discontinuity,which makes the solution non-smooth in time.We formally analyze the convergence order reduction problems of some popular numerical methods including finite difference and time-splitting for solving the DM-KdV equation,where a necessary constraint on the time step has been identified.Then,two exponential-type dispersionmap integrators up to second order accuracy are derived,which are efficiently incorporatedwith the Fourier pseudospectral discretization in space,and they can converge regardless the discontinuity and the step size.Numerical comparisons show the advantage of the proposed methods with the application to solitary wave dynamics and extension to the fast&strong dispersion-management regime.展开更多
In a previous paper, some particular multistep cosine methods were constructed which proved to be very efficient because of being able to integrate in a stable and explicit way linearly stiff problems of second-order ...In a previous paper, some particular multistep cosine methods were constructed which proved to be very efficient because of being able to integrate in a stable and explicit way linearly stiff problems of second-order in time. In the present paper, the conditions which guarantee stability for general methods of this type are given, as well as a thorough study of resonances and filtering for symmetric ones （which, in another paper, have been proved to behave very advantageously with respect to conservation of invariants in Hamiltonian wave equations）. What is given here is a systematic way to analyse and treat any of the methods of this type in the mentioned aspects.展开更多
A unified solution framework is proposed for efficiently solving conjugate fluid and solid heat transfer problems.The unified solution is solely governed by the compressible Navier-Stokes(N-S)equations in both fluid a...A unified solution framework is proposed for efficiently solving conjugate fluid and solid heat transfer problems.The unified solution is solely governed by the compressible Navier-Stokes(N-S)equations in both fluid and solid domains.Such method not only provides the computational capability for solid heat transfer simulations with existing successful N-S flow solvers,but also can relax time-stepping restrictions often imposed by the interface conditions for conjugate fluid and solid heat transfer.This paper serves as Part I of the proposed unified solution framework and addresses the handling of solid heat conduction with the nondimensional N-S equations.Specially,a parallel,adaptive high-order discontinuous Galerkin unified solver has been developed and applied to solve solid heat transfer problems under various boundary conditions.展开更多
基金Supported by the National Natural Science Foundation of China(No.U22A20246)the Key Project of Natural Science Foundation of Hebei Province of China(Basic Research Base Project)(No.A2023210064)the Science and Technology Program of Hebei Province of China(Nos.246Z1904G and 225676162GH)。
文摘The vehicle-road coupling dynamics problem is a prominent issue in transportation,drawing significant attention in recent years.These dynamic equations are characterized by high-dimensionality,coupling,and time-varying dynamics,making the exact solutions challenging to obtain.As a result,numerical integration methods are typically employed.However,conventional methods often suffer from low computational efficiency.To address this,this paper explores the application of the parameter freezing precise exponential integrator to vehicle-road coupling models.The model accounts for road roughness irregularities,incorporating all terms unrelated to the linear part into the algorithm's inhomogeneous vector.The general construction process of the algorithm is detailed.The validity of numerical results is verified through approximate analytical solutions(AASs),and the advantages of this method over traditional numerical integration methods are demonstrated.Multiple parameter freezing precise exponential integrator schemes are constructed based on the Runge-Kutta framework,with the fourth-order four-stage scheme identified as the optimal one.The study indicates that this method can quickly and accurately capture the dynamic system's vibration response,offering a new,efficient approach for numerical studies of high-dimensional vehicle-road coupling systems.
基金supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University(IMSIU)(grant number IMSIU-DDRSP2503).
文摘In this study,we proposed a numerical technique for solving time-dependent partial differential equations that arise in the electro-osmotic flowofCarreau fluid across a stationary plate based on amodified exponential integrator.The scheme is comprised of two explicit stages.One is the exponential integrator type stage,and the second is the Runge-Kutta type stage.The spatial-dependent terms are discretized using the compact technique.The compact scheme can achieve fourth or sixth-order spatial accuracy,while the proposed scheme attains second-order temporal accuracy.Also,a mathematical model for the electro-osmotic flow of Carreau fluid over the stationary sheet is presented with heat and mass transfer effects.The governing equations are transformed into dimensionless partial differential equations and solved by the proposed scheme.Simulation results reveal that increasing the Helmholtz-Smoluchowski velocityUHS by 400%leads to a 60%-75%rise in peak flowvelocity,while the electro-osmotic parameter me enhances near-wall acceleration.Conversely,velocity decreases significantly with higher Weissenberg numbers,indicating the Carreau fluid’s elastic resistance and increased magnetic field strength due to improved Lorentz forces.Temperature rises with the thermal conductivity parameter 2,while higher reaction ratesγdiminish concentration and local Sherwood number values.The simulation findings show the scheme’s correctness and efficacy in capturing the complicated interactions in non-Newtonian electro-osmotic transport by revealing the notable impact of electrokinetic factors on flowbehaviour.Theproposedmodel is particularly relevant for BiologicalMicro-Electro-Mechanical Systems(BioMEMS)applications,where precise control of electro-thermal transport in non-Newtonian fluids is critical for lab-on-a-chip diagnostics,drug delivery,and micro-scale thermal management.
基金This project is supported by the National Science Foundation of China.
文摘In this paper, by Laplace transform version of the Trotter-Kato approximation theorem and the integrated C-semigroup introduced by Myadera, the authors obtained some Trotter-Kato approximation theorems on exponentially bounded C-semigroups, where the range of C (and so the domain of the generator) may not be dense. The authors deduced the corresponding results on exponentially bounded integrated semigroups with nondensely generators. The results of this paper extended and perfected the results given by Lizama, Park and Zheng.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘Exponential integral for real arguments is evaluated by employing a fast-converging power series originally developed for the resolution of Grandi’s paradox. Laguerre’s historic solution is first recapitulated and then the new solution method is described in detail. Numerical results obtained from the present series solution are compared with the tabulated values correct to nine decimal places. Finally, comments are made for the further use of the present approach for integrals involving definite functions in denominator.
基金supported by the NSFC(11901120)supported by the NSFC(12171356)the Science and Technology Program of Guangzhou,China(2024A04J4027)。
文摘In this work,we propose a low-regularity Fourier integrator with almost mass conservation to solve the Davey-StewartsonⅡsystem(hyperbolic-elliptic case).Arbitrary order mass convergence could be achieved by the suitable addition of correction terms,while keeping the first order accuracy in H~γ×H^(γ+1)for initial data in H^(γ+1)×H^(γ+1)withγ>1.The main theorem is that,up to some fixed time T,there exist constantsτ_(0)and C depending only on T and‖u‖_(L^(∞)((0,T);H^(γ+1)))such that,for any 0<τ≤τ_(0),we have that‖u(t_(n),·)-u^(n)‖H_γ≤C_(τ),‖v(t_(n),·)-v^(n)‖_(Hγ+1)≤C_(τ),where u^(n)and v^(n)denote the numerical solutions at t_(n)=nτ.Moreover,the mass of the numerical solution M(u^(n))satisfies that|M(u^(n))-M(u_0)|≤Cτ~5.
基金supported in part by the Natural Science Foundation of Hebei Province(Grant No.A2021205036).
文摘Abstract Recently,the numerical methods for long-time dynamics of PDEs with weak nonlinearity have received more and more attention.For the nonlinear Schrödinger equation(NLS)with wave operator(NLSW)and weak nonlinearity controlled by a small valueε∈(0,1],an exponential wave integrator Fourier pseudo-spectral(EWIFP)discretization has been developed(Guo et al.,2021)and proved to be uniformly accurate aboutεup to the time atΟ(1/ε^(2))However,the EWIFP method is not time symmetric and can not preserve the discrete energy.As we know,the time symmetry and energy-preservation are the important structural features of the true solution and we hope that this structure can be inherited along the numerical solution.In this work,we propose a time symmetric and energy-preserving exponential wave integrator Fourier pseudo-spectral(SEPEWIFP)method for the NLSW with periodic boundary conditions.Through rigorous error analysis,we establish uniform error bounds of the numerical solution atΟ(h^(mo)+ε^(2-βτ2))up to the time atΟ(1/ε^(β))forβ∈[0,2]where h andτare the mesh size and time step,respectively,and m0 depends on the regularity conditions.The tools for error analysis mainly include cut-off technique and the standard energy method.We also extend the results on error bounds,energy-preservation and time symmetry to the oscillatory NLSW with wavelength atΟ(1/ε^(2))in time which is equivalent to the NLSW with weak nonlinearity.Numerical experiments confirm that the theoretical results in this paper are correct.Our method is novel because that to the best of our knowledge there has not been any energy-preserving exponential wave integrator method for the NLSW.
基金Jie Du is supported by the National Natural Science Foundation of China under Grant Number NSFC 11801302Tsinghua University Initiative Scientific Research Program+1 种基金Eric Chung is supported by Hong Kong RGC General Research Fund(Projects 14304217 and 14302018)The third author is supported by the NSF grant DMS-1818467.
文摘In this paper, we study the classical Allen-Cahn equations and investigate the maximum- principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen- Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to dem-onstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use rela-tively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.
基金supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University(IMSIU)(grant number IMSIU-DDRSP2503).
文摘This study develops a high-order computational scheme for analyzing unsteady tangent hyperbolic fluid flow with variable thermal conductivity,thermal radiation,and coupled heat andmass transfer effects.Amodified twostage Exponential Time Integrator is introduced for temporal discretization,providing second-order accuracy in time.A compact finite difference method is employed for spatial discretization,yielding sixth-order accuracy at most grid points.The proposed framework ensures numerical stability and convergence when solving stiff,nonlinear parabolic systems arising in fluid flow and heat transfer problems.The novelty of the work lies in combining exponential integrator schemes with compact high-order spatial discretization,enabling accurate and efficient simulations of tangent hyperbolic fluids under complex boundary conditions,such as oscillatory plates and varying thermal conductivity.This approach addresses limitations of classical Euler,Runge–Kutta,and spectral methods by significantly reducing numerical errors(up to 45%)and computational cost.Comprehensive parametric studies demonstrate how viscous dissipation,chemical reactions,the Weissenberg number,and the Hartmann number influence flow behaviour,heat transfer,and mass transfer.Notably,heat transfer rates increase by 18.6%with stronger viscous dissipation,while mass transfer rates rise by 21.3%with more intense chemical reactions.The real-world relevance of the study is underscored by its direct applications in polymer processing,heat exchanger design,radiative thermal management in aerospace,and biofluid transport in biomedical systems.The proposed scheme thus provides a robust numerical framework that not only advances the mathematical modelling of non-Newtonian fluid flows but also offers practical insights for engineering systems involving tangent hyperbolic fluids.
基金supported by the projects ADA(Grant No.ANR-19-CE40-0019-02)and SIMALIN(Grant No.ANR-19-CE40-0016)operated by the French National Research Agencysupported by the NSF of China(Grant Nos.11971488,12371417)supported by the China Scholarship Council(Grant No.202206370085)。
文摘We study a class of stochastic semilinear damped wave equations driven by additive Wiener noise.Owing to the damping term,under appropriate conditions on the nonlinearity,the solution admits a unique invariant distribution.We apply semi-discrete and fully-discrete methods in order to approximate this invariant distribution,using a spectral Galerkin method and an exponential Euler integrator for spatial and temporal discretization respectively.We prove that the considered numerical schemes also admit unique invariant distributions,and we prove error estimates between the approximate and exact invariant distributions,with identification of the orders of convergence.To the best of our knowledge this is the first result in the literature concerning numerical approximation of invariant distributions for stochastic damped wave equations.
文摘In this paper, a novel class of exponential Fourier collocation methods （EFCMs） is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation meth- ods. We discuss in detail the connections of EFCMs with trigonometric Fourier colloca- tion methods （TFCMs）, the well-known Hamiltonian Boundary Value Methods （HBVMs）, Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an es- sential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first- order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM（2,2） as an illustrative example. The numerical experiments using EFCM（2,2） are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM（2,2）.
文摘We study an explicit exponential scheme for the time discretisation of stochastic SchrS- dinger Equations Driven by additive or Multiplicative It6 Noise. The numerical scheme is shown to converge with strong order 1 if the noise is additive and with strong order 1/2 for multiplicative noise. In addition, if the noise is additive, we show that the exact solutions of the linear stochastic Sehr6dinger equations satisfy trace formulas for the expected mass, energy, and momentum （i. e., linear drifts in these quantities）. Furthermore, we inspect the behaviour of the numerical solutions with respect to these trace formulas. Several numerical simulations are presented and confirm our theoretical results.
文摘An efficient and accurate exponential wave integrator Fourier pseudospectral （EWI-FP） method is proposed and analyzed for solving the symmetric regularized-long-wave （SRLW） equation, which is used for modeling the weakly nonlinear ion acoustic and space-charge waves. The numerical method here is based on a Gautschi-type exponential wave integrator for temporal approximation and the Fourier pseudospectral method for spatial discretization. The scheme is fully explicit and efficient due to the fast Fourier transform. Numerical analysis of the proposed EWI-FP method is carried out and rigorous error estimates are established without CFL-type condition by means of the mathematical induction. The error bound shows that EWI-FP has second order accuracy in time and spectral accuracy in space. Numerical results are reported to confirm the theoretical studies and indicate that the error bound here is optimal.
基金Supported by National Natural Science Foundation of China(Grant Nos.12071003,12201294)Natural Science Foundation of Jiangsu Province,China(Grant No.BK20220865)。
文摘In this paper,we consider the derivatives of intersection local time for two independent d-dimensional symmetricα-stable processes X^(α) and X^(α)with respective indices α and α.We first study the sufficient condition for the existence of the derivatives,which makes us obtain the exponential integrability and H?lder continuity.Then we show that this condition is also necessary for the existence of derivatives of intersection local time at the origin.Moreover,we also study the power variation of the derivatives.
基金supported by the Ministry of Education of Singapore(Grant No.R146-000-196-112)National Natural Science Foundation of China(Grant No.91430103)
文摘We present several numerical methods and establish their error estimates for the discretization of the nonlinear Dirac equation （NLDE） in the nonrelativistic limit regime, involving a small dimensionless parameter 0 〈 ε〈〈1 which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e., there are propagating waves with wavelength O（ ε^2） and O（1） in time and space, respectively. We begin with the conservative Crank-Nicolson finite difference （CNFD） method and establish rigorously its error estimate which depends explicitly on the mesh size h and time step τ- as well as the small parameter 0 〈 ε≤1 Based on the error bound, in order to obtain ＇correct＇ numerical solutions in the nonrelativistic limit regime, i.e., 0 〈 ε≤1 , the CNFD method requests the ε-scalability： τ- = O（ε3） and h = O（√ε）. Then we propose and analyze two numerical methods for the discretization of NLDE by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time- splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their ε-scalability is improved to τ = O（ε2） and h = O（1） when 0 〈 ε 〈〈 1. Extensive numerical results are reported to confirm our error estimates.
基金supported by the National Key Research and Development Program of China(No.2020YFA0714200)the Natural Science Foundation of Hubei Province No.2019CFA007,the NSFC 11901440。
文摘In this paper,we consider the numerics of the dispersion-managed Kortewegde Vries(DM-KdV)equation for describingwave propagations in inhomogeneous media.The DM-KdV equation contains a variable dispersion map with discontinuity,which makes the solution non-smooth in time.We formally analyze the convergence order reduction problems of some popular numerical methods including finite difference and time-splitting for solving the DM-KdV equation,where a necessary constraint on the time step has been identified.Then,two exponential-type dispersionmap integrators up to second order accuracy are derived,which are efficiently incorporatedwith the Fourier pseudospectral discretization in space,and they can converge regardless the discontinuity and the step size.Numerical comparisons show the advantage of the proposed methods with the application to solitary wave dynamics and extension to the fast&strong dispersion-management regime.
文摘In a previous paper, some particular multistep cosine methods were constructed which proved to be very efficient because of being able to integrate in a stable and explicit way linearly stiff problems of second-order in time. In the present paper, the conditions which guarantee stability for general methods of this type are given, as well as a thorough study of resonances and filtering for symmetric ones （which, in another paper, have been proved to behave very advantageously with respect to conservation of invariants in Hamiltonian wave equations）. What is given here is a systematic way to analyse and treat any of the methods of this type in the mentioned aspects.
基金S.Li acknowledges the support from the National Natural Science Foundation of China(NSFC)under grant No.U1930402L.Ju’s work is partially supported by U.S.National Science Foundation DMS-2109633.
文摘A unified solution framework is proposed for efficiently solving conjugate fluid and solid heat transfer problems.The unified solution is solely governed by the compressible Navier-Stokes(N-S)equations in both fluid and solid domains.Such method not only provides the computational capability for solid heat transfer simulations with existing successful N-S flow solvers,but also can relax time-stepping restrictions often imposed by the interface conditions for conjugate fluid and solid heat transfer.This paper serves as Part I of the proposed unified solution framework and addresses the handling of solid heat conduction with the nondimensional N-S equations.Specially,a parallel,adaptive high-order discontinuous Galerkin unified solver has been developed and applied to solve solid heat transfer problems under various boundary conditions.
