小轿车绕流场数值计算和分析

本文求解了雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流的模拟采用了标准的两方程模型,修正的两方程模型,亚格子模型的大涡模拟三种要用有限体积方法离散求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。对流项的计算采用三阶迎风格式。为了提高求解过程的稳定性,采用显隐式混合格式来获得主对角元素区域占优。     

高空台扩压器的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发,利用时间相关法,采用ＭａｃＣｏｒｍａｃｋ两步显格式,结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,数值模拟了高空台扩压器的内部流场,并计算了扩压器的性能。结果表明,采用求解Ｎ－Ｓ方程的方法计算高空台扩压器性能是一种有效的方法。     

平板绕流的抛物化Navier－Stokes方程分析

对于绕有限长平板不可压粘性层流问题,循抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＰＮＳ）方程,进行ＰＬＫ坐标摄动分析,推得与Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＮＳ）方程结果完全一致的摩擦阻力规律。该规律在１５＜Ｒｅ＜１０＿４范围内成立,证实ＰＮＳ方程适应中等Ｒｅｙｎｏｌｄｓ数流动。     

再入飞行器迎风面边界层转捩的数值预估

基于求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,首先数值预估了高超声速三维流动的边界层特性;采用已有的高超声速三维边界层转捩准则,即在转捩位置Ｒｅθ／Ｍｅ＝常数,以及通过修正由ＯＮＥＲＡ提出的代数转捩模型,建立了预估三维边界层转捩的半经验数值技术并应用于三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的求解。利用ＨＡＬＩＳ模型及相关的实验数据,首次进行了在Ｍ∞＝６时不同雷诺数及迎角的验证试验。预估结果表明：在不考虑     

旋转矩形通道中道中粘性流动的数值模拟

通过有限差分方法求解时间相关的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，对旋转矩形通道中粘性流动进行了数值研究，旋转轴方向与主流方向垂直，给出了流分发展时矩形通道中的二次流结构及主流速度分布，得到的数值解与实验吻合良好。     

Navier—Stokes方程粘性项的有限体积离散法研究

Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程粘性项的有限体积离散法研究封建湖西北工业大学应用数学系［摘要］：三维Ｎ－Ｓ方程是描述真实气体流动的非线性偏微分方程组，如何快速准确地计算其粘性项至关重要。这不仅因为粘性项的计算量非常巨大，而且粘性项对急剧变化的流动梯度的...     

求解Euler／Navier—Stokes方程的有限体积龙格库塔方法

提出了一个求解定常Ｅｕｌｅｒ／Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的计算方法,求解的主管方程是积分形式的Ｅｕｌｅｒ方程或雷诺平均的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,计算方法是有限体积多步龙格库塔方法,通过在格式中加入自适应的二阶加四阶人工粘性项以保证数值计算的稳定性,用空间变时间步长、焓阻尼来加速定常问题的收敛速度。     

用Navier—Stokes方程解提供的开孔壁边界条件求洞壁干扰修正

将一个孔及其邻近区域改为开孔壁流态的基本单元，对其求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，找出穿透流量对附近诸如压力、速度等流行参数的依赖关系。根据这种透气规律提出的开孔壁边界条件改善了高速风洞洞壁干扰修正的准确度。用椭圆型差分方程组生成了贴体座标网格，在有限体积构架下，甩隐式的ＬＵ－ＡＤＩ因子分解对角化方法和显式的多步龙格－摩塔时间推进方法求解完全的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并分别试用了Ｋ－ε     

N—S方程数值解用于SB—03模型的开孔洞壁干扰修正

将一个孔及其邻近区域做为开孔壁流态的基本单元，对其求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，找出穿透流量对附近诸如压力、速度等流动参数的依赖关系。根据这种透气规律提出的开孔壁边界条件改善了高速风洞洞壁干扰修正的准确度。用椭圆型差分方程组生成了贴体座标网格。在有限体积构架下，用隐式的ＬＵ－ＡＤＩ因子分解对角化方法、显式和隐式ＴＶＤ方法、以及显式的多步龙格－库塔时间推进方法求解完全的Ｎａｖｉｅｒ－Ｓｔｏｋｅ     

跨声速机翼流场的N—S方程计算

