BOUNDARY LAYER ASYMPTOTIC BEHAVIOR OF INCOMPRESSIBLE NAVIER-STOKES EQUATION IN A CYLINDER WITH SMALL VISCOSITY

The object of this article is to study the boundary layer appearing at large Reynolds number （small viscosity ε） incompressible Navier Stokes Equation in a cylinder in space dimension three. These are Navier-Stokes equations linearized around a fixed velocity flow： the authors study the convergence as ε →0 to the inviscid type equations, the authors define the correctors needed to resolve the boundary layer and obtain convergence results valid up to the boundary and the authors also study the behavior of the boundary layer when, simultaneously, time and the Reynolds number tend to infinity, in which case the boundary layer tends to pervade the whole domain.

Applicability of Temperature Discrete Equation to NMRF Boundary Layer Scheme in GRAPES Model

By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast(NMRF)boundary layer scheme in the GRAPES model,the adjusted discrete equation for temperature is obviously different from the original equation under the background of hydrostatic equilibrium and adiabatic hypothesis.In the present research,three discrete equations for temperature in the NMRF boundary layer scheme are applied,namely the original(hereafter NMRF),the adjustment(hereafter NMRF-gocp),and the one in the YSU boundary-layer scheme(hereafter NMRF-TZ).The results show that the deviations of height,temperature,U and V wind in the boundary layer in the NMRF-gocp and NMRF-TZ experiments are smaller than those in the NMRF experiment and the deviations in the NMRF-gocp experiment are the smallest.The deviations of humidity are complex for the different forecasting lead time in the three experiments.Moreover,there are obvious diurnal variations of deviations from these variables,where the diurnal variations of deviations from height and temperature are similar and those from U and V wind are also similar.However,the diurnal variation of humidity is relatively complicated.The root means square errors of 2m temperature(T2m)and 10m speed(V10m)from the three experiments show that the error of NMRF-gocp is the smallest and that of NMRF is the biggest.There is also a diurnal variation of T2m and V10m,where T2m has double peaks and V10m has only one peak.Comparison of the discrete equations between NMRF and NMRF-gocp experiments shows that the deviation of temperature is likely to be caused by the calculation of vertical eddy diffusive coefficients of heating,which also leads to the deviations of other elements.

Analytic solution for magnetohydrodynamic boundary layer flow of Casson fluid over a stretching/shrinking sheet with wall mass transfer

In this analysis,the magnetohydrodynamic boundary layer flow of Casson fluid over a permeable stretching/shrinking sheet in the presence of wall mass transfer is studied.Using similarity transformations,the governing equations are converted to an ordinary differential equation and then solved analytically.The introduction of a magnetic field changes the behavior of the entire flow dynamics in the shrinking sheet case and also has a major impact in the stretching sheet case.The similarity solution is always unique in the stretching case,and in the shrinking case the solution shows dual nature for certain values of the parameters.For stronger magnetic field,the similarity solution for the shrinking sheet case becomes unique.

A STUDY ON THE BOUNDARY LAYER FOR THE PLANAR MAGNETOHYDRODYNAMICS SYSTEM

The paper aims to estimate the thickness of the boundary layer for the planar MHD system with vanishing shear viscosity μ. Under some conditions on the initial and boundary data, we show that the thickness is of the order √μ｜lnμ｜. Note that this estimate holds also for the Navier-Stokes system so that it extends the previous works even without the magnetic effect.

Improvement for expansion of parabolized stability equation method in boundary layer stability analysis

An improved expansion of the parabolized stability equation(iEPSE) method is proposed for the accurate linear instability prediction in boundary layers. It is a local eigenvalue problem, and the streamwise wavenumber α and its streamwise gradient dα/dx are unknown variables. This eigenvalue problem is solved for the eigenvalue dα/dx with an initial α, and the correction of α is performed with the conservation relation used in the PSE. The i EPSE is validated in several compressible and incompressible boundary layers. The computational results show that the prediction accuracy of the i EPSE is significantly higher than that of the ESPE, and it is in excellent agreement with the PSE which is regarded as the baseline for comparison. In addition, the unphysical multiple eigenmode problem in the EPSE is solved by using the i EPSE. As a local non-parallel stability analysis tool, the i EPSE has great potential application in the eNtransition prediction in general three-dimensional boundary layers.

Self-consistent parabolized stability equation(PSE) method for compressible boundary layer

The parabolized stability equation （PSE） method has been proven to be a useful and convenient tool for the investigation of the stability and transition problems of boundary layers. However, in its original formulation, for nonlinear problems, the complex wave number of each Fourier mode is determined by the so-called phase-locked rule, which results in non-self-consistency in the wave numbers. In this paper, a modification is proposed to make it self-consistent. The main idea is that, instead of allowing wave numbers to be complex, all wave numbers are kept real, and the growth or decay of each mode is simply manifested in the growth or decay of the modulus of its shape function. The validity of the new formulation is illustrated by comparing the results with those from the corresponding direct numerical simulation （DNS） as applied to a problem of compressible boundary layer with Mach number 6.

