Search

In this research, the three-dimensional (3D) steady and incompressible laminar Homann stagnation point nanofluid flow over a porous moving surface is addressed. The disturbance in the porous medium has been characterized by the Darcy-Forchheimer relation. The slip for viscous fluid is considered. The energy equation is organized in view of radiative heat flux which plays an important role in the heat transfer rate. The governing flow expressions are first altered into first-order ordinary ones and then solved numerically by the shooting method. Dual solutions are obtained for the velocity, skin friction coefficient, temperature, and Nusselt number subject to sundry flow parameters, magnetic parameter, Darcy-Forchheimer number, thermal radiation parameter, suction parameter, and dimensionless slip parameter. In this research, the main consideration is given to the engineering interest like skin friction coefficient (velocity gradient or surface drag force) and Nusselt number (temperature gradient or heat transfer rate) and discussed numerically through tables. In conclusion, it is noticed from the stability results that the upper branch solution (UBS) is more reliable and physically stable than the lower branch solution (LBS).

The present work concerns the momentum and heat transmission of the electro-magnetohydrodynamic (E-MHD) boundary layer Darcy-Forchheimer flow of a Sutterby fluid over a linear stretching sheet with slip. The nonlinear equations for the proposed model are analyzed numerically. Suitable techniques are used to transform the coupled nonlinear partial differential equations (PDEs) conforming to the forced balance law, energy, and concentration equations into a nonlinear coupled system of ordinary differential equations (ODEs). Numerical solutions of the transformed nonlinear system are obtained using a shooting method, improved by the Cash and Carp coefficients. The influence of important physical variables on the velocity, the temperature, the heat flux coefficient, and the skin-friction coefficient is verified and analyzed through graphs and tables. From the comprehensive analysis of the present work, it is concluded that by intensifying the magnitude of the Hartmann number, the momentum distribution decays, whereas the thermal profile of fluid increases. Furthermore, it is also shown that by aug- menting the values of the momentum slip parameter, the velocity profile diminishes. It is found that the Sutterby fluid model shows shear thickening and shear thinning behaviors. The momentum profile shows that the magnitude of velocity for the shear thickening case is dominant as compared with the shear thinning case. It is also demonstrated that the Sutterby fluid model reduces to a Newtonian model by fixing the fluid parameter to zero. In view of the limiting case, it is established that the surface drag in the case of the Sutterby model shows a trifling pattern as compared with the classical case.

The present article explores the entropy generation of radiating viscoelastic second grade nanofluid in a porous channel confined between two parallel plates. The boundaries of the plates are maintained at distinct temperatures and concentrations while the fluid is being sucked and injected periodically through upper and lower plates. The buoyancy forces, thermophoresis and Brownian motion are also considered due to the temperature and concentration differences across the channel. The system of governing partial differential equations has been transferred into a system of ordinary differential equations (ODEs) by appropriate similarity relations, and a shooting method with the fourth-order Runge-Kutta scheme is used for the solutions. The results are analyzed in detail for dimensionless velocity components. The temperature, concentration distributions, the entropy generation number, and the Bejan number corresponding to various fluid and geometric parameters are shown graphically. The skin friction, heat and mass transfer rates are presented in the form of tables. It is noticed that the temperature profile of the fluid is enhanced with the Brownian motion, whereas the concentration profile of the fluid is decreased with the thermophoresis parameter, and the entropy and Bejan numbers exhibit the opposite trend for the suction and injection ratio.

A mathematical model is proposed to execute the features of the non-uniform heat source or sink in the chemically reacting magnetohydrodynamic (MHD) Casson fluid across a slendering sheet in the presence of microorganisms and Cattaneo-Christov heat flux. Multiple slips (diffusion, thermal, and momentum slips) are applied in the modeling of the heat and mass transport processes. The Runge-Kutta based shooting method is used to find the solutions. Numerical simulation is carried out for various values of the physical constraints when the Casson index parameter is positive, negative, or infinite with the aid of plots. The coefficients of the skin factors, the local Nusselt number, and the Sherwood number are estimated for different parameters, and discussed for engineering interest. It is found that the gyrotactic microorganisms are greatly encouraged when the dimensionless parameters increase, especially when the Casson fluid parameter is negative. It is worth mentioning that the velocity profiles when the Casson fluid parameter is positive are higher than those when the Casson fluid parameter is negative or infinite, whereas the temperature and concentration fields show exactly opposite phenomena.

