Two Finite Difference Schemes for Multi-Dimensional Fractional Wave Equations with Weakly Singular Solutions

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Jinye Shen ◽  
Martin Stynes ◽  
Zhi-Zhong Sun
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
Mesh Size ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Difference Schemes ◽  
Stability And Convergence ◽  
Initial Time ◽  
Finite Difference Schemes

Abstract A time-fractional initial-boundary value problem of wave type is considered, where the spatial domain is ( 0 , 1 ) d (0,1)^{d} for some d ∈ { 1 , 2 , 3 } d\in\{1,2,3\} . Regularity of the solution 𝑢 is discussed in detail. Typical solutions have a weak singularity at the initial time t = 0 t=0 : while 𝑢 and u t u_{t} are continuous at t = 0 t=0 , the second-order derivative u t ⁢ t u_{tt} blows up at t = 0 t=0 . To solve the problem numerically, a finite difference scheme is used on a mesh that is graded in time and uniform in space with the same mesh size ℎ in each coordinate direction. This scheme is generated through order reduction: one rewrites the differential equation as a system of two equations using the new variable v := u t v:=u_{t} ; then one uses a modified L1 scheme of Crank–Nicolson type for the driving equation. A fast variant of this finite difference scheme is also considered, using a sum-of-exponentials (SOE) approximation for the kernel function in the Caputo derivative. The stability and convergence of both difference schemes are analysed in detail. At each time level, the system of linear equations generated by the difference schemes is solved by a fast Poisson solver, thereby taking advantage of the fast difference scheme. Finally, numerical examples are presented to demonstrate the accuracy and efficiency of both numerical methods.

Stability of Finite-difference Schemes for Semilinear Multidimensional Parabolic Equations

2012 ◽  
Vol 12 (3) ◽  
pp. 289-305 ◽  
Author(s):  
Bosko Jovanovic ◽  
Magdalena Lapinska-Chrzczonowicz ◽  
Aleh Matus ◽  
Piotr Matus
Keyword(s):  
Finite Difference ◽  
Parabolic Equations ◽  
Blow Up ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Difference Schemes ◽  
Finite Difference Schemes ◽  
Nonlinear Odes ◽  
The Stability ◽  
Initial Boundary Value

Abstract Abstract — We have studied the stability of finite-difference schemes approximating initial-boundary value problem (IBVP) for multidimensional parabolic equations with a nonlinear source of a power type. We have obtained simple sufficient input data conditions, in which the solutions of differential and difference problems are globally bounded for all t. It is shown that if these conditions are not satisfied, then the solution can blow-up (go to infinity) in finite time. The lower bound of the blow-up time has been determined. The stability of the difference solution has been proven. In all cases, we used the method of energy inequalities based on the application of the Chaplygin comparison theorem for nonlinear ODEs, Bihari-type inequalities and their discrete analogs.

scholarly journals Conservative Linear Difference Scheme for Rosenau-KdV Equation

2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Jinsong Hu ◽  
Youcai Xu ◽  
Bing Hu
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
Kdv Equation ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
The Difference ◽  
Second Order Convergence ◽  
Initial Boundary Value ◽  
Theoretical Results

A conservative three-level linear finite difference scheme for the numerical solution of the initial-boundary value problem of Rosenau-KdV equation is proposed. The difference scheme simulates two conservative quantities of the problem well. The existence and uniqueness of the difference solution are proved. It is shown that the finite difference scheme is of second-order convergence and unconditionally stable. Numerical experiments verify the theoretical results.

scholarly journals A Family of Finite-Difference Schemes with Discrete Transparent Boundary Conditions for a Parabolic Equation on the Half-Axis

2013 ◽  
Vol 13 (2) ◽  
pp. 119-138
Author(s):  
Alexander Zlotnik ◽  
Natalya Koltsova
Keyword(s):  
Parabolic Equation ◽  
Finite Difference ◽  
Boundary Conditions ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Difference Schemes ◽  
Finite Difference Schemes ◽  
Transparent Boundary Conditions ◽  
Half Axis ◽  
Two Parameters

Abstract. An initial-boundary value problem for the 1D self-adjoint parabolic equation on the half-axis is solved. We study a broad family of two-level finite-difference schemes with two parameters related to averages both in time and space. Stability in two norms is proved by the energy method. Also discrete transparent boundary conditions are rigorously derived for schemes by applying the method of reproducing functions. Results of numerical experiments are included as well.

COMPARISON OF FINITE‐DIFFERENCE SCHEMES FOR THE GROSS‐PITAEVSKII EQUATION

2009 ◽  
Vol 14 (1) ◽  
pp. 109-126 ◽  
Author(s):  
Vyacheslav A. Trofimov ◽  
Nikolai Peskov
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Conservation Laws ◽  
Finite Difference Scheme ◽  
Numerical Solutions ◽  
Difference Schemes ◽  
Finite Difference Schemes ◽  
Pitaevskii Equation ◽  
Conservative Scheme ◽  
Conservative Finite Difference Scheme

