Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBackward diffusion-wave problem: stability, regularization and approximation
85
0
0.0
(
0
)
Added by
Zhi Zhou
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We aim at the development and analysis of the numerical schemes for approximately solving the backward diffusion-wave problem, which involves a fractional derivative in time with order $alphain(1,2)$. From terminal observations at two time levels, i.e., $u(T_1)$ and $u(T_2)$, we simultaneously recover two initial data $u(0)$ and $u_t(0)$ and hence the solution $u(t)$ for all $t > 0$. First of all, existence, uniqueness and Lipschitz stability of the backward diffusion-wave problem were established under some conditions about $T_1$ and $T_2$. Moreover, for noisy data, we propose a quasi-boundary value scheme to regularize the mildly ill-posed problem, and show the convergence of the regularized solution. Next, to numerically solve the regularized problem, a fully discrete scheme is proposed by applying finite element method in space and convolution quadrature in time. We establish error bounds of the discrete solution in both cases of smooth and nonsmooth data. The error estimate is very useful in practice since it indicates the way to choose discretization parameters and regularization parameter, according to the noise level. The theoretical results are supported by numerical experiments.
rate research
Read More
A thorough backward stability analysis of Hotellings deflation, an explicit external deflation procedure through low-rank updates for computing many eigenpairs of a symmetric matrix, is presented. Computable upper bounds of the loss of the orthogonality of the computed eigenvectors and the symmetric backward error norm of the computed eigenpairs are derived. Sufficient conditions for the backward stability of the explicit external deflation procedure are revealed. Based on these theoretical results, the strategy for achieving numerical backward stability by dynamically selecting the shifts is proposed. Numerical results are presented to corroborate the theoretical analysis and to demonstrate the stability of the procedure for computing many eigenpairs of large symmetric matrices arising from applications.
Irreversible drift-diffusion processes are very common in biochemical reactions. They have a non-equilibrium stationary state (invariant measure) which does not satisfy detailed balance. For the corresponding Fokker-Planck equation on a closed manifold, via Voronoi tessellation, we propose two upwind finite volume schemes with or without the information of the invariant measure. Both two schemes enjoy stochastic $Q$-matrix structures and can be decomposed as a gradient flow part and a Hamiltonian flow part, which enable us to prove unconditional stability, ergodicity and error estimates. Based on two upwind schemes, several numerical examples - including sampling accelerated by a mixture flow, image transformations and simulations for stochastic model of chaotic system - are conducted. These two structure-preserving schemes also give a natural random walk approximation for a generic irreversible drift-diffusion process on a manifold. Thus they can be adapted to manifold-related computations induced from high dimensional molecular dynamics.
In this work, we study the numerical approximation of a class of singular fully coupled forward backward stochastic differential equations. These equations have a degenerate forward component and non-smooth terminal condition. They are used, for example, in the modeling of carbon market[9] and are linked to scalar conservation law perturbed by a diffusion. Classical FBSDEs methods fail to capture the correct entropy solution to the associated quasi-linear PDE. We introduce a splitting approach that circumvent this difficulty by treating differently the numerical approximation of the diffusion part and the non-linear transport part. Under the structural condition guaranteeing the well-posedness of the singular FBSDEs [8], we show that the splitting method is convergent with a rate $1/2$. We implement the splitting scheme combining non-linear regression based on deep neural networks and conservative finite difference schemes. The numerical tests show very good results in possibly high dimensional framework.
This paper is concerned with the introduction of Tikhonov regularization into least squares approximation scheme on $[-1,1]$ by orthonormal polynomials, in order to handle noisy data. This scheme includes interpolation and hyperinterpolation as special cases. With Gauss quadrature points employed as nodes, coefficients of the approximation polynomial with respect to given basis are derived in an entry-wise closed form. Under interpolatory conditions, the solution to the regularized approximation problem is rewritten in forms of two kinds of barycentric interpolation formulae, by introducing only a multiplicative correction factor into both classical barycentric formulae. An $L_2$ error bound and a uniform error bound are derived, providing similar information that Tikhonov regularization is able to reduce the operator norm (Lebesgue constant) and the error term related to the level of noise, both by multiplying a correction factor which is less than one. Numerical examples show the benefits of Tikhonov regularization when data is noisy or data size is relatively small.
In this paper we consider the numerical solution of Boussinesq-Peregrine type systems by the application of the Galerkin finite element method. The structure of the Boussinesq systems is explained and certain alternative nonlinear and dispersive terms are compared. A detailed study of the convergence properties of the standard Galerkin method, using various finite element spaces on unstructured triangular grids, is presented. Along with the study of the Peregrine system, a new Boussinesq system of BBM-BBM type is derived. The new system has the same structure in its momentum equation but differs slightly in the mass conservation equation compared to the Peregrine system. Further, the finite element method applied to the new system has better convergence properties, when used for its numerical approximation. Due to the lack of analytical formulas for solitary wave solutions for the systems under consideration, a Galerkin finite element method combined with the Petviashvili iteration is proposed for the numerical generation of accurate approximations of line solitary waves. Various numerical experiments related to the propagation of solitary and periodic waves over variable bottom topography and their interaction with the boundaries of the domains are presented. We conclude that both systems have similar accuracy when approximate long waves of small amplitude while the Galerkin finite element method is more effective when applied to BBM-BBM type systems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
