A RESIDUAL-BASED UNRESOLVED-SCALE FINITE ELEMENT MODELLING FOR IMPLICT LARGE EDDY SIMULATION 1)

对应于湍流的大尺度与小尺度流场信息, 本文在有限元的框架下, 假设Navier-Stokes方程的解的形函数由大尺度和不可解尺度形函数叠加组成, 引入对应的权函数, 将Navier-Stokes方程的有限元变分形式分解为大尺度和不可解尺度系统. 根据不可解尺度系统, 构建基于Navier-Stokes大尺度方程残差的不可解尺度模型, 将其代入Navier-Stokes方程的大尺度系统, 进而数值求解大尺度系统得到Navier-Stokes方程的大尺度解. 该方法无需像传统的大涡模拟方法那样对方程的解进行过滤, 通过对形函数进行尺度分解实现解的尺度分解. 本文使用该方法的自编程序代码开展了槽道湍流的数值模拟. 通过与有限差分大涡模拟、DNS计算结果的比较, 发现在使用较少网格情况下该方法预测的平均流向速度在近壁区与DNS数据吻合, 在黏性外层略偏高; 该方法对雷诺应力预测偏低导致从流向向垂向方向上湍动能输运略偏低. 流向速度等值面图显示该方法有效捕捉到了大尺度旋涡结构; 同时在近壁区可以观察到明显的低速条带结构.     

槽道湍流近壁区多尺度输运特性研究

利用槽道湍流直接数值模拟的数据库和离散正交子波，对近壁湍流的多尺度输运特性进 行了研究. 通过在流向和展向分别进行子波多尺度分解，得到了近壁区湍动能在流向和展向 多尺度传输的不同性质，发现流向传输以能量的反传为主，而在展向能量存在明显的正传， 并且当过滤尺度较大时以正传为主. 近壁湍流能量传输的各向异性为进一步构造各向异 性大涡模拟亚格子模式提供了必要的参考.     

壁面加热湍流相干结构的子波分析实验研究

通过在平板湍流边界层沿流向固壁表面平行放置若干条通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的大尺度流向涡结构,改变了平板湍流边界层中不同尺度结构及其能量分布。采用对壁湍流多尺度结构的子波分析表明,在湍流边界层近壁区域产生规则的流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,抑制了壁湍流各种尺度湍涡结构脉动,特别抑制了能量最大尺度湍涡结构的脉动,减小由于湍流脉动引起的在湍流边界层法向和展向的动量和能量损耗,从而减小了湍流的阻力。     

旋转圆管湍流的大涡模拟数值研究

采用大涡模拟（LES）方法,并结合动力学亚格子尺度应力（SGS）模型,通过数值求解柱坐标系下的滤波Navier-Stokes方程,研究了绕管轴旋转圆管内的湍流流动特性.为验证计算的可靠性,以及动力学SGS模型对于旋转湍流的适用性,将大涡模拟计算所得的结果,与相应的直接模拟（DNS）结果和实验数据进行了对比验证,吻合良好.进一步对旋转圆管湍流的物理机理进行了探讨,研究了湍流特性随旋转速率的变化规律.当旋转速率增加时,湍流流动有层流化的发展趋势.基于湍动能变化的关系,分析了旋转效应对湍流脉动生成的抑制作用.     

可解尺度各向同性湍流的标度律

针对大涡模拟, 首先利用 EDQNM 能谱和传输谱理论, 在3种不同的过滤器下分别计算可解尺度各向同性湍流的二阶、三阶结构函数标度律, 并通过速度差扭率验证了可解尺度湍流的 ESS 理论. 研究了该可解尺度湍流标度律及速度差扭率与多个因素的关系, 这些因素包括: 两点距离与过滤尺度的比、大涡模拟雷诺数和过滤器类型. 结果显示, 当两点距离位于过滤尺度量级时或大涡模拟雷诺数较小时, 可解尺度标度律与未过滤流场的标度律相差很大进而必须加以修正, 而且可解尺度流场也不再总是满足 ESS 理论进而对应的速度差扭率也需要修正. 然后通过3个例子介绍了这些结果在大涡模拟亚格子模型中的应用.      

槽道湍流展向振荡电磁力减阻的机理研究

