二维平板可压缩边界层的二次稳定性分析

本文在二维可压缩边界层中应用Floquet分析,建立了控制次谐波稳定性的方程组,研究在二维可压缩边界层转捩过程中二维有限振幅的T-S波对三维线性小扰动的作用,并计算了来流马赫数对次谐波的产生和发展情况的影响,从中可以看出二维和三维扰动波相互作用对二维可压缩流动边界层的发展过程所产生的影响.     

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.     

Analysis on linear stability of oblique shock waves in steady supersonic flow

An attached oblique shock wave is generated when a sharp solid projectile flies supersonically in the air. We study the linear stability of oblique shock waves in steady supersonic flow under three dimensional perturbation in the incoming flow. Euler system of equations for isentropic gas model is used. The linear stability is established for shock front with supersonic downstream flow, in addition to the usual entropy condition.     

Stability of transonic shocks for the full Euler system in supersonic flow past a wedge

We study the stability of transonic shocks in steady supersonic flow past a wedge. It is known that in generic case such a problem admits two possible locations of the shock front, connecting the flow ahead of it and behind it. They can be distinguished as supersonic–supersonic shock and supersonic–subsonic shock (or transonic shock). Both these possible shocks satisfy the Rankine–Hugoniot conditions and the entropy condition. We prove that the transonic shock is conditionally stable under perturbation of the upstream flow or perturbation of wedge boundary. Copyright © 2005 John Wiley & Sons, Ltd.     

Wave motions in a spatial boundary layer

The propagation of perturbations in a boundary layer under conditions when the velocity of the approaching stream may be both subsonic and supersonic is considered. With regard to the initial flow in the boundary layer it is assumed that it is stationary and possesses a spatial character which is caused by the external pressure gradient and not by the curvature of the body around which the flow occurs (boundary layers of this kind are extensively used in experiments at the present time). The linearized equations describing waves of vanishingly small amplitude are studied in detail. An analysis of the dispersion relation which links the frequency of the free oscillations with the components of the wave vector reveals a number of special features which are only present in motions with a three-dimensional velocity field. In particular, it is established that the Cauchy problem for the system of linear equations is ill-posed.     

抛物化稳定性方程在可压缩边界层中应用的检验

用抛物化稳定性方程（PSE）,研究了可压缩边界层中扰动的演化,并与由直接数值模拟（DNS）所得进行比较．目的在检验PSE方法用于研究可压缩边界层中扰动演化的可靠性．结果显示,无论是亚音速还是超音速边界层,由PSE方法和由DNS方法所得结果都基本一致,而温度比速度吻合得更好．对超音速边界层,还计算了小扰动的中性曲线．与线性稳定性理论（LST）的结果相比,二者的关系和不可压边界层的情况相似．     

Multidimensional transonic shocks and free boundary problems for nonlinear equations of mixed type

We establish the existence and stability of multidimensional transonic shocks for the Euler equations for steady potential compressible fluids. The Euler equations, consisting of the conservation law of mass and the Bernoulli law for the velocity, can be written as a second-order, nonlinear equation of mixed elliptic-hyperbolic type for the velocity potential. The transonic shock problem can be formulated into the following free boundary problem: The free boundary is the location of the transonic shock which divides the two regions of smooth flow, and the equation is hyperbolic in the upstream region where the smooth perturbed flow is supersonic. We develop a nonlinear approach to deal with such a free boundary problem in order to solve the transonic shock problem. Our results indicate that there exists a unique solution of the free boundary problem such that the equation is always elliptic in the downstream region and the free boundary is smooth, provided that the hyperbolic phase is close to a uniform flow. We prove that the free boundary is stable under the steady perturbation of the hyperbolic phase. We also establish the existence and stability of multidimensional transonic shocks near spherical or circular transonic shocks.

     

Inflow problem for the one-dimensional compressible Navier–Stokes equations under large initial perturbation

