具势亚声速小扰动流动的格子Boltzmann模拟

给出一种新的用于模拟亚声速流小扰动下的格子Boltzmann模型, 通过使用Chapman Enskog展开和多尺度技术, 得到一系列的格子Boltzmann方程,给出了满足亚声速流小扰动线化方程的平衡态分布函数表达式. 数值结果表明, 这种方法能很好地再现经典理论结果.     

基函数法在黏性可压缩流动中的研究与应用

研究基函数法在二维和三维轴对称情形下可压缩黏性流动中的应用, 并利用此法数值地计算了二维圆柱的超音速黏性绕流和三维球头的轴对称超音速黏性绕流两个算例。取一阶三角函数为基函数, 构造出导数的中心格式和迎风格式。对于可压缩N-S 方程中的对流项, 采用通量分裂法及中心格式与迎风格式相结合的技术;对于可压缩N-S方程中的黏性项, 则采用中心格式进行处理。由此, 构造出了数值求解黏性可压缩流动的一阶三角函数类型的基函数格式。两个算例的数值计算结果表明, 基函数法不仅在处理无黏可压缩流动时, 是一种高精度、高分辨率的新型计算方法, 而且在处理黏性可压缩流动问题时也一样行之有效。     

并排方柱绕流的大涡数值模拟及显示

采用大涡模拟的方法对雷诺数为7.5×105并排双方柱绕流流场进行了数值模拟,并对所采用的模拟方法进行了实验验证.利用有限体积法和SIMPLE算法计算程式,对双方柱三维物理模型求解不可压缩的N-S方程.计算得到了方柱绕流的速度场和涡量场以及不同时刻的速度分量分布情况.对方柱后速度场、涡旋形成、脱落以及波动性等问题进行了分析,结果符合物理学规律.     

低速湍流中气动噪声预测方法

 为了研究气动噪声的产生机理和传播过程,在考虑介质黏性的影响情况下,采用分解法并结合湍流模型,对低速湍流流动中气动噪声问题在时域上进行数值计算。基于非结构化同位网格和有限体积法,把可压缩N-S方程分解成不可压N-S方程、含黏性项的声扰动方程;为了考虑质点振速和声压耦合,采用SimpleC算法来同步求解不可压N-S方程和声扰动方程。进出口远场边界采用以渐近解为基础的无反射边界条件,并采用与内部区域相对应的有限体积法、时间隐式格式对其进行求解。利用所编制的程序进行层流状态下圆柱绕流气动噪声仿真验证,并与文献结果进行对比,检验本方法的正确性;并结合湍流模型将数值解法推广到湍流状态下气动噪声数值模拟中。结果表明该方法能够很好地反映流场和声场的形态,无反射边界能很好地抑制声波在边界处的反射,适合低速流气动噪声问题模拟,为实际工程中的降噪工作提供预测信息。     

展向周期变化槽道流中全局稳定性的数值模拟

用数值模拟的方法研究了基本流在展向周期变化的槽道流的全局稳定性问题.首先用线性稳定性理论求得的扰动特征函数作为初始值,计算扰动随时间的演化.经过充分的发展后,得到了全局稳定性的特征解.给出了不同展向变化幅值下的计算结果和全局稳定性下的中性曲线,分析了全局中性时雷诺数与展向不均匀程度的关系以及特征函数的性质.     

高速列车通过带有套衬结构隧道气动效应分析

采用三维、可压缩、非定常N-S方程的数值计算方法,对8辆编组的高速列车以300 km/h速度通过带有套衬结构隧道时车体表面及隧道壁面的瞬变压力进行分析。研究结果表明:数值计算结果与动模型实验结果较吻合,2种方法得到的压力曲线变化规律一致,幅值误差在5%以内;列车通过隧道时,车体头、尾处测点压力差别较大,中部测点压力差异较小;沿列车车身方向,测点正压幅值逐渐减小,负压幅值逐渐增大;隧道壁面测点压力峰峰值在隧道进、出口附近较小,而在靠近隧道中部时较大;隧道内安装套衬对于高速铁路双线隧道气动效应影响很小,加装套衬前后,测点压力幅值差异在2%以内。因此,建议在对高速铁路隧道病害整治中,考虑使用套衬技术。     

小攻角高超音速尖锥边界层的转捩预测及e~N法的改进

