零攻角小钝头钝锥高超音速绕流边界层的稳定性分析和转捩预报

研究了零攻角小钝头圆锥高超音速边界层的稳定性及转捩预测问题．小钝头的球头半径为0.5 mm,锥的半锥角为5°,来流马赫数为6．采用直接数值模拟方法得到了钝锥的基本流场,利用线性稳定性理论分析了等温壁面和绝热壁面条件下的第一、第二模态不稳定波,并用“e-N”方法对转捩位置进行了预测．在没有实验给出N值的情况下,暂取N为10．研究发现,壁面温度条件对于转捩位置有较大影响．绝热边界层的转捩位置比等温边界层的靠后．且尽管高马赫数下第二模态波的最大增长率远大于第一模态波的最大增长率,但绝热边界层的转捩位置是由第一模态不稳定波决定的．研究方法应能推广到有攻角的三维边界层流动的转捩预测．     

自由来流湍流与三维壁面局部粗糙诱导平板边界层不稳定T-S波的数值研究

采用直接数值模拟（DNS）方法,研究了在自由来流湍流与三维壁面局部粗糙作用下平板边界层内诱导产生不稳定T-S波的物理问题.数值结果可知,在平板边界层内发现了二维和三维T-S波组成的波包空间序列以及求得了波包向前传播的群速度大小,从而证明了自由来流湍流与三维壁面局部粗糙作用是激励平板边界层内诱导产生不稳定T-S波的一种机制.随后,建立了平板边界层内被激发的二维和三维T S波的初始幅值与自由来流湍流度,三维壁面局部粗糙的流向长度、展向宽度及法向高度之间的关系.这一问题的深入研究,进一步完善了流动稳定性与湍流理论.     

PSE在超音速边界层二次失稳问题中的应用

用抛物化稳定性方程（PSE）研究超音速边界层中的二次失稳问题．结果显示无论二维基本扰动是第一模态还是第二模态的T-S波,二次失稳机制都起作用．三维亚谐波的放大率随其展向波数和二维基本波幅值的变化关系与不可压缩边界层中所得类似．但是,即使二维波的幅值大到2％的量级,三维亚谐波的最大放大率仍远小于最不稳定的第二模态二维T-S波的放大率．因此,二次失稳应该不是导致超音速边界层转捩的主要因素．     

Resonant three-wave interactions in non-linear hyper-elastic fluid-filled tubes

A multiple-scales method is used to derive the Three-Wave Interaction (TWI) equations describing resonantly interacting triads in nonlinear hyper-elastic fluid-filled tubes. The tube wall is assumed to be an axially-tethered nonlinear membraneous cylindrical shell for which the resultant stresses can be determined by a strain-energy functional. The fluid within the tube is assumed to be two-dimensional, axi-symmetric and inviscid. We show that small-but-finite amplitude strongly dispersive pressure wave packets can continuously exchange energy in a resonant triad while conserving total energy. For a Mooney-Rivlin shell wall the theory presented predicts a short wavelength cutoff on the order of the tube radius. Thus pressure pulses containing wavelengths on the order of the tube radius and longer may contain resonantly interacting modes. Special solutions are presented: temporally developing modes, pump-wave approximations and explosively unstable steadily-traveling wave packets.     

Three-dimensional Tollmien-Schlichting waves generated by sound in the boundary layer on an elastic surface at transonic free-stream velocities

Within the framework of the triple-deck theory, the effect of surface elasticity on three-dimensional packets of Tollmien-Schlichting waves generated by acoustic disturbances induced near the boundary layer at transonic free-stream velocities is investigated. It is shown that the elasticity of the surface considerably weakens the most unstable oblique waves but does not change the characteristic horseshoe shape of wave packets with two disturbance peaks propagating at an angle to the incoming flow.     

On the dissipation due to wave ringing in nonelliptic elastic materials

