振动翼型和襟翼的绕流特性研究

用数值模拟手段详细地研究了振动翼型和襟翼的绕流问题,数值模拟的出发方程为Euler和N-S方程,格式为Bcam-Warming格式的改进型。数值实验主要针对流场的二大特性进行的,即振动对激波的影响和振动对分离的抑制作用,结果表明:(1)随翼型或襟翼的振动激波强度和位置也相应地变化但这一变化滞后于攻角的变化;(2)振幅加大激波强度的变化和激波运动范围也加大;(3)振动频率越高对激波的影响反而较低频时要小;(4)流动条件的不同可使升力回线的走向发生变化;(5)振动对分离有明显的抑制作用。     

三角翼大迎角不可压粘流的数值模拟

研究了人工压缩法拟压缩性系数β的选取，采用函数形式的β有效地加速了收敛过程．采用求解不可压Ｎ－Ｓ方程，对三角翼大迎角绕流进行了数值模拟，得到了与实验吻合很好的结果．分析和讨论了大迎角旋涡流动的复杂物理现象     

横向振荡圆柱绕流的格子Boltzmann方法模拟

基于格子Boltzmann方法(LBM)对不可压横向振荡圆柱绕流问题进行了数值研究. 与传统的求解宏观的N-S方程的数值方法不同, LBM求解此类问题不需要采用动网格, 而且不需要对网格进行特殊处理, 从而节约了计算成本. 结果显示, 当振荡频率增加到相应的静止圆柱绕流的自然涡脱落频率附近时, 圆柱后最新形成的集中涡距离柱体越来越近, 直到达到一个极限位置. 随后, 集中涡突然转向圆柱体另一侧脱落. 当振荡频率接近于静止圆柱的自然涡脱落频率时, 发生频率同步的现象. 随着振荡频率远离自然涡脱落频率, 同步现象消失. 在几种次谐振荡和超谐振荡下, 尾流区的涡脱落频率仍为相应的静止圆柱绕流的自然涡脱落频率.      

钝锥三维粘性绕流背风面分离的数值模拟

本文将作者在文献[1]中提出的方法推广应用于求解三维可压缩 N-S 方程和简化 N-S 方程,并对近似因式分解法应用于三维问题的稳定性进行了分析。指出,对二维问题原无条件稳定的格式,经近似因式分解后仍是无条件稳定的;对于三维问题,原无条件稳定的格式经普通近似因式分解后所得到的格式可能是不稳定的或条件稳定的。利用系数矩阵分裂法所得到的近似因式分解格式可仍是无条件稳定的,只要适当加大分裂后的系数反差即可。 文中给出了钝锥超音速三维粘性绕流结果。得到了背风面分离的流动图像,物面压力值与实验值吻合。     

翼型大攻角状态下表面吸气驻涡增升的数值模拟实验

用数值模拟方法给出了翼型大攻角状态表面吸气后绕翼型流动的某些新现象并对流场的特性进行了机理性研究,其中包括吸气对翼型背风面分离涡的驻涡增升作用;吸气孔位置对流场的影响;不同吸气强度以及间歇式吸气的增升效应。数值模拟的出发方程为N-S方程,差分格式为Beam-Warming格式。数值实验表明:(1)吸气可有效地提高翼型大攻角状态下的升力;(2)在一定吸气强度下吸气可使翼型背风面上涡的非定常脱落现象消失从而起到驻涡作用;(3)吸气孔位置在翼面的中部附近增升效果较好;(4)在一定范围内吸气强度越强其升力越高;(5)间歇式吸气也可提高平均升力,但引起升力的波动。     

超声速大攻角旋成体迎风激波数值模拟

旋成弹体在大攻角超声速流动中其迎风面上有时会出现一种迎风激波．采用非定常雷诺平均的NavierStokes方程以及两方程kω—sst和代数B—L湍流模型成功地模拟了该现象,验证了该迎风激波是由于流场主旋涡导致超声速流动偏移造成的涡诱导激波的一部分的结论．     

脱体涡流型的数值计算

本文在回顾空气动力体大迎角脱体涡流场数值计算的发展过程之后,详细讨论了各种势流计算模型。      

