分区LES/DES混合方法及缝翼三维流场模拟

提出一种分区大涡/脱体涡模拟混合方法,采用高阶空间有限差分方法和二阶隐式LU-SGS时间推进方法,对有限展长缝翼的三维流场进行了数值计算。该混合方法吸收了脱体涡模拟方法和大涡模拟方法各自的优点,在非核心区域使用脱体涡方法进行计算,相对于大涡模拟方法具有较高的计算效率;相对于分区脱体涡模拟中采用的雷诺平均方法,在计算量没有明显增加的条件下提高了对流场各尺度流动的模拟能力。计算结果与风洞实验数据吻合良好,同时缝翼流场主要观测点的压力脉动数据为后续缝翼噪声分析和低噪声优化提供了基础。     

模型燃烧室两相燃烧大涡模拟的并行计算

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对模型燃烧室气液两相瞬态喷雾燃烧过程进行大涡模拟并行计算研究。通过大涡模拟计算结果与雷诺平均计算结果对比表明大涡模拟方法能更好的模拟流场细观结构。同时分别采用多个处理器和动态内存分配对大涡模拟并行计算程序进行数值模拟,并行计算结果表明:多个处理器之间并行计算结果与PIV测量的瞬态速度场以及出口温度分布实验数据相互吻合,表明采用并行计算和动态内存分配在保证计算结果正确的前提下,大幅度降低计算机内存和大涡模拟计算时间。     

地面效应对尾涡消散的影响研究

采用大涡模拟方法可以分析和计算近地阶段尾涡消散和运动规律,但模型繁琐、计算复杂。本研究基于镜像涡方法建立了形式相对简单、计算迅速的近地尾涡运动模型和消散模型,计算结果与激光雷达测量数据和大涡模拟数据的偏差都不超过5％,满足动态尾涡间隔研究的需要。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

舰船空气流场数值模拟方法研究

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法,兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势。论文分别采用RANS和DES方法对舰船空气流场进行数值模拟,通过所得流场结果与实测数据的对比分析,说明DES方法对大分离流动具有更强的模拟能力,更加适用于舰船空气流场预报。     

特征无反射边界条件应用于二维圆柱绕流模拟

采用特征无反射边界条件和吸收层技术相结合的边界处理方法,研究了二维圆柱绕流模拟的边界反射问题。近年来在计算气动声学中提出了将大涡模拟和声学传播方程结合的混合方法求解声场,基于计算量的考虑,须在较小的计算区域内进行大涡模拟,这就需要在计算区域较小时,边界的反射尽可能小。作为大涡模拟研究的技术储备,将该方法应用于二维圆柱绕流的模拟中。数值结果表明这种方法在计算区域较小时,能起到很好的无反射效果,计算结果与前人在计算区域取得很大时的结果比较吻合。     

用确定性涡方法模拟二维钝体绕流

以离散涡方法为基础,采用粒子强度交换法模拟涡量扩散方程,对粒子场分布产生扭曲变形时采用粒子强度重分配,数值模拟了突然起动圆柱、风攻角为的方形截面和大海带东桥的流场,计算结果与解析解及其他数值方法计算的结果均吻合的很好,数值结果显示确定性涡方法的高精度的特性使其在模拟钝体扰流问题方面具有很大的优势.     

超音速底部流动的大涡模拟研究

超音速底部流动是CFD的难点问题之一。采用基于可压缩流动的大涡模拟方法求解N-S方程,对超音速底部流动进行了模拟研究,分析了其流动结构。计算结果表明：得到的速度型、底面压力分布与实验值吻合,优于雷诺平均方法的结果。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

固体火箭内流场的时均和大涡模拟技术对比研究

分别采用时均RNG湍流模拟和大涡模拟(LES)方法对某型固体火箭发动机内流场进行了瞬态数值模拟研究,建模过程中针对推进剂燃烧以及燃面边界采用了修正壁面通量的方法,并使用了动网格技术。结果表明:大涡模拟技术可以得到发动机内湍流流场中详细拟序结构,时均湍流模拟对于火箭发动机内流场时均量的模拟仍有较好的工程精度。     