A conservative finite‐difference scheme for numerical solution of the Gross‐Pitaevskii equation is proposed. The scheme preserves three invariants of the problem: the L 2 norm of the solution, the impulse functional, and the energy functional. The advantages of the scheme are demonstrated via several numerical examples in comparison with some other well‐known and widely used methods. The paper is organized as follows. In Section 2 we consider three main conservation laws of GPE and derive the evolution equations for first and second moments of a solution of GPE. In Section 3 we define the conservative finite‐difference scheme and prove the discrete analogs of conservation laws. The remainder of Section 3 consists of a brief description of other finite‐difference schemes, which will be compared with the conservative scheme. Section 4 presents the results of numerical solutions of three typical problems related to GPE, obtained by different methods. Comparison of the results confirms the advantages of conservative scheme. And finally we summarize our conclusions in Section 5.

scholarly journals An Average Linear Difference Scheme for the Generalized Rosenau-KdV Equation

2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
Maobo Zheng ◽  
Jun Zhou
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
Kdv Equation ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Discrete Energy ◽  
The Difference ◽  
Initial Boundary Value ◽  
Theoretical Results

An average linear finite difference scheme for the numerical solution of the initial-boundary value problem of Generalized Rosenau-KdV equation is proposed. The existence, uniqueness, and conservation for energy of the difference solution are proved by the discrete energy norm method. It is shown that the finite difference scheme is 2nd-order convergent and unconditionally stable. Numerical experiments verify that the theoretical results are right and the numerical method is efficient and reliable.

The Finite Difference Scheme for 3d Mathematical Modeling of a Wood Drying Process

2001 ◽  
Vol 1 (2) ◽  
pp. 125-137 ◽  
Author(s):  
Raimondas Čiegis ◽  
Vadimas Starikovičius
Keyword(s):  
Finite Difference ◽  
Boundary Conditions ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
Numerical Experiments ◽  
Stability And Convergence ◽  
Finite Difference Schemes ◽  
Drying Process ◽  
The Stability ◽  
Robin Boundary

AbstractThis work discusses issues on the design and analysis of finite difference schemes for 3D modeling the process of moisture motion in the wood. A new finite difference scheme is proposed. The stability and convergence in the maximum norm are proved for Robin boundary conditions. The influence of boundary conditions is investigated, and results of numerical experiments are presented.

scholarly journals On the Computation of the Mountain Torque from Gridded Global Datasets

2008 ◽  
Vol 136 (10) ◽  
pp. 4005-4009 ◽  
Author(s):  
Huei-Ping Huang ◽  
Klaus M. Weickmann
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Spectral Method ◽  
Surface Pressure ◽  
Finite Difference Scheme ◽  
Difference Schemes ◽  
Finite Difference Schemes ◽  
Axial Component ◽  
Practical Applications ◽  
Benchmark Calculation

This note evaluates the numerical schemes used for computing the axial component of the mountain torque from gridded global surface pressure and topography datasets. It is shown that the two formulas of the mountain torque based on (i) an integral of the product of the surface pressure and the gradient of topography, and (ii) an integral of the product of the topography and the surface pressure gradient, should produce identical results if a centered even-ordered finite-difference scheme or the spectral method is used to evaluate the integrand. Noncentered finite-difference schemes are not recommended not only because they produce extremely large errors but also because they produce different results for the two formulas. When compared with the benchmark calculation using the spectral method, it is found that the centered fourth-order finite-difference scheme is an efficient and generally accurate approximation for practical applications. Using the data from NCEP–NCAR reanalysis, the finite-difference schemes generally underestimate the global mountain torque compared to the benchmark. This negative error is interpreted as due to the asymmetry in the distribution of surface pressure and in the steepness of the topography between the western and eastern slopes of the mountains.

scholarly journals A linear finite difference scheme for the generalized dissipative symetric regularized long wave equation with damping

2019 ◽  
Vol 23 (Suppl. 3) ◽  
pp. 719-726 ◽  
Author(s):  
Xi Wang ◽  
Jin-Song Hu ◽  
Hong Zhang
Keyword(s):  
Wave Equation ◽  
Finite Difference ◽  
Difference Scheme ◽  
Finite Difference Scheme ◽  
Numerical Solutions ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Long Wave ◽  
Regularized Long Wave Equation ◽  
Initial Boundary Value

In this paper, we study and analyze a three-level linear finite difference scheme for the initial boundary value problem of the symmetric regularized long wave equation with damping. The proposed scheme has the second accuracy both for the spatial and temporal discretization. The convergence and stability of the numerical solutions are proved by the mathematical induction and the discrete functional analysis. Numerical results are given to verify the accuracy and the efficiency of proposed algorithm.

scholarly journals Factorized difference scheme for two-dimensional subdiffusion equation in nonhomogeneous media

2016 ◽  
Vol 99 (113) ◽  
pp. 1-13 ◽  
Author(s):  
Aleksandra Delic ◽  
Sandra Hodzic ◽  
Bosko Jovanovic
Keyword(s):  
Finite Difference ◽  
Difference Scheme ◽  
Numerical Approximation ◽  
Initial Boundary ◽  
Initial Boundary Value Problem ◽  
Stability And Convergence ◽  
Two Dimensional ◽  
Nonhomogeneous Media ◽  
Subdiffusion Equation ◽  
Initial Boundary Value

A factorized finite-difference scheme for numerical approximation of initial-boundary value problem for two-dimensional subdiffusion equation in nonhomogeneous media is proposed. Its stability and convergence are investigated. The corresponding error bounds are obtained.

Sign in / Sign up

Export Citation Format

Share Document