对槽道湍流的展向振荡电磁力控制进行了实验和数值研究. 实验通过PIV系统和浮动床阻力测试系统记录近壁区的条带变化和壁面阻力变化. 计算时, 利用谱方法直接模拟电磁力控制下的近壁流场. 实验和计算结果定性一致, 皆表明展向振荡电磁力可以减少壁面阻力, 并使条带倾斜. 计算结果还进一步揭示了电磁力减阻的机理. 电磁力诱导产生的流向涡与壁湍流的相互作用, 在近壁处形成负的脉动展向涡, 该涡将导致流向涡的倾斜和振荡, 从而抑制湍流, 减少壁面阻力.      

基于人工神经网络的湍流大涡模拟方法

大涡模拟方法(LES)是研究复杂湍流问题的重要工具,在航空航天、湍流燃烧、气动声学、大气边界层等众多工程领域中具有广泛的应用前景.大涡模拟方法采用粗网格计算大尺度上的湍流结构,并用亚格子(SGS)模型近似表达滤波尺度以下的流动结构对大尺度流场的作用.传统的亚格子模型由于只利用了单点流场信息和简单的函数关系,在先验验证中...     

用基于M-SST模型的DES数值模拟喷流流场

脱体涡数值模拟方法(dettached eddy simulation,DES)是把雷诺平均Navier-Stokes方程(RANS)方法及大涡模拟方法(LES)结合起来模拟有脱体涡的湍流流场的数值模拟方法,其主要思想是在物面附近解雷诺平均Navier-Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文在剪切应力传输(SST)湍流模型的基础上用DES及混合非结构网格数值模拟具有横向喷流的湍流流场,算法采用Osher逆风格式,利用该套程序(包括网格生成及算法),对导弹在不同马赫数下的喷流流场进行了数值模拟,并与同时开展的实验研究的结果进行了对比,结果表明用该方法处理这类问题是较准确的。     

金属丝壁面加热湍流标度律的实验研究

通过在固壁表面的平板湍流边界层沿流向平行放置若干通电加热的金属细丝,在平板表面形成沿展向周期性分布的温度场,利用该温度场引起的空气热对流,在湍流边界层近壁区域产生一组沿湍流边界层展向周期分布的流向涡结构。对壁湍流小尺度结构标度律统计特性的研究表明,金属丝加热后形成的规则流向涡结构将壁湍流各种尺度湍涡结构不规则的脉动有序地组织起来,增强了湍流小尺度结构的层次结构相似性,减小了壁湍流中小尺度结构的间歇性和奇异性,抑制了壁湍流中奇异的湍涡结构。     

基于卡门尺度和滤波方法的SST方程改进

基于滤波方法和卡门尺度对原始剪切应力输运(shear stress transport, SST)湍流模型进行了改进,提出了一种卡门尺度修正的滤波SST 方法. 湍流多尺度效应必须在分离流场模拟中给予反映,该方法减弱了雷诺平均(Reynolds averaged Navier-Stokes, RANS)方法时间平均特性对于流场脉动量的压迫作用,在流场中引入了大涡模拟(large eddy simulation, LES)方法的亚格子模型,形成一种新型的脱体涡模拟方法(detached eddy simulation,DES)方法;同时,为了降低原始DES方法在网格加密过程中产生网格诱发的雷诺应力损耗,利用卡门尺度对滤波因子进行修正. 平板边界层算例中,卡门尺度对于RANS方法的跟随性远远强于DES方法,在边界层内的速度型和RANS方法吻合很好,而DES方法在加密过程中速度型的鲁棒性较差,说明卡门尺度在有效地保护了边界层内使用RANS求解,降低速度型偏离对数率现象的产生;HGR-01翼型算例证明BY-SST方法可以有效的避免网格诱导分离现象的产生;证明BY-SST方法在分离流动中的精度高于DES方法.     

A Gas-Kinetic Scheme for Turbulent Flow

