Calculation of the boundary layer on the electrically conducting wall of a plane channel

There are presently available quite a large number of works devoted to the study of the motion of an electrically conducting fluid in boundary layers formed on electrodes or on the nonconducting walls of various MHD devices. However, the methods of solving the boundary layer equations in these studies are based on various simplifying assumptions which allow the problem to be reduced to the solution of a system of ordinary differential equations. Thus, in [1] there is imposed on the flow the special magnetic fieldH1/x, which enables the problem to be reduced to the self-similar form, while in the studies of other authors [2, 3] either the solution is sought in the form of expansions in x, or it is assumed that the problem is locally self-similar [4]. In the present paper we construct the solution of the MHD boundary layer equations which is obtained by one of the numerical methods which has long been used for solving the boundary layer equations for a nonconducting fluid.     

The solution of large deflection problem of thin circular plate by the method of composite expansion

In this paper, the method of composite expansion in perturbation theory is used for the solution of large deflection problem of thin circular plate. In this method. the outer field solution and the inner boundary layer solution are combined together to satisfy all the boundary conditions. In this paper, Hencky’s membrane solution is used for the first approximation in outer field solution, and then the second approximate solution is obtained. The inner boundary layer solution is found on the bases of boundary layer coondinate. In this paper, the reciprocal ratio of maximum deflection and thickness of the plate is used as the small parameter. The results of this paper improves quite a bit in comparison with the results obtained in 1948 by Chien Wei-zang.     

Exact solution of conductive heat transfer in cylindrical composite laminate

This paper presents an exact solution for steady-state conduction heat transfer in cylindrical composite laminates. This laminate is cylindrical shape and in each lamina, fibers have been wound around the cylinder. In this article heat transfer in composite laminates is being investigated, by using separation of variables method and an analytical relation for temperature distribution in these laminates has been obtained under specific boundary conditions. Also Fourier coefficients in each layer obtain by solving set of equations that related to thermal boundary layer conditions at inside and outside of the cylinder also thermal continuity and heat flux continuity between each layer is considered. In this research LU factorization method has been used to solve the set of equations.     

Approximate solution of the problem of the three-dimensional boundary layer in an incompressible fluid

A method of successive approximations is proposed for the solution of the equations of the three-dimensional incompressible boundary layer on bodies of arbitrary shape. A coordinate system connected with the streamlines of the external nonviscous flow is used. It is assumed that the velocity across the external streamlines is small. When the intensity of secondary flow is low the equations describing the boundary layer in an incompressible fluid are reduced to a form analogous to the equations for the boundary layer on axially symmetrical bodies. An approximate analytical solution is obtained for the velocity and for the friction in the form of equations which can be used for any problems of a three-dimensional incompressible boundary layer. The method developed was applied to the problem of the three-dimensional boundary layer at a plate with a cylindrical obstacle in the presence of a slip angle.     

The laminar boundary layer near a body corner point

A method is developed for calculating the characteristics of a laminar boundary layer near a body contour corner point, in the vicinity of which the outer supersonic stream passes through a rarefaction flow. In the study we use the asymptotic solution of the Navier-Stokes equations in the region with large longitudinal gradients of the flow functions for large values of the Reynolds number, the general form of which was used in [1].The pressure, heat flux, and friction distributions along the body surface are obtained. For small pressure differentials near the corner the solution of the corresponding equations for small disturbances is obtained in analytic form.The conventional method for studying viscous gas flow near body surfaces for large values of the Reynolds number is the use of the Prandtl boundary layer theory. Far from the body the asymptotic solution of the Navier-Stokes equations in the first approximation reduces to the solution of the Euler equations, while near the body it reduces to the solution of the Prandtl boundary layer equations. The characteristic feature of the boundary layer region is the small variation of the flow functions in the longitudinal direction in comparison with their variation in the transverse direction. However, in many cases this condition is violated.The necessity arises for constructing additional asymptotic expansions for the region in which the longitudinal and transverse variations of the flow functions are quantities of the same order. The general method for constructing asymptotic solutions for such flows with the use of the known method of outer and inner expansions is presented in [1].In the following we consider the flow in a laminar boundary layer for the case of a viscous supersonic gas stream in the vicinity of a body corner point. Behind the corner the flow separates from the body surface and flows around a stagnant zone, in which the pressure differs by a specified amount from the pressure in the undisturbed flow ahead of the point of separation. A pressure (rarefaction) disturbance propagates in the subsonic portion of the boundary layer upstream for a distance which in order of magnitude is equal to several boundary layer thicknesses. In the disturbed region of the boundary layer the longitudinal and transverse pressure and velocity disturbances are quantities of the same order. In this study we construct additional asymptotic expansions in the first approximation and calculate the distributions of the pressure, friction stress, and thermal flux along the body surface.     