This paper is concerned with the inflow problem for the one-dimensional compressible Navier–Stokes equations. For such a problem, Matsumura and Nishihara showed in [10] that there exists boundary layer solution to the inflow problem, and that both the boundary layer solution, the rarefaction wave, and the superposition of boundary layer solution and rarefaction wave are nonlinear stable under small initial perturbation. The main purpose of this paper is to show that similar stability results for the boundary layer solution and the supersonic rarefaction wave still hold for a class of large initial perturbation which can allow the initial density to have large oscillation. The proofs are given by an elementary energy method and the key point is to deduce the desired lower and upper bounds on the density function.     

Transonic shock in a nozzle I: 2D case

In this paper we establish the existence and uniqueness of a transonic shock for the steady flow through a general two‐dimensional nozzle with variable sections. The flow is governed by the inviscid potential equation, and is supersonic upstream, has no‐flow boundary conditions on the nozzle walls, and a given pressure at the exit of the exhaust section. The transonic shock is a free boundary dividing two regions of C flow in the nozzle. The potential equation is hyperbolic upstream where the flow is supersonic, and elliptic in the downstream subsonic region. In particular, our results show that there exists a solution to the corresponding free boundary problem such that the equation is always subsonic in the downstream region of the nozzle when the pressure in the exit of the exhaustion section is appropriately larger than that in the entry. This confirms exactly the conjecture of Courant and Friedrichs on the transonic phenomena in a nozzle [10]. Furthermore, the stability of the transonic shock is also proved when the upstream supersonic flow is a small steady perturbation for the uniform supersonic flow or the pressure at the exit has a small perturbation. The main ingredients of our analysis are a generalized hodograph transformation and multiplier methods for elliptic equation with mixed boundary conditions and corner singularities. © 2004 Wiley Periodicals, Inc.     

Analyzing the longitudinal effect of hypersonic flow past a conical cone via the perturbation method

To analyze the hypersonic flow past a conical cone, the variations of gasdynamic properties subjected to the longitudinal curvature effect by using the perturbation method. An outer perturbation expansion has been carried out by recent researchers, but a problem occurred, the outer expansion solutions are not uniformly valid in the shock layer, however, the outcome near the conical body surface called vortical layer remains deflective. This study intends to discover uniformly valid analytical solutions in the shock layer by applying the inner perturbation expansions matching with the out expansions to analyze the characteristics in the whole region including shock layer and vortical layer. Starting from the zero-order approximate solutions for hypersonic conical flow is then applied as the basic solutions for the outer perturbation expansions of a flow field. The governing equations and boundary conditions are also expanded via outer perturbations. Using an approximate analytical scheme in the derivation process, first-order perturbation equations can be simplified and the approximate closed-form solutions are obtained; furthermore, the various flow field quantities, including the normal force coefficient on the cone surface, have been calculated. According to the variations of gasdynamic properties, the longitudinal curvature effect for the hypersonic flow past a conical cone can be determined. Thicknesses of shock layer and vortical layer can be predicted as well. The physical phenomena inside both layers can be investigated carefully, the conditions for an elliptic cone with longitudinal curvature, m = 1 and n = 2 and other conditions of parameters; the perturbation parameter, ε_m2 = 0.1, semi-vertex angle of the unperturbed cone, δ = 10°, and hypersonic similarity parameter, K_δ = M_∞δ = 1.0, the thickness of vortical layer, η_VL, can be calculated at the position angle of conical cone body, ? = 30° was demonstrated here. Results show how very thin the vortical layer is approximately only 10% of the shock layer close to the body, the pressure in the whole shock layer is verified to be uniformly valid which agrees with previous studies. Large gradient changes in entropy and density are found when the flow approaches the cone surface, the most important is, this method provides a benchmark solution to the hypersonic flow past a conical cone and to assist the grids and numerics for numerical computation should be fashioned to accommodate the whole flow field region including the vortical layer of rapid adjustment, and let the analysis become more effective and low cost.     

Supersonic flow past a flat lattice of cylindrical rods

Two-dimensional supersonic laminar ideal gas flows past a regular flat lattice of identical circular cylinders lying in a plane perpendicular to the free-stream velocity are numerically simulated. The flows are computed by applying a multiblock numerical technique with local boundary-fitted curvilinear grids that have finite regions overlapping the global rectangular grid covering the entire computational domain. Viscous boundary layers are resolved on the local grids by applying the Navier–Stokes equations, while the aerodynamic interference of shock wave structures occurring between the lattice elements is described by the Euler equations. In the overlapping grid regions, the functions are interpolated to the grid interfaces. The regimes of supersonic lattice flow are classified. The parameter ranges in which the steady flow around the lattice is not unique are detected, and the mechanisms of hysteresis phenomena are examined.     

