Lattice Boltzmann simulation for high-speed compressible viscous flows with a boundary layer

In this study, the lattice Boltzmann method is employed for simulating high-speed compressible viscous flows with a boundary layer. The coupled double-distribution-function lattice Boltzmann method proposed by Li et al. (2007) is employed because of its good numerical stability and non-free-parameter feature. The non-uniform mesh construction near the wall boundary in fine grids is combined with an appropriate wall boundary treatment for the finite difference method in order to obtain accurate spatial resolution in the boundary layer problem. Three typical problems in high-speed viscous flows are solved in the lattice Boltzmann simulation, i.e., the compressible boundary layer problem, shock wave problem, and shock boundary layer interaction problem. In addition, in-depth comparisons are made with the non-oscillatory and non-free-parameter dissipation (NND) scheme and second order upwind scheme in the present lattice Boltzmann model. Our simulation results indicate the great potential of the lattice Boltzmann method for simulating high-speed compressible viscous flows with a boundary layer. Further research is needed (e.g., better numerical models and appropriate finite difference schemes) because the lattice Boltzmann method is still immature for high-speed compressible viscous flow applications.     

Well-posedness of thermal boundary layer equation in two-dimensional incompressible heat conducting flow with analytic datum

In this paper, we study the well-posedness of the thermal boundary layer equation in two-dimensional incompressible heat conducting flow. The thermal boundary layer equation describes the behavior of thermal layer and viscous layer for the two-dimensional incompressible viscous flow with heat conduction in the small viscosity and heat conductivity limit. When the initial datum are analytic, with respect to the tangential variable of the boundary, and without the monotonicity condition of the tangential velocity, by using the Littlewood-Paley theory, we obtain the local-in-time existence and uniqueness of solution to this thermal boundary layer problem.     

The applicability of continuum models in the transitional regime of hypersonic flow over blunt bodies

Hypersonic rarefied gas flow over blunt bodies in the transitional flow regime (from continuum to free-molecule) is investigated. Asymptotically correct boundary conditions on the body surface are derived for the full and thin viscous shock layer models. The effect of taking into account the slip velocity and the temperature jump in the boundary condition along the surface on the extension of the limits of applicability of continuum models to high free-stream Knudsen numbers is investigated. Analytic relations are obtained, by an asymptotic method, for the heat transfer coefficient, the skin friction coefficient and the pressure as functions of the free-stream parameters and the geometry of the body in the flow field at low Reynolds number; the values of these coefficients approach their values in free-molecule flow (for unit accommodation coefficient) as the Reynolds number approaches zero. Numerical solutions of the thin viscous shock layer and full viscous shock layer equations, both with the no-slip boundary conditions and with boundary conditions taking into account the effects slip on the surface are obtained by the implicit finite-difference marching method of high accuracy of approximation. The asymptotic and numerical solutions are compared with the results of calculations by the Direct Simulation Monte Carlo method for flow over bodies of different shape and for the free-stream conditions corresponding to altitudes of 75–150 km of the trajectory of the Space Shuttle, and also with the known solutions for the free-molecule flow regine. The areas of applicability of the thin and full viscous shock layer models for calculating the pressure, skin friction and heat transfer on blunt bodies, in the hypersonic gas flow are estimated for various free-stream Knudsen numbers.     

Interactive computing methods for aeroelastic analysis of turbomachinery

An interaction viscous‐inviscid method for efficiently computing steady and unsteady viscous flows is presented. The inviscid domain is modeled using a finite element discretization of the full potential equation. The viscous region is modeled using a finite difference boundary layer technique. The two regions are simultaneously coupled using the transpiration approach. A time linearization technique is applied to this interactive method. For unsteady flows, the fluid is assumed to be composed of a mean or steady flow plus a harmonically varying small unsteady disturbance. Numerically exact nonreflecting boundary conditions are used for the far field conditions. Results for some steady and unsteady, laminar and turbulent flow problems are compared to linearized Navier‐Stokes or time‐marching boundary layer methods. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

Hele Shaw Convection with Imposed Shear Flows: Boundary-Layer Formulation

