Solution of third-order Emden–Fowler-type equations using wavelet methods

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Arshad Khan ◽  
Mo Faheem ◽  
Akmal Raza

Purpose The numerical solution of third-order boundary value problems (BVPs) has a great importance because of their applications in fluid dynamics, aerodynamics, astrophysics, nuclear reactions, rocket science etc. The purpose of this paper is to develop two computational methods based on Hermite wavelet and Bernoulli wavelet for the solution of third-order initial/BVPs. Design/methodology/approach Because of the presence of singularity and the strong nonlinear nature, most of third-order BVPs do not occupy exact solution. Therefore, numerical techniques play an important role for the solution of such type of third-order BVPs. The proposed methods convert third-order BVPs into a system of algebraic equations, and on solving them, approximate solution is obtained. Finally, the numerical simulation has been done to validate the reliability and accuracy of developed methods. Findings This paper discussed the solution of linear, nonlinear, nonlinear singular (Emden–Fowler type) and self-adjoint singularly perturbed singular (generalized Emden–Fowler type) third-order BVPs using wavelets. A comparison of the results of proposed methods with the results of existing methods has been given. The proposed methods give the accuracy up to 19 decimal places as the resolution level is increased. Originality/value This paper is one of the first in the literature that investigates the solution of third-order Emden–Fowler-type equations using Bernoulli and Hermite wavelets. This paper also discusses the error bounds of the proposed methods for the stability of approximate solutions.

Author(s):  
Chengdong Yuan ◽  
Siyang Hu ◽  
Tamara Bechtold

Purpose Based on the framework of Krylov subspace-based model order reduction (MOR), compact models of the piezoelectric energy harvester devices can be generated. However, the stability of reduced piezoelectric model often cannot be preserved. In previous research studies, “MOR after Schur,” “Schur after MOR” and “multiphysics structure preserving MOR” methods have proven successful in obtaining stable reduced piezoelectric energy harvester models. Though the stability preservation of “MOR after Schur” and “Schur after MOR” methods has already been mathematically proven, the “multiphysics structure preserving MOR” method was not. This paper aims to provide the missing mathematical proof of “multiphysics structure preserving MOR.” Design/methodology/approach Piezoelectric energy harvesters can be represented by system of differential-algebraic equations obtained by the finite element method. According to the block structure of its system matrices, “MOR after Schur” and “Schur after MOR” both perform Schur complement transformations either before or after the MOR process. For the “multiphysics structure preserving MOR” method, the original block structure of the system matrices is preserved during MOR.  Findings This contribution shows that, in comparison to “MOR after Schur” and “Schur after MOR” methods, “multiphysics structure preserving MOR” method performs the Schur complement transformation implicitly, and therefore, stabilizes the reduced piezoelectric model. Originality/value The stability preservation of the reduced piezoelectric energy harvester model obtained through “multiphysics structure preserving MOR” method is proven mathematically and further validated by numerical experiments on two different piezoelectric energy harvester devices.


2019 ◽  
Vol 14 (9) ◽  
pp. 15
Author(s):  
Nguyen Tien Dung ◽  
Nguyen Anh Tra

In this paper, finite difference schemes are proposed to approximate solutions of stochastic advection-diffusion equations. We used central-difference formula of third-order to approximate spatial derivatives. The stability, consistency and convergence of the scheme are analysed and established. A numerical result is also given to demonstrate the computational efficiency of the stochastic schemes.


2017 ◽  
Vol 27 (8) ◽  
pp. 1687-1705 ◽  
Author(s):  
Velinda Calvert ◽  
Mohsen Razzaghi

Purpose This paper aims to propose a new numerical method for the solution of the Blasius and magnetohydrodynamic (MHD) Falkner-Skan boundary-layer equations. The Blasius and MHD Falkner-Skan equations are third-order nonlinear boundary value problems on the semi-infinite domain. Design/methodology/approach The approach is based upon modified rational Bernoulli functions. The operational matrices of derivative and product of modified rational Bernoulli functions are presented. These matrices together with the collocation method are then utilized to reduce the solution of the Blasius and MHD Falkner-Skan boundary-layer equations to the solution of a system of algebraic equations. Findings The method is computationally very attractive and gives very accurate results. Originality/value Many problems in science and engineering are set in unbounded domains. One approach to solve these problems is based on rational functions. In this work, a new rational function is used to find solutions of the Blasius and MHD Falkner-Skan boundary-layer equations.


Author(s):  
William P. Walters ◽  
Cyril L. Williams

The Alekseevski-Tate equations have long been used to predict the penetration, penetration velocity, rod velocity, and rod erosion of long-rod projectiles or kinetic-energy penetrators [1]. These nonlinear equations were originally solved numerically, then by the exact analytical solution of Walters and Segletes [2, 3]. However, due to the nonlinear nature of the equations, the penetration was obtained implicitly as a function of time, so that an explicit functional dependence of the penetration on material properties was not obtained. Walters and Williams [4, 5, 6] obtained the velocities, length, and penetration as an explicit function of time by employing a perturbation solution of the non-dimensional Alekseevski-Tate equations. Algebraic equations were obtained for a third-order perturbation solution which showed excellent agreement with the exact solution of the Tate equations for tungsten heavy alloy rods penetrating a semi-infinite armor plate. The current paper employs this model to rapidly assess the effect of increasing the impact velocity of the penetrator and increasing the armor material properties (density and target resistance) on penetration. This study is applicable to the design of hardened targets.


