流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

聚合物槽道湍流的直接数值模拟

应用直接数值模拟方法研究了聚合物对槽道湍流的减阻特性.数值计算方法采用谱方法,时间积分采用二阶精度的时间分裂格式.应用FENE-P模型对聚合物槽道湍流进行定压降的直接数值模拟计算,揭示了聚合物湍流减阻的启动(onset)现象;应用O-B模型进行定流量的直接数值模拟计算,分析了其流动特征.计算结果表明:在低减阻形式中聚合物对湍流的影响主要发生在壁面区,而高减阻形式中聚合物对湍流的影响则遍及全流场.     

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.      

槽道内涡波流场展向涡的分布特征

通过2DPIV测量了槽道内涡波流场的瞬态速度矢量场,利用雷诺分解和旋转强度的λci准则进行了展向涡的识别;并从识别的逆时针和顺时针展向涡的尺度属性、力学属性和运动属性方面分析了展向涡的分布特征。结果表明,流场中的顺时针展向涡的数量大于逆时针展向涡的数量;在法向位置0.05~0.45范围内,逆时针展向涡的旋转强度呈抛物线形状变化,顺时针展向涡整体呈下降趋势;逆时针展向涡的平均直径和椭圆平均长轴长度随着法向位置的增加而缓慢增加,而顺时针展向涡趋于下降;逆时针展向涡的椭圆长轴倾角最大值为67.28°,顺时针展向涡的等效椭圆的离心率呈下降趋势,似圆性较好;逆时针展向涡对全局平均的〈vωz〉分量在法向位置0.15处出现最大值61%;展向涡对流场雷诺应力的影响较小;逆时针展向涡的数量密度在0~0.08法向位置范围内呈急剧上升趋势,并达到最大值0.36,顺时针展向涡数量密度整体呈倒U形分布;不同旋向的展向涡占总展向涡数量密度的比值都大于0.3;展向涡的流向对流速度整体上小于流场的平均速度,且在槽道中心区域差值较大;展向涡的法向对流速度随着法向距离增加,从负值上升到正值,在中心主流两侧分别发生喷射和扫掠2种现象。     

微重力量值对气液两相流相分布及液相湍流统计量影响的数值研究

气液两相流在空间领域具有广阔的应用前景, 深入理解微重力量值(微重力大小)对相分布和液相湍流的影响十分必要. 采用欧拉elax-elax拉格朗日双向耦合模型深入研究了不同微重力量值对相分布和液相湍流的影响. 液相速度场通过直接数值模拟求解, 气泡的运动轨迹由牛顿运动方程跟踪. 研究表明, 气泡分布和液相湍流与微重力量值均具有直接联系. 在低微重力量值下, 气泡近似均匀分布在槽道内, 且对液相湍流统计量几乎没有影响; 然而当微重力量值较高时, 大量气泡聚集在壁面附近, 液相湍流由于气泡的注入受到极大调制.      

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

基于非湍流脉动动能方程的高超声速转捩预测

通过求解基于稳定性理论的非湍流脉动动能输运方程,预测了来流Ma为5.91的裙锥绕流流动的层流-湍流转捩点,并采用SST(Shear-Stress Transport)湍流模型结合代数间歇因子模型对转捩流场进行了数值模拟.结果表明非湍流脉动动能输运方程模型可以很好地捕捉逆压梯度对转捩前期边界层内的不稳定波动频率的影响.等温壁面条件下计算得到的转捩位置与静风洞实验结果基本一致.壁面冷却促进第2模态不稳定波动的增长,使转捩提前发生.代数间歇因子模型模拟高马赫数流动的转捩区长度较短,温度峰值偏低.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

近壁区Reynolds应力输运方程的理论模型

在湍流边界层的近壁区,采用对称共振三波的理论模型描述相干结构,根据理论模型对Reynolds应力输运方程逐项进行计算和分析,结果与直接数值模拟符合很好.这不仅在理论上有益于对湍流物理机制的了解,而且展现了一种可能性,即根据相干结构的理论知识来改造湍流模型,使之具有更清晰的物理内涵以提高近壁区的预估精度.     

定常压力梯度边界层相干结构的直接数值模拟

根据流动稳定性理论,将不稳定波的一个周期作为相干结构的初值,采用直接数值模拟方法对有压力梯度湍流边界层中相干结构的演化进行了研究,得出其各种特性的变化与实验观测到的结果一致.     

逆压梯度湍流边界层中的相干结构

利用流动稳定性理论中一般共振三波的概念对4种定常逆压梯度以及两种变逆压梯度湍流边界层的相干结构进行研究.其中2种变逆压梯度分别为压力梯度由小变大和由大变小,变化的终了具有相同的压力梯度.通过计算得出逆压梯度能够激励相干结构的生成和发展,而且逆压梯度越大作用越强.在变压力梯度流动中,压力梯度由大变小的情况对相干结构的激励最强.此外,逆压梯度使相干结构的流向尺度变小,而展向间距基本在100个粘性长度左右,这些变化趋势与实验以及直接数值模拟结果一致.     

近壁区理论耗散率模型的改进

提出了在边界层的近壁区,采用共振三波的理论模型描述湍流边界层相干结构,由此解决了对称三维波理论模型不能处理三阶相关量的问题,根据新的理论模型对ε方程逐项进行计算和分析.与原理论模型计算结果相比,本文理论结果与直接数值模拟(DNS)符合更好.     

旋转方位角对斜肋通道内部流动与换热的影响

数值模拟了静止及旋转状态下径向出流60°斜肋通道的内流动与换热分布,研究了3种旋转方位角:0°,-45°及+45°.分析了哥氏力与肋片对旋转通道内流及换热的耦合作用机理.计算结果表明:旋转方位角对通道的内流与换热影响显著.在0°旋转方位角下,通道前后缘的换热差异最大;在-45°方位角下,通道整体换热性能最差;在+45°方位角下,通道具有最佳的整体换热性能.