本文给出了三维机翼粘性流场的数值模拟过程,控制方程为时均化的三维可压缩薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ两层代数湍流模型,空间离散采用中心有限体积格式,时间采用Ｒｕｎｇｅ－Ｋｕｔｔａ多齿格式进行式积分。     

二维引射—混合器流场的数值研究与验证

本文运用一维引射特性方程和不可压缩流动Navier-Stokes方程,对二维引射—混合流场进行了数值研究。首先用一维引射特性方程估算理想的引射气流流量;然后用N-S方程对混合管内的气流混合流动进行计算,用计算得到的混合管速度分布和壁面粘性剪切力来修正一维引射特性方程中的经验常数,逐次逼近真实的引射流量比。计算结果表明,本文的数值方法可以有效地预测引射流量比和混合流场。     

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

有二次涡的激波边界层干扰流动的DPNS和NS方程计算

本文针对具有二次涡复杂分离再附现象的激波边界层干扰流动,数值地考察了扩散抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＤＰＮＳ）方程组的适用情况。壁面摩阻和压力、主涡和二次涡的涡高和涡长、分离再附位置以及流线图等特性的计算表明：ＤＰＮＳ方程组的数值结果均与ＮＳ方程组的数值结果很好相符。     

跨音速操纵面嗡鸣Hopf分叉分析及结构参数对嗡鸣特性影响的研究

采用Ｈｏｐｆ分叉分析方法,对跨音速操纵面嗡鸣问题进行了研究。用二维守恒型非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ 方程计算了舵面振动时的非定常气动力,计算了出现Ｈｏｐｆ分叉的Ｍａ 数,研究了一些参数对嗡鸣的影响。将Ｈｏｐｆ分叉的计算结果与嗡鸣的时间历程计算做了对比,两者相当吻合,与飞行试验结果相比,也比较接近。     

桥梁断面的气动导数和颤振临界风速的数值计算

选取有代表性的大海带东桥、南京二桥流线型闭口箱梁断面和荆沙桥钝头开口板梁断面,用动网格法考虑流体和结构的耦合作用,用有限单元法求解原始变量二维不可压粘性流体的Ｎ－Ｓ方程,计算其气动导数,并用半逆解法确定颤振临界风速。计算结果和实验值相当吻合。     

直升机旋翼非定常气动力及结构动力响应计算

以二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为模型方程,采用Ｊａｍｅｓｏｎ有限体积并引入变系数隐式残值光顺及固接网格计算了直升机旋翼不同展向位置处翼剖面的非定常气动力,进而采用插值的方法模拟直升机在悬停及前飞状态下,旋翼的桨叶所受的非常气动力载荷,并将之与桨叶的结构振动方程相耦合,求解其时域内的结构响应,非定常下气动力计算结果与实验结果做了对比,两者比较吻合,旋翼结构响应计算结果合理。     

有侧壁干扰的机翼半模型风洞实验的N-S方程数值模拟

应用三维可压、雷诺平均Navier-Stokes方程数值模拟了机翼半模实验风洞侧壁干扰和三维机翼半模与安装侧壁结合部流场。计算采用中心有限体积多步Runge-Kutta时间步长格式,将Baldwin-Lomax两层代数紊流模型扩展应用于三维拐角区流动。用本文方法计算了3个算例 ,并与国外的有关实验数据进行比较,计算和实验吻合良好。     

高速粘性绕流的有限体积TVD格式

本文通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,解决二维超声速及高速声速流动问题。通过采用坐标交换技术,可以方便地将有限体积ＴＶＤ格式应用到各种形状的流动问题中,各种典型算例的计算表明该方法数值计算稳定,捕捉激波分辩率高,能很好地模拟超声速和高速超声流动中激波反射及激波与边界层的干扰等各种复杂流动现象。湍流计算中采用ＢａｌｄｗｉｎＬｏｍａｘ湍流。数值结果与实验结果及理论解的对比令人满意。     

多级风扇三维粘性流场的数值模拟

 进行了多级跨音速风扇粘性流动的数值模拟研究。求解全三维N-S方程,通过高分辨率格式提高跨音速流场计算的准确性,利用LU-SGS隐式解法加快计算收敛速度。给出了一双级高负荷跨音速风扇设计和非设计工况下的详细计算结果,各工况下的计算与实验数据均吻合较好。