scholarly journals Stable finite difference method for fractional reaction–diffusion equations by compact implicit integration factor methods

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Rongpei Zhang ◽  
Mingjun Li ◽  
Bo Chen ◽  
Liwei Zhang
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Difference Method ◽  
Reaction Diffusion ◽  
Second Order ◽  
Reaction Diffusion Equations ◽  
Implicit Integration ◽  
Fractional Reaction

AbstractIn this paper we propose a stable finite difference method to solve the fractional reaction–diffusion systems in a two-dimensional domain. The space discretization is implemented by the weighted shifted Grünwald difference (WSGD) which results in a stiff system of nonlinear ordinary differential equations (ODEs). This system of ordinary differential equations is solved by an efficient compact implicit integration factor (cIIF) method. The stability of the second order cIIF scheme is proved in the discrete $L^{2}$ L 2 -norm. We also prove the second-order convergence of the proposed scheme. Numerical examples are given to demonstrate the accuracy, efficiency, and robustness of the method.

scholarly journals A New Approach and Solution Technique to Solve Time Fractional Nonlinear Reaction-Diffusion Equations

2015 ◽  
Vol 2015 ◽  
pp. 1-13 ◽  
Author(s):  
Inci Cilingir Sungu ◽  
Huseyin Demir
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Difference Method ◽  
Reaction Diffusion ◽  
Reaction Diffusion Equations ◽  
Diffusion Equations ◽  
Time Stepping ◽  
Nonlinear Reaction ◽  
Differential Transform ◽  
Generalized Differential

A new application of the hybrid generalized differential transform and finite difference method is proposed by solving time fractional nonlinear reaction-diffusion equations. This method is a combination of the multi-time-stepping temporal generalized differential transform and the spatial finite difference methods. The procedure first converts the time-evolutionary equations into Poisson equations which are then solved using the central difference method. The temporal differential transform method as used in the paper takes care of stability and the finite difference method on the resulting equation results in a system of diagonally dominant linear algebraic equations. The Gauss-Seidel iterative procedure then used to solve the linear system thus has assured convergence. To have optimized convergence rate, numerical experiments were done by using a combination of factors involving multi-time-stepping, spatial step size, and degree of the polynomial fit in time. It is shown that the hybrid technique is reliable, accurate, and easy to apply.

scholarly journals Timoshenko Beam Theory:First-Order Analysis, Second-Order Analysis, Stability, and Vibration Analysis Using the Finite Difference Method

2021 ◽  
Author(s):  
Valentin Fogang
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Timoshenko Beam ◽  
Difference Method ◽  
Vibration Analysis ◽  
Beam Theory ◽  
Second Order ◽  
Order Analysis ◽  
The Finite Difference Method

This paper presents an approach to the Timoshenko beam theory (TBT) using the finite difference method (FDM). The TBT covers cases associated with small deflections based on shear deformation considerations, whereas the Euler–Bernoulli beam theory neglects shear deformations. The FDM is an approximate method for solving problems described with differential or partial differential equations. It does not involve solving differential equations; equations are formulated with values at selected points of the structure. The model developed in this paper consists of formulating partial differential equations with finite differences and introducing new points (additional or imaginary points) at boundaries and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness). The introduction of additional points allows satisfying boundary and continuity conditions. First-order, second-order, and vibration analyses of structures were conducted with this model. Efforts, displacements, stiffness matrices, buckling loads, and vibration frequencies were determined. In addition, tapered beams were analyzed (e.g., element stiffness matrix, second-order analysis, and vibration analysis). Finally, the direct time integration method (DTIM) was presented. The FDM-based DTIM enabled the analysis of forced vibration of structures, considering the damping. The efforts and displacements could be determined at any time.

scholarly journals CALCULATION OF A THREE-LAYER CYLINDRICAL SHELL TAKING THE CREEP INTO ACCOUNT