研究小攻角高超音速尖锥边界层的转捩预测问题.来流马赫数为6,锥的半锥角为5°,攻角为1°.采用传统的eN方法得到的转捩位置很不令人满意.原因在于,在传统的做法中,只考虑了增长中的扰动波,而事实上对于不同的子午面,扰动在开始增长前的衰减过程是各不相同的.根据我们以前的一些工作,对eN方法进行了新的诠释,并且提出实质上的改进.不仅考虑了扰动的增长过程还考虑了扰动的衰减过程.由线性理论计算得到的扰动从初始幅值增长到1%的位置被看作是转捩位置.这样得到的新的转捩预测结果相当令人满意.此外,对于基本流的计算,研究表明边界层方程可以用于小攻角的情况,其计算量远远小于直接数值模拟.     

小攻角高超音速尖锥边界层的转捩预测及e^N法的改进

研究小攻角高超音速尖锥边界层的转捩预测问题．来流马赫数为6,锥的半锥角为5°,攻角为1°．采用传统的e^N方法得到的转捩位置很不令人满意．原因在于,在传统的做法中,只考虑了增长中的扰动波,而事实上对于不同的子午面,扰动在开始增长前的衰减过程是各不相同的．根据我们以前的一些工作,对e^N方法进行了新的诠释,并且提出实质上的改进．不仅考虑了扰动的增长过程还考虑了扰动的衰减过程．由线性理论计算得到的扰动从初始幅值增长到1％的位置被看作是转捩位置．这样得到的新的转捩预测结果相当令人满意．此外,对于基本流的计算,研究表明边界层方程可以用于小攻角的情况,其计算量远远小于直接数值模拟．     

基于康达效应的超微粉气流分级数值模拟研究

基于不可压缩粘性流体的N-S方程和κ-ε两方程模型，结合固体微粒运动方程，采用有限容积法对射流分级机中微粒的运动进行了数值模拟，得到射流分级机内的压力场、速度场和颗粒运动轨迹，计算得到的粒径分布规律与分级实验结果相吻合。     

色散和Kerr系数扰动影响下光孤子传输演化的研究

从描述色散位移光纤中Kerr系数随机扰动影响下皮秒光孤子的修正非线性薛定锷（NLS）方程出发，利用 对称分步Fourier方法，对光孤子Kerr系数和群速度色散随机扰动情况下的传输演化特性进行了数值研究。结 果表明Kerr系数的随机扰动使单孤子幅值随传输距离产生了微小波动，对两孤子传输的影响主要表现在孤子间 相互作用被加强，同时还表明，群速度色散的随机变化没有加强Kerr系数随机扰动对孤子间相互作用的影响，而 是部分地互相抵消。     

Adaptive Lattice Boltzmann Model for Compressible Flows

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｔｈｅＬａｔｔｉｃｅＢｏｌｔｚｍａｎｎ(ＬＢ)ｍｅｔｈｏｄｈａｓａｔｔｒａｃｔｅｄｍｕｃｈａｔｔｅｎｔｉｏｎａｓａｎｏｖｅｌａｌｔｅｒｎａｔｉｖｅｔｏｔｒａｄｉｔｉｏｎａｌｍｅｔｈｏｄｓｆｏｒｎｕｍｅｒｉｃａｌｌｙｓｏｌｖｉｎｇｔｈｅＮＳｅｑｕａｔｉｏｎｓ[1].ＴｈｉｓｍｅｔｈｏｄｏｒｉｇｉｎａｔｅｄｆｒｏｍａＢｏｏｌｅａｎｍｏｄｅｌｋｎｏｗｎａｓＬａｔｔｉｃｅＧａｓＡｕｔｏｍａｔａ(ＬＧＡ)[2].ＴｈｅｓｔａｎｄａｒｄＬＧＡｍｏｄｅｌｉｍｐｏｓｅｓ,ｆ…     

三维非定常不可压N-S方程一个高精度求解格式

对三维、非定常、不可压Navier-Stokes方程提出了一种新的数值计算方法.在时间离散上应用3阶精度混合显-隐相结合的分裂格式,空间离散在x及y轴向采用非等间距网格的紧致有限差分格式与z向应用Fourier谱展开相杂交的数值方法逼近.经平板边界层流的验证表明,该算法具有计算精度高、稳定性好、收敛速度快等特点.同时也研究了三维、非定常流体运动下游边界问题,提出了无反射出流边界条件,以减少在有限计算区域内人工出流边界反射引起的数值误差,保证直接数值模拟的精度和准确性.该算法的提出对于求解边界层、射流及混合层等流动中的转捩与演化问题具有重要的理论意义.     

