Search

Researchers are using different types of nanomaterials for the enhancement of the thermal performance of regular fluids such as water, kerosene oil, etc. However, these days, the researchers are more interested in hybrid nanomaterials. The purpose of this communication is to examine the stability analysis of Cu–Al2O3/water hybrid nanofluid over a non-linear shrinking sheet. The hybrid nanomaterials are composed of Cu and Al2O3. These hybridized nanomaterials are then dissolved in water taken as base fluid to form Cu–Al2O3//water hybrid nanofluid. Mathematical analysis and modeling have been attended in the presence of viscous dissipation and suction/injection effects. The governing equations of mathematical models are transformed into self-similar solutions in the form of ODEs by using similarity transformation. Solutions of the non-linear ODEs are created by employing of three-stage Lobatto IIIa formula which is built-in BVP4C function in the MATLAB software. A comparison of the current study has been done with the preceding published literature. The distributions of velocity, temperature profiles, coefficient of skin friction and heat transfer rate are presented graphically and conferred for numerous significant parameters entering into the problem. Results revealed the existence of dual solutions for a certain range of the suction/blowing parameter. Stability analysis is also done in order to obtain dual solutions stability. The smallest eigenvalues suggest that the first solution is stable from the second solution. Hybrid nanomaterials have a high scope toward nurturing our day-to-day life. Keywords: Hybrid nanomaterials, Nanofluid, Dual solutions, Stability analysis, Viscous dissipation, Shrinking surface

Block methods for the solution of fourth order ordinary differential equations are conventionally developed with specific grid or off-grid points. An implicit five-step block method with generalised equidistant points for solving fourth order initial value problems (IVPs) is presented in this article. The generalised points allow for flexibility in choice of points and the strategy adopted for the developing the block method considers the general block form with integration of Taylor series expansion. This results in a family of schemes implemented as simultaneous integrators of the fourth order IVPs at all grid points. The basic properties of the block method are investigated and the block method is seen to satisfy the property of convergence which is evident in the numerical solutions. In addition, the numerical results displayed the five-step block method performing better than existing methods. The block method is also applied to a ship dynamics fourth order model. Keywords: Block method, Generalised points, Fourth order, Linear, Non-linear, Initial value problems

his article presents a new generalized algorithm for developing k-step (m+1)^th derivative block methods for solving m^th order ordinary differential equations. This new algorithm utilizes the concept from the conventional Taylor series approach of developing linear multistep methods. Certain examples are given to show the simplicity involved in the usage of this new generalized algorithm.

The existing hybrid methods for solving ordinary differential equations were only derived using specific off-step points. Thus, this paper proposes a new four-step block method with three generalized off-step points for solving initial value problems of fourth order ordinary differential equations directly. The strategy employed to develop this method is interpolating the basis function at yn+j,j=0(1)3 and collocating the fourth derivative of the basis function at all points within the selected interval. The implementation of this method in a block-by-block fashion can overcome the setbacks of applying starting values and predictors which are created in predictor-corrector approach. The convergence analysis of the developed method is performed and the accuracy of the method is tested on several problems. The numerical results indicate that the new method outperforms the existing ones in terms of errors. In addition, the new method does not require much computation when compared with previous methods.

In this paper, we explore dual solutions of MHD flow, heat and mass transfer of micropolar nanofluid over a linear vertical shrinking surface with buoyancy effects, which was not considered in the previous works. The governing fluid flow equations of this problem are transformed into nonlinear boundary value problems (BVPs) of ordinary differential equations (ODEs) by applying similarity variables. The resultant BVPs are converted into initial value problems (IVPs) by using shooting method which then resolved by employing Runge Kutta of order four. The impacts of the governing parameters, such as suction parameter, material parameter, Richardson number, magnetic parameter, Prandtl number, thermophoresis and Brownian motion parameters on velocity, angular velocity, temperature, and concentration are illustrated graphically. The results indicate that the existence of a range of dual solutions and no-solutions. When Richardson number (δ) is increased, the reduction of the velocity of micropolar nanofluid has occurred in the second solution. The stability analysis on dual solutions, however, reveals that only the first solution is stable.

This study investigates the numerical solutions of MHD boundary layer and heat transfer of the Williamson fluid flow on the exponentially vertical shrinking sheet, having variable thickness and thermal conductivity under effects of the velocity and thermal slip parameters. It is also assumed that shrinking/stretching velocity, as well as the wall temperature, has the exponential function form. In this study, the continuity, momentum and energy equations with buoyancy parameter and Hartmann number are incorporated especially in the Williamson fluid flow case. Similarity transformation variables have been employed to formulate the ordinary differential equations (ODEs) from partial differential equations (PDEs). The resultant ODEs are solved by shooting method with Runge Kutta of fourth order method in Maple software. The effects of the different applied non-dimensional physical parameters on the boundary layer and heat transfer flow problems are presented in graphs. The effects of Williamson parameter, Prandtl number, and slip parameters on velocity and temperature profiles have been thoroughly demonstrated and discussed. The numerical results show that the buoyancy force and the slip parameters contribute to the occurrence of the dual solutions on the boundary layer and heat transfer flow problems. Furthermore, the stability analysis suggests that the first solution is stable and physically possible.