Steady transonic shocks and free boundary problems in infinite cylinders for the Euler equations

We establish the existence and stability of multidimensional transonic shocks (hyperbolic‐elliptic shocks) for the Euler equations for steady compressible potential fluids in infinite cylinders. The Euler equations, consisting of the conservation law of mass and the Bernoulli law for velocity, can be written as a second order nonlinear equation of mixed elliptic‐hyperbolic type for the velocity potential. The transonic shock problem in an infinite cylinder can be formulated into the following free boundary problem: The free boundary is the location of the multidimensional transonic shock which divides two regions of C^1,α flow in the infinite cylinder, and the equation is hyperbolic in the upstream region where the C^1,α perturbed flow is supersonic. We develop a nonlinear approach to deal with such a free boundary problem in order to solve the transonic shock problem in unbounded domains. Our results indicate that there exists a solution of the free boundary problem such that the equation is always elliptic in the unbounded downstream region, the uniform velocity state at infinity in the downstream direction is uniquely determined by the given hyperbolic phase, and the free boundary is C^1,α, provided that the hyperbolic phase is close in C^1,α to a uniform flow. We further prove that, if the steady perturbation of the hyperbolic phase is C^2,α, the free boundary is C^2,α and stable under the steady perturbation. © 2003 Wiley Periodicals Inc.     

Stability of a Mach configuration

We prove the stability of a Mach configuration, which occurs in shock reflection off an obstacle or shock interaction in compressible flow. The compressible flow is described by a full, steady Euler system of gas dynamics. The unperturbed Mach configuration is composed of three straight shock lines and a slip line carrying contact discontinuity. Among four regions divided by these four lines in the neighborhood of the intersection, two are supersonic regions, and other two are subsonic regions. We prove that if the constant states in the supersonic regions are slightly perturbed, then the structure of the whole configuration holds, while the other two shock fronts and the slip line, as well as the flow field in the subsonic regions, are also slightly perturbed. Such a conclusion asserts the existence and stability of the general Mach configuration in shock theory. In order to prove the result, we reduce the problem to a free boundary value problem, where two unknown shock fronts are free boundaries, while the slip line is transformed to a fixed line by a Lagrange transformation. In the region where the solution is to be determined, we have to deal with an elliptic‐hyperbolic composed system. By decoupling this system and combining the technique for both hyperbolic equations and elliptic equations, we establish the required estimates, which are crucial in the proof of the existence of a solution to the free boundary value problem. © 2005 Wiley Periodicals, Inc.     

The antiplane problem of the motion of a point shear load over a two-layer elastic base at a constant velocity

The problem of modelling the motion of a force disturbance in an elastic medium that is heterogeneous over its depth is investigated. It is in an antiplane formulation in a moving system of coordinates that all possible versions of the ratio of the velocity of motion of the surface point shear load to the velocities of the shear waves in the layers of the two-layer elastic base are examined. Cases of a subsonic regime (SBR) in the upper and lower layers, of a supersonic regime (SPR) in the upper layer and an SBR in the lower layer, and of an SBR in the upper layer and an SPR in the lower layer are studied using the Fourier transform and the theory of residues. The last two cases are extremely interesting from the mathematical point of view, as here, on the boundary between the layers, the solutions of elliptic and hyperbolic equations meet, and previously unknown features arise in the displacements that,it seems, should also occur in the solution of the corresponding plane problem. The case of an SPR in the upper and lower layers is investigated using a special method for successive allowance for the incident, reflected and refracted shock wave fronts. In all cases, expressions are obtained for the displacements in the layers, and their characteristic features are investigated.     

Oblique Shock‐Vortex Interaction

