scholarly journals Numerical Solutions of the Mathematical Models on the Digestive System and COVID-19 Pandemic by Hermite Wavelet Technique

Symmetry ◽  
2021 ◽  
Vol 13 (12) ◽  
pp. 2428
Author(s):  
Kumbinarasaiah Srinivasa ◽  
Haci Mehmet Baskonus ◽  
Yolanda Guerrero Sánchez

This article developed a functional integration matrix via the Hermite wavelets and proposed a novel technique called the Hermite wavelet collocation method (HWM). Here, we studied two models: the coupled system of an ordinary differential equation (ODE) is modeled on the digestive system by considering different parameters such as sleep factor, tension, food rate, death rate, and medicine. Here, we discussed how these parameters influence the digestive system and showed them through figures and tables. Another fractional model is used on the COVID-19 pandemic. This model is defined by a system of fractional-ODEs including five variables, called S (susceptible), E (exposed), I (infected), Q (quarantined), and R (recovered). The proposed wavelet technique investigates these two models. Here, we express the modeled equation in terms of the Hermite wavelets along with the collocation scheme. Then, using the properties of wavelets, we convert the modeled equation into a system of algebraic equations. We use the Newton–Raphson method to solve these nonlinear algebraic equations. The obtained results are compared with numerical solutions and the Runge–Kutta method (R–K method), which is expressed through tables and graphs. The HWM computational time (consumes less time) is better than that of the R–K method.

2020 ◽  
Vol 31 (07) ◽  
pp. 2050093
Author(s):  
M. M. Khader ◽  
Mohammed M. Babatin

This paper is devoted to introduce an efficient solver using the Hermite collocation technique (HCT), of the coupled system of fractional differential equations (FDEs). The given systems are of basic importance in modeling various phenomena like Cascades and Compartment Analysis, Pond Pollution, Home Heating, Chemostats, and Microorganism Culturing, Nutrient Flow in an Aquarium, Biomass Transfer, Forecasting Prices, Electrical Network, Earthquake Effects on Buildings. The proposed method reduces the system of FDEs to a system of algebraic equations in the coefficients of the expansion using the Hermite polynomials. The introduced method is computer oriented and provides highly accurate solution. To demonstrate the efficiency of the method, two examples are solved and the results are displayed graphically. Finally, we convert the presented coupled systems from the case of its standard form to a first-order ordinary differential equations to compare the obtained numerical solutions with those solutions using the fourth-order Runge–Kutta method (RK4).


Author(s):  
S.C. Shiralashetti ◽  
A.B. Deshi

In this paper, numerical solutions of Riccati and fractional Riccati differential equations are obtained by the Haar wavelet collocation method. An operational matrix of integration based on the Haar wavelet is established, and the procedure for applying the matrix to solve these equations. The fundamental idea of Haar wavelet method is to convert the proposed differential equations into a group of non-linear algebraic equations. The accuracy of approximate solution can be further improved by increasing the level of resolution and an error analysis is computed. The examples are given to demonstrate the fast and flexibility of the method. The results obtained are in good agreement with the exact in comparison with existing ones and it is shown that the technique introduced here is robust, easy to apply and is not only enough accurate but also quite stable.


2018 ◽  
Vol 846 ◽  
pp. 166-189 ◽  
Author(s):  
Nicholas R. Buttle ◽  
Ravindra Pethiyagoda ◽  
Timothy J. Moroney ◽  
Scott W. McCue

We consider steady nonlinear free surface flow past an arbitrary bottom topography in three dimensions, concentrating on the shape of the wave pattern that forms on the surface of the fluid. Assuming ideal fluid flow, the problem is formulated using a boundary integral method and discretised to produce a nonlinear system of algebraic equations. The Jacobian of this system is dense due to integrals being evaluated over the entire free surface. To overcome the computational difficulty and large memory requirements, a Jacobian-free Newton–Krylov (JFNK) method is utilised. Using a block-banded approximation of the Jacobian from the linearised system as a preconditioner for the JFNK scheme, we find significant reductions in computational time and memory required for generating numerical solutions. These improvements also allow for a larger number of mesh points over the free surface and the bottom topography. We present a range of numerical solutions for both subcritical and supercritical regimes, and for a variety of bottom configurations. We discuss nonlinear features of the wave patterns as well as their relationship to ship wakes.


2017 ◽  
Vol 13 (2) ◽  
pp. 4657-4670
Author(s):  
W. S. Amer

This work touches two important cases for the motion of a pendulum called Sub and Ultra-harmonic cases. The small parameter method is used to obtain the approximate analytic periodic solutions of the equation of motion when the pivot point of the pendulum moves in an elliptic path. Moreover, the fourth order Runge-Kutta method is used to investigate the numerical solutions of the considered model. The comparison between both the analytical solution and the numerical ones shows high consistency between them.


Author(s):  
V. F. Edneral ◽  
O. D. Timofeevskaya

Introduction:The method of resonant normal form is based on reducing a system of nonlinear ordinary differential equations to a simpler form, easier to explore. Moreover, for a number of autonomous nonlinear problems, it is possible to obtain explicit formulas which approximate numerical calculations of families of their periodic solutions. Replacing numerical calculations with their precalculated formulas leads to significant savings in computational time. Similar calculations were made earlier, but their accuracy was insufficient, and their complexity was very high.Purpose:Application of the resonant normal form method and a software package developed for these purposes to fourth-order systems in order to increase the calculation speed.Results:It has been shown that with the help of a single algorithm it is possible to study equations of high orders (4th and higher). Comparing the tabulation of the obtained formulas with the numerical solutions of the corresponding equations shows good quantitative agreement. Moreover, the speed of calculation by prepared approximating formulas is orders of magnitude greater than the numerical calculation speed. The obtained approximations can also be successfully applied to unstable solutions. For example, in the Henon — Heyles system, periodic solutions are surrounded by chaotic solutions and, when numerically integrated, the algorithms are often unstable on them.Practical relevance:The developed approach can be used in the simulation of physical and biological systems.


