BZPEERReading, trust and collaboration for research
Menu

Discover

SearchFind papers, people, topics and opportunities from one place.GraphOpen maps around papers, topics and saved views.

Workspaces

HomeOpen the ranked research feed.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Source author record

Jialin Hong

Jialin Hong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.PR

Catalog footprint

What is connected

30works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

Central limit theorem for full discretization of parabolic SPDE

In order to characterize the fluctuation between the ergodic limit and the time-averaging estimator of a full discretization in a quantitative way, we establish a central limit theorem for the full discretization of the parabolic stochastic partial differential equation. The theorem shows that the normalized time-averaging estimator converges to a normal distribution with the variance being the same as that of the continuous case, where the scale used for the normalization corresponds to the temporal strong convergence order of the considered full discretization. A key ingredient in the proof is to extract an appropriate martingale difference series sum from the normalized time-averaging estimator so that the convergence to the normal distribution of such a sum and the convergence to zero in probability of the remainder are well balanced. The main novelty of our method to balance the convergence lies in proposing an appropriately modified Poisson equation so as to possess the space-independent regularity estimates. As a byproduct, the full discretization is shown to fulfill the weak law of large numbers, namely, the time-averaging estimator converges to the ergodic limit in probability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal Hölder Continuity and Hitting Probabilities for SPDEs with Rough Fractional Noises

We investigate the optimal Hölder continuity and hitting probabilities for systems of stochastic heat equations and stochastic wave equations driven by an additive fractional Brownian sheet with temporal index $1/2$ and spatial index $H\le1/2$. Using stochastic calculus for fractional Brownian motion, we prove that these systems are well-posed and the solutions are Hölder continuous. Furthermore, the optimal Hölder exponents are obtained, which is the first result, as far as we knew, on the optimal Hölder continuity of SHEs and SWEs driven by fractional Brownian sheet that is rough in space. Based on this sharp regularity, we obtain lower and upper bounds of hitting probabilities of the solutions in terms of Bessel--Riesz capacity and Hausdorff measure, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Semi-implicit energy-preserving numerical schemes for stochastic wave equation via SAV approach

In this paper, we propose and analyze semi-implicit numerical schemes for the stochastic wave equation (SWE) with general nonlinearity and multiplicative noise. These numerical schemes, called stochastic scalar auxiliary variable (SAV) schemes, are constructed by transforming the considered SWE into a higher dimensional stochastic system with a stochastic SAV. We prove that they can be solved explicitly and preserve the modified energy evolution law and the regularity structure of the original system. These structure-preserving properties are the keys to overcoming the mutual effect of the noise and nonlinearity. By proving new regularity estimates of the introduced SAV, we establish the strong convergence rate of stochastic SAV schemes and the further fully-discrete schemes with the finite element method in spatial direction. To the best of our knowledge, this is the first result on the construction and strong convergence of semi-implicit energy-preserving schemes for nonlinear SWE.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three kinds of novel multi-symplectic methods for stochastic Hamiltonian partial differential equations

Stochastic Hamiltonian partial differential equations, which possess the multi-symplectic conservation law, are an important and fairly large class of systems. The multi-symplectic methods inheriting the geometric features of stochastic Hamiltonian partial differential equations provide numerical approximations with better numerical stability, and are of vital significance for obtaining correct numerical results. In this paper, we propose three novel multi-symplectic methods for stochastic Hamiltonian partial differential equations based on the local radial basis function collocation method, the splitting technique, and the partitioned Runge-Kutta method. Concrete numerical methods are presented for nonlinear stochastic wave equations, stochastic nonlinear Schrödinger equations, stochastic Korteweg-de Vries equations and stochastic Maxwell equations. We take stochastic wave equations as examples to perform numerical experiments, which indicate the validity of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A new efficient operator splitting method for stochastic Maxwell equations

