Richardson Extrapolation: An Info-Gap Analysis of Numerical Uncertainty

Author(s):  
Yakov Ben-Haim ◽  
François Hemez

Abstract Computational modeling and simulation is a central tool in science and engineering, directed at solving partial differential equations for which analytical solutions are unavailable. The continuous equations are generally discretized in time, space, energy, etc., to obtain approximate solutions using a numerical method. The aspiration is for the numerical solutions to asymptotically converge to the exact-but-unknown solution as the discretization size approaches zero. A generally applicable procedure to assure convergence is unavailable. The Richardson extrapolation is the main method for dealing with this challenge, but its assumptions introduce uncertainty to the resulting approximation. We use info-gap decision theory to model and manage its main uncertainty, namely, in the rate of convergence of numerical solutions. The theory is illustrated with a numerical application to Hertz contact in solid mechanics.

Symmetry ◽  
2021 ◽  
Vol 13 (12) ◽  
pp. 2296
Author(s):  
Mariam Sultana ◽  
Uroosa Arshad ◽  
Md. Nur Alam ◽  
Omar Bazighifan ◽  
Sameh Askar ◽  
...  

Symmetry performs an essential function in finding the correct techniques for solutions to time space fractional differential equations (TSFDEs). In this article, we present the Novel Analytic Method (NAM) for approximate solutions of the linear and non-linear KdV equation for TSFDs. To enunciate the non-integer derivative for the aforementioned equation, the Caputo operator is manipulated. Furthermore, the formula implemented is a numerical way that is postulated from Taylor’s series, which confirms an analytical answer in the form of a convergent series. For delineation of the efficiency and functionality of the method in question, four applications are exemplified along with graphical interpretation and numerical solutions to finitely illustrate the behavior of the solution to this equation. Moreover, the 3D graphs of some of these numerical examples are plotted with specific values. Observing the effectiveness of this process, we can easily decide that this process can be implemented to other TSFDEs applied in the mathematical modeling of a real-world aspect.


2020 ◽  
Vol 9 (3) ◽  
pp. 633-644
Author(s):  
A. K. Mittal

Abstract In this paper, a new numerical technique implements on the time-space pseudospectral method to approximate the numerical solutions of nonlinear time- and space-fractional coupled Burgers’ equation. This technique is based on orthogonal Chebyshev polynomial function and discretizes using Chebyshev–Gauss–Lobbato (CGL) points. Caputo–Riemann–Liouville fractional derivative formula is used to illustrate the fractional derivatives matrix at CGL points. Using the derivatives matrices, the given problem is reduced to a system of nonlinear algebraic equations. These equations can be solved using Newton–Raphson method. Two model examples of time- and space-fractional coupled Burgers’ equation are tested for a set of fractional space and time derivative order. The figures and tables show the significant features, effectiveness, and good accuracy of the proposed method.


2020 ◽  
Vol 28 (3) ◽  
pp. 209-216
Author(s):  
S. Singh ◽  
S. Saha Ray

AbstractIn this article, hybrid Legendre block-pulse functions are implemented in determining the approximate solutions for multi-dimensional stochastic Itô–Volterra integral equations. The block-pulse function and the proposed scheme are used for deriving a methodology to obtain the stochastic operational matrix. Error and convergence analysis of the scheme is discussed. A brief discussion including numerical examples has been provided to justify the efficiency of the mentioned method.


2017 ◽  
Vol 72 (1) ◽  
pp. 59-69 ◽  
Author(s):  
M.M. Fatih Karahan ◽  
Mehmet Pakdemirli

AbstractStrongly nonlinear cubic-quintic Duffing oscillatoris considered. Approximate solutions are derived using the multiple scales Lindstedt Poincare method (MSLP), a relatively new method developed for strongly nonlinear oscillators. The free undamped oscillator is considered first. Approximate analytical solutions of the MSLP are contrasted with the classical multiple scales (MS) method and numerical simulations. It is found that contrary to the classical MS method, the MSLP can provide acceptable solutions for the case of strong nonlinearities. Next, the forced and damped case is treated. Frequency response curves of both the MS and MSLP methods are obtained and contrasted with the numerical solutions. The MSLP method and numerical simulations are in good agreement while there are discrepancies between the MS and numerical solutions.


Author(s):  
Vladimir P. Gerdt ◽  
Mikhail D. Malykh ◽  
Leonid A. Sevastianov ◽  
Yu Ying

The article considers the midpoint scheme as a finite-difference scheme for a dynamical system of the form ̇ = (). This scheme is remarkable because according to Cooper’s theorem, it preserves all quadratic integrals of motion, moreover, it is the simplest scheme among symplectic Runge-Kutta schemes possessing this property. The properties of approximate solutions were studied in the framework of numerical experiments with linear and nonlinear oscillators, as well as with a system of several coupled oscillators. It is shown that in addition to the conservation of all integrals of motion, approximate solutions inherit the periodicity of motion. At the same time, attention is paid to the discussion of introducing the concept of periodicity of an approximate solution found by the difference scheme. In the case of a nonlinear oscillator, each step requires solving a system of nonlinear algebraic equations. The issues of organizing computations using such schemes are discussed. Comparison with other schemes, including those symmetric with respect to permutation of and .̂


Author(s):  
Liviu Ciupitu

Present paper is dealing with the adaptive static balancing of robot or other mechatronic arms that are moving in vertical plane and whose static loads are variable, by using counterweights and springs. Some simple passive and approximate solutions are proposed and an example is shown. The active and exact solutions by using adaptive real time control in the case of unknown variation of static loads are simulated on VIPRO platform developed at Institute of Solid Mechanics of Romanian Academy.


Symmetry ◽  
2020 ◽  
Vol 12 (9) ◽  
pp. 1446
Author(s):  
Weidong Yang ◽  
Xuehui Chen ◽  
Yuan Meng ◽  
Xinru Zhang ◽  
Shiyun Mi

In this paper, we develop a new time-space fractional constitution relation to study the unsteady boundary layer flow over a stretching sheet. For the convenience of calculation, the boundary layer flow is simulated as a symmetrical rectangular area. The implicit difference method combined with an L1-algorithm and shift Grünwald scheme is used to obtain the numerical solutions of the fractional governing equation. The validity and solvability of the present numerical method are analyzed systematically. The numerical results show that the thickness of the velocity boundary layer increases with an increase in the space fractional parameter γ. For a different stress fractional parameter α, the viscoelastic fluid will exhibit viscous or elastic behavior, respectively. Furthermore, the numerical method in this study is validated and can be extended to other time-space fractional boundary layer models.


2003 ◽  
Vol 2003 (2) ◽  
pp. 87-114 ◽  
Author(s):  
J. R. Fernández ◽  
M. Sofonea

We consider the quasistatic Signorini′s contact problem with damage for elastic-viscoplastic bodies. The mechanical damage of the material, caused by excessive stress or strain, is described by a damage function whose evolution is modeled by an inclusion of parabolic type. We provide a variational formulation for the mechanical problem and sketch a proof of the existence of a unique weak solution of the model. We then introduce and study a fully discrete scheme for the numerical solutions of the problem. An optimal order error estimate is derived for the approximate solutions under suitable solution regularity. Numerical examples are presented to show the performance of the method.


Sign in / Sign up

Export Citation Format

Share Document