The large deflection of an axially extensible curved beam with a rectangular cross-section is investigated. The elastic beam is assumed to satisfy the Euler-Bernoulli postulation and be made of the Ludwick type material. Through reasonably simplified integration, the strain and curvature of the axis of the beam are presented in implicit formulations. The governing equations involving both geometric and material nonlinearities of the curved beam are derived and solved by the shooting method. When the initial curvature of the beam is zero, the curved beam is degenerated into a straight beam, and the predicted results obtained by the present model are consistent with those in the open literature. Numerical examples are further given for curved cantilever and simply supported beams, and the couplings between elongation and bending are found for the curved beams.

The capillary interactions, including the capillary force and capillary suction, between two unequal-sized particles with a fixed liquid volume are investigated. The capillary interaction model is used within the Young-Laplace framework. With the profile of the meridian of the liquid bridge, the capillary suction, and the liquid volume as state variables, the governing equations with two-fixed-point boundary are first derived using a variable substitution technique, in which the gravity effects are neglected. The capillary suction and geometry of the liquid bridge with a fixed volume are solved with a shooting method. In modeling the capillary force, the Gorge method is applied. The effects of various parameters including the distance between two particles, the ratio of particle radii, and the liquid-solid contact angles are discussed.

Free vibration response of functionally graded material (FGM) beams is studied based on the Levinson beam theory (LBT). Equations of motion of an FGM beam are derived by directly integrating the stress-form equations of elasticity along the beam depth with the inertial resultant forces related to the included coupling and higherorder shear strain. Assuming harmonic response, governing equations of the free vibration of the FGM beam are reduced to a standard system of second-order ordinary differential equations associated with boundary conditions in terms of shape functions related to axial and transverse displacements and the rotational angle. By a shooting method to solve the two-point boundary value problem of the three coupled ordinary differential equations, free vibration response of thick FGM beams is obtained numerically. Particularly, for a beam with simply supported edges, the natural frequency of an FGM Levinson beam is analytically derived in terms of the natural frequency of a corresponding homogenous Euler-Bernoulli beam. As the material properties are assumed to vary through the depth according to the power-law functions, the numerical results of frequencies are presented to examine the effects of the material gradient parameter, the length-to-depth ratio, and the boundary conditions on the vibration response.

The present article investigates the dual nature of the solution of the magneto-hydrodynamic (MHD) stagnation-point flow of a Prandtl fluid model towards a shrinking surface. The self-similar nonlinear ordinary differential equations are solved numerically by the shooting method. It is found that the dual solutions of the flow exist for certain values of the velocity ratio parameter. The special case of the first branch solutions (the classical Newtonian fluid model) is compared with the present numerical results of stretching flow. The results are found to be in good agreement. It is also shown that the boundary layer thickness for the second solution is thicker than that for the first solution.

This paper studies the thermal-diffusion and diffusion thermo-effects in the hydro-magnetic unsteady flow by a mixed convection boundary layer past an impermeable vertical stretching sheet in a porous medium in the presence of chemical reaction. The velocity of the stretching surface, the surface temperature, and the concentration are assumed to vary linearly with the distance along the surface. The governing partial differential equations are transformed into self-similar unsteady equations using similarity transformations and solved numerically by the Runge-Kutta fourth order scheme in association with the shooting method for the whole transient domain from the initial state to the final steady state flow. Numerical results for the velocity, the temperature, the concentration, the skin friction, and the Nusselt and Sherwood numbers are shown graphically for various flow parameters. The results reveal that there is a smooth transition of flow from unsteady state to the final steady state. A special case of our results is in good agreement with an earlier published work.

The magnetohydrodynamics (MHD) convection flow and heat transfer of an incompressible viscous nanofluid past a semi-infinite vertical stretching sheet in the presence of thermal stratification are examined. The partial differential equations governing the problem under consideration are transformed by a special form of the Lie symmetry group transformations, i.e., a one-parameter group of transformations into a system of ordinary differential equations which are numerically solved using the Runge-Kutta-Gill-based shooting method. It is concluded that the flow field, temperature, and nanoparticle volume fraction profiles are significantly influenced by the thermal stratification and the magnetic field.

Based on the assumption of finite deformation, the Hamilton variational principle is extended to a nonlinear elastic Euler-type beam-column structure located on a nonlinear elastic foundation. The corresponding three-dimensional (3D) mathematical model for analyzing the nonlinear mechanical behaviors of structures is established, in which the effects of the rotation inertia and the nonlinearity of material and geometry are considered. As an application, the nonlinear stability and the post-buckling for a linear elastic beam with the equal cross-section located on an elastic foundation are analyzed. One end of the beam is fully fixed, and the other end is partially fixed and subjected to an axial force. A new numerical technique is proposed to calculate the trivial solution, bifurcation points, and bifurcation solutions by the shooting method and the Newton-Raphson iterative method. The first and second bifurcation points and the corresponding bifurcation solutions are calculated successfully. The effects of the foundation resistances and the inertia moments on the bifurcation points are considered.