亚临界雷诺数下圆柱绕流的大涡模拟

本文应用Smagorinsky涡粘性模式和二阶精度的有限体积法对圆柱绕流的流场进行大涡模拟．求解了非正交曲线坐标系下的N-S方程,对雷诺数为100和20000的工况进行了计算．计算结果与实验及动力涡粘性模式的结果进行了比较,表明计算对于层流及高亚临界雷诺数的湍流流动是合理的     

流向振荡圆柱绕流的格子Boltzmann方法模拟

用一种新近发展起来的格子Boltzmann方法(LBM)在相对较小的雷诺数(Re \le 200)条件下模拟了不可压缩的流向振荡圆柱绕流问题, 考查了涡脱落模态和升阻力特性. 通过模拟, 在近尾流区发现了实验研究中已经发现的对称/反对称的涡脱落模态, 包括有些传统数值方法未发现的模态. 研究了频率锁定区域的范围及其与振幅的关系, 发现振幅越大, 发生锁定的频率区域越宽. 此外还对升阻力进行了定量意义的模拟,研究了振荡频率和振幅与升阻力的关系.      

振荡潜流带沉积层?水界面污染物输运的研究

潜流带中污染物质交换与输运特性是影响水资源环境的重要问题之一. 本文对底部为高渗透沉积层的三维槽道振荡流高Schmidt数传质问题进行了大涡模拟研究. 采用动力学亚格子模型来封闭滤波后的三维不可压缩Navier-Stokes方程以及污染物输运方程, 同时采用修正的Darcy-Brinkman-Forcheimer模型来描述沉淀有锌离子污染溶质的可渗透沉积层. 通过对沉积层内外流场和浓度场的统计特性以及瞬态结构的分析, 探究了上覆水体中振荡流驱动作用对污染物输运的动力学影响以及扩散率随振荡周期和振荡角的变化规律. 研究结果表明, 浓度通量中的湍流浓度分量在垂向物质交换中起主导作用, 流向、展向速度, 湍流强度和污染物浓度的波动跟随振荡驱动力呈现准周期变化, 同时发现沉积层?水交界面处的湍流浓度通量与法向湍流强度两者之间的变化具有明确的相关性. 并且在较大振荡角和低频振荡的情况下, 沉积层?水交界面处的有效扩散率增大, 这主要是来自于沉积层?水交界面处流体的猝发行为有效促进了湍流混合和物质交换, 进而增强了污染物的垂向输运.      

基于波尔兹曼方法研究幂律流体的椭圆柱绕流

基于格子波尔兹曼方法(lattice Boltzmann method) 和幂律流体本构方程, 建立二维流动模型, 将充分发展的速度分布与理论解进行对比, 吻合良好. 对幂律流体的圆柱绕流进行模拟, 采用了反弹格式的无滑移圆柱边界, 并使用应力积分法计算阻力系数, 分析了稳态圆柱绕流时, 阻力系数随幂律指数n 以及雷诺数Re 的变化规律. 分析了椭圆横轴/纵轴长度比和幂律指数n, 对压力系数Cp 和黏度系数Cv 的影响. 得到的变化规律与有限元方法规律一致, 验证了格子波尔兹曼模型的可行性.     

跨音速大迎角绕流的Euler方程数值模拟

采用保角变换与代数方法相结合，生成全场统一的贴体、正交Ｏ－Ｈ型网格．采用有限体积法求解Ｅｕｌｅｒ方程，模拟具有歼击机外形的全机及翼身组合体大迎角跨音速绕流．计算表明，法向力系数、气动中心位置及压力分布的计算结果与实验结果吻合良好     

大涡模拟研究展向旋转槽道湍流流动

本文用大涡模拟预测了以不同转速做展向旋转的槽道湍流流动,统计平均的流向速度型在壁面附近与已有实验数据符合很好,在通道中部的预测差异也能给出合理解释,对比不同转速的计算结果,表明展向旋转通道的湍流应力和壁面摩擦力在压力面附近提高、在吸力面附近降低,这些高阶湍流统计量的变化规律可以结合湍流应力输运方程加以解释,漩涡识别技术显示了近壁条带结构,其形态和猝发率受旋转附加力的影响发生改变,进而影响壁面摩擦速度的数值和分布,进一步考察垂直流动方向的截面内速度分布,发现旋转引起了垂直壁面方向的流动,形成正负相间排列的流向涡对,并随着转速的增加向压力面靠近。      

