Analytical solution to stagnation-point flow and heat transfer over a stretching sheet based on homotopy analysis

This paper is concerned with two-dimensional stagnation-point steady flow of an incompressible viscous fluid towards a stretching sheet whose velocity is proportional to the distance from the slit. The governing system of partial differential equations is first transformed into a system of dimensionless ordinary differential equations. Analytical solutions of the velocity distribution and dimensionless temperature profiles are obtained for different ratios of free stream velocity and stretching velocity, Prandtl number, Eckert number and dimensionality index in series forms using homotopy analysis method（HAM）. It is shown that a boundary layer is formed when the free stream velocity exceeds the stretching velocity, and an inverted boundary layer is formed when the free stream velocity is less than the stretching velocity. Graphs are presented to show the effects of different parameters.     

ANAI.YTIC SOLUTIONS FOR A FINITE CRACK RUPTURING WITH TYPE-II AT ONE TIP WITH EXCITING AND DECAYING RPOCESSES

Under the action of Rayleigh damping,when the shear stressexerts at the boundary of the crack and causes one tip ofthe crack to rupture with varying velocity,by the singularperturbation method,we reduce the governing nonlinearpartial differential equations to a system of linear onesand solve them by using generalized Fourier series.     

Differential transformation method for studying flow and heat transfer due to stretching sheet embedded in porous medium with variable thickness,variable thermal conductivity,and thermal radiation

This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.     

MHD free convective flow past semi-infinite vertical permeable wall

In this paper,the basic equations governing the flow and heat transfer of an incompressible viscous and electrically conducting fluid past a semi-infinite vertical permeable plate in the form of partial differential equations are reduced to a set of non-linear ordinary differential equations by applying a suitable similarity transformation.Approximate solutions of the transformed equations are obtained by employing the perturbation method for two cases,i.e.,small and large values of the suction parameter.From the numerical evaluations of the solution,it can be seen that the velocity field at any point decreases as the values of the magnetic and suction parameters increase.The effect of the magnetic parameter is to increase the thermal boundary layer.It is also found that the velocity and temperature fields decrease with the increase in the sink parameter.     

Non-uniform slot suction/injection into mixed convection boundary layer flow over vertical cone

The aim of this work is to study the effect of non-uniform single and double slot suction/injection into a steady mixed convection boundary layer flow over a vertical cone, while the axis of the cone is inline with the flow. The governing boundary layer equations are transformed into a non-dimensional form by a group of non-similar trans- formations. The resulting coupled non-linear partial differential equations are solved nu- merically by employing the quasi-linearization technique and an implicit finite-difference scheme. Numerical computations are performed for different values of the dimensionless parameters to display the velocity and temperature profiles graphically. Also, numerical results are presented for the skin friction and heat transfer coefficients. Results indicate that the skin friction and heat transfer coefficients increase with non-uniform slot suction, but the effect of non-uniform slot injection is just opposite.     

INCREMENTAL HARMONIC BALANCE METHOD FOR AIRFOIL FLUTTER WITH MULTIPLE STRONG NONLINEARITIES

The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.     

Critical velocity of sandwich cylindrical shell under moving internal pressure

Critical velocity of an infinite long sandwich shell under moving internal pressure is studied using the sandwich shell theory and elastodynamics theory. Propagation of axisymmetric free harmonic waves in the sandwich shell is studied using the sandwich shell theory by considering compressibility and transverse shear deformation of the core, and transverse shear deformation of face sheets. Based on the elastodynamics theory, displacement components expanded by Legendre polynomials, and position-dependent elastic constants and densities are introduced into the equations of motion. Critical velocity is the minimum phase velocity on the desperation relation curve obtained by using the two methods. Numerical examples and the finite element （FE） simulations are presented. The results show that the two critical velocities agree well with each other, and two desperation relation curves agree well with each other when the wave number k is relatively small. However, two limit phase velocities approach to the shear wave velocities of the face sheet and the core respectively when k limits to infinite. The two methods are efficient in the investigation of wave propagation in a sandwich cylindrical shell when k is relatively small. The critical velocity predicted in the FE simulations agrees with theoretical prediction.     

AN APPROXIMATE SOLUTION CONSIDERING FLOW INERTIA BETWEEN CONICAL SURFACES