RANS simulations may not provide accurate results for all flow conditions. The interaction between a shock wave and a turbulent boundary layer is an example which may still be difficult to simulate accurately. Beside the inability to reproduce physical phenomena such as shock unsteadiness, the argument is put forward that the conventional numerical schemes, based on the Navier-Stokes equations, may be unable to generate a physically consistent turbulent stress tensor in the presence of large unresolved scales of motion. A large ratio between unresolved and resolved scales of motion, a sort of Knudsen number based on turbulent fluctuations, might introduce inaccuracies for which the turbulence model is not accountable. In order to improve the accuracy of RANS simulations, researchers have suggested various ad-hoc modifications to standard turbulence models which limit eddy viscosity or the turbulent stress tensor in the presence of strong gradients. Gas-kinetic schemes might be able to improve RANS predictions in shocklayers by removing or limiting the errors caused by the large scales ratio. These schemes are a class of their own; in the framework of a finite-volume or finite-elements discretizations, they model the numerical fluxes on the basis of the Boltzmann equation instead of the Navier-Stokes equations as is conventionally done. In practical terms, these schemes provide a higher accuracy and, more importantly, an in-built “multiscalar” mechanism, i.e. the ability to adjust to the size of unresolved scales of motion. This property makes them suitable for shock-capturing and rarefied flow. Gas-kinetic scheme may be coupled to a conventional RANS turbulence model; it is shown that the turbulent stress tensor is naturally adjusted as a function of the unresolved-to-resolved scales ratios and achieves a higher physical consistency than conventional schemes. The simulations shown - well-known benchmark cases with strong shock-boundary layer interactions - have been obtained with a standard two-equation turbulence model (k- ω). It is shown that the gas-kinetic scheme provides good quality predictions, where conventional schemes with the same turbulence model are known to fail.     

Interaction between longitudinal vortices and flat plate boundary layer

The interaction between longitudinal vortices and flat plate boundary layer has been studied numerically for both laminar and turbulent flow situations. The vortices are assumed to be placed in an otherwise two-dimensional boundary layer flow. The flow is assumed to be incompressible and steady. Considering the fact that the velocity, vorticity and temperature gradients in the transverse directions are much larger than the longitudinal (streamwise) gradients for these flows, the original Navier Stokes equations are parabolized in the streamwise direction. A simple model, based on Boussinesq hypothesis, is used for turbulent flow. The discretized equations are then solved step by step in the streamwise direction, using an iterative procedure at each station. Numerical solutions have been obtained for different parameters, such as the Reynolds number, the circulation and the initial position of the vortices. The computed flow patterns and the skin friction coefficient and Stanton number are found to be qualitatively consistent with available experimental results. It is shown that the interaction between the vortices and the boundary layer may severely disturb the boundary layer flow field and thus considerably increase the local skin friction and heat transfer rate on surface of an aircraft.     

Projection‐based variational multiscale method for incompressible Navier–Stokes equations in time‐dependent domains

A variational multiscale method for computations of incompressible Navier–Stokes equations in time‐dependent domains is presented. The proposed scheme is a three‐scale variational multiscale method with a projection‐based scale separation that uses an additional tensor valued space for the large scales. The resolved large and small scales are computed in a coupled way with the effects of unresolved scales confined to the resolved small scales. In particular, the Smagorinsky eddy viscosity model is used to model the effects of unresolved scales. The deforming domain is handled by the arbitrary Lagrangian–Eulerian approach and by using an elastic mesh update technique with a mesh‐dependent stiffness. Further, the choice of orthogonal finite element basis function for the resolved large scale leads to a computationally efficient scheme. Simulations of flow around a static beam attached to a square base, around an oscillating beam and around a plunging aerofoil are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

A quasi-similarity analysis of the turbulent boundary layer on a flat plate

The partial differential equation of the boundary layer on a flat plate are simplified by using the universal variables for turbulent flow. For laminar flow this gives boundary layer having a finite thickness and a friction coefficient differing by a few percent from the Blasius value. For a turbulent flow a differential equation for the velocity distribution is obtained with a parameter which varies slowly with the streamwise coordinate. The numerical value of this parameter is determined as an eigenvalue of the differential equations giving a velocity profile which evolves as the boundary layer thickens. Numerical calculations using a simple eddy viscosity model gave results in very good agreement with experiment.     

