A high-order uniform Cartesian grid compact finite difference scheme for the Goursat problem is developed. The basic idea of high-order compact schemes is to find the compact approximations to the derivatives terms by differentiating centrally the governing equations. Our compact scheme will approximate the derivative terms by involving the higher terms and reducing the number of grid points. The compact finite difference scheme is given for general form of the Goursat problem in uniform domain and illustrates the performance by applying a linear problem. Numerical experiments have been conducted with the new scheme and encouraging results have been obtained. In this paper we present the compact finite difference scheme for the Goursat problem. With the aid of computational software the scheme was programmed for determining the relative errors of linear Goursat problem.
In this paper, the optimal homotopy asymptotic method (OHAM) is applied to obtain an approximate solution of the nonlinear Riccati differential equation. The method is tested on several types of Riccati differential equations and comparisons that were made with numerical results showed the effectiveness and accuracy of this method.
We study and discuss the effect of thermal slip on steady free convection flow of a viscous, incompressible micropolar fluid past a vertical moving plate in a saturated porous medium. The effect of viscous dissipation is incorporated in the energy equation. The associated partial differential equations are transformed into a system of ordinary differential equations using similarity transformations generated by a group method and this system is then solved numerically. The effect of controlling parameters on the dimensionless velocity, angular velocity and temperature as well as friction factor, couple stress factor and heat transfer rate are shown graphically and discussed in detail. It is found that the dimensional velocity and angular velocity decrease whilst the temperature increases with velocity slip parameter. It is further found that thermal slip decreases the dimensional velocity and temperature but increases the dimensional angular velocity. Data from published work and our results are found to be in good agreement.
In this paper, the combined influences of biotic interactions, environmental components and harvesting strategy on the spread of Hantavirus are investigated. By employing a multi-species model consisting of (susceptible and infected) rodents and alien species, we show that interspecific competition from alien species has an effect in reducing the spread of infection, and this species could be employed as a potential biocontrol agent. Our analysis using numerical continuation and simulation also reveals the conditions under which Hantavirus infection occurs and disappears as the environmental conditions and the intensity of harvesting change. Without harvesting, infection emerges when environments are conducive. Inclusion of moderate harvesting in favourable environments can lead to disappearance of infection among rodent species. However, as the intensity of harvesting increases, this situation can cause extinction of all rodents species and consequently, jeopardise biodiversity. Overall, our results demonstrate how the interplay of different factors can combine to determine the spread of infectious diseases.
Homotopy continuation methods (HCMs) can be used to find the solutions
of polynomial equations. The advantages of HCMs over classical methods such as the
Newton and bisection methods are that HCMs are able to resolve divergence and starting
value problems. In this paper, we develop Super Ostrowski-HCM as a technique to
overcome the starting value problem. We compare the performance of this proposed
method with Ostrowski-HCM. The results provide evidence of the superiority of Super
Ostrowski-HCM.
Second order linear two-point boundary value problems were solved using extended cubic B-spline interpolation method. Extended cubic B-spline is an extension of cubic B-spline consisting of one shape parameter, called λ. The resulting approximated analytical solution for the problems would be a function of λ. Optimization of λ was carried out to find the best value of λ that generates the closest fit to the differential equations in the problems. This method approximated the solutions for the problems much more accurately compared to finite difference, finite element, finite volume and cubic B-spline interpolation methods.
Steady laminar mixed convection boundary layer flow past a horizontal circular cylinder with constant wall heat flux, immersed in a viscous and incompressible fluid of temperature-dependent viscosity is considered in this study. The governing partial differential equations were transformed using non-similar transformation and then solved numerically by an implicit finite-difference scheme known as the Keller-box method. The effects of temperature-dependent viscosity parameter θr on the flow and heat transfer characteristics were examined for various values of Prandtl number, Pr and the mixed convection parameter, λ. It was found that for both assisting and opposing flows, as θr increases, the local skin friction coefficient increases while the wall temperature decreases for air but for water, the local skin friction coefficient decreases then slightly increases while temperature decreases.
A combined similarity-numerical solution of the magnetohydrodynamic boundary layer slip flow of an electrically conducting non-Newtonian power-law nanofluid along a heated radiating moving vertical plate is explored. Our nanofluid model incorporates the influences of the thermophoresis and the Brownian motion. The basic transport equations are made dimensionless first and then suitable similarity transformations are applied to reduce them into a set of nonlinear ordinary differential equations with the associated boundary conditions. The reduced equations are then solved numerically. Graphical results for the non-dimensional flow velocity, the temperature and the nanoparticles volume fraction profiles as well as for the friction factor, the local Nusselt and the Sherwood numbers are exhibited and examined for various values of the controlling parameters to display the interesting aspects of the solutions. It was found that the friction factor increases with the increase of the magnetic field (M), whilst it is decreased with the linear momentum slip parameter (a). The linear momentum slip parameter (a) reduces the heat transfer rates and the nanoparticles volume fraction rates. Our results are compatible with the existing results for a special case.