In this paper, analytical calculation expressions of the pressure distribution, velocity distribution and the rate of the flow between conical surfaces are found by using the method of iterative approximate solution when the inertia terms of the Navier-Stokes equations in conical coordinates are taken into account. Furthermore, we compare the centrifugal flow with the centripetal flow of axisymmetrical passing flow.     

Flow of variable thermal conductivity fluid due to inclined stretching cylinder with viscous dissipation and thermal radiation

The aim of the present study is to investigate the flow of the Casson fluid by an inclined stretching cylinder. A heat transfer analysis is carried out in the presence of thermal radiation and viscous dissipation effects. The temperature dependent thermal conductivity of the Casson fluid is considered. The relevant equations are first simplified under usual boundary layer assumptions, and then transformed into ordinary differential equations by suitable transformations. The transformed ordinary differential equations are computed for the series solutions of velocity and temperature. A convergence analysis is shown explicitly. Velocity and temperature fields are discussed for different physical parameters by graphs and numerical values. It is found that the velocity decreases with the increase in the angle of inclination while increases with the increase in the mixed convection parameter. The enhancement in the thermal conductivity and radiation effects corresponds to a higher fluid temperature. It is also found that heat transfer is more pronounced in a cylinder when it is compared with a flat plate. The thermal boundary layer thickness increases with the increase in the Eckert number. The radiation and variable thermal conductivity decreases the heat transfer rate at the surface.     

Effect of double dispersion on non-Darcy mixed convective flow over vertical surface embedded in porous medium

A numerical study of a non-Darcy mixed convective heat and mass transfer flow over a vertical surface embedded in a dispersion, melting, and thermal radiation is porous medium under the effects of double investigated. The set of governing boundary layer equations and the boundary conditions is transformed into a set of coupled nonlinear ordinary differential equations with the relevant boundary conditions. The transformed equations are solved numerically by using the Chebyshev pseudospectral method. Comparisons of the present results with the existing results in the literature are made, and good agreement is found. Numerical results for the velocity, temperature, concentration profiles, and local Nusselt and Sherwood numbers are discussed for various values of physical parameters.     

Reduced projection augmented Lagrange bi-conjugate gradient method for contact and impact problems

Based on the numerical governing formulation and non-linear complementary conditions of contact and impact problems,a reduced projection augmented Lagrange bi- conjugate gradient method is proposed for contact and impact problems by translating non-linear complementary conditions into equivalent formulation of non-linear program- ming.For contact-impact problems,a larger time-step can be adopted arriving at numer- ical convergence compared with penalty method.By establishment of the impact-contact formulations which are equivalent with original non-linear complementary conditions,a reduced projection augmented Lagrange bi-conjugate gradient method is deduced to im- prove precision and efficiency of numerical solutions.A numerical example shows that the algorithm we suggested is valid and exact.     

MHD axisymmetric flow of third grade fluid between porous disks with heat transfer

The magnetohydrodynamic(MHD) flow of the third grade fluid between two permeable disks with heat transfer is investigated.The governing partial differential equations are converted into the ordinary differential equations by suitable transformations.The transformed equations are solved by the homotopy analysis method(HAM).The expressions for square residual errors are defined,and the optimal values of convergencecontrol parameters are selected.The dimensionless velocity and temperature fields are examined for various dimensionless parameters.The skin friction coefficient and the Nusselt number are tabulated to analyze the effects of dimensionless parameters.     

Nonstationary plane contact problem in theory of elasticity for conformal cylindrical surfaces

A numerical–analytical approach is described to investigate the process of impact interaction between a long smooth rigid body and the surface of a circular cylindrical cavity in elastic space. A non-stationary mixed initial boundary value problem is formulated with a priori unknown boundaries moving with variable velocity. The problem is solved using the methods of the theory of integral transforms, expansion of desired variables into a Fourier series, and the quadrature method to reduce the problem to solving a system of linear algebraic equations at each time step. Some concrete numerical computations are presented.The cylindrical body mass and radius impact on the proile of the transient process of contact interaction has been analysed.     

Iterative and adjusting method for computing stream function and velocity potential in limited domains and convergence analysis

