Approximate solutions to MHD Falkner-Skan flow over permeable wall

The magnetohydrodynamic (MHD) Falkner-Skan boundary layer flow over a permeable wall in the presence of a transverse magnetic field is examined. The approximate solutions and skin friction coefficients of the MHD boundary layer flow are obtained by using a method that couples the differential transform method (DTM) with the Padé approximation called DTM-Padé. The approximate solutions are expressed in the form of a power series that can be easily computed with an iterative procedure. The approximate solutions are tabulated, plotted for the values of different parameters and compared with the numerical ones obtained by employing the shooting technique. It is found that the approximate solution agrees very well with the numerical solution, showing the reliability and validity of the present work. Moreover, the effects of various physical parameters on the boundary layer flow are presented graphically and discussed.     

Approximate analytical solutions and approximate value of skin friction coefficient for boundary layer of power law fluids

This paper presents a theoretical analysis for laminar boundary layer flow in a power law non-Newtonian fluids.The Adomian analytical decomposition technique is presented and an approximate analytical solution is obtained.The approximate analytical solution can be expressed in terms of a rapid convergent power series with easily computable terms.Reliability and efficiency of the approximate solution are verified by comparing with numerical solutions in the literature.Moreover,the approximate solution can be successfully applied to provide values for the skin friction coefficient of the laminar boundary layer flow in power law non-Newtonian fluids.     

Laminar boundary layer between two planes perpendicular to each other

In this paper, we obtain a third-order approximate solution for the laminar boundary layer between two planes perpendicular to each other.In boundary layer equations, the viscous and the inertial terms have the same quantity step. In this paper, at first, supposing that the inertial terms are bigger than the viscous terms, we solve the boundary layer equations, and then we suppose that the viscous terms are bigger than the inertial terms. At last, we take the mean value as the valid solution of the boundary layer equations.The first- and the second-order approximate solutions obtained in this paper coincide with the results in ref. [1], while the third-order solution obtained in this paper is better than that in ref. [1].     

Investigation of three-dimensional boundary layers on blunt bodies with permeable surface

Numerical and approximate analytic methods are used to investigate three-dimensional laminar boundary layers on blunt bodies with permeable surface in a supersonic gas stream. In the first approximation of the integral method of successive approximation an analytic solution is obtained to the problem for an impermeable surface, small values of the blowing parameter, and arbitrary suction. For large parameters of the blowing (or suction), whose velocity vector in the general case is directed at a certain angle to the vector of the outer normal to the body, asymptotic expressions are derived for the components of the frictional stress and the heat flux. A numerical solution is obtained to the equations of the three-dimensional boundary layer in a wide range of variation of the blowing (or suction) parameter. The accuracy and region of applicability of the analytic solutions is estimated by comparison with the numerical solutions. On the basis of the solutions obtained in the present paper and the work of other authors an expression is proposed for calculating the heat fluxes to a perfectly catalytic surface of a body in a three-dimensional supersonic flow of dissociated or ionized air. The present paper continues earlier work of the authors [1, 2] on boundary layers in the neighborhood of a symmetry plane and on sweptback wings of infinite span.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 49–58, May–June, 1982.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

MHD effect of mixed convection boundary-layer flow of Powell-Eyring fluid past nonlinear stretching surface

Sufficient conditions are found for the existence of similar solutions of the mixed convection flow of a Powell-Eyring fluid over a nonlinear stretching permeable sur- face in the presence of magnetic field. To achieve this, one parameter linear group trans- formation is applied. The governing momentum and energy equations are transformed to nonlinear ordinary differential equations by use of a similarity transformation. These equations are solved by the homotopy analysis method （HAM） to obtain the approximate solutions. The effects of magnetic field, suction, and buoyancy on the Powell-Eyring fluid flow with heat transfer inside the boundary layer are analyzed. The effects of the non- Newtonian fluid （Powell-Eyring model） parameters ε and δon the skin friction and local heat transfer coefficients for the cases of aiding and opposite flows are investigated and discussed. It is observed that the momentum boundary layer thickness increases and the thermal boundary layer thickness decreases with the increase in ε whereas the momentum boundary layer thickness decreases and thermal boundary layer thickness increases with the increase in δ for both the aiding and opposing mixed convection flows.     

Fluctuating flow of a second-order fluid near a stagnation point