MHD rotating flow of a viscous fluid over a shrinking surface

This study is concerned with the magnetohydrodynamic (MHD) rotating boundary layer flow of a viscous fluid caused by the shrinking surface. Homotopy analysis method (HAM) is employed for the analytic solution. The similarity transformations have been used for reducing the partial differential equations into a system of two coupled ordinary differential equations. The series solution of the obtained system is developed and convergence of the results are explicitly given. The effects of the parameters M, s and λ on the velocity fields are presented graphically and discussed. It is worth mentioning here that for the shrinking surface the stable and convergent solutions are possible only for MHD flows.     

Free convection from a vertical permeable circular cone with non-uniform surface heat flux

 In this paper, the problem of laminar free convection from a vertical permeable circular cone maintained with non-uniform surface heat flux is considered. The governing boundary layer equations are reduced non-similar boundary layer equations with surface heat flux proportional to x ⁿ (where x is the distance measured from the leading edge). The solutions of the reduced equations are obtained by using three distinct solution methodologies; namely, (i) perturbation solution for small transpiration parameter, ξ, (ii) asymptotic solution for large ξ, and (iii) the finite difference solutions for all ξ. The solutions are presented in terms of local skin-friction and local Nusselt number for smaller values of Prandtl number and heat flux gradient and are displayed in tabular form as well as graphically. Effects of pertinent parameters on velocity and temperature profiles are also shown graphically. Solutions obtained by finite difference method are also compared with the perturbation solutions for small and large ξ and found to be in excellent agreement. Received on 1 October 1999     

Unsteady three‐dimensional boundary layer flow due to a stretching surface in a micropolar fluid

In this paper, the unsteady three‐dimensional boundary layer flow due to a stretching surface in a viscous and incompressible micropolar fluid is considered. The partial differential equations governing the unsteady laminar boundary layer flow are solved numerically using an implicit finite‐difference scheme. The numerical solutions are obtained which are uniformly valid for all dimensionless time from initial unsteady‐state flow to final steady‐state flow in the whole spatial region. The equations for the initial unsteady‐state flow are also solved analytically. It is found that there is a smooth transition from the small‐time solution to the large‐time solution. The features of the flow for different values of the governing parameters are analyzed and discussed. The solutions of interest for the skin friction coefficient with various values of the stretching parameter c and material parameter K are presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Boundary layer flow of Reiner-Philippoff fluids

An analysis is made of the boundary layer flow of Reiner-Philippoff fluids. This work is an extension of a previous analysis by Hansen and Na [A.G. Hansen and T.Y. Na, Similarity solutions of laminar, incompressible boundary layer equations of non-Newtonian fluids. ASME 67-WA/FE-2, presented at the ASME Winter Annual Meeting, November (1967)], where the existence of similar solutions of the boundary layer equations of a class of general non-Newtonian fluids were investigated. It was found that similarity solutions exist only for the case of flow over a 90° wedge and, being similar, the solution of the non-linear boundary layer equations can be reduced to the solution of non-linear ordinary differential equations. In this paper, the more general case of the boundary layer flow of Reiner-Philippoff fluids over other body shapes will be considered. A general formulation is given which makes it possible to solve the boundary layer equations for any body shape by a finite-difference technique. As an example, the classical solution of the boundary layer flow over a flat plate, known as the Blasius solution, will be considered. Numerical results are generated for a series of values of the parameters in the Reiner-Philippoff model.     

On automodel numerical and asymptotic solutions of boundary-layer equations for high rates of injection

Automodel solutions of the equations of a laminar, multicomponent, isothermal boundary layer are considered for high rates of injection. The asymptotic velocity profiles and the thickness of the boundary layer are given for various negative pressure gradients (>0), A numerical solution is presented for the boundary-layer equations when injection involves the flow of a gas mixture comprising hydrogen, nitrogen, and carbon dioxide around the surface. The asymptotic solution is compared with the numerical solution, and its ranges of applicability are established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 49–52, May–June, 1971.     

MHD boundary-layer flow of an upper-convected Maxwell fluid in a porous channel