This paper proposes and analyzes a new operator splitting method for stochastic Maxwell equations driven by additive noise, which not only decomposes the original multi-dimensional system into some local one-dimensional subsystems, but also separates the deterministic and stochastic parts. This method is numerically efficient, and preserves the symplecticity, the multi-symplecticity as well as the growth rate of the averaged energy. A detailed $H^2$-regularity analysis of stochastic Maxwell equations is obtained, which is a crucial prerequisite of the error analysis. Under the regularity assumptions of the initial data and the noise, the convergence order one in mean square sense of the operator splitting method is established.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Energy-preserving fully-discrete schemes for nonlinear stochastic wave equations with multiplicative noise

In this paper, we focus on constructing numerical schemes preserving the averaged energy evolution law for nonlinear stochastic wave equations driven by multiplicative noise. We first apply the compact finite difference method and the interior penalty discontinuous Galerkin finite element method to discretize space variable and present two semi-discrete schemes, respectively. Then we make use of the discrete gradient method and the Padé approximation to propose efficient fully-discrete schemes. These semi-discrete and fully-discrete schemes are proved to preserve the discrete averaged energy evolution law. In particular, we also prove that the proposed fully-discrete schemes exactly inherit the averaged energy evolution law almost surely if the considered model is driven by additive noise. Numerical experiments are given to confirm theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Weak intermittency and second moment bound of a fully discrete scheme for stochastic heat equation

In this paper, we first prove the weak intermittency, and in particular the sharp exponential order $Cλ^4t$ of the second moment of the exact solution of the stochastic heat equation with multiplicative noise and periodic boundary condition, where $λ>0$ denotes the level of the noise. In order to inherit numerically these intrinsic properties of the original equation, we introduce a fully discrete scheme, whose spatial direction is based on the finite difference method and temporal direction is based on the $θ$-scheme. We prove that the second moment of numerical solutions of both spatially semi-discrete and fully discrete schemes grows at least as $\exp\{Cλ^2t\}$ and at most as $\exp\{Cλ^4t\}$ for large $t$ under natural conditions, which implies the weak intermittency of these numerical solutions. Moreover, a renewal approach is applied to show that both of the numerical schemes could preserve the sharp exponential order $Cλ^4t$ of the second moment of the exact solution for large spatial partition number.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absolute continuity and numerical approximation of stochastic Cahn--Hilliard equation with unbounded noise diffusion

In this article, we develop and analyze a full discretization, based on the spatial spectral Galerkin method and the temporal drift implicit Euler scheme, for the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise. By introducing an appropriate decomposition of the numerical approximation, we first use the factorization method to deduce the a priori estimate and regularity estimate of the proposed full discretization. With the help of the variation approach, we then obtain the sharp spatial and temporal convergence rate in negative Sobolev space in mean square sense. Furthermore, the sharp mean square convergence rates in both time and space are derived via the Sobolev interpolation inequality and semigroup theory. To the best of our knowledge, this is the first result on the convergence rate of temporally and fully discrete numerical methods for the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical analysis of a full discretization for stochastic Cahn--Hilliard equation driven by additive noise

In this article, we consider the stochastic Cahn--Hilliard equation driven by space-time white noise. We discretize this equation by using a spatial spectral Galerkin method and a temporal accelerated implicit Euler method. The optimal regularity properties and uniform moment bounds of the exact and numerical solutions are shown. Then we prove that the proposed numerical method is strongly convergent with the sharp convergence rate in a negative Sobolev space. By using an interpolation approach, we deduce the spatial optimal convergence rate and the temporal super-convergence rate of the proposed numerical method in strong convergence sense. To the best of our knowledge, this is the first result on the strong convergence rates of numerical methods for the stochastic Cahn--Hilliard equation driven by space-time white noise. This interpolation approach is also applied to the general noise and high dimension cases, and strong convergence rate results of the proposed scheme are given.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure-preserving numerical methods for stochastic Poisson systems

