Geometric grid network and third-order compact scheme for solving nonlinear variable coefficients 3D elliptic PDEs

Author(s):  
Navnit Jha ◽  
Venu Gopal ◽  
Bhagat Singh

By using nonuniform (geometric) grid network, a new high-order finite-difference compact scheme has been obtained for the numerical solution of three-space dimensions partial differential equations of elliptic type. Single cell discretization to the elliptic equation makes it easier to compute and exhibit stability of the numerical solutions. The monotone and irreducible property of the Jacobian matrix to the system of difference equations analyses the converging behavior of the numerical solution values. As an experiment, applications of the compact scheme to Schrödinger equations, sine-Gordon equations, elliptic Allen–Cahn equation and Poisson’s equation have been presented with root mean squared errors of exact and approximate solution values. The results corroborate the reliability and efficiency of the scheme.

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Weiqiu Pan ◽  
Tianzeng Li ◽  
Safdar Ali

AbstractThe Ebola outbreak in 2014 caused many infections and deaths. Some literature works have proposed some models to study Ebola virus, such as SIR, SIS, SEIR, etc. It is proved that the fractional order model can describe epidemic dynamics better than the integer order model. In this paper, we propose a fractional order Ebola system and analyze the nonnegative solution, the basic reproduction number $R_{0}$ R 0 , and the stabilities of equilibrium points for the system firstly. In many studies, the numerical solutions of some models cannot fit very well with the real data. Thus, to show the dynamics of the Ebola epidemic, the Gorenflo–Mainardi–Moretti–Paradisi scheme (GMMP) is taken to get the numerical solution of the SEIR fractional order Ebola system and the modified grid approximation method (MGAM) is used to acquire the parameters of the SEIR fractional order Ebola system. We consider that the GMMP method may lead to absurd numerical solutions, so its stability and convergence are given. Then, the new fractional orders, parameters, and the root-mean-square relative error $g(U^{*})=0.4146$ g ( U ∗ ) = 0.4146 are obtained. With the new fractional orders and parameters, the numerical solution of the SEIR fractional order Ebola system is closer to the real data than those models in other literature works. Meanwhile, we find that most of the fractional order Ebola systems have the same order. Hence, the fractional order Ebola system with different orders using the Caputo derivatives is also studied. We also adopt the MGAM algorithm to obtain the new orders, parameters, and the root-mean-square relative error which is $g(U^{*})=0.2744$ g ( U ∗ ) = 0.2744 . With the new parameters and orders, the fractional order Ebola systems with different orders fit very well with the real data.


2018 ◽  
Vol 35 (3) ◽  
pp. 1081-1102 ◽  
Author(s):  
Nadeem Haider ◽  
Imran Aziz ◽  
Siraj-ul-Islam

2015 ◽  
Vol 2015 ◽  
pp. 1-10
Author(s):  
M. Fakharany ◽  
R. Company ◽  
L. Jódar

This paper is concerned with the numerical solution of partial integrodifferential equation for option pricing models under a tempered stable process known as CGMY model. A double discretization finite difference scheme is used for the treatment of the unbounded nonlocal integral term. We also introduce in the scheme the Patankar-trick to guarantee unconditional nonnegative numerical solutions. Integration formula of open type is used in order to improve the accuracy of the approximation of the integral part. Stability and consistency are also studied. Illustrative examples are included.


2021 ◽  
Author(s):  
Ping-Cheng Hsieh ◽  
Tzu-Ting Huang

Abstract. This study discussed water storage in aquifers of hillslopes under temporally varied rainfall recharge by employing a hillslope-storage equation to simulate groundwater flow. The hillslope width was assumed to vary exponentially to denote the following complex hillslope types: uniform, convergent, and divergent. Both analytical and numerical solutions were acquired for the storage equation with a recharge source. The analytical solution was obtained using an integral transform technique. The numerical solution was obtained using a finite difference method in which the upwind scheme was used for space derivatives and the third-order Runge–Kutta scheme was used for time discretization. The results revealed that hillslope type significantly influences the drains of hillslope storage. Drainage was the fastest for divergent hillslopes and the slowest for convergent hillslopes. The results obtained from analytical solutions require the tuning of a fitting parameter to better describe the groundwater flow. However, a gap existed between the analytical and numerical solutions under the same scenario owing to the different versions of the hillslope-storage equation. The study findings implied that numerical solutions are superior to analytical solutions for the nonlinear hillslope-storage equation, whereas the analytical solutions are better for the linearized hillslope-storage equation. The findings thus can benefit research on and have application in soil and water conservation.


2020 ◽  
Author(s):  
◽  
Uriel Jacket Tresor Demby's

