Analytical approximate solutions for two-dimensional incompressible Navier-Stokes equations

Analytical approximate solutions of the two-dimensional incompressible Navier-Stokes equations by means of Adomian decomposition method are presented. The power of this manageable method is confirmed by applying it for two selected  flow problems: The first is the Taylor decaying vortices, and the second is the flow behind a grid, comparison with High-order upwind compact finite-difference method is made. The numerical results that are obtained for two incompressible flow problems  showed that the proposed method is less time consuming, quite accurate and easily implemented. In addition, we prove the convergence of this method when it is applied to the flow problems, which are describing them by  unsteady two-dimensional incompressible Navier-Stokes equations.   Keywords: Navier-Stokes equations, Adomian decomposition, upwind compact difference, Accuracy, Convergence analysis,Taylor's decay vortices, flow behind a grid

Fractional model of cancer immunotherapy and its optimal control

Cancer is one of the most serious illnesses in all of the world. Although most of the cancer patients are treated with chemotherapy, radiotherapy and surgery, wide research is conducted related to experimental and theoretical immunology. In recent years, the research on cancer immunotherapy has led to major medical advances. Cancer immunotherapy refers to the stimulation of immune system to deal with cancer cells. In medical practice, it is mainly achieved by using effector cells such as activated T-cells and Interleukin-2 (IL-2), which is the main cytokine responsible for lymphocyte activation, growth and differentiation. A well-known mathematical model, named as Kirschner-Panetta (KP) model, represents richly the dynamics of the interaction between cancer cells, IL-2 and the effector cells. The dynamics of the KP model is described and the solution to which is approximated by using polynomial approximation based methods such as Adomian decomposition method and differential transform method. The rich nonlinearity of the KP model causes these approaches to become so complicated in order to deal with the representation of polynomial approximations. It is illustrated that the approximated polynomials are in good agreement with the solution obtained by common numerical approaches. In the KP model, the growth of the tumour cells can be expressed by a linear function or any limited-growth function such as logistic equation, in which the cancer population possesses an upper bound mentioned as carrying capacity. Effector cells and IL-2 construct two external sources of medical treatment to stimulate immune system to eradicate cancer cells. Since the main goal in immunotherapy is to remove the tumour cells with the least probable medication side effects, an advanced version of the model may include a time dependent external sources of medical treatment, meaning that the external sources of medical treatment could be considered as control functions of time and therefore the optimum use of medical sources can be evaluated in order to achieve the optimal measure of an objective function. With this sense of direction, two distinct strategies are explored. The first one is to only consider the external source of effector cells as the control function to formulate an optimal control problem. It is shown under which circumstances, the tumour is eliminated. The approach in the formulation of the optimal control is the Pontryagin maximum principal. Furthermore the optimal control problem will be dealt with using particle swarm optimization (PSO). It is shown that the obtained results are significantly better than those obtained by previous researchers. The second strategy is to formulate an optimal control problem by considering both the two external sources as the controls. To our knowledge, it is the first time to present a multiple therapeutic protocol for the KP model. Some MATLAB routines are develop to solve the optimal control problems based on Pontryagin maximum principal and also the PSO. As known, fractional differential equations are more appropriate to describe the persistent memory of physical phenomena. Thus, the fractional KP model is defined in the sense of Caputo differentiation operator. An effective method for numerical treatment of the model is described, namely Predictor-Corrector method of Adams-Bashforth-Moulton type. A robust MATLAB routine is coded based on the mentioned approach and the solution obtained will be compared with those of the classical KP model. The code is prepared in such a way to be able to deal with systems of fractional differential equations, in which each equation has its own fractional order (i.e. multi-order systems of fractional differential equations). The theorems for existence of solutions and the stability analysis of the fractional KP model are represented. In this regard, a frequently used method of solving fractional differential equations (FDEs) is described in details, namely multi-step generalized differential transform method (MSGDTM), then it is illustrated that the method neglects the persistent memory property and takes the incorrect approach in dealing with numerical solutions of FDEs and therefore it is unfit to be used in differential equations governed by fractional differentiation operators. The sigmoidal behavior of the solution to the logistic equation caused it to be one of the most versatile models in natural sciences and therefore the fractional logistic equation would be a relevant problem to be dealt with. Thus, a power series of Mittag-Leffer functions is introduced, the behaviour of which is in good agreement with the solution to fractional logistic equation (FLE), and then a fractional integro-differential equation is represented and proved to be satisfied with the power series of Mittag-Leffler function. The obtained fractional integro-differential equation is named as modified fractional differential equation (MFDL) and possesses a nonlinear additive term related to the solution of the logistic equation (LE). The method utilized in the thesis, may be appropriately applied to the analysis of solutions to nonlinear fractional differential equations of mathematical physics. Inverse problems to FDEs occur in many branches of science. Such problems have been investigated, for instance, in fractional diffusion equation and inverse boundary value problem for semi- linear fractional telegraph equation. The determination of the order of fractional differential equations is an issue, which has been analyzed and discussed in, for instance, fractional diffusion equations. Thus, fractional order estimation has been conducted for some classes of linear fractional differential equations, by introducing the relationship between the fractional order and the asymptotic behaviour of the solutions to linear fractional differential equations. Fractional optimal control problems, in which the system and (or) the objective function are described based on fractional derivatives, are much more complicated to be solved by using a robust and reliable numerical approach. Thus, a MATLAB routine is provided to solve the optimal control for fractional KP model and the obtained solutions are compared with those of classical KP model. It is shown that the results for fractional optimal control problems are better than classical optimal control problem in the sense of the amount of drug administration