We propose a class of numerical integration methods for stochastic Poisson systems (SPSs) of arbitrary dimensions. Based on the Darboux-Lie theorem, we transform the SPSs to their canonical form, the generalized stochastic Hamiltonian systems (SHSs), via canonical coordinate transformations found by solving certain PDEs defined by the Poisson brackets of the SPSs. An a-generating function approach with α\in [0,1] is then used to create symplectic discretizations of the SHSs, which are then transformed back by the inverse coordinate transformation to numerical integrators for the SPSs. These integrators are proved to preserve both the Poisson structure and the Casimir functions of the SPSs. Applications to a three-dimensional stochastic rigid body system and a three-dimensional stochastic Lotka-Volterra system show efficiency of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Super-convergence analysis on exponential integrator for stochastic heat equation driven by additive fractional Brownian motion

In this paper, we consider the strong convergence order of the exponential integrator for the stochastic heat equation driven by an additive fractional Brownian motion with Hurst parameter $H\in(\frac12,1)$. By showing the strong order one of accuracy of the exponential integrator under appropriote assumptions, we present the first super-convergence result in temporal direction on full discretizations for stochastic partial differential equations driven by infinite dimensional fractional Brownian motions with Hurst parameter $H\in(\frac12,1)$. The proof is a combination of Malliavin calculus, the $L^p(Ω)$-estimate of the Skorohod integral and the smoothing effect of the Laplacian operator.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Well-posedness and Finite Element Approximations for Elliptic SPDEs with Gaussian Noises

The paper studies the well-posedness and optimal error estimates of spectral finite element approximations for the boundary value problems of semi-linear elliptic SPDEs driven by white or colored Gaussian noises. The noise term is approximated through the spectral projection of the covariance operator, which is not required to be commutative with the Laplacian operator. Through the convergence analysis of SPDEs with the noise terms replaced by the projected noises, the well-posedness of the SPDE is established under certain covariance operator-dependent conditions. These SPDEs with projected noises are then numerically approximated with the finite element method. A general error estimate framework is established for the finite element approximations. Based on this framework, optimal error estimates of finite element approximations for elliptic SPDEs driven by power-law noises are obtained. It is shown that with the proposed approach, convergence order of white noise driven SPDEs is improved by half for one-dimensional problems, and by an infinitesimal factor for higher-dimensional problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wellposedness and regularity estimate for stochastic Cahn--Hilliard equation with unbounded noise diffusion

In this article, we consider the stochastic Cahn--Hilliard equation driven by multiplicative space-time white noise with diffusion coefficient of sublinear growth. By introducing the spectral Galerkin method, we first obtain the well-posedness of the approximated equation in finite dimension. Then with the help of the semigroup theory and the factorization method, the approximation processes is shown to possess many desirable properties. Further, we show that the approximation process is strongly convergent in certain Banach space via the interpolation inequality and variational approach. Finally, the global existence and regularity estimate of the unique solution process are proven by means of the strong convergence of the approximation process.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exponential Integrators for Stochastic Maxwell's Equations Driven by Itô Noise

This article presents explicit exponential integrators for stochastic Maxwell's equations driven by both multiplicative and additive noises. By utilizing the regularity estimate of the mild solution, we first prove that the strong order of the numerical approximation is $\frac 12$ for general multiplicative noise. Combing a proper decomposition with the stochastic Fubini's theorem, the strong order of the proposed scheme is shown to be $1$ for additive noise. Moreover, for linear stochastic Maxwell's equation with additive noise, the proposed time integrator is shown to preserve exactly the symplectic structure, the evolution of the energy as well as the evolution of the divergence in the sense of expectation. Several numerical experiments are presented in order to verify our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimal Strong Convergence Rate of a Backward Euler Type Scheme for the Cox--Ingersoll--Ross Model Driven by Fractional Brownian Motion