The nonlinear convection forced by the boundaries of a Hele Shaw cell to align perpendicular to an imposed shear flow was analytically investigated by the boundary-layer method. The imposed shear flow may be a Couette flow that extends throughout the convecting layer or flow confined to a boundary, depending on the geometry of the Hele Shaw cell. This study examined the case in which the imposed shear flow has a boundary-layer structure and its interaction with the convecting interior. Analytical solutions for both the boundary layer and interior were obtained. The study revealed the following.For large aspect ratio A , the interaction of the imposed shear flow and convection is confined to the boundary layer. The boundary layer is a viscous rather than a thermal layer. The results showed that the range of validity of the Hele Shaw equations used in the literature is of order 1/ A ². For an asymptotically large aspect ratio A up to order 1/ A ², the velocity in the y -direction must be zero. The velocity in the x -direction and the z -direction has a parabolic dependence on y , but the temperature perturbation does not depend on y . These results may have implication for convection in porous media.     

HAM solutions for boundary layer flow in the region of the stagnation point towards a stretching sheet

In the present study, we have described the stagnation point flow of a viscous fluid towards a stretching sheet. The complete analytical solution of the boundary layer equation has been obtained by homotopy analysis method (HAM). The solutions are compared with the available numerical results obtained by Nazar et al. [Nazar R, Amin N, Filip D, Pop I. Unsteady boundary layer flow in the region of the stagnation point on a stretching sheet. Int J Eng Sci 2004;42:1241–53] and a good agreement is found. The convergence region is also computed which shows the validity of the HAM solution.     

A numerical study of boundary layer flow of viscous incompressible fluid past an inclined stretching sheet and heat transfer using nonpolynomial spline method

In the present paper, we study the boundary layer flow of viscous incompressible fluid over an inclined stretching sheet with body force and heat transfer. Considering the stream function, we convert the boundary layer equation into nonlinear third-order ordinary differential equation together with appropriate boundary conditions in an infinite domain. The nonlinear boundary value problem has been linearized by using the quasilinearization technique. Then, we develop a nonpolynomial spline method, which is used to solve the flow problem. The convergence analysis of the method is also discussed. We study the velocity function for different angles of inclination and Froude number with the help of various graphs and tables. Then using these in heat convection flow, we obtain the expression for temperature field. Skin friction is also calculated. The various results have been given in tables. At last, we calculated the Nusselt number.     

Nonlinear Stability of Nonstationary Cross-Flow Vortices in Compressible Boundary Layers

The nonlinear evolution of long-wavelength non stationary cross-flow vortices in a compressible boundary layer is investigated; the work extends that of Gajjar [1] to flows involving multiple critical layers. The basic flow profile considered in this paper is that appropriate for a fully three-dimensional boundary layer with O(1) Mach number and with wall heating or cooling. The governing equations for the evolution of the cross-flow vortex are obtained, and some special cases are discussed. One special case includes linear theory, where exact analytic expressions for the growth rate of the vortices are obtained. Another special case is a generalization of the Bassom and Gajjar [2] results for neutral waves to compressible flows. The viscous correction to the growth rate is derived, and it is shown how the unsteady nonlinear critical layer structure merges with that for a Haberman type of viscous critical layer.     

Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Effect of injection on the flow of Oldroyd fluid in the inlet region of a channel

The effect of injection on the flow of Oldroyd fluid in the inlet region of a channel has been investigated using the moment and energy integrals, taking into account the loss of energy due to viscous dissipation in the boundary layer. Analytical expression for boundary layer development has been presented.     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

双曲拉伸面上的流动及其热交换

研究不可压缩粘性流体,在双曲拉伸面上的边界层流动及其热传导．分别使用级数展开法和局部非相似（LNS）法,得到解析结果和数值结果,给出了表面摩擦和Nusselt数的解析结果和数值结果,并进行了互相比较．同时发现动量和热边界层厚度,随着离前缘距离的增加而减小．众所周知,线性拉伸项方程的解,可以作为双曲拉伸首次项方程的解．     

First stages of drop impact on a dry surface: asymptotic model

After impact of a viscous liquid drop on a dry wall surrounded by a gas, the drop surface is highly deformed, leading to the formation of an axisymmetrical lateral lamella along the wall. A local asymptotic model for the potential flow and unsteady boundary layer flow is developed to describe the lamella dynamics at early stages after impact. The second-order potential flow displaced by the unsteady boundary layer is taken into account. The lamella shape, its velocity and pressure are calculated with this model in parametrical forms. The three model parameters are evaluated here by fitting with recent experimental findings.     

On multiscale homogenization problems in boundary layer theory

