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.

scholarly journals A comparative analysis of methods: mimetics, finite differences and finite elements for 1-dimensional stationary problems

2021 ◽  
Vol 8 (1) ◽  
pp. 1-11
Author(s):  
Abdul Abner Lugo Jiménez ◽  
Guelvis Enrique Mata Díaz ◽  
Bladismir Ruiz
Keyword(s):  
Numerical Methods ◽  
Finite Difference ◽  
Finite Difference Method ◽  
Differential Equations ◽  
Difference Method ◽  
Initial Conditions ◽  
Stationary Regime ◽  
Heat Transfer Fluid ◽  
Mimetic Finite Difference Method ◽  
Linear Differential

Numerical methods are useful for solving differential equations that model physical problems, for example, heat transfer, fluid dynamics, wave propagation, among others; especially when these cannot be solved by means of exact analysis techniques, since such problems present complex geometries, boundary or initial conditions, or involve non-linear differential equations. Currently, the number of problems that are modeled with partial differential equations are diverse and these must be addressed numerically, so that the results obtained are more in line with reality. In this work, a comparison of the classical numerical methods such as: the finite difference method (FDM) and the finite element method (FEM), with a modern technique of discretization called the mimetic method (MIM), or mimetic finite difference method or compatible method, is approached. With this comparison we try to conclude about the efficiency, order of convergence of these methods. Our analysis is based on a model problem with a one-dimensional boundary value, that is, we will study convection-diffusion equations in a stationary regime, with different variations in the gradient, diffusive coefficient and convective velocity.

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.

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