Breakdown of a slender vortex caused by an oblique (OSVI) shock is studied using numerical solutions of the Euler equations for unsteady three‐dimensional flow. A Burgers vortex with a given circulation and axial velocity distribution is prescribed at the inflow boundary. The calculations show that like in incompressible flows supersonic breakdown is primarily controlled by pressure forces. The shock is deformed into an ‘s’‐shaped part near the vortex core where the shock becomes normal. The results indicate that initiation of breakdown is more sensitive to variations in the axial velocity than in the circulation, and that the .ow structure is clearly time‐dependent.     

Hamiltonian long‐wave expansions for free surfaces and interfaces

The theory of internal waves between two bodies of immiscible fluid is important both for its interest to ocean engineering and as a source of numerous interesting mathematical model equations that exhibit nonlinearity and dispersion. In this paper we derive a Hamiltonian formulation of the problem of a dynamic free interface (with rigid lid upper boundary conditions), and of a free surface and a free interface, this latter situation occurring more commonly in experiment and in nature. From the formulation, we develop a Hamiltonian perturbation theory for the long‐wave limits, and we carry out a systematic analysis of the principal long‐wave scaling regimes. This analysis provides a uniform treatment of the classical works of Peters and Stoker (28), Benjamin (3, 4), Ono (26), and many others. Our considerations include the Boussinesq and Korteweg–de Vries (KdV) regimes over finite‐depth fluids, the Benjamin‐Ono regimes in the situation in which one fluid layer is infinitely deep, and the intermediate long‐wave regimes. In addition, we describe a novel class of scaling regimes of the problem, in which the amplitude of the interface disturbance is of the same order as the mean fluid depth, and the characteristic small parameter corresponds to the slope of the interface. Our principal results are that we highlight the discrepancies between the case of rigid lid and of free surface upper boundary conditions, which in some circumstances can be significant. Motivated by the recent results of Choi and Camassa (6, 7), we also derive novel systems of nonlinear dispersive long‐wave equations in the large‐amplitude, small‐slope regime. Our formulation of the dynamical free‐surface, free‐interface problem is shown to be very effective for perturbation calculations; in addition, it holds promise as a basis for numerical simulations. © 2005 Wiley Periodicals, Inc.     

Asymptotic Approach to the Problem of Boundary Layer Instability in Transonic Flow

Tollmien–Schlichting waves can be analyzed using the Prandtl equations involving selfinduced pressure. This circumstance was used as a starting point to examine the properties of the dispersion relation and the eigenmode spectrum, which includes modes with amplitudes increasing with time. The fact that the asymptotic equations for a nonclassical boundary layer (near the lower branch of the neutral curve) have unstable fluctuation solutions is well known in the case of subsonic and transonic flows. At the same time, similar solutions for supersonic external flows do not contain unstable modes. The bifurcation pattern of the behavior of dispersion curves in complex domains gives a mathematical explanation of the sharp change in the stability properties occurring in the transonic range.     

Transonic shocks in 3-D compressible flow passing a duct with a general section for Euler systems

This paper is devoted to the study of a transonic shock in three-dimensional steady compressible flow passing a duct with a general section. The flow is described by the steady full Euler system, which is purely hyperbolic in the supersonic region and is of elliptic-hyperbolic type in the subsonic region. The upstream flow at the entrance of the duct is a uniform supersonic one adding a three-dimensional perturbation, while the pressure of the downstream flow at the exit of the duct is assigned apart from a constant difference. The problem to determine the transonic shock and the flow behind the shock is reduced to a free boundary value problem of an elliptic-hyperbolic system. The new ingredients of our paper contain the decomposition of the elliptic-hyperbolic system, the determination of the shock front by a pair of partial differential equations coupled with the three-dimensional Euler system, and the regularity analysis of solutions to the boundary value problems introduced in our discussion.

     

Uniform convergence of a finite difference scheme for a system of coupledreaction‐diffusion equations

We study a system of coupled reaction‐diffusion equations. The equations have diffusion parameters of different magnitudes associated with them. Near each boundary, their solution exhibit two overlapping layers. A difference scheme on layer‐adapted piecewise uniform meshes is used to solve the system numerically. We show that the scheme is almost second‐order convergent, uniformly in both perturbation parameters, thus improving previous results [3].