NAMar 4, 2016
Energy-conserving method for Stochastic Maxwell Equations with Multiplicative NoiseJialin Hong, Lihai Ji, Liying Zhang et al.
In this paper, it is shown that three-dimensional stochastic Maxwell equations with multiplicative noise are stochastic Hamiltonian partial differential equations possessing a geometric structure (i.e. stochastic mutli-symplectic conservation law), and the energy of system is a conservative quantity almost surely. We propose a stochastic multi-symplectic energy-conserving method for the equations by using the wavelet collocation method in space and stochastic symplectic method in time. Numerical experiments are performed to verify the excellent abilities of the proposed method in providing accurate solution and preserving energy. The mean square convergence result of the method in temporal direction is tested numerically, and numerical comparisons with finite difference method are also investigated.
NAJan 14, 2018
Optimal Error Estimates of Conservative Local Discontinuous Galerkin Method for Nonlinear Schrödinger EquationJialin Hong, Lihai Ji, Zhihui Liu
In this paper, we propose a conservative local discontinuous Galerkin method for one-dimensional nonlinear Schrödinger equation. By using special upwind-biased numerical fluxes, we establish the optimal rate of convergence $\mathcal O(h^{k+1})$, with polynomial of degree $k$ and grid size $h$. Meanwhile, we show that this method preserves the charge conservation law and thus we call it a conservative local discontinuous Galerkin method. Numerical experiments verify our theoretical result.
NAJun 6, 2018
Runge-Kutta semidiscretizations for stochastic Maxwell equations with additive noiseChuchu Chen, Jialin Hong, Lihai Ji
The paper concerns semidiscretizations in time of stochastic Maxwell equations driven by additive noise. We show that the equations admit physical properties and mathematical structures, including regularity, energy and divergence evolution laws, and stochastic symplecticity, etc. In order to inherit the intrinsic properties of the original system, we introduce a general class of stochastic Runge-Kutta methods, and deduce the condition of symplecticity-preserving. By utilizing a priori estimates on numerical approximations and semigroup approach, we show that the methods, which are algebraically stable and coercive, are well-posed and convergent with order one in mean-square sense, which answers an open problem in [Chen and Hong, SIAM J. Numer. Anal., 2016] for stochastic Maxwell equations driven by additive noise.
SGMar 29, 2018
Stochastic conformal multi-symplectic method for damped stochastic nonlinear Schrodinger equationChuchu Chen, Jialin Hong, Lihai Ji
In this paper, we propose a stochastic conformal multi-symplectic method for a class of damped stochastic Hamiltonian partial differential equations in order to inherit the intrinsic properties, and apply the numerical method to solve a kind of damped stochastic nonlinear Schrodinger equation with multiplicative noise. It is shown that the stochastic conformal multi-symplectic method preserves the discrete stochastic conformal multi-symplectic conservation law, the discrete charge exponential dissipation law almost surely, and we also deduce the recurrence relation of the discrete global energy. Numerical experiments are preformed to verify the good performance of the proposed stochastic conformal multi-symplectic method, compared with a Crank-Nicolson type method. Finally, we present the mean square convergence result of the proposed numerical method in temporal direction numerically.
NANov 9, 2017
Stochastic K-symplectic integrators for stochastic non-canonical Hamiltonian systems and applications to the Lotka--Volterra modelJialin Hong, Lihai Ji, Xu Wang et al.
We give a theoretical framework of stochastic non-canonical Hamiltonian systems as well as their modified symplectic structure which is named stochastic K-symplectic structure. The framework can be applied to the study of the Lotka--Volterra model perturbed by external noises. In terms of internal properties of the stochastic Lotka--Volterra model, we propose different kinds of stochastic K-symplectic integrators which could inherit the positivity of the solution. The K-symplectic conditions are also obtained to ensure that the proposed schemes admit the same geometric structure as the original system. Besides, the first-order condition of the proposed schemes in $L^1(Ω)$ sense are given based on the uniform boundedness of both the exact solution and the numerical one. Several numerical examples are illustrated to verify above properties of proposed schemes compared with non-K-symplectic ones.
NANov 27, 2016
Convergence in probability of an ergodic and conformal multi-symplectic numerical scheme for a damped stochastic NLS equationJialin Hong, Lihai Ji, Xu Wang
In this paper, we investigate the convergence order in probability of a novel ergodic numerical scheme for damped stochastic nonlinear Schrödinger equation with an additive noise. Theoretical analysis shows that our scheme is of order one in probability under appropriate assumptions for the initial value and noise. Meanwhile, we show that our scheme possesses the unique ergodicity and preserves the discrete conformal multi-symplectic conservation law. Numerical experiments are given to show the longtime behavior of the discrete charge and the time average of the numerical solution, and to test the convergence order, which verify our theoretical results.
NAFeb 28, 2018
Mean-square convergence of a semi-discrete scheme for stochastic nonlinear Maxwell equationsChuchu Chen, Jialin Hong, Lihai Ji
In this paper, we propose a semi-implicit Euler scheme to discretize the stochastic nonlinear Maxwell equations with multiplicative Ito noise, which is implicit in the drift term and explicit in the diffusion term of the equations, in order to suited to Ito product. Uniform bounds with high regularities of solutions for both the continuous and the discrete problems are obtained, which are crucial properties to derive the mean-square convergence with certain order. Allowing sufficient spatial regularity and utilizing the energy estimate technique, the convergence order 1/2 in mean-square sense is obtained.
NAMay 23, 2017
Parareal algorithms applied to stochastic differential equations with conserved quantitiesLiying Zhang, Weien Zhou, Lihai ji
In this papers, we couple the parareal algorithm with projection methods of the trajectory on a specific manifold, defined by the preservation of some conserved quantities of the differential equations. First, projection methods are introduced as the coarse and fine propagators. Second, we also apply the projection methods for systems with conserved quantities in the correction step of original parareal algorithm. Finally, three numerical experiments are performed by different kinds of algorithms to show the property of convergence in iteration, and preservation in conserved quantities of model systems.
NAMar 24, 2015
Mean-square Convergence of a Symplectic Local Discontinuous Galerkin Method Applied to Stochastic Linear Schroedinger EquationChuchu Chen, Jialin Hong, Lihai Ji
In this paper, we investigate the mean-square convergence of a novel symplectic local discontinuous Galerkin method in L^2-norm for stochastic linear Schroedinger equation with multiplicative noise. It is shown that the mean-square error is bounded not only by the temporal and spatial step-sizes, but also by their ratio. The mean-square convergence rate with respect to the computational cost is derived under appropriate assumptions for initial data and noise. Meanwhile, we show that the method preserves the discrete charge conservation law which implies an L^2-stability
NAOct 27, 2014
Energy evolution of multi-symplectic methods for Maxwell equations with perfectly matched layer boundaryJialin Hong, Lihai Ji
In this paper, we consider the energy evolution of multi-symplectic methods for three-dimensional (3D) Maxwell equations with perfectly matched layer boundary, and present the energy evolution laws of Maxwell equations under the discretization of multi-symplectic Yee method and general multi-symplectic Runge-Kutta methods.