翼吊式双发民机机体/动力装置一体化数值分析

介绍了多块网格技术与流场分区求解方法在翼吊式双发民机机体／动力装置一体化研究中的应用。数值求解Ｅｕｌｅｒ方程模拟复杂组合体绕流。采用边界层方程／Ｅｕｌｅｒ方程耦合迭代技术进行翼面粘性修正。为保持Ｅｕｌｅｒ注解中计算网格固定,粘流／无粘流耦合迭采用表面源模型。该方法对某民用飞机模型跨音速绕流流场进行了数值模拟,机翼表面计算压力分布与实验吻合良好。     

紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用

利用紧致插值曲线(constrained interpolation profile method in Zhejiang University, CIP-ZJU) 数学模型, 对低科勒冈-卡朋特(Keulegan–Carpenter) 数KC 静止流体中振荡圆柱以及雷诺数Re = 200 时流向强迫振荡圆柱绕流进行了数值模拟. 模型在直角坐标系统下建立, 采用紧致插值曲线方法作为流场的基本求解器离散了纳维-斯托克斯方程, 基于多相流的理论实现流固耦合同步求解, 利用浸入边界方法处理固体边界. 模拟结果与现有文献结果进行比较, 二者吻合情况较好, 验证了此方法对于计算复杂流动问题的可靠性.      

自主推进俯仰震荡翼型的数值模拟研究

运用自适应多重网格法和内置边条法研究了俯仰震荡翼型的运动.通过研究翼型的原地摆动与自由游动,提出了一种确定自主推进俯仰震荡翼型推力的方法,得到了推力系数、功率系数和推进效率与Strouhal数的关系.与以往研究不同的是,我们还得到了Strouhal数随雷诺数的变化规律.此外,从自由游动翼型诱导出的涡量场中可以清楚的观察到旋涡的合并,这与实验研究非常吻合.     

??????????????????о?

应用一维水沙耦合数学模型研究了可冲刷坡面上滚波的水动力学特性. 模型的基本控制方程 采用完整的、基于守恒定律的一维浅水动力学方程,运用能够捕捉激波和泥沙运动不连续性 的WAF TVD二阶数值格式离散控制方程. 通过复演定床滚波的运动特点,对模型进行了验证. 可冲刷坡面的滚波数值模拟结果表明,床面形态对滚波水动力学特性有显著的影响.     

来流脉动对自激振荡脉冲流的影响

基于瞬变流和流体网络理论建立了由往复柱塞泵,管路和自激振荡呀嘴组成的装置系统的数学模型,分析了泵源脉动对自激振荡脉冲射流振荡特性的影响,表明当来流脉动频率小于或等于嘴嘴装置的固有频率时,能提高射流的振荡,冲蚀效果,理论分析与实验结果吻合,其结论对合理设计自激振荡射流喷嘴有一定的指导意义。     

透水床面明渠湍流时空平均特征的大涡模拟与分析

基于大涡模拟数据,研究了理想粗糙透水床面明渠湍流的时空平均特性. 考虑到空间异构性,对比分析了不同位置的时空平均流速、雷诺剪应力、构造剪应力、脉动幅度的垂线分布. 结果表明:第一,顶层床面之上,空间异构性的影响较小,不同位置的双平均流速符合类似的对数分布,但由于透水床面影响,卡门常数较不透水床面小;在床面附近,空间异构性影响较大,不同位置的双平均流速分别符合线性分布与多项式分布;在透水河床内部,靠近底层球孔的双平均流速为上部球孔双平均流速的1.55 倍. 第二,床面之上,雷诺剪应力占总剪应力的95% 以上,占有主体地位;床面附近,紊动较大,构造剪应力不能忽略,其值大约占总剪应力的15%.由于流场的各向异性,纵向与垂向的脉动幅度有所差异.