The adomian decomposition method applied to blood flow through arteries in the presence of a magnetic field

A dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, in fulfilment of requirements for the degree of Master of Science. February 16, 2015.The Adomian decomposition method is an effective procedure for the analytical solution of a wide class of dynamical systems without linearization or weak nonlinearity assumptions, closure approximations, perturbation theory, or restrictive assumptions on stochasticity. Our aim here is to apply the Adomian decomposition method to steady two-dimensional blood flow through a constricted artery in the presence of a uniform transverse magnetic field. Blood flow is the study of measuring blood pressure and determining flow through arteries. Blood flow is assumed to be Newtonian and is governed by the equation of continuity and the momentum balanced equation (which are known as the Navier-Stokes equations). This model is consistent with the principles of ferro-hydrodynamics and magnetohydrodynamics and takes into account both magnetization and electrical conductivity of blood. We apply the Adomian decomposition method to the equations governing blood flow through arteries in the presence of an external transverse magnetic field. The results show that the e ect of a uniform external transverse magnetic field applied to blood flow through arteries favors the physiological condition of blood. The motion of blood in stenosed arteries can be regulated by applying a magnetic field externally and increasing/decreasing the intensity of the applied field

A finite volume approach for the numerical analysis and solution of the Buckley-Leverett equation including capillary pressure

The study of petroleum recovery is significant for reservoir engineers. Mathematical models of the immiscible displacement process contain various assumptions and parameters, resulting in nonlinear governing equations which are tough to solve. The Buckley-Leverett equation is one such model, where controlling forces like gravity and capillary forces directly act on saturation profiles. These saturation profiles have important features during oil recovery. In this thesis, the Buckley-Leverett equation is solved through a finite volume scheme, and capillary forces are considered during this calculation. The detailed derivation and calculation are also illustrated here. First, the method of characteristics is used to calculate the shock speed and characteristics curve behaviour of the Buckley-Leverett equation without capillary forces. After that, the local Lax-Friedrichs finite-volume scheme is applied to the governing equation (assuming there are no capillary and gravity forces). This mathematical formulation is used for the next calculation, where the cell-centred finite volume scheme is applied to the Buckley- Leverett equation including capillary forces. All calculations are performed in MATLAB. The fidelity is also checked when the finite-volume scheme is computed in the case where an analytical solution is known. Without capillary pressure, all numerical solutions are calculated using explicit methods and smaller time steps are used for stability. Later, the fixed-point iteration method is followed to enable the stability of the local Lax-Friedrichs and Cell-centred finite volume schemes using an implicit formulation. Here, we capture the number of iterations per time-steps (including maximum and average iterations per time-step) to get the solution of water saturation for a new time-step and obtain the saturation profile. The cumulative oil production is calculated for this study and illustrates capillary effects. The influence of viscosity ratio and permeability in capillary effects is also tested in this study. Finally, we run a case study with valid field data and check every calculation to highlight that our proposed numerical schemes can capture capillary pressure effects by generating shock waves and providing single-valued saturation at each position. These saturation profiles help find the amount of water needed in an injection well to displace oil through a production well and obtains good recovery using the water flooding technique