The problem of the steady magnetohydrodynamic (MHD) stagnation-point flow of an incompressible viscous fluid over a stretching sheet is studied. The effect of an induced magnetic field is taken into account. The nonlinear partial differential equations are transformed into ordinary differential equations via the similarity transformation. The transformed boundary layer equations are solved numerically using the shooting method. Numerical results are obtained for various magnetic parameters and Prandtl numbers. The effects of the induced magnetic field on the skin friction coefficient, the local Nusselt number, the velocity, and the temperature profiles are presented graphically and discussed in detail.

Combined heat and mass transfer on free, forced, and mixed convection flow along a porous wedge with magnetic effect in the presence of chemical reaction is investigated. The flow field characteristics are analyzed by the Runge-Kutta-Gill scheme with the shooting method as well as the local non-similarity method up to the third level of truncation, which are used to reduce the governing partial differential equations into nine ordinary differential equations. The governing boundary layer equations are converted to a dimensionless form by Falkner-Skan transformations. Because of the effect of suction/injection on the wall of the wedge with buoyancy force and variable wall temperature, the flow field is locally non-similar. Numerical calculations up to the third order level of truncation are carried out as a special case for different values of dimensionless parameters. Effects of the magnetic field strength in the presence of chemical reaction with variable wall temperature and concentration on the dimensionless velocity, temperature and concentration profiles are shown graphically.

This study is devoted to the investigation of thermal criticality for a reactive gravity driven thin film flow of a third-grade fluid with adiabatic free surface down an inclined isothermal plane. It is assumed that the reaction is exothermic under Arrhenius kinetics, neglecting the consumption of the material. The governing non-linear equations for conservation of momentum and energy are obtained and solved by using a new computational approach based on a special type of Hermite-Pade approximation technique implemented in MAPLE. This semi-numerical scheme offers some advantages over solutions obtained with traditional methods such as finite differences, spectral method, and shooting method. It reveals the analytical structure of the solution function. Important properties of overall flow structure including velocity field, temperature field, thermal criticality, and bifurcations are discussed.

The existence of solutions for the 2n-order m-point boundary value problem at resonance is obtained by using the coincidence degree theory of Mawhin. We give an example to demonstrate our result. The interest is that the nonlinear term may be noncontinuous.

The non-linear governing differential equations of immovably simply supported functionally graded material (FGM) rod subjected to thermal loads were derived. The thermal post-buckling behaviors of FGM rod made of ZrO2 and Ti-6A1-4V were analyzed by shooting method. Firstly, the thermal post-buckling equilibrium paths of the FGM rod with different gradient index in the uniform temperature field were plotted, and compared with the behaviors of the homogeneous rods made of ZrO2 and Ti-6A1-4V materials, respectively. For given value of end rotation angles, the influence of gradient index on the thermal post-buckling behaviors of FGM rod was discussed. Secondly, the thermal post-buckling characteristics of the FGM rod were analyzed when the temperature difference parameter is changed while the bottom temperature parameter remains constant, and when the bottom temperature parameter is changed while the temperature difference parameter remains constant, and compared with the characteristics of the two homogeneous material rods.

Analysis of thermal post-buckling of FGM (Functionally Graded Material) Timoshenko beams subjected to transversely non-uniform temperature rise is presented. By accurately considering the axial extension and transverse shear deformation in the sense of theory of Timoshenko beam, geometrical nonlinear governing equations including seven basic unknown functions for functionally graded beams subjected to mechanical and thermal loads were formulated. In the analysis, it was assumed that the material properties of the beam vary continuously as a power function of the thickness coordinate. By using a shooting method, the obtained nonlinear boundary value problem was numerically solved and thermal buckling and post-buckling response of transversely non-uniformly heated FGM Timoshenko beams with fixed-fixed edges were obtained. Characteristic curves of the buckling deformation of the beam varying with thermal load and the power law index are plotted. The effects of material gradient property on the buckling deformation and critical temperature of beam were discussed in details. The results show that there exists the tension-bend coupling deformation in the uniformly heated beam because of the transversely non-uniform characteristic of materials.

The well-posedness of the initial value problem of the Euler equations was mainly discussed based on the stratification theory, and the necessary and sufficient conditions of well-posedness are presented for some representative initial or boundary value problem, thus the structure of solution space for local (exact) solution of the Euler equations is determined. Moreover the computation formulas of the analytical solution of the well-posed problem are also given.

The influence of labyrinth seal on the stability of unbalanced rotor system was presented. Under the periodic excitation of rotor unbalance, the whirling vibration of rotor is synchronous if the rotation speed is below stability threshold, whereas the vibration becomes severe and asynchronous which is defined as unstable if the rotation speed exceeds threshold. The Muszynska model of seal force and shooting method were used to investigate synchronous solution of the dynamic equation of rotor system. Then, based on Floquet theory the stability of synchronous solution and unstable dynamic characteristic of system were analyzed.