This article aims to study numerically the rotating, steady, and three-dimensional (3D) flow of a hybrid nanofluid over an exponentially shrinking sheet with the suction effect. We considered water as base fluid and alumina (Al2O3), and copper (Cu) as solid nanoparticles. The system of governing partial differential equations (PDEs) was transformed by an exponential similarity variable into the equivalent system of ordinary differential equations (ODEs). By applying a three-stage Labatto III-A method that is available in bvp4c solver in the Matlab software, the resultant system of ODEs was solved numerically. In the case of the hybrid nanofluid, the heat transfer rate improves relative to the viscous fluid and regular nanofluid. Two branches were obtained in certain ranges of the involved parameters. The results of the stability analysis revealed that the upper branch is stable. Moreover, the results also indicated that the equations of the hybrid nanofluid have a symmetrical solution for different values of the rotation parameter (Ω).

A new numerical method for solving fourth order ordinary differential equations directly is proposed in this paper. Interpolation and collocation were employed in developing this method using seven steps. The use of the approximated power series as an interpolation equation was adopted in deriving the method. The basic properties of the new method such as zero-stability, consistency, convergence and order are established. The numerical results indicate that the new method gives better accuracy than the existing methods when it is applied to fourth ordinary differential equations.

In this analysis, we aim to examine the heat transfer and flow characteristics of a copper-aluminum/water hybrid nanofluid in the presence of viscous dissipation, magnetohydrodynamic (MHD), and porous medium effect over the shrinking sheet. The governing equations of the fluid model have been acquired by employment of the model of Tiwari and Das, with additional properties of the hybrid nanofluid. The system of partial differential equations (PDEs) has been converted into ordinary differential equations (ODEs) by adopting the exponential similarity transformation. Similarity transformation is an essential class of phenomenon where the symmetry of the scale helps to reduce the number of independent variables. Note that ODE solutions demonstrate the PDEs symmetrical behavior for the velocity and temperature profiles. With BVP4C solver in the MATLAB program, the system of resulting equations has been solved. We have compared the present results with the published results and found in excellent agreements. The findings of the analysis are also displayed and discussed in depth graphically and numerically. It is discovered that two solutions occur in definite ranges of suction and magnetic parameters. Dual (no) similarity solutions can be found in the range of Sc≤S and Mc≤M (Sc>S and Mc>M). By performing stability analysis, the smallest values of eigenvalue are obtained, suggesting that a stable solution is the first one. Furthermore, the graph of the smallest eigenvalue shows symmetrical behavior. By enhancing the Eckert number values the temperature of the fluid is raised.

In this research, we intend to develop a dynamical system for the magnetohydrodynamic (MHD) flow of an electrically conducting Casson nanofluid on exponentially shrinking and stretching surfaces, in the presence of a velocity and concertation slip effect, with convective boundary conditions. There are three main objectives of this article, specifically, to discuss the heat characteristics of flow, to find multiple solutions on both surfaces, and to do stability analyses. The main equations of flow are governed by the Brownian motion, the Prandtl number, and the thermophoresis parameters, the Schmid and Biot numbers. The shooting method and three-stage Lobatto IIIa formula have been employed to solve the equations. The ranges of the dual solutions are f w c 1 ≤ f w and λ c ≤ λ , while the no solution ranges are f w c 1 > f w and λ c > λ . The results reveal that the temperature of the fluid increases with the extended values of the thermophoresis parameter, the Brownian motion parameter, and the Hartmann and Biot numbers, for both solutions. The presence of dual solutions depends on the suction parameter. In order to indicate that the first solution is physically relevant and stable, a stability analysis has been performed.