This paper is concerned with the homogenization of the equations describing a magnetohydrodynamic boundary layer flow past a flat plate, the flow being subjected to velocities caused by injection and suction. The fluid is assumed incompressible, viscous and electrically conducting with a magnetic field applied transversally to the direction of the flow. The velocities of injection and suction and the applied magnetic field are represented by rapidly oscillating functions according to several scales. We derive the homogenized equations, prove convergence results and establish error estimates in a weighted Sobolev norm and in C ⁰-norm. We also examine the asymptotic behavior of the solutions of the equations governing a boundary layer flow past a rough plate with a locally periodic oscillating structure.     

Combined effects of non-uniform heat source/sink and thermal radiation on heat transfer over an unsteady stretching permeable surface

The present paper is concerned with the study of flow and heat transfer characteristics in the unsteady laminar boundary layer flow of an incompressible viscous fluid over continuously stretching permeable surface in the presence of a non-uniform heat source/sink and thermal radiation. The unsteadiness in the flow and temperature fields is because of the time-dependent stretching velocity and surface temperature. Similarity transformations are used to convert the governing time-dependent nonlinear boundary layer equations for momentum and thermal energy are reduced to a system of nonlinear ordinary differential equations containing Prandtl number, non-uniform heat source/sink parameter, thermal radiation and unsteadiness parameter with appropriate boundary conditions. These equations are solved numerically by applying shooting method using Runge–Kutta–Fehlberg method. Comparison of numerical results is made with the earlier published results under limiting cases. The effects of the unsteadiness parameter, thermal radiation, suction/injection parameter, non-uniform heat source/sink parameter on flow and heat transfer characteristics as well as on the local Nusselt number are shown graphically.     

两互相垂直的平面间的层流边界层

本文得到了两互相垂直的平面间的层流边界层的三级近似解.在边界层中,边界层方程中的粘性项和惯性项具有相同的数量级[3].本文则首先假定惯性项大于粘性项去求解边界层方程;然后,令粘性项大于贯性项.最后,取二者的平均值作为边界层方程的真实解.本文所得一级及二级近似解和文献[1]的结果相同.本文的三级近似解则较[1]的结果更精确.     

Classroom Notes

This note gives a simple‐minded approach to the two‐dimensional boundary layer equations. The pressure is eliminated from the equations of motion and the resulting equation is simplified by assuming that certain derivatives in the direction of the boundary are small compared with those at right angles to it. The simplified equation is then integrated to give a single boundary layer equation which, together with the stress rate of strain law and the continuity equation, is sufficient (in theory at least) to predict the flow.

The boundary layer equation as given does not depend on a particular form for the stress rate of strain law and could possibly form the basis for a non‐Newtonian investigation. The viscous boundary layer is given as a special case.     

Heat transfer in a Walter’s liquid B fluid over an impermeable stretching sheet with non-uniform heat source/sink and elastic deformation

In this present article an analysis is carried out to study the boundary layer flow behavior and heat transfer characteristics in Walter’s liquid B fluid flow. The stretching sheet is assumed to be impermeable, the effects of viscous dissipation, non-uniform heat source/sink in the presence and in the absence of elastic deformation (which was escaped from attention of researchers while formulating the viscoelastic boundary layer flow problems)on heat transfer are addressed. The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. Analytical solutions are obtained for the resulting boundary value problems. The effects of viscous dissipation, Prandtl number, Eckert number and non-uniform heat source/sink on heat transfer (in the presence and in the absence of elastic deformation) are shown in several plots and discussed. Analytical expressions for the wall frictional drag coefficient, non-dimensional wall temperature gradient and non-dimensional wall temperature are obtained and are tabulated for various values of the governing parameters. The present study reveals that, the presence of work done by deformation in the energy equation yields an augment in the fluid’s temperature.     

Effects of thermal radiation and magnetic field on unsteady mixed convection flow and heat transfer over an exponentially stretching surface with suction in the presence of internal heat generation/absorption

In this paper, the problem of unsteady laminar two-dimensional boundary layer flow and heat transfer of an incompressible viscous fluid in the presence of thermal radiation, internal heat generation or absorption, and magnetic field over an exponentially stretching surface subjected to suction with an exponential temperature distribution is discussed numerically. The governing boundary layer equations are reduced to a system of ordinary differential equations. New numerical method using Mathematica has been used to solve such system after obtaining the missed initial conditions. Comparison of obtained numerical results is made with previously published results in some special cases, and found to be in a good agreement.