Two-dimensional magnetohydrodynamic (MHD) boundary layer flow of an upper-convected Maxwell fluid is investigated in a channel. The walls of the channel are taken as porous. Using the similarity transformations and boundary layer approximations, the nonlinear partial differential equations are reduced to an ordinary differential equation. The developed nonlinear equation is solved analytically using the homotopy analysis method. An expression for the analytic solution is derived in the form of a series. The convergence of the obtained series is shown. The effects of the Reynolds number Re, Deborah number De and Hartman number M are shown through graphs and discussed for both the suction and injection cases.     

Conjugate free convection over a vertical slender hollow cylinder embedded in a porous medium

A numerical study of the steady conjugate free convection over a vertical slender, hollow circular cylinder with the inner surface at a constant temperature and embedded in a porous medium is reported. The governing boundary layer equations for the fluid-saturated porous medium over the cylinder along with the one-dimensional heat conduction equation for the cylinder are cast into dimensionless form, by using a non-similarity transformation. The resulting non-similarity equations with their corresponding boundary conditions are solved by using the Keller box method. Emphasis is placed on the effects caused by the wall conduction parameter, p, and calculations have covered a wide range of this parameter. Heat transfer results including the temperature profiles, the interface temperature profiles and the local Nusselt number are presented. Received on 17 November 1997     

Temperature distribution in a Newtonian fluid injected between two semi-infinite plates

The analysis of the temperature distribution in time and place of a hot heat-conducting Newtonian fluid injected between two cooled parallel plates is presented. The 2-dimensional flow has a free flow front moving with constant velocity. The kernel of the fluid remains almost at the inlet temperature, but at the walls boundary layers occur with steeply descending temperature. The inner solutions inside these boundary layers are determined. To this end, the total region is divided into three distinct regions: the region G_I far behind the flow front, the flow front region G_II, and the intermediate region G_III between G_I and G_II. The asymptotics owing to each region are presented. The fundamental small parameter here is the thickness-to-length ratio of the 2-dimensional flow region. In most of the cases, similarity solutions are found. In the flow front region, for the formulation of the inner solution a Wiener–Hopf technique is used. Via matching procedures, the separate boundary layers are linked to each other to form one global boundary layer for the whole front region. All calculations in this paper are performed by analytical means, and all results are in analytical form. Comparison of our results with numerical solutions shows perfect agreement.     

Use of the generalized similarity method for calculating the inner zone of a turbulent boundary layer

The use of the generalized similarity method for calculating laminar boundary layers has been fully justified (see [1, §113, 114, 148]). The replacement of the partial differential equations by ordinary differential equations, their universality and the possibility of physically interpreting the solutions in the first, parametric stage of the calculations, which distinguish the generalized similarity method from direct numerical integration methods, are preserved in the case of a turbulent boundary layer also. A comparison of the calculated and experimental velocity profiles in the inner zone of the turbulent boundary layer suggests that the generalized similarity method could be used for calculating the turbulent layer as a whole.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 25–34, September–October, 1990.     

Forced convection from a microstructure on a flat plate

In this work, a cooling of a flat microelectronic structure with single-phase forced convection is investigated. The axial conduction, usually neglected in boundary layer theory, is considered here since the length of the heated element is in the same order of magnitude as the thickness of the boundary layer. The microstructure represents a package of chips mounted flush with the surface of the plate, and uniformly heated with a constant heat flux. The differential method is used to reduce the governing partial differential equations to ordinary differential ones, which are solved numerically by the use of a computational code developed by the authors. This code is based on Keller–Box method. The temperature profiles and Nusselt numbers are plotted at several locations on the heated element and are given as functions of the Reynolds number at the beginning of heated microstructure and of the ratio of unheated to heated length. Furthermore, the average Nusselt numbers on the heated length are computed for Prandtl numbers in the range 0.7≤Pr≤7,000. The results are compared to the boundary layer solution of unheated starting length problem. The results will be used as a baseline for successively more complex situations of cooling in electronics.     

Uniqueness of the self-similar solutions of three-dimensional laminar boundary-layer equations