In the context of articulated robotic manipulators, the Forward Kinematics (FK) is a highly non-linear function that maps joint configurations of the robot to poses of its endeffector. Furthermore, while in the most useful cases these functions are neither injective (one-to-one) nor surjective (onto), depending on the robot configuration -- i.e. the sequence of prismatic versus revolute joints, and the number of Degrees of Freedom (DoF) -- the associated Inverse Kinematics (IK) problem may be practically or even theoretically impossible to be solved analytically. Therefore, in the past decades, several approximate methods have been developed for many instances of IK problems. The approximate methods can be divided into two distinct categories: data-driven and numerical approaches. In the first case, data-driven approaches have been successfully used for small workspace domains (e.g., task-driven applications), but not fully explored for large ones, i.e. in task-independent applications where a more general IK is required. Similarly, and despite many successful implementations over the years, numerical solutions may fail if an improper matrix inverse is employed (e.g., Moore-Penrose generalized inverse). In this research, we propose a systematic, robust and accurate numerical solution for the IK problem using the Unit-Consistent (UC) and the Mixed (MX) Inverse methods to invert the Jacobians derived from the Denavit-Hartenberg (D-H) representation of the FK for any robot. As we demonstrate, this approach is robust to whether the system is underdetermined (less than 6 DoF) or overdetermined (more than 6 DoF). We compare the proposed numerical solution to data driven solutions using different robots -- with DoF varying from 3 to 7. We conclude that numerical solutions are easier to implement, faster, and more accurate than most data-driven approaches in the literature, specially for large workspaces as in task-independent applications. We particularly compared the proposed numerical approach against two data-driven approaches: Multi-Layer Perceptron (MLP) and Adaptive Neuro-Fuzzy Inference System (ANFIS), while exploring various architectures of these Neural Networks (NN): i.e. number of inputs, number of outputs, depth, and number of nodes in the hidden layers.


1990 ◽  
Vol 17 (2) ◽  
pp. 252-261 ◽  
Author(s):  
Kevin R. Hall

The interaction of a wave with a rubblemound breakwater results in a complex flow field which is both nonlinear and turbulent, particularly within a region close to the surface of the structure. Numerical models describing internal flow in a rubblemound breakwater are becoming increasingly important, particularly as the influence of scale effects on internal flow in physical hydraulic models are becoming understood as important. A number of numerical models to predict the internal breakwater flow kinematics have been produced in the past two decades. This paper provides a review of the state-of-the-art of numerical modelling of wave interaction with rubblemound breakwaters. Details of the theoretical development and the resulting numerical solution techniques are presented. Methods for incorporating secondary effects such as two-phase (air–water) flow, inertia, and unbalanced boundary conditions are discussed. Limitations of the models resulting from the validity of the assumptions made in order to effect a numerical solution are discussed. Key words: breakwaters, internal flow, porous media flow, numerical modelling, rubblemound breakwaters.


1998 ◽  
Vol 120 (1) ◽  
pp. 46-51 ◽  
Author(s):  
L. N. Srinivasan ◽  
Q. Jeffrey Ge

This paper presents two algorithms for fine-tuning rational B-spline motions suitable for Computer Aided Design. The problem of fine-tuning of rational motions is studied as that of fine-tuning rational curves in a projective dual three-space, called the image curves. The path-smoothing algorithm automatically detects and smoothes out the third order geometric discontinuities in the path of a cubic rational B-spline image curve. The speed-smoothing algorithm uses a quintic rational spline image curve to obtain a second-order geometric approximation of the path of a cubic rational B-spline image curve while allowing specification of the speed and the rate of change of speed at the key points to obtain a near constant kinetic energy parameterization. The results have applications in Cartesian trajectory planning in robotics, spatial navigation in visualization and virtual reality systems, as well as mechanical system simulation.


1976 ◽  
Vol 20 (03) ◽  
pp. 160-170
Author(s):  
Nils Salvesen ◽  
C. von Kerczek

Numerical solutions of the nonlinear problem of the steady two-dimensional potential flow past a submerged line vortex are obtained using the finite-difference iterative technique previously presented by the authors. These solutions are compared in detail with third-order perturbation theory solutions. It is found that very good agreement is obtained for cases of positive circulation of the vortex with strength large enough to produce downstream waves whose steepness is within 15 percent of the maximum possible steepness of irrotational free waves. These computed waves are as steep as the steepest waves obtained in a certain experiment involving the flow past a two-dimensional hydrofoil. For negative circulation, there is substantial difference between the numerical results and third-order perturbation theory. The failure of the perturbation theory is discussed. Details of the far-downstream wave system obtained by the numerical method are compared with other numerical solutions and very high-order perturbation theory solutions of the free-wave problem. Very good agreement is obtained in most cases.


2021 ◽  
Vol 28 (28) ◽  
pp. 85-93
Author(s):  
I. VOLINSKY VOLINSKY ◽  
M. BERSHADSKY BERSHADSKY

In this paper we deal the model of infection diseases, constructed by G.I.Marchuk. We add distributed control to one of the equations describing this model, and obtain numerical solutions of this model. Several numerical solutions demonstrate advances of distributed control.


2014 ◽  
Vol 2014 ◽  
pp. 1-10
Author(s):  
Suheel Abdullah Malik ◽  
Ijaz Mansoor Qureshi ◽  
Muhammad Amir ◽  
Ihsanul Haq

We present a hybrid heuristic computing method for the numerical solution of nonlinear singular boundary value problems arising in physiology. The approximate solution is deduced as a linear combination of some log sigmoid basis functions. A fitness function representing the sum of the mean square error of the given nonlinear ordinary differential equation (ODE) and its boundary conditions is formulated. The optimization of the unknown adjustable parameters contained in the fitness function is performed by the hybrid heuristic computation algorithm based on genetic algorithm (GA), interior point algorithm (IPA), and active set algorithm (ASA). The efficiency and the viability of the proposed method are confirmed by solving three examples from physiology. The obtained approximate solutions are found in excellent agreement with the exact solutions as well as some conventional numerical solutions.


Sign in / Sign up

Export Citation Format

Share Document