Summary Initial-boundary value problems describing the mechanics of nonelliptic elastic materials give rise to solutions that involve phase boundaries, the motion of which can dissipate mechanical energy. We investigate whether this dissipation, acting alone, can drive such a system toward equilibrium. Moving phase boundaries are regarded as a localized dissipative mechanism, and we consider a model which specifically excludes dissipation away from a phase boundary (such as that due to viscoelastic damping). In the problem under consideration, wave packets reverberate between the fixed external boundary and a single internal phase boundary. The phase boundary remains stationary unless it is acted upon by one of these wave packets, and each such interaction dissipates a finite amount of energy while causing the initiating wave packet to split into a reflected wave packet and a transmitted wave packet. Consequently, the number of wave packets increases in a geometric fashion. Each individual interaction of a wave packet with the phase boundary is, in a certain sense, mechanically underdetermined, and we augment the mechanical theory with two alternative energy criteria, each of which determines a different interaction dynamics. These alternative energy criteria are motivated by considerations of maximizing the energy dissipation in the system. We treat a system that is perturbed out of an initial minimum energy equilibrium state by a disturbance at the external boundary. A framework is developed for treating the resulting wave reverberations and calculating the energy dissipation for large time. Numerical computation indicates that the total energy dissipated in both versions of the dynamical problem is that which is necessary to settle into a new energy-minimal equilibrium state. We then establish the same result analytically for a meaningful limit involving a vanishingly small dynamical perturbation.     

Internal Kelvin waves in a two-layer liquid model

The subinertial internal Kelvin wave solutions of a linearized system of the ocean dynamics equations for a semi-infinite two-layer f-plane model basin of constant depth bordering a straight, vertical coast are imposed. A rigid lid surface condition and no-slip wall boundary condition are imposed. Some trapped wave equations are presented and approximate solutions using an asymptotic method are constructed. In the absence of bottom friction, the solution consists of a frictionally modified Kelvin wave and a vertical viscous boundary layer. With a no-slip bottom boundary condition, the solution consists of a modified Kelvin wave, two vertical viscous boundary layers, and a large cross-section scale component. The numerical solutions for Kelvin waves are obtained for model parameters that take account of a joint effect of lateral viscosity, bottom friction, and friction between the layers.     

Logarithmic lower bounds for Néel walls

Most mathematical models for interfaces and transition layers in materials science exhibit sharply localized and rapidly decaying transition profiles. We show that this behavior can largely change when non-local interactions dominate and internal length scales fail to be determined by dimensional analysis: we consider a reduced model for Néel walls, micromagnetic transition layers which are observed in a suitable thin-film regime. The typical phenomenon associated with this wall type is the very long logarithmic tail of transition profiles. Recently, we derived logarithmic upper bounds. Here, we prove that the latter result is indeed optimal. In particular, we show that Néel wall profiles are supported by explicitly known comparison profiles that minimize relaxed variational principles and exhibit logarithmic decay behavior. This lower bound is established by a comparison argument based on a global maximum principle for the non-local field operator and the qualitative decay behavior of comparison profiles.Received: 17 June 2003, Accepted: 18 November 2003, Published online: 25 February 2004Mathematics Subject Classification (2000): 78A30, 49S05, 45G15, 35B25     

The behavior of plasma with an arbitrary degree of degeneracy of electron gas in the conductive layer

We obtain an analytic solution of the boundary problem for the behavior (fluctuations) of an electron plasma with an arbitrary degree of degeneracy of the electron gas in the conductive layer in an external electric field. We use the kinetic Vlasov–Boltzmann equation with the Bhatnagar–Gross–Krook collision integral and the Maxwell equation for the electric field. We use the mirror boundary conditions for the reflections of electrons from the layer boundary. The boundary problem reduces to a one-dimensional problem with a single velocity. For this, we use the method of consecutive approximations, linearization of the equations with respect to the absolute distribution of the Fermi–Dirac electrons, and the conservation law for the number of particles. Separation of variables then helps reduce the problem equations to a characteristic system of equations. In the space of generalized functions, we find the eigensolutions of the initial system, which correspond to the continuous spectrum (Van Kampen mode). Solving the dispersion equation, we then find the eigensolutions corresponding to the adjoint and discrete spectra (Drude and Debye modes). We then construct the general solution of the boundary problem by decomposing it into the eigensolutions. The coefficients of the decomposition are given by the boundary conditions. This allows obtaining the decompositions of the distribution function and the electric field in explicit form.     

The Nonparallel Evolution of Nonlinear Short Waves in Buoyant Boundary Layers

Buoyant boundary-layer flows, typified by the flow over a heated flat plate, have the curious property that they can exhibit regions of "overshoot" in which the streamwise velocity exceeds its free-stream value. A consequence of this is the streamwise velocity develops a local maximum and is inflectional in nature. It is therefore inviscidly unstable, and the fastest growing wave mode is known to be one whose wavelength is short compared to the boundary-layer thickness. In this work we consider the nonparallel evolution of these short waves and show that they can be described in terms of the solution of a system of ordinary differential equations. Numerical and asymptotic studies enable us to explain the ultimate fate of the wave and show, depending on a key parameter which is a function of the underlying boundary layer, that two possibilities can arise. Nonparallelism may be sufficiently stabilizing so as to extinguish the linearly unstable waves or, in other cases, the mode may intensify but concentrate itself in a very thin zone surrounding the maximum in the streamwise velocity. These findings enable us to give some indication of the part these modes play in the transition to turbulence in buoyant boundary layers.     