In the current study, the characteristics of heat transfer of a steady, two-dimensional, stagnation point, and magnetohydrodynamic (MHD) flow of shear thickening Casson fluid on an exponentially vertical shrinking/stretching surface are examined in attendance of convective boundary conditions. The impact of the suction parameter is also considered. The system of governing partial differential equations (PDEs) and boundary conditions is converted into ordinary differential equations (ODEs) with the suitable exponential similarity variables of transformations and then solved using the shooting method with the fourth order Runge–Kutta method. Similarity transformation is an important class of phenomena in which scale symmetry allows one to reduce the number of independent variables of the problem. It should be noted that solutions of the ODEs show the symmetrical behavior of the PDES for the profiles of velocity and temperature. Similarity solutions are obtained for the case of stretching/shrinking and suction parameters. It is revealed that there exist two ranges of the solutions in the specific ranges of the physical parameters, three solutions depend on the opposing flow case where stagnation point (A) should be equal to 0.1, two solutions exist when λ1 = 0 where λ1 is a mixed convection parameter and A > 0.1, and a single solution exists when λ1 > 0. Moreover, the effects of numerous applied parameters on velocity, temperature distributions, skin friction, and local Nusselt number are examined and given through tables and graphs for both shrinking and stretching surfaces.

In this article, we reconsidered the problem of Aurangzaib et al., and reproduced the results for triple solutions. The system of governing equations has been transformed into the system of non-linear ordinary differential equations (ODEs) by using exponential similarity transformation. The system of ODEs is reduced to initial value problems (IVPs) by employing the shooting method before solving IVPs by the Runge Kutta method. The results reveal that there are ranges of multiple solutions, triple solutions, and a single solution. However, Aurangzaib et al., only found dual solutions. The effect of the micropolar parameter, suction parameter, and Prandtl number on velocity, angular velocity, and temperature profiles have been taken into account. Stability analysis of triple solutions is performed and found that a physically possible stable solution is the first one, while all leftover solutions are not stable and cannot be experimentally seen.

In this article, the magnetohydrodynamic (MHD) flow of Casson nanofluid with thermal radiation over an unsteady shrinking surface is investigated. The equation of momentum is derived from the Navier−Stokes model for non-Newtonian fluid where components of the viscous terms are symmetric. The effect of Stefan blowing with partial slip conditions of velocity, concentration, and temperature on the velocity, concentration, and temperature distributions is also taken into account. The modeled equations of partial differential equations (PDEs) are transformed into the equivalent boundary value problems (BVPs) of ordinary differential equations (ODEs) by employing similarity transformations. These similarity transformations can be obtained by using symmetry analysis. The resultant BVPs are reduced into initial value problems (IVPs) by using the shooting method and then solved by using the fourth-order Runge−Kutta (RK) technique. The numerical results reveal that dual solutions exist in some ranges of different physical parameters such as unsteadiness and suction/injection parameters. The thickness of the velocity boundary layer is enhanced in the second solution by increasing the magnetic and velocity slip factor effect in the boundary layer. Increment in the Prandtl number and Brownian motion parameter is caused by a reduction of the thickness of the thermal boundary layer and temperature. Moreover, stability analysis performed by employing the three-stage Lobatto IIIA formula in the BVP4C solver with the help of MATLAB software reveals that only the first solution is stable and physically realizable.

In this paper, the unsteady magnetohydrodynamic (MHD) flow of hybrid nanofluid (HNF) composed of C u − A l 2 O 3 /water in the presence of a thermal radiation effect over the stretching/shrinking sheet is investigated. Using similarity transformation, the governing partial differential equations (PDEs) are transformed into a system of ordinary differential equations (ODEs), which are then solved by using a shooting method. In order to validate the obtained numerical results, the comparison of the results with the published literature is made numerically as well as graphically and is found in good agreements. In addition, the effects of many emerging physical governing parameters on the profiles of velocity, temperature, skin friction coefficient, and heat transfer rate are demonstrated graphically and are elucidated theoretically. Based on the numerical results, dual solutions exist in a specific range of magnetic, suction, and unsteadiness parameters. It was also found that the values of f ″ ( 0 ) rise in the first solution and reduce in the second solution when the solid volume fraction ϕ C u is increased. Finally, the temporal stability analysis of the solutions is conducted, and it is concluded that only the first solution is stable.

A numerical study was carried out to examine the magnetohydrodynamic (MHD) flow of micropolar fluid on a shrinking surface in the presence of both Joule heating and viscous dissipation effects. The governing system of non-linear ordinary differential equations (ODEs) was obtained from the system of partial differential equations (PDEs) by employing exponential transformations. The resultant equations were transformed into initial value problems (IVPs) by shooting technique and then solved by the Runge−Kutta (RK) method. The effects of different parameters on velocity, angular velocity, temperature profiles, skin friction coefficient, and Nusselt number were obtained and demonstrated graphically. We observed that multiple solutions occurred in certain assortments of the parameters for suction on a surface. The stability analysis of solutions was performed, and we noted that the first solution was stable while the remaining two solutions were not. The results also showed that the velocity of the fluid increased as the non-Newtonian parameter rose in all solutions. Furthermore, it was detected that the temperature of fluid rose at higher values of the Eckert number in all solutions.

