Volume 7, Issue 3

Mathematical model for an adiabatic tubular chemical reactor which processes an irreversible exothermic chemical reaction has been considered. For steady state solution for an adiabatic tubular chemical reactor, the model can be reduced to ordinary differential equation with a parameter in the boundary conditions. Again the ordinary differential equation has been converted into a Hammerstein integral equation which can be solved numerically. B-spline wavelet method has been developed to approximate the solution of Hammerstein integral equation. This method reduces the integral equation to a system of algebraic equations. The numerical results obtained by the present method have been compared with the available results.

We introduce new numerical techniques for solving nonlinear unsteady Burgers equation. The numerical technique involves discretization of all variables except the time variable which converts nonlinear PDE into nonlinear ODE system. Stability of the nonlinear system is verified using Lyapunov’s stability criteria. Implicit stiff solvers backward differentiation formula of order one, two and three are used to solve the nonlinear ODE system. Four test problems are included to show the applicability of introduced numerical techniques. Numerical solutions so obtained are compared with solutions of existing schemes in literature. The proposed numerical schemes are found to be simple, accurate, fast, practical and superior to some existing methods.

In this study, the peristaltic flow of a Casson fluid in a channel is considered in the presence of an applied magnetic field. Flow is considered in the moving frame of reference with constant velocity along the wave. The developed mathematical model is presented by a set of partial differential equations. A numerical algorithm based on finite element method is implemented to evaluate the numerical solution of the governing partial differential equations in the stream-vorticity formulation. The obtained results are independent of low Reynolds number and long wavelength assumptions, so the effects of non-zero moderate Reynolds number are presented. The expression for the pressure is also calculated implicitly and discussed through graphs. Computed solutions are presented in the form of contours of streamlines and vorticity. Velocity profile and pressure distribution for variation of different involved parameters are also presented through graphs. The investigation shows that the strength of circulation for stream function increases by increasing the Reynolds and Hartmann numbers. Enhancement in longitudinal velocity is noted by increasing the Reynolds number and Casson parameter while increasing Hartmann number reduces the longitudinal velocity. Comparison of the present results with the available results in literature is also included and found in good agreement.

The laminar boundary layer MHD three-dimensional mixed convective flow of Maxwell nanofluid towards a bidirectional stretching sheet with non-linear radiation is analyzed. A constant magnetic field is implemented normal to the fluid flow direction. A numerical technique of Runge-Kutta-Fehlberg (RFK45) is utilized to obtain the numerical solution of the dimensionless coupled ODEs with associated boundary conditions. The various pertinent dimensionless parameters on the flow are examined with the help of graphs and tables. Results shows that, nonlinear thermal radiation is more influential o on temperature profile when compared to linear thermal radiation.

In quest to contain and subsequently eradication Human Immunodeficiency virus (HIV) in the society, mathematical modelling remains an important research tool. In this paper, we formulated a mathematical model to study the effects of cortisol on immune response to HIV capturing the roles played by dendritic cells, T helper cells, regulatory T cells and cytotoxic T cells in the virus replication dynamics. The primary source of concentration of cortisol in this work is through psychological stress. Numerical experiments are performed to examine the effect of cortisol on selective inhibition of antigen presentation activities and up-regulation of naive cytotoxic T cells activation in the case of acute and persistent stressful conditions.

In present work, nonlinear fractional partial differential equations namely transport equation and Fokker-Planck equation involving local fractional differential operators, are investigated by means of the local fractional homotopy perturbation Sumudu transform method. The proposed method is a coupling of homotopy perturbation method with local fractional Sumudu transform and is used to describe the non-differentiable problems. Numerical simulation results are projected to show the efficiency of the proposed technique.

In the present paper, the flow of an incompressible, electrically conducting dusty fluid over a stretching sheet is considered. The Cattaneo- Christov heat flux theory is employed to control the thermal boundary layer. The flow equations are transformed into nonlinear ordinary differential equations (NODEs) and which are solved with help of Runge-Kutta 4 th order method. Flow equations are examined with respect to boundary conditions namely prescribed wall temperature (PWT) and prescribed heat flux (PHF) cases. In general PWT and PHF boundary conditions are very useful in the industrial as well as manufacturing up and down processes. Impact of the emerging parameters on the dimensionless velocity and temperature as well as friction coefficient and local Nusselt number are examined. We also validated my results with already available literature. It is found that the heat transfer rate of the flow in PWT case is higher than that of PHF case. These results can help us to conclude that for higher heating processes (Heating industries) PWT case and lesser heating processes (Cooling industries) PHF boundary condition is useful.

In present study, some sufficient conditions for the exponential stability of impulsive delay differential equations are obtained by introducing weight function in the norm and applying the concept of Lyapunov functions and Razumikhin techniques. The function ψ plays the role of weight and hence increases the rate of convergence towards stability. The obtained results are demonstrated with examples.