2021 ◽  
Vol 5 (1) ◽  
pp. 8
Author(s):  
Cundi Han ◽  
Yiming Chen ◽  
Da-Yan Liu ◽  
Driss Boutat

This paper applies a numerical method of polynomial function approximation to the numerical analysis of variable fractional order viscoelastic rotating beam. First, the governing equation of the viscoelastic rotating beam is established based on the variable fractional model of the viscoelastic material. Second, shifted Bernstein polynomials and Legendre polynomials are used as basis functions to approximate the governing equation and the original equation is converted to matrix product form. Based on the configuration method, the matrix equation is further transformed into algebraic equations and numerical solutions of the governing equation are obtained directly in the time domain. Finally, the efficiency of the proposed algorithm is proved by analyzing the numerical solutions of the displacement of rotating beam under different loads.


2021 ◽  
Vol 2 (3) ◽  
pp. 431-441
Author(s):  
Odysseas Kosmas

In previous works we developed a methodology of deriving variational integrators to provide numerical solutions of systems having oscillatory behavior. These schemes use exponential functions to approximate the intermediate configurations and velocities, which are then placed into the discrete Lagrangian function characterizing the physical system. We afterwards proved that, higher order schemes can be obtained through the corresponding discrete Euler–Lagrange equations and the definition of a weighted sum of “continuous intermediate Lagrangians” each of them evaluated at an intermediate time node. In the present article, we extend these methods so as to include Lagrangians of split potential systems, namely, to address cases when the potential function can be decomposed into several components. Rather than using many intermediate points for the complete Lagrangian, in this work we introduce different numbers of intermediate points, resulting within the context of various reliable quadrature rules, for the various potentials. Finally, we assess the accuracy, convergence and computational time of the proposed technique by testing and comparing them with well known standards.


2014 ◽  
Vol 2014 ◽  
pp. 1-9 ◽  
Author(s):  
E. H. Doha ◽  
D. Baleanu ◽  
A. H. Bhrawy ◽  
R. M. Hafez

A new Legendre rational pseudospectral scheme is proposed and developed for solving numerically systems of linear and nonlinear multipantograph equations on a semi-infinite interval. A Legendre rational collocation method based on Legendre rational-Gauss quadrature points is utilized to reduce the solution of such systems to systems of linear and nonlinear algebraic equations. In addition, accurate approximations are achieved by selecting few Legendre rational-Gauss collocation points. The numerical results obtained by this method have been compared with various exact solutions in order to demonstrate the accuracy and efficiency of the proposed method. Indeed, for relatively limited nodes used, the absolute error in our numerical solutions is sufficiently small.


Author(s):  
Miguel Moreira ◽  
Jose´ Antunes

Fluid-coupling effects lead to a complex dynamical behavior of immersed spent fuel assembly storage racks. Predicting their responses under strong earthquakes is of prime importance for the safety of nuclear plant facilities. In the near-past we introduced a simplified linearized model for the vibrations of such systems, in which gap-averaged velocity and pressure fields were described analytically in terms of a single space-coordinate for each fluid inter-rack channel. Using such approach it was possible to generate and assemble a complete set of differential-algebraic equations describing the multi-rack fluid coupled system dynamics. Because of the linearization assumptions, we achieved computation of the flow-structure coupled modes, but also time-domain simulations of the system responses. However, nonlinear squeeze-film and dissipative flow effects, connected with very large amplitude responses and/or relatively small water gaps, cannot be properly accounted unless the linearization assumption is relaxed. Such is the aim of the present paper. Here, using a similar approach, we generalize our theoretical model to deal with nonlinear flow effects. Besides that the proposed methodology can be automatically implemented in a symbolic computational environment, it is much less computer-intensive than finite element formulations. Using the proposed technique, computations of basic flow-coupled rack configurations subjected to impulse excitations are presented, in order to highlight the essential features of such systems as well as the relevance of squeeze-film and dissipative effects. Finally, more realistic simulations of complex system responses to strong seismic excitations are presented and discussed.


Author(s):  
S. C. Sinha ◽  
Der-Ho Wu ◽  
Vikas Juneja ◽  
Paul Joseph

Abstract In this paper a general method for the analysis of multidimensional second-order dynamic systems with periodically varying parameters is presented. The state vector and the periodic matrices appearing in the equations are expanded in Chebyshev polynomials over the principal period and the original differential problem is reduced to a set of linear algebraic equations. The technique is suitable for constructing either numerical or approximate analytical solutions. As an illustrative example, approximate analytical expressions for the Floquet characteristic exponents of Mathieu’s equation are obtained. Stability charts are drawn to compare the results the proposed method with those obtained by Runge-Kutta and perturbation methods. Numerical solutions for the flap-lag motion of a three blade helicopter rotor are constructed in the next example. The numerical accuracy and efficiency of the proposed technique is compared with standard numerical codes based on Runge-Kutta, Adams-Moulton and Gear algorithms. The results obtained in the both examples indicate that the suggested approach extremely accurate and is by far the most efficient one.


Sign in / Sign up

Export Citation Format

Share Document