547 research outputs found
Analytical Solution of Biological Population of Fractional Differential Equations by Reconstruction of Variational Iteration Method
This article presents a brand-new approximation analytical technique we refer to as the reconstruction of variational iteration method. For the goal of solving fractional biological population option pricing equations, this methodology was created. In certain circumstances, you may actually use the well-known Mittag-Leffer function to get an explicit response. The usage of the three examples below demonstrates the precision and effectiveness of the suggested method. The results show that the RVIM is not only quite straightforward but also very successful at resolving non-linear problems
A New Adjustment of Laplace Transform for Fractional Bloch Equation in NMR Flow
This work purpose suggest a new analytical technique called the fractional homotopy analysis transform method (FHATM) for solving time fractional Bloch NMR (nuclear magnetic resonance) flow equations, which are a set of macroscopic equations that are used for modeling nuclear magnetization as a function of time. The true beauty of this article is the coupling of the homotopy analysis method and the Laplace transform method for systems of fractional differential equations. The solutions obtained by the proposed method indicate that the approach is easy to implement and computationally very attractive
Optimal Perturbation Iteration Method for Solving Fractional Model of Damped Burgers’ Equation
The newly constructed optimal perturbation iteration procedure with Laplace transform is applied to obtain the new approximate semi-analytical solutions of the fractional type of damped Burgers’ equation. The classical damped Burgers’ equation is remodeled to fractional differential form via the Atangana–Baleanu fractional derivatives described with the help of the Mittag–Leffler function. To display the efficiency of the proposed optimal perturbation iteration technique, an extended example is deeply analyzed.This work was supported in part by the Basque Government, through project IT1207-19
On new and improved semi-numerical techniques for solving nonlinear fluid flow problems.
Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2012.Most real world phenomena is modeled by ordinary and/or partial differential equations.
Most of these equations are highly nonlinear and exact solutions are not always possible.
Exact solutions always give a good account of the physical nature of the phenomena modeled.
However, existing analytical methods can only handle a limited range of these equations.
Semi-numerical and numerical methods give approximate solutions where exact solutions are
impossible to find. However, some common numerical methods give low accuracy and may lack
stability. In general, the character and qualitative behaviour of the solutions may not always
be fully revealed by numerical approximations, hence the need for improved semi-numerical
methods that are accurate, computational efficient and robust.
In this study we introduce innovative techniques for finding solutions of highly nonlinear
coupled boundary value problems. These techniques aim to combine the strengths of both
analytical and numerical methods to produce efficient hybrid algorithms. In this work, the
homotopy analysis method is blended with spectral methods to improve its accuracy. Spectral
methods are well known for their high levels of accuracy. The new spectral homotopy analysis
method is further improved by using a more accurate initial approximation to accelerate
convergence. Furthermore, a quasi-linearisation technique is introduced in which spectral
methods are used to solve the linearised equations. The new techniques were used to solve
mathematical models in fluid dynamics.
The thesis comprises of an introductory Chapter that gives an overview of common numerical
methods currently in use. In Chapter 2 we give an overview of the methods used in this
work. The methods are used in Chapter 3 to solve the nonlinear equation governing two-dimensional
squeezing flow of a viscous fluid between two approaching parallel plates and the
steady laminar flow of a third grade fluid with heat transfer through a flat channel. In Chapter
4 the methods were used to find solutions of the laminar heat transfer problem in a rotating
disk, the steady flow of a Reiner-Rivlin fluid with Joule heating and viscous dissipation and
the classical von Kάrmάn equations for boundary layer flow induced by a rotating disk. In
Chapter 5 solutions of steady two-dimensional flow of a viscous incompressible fluid in a
rectangular domain bounded by two permeable surfaces and the MHD viscous flow problem
due to a shrinking sheet with a chemical reaction, were solved using the new methods
An application of modern analytical solution techniques to nonlinear partial differential equations.
Thesis (M.Sc.)-University of KwaZulu-Natal, Pietermaritzburg, 2013.Many physics and engineering problems are modeled by differential equations. In
many instances these equations are nonlinear and exact solutions are difficult to
obtain. Numerical schemes are often used to find approximate solutions. However,
numerical solutions do not describe the qualitative behaviour of mechanical systems
and are insufficient in determining the general properties of certain systems of
equations. The need for analytical methods is self-evident and major developments
were seen in the 1990’s. With the aid of faster processing equipment today, we are
able to compute analytical solutions to highly nonlinear equations that are more
accurate than numerical solutions.
In this study we discuss solutions to nonlinear partial differential equations with
focus on non-perturbation analytical methods. The non-perturbation methods of
choice are the homotopy analysis method (HAM) developed by Shijun Liao and the
variational iteration method (VIM) developed by Ji-Huan He. The aim is to compare the solutions obtained by these modern day analytical methods against each other
focusing on accuracy, convergence and computational efficiency.
The methods were applied to three test problems, namely, the heat equation, Burgers
equation and the Bratu equation. The solutions were compared against both the exact
results as well as solutions generated using the finite difference method, in some cases.
The results obtained show that the HAM successfully produces solutions which are
accurate, faster converging and requires less computational resources than the VIM.
However, the VIM still provides accurate solutions that are also in good agreement
with the closed form solutions of the test problems. The FDM also produced good
results which were used as a further comparison to the analytical solutions. The
findings of this study is in agreement with those published in the literature
Modeling and inversion of seismic data using multiple scattering, renormalization and homotopy methods
Seismic scattering theory plays an important role in seismic forward modeling and is the theoretical foundation for various seismic imaging methods. Full waveform inversion is a powerful technique for obtaining a high-resolution model of the subsurface. One objective of this thesis is to develop convergent scattering series solutions of the Lippmann-Schwinger equation in strongly scattering media using renormalization and homotopy methods. Other objectives of this thesis are to develop efficient full waveform inversion methods of time-lapse seismic data and, to investigate uncertainty quantification in full waveform inversion for anisotropic elastic media based on integral equation approaches and the iterated extended Kalman filter. The conventional Born scattering series is obtained by expanding the Lippmann-Schwinger equation in terms of an iterative solution based on perturbation theory. Such an expansion assumes weak scattering and may have the problems of convergence in strongly scattering media. This thesis presents two scattering series, referred to as convergent Born series (CBS) and homotopy analysis method (HAM) scattering series for frequency-domain seismic wave modeling. For the convergent Born series, a physical interpretation from the renormalization prospective is given. The homotopy scattering series is derived by using homotopy analysis method, which is based on a convergence control parameter and a convergence control operator that one can use to ensure convergence for strongly scattering media. The homotopy scattering scattering series solutions of the Lippmann-Schwinger equation, which is convergent in strongly scattering media. The homotopy scattering series is a kind of unified scattering series theory that includes the conventional and convergent Born series as special cases. The Fast Fourier Transform (FFT) is employed for efficient implementation of matrix-vector multiplication for the convergent Born series and the homotopy scattering series. This thesis presents homotopy methods for ray based seismic modeling in strongly anisotropic media. To overcome several limitations of small perturbations and weak anisotropy in obtaining the traveltime approximations in anisotropic media by expanding the anisotropic eikonal equation in terms of the anisotropic parameters and the elliptically anisotropic eikonal equation based on perturbation theory, this study applies the homotopy analysis method to the eikonal equation. Then this thesis presents a retrieved zero-order deformation equation that creates a map from the anisotropic eikonal equation to a linearized partial differential equation system. The new traveltime approximations are derived by using the linear and nonlinear operators in the retrieved zero-order deformation equation. Flexibility on variable anisotropy parameters is naturally incorporated into the linear differential equations, allowing a medium of arbitrarily anisotropy. This thesis investigates efficient target-oriented inversion strategies for improving full waveform inversion of time-lapse seismic data based on extending the distorted Born iterative T-matrix inverse scattering to a local inversion of a small region of interest (e. g. reservoir under production). The target-oriented approach is more efficient for inverting the monitor data. The target-oriented inversion strategy requires properly specifying the wavefield extrapolation operators in the integral equation formulation. By employing the T-matrix and the Gaussian beam based Green’s function, the wavefield extrapolation for the time-lapse inversion is performed in the baseline model from the survey surface to the target region. I demonstrate the method by presenting numerical examples illustrating the sequential and double difference strategies. To quantify the uncertainty and multiparameter trade-off in the full waveform inversion for anisotropic elastic media, this study applies the iterated extended Kalman filter to anisotropic elastic full waveform inversion based on the integral equation method. The sensitivity matrix is an explicit representation with Green’s functions based on the nonlinear inverse scattering theory. Taking the similarity of sequential strategy between the multi-scale frequency domain full waveform inversion and data assimilation with an iterated extended Kalman filter, this study applies the explicit representation of sensitivity matrix to the the framework of Bayesian inference and then estimate the uncertainties in the full waveform inversion. This thesis gives results of numerical tests with examples for anisotropic elastic media. They show that the proposed Bayesian inversion method can provide reasonable reconstruction results for the elastic coefficients of the stiffness tensor and the framework is suitable for accessing the uncertainties and analysis of parameter trade-offs
- …