The stream function and the velocity potential can be easily computed by solving the Poisson equations in a unique way for the global domain.Because of the various assumptions for handling the boundary conditions,the solution is not unique when a limited domain is concerned.Therefore,it is very important to reduce or eliminate the effects caused by the uncertain boundary condition.In this paper,an iterative and adjusting method based on the Endlich iteration method is presented to compute the stream function and the velocity potential in limited domains.This method does not need an explicitly specifying boundary condition when used to obtain the effective solution,and it is proved to be successful in decomposing and reconstructing the horizontal wind field with very small errors.The convergence of the method depends on the relative value for the distances of grids in two different directions and the value of the adjusting factor.It is shown that applying the method in Arakawa grids and irregular domains can obtain the accurate vorticity and divergence and accurately decompose and reconstruct the original wind field.Hence,the iterative and adjusting method is accurate and reliable.     

AN APPROXIMATE SOLUTION CONSIDERING FLOW INERTIA BETWEEN SPHERICAL SURFACES

In this paper, the analytical expressions of the pressure distribution, velocity distribution and discharge of the flow between spherical surfaces are found by using the method of iterative approximate solution when the inertia terms of Navier-Stokes equations in spherical coordinates are taken into consideration. Furthermore, using these expressions, we can directly obtain the corresponding analytical expressions of the laminar radial flow between parallel disks, which are fully identical with corresponding results presented by refs.     

Local non-similarity solution to impact of chemical reaction on MHD mixed convection heat and mass transfer flow over porous wedge in the presence of suction / injection

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.     

Nonlinear coupled dynamics of liquid-filled spherical container in microgravity

Nonlinear coupled dynamics of a liquid-filled spherical container in microgravity are investigated. The governing equations of the low-gravity liquid sloshing in a convex axisymmetrical container subjected to lateral excitation is obtained by the variational principle and solved with a modal analysis method. The variational formulas are transformed into a frequency equation in the form of a standard eigenvalue problem by the Galerkin method, in which admissible functions for the velocity potential and the liquid flee surface displacement are determined analytically in terms of the Gaussian hypergeometric series. The coupled dynamic equations of the liquid-filled container are derived using the Lagrange＇s method and are numerically solved. The time histories of the modal solutions are obtained in numerical simulations.     

LARGE RANGE ANALYSIS FOR NONLINEAR DYNAMIC SYSTEMS——ELEMENT MAPPING METHOD

This paper presents a new method for global analysis of nonlinear system. By means of transforming the nonlinear dynamic problems into point mapping forms which are single-valued and continuous, the state space can be regularly divided into a certain number of finitely small triangle elements on which the non-linear mapping can be approximately substituted by the linear mapping given by definition. Hence, the large range distributed problem of the mapping fixed points will be simplified as a process for solving a set of linear equations. Still further, the exact position of the fixed points can be found by the iterative technique. It is convenient to judge the stability of fixed points and the shrinkage zone in the state space by using the deformation matrix of linear mapping. In this paper, the attractive kernel for the stationary fixed points is defined, which makes great advantage for describing the attractive domains of the fixed points. The new method is more convenient and effective than the cell     

Thermal criticality for a reactive gravity driven thin film flow of a third-grade fluid with adiabatic free surface down an inclined plane

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-Padé 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.     

Modeling the interaction of solitary waves and semi-circular breakwaters by using unsteady reynolds equations

A vertical 2 -D numerical wave model was developed based on unsteady Reynolds equations. In this model, the k-epsilon models were used to close the Reynolds equations, and volume of fluid( VOF) method was used to reconstruct the free surface. The model was verified by experimental data. Then the model was used to simulate solitary wave interaction with submerged, alternative submerged and emerged semi-circular breakwaters. The process of velocity field, pressure field and the wave surface near the breakwaters was obtained. It is found that when the semi-circular breakwater is submerged, a large vortex will be generated at the bottom of the lee side wall of the breakwater ; when the still water depth is equal to the radius of the semi-circular breakwater, a pair of large vortices will be generated near the shoreward wall of the semi-circular breakwater due to wave impacting, but the velocity near the bottom of the lee side wall of the breakwater is always relatively small. When the semi-circular breakwater is emerged, and solitary wave cannot overtop it, the solitary wave surface will run up and down secondarily during reflecting from the breakwater. It can be further used to estate the diffusing and transportation of the contamination and transportation of suspended sediment.