In this paper, we investigate the optimal strong convergence rate of numerical approximations for the Cox--Ingersoll--Ross model driven by fractional Brownian motion with Hurst parameter $H\in(1/2,1)$. To deal with the difficulties caused by the unbounded diffusion coefficient, we study an auxiliary equation based on Lamperti transformation. By means of Malliavin calculus, we prove that the backward Euler scheme applied to this auxiliary equation ensures the positivity of the numerical solution, and is of strong order one. Furthermore, a numerical approximation for the original model is obtained and converges with the same order.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Approximating Stochastic Evolution Equations with Additive White and Rough Noises

In this paper, we analyze Galerkin approximations for stochastic evolution equations driven by an additive Gaussian noise which is temporally white and spatially fractional with Hurst index less than or equal to $1/2$. First we regularize the noise by the Wong-Zakai approximation and obtain its optimal order of convergence. Then we apply the Galerkin method to discretize the stochastic evolution equations with regularized noises. Optimal error estimates are obtained for the Galerkin approximations. In particular, our error estimates remove an infinitesimal factor which appears in the error estimates of various numerical methods for stochastic evolution equations in existing literatures.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Finite element approximations for second order stochastic differential equation driven by fractional Brownian motion

We consider finite element approximations for a one dimensional second order stochastic differential equation of boundary value type driven by a fractional Brownian motion with Hurst index $H\le 1/2$. We make use of a sequence of approximate solutions with the fractional noise replaced by its piecewise con- stant approximations to construct the finite element approximations for the equation. The error estimate of the approximations is derived through rigorous convergence analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Approximation of Invariant Measure for Damped Stochastic Nonlinear Schrödinger Equation via an Ergodic Numerical Scheme

In order to inherit numerically the ergodicity of the damped stochastic nonlinear Schrödinger equation with additive noise, we propose a fully discrete scheme, whose spatial direction is based on spectral Galerkin method and temporal direction is based on a modification of the implicit Euler scheme. We not only prove the unique ergodicity of the numerical solutions of both spatial semi-discretization and full discretization, but also present error estimations on invariant measures, which gives order $2$ in spatial direction and order ${\frac12}$ in temporal direction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence in probability of an ergodic and conformal multi-symplectic numerical scheme for a damped stochastic NLS equation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy-conserving method for Stochastic Maxwell Equations with Multiplicative Noise

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Analysis on Ergodic Limit of Approximations for Stochastic NLS Equation via Multi-symplectic Scheme

We consider a finite dimensional approximation of the stochastic nonlinear Schrödinger equation driven by multiplicative noise, which is derived by applying a symplectic method to the original equation in spatial direction. Both the unique ergodicity and the charge conservation law for this finite dimensional approximation are obtained on the unit sphere. To simulate the ergodic limit over long time for the finite dimensional approximation, we discretize it further in temporal direction to obtain a fully discrete scheme, which inherits not only the stochastic multi-symplecticity and charge conservation law of the original equation but also the unique ergodicity of the finite dimensional approximation. The temporal average of the fully discrete numerical solution is proved to converge to the ergodic limit with order one with respect to the time step for a fixed spatial step. Numerical experiments verify our theoretical results on charge conservation, ergodicity and weak convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Projection methods for stochastic differential equations with conserved quantities

In this paper, we consider the numerical methods preserving single or multiple conserved quantities, and these methods are able to reach high order of strong convergence simultaneously based on some kinds of projection methods. The mean-square convergence orders of these methods under certain conditions are given, which can reach order 1.5 or even 2 according to the supporting methods embedded in the projection step. Finally, three numerical experiments are taken into account to show the superiority of the projection methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Symplectic Runge-Kutta Semi-discretization for Stochastic Schrödinger Equation