Mathematics ◽  
2019 ◽  
Vol 7 (10) ◽  
pp. 923 ◽  
Author(s):  
Abdul Ghafoor ◽  
Sirajul Haq ◽  
Manzoor Hussain ◽  
Poom Kumam ◽  
Muhammad Asif Jan

In this paper, a wavelet based collocation method is formulated for an approximate solution of (1 + 1)- and (1 + 2)-dimensional time fractional diffusion wave equations. The main objective of this study is to combine the finite difference method with Haar wavelets. One and two dimensional Haar wavelets are used for the discretization of a spatial operator while time fractional derivative is approximated using second order finite difference and quadrature rule. The scheme has an excellent feature that converts a time fractional partial differential equation to a system of algebraic equations which can be solved easily. The suggested technique is applied to solve some test problems. The obtained results have been compared with existing results in the literature. Also, the accuracy of the scheme has been checked by computing L 2 and L ∞ error norms. Computations validate that the proposed method produces good results, which are comparable with exact solutions and those presented before.


2020 ◽  
Vol 25 (2) ◽  
pp. 29
Author(s):  
Desmond Adair ◽  
Aigul Nagimova ◽  
Martin Jaeger

The vibration characteristics of a nonuniform, flexible and free-flying slender rocket experiencing constant thrust is investigated. The rocket is idealized as a classic nonuniform beam with a constant one-dimensional follower force and with free-free boundary conditions. The equations of motion are derived by applying the extended Hamilton’s principle for non-conservative systems. Natural frequencies and associated mode shapes of the rocket are determined using the relatively efficient and accurate Adomian modified decomposition method (AMDM) with the solutions obtained by solving a set of algebraic equations with only three unknown parameters. The method can easily be extended to obtain approximate solutions to vibration problems for any type of nonuniform beam.


2012 ◽  
Vol 67 (12) ◽  
pp. 665-673 ◽  
Author(s):  
Kourosh Parand ◽  
Mehran Nikarya ◽  
Jamal Amani Rad ◽  
Fatemeh Baharifard

In this paper, a new numerical algorithm is introduced to solve the Blasius equation, which is a third-order nonlinear ordinary differential equation arising in the problem of two-dimensional steady state laminar viscous flow over a semi-infinite flat plate. The proposed approach is based on the first kind of Bessel functions collocation method. The first kind of Bessel function is an infinite series, defined on ℝ and is convergent for any x ∊ℝ. In this work, we solve the problem on semi-infinite domain without any domain truncation, variable transformation basis functions or transformation of the domain of the problem to a finite domain. This method reduces the solution of a nonlinear problem to the solution of a system of nonlinear algebraic equations. To illustrate the reliability of this method, we compare the numerical results of the present method with some well-known results in order to show the applicability and efficiency of our method.


2016 ◽  
Vol 53 (9) ◽  
pp. 1522-1532 ◽  
Author(s):  
Farshid Vahedifard ◽  
Shahriar Shahrokhabadi ◽  
Dov Leshchinsky

This study presents a methodology to determine the stability and optimal profile for slopes with concave cross section under static and seismic conditions. Concave profiles are observed in some natural slopes suggesting that such geometry is a more stable configuration. In this study, the profile of a concave slope was idealized by a circular arc defined by a single variable, the mid-chord offset (MCO). The proposed concave profile formulation was incorporated into a limit equilibrium–based log spiral slope stability method. Stability charts are presented to show the stability number, MCO, and mode of failure for homogeneous slopes corresponding to the most stable configuration under static and pseudostatic conditions. It is shown that concave profiles can significantly improve the stability of slopes. Under seismic conditions, the impact of concavity is most pronounced. Good agreement was demonstrated upon comparison of the results from the proposed method against those attended from a rigorous upper bound limit analysis. The proposed methodology, along with recent advances in construction technology, can be employed to use concave profiles in trenches, open mine excavations, earth retaining systems, and naturally cemented and stabilized soil slopes. The results presented provide a useful tool for preliminary evaluation for adopting such concave profiles in practice.


1968 ◽  
Vol 78 (1) ◽  
pp. 91-103 ◽  
Author(s):  
G. P. Szegö ◽  
C. Olech ◽  
A. Cellina

1989 ◽  
Vol 111 (2) ◽  
pp. 187-193 ◽  
Author(s):  
C. Nataraj ◽  
H. D. Nelson

A new quantitative method of estimating steady state periodic behavior in nonlinear systems, based on the trigonometric collocation method, is outlined. A procedure is developed to analyze large rotor dynamic systems with nonlinear supports by the use of the above method in conjunction with Component Mode Synthesis. The algorithm discussed is seen to reduce the original problem to solving nonlinear algebraic equations in terms of only the coordinates associated with the nonlinear supports and is a big improvement over commonly used integration methods. The feasibility and advantages of the procedure so developed are illustrated with the help of an example of a typical rotor dynamic system with an uncentered squeeze film damper. Future work on the investigation of the stability of the periodic response so obtained is outlined.


Sign in / Sign up

Export Citation Format

Share Document