Stable boundary element domain decomposition methods for the Helmholtz equation

In this paper, we present a stable boundary element domain decomposition method to solve boundary value problems of the Helmholtz equation via a tearing and interconnecting approach. A possible non-uniqueness of the solution of local boundary value problems due to the appearance of local eigensolutions is resolved by using modified interface conditions of Robin type, which results in a Galerkin boundary element discretization which is robust for all local wave numbers. Numerical examples confirm the stability of the proposed approach.     

On the Generation of Mean Flows by the Interaction of Görtler Vortices and Tollmien-Schlichting Waves in Curved Channel Flows

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. Here the interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmien-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; here the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it is found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance the mean state driven by a vortex of short wavelength in the absence of a Tollmien-Schlichting wave is considered. It is shown that for a given channel curvature and vortex wavelength there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When Tollmien-Schlichting waves are present then the nonlinear differential equation to determine the mean state is modified. At sufficiently high Tollmien-Schlichting amplitudes it is found that the vortex flows are destroyed, but there is a range of amplitudes where a fully nonlinear mixed vortex-wave state exists and indeed drives a mean state having little similarity with the flow which occurs without the instability modes. The vortex and Tollmien-Schlichting wave structure in the nonlinear regime has viscous wall layers and internal shear layers; the thickness of the internal layers is found to be a function of the Tollmien-Schlichting wave amplitude.     

Nonlinear evolution of a first mode oblique wave in a compressible boundary layer: I. Heated/cooled walls

The nonlinear stability of an oblique mode propagating in atwo-dimensional compressible boundary layer is considered underthe long wavelength approximation. The growth rate of the waveis assumed to be small so that the ideas of unsteady nonlinearcritical layers can be applied. It is shown that the spatial/temporalevolution of the mode is governed by a pair of coupled unsteadynonlinear equations for the disturbance vorticity and density.Expressions for the linear growth rate show clearly the effectsof wall heating and cooling, and in particular how heating destabilizesthe boundary layer for these long wavelength inviscid modesat O(1) Mach numbers. A generalized expression for the lineargrowth rate is obtained and is shown to compare very well fora range of frequencies and wave angles at moderate Mach numberswith full numerical solutions of the linear stability problem.The numerical solution of the nonlinear unsteady critical layerproblem using a novel method based on Fourier decompositionand Chebyshev collocation is discussed and some results arepresented.     

Implementation of near-wall boundary conditions for modeling boundary layers with free-stream turbulence

The paper is devoted to the extension of the near-wall domain decomposition, earlier developed in some previous works by the authors, to modeling flat-plate boundary layers undergoing laminar-to-turbulent bypass transition. The steady-state wall boundary layers at high-intensity free-stream turbulence are studied on the basis of differential turbulence models with the use of non-overlapping domain decomposition. In the approach the near-wall resolution is replaced by the interface boundary conditions of Robin type. In contrast to the previous studies, the main attention is paid to the laminar–turbulent transitional regime. With the use of modified turbulence models we study an effect of free-stream parameters on the development of dynamic processes in the boundary layer including a transitional regime and fully developed turbulent flow. In addition, for the first time a full scale domain decomposition is realized via iterations between the inner and outer subregions until a convergence. The computational profiles of the velocity and intensity of the turbulence kinetic energy are compared with experimental data. A possible range of location of the near-wall interface boundary is found.     

The Propagation of Finite-Amplitude Waves in a Model Boundary Layer

Finite-amplitude wave propagation is considered in flows of boundary-layer type when the wavelength is long compared to the boundary layer thickness. In this limit, the evolution of the amplitude is governed by the Benjamin-Ono equation and we have computed the coefficients of its nonlinear and dispersive terms for the specific case of Tietjens's model. The propagation of wave packets is also considered, and it is found that for packets centered about an O(1) wavenumber questions again arise relative to long waves, except that now the packet-induced mean flow is the “long wave.” By introducing an appropriate scaling for the far field and employing multiple scales in the direction transverse to the flow, it is shown how the mean-flow distortion can be made to vanish at infinity.     

Periodic solutions for the Schrödinger equation with nonlocal smoothing nonlinearities in higher dimension