Improved nonlinear parabolized stability equations approach for hypersonic boundary layers

The nonlinear parabolized stability equations(NPSEs)approach is widely used to study the evolution of disturbances in hypersonic boundary layers owing to its high computational efficiency.However,divergence of the NPSEs will occur when disturbances imposed at the inlet no longer play a leading role or when the nonlinear effect becomes very strong.Two major improvements are proposed here to deal with the divergence of the NPSEs.First,all disturbances are divided into two types:dominant waves and non-dominant waves.Disturbances imposed at the inlet or playing a leading role are defined as dominant waves,with all others being defined as non-dominant waves.Second,the streamwise wavenumbers of the non-dominant waves are obtained using the phase-locked method,while those of the dominant waves are obtained using an iterative method.Two reference wavenumbers are introduced in the phase-locked method,and methods for calculating them for different numbers of dominant waves are discussed.Direct numerical simulation(DNS)is performed to verify and validate the predictions of the improved NPSEs in a hypersonic boundary layer on an isothermal swept blunt plate.The results from the improved NPSEs approach are in good agreement with those of DNS,whereas the traditional NPSEs approach is subject to divergence,indicating that the improved NPSEs approach exhibits greater robustness.

DIFFERENCE SCHEME FOR AN INITIAL-BOUNDARY VALUE PROBLEM FOR LINEAR COEFFICIENT-VARIED PARABOLIC DIFFERENTIAL EQUATION WITH A NONSMOOTH BOUNDARY LAYER FUNCTION

In this paper, using nonuniform mesh and exponentially fitted difference method, a uniformly convergent difference scheme for an initial-boundary value problem of linear parabolic differential equation with the nonsmooth boundary layer function with respect to small parameter e is given, and error estimate and numerical result are also given.

BOUNDARY LAYER AND VANISHING DIFFUSION LIMIT FOR NONLINEAR EVOLUTION EQUATIONS

In this paper, we consider an initial-boundary value problem for some nonlinear evolution equations with damping and diffusion. The main purpose is to investigate the boundary layer effect and the convergence rates as the diffusion parameter α goes to zero.

BOUNDARY LAYER ESTIMATION OF THE QUASILINEAR STOKES EQUATIONS

We first prove the existence and uniqueness of solution of quasilinear Stokes equations. Then it is shown when the viscosity vanishes, the solution of the quasilinear Stokes equations tends to the solution of the degenerate equations, in which the viscous term is omitted from the quasilinear Stokes equations and the boundary condition is weakened. In the end, we obtain the boundary layer estimation, Our result shows that the thickness of the boundary layer is proportional to epsilon(1/4).

DEGENERATE BOUNDARY LAYER SOLUTIONS TO THE GENERALIZED BENJAMIN-BONAMAHONY-BURGERS EQUATION

This paper is concerned with the convergence rates of the global solutions of the generalized Benjamin-Bona-Mahony-Burgers（BBM-Burgers） equation to the corresponding degenerate boundary layer solutions in the half-space.It is shown that the convergence rate is t-（α/4） as t →∞ provided that the initial perturbation lies in H α 1 for α 〈 α（q）：= 3 ＋（2/q）,where q is the degeneracy exponent of the flux function.Our analysis is based on the space-time weighted energy method combined with a Hardy type inequality with the best possible constant introduced in [1]

MULTIPLE BOUNDARY LAYER PHENOMENA OF THE SOLUTION OF NONLINEAR HIGHER ORDER ELLIPTIC EQUATIONS WITH PERTURBATION BOTH IN BOUNDARY AND IN OPERATOR

In this paper, we reveal the multiple boundary layer phenomena of the solution of nonlinear higher order elliptic equations with perturbation both in boundary and in operator, and provide a method to find uniformly valid asymptotic solution of arbitrary order for these types of problems.

SIMILARITY SOLUTIONS OF BOUNDARY LAYER EQUATIONS FOR A SPECIAL NON-NEWTONIAN FLUID IN A SPECIAL COORDINATE SYSTME

