scholarly journals Stability Analysis of the Explicit Difference Scheme for Richards Equation

Entropy ◽  
2020 ◽  
Vol 22 (3) ◽  
pp. 352
Author(s):  
Fengnan Liu ◽  
Yasuhide Fukumoto ◽  
Xiaopeng Zhao

A stable explicit difference scheme, which is based on forward Euler format, is proposed for the Richards equation. To avoid the degeneracy of the Richards equation, we add a perturbation to the functional coefficient of the parabolic term. In addition, we introduce an extra term in the difference scheme which is used to relax the time step restriction for improving the stability condition. With the augmented terms, we prove the stability using the induction method. Numerical experiments show the validity and the accuracy of the scheme, along with its efficiency.

2020 ◽  
Vol 25 (6) ◽  
pp. 997-1014
Author(s):  
Ozgur Yildirim ◽  
Meltem Uzun

In this paper, we study the existence and uniqueness of weak solution for the system of finite difference schemes for coupled sine-Gordon equations. A novel first order of accuracy unconditionally stable difference scheme is considered. The variational method also known as the energy method is applied to prove unique weak solvability.We also present a new unified numerical method for the approximate solution of this problem by combining the difference scheme and the fixed point iteration. A test problem is considered, and results of numerical experiments are presented with error analysis to verify the accuracy of the proposed numerical method.


2013 ◽  
Vol 732-733 ◽  
pp. 487-490
Author(s):  
Zhen Ye Wang ◽  
Jiang Fei Li ◽  
Lian Yuan ◽  
Zhi Zhong Fu ◽  
Bo Li ◽  
...  

In this paper, explicit difference scheme, implicit difference scheme and characteristics method are separately used to simulate the transient flow in products pipeline. The simulation result can be used to prevent water hammer in the pipeline of unsteady situation and to improve the efficiency and safety in oil transmission systems. And then, the stability and accuracy of the three methods are compared by adopting different time steps. For explicit difference method, large fluctuation may occur in case of large time step. For implicit method, the result is weakly affected by time step, only if the relaxation factor selected is reasonable. For characteristics method, the results have a high convergence speed and precision. The results show that, in the situation of valve shut down in terminal, it takes about 1.1×104 seconds to return to a new steady state.


2007 ◽  
Vol 7 (3) ◽  
pp. 208-220 ◽  
Author(s):  
P. Matus ◽  
O. Korolyova ◽  
M. Chuiko

Abstract A priory estimates of the stability in the sense of the initial data of the difference scheme approximating weakly compressible liquid equations in the Riemann invariants have been obtained. These estimates have been proved without any assumptions about the properties of the solution of the differential problem and depend only on the behavior of the initial conditions. As distinct from linear problems, the obtained estimates of stability in the general case exist only for a finite instant of time t≤t_0. In particular, this is confirmed by the fact, that nonfulfilment of these stability conditions lead to the appearance of supersonic flows or domains with large gradients. The questions of uniqueness and convergence of the difference solution are considered also. The results of the computating experiment confirming the theoretical conclusions are given.


2016 ◽  
Vol 26 (3) ◽  
pp. 429-435 ◽  
Author(s):  
Roman I. Parovik

Abstract The paper deals with the model of variable-order nonlinear hereditary oscillator based on a numerical finite-difference scheme. Numerical experiments have been carried out to evaluate the stability and convergence of the difference scheme. It is argued that the approximation, stability and convergence are of the first order, while the scheme is stable and converges to the exact solution.


2012 ◽  
Vol 2012 ◽  
pp. 1-18 ◽  
Author(s):  
Maksat Ashyraliyev

The integral-differential equation of the parabolic type in a Banach space is considered. The unique solvability of this equation is established. The stability estimates for the solution of this equation are obtained. The difference scheme approximately solving this equation is presented. The stability estimates for the solution of this difference scheme are obtained.


2001 ◽  
Vol 1 (2) ◽  
pp. 125-137 ◽  
Author(s):  
Raimondas Čiegis ◽  
Vadimas Starikovičius

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.


2013 ◽  
Vol 2013 ◽  
pp. 1-8 ◽  
Author(s):  
Allaberen Ashyralyev ◽  
Mesut Urun

The first order of accuracy difference scheme for the numerical solution of the boundary value problem for the differential equation with parameterp,i(du(t)/dt)+Au(t)+iu(t)=f(t)+p,0<t<T,u(0)=φ,u(T)=ψ, in a Hilbert spaceHwith self-adjoint positive definite operatorAis constructed. The well-posedness of this difference scheme is established. The stability inequalities for the solution of difference schemes for three different types of control parameter problems for the Schrödinger equation are obtained.


2021 ◽  
Vol 4 (4(112)) ◽  
pp. 47-56
Author(s):  
Rakhmatillo Aloev ◽  
Abdumauvlen Berdyshev ◽  
Aziza Akbarova ◽  
Zharasbek Baishemirov

The problem of numerical determination of Lyapunov-stable (exponential stability) solutions of the Saint-Venant equations system has remained open until now. The authors of this paper previously proposed an implicit upwind difference splitting scheme, but its practical applicability was not indicated there. In this paper, the problem is solved successfully, namely, an algorithm for calculating Lyapunov-stable solutions of the Saint-Venant equations system is developed and implemented using an upwind implicit difference splitting scheme on the example of the Big Almaty Canal (hereinafter BAC). As a result of the proposed algorithm application, it was established that: 1) we were able to perform a computational calculation of the numerical determination problem of the water level and velocity on a part of the BAC (10,000 meters) located in the Almaty region; 2) the numerical values of the water level height and horizontal velocity are consistent with the actual measurements of the parameters of the water flow in the BAC; 3) the proposed computational algorithm is stable; 4) the numerical stationary solution of the system of Saint-Venant equations on the example of the BAC is Lyapunov-stable (exponentially stable); 5) the obtained results (according to the BAC) show the efficiency of the developed algorithm based on an implicit upwind difference scheme according to the calculated time. Since we managed to increase the values of the difference grid time step up to 0.8 for calculating the numerical solution according to the proposed implicit scheme.


Sign in / Sign up

Export Citation Format

Share Document