We consider the nonlinear Schrödinger equation in higher dimension with Dirichlet boundary conditions and with a nonlocal smoothing nonlinearity. We prove the existence of small amplitude periodic solutions. In the fully resonant case we find solutions which at leading order are wave packets, in the sense that they continue linear solutions with an arbitrarily large number of resonant modes. The main difficulty in the proof consists in a “small divisor problem” which we solve by using a renormalisation group approach.     

Spatial Stability for Mixed Convection Boundary Layer over a Heated Horizontal Plate

The spatial stability properties of a mixed convection boundary layer developing over a heated horizontal plate is studied here under linear and quasi-parallel flow assumption. The main aim of the present work is to find out if there is a critical buoyancy parameter that would indicate the importance of heat transfer in destabilizing mixed convection boundary layers, when the buoyancy effect is given by Boussinesq approximation. The undisturbed flow used here is that given by the similarity solution of [ 1 ] that implies the wall temperature to vary as the inverse square root of the distance from the leading edge of the plate. The stability of this flow has been investigated by using the compound matrix method (CMM)—that allows finding all the modes in the chosen range in the complex wave number plane for spatial stability analysis. Presented neutral curves for mixed convection boundary layer show the existence of two types of disturbances present simultaneously, for large buoyancy parameter. One notices very unstable high-frequency mode when the buoyancy parameter exceeds the above-mentioned critical value. This unstable thermal mode is in addition to the hydrodynamic mode of isothermal flow given by corresponding similarity profile. The calculated critical buoyancy parameter is shown to qualitatively match with experimental results.     

Investigation of Dispersion Characteristics of Composites

This numerical study of dispersion characteristics of some symmetric and asymmetric composites modelled as multilayered packets of layers with arbitrary anisotropy of each layer is carried out. The authors introduce a subsidiary boundary problem of three-dimensional theory of elasticity for the system of partial differential equations describing the harmonic oscillations of the composite excited by a surface load. The authors obtained the curves and surfaces describing slowness of the wave front of Lamb waves propagating with respect to the direction of propagation. The method suggested by the authors of the present article does not impose any restrictions on the type of elastic anisotropy of the layers, their orientations and mutual positions. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rogue waves in baroclinic flows

We investigate an AB system, which can be used to describe marginally unstable baroclinic wave packets in a geophysical fluid. Using the generalized Darboux transformation, we obtain higher-order rogue wave solutions and analyze rogue wave propagation and interaction. We obtain bright rogue waves with one and two peaks. For the wave packet amplitude and the mean-flow correction resulting from the self-rectification of the nonlinear wave, the positions and values of the wave crests and troughs are expressed in terms of a parameter describing the state of the basic flow, in terms of a parameter responsible for the interaction of the wave packet and the mean flow, and in terms of the group velocity. We show that the interaction of the wave packet and mean flow and also the group velocity affect the propagation and interaction of the amplitude of the wave packet and the self-rectification of the nonlinear wave.     

Influence of Finite Amplitude Disturbances on the Nonstationary Modes of a Compressible Boundary Layer Flow

A weakly nonlinear stability analysis is performed to search for the effects of compressibility on a mode of instability of the three-dimensional boundary layer flow due to a rotating disk. The motivation is to extend the stationary work of [ 1 ] (hereafter referred to as S90) to incorporate into the nonstationary mode so that it will be investigated whether the finite amplitude destabilization of the boundary layer is owing to this mode or the mode of S90. Therefore, the basic compressible flow obtained in the large Reynolds number limit is perturbed by disturbances that are nonlinear and also time dependent. In this connection, the effects of nonlinearity are explored allowing the finite amplitude growth of a disturbance close to the neutral location and thus, a finite amplitude equation governing the evolution of the nonlinear lower branch modes is obtained. The coefficients of this evolution equation clearly demonstrate that the nonlinearity is destabilizing for all the modes, the effect of which is higher for the nonstationary waves as compared to the stationary waves. Some modes particularly having positive frequency, regardless of the adiabatic or wall heating/cooling conditions, are always found to be unstable, which are apparently more important than those stationary modes determined in S90. The solution of the asymptotic amplitude equation reveals that compressibility as the local Mach number increases, has the influence of stabilization by requiring smaller initial amplitude of the disturbance for the laminar rotating disk boundary layer flow to become unstable. Apart from the already unstable positive frequency waves, perturbations with positive frequency are always seen to compete to lead the solution to unstable state before the negative frequency waves do. Also, cooling the surface of the disk will be apparently ineffective to suppress the instability mechanisms operating in this boundary layer flow.