大气边界层内羽流扩散研究

主要研究了高架点源的污染,即``羽流扩散'. 由于大气边界层内的湍流运动,引 起污染物扩散的主要因素是湍流扩散, 研究湍流扩散有两种基本方法：统计理论和梯度理论. 采用拉格朗日粒子扩散模型 (Lagrangian particle dispersion model, LPDM)来评价羽流扩散的平均浓度,分别计算了表征扩散程度的3个参数：沿烟轴方向的 着地浓度、水平方向的羽流扩散宽度和垂直方向的羽流扩散宽度. 从计算结果和风洞实验的 数据对比来看,用该模式可以较好地模拟大气边界层内复杂气候条件下的羽流扩散.     

基于雷诺应力模型的平衡大气边界层模拟

基于雷诺应力湍流模型（简称RSM模型）,研究了平衡大气边界层风场数值模拟问题．假设流体不可压,且不计雷诺应力输运方程中的对流项、浮力产生项、系统旋转产生项和扩散项,在准各向同性的条件下,推导出RSM模型湍动能k的表达式是标准k-ε模型k常数表达式的0.893倍．考虑k沿高度变化的修正,根据在标准k-ε模型中满足水平均匀性的湍流来流边界条件,提出在RSM模型中产生平衡大气边界层的湍流来流边界条件．基于空风洞的数值模拟结果表明,与工程上常用的湍流来流边界条件相比,基于本文提出的湍流来流边界条件得到的风场水平均匀性更优,且在整个流域内,得到的雷诺应力剖面更合适．从而验证了该湍流来流边界条件的适用性．     

大气边界层的一些空气动力学特征

本文在定常和中性温度层结条件下,首先介绍了均匀平坦地面上大气边界层的一些空气动力学特征,并将大气边界层与一般空气动力学边界层作了比较。然后介绍了地面粗糙度分布有阶跃变化以及非平坦地形上的大气边界层特征。文中说明,气流的平均流场结构以及湍流特性与均匀平坦地面的情形相比有显著不同。最后,提出了一些有待今后进一步研究的问题。      

化学非平衡边界层壁面边界条件研究

通过理论分析导出了化学非平衡边界层，化学反应壁面边界条件的一般性提法，并对其具体应用进行了研究．在壁面完全催化时，可利用传递系数法使烧蚀壁面条件计算与边界层计算解耦，给出了一般计算方法．     

在1.4米航空风洞中模拟大气边界层

在风洞中正确模拟大气边界层的流动特性是风工程风洞试验结果可信的必要条件,本次试验研究的目的是在短试验段的航空风洞中建立大比例的大气边界层模拟流场。通过适当的方式延长1.4m×1.4m航空风洞的试验段长度,并利用尖塔、粗糙元等边界层发生装置,在该风洞中建立了边界层流场,测量了流场的平均风速剖面、湍流强度剖面、脉动风速的自相关系数、风谱等参数,讨论了湍流积分尺度的处理和大气边界层几何模拟比例的确定,用谱拟合法和自相关系数积分法求出了湍流积分尺度。结果分析表明:试验所得流场是合理的大气边界层模拟流场,其平均风速剖面幂指数α=0.3,大气边界层模拟比例为1∶500,为后续的建筑物模型动态风荷载试验提供了前提条件。     

流过泰氟龙烧蚀表面的化学非平衡边界层数值研究

对流过泰氟龙烧蚀表面的化学非平衡边界层进行了数值分析．研究了气体模型、化学反应速率常数和壁画催化对边界层特性的影响．研究表明，气体模型对烧蚀产物组元浓度大小的排序有很大影响，但不同气体模型造成的电子数密度差最大仅一个数量级；不同的氟碳反应速率常数，虽在非催化壁条件下，对组元浓度剖面有强烈影响，但对电子数密度剖面，不管壁面催化特性如何，几乎都没有影响；表面烧蚀可能引起边界层分离．     

大气边界层的风洞模拟

本文通过比较大气边界层与不同比尺的模型实验，研究了风洞模型实验中的比尺效应。结果表明，大气边界层粗糙高度的比尺应当和边界层厚度比尺相等。边界层厚度的比尺在２００～４００之间时可以得到最好的速度分布的相似性，表面压强分布主要依赖于边界层比尺而与建筑物几何比尺关系不大。如果正确选择了粗糙高度比尺，扩散系数比尺和边界层厚度比尺相同。     

流过泰氟龙烧蚀表面的化学非平衡边界层数值研究