The present study accentuates the heat transfer characteristics of a convective condition of micropolar nanofluid on a permeable shrinking/stretching inclined surface. Brownian and thermophoresis effects are also involved to incorporate energy and concentration equations. Moreover, linear similarity transformation has been used to transform the system of governing partial differential equations (PDEs) into a set of nonlinear ordinary differential equations (ODEs). The numerical comparison has been done with the previously published results and found in good agreement graphically and tabular form by using the shooting method in MAPLE software. Dual solutions have been found in the specific range of shrinking/stretching surface parameters and the mass suction parameter for the opposing flow case. Moreover, the skin friction coefficient, the heat transfer coefficient, the couple stress coefficient, and the concentration transfer rate decelerate in both solutions against the mass suction parameter for the augmentation of the micropolar parameter respectively. The first (second) solution is the stable (unstable) solution and can (not) be considered as a real solution as the values of the smallest eigenvalues are positive (negative).

A new one-step block method with generalized three hybrid points for solving initial value problems of second-order ordinary differential equations directly is proposed. In deriving this method, a power series approximate function is interpolated at {xn,xn+r} while its second and third derivatives are collocated at all points {xn,xn+r,xn+s,xn+t,xn+1} in the given interval. The proposed method is then tested on initial value problems of second-order ordinary differential equations solved by other methods previously. The numerical results confirm the superiority of the new method to the existing methods in terms of accuracy.

Nonlinear boundary value problems (BVPs) are more tedious to solve than their linear counterparts. This is observed in the extra computation required when determining the missing conditions in transforming BVPs to initial value problems. Although a number of numerical approaches are already existent in literature to solve nonlinear BVPs, this article presents a new block method with improved accuracy to solve nonlinear BVPs. A m+1th-step block method is developed using a modified Taylor series approach to directly solve fourth-order nonlinear boundary value problems (BVPs) where m is the order of the differential equation under consideration. The schemes obtained were combined to simultaneously produce solution to the fourth-order nonlinear BVPs at m+1 points iteratively. The derived block method showed improved accuracy in comparison to previously existing authors when solving the same problems. In addition, the suitability of the m+1th-step block method was displayed in the solution for magnetohydrodynamic squeezing flow in porous medium.

A general one-step three-hybrid (off-step) points block method is proposed for solving fourth-order initial value problems of ordinary differential equations directly. A power series approximate function is employed for deriving this method. The approximate function is interpolated at xn,xn+r,xn+s,xn+t while its fourth and fifth derivatives are collocated at all points xi, i=0,r,s,t,1, in the interval of approximation. Several fourth-order initial value problems of ordinary differential equations are then solved to compare the performance of the proposed method with the derived methods. The analysis of the method reveals that the method is consistent and zero stable concluding that the method is also convergent. The numerical results demonstrate the superiority of the new method over the existing ones in terms of error.

In this paper, steady two-dimensional laminar incompressible magnetohydrodynamic flow over an exponentially shrinking sheet with the effects of slip conditions and viscous dissipation is examined. An extended Darcy Forchheimer model was considered to observe the porous medium embedded in a non-Newtonian-Casson-type nanofluid. The governing equations were converted into nonlinear ordinary differential equations using an exponential similarity transformation. The resultant equations for the boundary values problem (BVPs) were reduced to initial values problems (IVPs) and then shooting and Fourth Order Runge-Kutta method (RK-4th method) were applied to obtain numerical solutions. The results reveal that multiple solutions occur only for the high suction case. The results of the stability analysis showed that the first (second) solution is physically reliable (unreliable) and stable (unstable).

A numerical investigation is made to determine the occurrence of the multiple solutions of MHD Casson fluid in a porous channel. Governing partial differential equation of the proposed problem converted into nonlinear ordinary differential equations by using similarity transformation. Numerical technique known as shooting method is used to investigate the existence of the multiple solutions for the variations of different parameters. Effects of physical parameters on velocity profile, temperature, concentration, and skin friction are presented in pictorial and tabulation representation.

Ordinary differential equations are commonly used for mathematical modeling in many diverse fields such as engineering, industrial mathematics, operation research, artificial intelligence, management, sociology and behavioural sciences. Numerous problems encountered in these fields require lengthy computation and immediate solution. In this paper, a new method called parallel 3-point implicit block method for solving second order ODES is developed. This method takes full advantage of parallel computers because the numerical solution can be computed at three points simultaneously. As a result, the solution can be obtained faster if compared to the conventional methods where the numerical solution is computed at one point at a time. Computational advantages are presented comparing the results obtained by the new method with that of 1-point and 2-point implicit block methods. The numerical results show that parallel 3-point implicit block method reduces the total number of steps and execution time without sacrificing the accuracy.