Summary The boundary layer equations for axisymmetric flow of an incompressible second-order fluid have been deduced. The flow of such a fluid near a stagnation point when the main stream outside the boundary layer fluctuates in magnitude but not in direction has been discussed. The velocity distribution is found for various values of the steady mean in two limiting cases of small and large values of the frequency of the oscillation of the main stream. The frequency for which two approximate solutions overlap has been calculated in each case.     

Falkner-Skan方程的近似解析解

研究了粘性流体绕流楔型物体的Falkner-Skan边界层方程求解问题.利用Adomian拆分方法,通过引入Crocco变量变换将无穷区间的边界值问题转为初值问题并利用Padé逼近技巧确定初值,给出了一种有效的解析分解方法.进一步,本文设计了一种数值解法,将本文得到的近似解析解及数值结果与早期研究者Hartree等人的结果进行了比较,证明了本文提出的解法的有效性和可靠性.     

MHD flow of power-law fluid on moving surface with power-law velocity and special injection/blowing

The problem of magnetohydrodynamic （MHD） flow on a moving surface with the power-law velocity and special injection/blowing is investigated. A scaling group transformation is used to reduce the governing equations to a system of ordinary differen- tial equations. The skin friction coefficients of the MHD boundary layer flow are derived, and the approximate solutions of the flow characteristics are obtained with the homotopy analysis method （HAM）. The approximate solutions are easily computed by use of a high order iterative procedure, and the effects of the power-law index, the magnetic parameter, and the special suction/blowing parameter on the dynamics are analyzed. The obtained results are compared with the numerical results published in the literature, verifying the reliability of the approximate solutions.     

Approximate dynamic boundary conditions for a thin piezoelectric layer

Approximate dynamic boundary conditions of different orders are derived for the case of a thin piezoelectric coating layer bonded to an elastic material. The approximate boundary conditions are derived using series expansions of the elastic displacements and the electric potential in the thickness coordinate of the layer. All the expansion functions are then eliminated with the aid of the equations of motion and boundary/interface conditions of the layer. This results in boundary conditions on the elastic material that may be truncated to different orders in the thickness of the layer to obtain approximate boundary conditions. The approximate boundary conditions may be used as a replacement for the piezoelectric layer and thus simplify the analysis significantly. Numerical examples show that the approximate boundary conditions give good results for low frequencies and/or thin piezoelectric layers.     

坐标变换法数值求解微通道行波电场电渗流

采用坐标变换法数值求解了耦合的Poisson-Nernst-Planck (PNP)方程和Navier-Stokes(NS)方程, 研究二维狭窄微通道行波电场电渗流数值解. 数值结果表明,坐标变换法能有效降低电渗流解数值解在双电层的高梯度, 有效改善数值解的收敛性和稳定性. 坐标变换的电渗流数值解和原始坐标下的数值解完全一致. 坐标变换后采用简单的网格也能得到和原始坐标下复杂网格相同的解. 给出了滑移边界的近似解与完整的PNP-NS数值解的比较. 在双电层厚度与微通道深度比值(λ/H)很小的情况下(相对深通道), 两者的解基本一致. 但在λ/H较大时(相对浅通道)滑移边界的解高于电渗流速度.      

On the polynomial approximation of boundary-layer flow profiles

The behaviour of the polynomial approximation to the boundary layer velocity profile is investigated. Various orders of polynomials and 4 different schemes of “reasonable” boundary conditions are examined for applicability as approximate solutions to the Blasius flow over a flat plate. A variational formulation, based upon the local potential is used to obtain the solution. It is found that the best and most consistent results are obtained when a symmetric distribution of auxiliary boundary conditions on the wall and outer edge of the boundary layer is used. The 6th order polynomial of this type, for example, already gives a wall friction factor within 0.5% of the exact solution.     

Numerical method for unsteady compressible boundary layers driven by compression or expansion wave

A numerical analysis is presented for the unsteady compressible laminar boundary layer driven by a compression or expansion wave. Approximate or series expansion methods have been used for the problems because of the characteristics of the governing equations, such as non-linearity, coupling with the thermal boundary layer equation and initial conditions. Here a transformation of the governing equations and the numerical linearization technique are introduced to deal with the difficulties. First, the governing equations are transformed for the initial conditions by Howarth and semisimilarity variables. These transformations reduce the number of independent variables from three to two and the governing equations from partial to ordinary differential equations at the initial point. Next, the numerical linearization technique is introduced for the non-linearity and the coupling with the thermal boundary layer equation. Because the non-linear terms are linearized without sacrifice of numerical accuracy, the solutions can be obtained without numerical iterations. Therefore the exact numerical solution, not approximate or series expansion, can be obtained. Compared with the approximate or series expansion method, this method is much improved. Results are compared with the series expansion solutions.     