The equations of the three-dimensional laminar boundary layer on lines of flow outflow and inflow are studied for conical outer flow under the assumption that the Prandtl number and the productρμ are constant. It is shown that in the case of a positive velocity gradient of the secondary flow (α₁>0) the additional conditions which result from the physical flow pattern determine a unique solution of the system of boundary-layer equations. For a negative velocity gradient of the secondary flow (α₁≤0) these conditions are satisfied by two solutions. An approximate solution is obtained for the boundary layer equations which is in rather good agreement with the numerical integration results. Compressible gas flow in a three-dimensional laminar boundary layer is described by a system of nonlinear differential equations whose solution is not unique for given boundary conditions. Therefore additional conditions resulting from the physical pattern of the gas flow are imposed on the resulting solution. In the solution of problems with a negative pressure gradient these additional conditions are sufficient for a unique selection of the solution of the boundary-layer equations. However, in the case of a positive pressure gradient the solution of the boundary-layer equations satisfying the boundary and additional conditions may not be unique. In particular, in [1] in a study of a three-dimensional laminar boundary layer in the vicinity of the stagnation point it was shown that for $$c = {{\frac{{\partial v_e }}{{\partial y}}} \mathord{\left/ {\vphantom {{\frac{{\partial v_e }}{{\partial y}}} {\frac{{\partial u_e }}{{\partial x}}}}} \right. \kern-\nulldelimiterspace} {\frac{{\partial u_e }}{{\partial x}}}} > 0$$ the solution is unique, while for c<0 there are two solutions. In the present paper we study the question of the uniqueness of the self-similar solution of the three-dimensional laminar boundary-layer equations on lines of flow outflow and inflow for a conical outer flow.     

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].     

Axisymmetric thermal stress in a thin circular disk by the method of matched asymptotic expansions

Using the method of matched asymptotic expansions, a solution is constructed for thin disks with stress-free edges which consists of an analytical solution in the interior and a boundary layer correction. The analytical solution is shown to be accurate to order (h/a)² compared to unity, and is that given by the classical theory of thin plates. There is no boundary layer in the case of a clamped disk, so that the thin plate solution is valid throughout. The solution of the boundary layer equations was obtained numerically by solving the finite difference form of the governing equations using an iterative scheme called “Dynamic Relaxation”.     

水泵叶片表面三维边界层计算

本文用有限差分法计算混流式可逆水力机械水泵工况叶片表面的三维边界层。水泵叶轮中主流区的三维势流由直接边界元法计算。对于叶片面附近的粘性流动。用三维半正交坐标系中的边界层方程表示。为了提高计算精度采用贴体坐标技术生成边界层区域的计算网格。并利用Cebeci等变换函数及Keller差分格式离散方程。用分块解法求解。计算叶轮叶片表面的压力分布与相应试验结果进行了对比。     

Boundary layer development in a gas flow with stagnation enthalpy varying across the streamlines

We consider a laminar boundary layer for which the stagnation enthalpy specified in the initial section is variable with height. Such problems arise, for example, for bodies located in the wake behind another body, for hypersonic flow past slender blunted bodies (as a result of the large transverse entropy gradients in the highentropy layer), for stepwise variation of the temperature of a surface on which there is an already developed boundary layer, for sudden expansion of the boundary layer as a result of its flow past a corner of the surface, etc.Strictly, we should in such cases solve the boundary layer equations (if the longitudinal gradients are much smaller than the transverse) with the specified initial distribution of the quantities. However, from the physical point of view, the distributed region may be broken down into two regions, the near-wall boundary layer and an outer region which is a gas flow with constant velocity and the specified initial temperature profile, whose calculation yields the edge conditions for the boundary layer. The boundary between the regions is determined from the condition of adequately smooth matching of the solutions. This approach is much preferable to the first, since it permits avoiding (within the framework of boundary layer theory) the difficulties associated with the presence of a possible singularity at the initial point of the surface due to the discontinuity of the boundary conditions at this point, and also permits using conventional boundary layer theory if the effect of the viscosity in the outer region is not significant. However, this partition requires additional justifications of the possibility of independent determination of the solution in the outer region and the determination of the edge of the boundary layer, considered as the region of influence of the wetted surface. The boundary layer in a nonuniform flow has been considered in several works for a linear initial velocity or temperature profile [1–3].It should be noted that the linear initial enthalpy or velocity profiles for constant gas properties do not undergo changes under the influence of viscosity or thermal conductivity. Thus the fundamental characteristic features noted above which are associated with the presence of the two regions and their interaction in essence cannot be investigated using these examples.In this study we obtain and analyze the exact solutions of the equations of the compressible boundary layer for a power-law variation of the initial stagnation enthalpy profile as a function of the stream function for a constant initial velocity. Here it is shown that the influence of the boundary conditions at the wall are actually localized in the near-wall boundary layer, which is similar in dimensions to the conventional velocity or thermal boundary layers. In the region which is external with relation to this layer, in accordance with the physical picture described above, the solution coincides with the solution of the Cauchy problem for the heat conduction equation, which describes the development of the initial temperature profile in an infinite steady-state flow with constant velocity.It is shown that for the sufficiently smooth initial profiles which are of interest in practice the outer flow undergoes practically no changes until we reach the inner boundary layer, and it may be calculated using the perfect gas laws.