Based on a variational principle with a stochastic forcing, we indicate that the stochastic Schrödinger equation in Stratonovich sense is an infinite-dimensional stochastic Hamiltonian system, whose phase flow preserves symplecticity. We propose a general class of stochastic symplectic Runge-Kutta methods in temporal direction to the stochastic Schrödinger equation in Stratonovich sense and show that the methods preserve the charge conservation law. We present a convergence theorem on the relationship between the mean-square convergence order of a semi-discrete method and its local accuracy order. Taking stochastic midpoint scheme as an example of stochastic symplectic Runge-Kutta methods in temporal direction, based on the theorem we show that the mean-square convergence order of the semi-discrete scheme is 1 under appropriate assumptions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of a $θ$-scheme to solve the stochastic nonlinear Schrödinger equation with Stratonovich noise

We propose a $θ$-scheme to discretize the $d$-dimensional stochastic cubic Schrödinger equation in Stratono\-vich sense. A uniform bound for the Hamiltonian of the discrete problem is obtained, which is a crucial property to verify the convergence in probability towards a mild solution. Furthermore, based on the uniform bounds of iterates in ${\mathbb H}^2(\mathcal{O})$ for $\mathcal{O}\subset\mathbb{R}^{1}$, the optimal convergence order 1 in strong local sense is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mean-square Convergence of a Symplectic Local Discontinuous Galerkin Method Applied to Stochastic Linear Schroedinger Equation

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Preservation of Physical Properties of Stochastic Maxwell Equations with Additive Noise via Stochastic Multi-symplectic Methods

Stochastic Maxwell equations with additive noise are a system of stochastic Hamiltonian partial differential equations intrinsically, possessing the stochastic multi-symplectic conservation law.It is shown that the averaged energy increases linearly with respect to the evolution of time and the flow of stochastic Maxwell equations with additive noise preserves the divergence in the sense of expectation. Moreover, we propose three novel stochastic multi-symplectic methods to discretize stochastic Maxwell equations in order to investigate the preservation of these properties numerically. We made theoretical discussions and comparisons on all of the three methods to observe that all of them preserve the corresponding discrete version of the averaged divergence. Meanwhile, we obtain the corresponding dissipative property of the discrete averaged energy satisfied by each method. Especially, the evolution rates of the averaged energies for all of the three methods are derived which are in accordance with the continuous case. Numerical experiments are performed to verify our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conservative methods for stochastic differential equations with a conserved quantity

This paper proposes a novel conservative method for numerical computation of general stochastic differential equations in the Stratonovich sense with a conserved quantity. We show that the mean-square order of the method is $1$ if noises are commutative and that the weak order is also $1$. Since the proposed method may need the computation of a deterministic integral, we analyse the effect of the use of quadrature formulas on the convergence orders. Furthermore, based on the splitting technique of stochastic vector fields, we construct conservative composition methods with similar orders as the above method. Finally, numerical experiments are presented to support our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modified equations for weak stochastic symplectic schemes via their generating functions

In this paper, a systematic approach of constructing modified equations for weak stochastic symplectic methods of stochastic Hamiltonian systems is given via using the generating functions of the stochastic symplectic methods. This approach is valid for stochastic Hamiltonian systems with either additive noises or half-multiplicative noises, and we prove that the modified equation of the weak stochastic symplectic methods are perturbed stochastic Hamiltonian systems of the original systems, which reveals in certain sense the reason for the good long time numerical behavior of stochastic symplectic methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Near Preservation of Quadratic Invariants by Stochastic Runge-Kutta Methods

Based on the combinatory theory of rooted colored trees, we investigate the conditions for the explicit stochastic Runge-Kutta (SRK) methods to preserve quadratic invariants (QI) up to certain orders of accuracy. These conditions can supply a practical approach of constructing explicit nearly conservative SRK methods. Meanwhile, we estimate errors in the preservation of QI resulting from iterative implementation of implicit conservative SRK methods with fixed-point and Newton's iterations. Finally, numerical experiments are performed to test the behavior of the methods in preserving QI.

0 citations0 reviewsNo code linkedOpen access