Based on the geometrically nonlinear theory of axially extensible elastic rods, the governing equations of post-buckling of a clamped-free rod with variable cross-sections, subjected to a combined load, a concentrated axial load P at the free end and a non-uniformly distributed axial load q, are established. By using shooting method, the strong nonlinear boundary value problems are numerically solved. The secondary equilibrium paths and the post-buckling configurations of the rod with linearly varied cross-sections are presented.

Based on the nonlinear geometric theory of axially extensible beams and by using the shooting method, the thermal post-buckling responses of an elastic beams, with immovably simply supported ends and subjected to a transversely non-uniformly distributed temperature rising, were investigated. Especially, the influences of the transverse temperature change on the thermal post-buckling deformations were examined and the corresponding characteristic curves were plotted. The numerical results show that the equilibrium paths of the beam are similar to what of an initially deformed beam because of the thermal bending moment produced in the beam by the transverse temperature change.

A fixed point theorem is used to study a singular second order three-point boundary value problem. The problem is more general. Combining the method of constructing Green functions with operators defined piecewise, the existence result of positive solutions to a singular second order three-point boundary value problem is established. The nonlinearity can be allowed to change sign.

With prior estimate method, the existence, uniqueness, stability and large time behavior of the solution of second initial-boundary value problem for a fast diffusion equation with nonlinear boundary conditions are investigated. The main results are: 1) there exists only one global weak solution which continuously depends on initial value; 2) when t

Based on the nonlinear geometric theory of extensible rods, an exact mathematical model of thermal post-buckling behavior of uniformly heated elastic rods with axially immovable ends is developed, in which the arc length s(x) of axial line and the longitudinal displacement u(x) are taken as the basic unknown functions. This is a two point boundary value problem of first order ordinary differential equations with strong nonlinearity. By using shooting method and analytical continuation, the nonlinear boundary value problems are numerically solved. The thermal post-buckled states of the rods with transversely simply supported and clamped ends are obtained respectively and the corresponding numerical data tables and characteristic curves are also given.

A new τ-method is presented for the two dimensional linear boundary value problems. Theoretical and numerical analyses are presented. There results indicate that our method works nicely and efficiently.

In this paper, the axisymmetric buckled states of an annular sandwich plate (Reissner-type sandwich plate) with the clamped inner edge which is subjected to a uniform radial compressive thrust at the clamped outer edge are studied. Firstly, the basic equation of the buckled problem is derived. Secondly, the critical loads for some parameters are obtained by using the shooting method. Finally, we discuss the existence of the buckled states in the vicinity of the critical loads and obtain the asymptotic expansions of the buckled states.

In this paper, using the differentiability of the solution with respect to the initial value and the parameter, we present a method which, different from Liapunov's direct method, will determine the stability of the non-stationary solution of the initial value problem when the non-stationary solution remains unknown.

On the basis of von Karman equations, the axisymmetric buckling and post-buckling of annular plates on anelastic foundation is systematically discussed byusing shooting methods.

In this paper, the problem of large axisymmetric deformation of hyperelastic membrane is reformulated on the basis of the general theory of finite elasticity. From the fundamental equations derived here, the tensile instability of toroidal rubberlike membrane with Mooney-Rivlin's model constitutive relation is solved by using the shooting method. The upper and lower limit loads and the response curve of the displacement to the load are given.

In this paper, we consider the initial value problem of a class of Hill's equation having a small parameter. Using the solvable condition of boundary value problem and the stretched parameter method in the perturbation techniques, we present the method which can be applied to obtain asymptotic periodic solution of the initial value problem. As an example, we consider Mathieu equation and present its computational result.

The convergence of approximate solutions of boundary value problem of fourth order elliptic differential equations with uncontinuous coefficients in a rectangular region is investigated in this paper. This is useful for certain bending problems of rectangular plate on elastic supports.

By using the adaptive steplength integration scheme with a shooting technique, a rather difficult singular perturbation problem of ordinary differential equations with boundary layers can be calculated effectively. Computing examples are given in this paper which show the convergence within one iteration of the method in the case of a linear problem, the efficiency of the method for many boundary layers and turning points, especially the convenience in calculating multiple solutions. A comparison with traditional difference method is given at the end of this paper.

In this paper. first we introduce the shooting method (including the method of adjoints and the method of quasilinearization) for two-point boundary value problems of ordinary differential equations with parameters, and give out an example of solving hydrodynamic lubrication equation by using the shooting method. The calculated results are satisfactory.