对流过泰氟龙烧蚀表面的化学非平衡边界层进行了数值分析．研究了气体模型、化学反应速率常数和壁画催化对边界层特性的影响．研究表明，气体模型对烧蚀产物组元浓度大小的排序有很大影响，但不同气体模型造成的电子数密度差最大仅一个数量级；不同的氟碳反应速率常数，虽在非催化壁条件下，对组元浓度剖面有强烈影响，但对电子数密度剖面，不管壁面催化特性如何，几乎都没有影响；表面烧蚀可能引起边界层分离．     

倾斜吹吸控制下湍流边界层减阻的直接数值模拟

雷诺切应力是壁湍流高摩擦阻力的重要来源, 有理论认为可以通过壁面生成负雷诺应力(数值上为正)的方式来削弱湍流流场中雷诺应力的分布, 以此获得流动减阻. 而通过对雷诺平均运动方程的法向二次积分, 可以发现壁面生成正雷诺应力(数值上为负)对壁面摩擦阻力系数才有负贡献. 文中在湍流边界层流动的控制区域下边界设置一系列倾斜狭缝, 利用该装置通过周期性吹吸的方法产生壁面生成正(负)雷诺应力, 并采用直接数值模拟方法考察和验证上文提到的减阻理论. 文中采用的湍流边界层流动模型, 其流动雷诺数(基于外流速度及动量损失厚度)从300 发展到860. 文中通过多组数值模拟算例, 考察了射流强度和频率对壁面摩擦阻力系数的影响, 并对比了壁面生成正或负雷诺应力对流动的影响. 研究表明, 壁面生成正雷诺应力控制的减阻率能达到3.26, 而壁面生成负雷诺应力控制的减阻效果较壁面生成正雷诺应力控制的要差; 壁面生成的正雷诺应力对壁面摩擦阻力有负贡献, 而壁面生成的负雷诺应力对壁面摩擦阻力有正贡献; 通过考察控制的收支比, 发现控制方案不能获得能量净收益.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

Stratified Atmospheric Boundary Layers and Breakdown of Models

The goal of this study is to assess complications in atmospheric stable boundary layers which are not included in numerical models of the stably stratified boundary layer and to provide a formulation of surface fluxes for use in numerical models. Based on an extensive interpretive literature survey and new eddy correlation data for the stable boundary layer, this study defines two prototype stable boundary layers: the weakly stable case and the very stable case. The weakly stable boundary layer is amenable to existing models. The very stable boundary layer eludes modeling attempts due to breakdown of existing formulations of turbulence and due to features found in the atmosphere which are not normally included in models. The latter includes clear-air radiative cooling, low-level jets, surface heterogeneity, gravity waves, meandering motions, and other mesoscale motions which propagate from outside the local domain. While these mechanisms are not essential to understanding idealized or laboratory versions of the stable boundary layer, they complicate comparisons of numerical models and theories with actual atmospheric boundary layers. Statistics which describe various features of the stable boundary layer are offered for future comparison with modeling results. Received 13 June 1997 and accepted 19 November 1997     

A Very Large Eddy Simulation Model Using a Reductionist Inlet Turbulence Generator and Wall Modeling for Stable Atmospheric Boundary Layers

Fluid Dynamics - Despite many advances in numerical simulation of stable boundary layers (SBL), most of the models developed are complex and computationally expensive. A computational fluid...     

Thermosolutal Marangoni Forced Convection Boundary Layers

An analysis is made for the forced convection thermal and solute concentration Marangoni boundary layers (thin dissipative layers) that can be formed along the surface, which separates two immiscible fluids in surface driven flows when the appropriately defined Reynolds number is large enough. Similarity equations for the case in which an external pressure gradient is imposed are derived. These equations are then solved numerically for some values of the involved parameters using very efficient numerical schemes known as Keller-box and superposition methods and the features of the flow and transport characteristics are analysed and discussed.     

Wall Boundary Layers for Maxwell Liquids

Flows of viscoelastic liquids at high Weissenberg number exhibit stress boundary layers near walls. These boundary layers are caused by the memory of the fluid: while particles at the wall remain in their position, particles at some distance from the wall move a long distance within one relaxation time if the Weissenberg number is high. Since the stresses depend on the flow history, this causes a steep boundary layer to form. A rescaling of the variables exploiting the thinness of this boundary layer can be used to derive a reduced set of boundary layer equations. This paper addresses the question of existence of solutions for these boundary layer equations. Using an implicit function argument, we prove the existence of a large class of solutions which arise from spatially periodic perturbations of uniform shear flow. The solutions we find can be characterized by the shear rate outside the boundary layer, which can be prescribed arbitrarily. Accepted: September 27, 1999     

A New Insight into the Convective Boundary Condition

