No Arabic abstract
This paper addresses the question whether there are numerical schemes for constant-coefficient advection problems that can yield convergent solutions for an infinite time horizon. The motivation is that such methods may serve as building blocks for long-time accurate solutions in more complex advection-dominated problems. After establishing a new notion of convergence in an infinite time limit of numerical methods, we first show that linear methods cannot meet this convergence criterion. Then we present a new numerical methodology, based on a nonlinear jet scheme framework. We show that these methods do satisfy the new convergence criterion, thus establishing that numerical methods exist that converge on an infinite time horizon, and demonstrate the long-time accuracy gains incurred by this property.
In this paper, we present a numerical verification method of solutions for nonlinear parabolic initial boundary value problems. Decomposing the problem into a nonlinear part and an initial value part, we apply Nakaos projection method, which is based on the full-discrete finite element method with constructive error estimates, to the nonlinear part and use the theoretical analysis for the heat equation to the initial value part, respectively. We show some verified examples for solutions of nonlinear problems from initial value to the neighborhood of the stationary solutions, which confirm us the actual effectiveness of our method.
This article is concerned with the derivation of numerical reconstruction schemes for the inverse moving source problem on determining source profiles in (time-fractional) evolution equations. As a continuation of the theoretical result on the uniqueness, we adopt a minimization procedure with regularization to construct iterative thresholding schemes for the reduced backward problems on recovering one or two unknown initial value(s). Moreover, an elliptic approach is proposed to solve a convection equation in the case of two profiles.
A study is presented on the convergence of the computation of coupled advection-diffusion-reaction equations. In the computation, the equations with different coefficients and even types are assigned in two subdomains, and Schwarz iteration is made between the equations when marching from a time level to the next one. The analysis starts with the linear systems resulting from the full discretization of the equations by explicit schemes. Conditions for convergence are derived, and its speedup and the effects of difference in the equations are discussed. Then, it proceeds to an implicit scheme, and a recursive expression for convergence speed is derived. An optimal interface condition for the Schwarz iteration is obtained, and it leads to perfect convergence, that is, convergence within two times of iteration. Furthermore, the methods and analyses are extended to the coupling of the viscous Burgers equations. Numerical experiments indicate that the conclusions, such as the perfect convergence, drawn in the linear situations may remain in the Burgers equations computation.
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation $w = - {mathcal L}^{-1} {mathcal F}(hat{u}) + {mathcal L}^{-1} {mathcal G}(w)$, where ${mathcal L}$ is a linearized operator, ${mathcal F}(hat{u})$ is a residual, and ${mathcal G}(w)$ is a local Lipschitz term. Therefore, the estimations of $| {mathcal L}^{-1} {mathcal F}(hat{u}) |$ and $| {mathcal L}^{-1}{mathcal G}(w) |$ play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement for matrix problems, we represent the inverse operator ${mathcal L}^{-1}$ as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as ${mathcal L}^{-1}$ are presented in the appendix.
In this paper, by employing the asymptotic method, we prove the existence and uniqueness of a smoothing solution for a singularly perturbed Partial Differential Equation (PDE) with a small parameter. As a by-product, we obtain a reduced PDE model with vanished high order derivative terms, which is close to the original PDE model in any order of this small parameter in the whole domain except a negligible transition layer. Based on this reduced forward model, we propose an efficient two step regularization algorithm for solving inverse source problems governed by the original PDE. Convergence rate results are studied for the proposed regularization algorithm, which shows that this simplification will not (asymptotically) decrease the accuracy of the inversion result when the measurement data contains noise. Numerical examples for both forward and inverse problems are given to show the efficiency of the proposed numerical approach.