列车进入隧道时产生涡流的数值模拟

给出了列车穿越隧道时在隧道入口形成涡流场的数值模拟过程．控制方程为三维黏性、可压缩、等熵、非定常流的N-S（Navier—Stokes）方程,空间离散采用了中心有限体积法格式,时间采用预处理二阶精度多步后差分格式进行离散,对列车与隧道之间的相对移动采用移动网格技术处理．计算结果与国外的试验结果基本一致．研究结果表明,在整个涡流的形成过程中,可以将这个过程分为连续的4个阶段,即涡的形成、发展、传播及破坏阶段．研究结果还表明,涡受到列车头部形状的影响要比受到列车速度的影响要大;喷射流速度与列车速度有关系．     

轿车外流场的数值计算和分析

认识汽车的外流场对研究汽车的空气动力特性及气动噪声具有重要的意义 本文采用k -ε二方程模型来求解N -S方程 ,用有限元法对轿车外流场进行数值模拟 ,计算结果反映了汽车外流场特性 ,较好地模拟了汽车拐角处的涡流分离流动     

基于图卷积循环网络的非定常周期性流体流动预测

基于纳维-斯托克斯（Navier-Stokes, N-S）方程的计算流体力学（computational fluid dynamics, CFD）数值模拟是研究非定常周期性流体流动问题的常用方法之一，但其计算需占用大量时间。针对N-S方程求解CFD过程费时的问题，本文提出一种基于图卷积循环网络的非定常周期性流体流动预测框架，实现流体流动状态的快速预测。文中输入历史流体流动数据，通过图卷积循环网络学习非定常周期性流体流动的物理过程，预测流场变量分布。结果表明，基于图卷积循环网络的预测模型可以准确预测流体力学规律，其在流速、涡旋、压力等流场变量预测方面均具有较好表现。相比于传统N-S方程求解方法，采用图卷积循环网络预测速度提高了5倍以上，为流场变量的预测提供了一种新方法。     

宽带噪声激励下含分数阶导数的Duffing-van del Pol振子的稳态响应

应用广义谐波平衡技术,将分数阶导数表示的回复力分解为幅值依赖的等效拟线性阻尼力和拟线性回复力.然后,应用基于广义谐和函数的随机平均法得到关于系统幅值的平均伊藤随机微分方程,系统幅值的稳态概率密度通过求解与之相应的平稳FPK方程得到.最后,用原方程的Monte Carlo模拟结果验证了近似解析解的正确性.     

三角翼的可压缩外部绕流流场数值模拟

据Spalart-Allmaras模型建立了三角翼的可压缩外部绕流流动的数学模型,并考虑了近壁面流动的处理方法。利用Fluent软件模拟某三角翼周围的二维湍流流场,用空气绕流阻力、升力和力矩系数来监测解的收敛性,得到三角翼周围的速度矢量图、壁面切应力图及边界层分离点范围,显示了边界层及尾迹区的复杂流动。结果表明,数值模拟方法能真实反映边界层及尾迹区的复杂流动,为解决边界层分离问题和机翼设计提供理论依据。     

可压多孔介质流的扩展混合元解法

讨论多孔介质中两种可压缩流体混溶驱动问题数值方法,假定介质是各向异性的,渗透率系数为张量形式。压力方程采用扩展混合元方法求解压力变量、梯度变量,以及速度变量;浓度方程采用标准有限元方法求解,这一方法对各向异性渗透率多孔介质流可以获得更可靠的数值解。构造了半离散数值格式,通过理论分析得到了压力、速度以及浓度等变量的最优L2模误差估计,对浓度变量获得了H1模最优误差估计。     

基于S-A模型的覆冰输电线绕流数值分析

基于Spalart-Allmaras湍流模型(S-A模型)封闭的雷诺平均N-S方程,采用沿均匀流线的三阶Runge-Kutta法和Galerkin法分别进行时间和空间离散,得到了湍流方程的有限元格式。采用自编Matlab程序,数值模拟了不同风攻角下覆冰输电线绕流问题,得到了覆冰输电线绕流瞬态和时均流场、气动力系数时程曲线和平均气动力系数随风攻角的变化关系,并与试验结果进行对比,验证了该算法的有效性。基于流场演化和平均压力分布,分析表明风攻角变化引起的边界层分离点转移是产生升力系数尖峰突跳现象的主要原因。