Two dimensional equations of steady motion for third order fluids are expressed in a special coordinate system generated by the potential flow corresponding to an inviscid fluid. For the inviscid flow around an arbitrary object, the streamlines are the phi-coordinates and velocity potential lines are psi-coordinates which form an orthogonal curvilinear set of coordinates. The outcome, boundary layer equations, is then shown to be independent of the body shape immersed into the flow. As a first approximation,assumption that second grade terms are negligible compared to viscous and third grade terms. Second grade terms spoil scaling transformation which is only transformation leading to similarity solutions for third grade fluid. By using Lie group methods,infinitesimal generators of boundary layer equations are calculated. The equations are transformed into an ordinary differential system. Numerical solutions of outcoming nonlinear differential equations are found by using combination of a Runge-Kutta algorithm and shooting technique.

FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL DIFFUSION-REACTION EQUATIONS OF BOUNDARY LAYER TYPE IN POROUS CATALYST PELLET

In this paper,finite element method(FEM)is used to solve two-dimensional diffu-sion-reaction equations of boundary layer type.This kind of equations are usually too complicatedand diffcult to be solved by applying the traditional methods used in chemical engineering becauseof the steep gradients of concentration and temperature.But,these difficulties are easy to be over-comed when the FEM is used.The integraded steps of solving this kind of problems by the FEMare presented in this paper.By applying the FEM to the two actual examples,the conclusion can bereached that the FEM has the advantages of simplicity and good accuracy.

A Modified Homotopy Analysis Method for Solving Boundary Layer Equations

A new modification of the Homotopy Analysis Method (HAM) is presented for highly nonlinear ODEs on a semi-infinite domain. The main advantage of the modified HAM is that the number of terms in the series solution can be greatly reduced;meanwhile the accuracy of the solution can be well retained. In this way, much less CPU is needed. Two typical examples are used to illustrate the efficiency of the proposed approach.

Unsteady Magnetohydrodynamic Boundary Layer Flow near the Stagnation Point towards a Shrinking Surface

The unsteady two-dimensional, laminar flow of a viscous, incompressible, electrically conducting fluid towards a shrinking surface in the presence of a uniform transverse magnetic field is studied. Taking suitable similarity variables, the governing boundary layer equations are transformed to ordinary differential equations and solved numerically by a perturbation technique for a small magnetic parameter. The effects of various parameters such as unsteadiness parameter, velocity parameter, magnetic parameter, Prandtl number and Eckert number for velocity and temperature distributions along with local Skin friction coefficient and local Nusselt number have been discussed in detail through numerical and graphical representations.

The Boundary Layer Equations and a Dimensional Split Method for Navier-Stokes Equations in Exterior Domain of a Spheroid and Ellipsoid

In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.

NONPARALLEL BOUNDARY LAYER STABILITY IN HIGH SPEED FLOWS

The parabolized stability equations （PSEs） for high speed flows, especially supersonic and hypersonic flows, are derived and used to analyze the nonparallel boundary layer stability. The proposed numerical techniques for solving PSE include the following contents： introducing the efficiently normal transformation of the boundary layer, improving the computational accuracy by using a high-order differential scheme near the wall, employing the predictor-corrector and iterative approach to satisfy the important normalization condition, and implementing the stable spatial marching. Since the second mode dominates the growth of the disturbance in high Mach number flows, it is used in the computation. The evolution and characteristics of the boundary layer stability in the high speed flow are demonstrated in the examples. The effects of the nonparallelizm, the compressibility and the cooling wall on the stability are analyzed. And computational results are in good agreement with the relevant data.

On Nonlinear Evolution of C-type Instability in Nonparallel Boundary Layers

The process of evolution, especially that of nonlinear evolution, of C-type instability of laminar-turbulent flow transition in nonparallel boundary layers are studied by means of a newly developed method called parabolic stability equations （PSE）. Initial conditions, which are very important for the nonlinear problem, are investigated by computing initial solution of the harmonic waves, modifying the mean-flow-distortion, and giving initial value of TS wave and its subharmonic waves at initial station by solving linear PSE. A numerical method with high-order accuracy are developed in the text, the key normalization conditions in the PSE are satisfied, and nonlinear PSE are solved efficiently and implemented stably by the spatial marching. It has been shown that the computed process of nonlinear evolution of C-type instability in Blasius boundary layer is in good agreement with the experimental results.

EVOLUTION ANALYSIS OF TS WAVE AND HIGH-ORDER HARMONIC WAVES IN BOUNDARY LAYERS

Linear and nonlinear evolutions of TS wave and high-order harmonic waves in boundary layers are studied based on the parabolic stability equation （PSE）. Initial conditions are derived by the local method with the Landau expansion. The evolution process and characteristics of the disturbance amplitude and the velocity profile, etc. , especially stronger nonlinear effects, are computed by an efficient numerical method. Effects and regulations of different initial amplitudes, frequencies and pressure gradients on the evolution of disturbances are explored, which are directly relative to the stability and the transition in boundary layers. Simulation results are in good agreement with the data of the accuracy direct numerical simulation （DNS） using full Navier-Stokes equations.