Taylor and Lagrange correlations in a turbulent free shear layer

 The objective of this research was to study the effect of various Lagrange-tracking correlation methods in estimating the eddy lifetime for a two-stream, turbulent, planar free shear layer. Zeroth-, first- and second-order Lagrange correlation methods were applied to the time-evolving velocity field data collected from a cinematic particle image velocimetry technique. A time scale associated with the eddy lifetime was obtained based on a 2/e correlation of either vorticity or streamwise velocity fluctuations. When based on vorticity, this time scale significantly increased as expected when the tracking was computed with a second-order Lagrangian tracking technique as compared to a (zeroth-order) Taylor hypothesis approach. However when based on streamwise velocity fluctuations, this time scale did not increase significantly for the higher order projection methods. The latter result is attributed to occurrences of “reverse correlation” of the instantaneous streamwise velocity fluctuations caused by eddy rotation. Received: 2 April 1997/Accepted: 3 September 1997     

Proper orthogonal decomposition analysis of coherent motions in a turbulent annular jet

A three-dimensional incompressible annular jet is simulated by the large eddy simulation(LES) method at a Reynolds number Re = 8 500. The time-averaged velocity field shows an asymmetric wake behind the central bluff-body although the flow geometry is symmetric. The proper orthogonal decomposition(POD) analysis of the velocity fluctuation vectors is conducted to study the flow dynamics of the wake flow.The distribution of turbulent kinetic energy across the three-dimensional POD modes shows that the first four eigenmodes each capture more than 1% of the turbulent kinetic energy, and hence their impact on the wake dynamics is studied. The results demonstrate that the asymmetric mean flow in the near-field of the annular jet is related to the first two POD modes which correspond to a radial shift of the stagnation point. The modes 3 and 4 involve the stretching or squeezing effects of the recirculation region in the radial direction. In addition, the spatial structure of these four POD eigenmodes also shows the counter-rotating vortices in the streamwise direction downstream of the flow reversal region.     

粘流-无粘干扰流动(IF)理论

对不可压缩层流二维干扰流动,本文提出一个干扰流动(IF)理论。IF理论要点为:1)干扰流动沿主流的法向被分为三层即粘性层、干扰层和无粘层,引进了法向动量交换为主导过程的干扰层概念。2)利用力学守恒律、三层匹配关系及文中引进的干扰模型,把三层的空间尺度及惯性-粘性诸力的数置级表示为单参数m的函数,m<1/2·3)导出描述各层流动的控制方程、导出描述全城流动的控制方程为简化Navie-Stokes(SNS)方程。IF理论适用于不存在分离的附着干扰流动以及存在分离的大范围干扰流动,经典边界层(CBL)理论和流动分离局部区域Triple-Deck(TD)理论分别是本文理论在参数m=O和1/4时的两个特例,本文理论容易推广到可压缩、三维及湍流流动。     

湍流边界层等动量区演化机理的实验研究

等动量区是瞬时流场中流体动量接近的局部区域,其生成和分布与相干结构密切相关. 对等动量区的研究有助于更深入认识湍流边界层相干结构,但目前对其演化过程还缺乏实验支持和机理分析. 设计并使用移动式高时间分辨率粒子图像测速技术(TRPIV)系统对光滑平板湍流边界层进行了跟踪测量,用滤波方式对数据进行降噪,结合对直接数值模拟数据的插值结果,获得脉动速度信号. 使用改进方法去掉非湍流的影响,检测边界层内的等动量区,得到其数量的时间序列,结合流向速度概率密度函数分布的变化,分析得出了等动量区的数量在大的时间尺度下从一个稳态到另一个稳态的阶梯状变化特点. 分解不同尺度的脉动速度,对大尺度和小尺度脉动信号进行条件平均,发现大尺度脉动对等动量区数量变化起主要作用,表现为不同速度流体通过发生不同猝发事件改变流向速度概率密度函数分布. 分析流向大尺度脉动空间分布的变化,发现等动量区内常含有多个大尺度脉动区域,不同区域的扩张、收缩、分裂、合并影响流向速度的集中程度,进而导致等动量区数量的变化.      

The response of a turbulent boundary layer to artificial disturbances

With periodic fluid injection through small slots, a turbulent boundary layer is artificially disturbed on scales that are of the order of those of the natural quasi-periodic events. The periodic phase-average of the streamwise fluid velocity is determined from hot-film measurements, and used to find the coherent velocity component as defined by the triple decomposition. It appears that, when a disturbance is active, the generated flow pattern is very similar to the one caused by the interaction of a crossflow and a jet. However, when it is terminated, the turbulent boundary layer returns to its undisturbed state.