The steady mixed convection boundary layer flows over a vertical surface adjacent to a Darcy porous medium and subject respectively to (i) a prescribed constant wall temperature, (ii) a prescribed variable heat flux, $q_\mathrm{w} =q_0 x^{-1/2}$ q w = q 0 x ? 1 / 2 , and (iii) a convective boundary condition are compared to each other in this article. It is shown that, in the characteristic plane spanned by the dimensionless flow velocity at the wall ${f}^{\prime }(0)\equiv \lambda $ f ′ ( 0 ) ≡ λ and the dimensionless wall shear stress $f^{\prime \prime }(0)\equiv S$ f ′ ′ ( 0 ) ≡ S , every solution $(\lambda , S)$ ( λ , S ) of one of these three flow problems at the same time is also a solution of the other two ones. There also turns out that with respect to the governing mixed convection and surface heat transfer parameters $\varepsilon $ ε and $\gamma $ γ , every solution $(\lambda , S)$ ( λ , S ) of the flow problem (iii) is infinitely degenerate. Specifically, to the very same flow solution $(\lambda , S)$ ( λ , S ) there corresponds a whole continuous set of values of $\varepsilon $ ε and $\gamma $ γ which satisfy the equation $S=-\gamma (1+\varepsilon -\lambda )$ S = ? γ ( 1 + ε ? λ ) . For the temperature solutions, however, the infinite degeneracy of the velocity solutions becomes lifted. These and further outstanding features of the convective problem (iii) are discussed in the article in some detail.     

Nonlinear Wave Motions in Stratified Boundary Layers

We consider nonlinear wave motions in strongly buoyant mixed forced–free convection boundary layer flows. In the natural limit of large Reynolds number the nonlinear evolution of a single monochromatic wave mode is shown to be governed by a novel wave/mean-flow interaction in which the wave amplitude and the wave induced mean-flow are of comparable size. A nonlinear integral equation describing the bifurcation to finite-amplitude travelling wave solutions is derived. Solutions of this equation are presented together with a discussion of their physical significance. Received 10 December 1996 and accepted 14 April 1997     

On Disturbance Propagation in Internal-Flow Boundary Layers

The unsteady interaction of plane-channel wall boundary layers with a supersonic inviscid flow is investigated. The flow regimes in which disturbances introduced by the boundary layer developing on one wall influence the boundary layer on the other wall are considered. The regime of relatively large pressure disturbance amplitudes generated near the nozzle outlet or by deforming the channel walls is studied. In these conditions, the interaction process is described by a system of Burgers equations with retarded arguments. Numerical solutions of this system are obtained for symmetric and antisymmetric perturbations of the channel walls.     

Bending Boundary Layers in a Moving Strip

The high-speed transport of thin-sheet materials occurs in avariety of industrial processes. In this paper the method of matchedasymptotic expansions is used to solve the low tension, large-deflectionshape of a convected strip hanging under gravity. To leading order theequations for the deflection in the boundary layer and in the centre ofthe span are nonlinear. Previous treatments of this problem have notaccurately solved the leading-order equation within the boundary layer.The composite analytic solution of the full nonlinear equation isobtained and compares well with numerical solutions.     

On the Secondary Instability of Three-Dimensional Boundary Layers

One of the possible transition scenarios in three-dimensional boundary layers, the saturation of stationary crossflow vortices and their secondary instability to high-frequency disturbances, is studied using the Parabolized Stability Equations (PSE) and Floquet theory. Starting from nonlinear PSE solutions, we investigate the region where a purely stationary crossflow disturbance saturates for its secondary instability characteristics utilizing global and local eigenvalue solvers that are based on the Implicitly Restarted Arnoldi Method and a Newton–Raphson technique, respectively. Results are presented for swept Hiemenz flow and the DLR swept flat plate experiment. The main focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain experimental observations of time-dependent occurrences of an explosive growth of traveling disturbances, on the origin of high-frequency disturbances, as well as on gaining more information about threshold amplitudes of primary disturbances necessary for the growth of secondary disturbances. Received 13 July 1998 and accepted 7 July 2000     

Optimal Control of Nonmodal Disturbances in Boundary Layers

The optimal control of infinitesimal flow disturbances experiencing the largest transient gain over a fixed time span, commonly termed “optimal perturbations,” is undertaken using a variational technique in two- and three-dimensional boundary layer flows. The cost function employed includes various energy metrics which can be weighted according to their perceived importance, simplifying the task of determining which terms are essential for a “good” control scheme. In the accelerated boundary layers investigated, disturbance kinetic energy can be typically reduced by about one order of magnitude. However, it seems impossible to suppress completely over the entire control interval; “good” control strategies still permit approximately an order magnitude growth over the initial energy at some point in the interval. It is shown that the control effort efficiently targets the physical mechanisms behind transient growth. Received 5 February 2001 and accepted 15 June 2001