2019 ◽  
Vol 6 (4) ◽  
pp. 14-18 ◽  
Author(s):  
Антон Чепурненко ◽  
Anton Chepurnenko ◽  
Батыр Языев ◽  
Batyr Yazyev ◽  
Анастасия Лапина ◽  
...  
Keyword(s):  
Cylindrical Shell ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Difference Method ◽  
Cylindrical Shells ◽  
Euler Method ◽  
Axisymmetric Loading ◽  
The Finite Difference Method

The article presents the derivation of the resolving equations for the calculation of three-layer cylindrical shells under axisymmetric loading, taking into account creep. The problem is reduced to a system of two ordinary differential equations. The solution is performed numerically using the finite difference method in combination with the Euler method.

A method for the estimation of errors propagated in the numerical solution of a system of ordinary differential equations

1966 ◽  
Vol 25 ◽  
pp. 281-287 ◽  
Author(s):  
P. E. Zadunaisky
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Asymptotic Behaviour ◽  
Numerical Solution ◽  
Ordinary Differential Equations ◽  
Difference Method ◽  
Initial Conditions ◽  
Current Theory ◽  
Numerical Process

Let bex′=f(t,x) a system of ordinary differential equations, with initial conditionsx(a) =s, which is integrated numerically by a finite difference method of orderpand constant steph.To estimate the truncation and round-off errors accumulated during the numerical process it is established a method based on the current theory of the asymptotic behaviour (whenh→0) of errors. This method should avoid the main difficulties that arise when the results of the theory must be applied to practical cases. The method has been successfully tested and applied to estimate the errors accumulated in a numerical computation of planetary perturbations on the orbit of a comet.

A new finite difference method for computing approximate solutions of boundary value problems including transition conditions

2021 ◽  
Vol 102 (2) ◽  
pp. 54-61
Author(s):  
S. Çavuşoğlu ◽  
◽  
O.Sh. Mukhtarov ◽  
◽  
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Numerical Solution ◽  
Ordinary Differential Equations ◽  
Boundary Value Problems ◽  
Difference Method ◽  
Numerical Solutions ◽  
Boundary Value ◽  
Approximate Solutions

This article is aimed at computing numerical solutions of new type of boundary value problems (BVPs) for two-linked ordinary differential equations. The problem studied here differs from the classical BVPs such that it contains additional conditions at the point of interaction, so-called transition conditions. Naturally, such type of problems is much more complicated to solve than classical problems. It is not clear how to apply the classical numerical methods to such type of boundary value transition problems (BVTPs). Based on the finite difference method (FDM) we have developed a new numerical algorithm for computing numerical solution of BVTPs for two-linked ordinary differential equations. To demonstrate the reliability and efficiency of the presented algorithm we obtained numerical solution of one BVTP and the results are compared with the corresponding exact solution. The maximum absolute errors (MAEs) are presented in a table.

Non-Linear Static Analysis of 2D Steep Wave Riser Under Current Load

2016 ◽  
Author(s):  
Hongdong Qiao ◽  
Weidong Ruan ◽  
Zhaohui Shang ◽  
Yong Bai
Keyword(s):  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Difference Method ◽  
Shooting Method ◽  
Non Linear ◽  
Newton Raphson ◽  
Raphson Method ◽  
Steep Wave

A new solution combining finite difference method and shooting method is developed to analyze the behavior of steep wave riser subjected to current loading. Based on the large deformation beam theory and mechanics equilibrium principle, a set of non-linear ordinary differential equations describing the motion of the steep wave riser are obtained. Then, finite difference method and shooting method are adopted and combined to solve the ordinary differential equations with zero moment boundary conditions at both the seabed end and surface end of the steep wave riser. The resulting non-linear finite difference formulations can be solved effectively by Newton-Raphson method. To improve iterative efficiency, shooting method is also employed to obtain the initial value for Newton-Raphson method. Results are compared with that of FEM by OrcaFlex, to verify the accuracy and reliability of the numerical method. Finally, a series of sensitivity analyses are also performed to highlight the influencing parameters in the steep wave riser.

Sign in / Sign up

Export Citation Format

Share Document