The coefficient of regeneration of the enthalpy in a three-dimensional laminar boundary layer

Self-similar solutions of the equations of a three-dimensional laminar boundary layer are of interest from two points of view. In the first place, they can be used to construct approximate calculating methods, making it possible to analyze several variants and to consider complex flows, in which it is impossible to neglect the interaction between the boundary layer and the external flow (for example, in the region of hypersonic interaction [1–3]). In the second place, the analysis of self-similar solutions permits clarifying the effect of individual parameters on one or another characteristic of the boundary layer and representing this effect in predictable form. One of the principal characteristics of a three-dimensional boundary layer, as also of a two-dimensional, is the coefficient of regeneration of the enthalpy. The value of this coefficient is needed for determining the temperature of a thermally insulated surface, as well as for finiing the real temperature (or enthalpy) head, which determines the value of the heat flux from a heated gas to the surface of the body around which the flow takes place. The article presents the results of calculations of the coefficient of regeneration of the enthalpy for locally self-similar solutions of the equations of a three-dimensional boundary layer, forming with flow around a cylindrical thermally insulated surface at an angle. It is clarified that the dependence of the coefficient of regeneration of the enthalpy on the determining parameters is not always continuous.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 60–63, January–February, 1973.     

Asymptotic investigation of turbulent buoyant jets

Results of a theoretical investigation are presented for the averaged characteristics of turbulent slightly swirling jets discharged vertically upward into a stagnant homogeneous medium. A semi-empirical model of turbulence is used to close the system of initial turbulent boundary layer equations. Exact and approximate analytical solutions are suggested. The results of calculations are compared with experimental and numerical data of other authors.     

A kinetic theory solution method for the Navier–Stokes equations

The kinetic-theory-based solution methods for the Euler equations proposed by Pullin and Reitz are here extended to provide new finite volume numerical methods for the solution of the unsteady Navier–Stokes equations. Two approaches have been taken. In the first, the equilibrium interface method (EIM), the forward- and backward-flowing molecular fluxes between two cells are assumed to come into kinetic equilibrium at the interface between the cells. Once the resulting equilibrium states at all cell interfaces are known, the evaluation of the Navier–Stokes fluxes is straightforward. In the second method, standard kinetic theory is used to evaluate the artificial dissipation terms which appear in Pullin's Euler solver. These terms are subtracted from the fluxes and the Navier–Stokes dissipative fluxes are added in. The new methods have been tested in a 1D steady flow to yield a solution for the interior structure of a shock wave and in a 2D unsteady boundary layer flow. The 1D solutions are shown to be remarkably accurate for cell sizes large compared to the length scale of the gradients in the flow and to converge to the exact solutions as the cell size is decreased. The steady-state solutions obtained with EIM agree with those of other methods, yet require a considerably reduced computational effort.     

A comparative assessment of fluid flow and heat transfer parameters in laminar flow around a circular cylinder and a sphere

Theoretical studies have been carried out for a comparative assessment of hydrodynamic boundary layer thickness, displacement thickness and shear stress at the wall for laminar flow around a circular cylinder and a sphere with the help of the approximate method due to Karman and Pohlhausen for two dimensional flow and the method as applied to bodies of revolution based on the work of F. W. Scholkemeier, respectively. Thermal boundary layer thickness and Nusselt number have been evaluated around the surface of the solids. Comparison is made with available solutions. The graphical presentation of the results depicts a concise and relative assessment of fluid flow and heat-transfer parameters for flow around cylinder and sphere.     

On the numerical solution of some two-dimensional boundary-contact delocalization problems

The elastic equilibrium of a multi-layer confocal elliptic ring is studied. The ring consists of steel, rubber and celluloid layers which differ in thickness and in the order in which they are placed relative to one another. Using the solutions of the considered problems, the following delocalization problem is solved: for a three-layer elliptic body, the external elliptic boundary of which is loaded by normal point force and the internal boundary is stress-free and the layers of which are in rigid or sliding contact with one another, by an appropriate choice of layer thickness and arrangement of the layers relative to one another we can obtain a sufficiently uniform distribution of normal displacements on the internal elliptic boundary. Numerical solutions are obtained by the boundary element method and the related graphs are constructed. For the two-layer ellipse, exact and approximate solutions of the same problem are obtained respectively by the method of separation of variables and by the boundary element method. The results obtained by both methods are compared and the conclusion as to the reliability of the numerical boundary element method is made.