A new hybrid method for solving nonlinear fractional differential equations

In this paper, numerical solution of initial and boundary value problems for nonlinear fractional differential equations is considered by pseudospectral method. In order to avoid solving systems of nonlinear equations resulting from the method, the residual function of the problem is constructed, as well as a suggested unconstrained optimization model solved by PSOGSA algorithm. Furthermore, the research inspects and discusses the spectral accuracy of Chebyshev polynomials in the approximation theory. The following scheme is tested for a number of prominent examples, and the obtained results demonstrate the accuracy and efficiency of the proposed method

Study of reactor constitutive model and analysis of nuclear reactor kinetics by fractional calculus approach

The diffusion theory model of neutron transport plays a crucial role in reactor theory since it is simple enough to allow scientific insight, and it is sufficiently realistic to study many important design problems. The neutrons are here characterized by a single energy or speed, and the model allows preliminary design estimates. The mathematical methods used to analyze such a model are the same as those applied in more sophisticated methods such as multi-group diffusion theory, and transport theory. The neutron diffusion and point kinetic equations are most vital models of nuclear engineering which are included to countless studies and applications under neutron dynamics. By the help of neutron diffusion concept, we understand the complex behavior of average neutron motion. The simplest group diffusion problems involve only, one group of neutrons, which for simplicity, are assumed to be all thermal neutrons. A more accurate procedure, particularly for thermal reactors, is to split the neutrons into two groups; in which case thermal neutrons are included in one group called the thermal or slow group and all the other are included in fast group. The neutrons within each group are lumped together and their diffusion, scattering, absorption and other interactions are described in terms of suitably average diffusion coefficients and cross-sections, which are collectively known as group constants. We have applied Variational Iteration Method and Modified Decomposition Method to obtain the analytical approximate solution of the Neutron Diffusion Equation with fixed source. The analytical methods like Homotopy Analysis Method and Adomian Decomposition Method have been used to obtain the analytical approximate solutions of neutron diffusion equation for both finite cylinders and bare hemisphere. In addition to these, the boundary conditions like zero flux as well as extrapolated boundary conditions are investigated. The explicit solution for critical radius and flux distributions are also calculated. The solution obtained in explicit form which is suitable for computer programming and other purposes such as analysis of flux distribution in a square critical reactor. The Homotopy Analysis Method is a very powerful and efficient technique which yields analytical solutions. With the help of this method we can solve many functional equations such as ordinary, partial differential equations, integral equations and so many other equations. It does not require enough memory space in computer, free from rounding off errors and discretization of space variables. By using the excellence of these methods, we obtained the solutions which have been shown graphically

Periodic Solution of Nonlinear Conservative Systems

Conservative systems represent a large number of naturally occurring and artificially designed scientific and engineering systems. A key consideration in the theory and application of nonlinear conservative systems is the solution of the governing nonlinear ordinary differential equation. This chapter surveys the recent approximate analytical schemes for the periodic solution of nonlinear conservative systems and presents a recently proposed approximate analytical algorithm called continuous piecewise linearization method (CPLM). The advantage of the CPLM over other analytical schemes is that it combines simplicity and accuracy for strong nonlinear and large-amplitude oscillations irrespective of the complexity of the nonlinear restoring force. Hence, CPLM solutions for typical nonlinear Hamiltonian systems are presented and discussed. Also, the CPLM solution for an example of a non-Hamiltonian conservative oscillator was presented. The chapter is aimed at showcasing the potential and benefits of the CPLM as a reliable and easily implementable scheme for the periodic solution of conservative systems

On Bernstein Polynomials Method to the System of Abel Integral Equations

This paper deals with a new implementation of the Bernstein polynomials method to the numerical solution of a special kind of singular system. For this aim, first the truncated Bernstein series polynomials of the solution functions are substituted in the given problem. Using some properties of these polynomials, the solution of the problem is reduced to solve a linear system of algebraic equations. In order to confirm the reliability and accuracy of the proposed method, some weakly Abel integral equations systems with comparisons are solved in detail as numerical examples