用DES数值模拟具有横向喷流的紊流流场

DES(DettachedEddySimulation)是把RANS(雷诺平均Navier Stokes方程)方法及LES(大涡模拟)方法结合起来的模拟有脱体涡的紊流流场的数值模拟方法。其主要思想是在物面附近解雷诺平均Navier Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文用DES及混合非结构网格数值模拟复杂紊流流场,算法采用Osher逆风格式,对超音速钝头体、三角翼、某型飞机及带喷流的导弹绕流用该程序进行了数值模拟。结果表明用DES模拟有分离的紊流流场是高效的。     

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

基于人工神经网络的可压缩湍流大涡模拟模型

在国家数值风洞（NNW）工程项目的指导下，空间人工神经网络（SANN）模型被用于强可压缩湍流大涡模拟（LES）研究，其中流场的湍流马赫数分别为0.6、0.8、1.0。基于湍流的多尺度空间结构特性和人工神经网络方法发展的高精度空间神经网络（SANN）模型适用于不可压缩湍流和弱可压缩湍流。对于强可压缩湍流，流场中会出现激波结构，给大涡模拟带来了挑战。本文的研究结果表明：SANN模型适用于强可压缩湍流的大涡模拟。在先验分析中，SANN模型预测的亚格子应力和亚格子热流的相关系数超过0.995，远远高于梯度模型和近似反卷积模型等传统模型；传统模型的相对误差大于30%，而SANN模型在这方面有很大的改进，相对误差低于11%。在后验分析中，与隐式大涡模拟（ILES）、动态Smagorinsky模型（DSM）、动态混合模型（DMM）相比，SANN模型能更精确地预测能谱、各类湍流统计特性以及瞬态流动结构。因此，基于湍流多尺度空间结构特性的人工神经网络模型加深了人们对强可压缩湍流亚格子建模的认识，同时可以服务于NNW工程的流体力学模型构造。     

间断Galerkin方法计算叶型转捩流动

使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.      

Large eddy simulation of turbulent diffusion flame combustion using a conserved scalar methodology

The present paper describes an LES prediction of turbulent diffusion flame combustion in a simplified axi-symmetric combustor geometry.The calculations are carried out using a well-tested finite volume incompressible LES code which has been modified to handle variable density and reacting flows.The basic mixture fraction conserved scalar method is used with the chemical state relationships described by fast chemistry.The turbulence-chemistry interaction is modelled by a sub-grid PDF method and the PDF is assumed to follow a Beta-function shape.The LES predictions have been time-averaged over 3.5 flow-through times to generate the mean radial profiles of mixture fraction,product mass fraction,temperature,axial velocity and axial rms.The agreement of the LES predictions with the experimental data is good for all the above quantities at four different axial positions with largest differences at the first measurement plane.The LES method also provides information on the unsteady nature of turbulent diffusion combustion. For turbulent reacting flows with large density ratio,it was found necessary to use a relaxation method in order to remove unphysical high-frequency fluctuations and to maintain numerical stability.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

厚覆冰导线升力突变机理及数值模拟研究

针对新月形厚覆冰导线的升力系数在风攻角15°附近存在突变的问题,分别采用基于k-ωSST湍流模型的雷诺时均法和大涡模拟(LES)的数值方法对新月形厚覆冰导线在风攻角10°~20°范围进行了模拟。通过对比两种数值方法计算得到的覆冰导线气动力系数、流场结构和表面风压,发现LES方法能够更好地捕捉新月形覆冰导线表面的小尺度涡结构,得到的覆冰导线气动力参数计算结果与风洞试验数据高度吻合;而k-ωSST湍流模型难以模拟壁面上小尺度涡,捕捉不到升力系数的突变。根据覆冰导线不同壁面区域的压力分布,发现上侧壁面处的涡结构影响整体流场,并在下侧壁面曲率、来流夹角和壁面切线方向共同作用下导致升力系数突变。LES的气动力参数模拟结果可为覆冰导线防舞提供参考。     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.