2004年 22卷 第4期
2004, 22(4): 377-383,.
摘要：
时间分裂算法很难获得带源项的守恒系统的定常解,本文中,考虑使用基本的隐式方法积分线性、二次和三次源项,研究了时间分裂算法不能收敛到数值定常态的原因.采用一个刚性参数的方式预测定常态的误差,刚性参数依赖于特定的源项.当源项是解的线性函数时,不存在定常态误差;当源项是解的非线形函数时,发现定常状态的误差是刚性参数的单调增加函数.定常态误差的分析将推广到标准k-epsilon双方程湍流模型计算的情况.
时间分裂算法很难获得带源项的守恒系统的定常解,本文中,考虑使用基本的隐式方法积分线性、二次和三次源项,研究了时间分裂算法不能收敛到数值定常态的原因.采用一个刚性参数的方式预测定常态的误差,刚性参数依赖于特定的源项.当源项是解的线性函数时,不存在定常态误差;当源项是解的非线形函数时,发现定常状态的误差是刚性参数的单调增加函数.定常态误差的分析将推广到标准k-epsilon双方程湍流模型计算的情况.
2004, 22(4): 384-388.
摘要：
在1m口径圆形截面低速风洞中用一组大小不同外形相似的前缘后掠角为70°的三角翼模型,测量了模型作俯仰运动时的上表面的动态压力分布,研究了洞壁对模型表面动态压力分布的影响.在3m×2.5m矩形截面低速风洞中,用另一组外形相似的前缘后掠角为70°的三角翼模型进行非定常测力试验,并测定了洞壁上有限测压点的动态压力,研究了模型大小对非定常壁压的影响和洞壁对模型非定常法向力系数的影响,并用混合法对大模型的非定常法向力进行了洞壁干扰修正,大模型的法向力系数修正结果与小模型的法向力系数基本一致.
在1m口径圆形截面低速风洞中用一组大小不同外形相似的前缘后掠角为70°的三角翼模型,测量了模型作俯仰运动时的上表面的动态压力分布,研究了洞壁对模型表面动态压力分布的影响.在3m×2.5m矩形截面低速风洞中,用另一组外形相似的前缘后掠角为70°的三角翼模型进行非定常测力试验,并测定了洞壁上有限测压点的动态压力,研究了模型大小对非定常壁压的影响和洞壁对模型非定常法向力系数的影响,并用混合法对大模型的非定常法向力进行了洞壁干扰修正,大模型的法向力系数修正结果与小模型的法向力系数基本一致.
2004, 22(4): 389-394,.
摘要：
本文运用工程方法设计了喷气Z箍缩装置的超声速喷管,通过数值求解Navier-Stokes方程模拟了气流从储气室到喷管出口的非定常流动过程,最后用DSMC方法计算了喷管出口气体向真空中的膨胀,得到了喷气负载的气流位形和质量线密度.计算结果表明,除了储气室的压力、喷管的喉道和出口尺寸等因素外,负载的特性参数较强地依赖于阀门开启和关闭时间、阀门开启的速度、开启的口径大小以及阀门后真空管道的几何形状.采用定常流动模型比非定常流动模型预测的负载质量线密度大10倍左右,而后者与实验结果较为一致.
本文运用工程方法设计了喷气Z箍缩装置的超声速喷管,通过数值求解Navier-Stokes方程模拟了气流从储气室到喷管出口的非定常流动过程,最后用DSMC方法计算了喷管出口气体向真空中的膨胀,得到了喷气负载的气流位形和质量线密度.计算结果表明,除了储气室的压力、喷管的喉道和出口尺寸等因素外,负载的特性参数较强地依赖于阀门开启和关闭时间、阀门开启的速度、开启的口径大小以及阀门后真空管道的几何形状.采用定常流动模型比非定常流动模型预测的负载质量线密度大10倍左右,而后者与实验结果较为一致.
2004, 22(4): 395-398.
摘要：
针对新近提出的双凹槽和预燃室结构的燃烧室概念,采用混合通量分离法及NND格式,对其中的冷态流场进行了数值模拟.结果表明:在没有喷流条件下,该结构是典型的开式流动,有利于混合,在工程应用范围内,可以用二维结果粗略地估算三维情况;但在有喷情况下,开式结构被打破,必须对三维流动作真实模拟.
针对新近提出的双凹槽和预燃室结构的燃烧室概念,采用混合通量分离法及NND格式,对其中的冷态流场进行了数值模拟.结果表明:在没有喷流条件下,该结构是典型的开式流动,有利于混合,在工程应用范围内,可以用二维结果粗略地估算三维情况;但在有喷情况下,开式结构被打破,必须对三维流动作真实模拟.
2004, 22(4): 399-403.
摘要：
侧向喷流是导弹在中高空进行姿态控制的重要手段.本文用数值方法求解N-S方程,模拟了二种带侧向喷流控制的喷流与绕流的干扰流场.通过定性流场结构和物面压力分布与实验结果的比较,数值模拟能较好地模拟侧向喷流与外流干扰产生的复杂流场特性.给出了合理的流场结构和干扰气动力特性,表明数值方法可以作为型号侧向喷流控制设计的研究手段.
侧向喷流是导弹在中高空进行姿态控制的重要手段.本文用数值方法求解N-S方程,模拟了二种带侧向喷流控制的喷流与绕流的干扰流场.通过定性流场结构和物面压力分布与实验结果的比较,数值模拟能较好地模拟侧向喷流与外流干扰产生的复杂流场特性.给出了合理的流场结构和干扰气动力特性,表明数值方法可以作为型号侧向喷流控制设计的研究手段.
2004, 22(4): 404-409.
摘要：
为了提高AUSM+格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM+格式的基础上经过数值研究分析,对AUSM+格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM+格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM+格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM+格式继承了原AUSM+格式优点,计算处理问题能力增强,可以
为了提高AUSM+格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM+格式的基础上经过数值研究分析,对AUSM+格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM+格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM+格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM+格式继承了原AUSM+格式优点,计算处理问题能力增强,可以
2004, 22(4): 416-421.
摘要：
为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.
为了利用被动控制增强可压缩混合层中混合效率,本文用数值模拟的方法研究了三维可压缩混合层中扰动向下游的演化.通过在入口处引入不同的不稳定T-S波,给出了扰动随空间的演化过程,及马赫数对扰动演化的影响.计算结果表明在马赫数较小时,扰动波增长很快,当马赫数增大时,扰动所形成的流向涡的强度会减弱.这一结论与流动稳定性分析结果是一致的.
2004, 22(4): 422-426.
摘要：
采用线性稳定性分析的方法对粘性气体中的粘性液体射流的分裂与雾化机理进行了分析,数值计算表明:液体射流分裂与雾化过程中存在一临界气体韦伯数We2c=1,We2＜We2c对应的是射流分裂过程,We2＞We2c对应的是射流雾化过程,射流分裂过程和雾化过程的机理有所不同.当We2＜1时,We2对射流分裂过程具有稳定性的作用;当We2＞1时,We2对射流雾化过程起着不稳定性的作用.液体Reynolds数Re1在整个射流过程中始终起着不稳定性的作用,气体Reynolds数Re2的作用却相反.气液密度比Q,即气动力在整
采用线性稳定性分析的方法对粘性气体中的粘性液体射流的分裂与雾化机理进行了分析,数值计算表明:液体射流分裂与雾化过程中存在一临界气体韦伯数We2c=1,We2＜We2c对应的是射流分裂过程,We2＞We2c对应的是射流雾化过程,射流分裂过程和雾化过程的机理有所不同.当We2＜1时,We2对射流分裂过程具有稳定性的作用;当We2＞1时,We2对射流雾化过程起着不稳定性的作用.液体Reynolds数Re1在整个射流过程中始终起着不稳定性的作用,气体Reynolds数Re2的作用却相反.气液密度比Q,即气动力在整
2004, 22(4): 427-431,.
摘要：
提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明
提出并介绍了基于直角切割网格的分区面搭接技术,采用变长宽比网格方法,成功地进行了存在外形间断(剪刀叉)的民机增升装置的网格生成和绕流Euler方程数值模拟.根据外形的特点,运用"根"网格分区算法,降低了整个网格的生成难度;通过基于外形的自适应网格加密技术,详细地描述了外形上的大量缝道和凹槽,提高了网格质量;在分区交界面上,利用重叠面积切割和面搭接算法实现了两侧网格间的流场信息传递,并保证了通量守恒;采用中心有限体积方法,结合双时间推进算法,完成流场的Euler方程数值模拟,计算结果与实验数据吻合良好,说明
2004, 22(4): 432-437.
摘要：
基于混合通量分裂的思想,通过应用Gauss-Seidel迭代求解差分方程,构造了一种时空二阶精度、无条件稳定的隐式迭代NND算法,并讨论了时间精度与稳定性,亚迭代收敛判则,几何守恒律的应用以及动壁边界条件等相关问题.通过引入一种简便易行的加权函数来综合刚性动网格生成技术和超限插值生成动网格这两种方法的优点,发展了一种计算量小、比较实用的加权动网格生成技术.作为应用实例,本文给出了多个动态物体绕流的数值模拟算例,计算结果表明了本文数值方法的成功.
基于混合通量分裂的思想,通过应用Gauss-Seidel迭代求解差分方程,构造了一种时空二阶精度、无条件稳定的隐式迭代NND算法,并讨论了时间精度与稳定性,亚迭代收敛判则,几何守恒律的应用以及动壁边界条件等相关问题.通过引入一种简便易行的加权函数来综合刚性动网格生成技术和超限插值生成动网格这两种方法的优点,发展了一种计算量小、比较实用的加权动网格生成技术.作为应用实例,本文给出了多个动态物体绕流的数值模拟算例,计算结果表明了本文数值方法的成功.
2004, 22(4): 438-442.
摘要：
对大坝渗流反问题进行了初步研究,基于渗流场的有限元计算建立了反演坝体渗流参数的计算方法.该方法将参数反演问题转化为参数优化问题,利用大坝运行水位变化过程中坝体内某些测点上的水头测量数据,用牛顿-拉夫逊算法优化出坝体的渗流参数.通过算例计算分析表明,该反演方法具有推导简洁、反演精度高、适应性好的优点,有较好的工程实用价值.
对大坝渗流反问题进行了初步研究,基于渗流场的有限元计算建立了反演坝体渗流参数的计算方法.该方法将参数反演问题转化为参数优化问题,利用大坝运行水位变化过程中坝体内某些测点上的水头测量数据,用牛顿-拉夫逊算法优化出坝体的渗流参数.通过算例计算分析表明,该反演方法具有推导简洁、反演精度高、适应性好的优点,有较好的工程实用价值.
2004, 22(4): 443-446,.
摘要：
反问题作为一种自由(未知)边界问题近年来得到了广泛的研究.变域变分通过把自由边界结合在变分泛函中,使其与求解流场的控制方程结合起来,从而自由边界求解和流场分析可以完全耦合进行.本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型,计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果.文中还采用了两种计算网格(H型和O型)来设计翼型,比较了它们的精度,计算表明本方法完全能用于翼型的反设计.
反问题作为一种自由(未知)边界问题近年来得到了广泛的研究.变域变分通过把自由边界结合在变分泛函中,使其与求解流场的控制方程结合起来,从而自由边界求解和流场分析可以完全耦合进行.本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型,计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果.文中还采用了两种计算网格(H型和O型)来设计翼型,比较了它们的精度,计算表明本方法完全能用于翼型的反设计.
2004, 22(4): 447-451.
摘要：
本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.
本文应用直接仿真蒙特卡罗(DSMC)方法,研究了Crookes热辐射计内由于叶片两表面的温度不同而引起的稀薄气体的流动,给出了热辐射计叶片高温面温度900K,低温表面温度为300 K,四种不同克努森数(Kn)下流速矢量场、温度场和叶片两面压力差的仿真结果.并与热辐射计在不同Kn数下的转速实验测量结果进行了比较,两者结果基本一致.
2004, 22(4): 452-456.
摘要：
本文介绍国内外男女用的几种品牌标枪的气动力特性、几何和物理参数测量.测量结果表明,同类标枪在总长、质量、重心、面心方面差异不大,一般在1%以内,但转动惯量相差较大,约为3～5%.标枪的气动阻力、升力随迎角增加而增加.标枪的最大升阻比在迎角20°左右,其值约为2.标枪重心前移后,其压力中心始终位在重心之后,因此,标枪气动力性能产生很大变化.
本文介绍国内外男女用的几种品牌标枪的气动力特性、几何和物理参数测量.测量结果表明,同类标枪在总长、质量、重心、面心方面差异不大,一般在1%以内,但转动惯量相差较大,约为3～5%.标枪的气动阻力、升力随迎角增加而增加.标枪的最大升阻比在迎角20°左右,其值约为2.标枪重心前移后,其压力中心始终位在重心之后,因此,标枪气动力性能产生很大变化.
2004, 22(4): 457-460.
摘要：
抖振是大跨度桥梁易发生的一种风振形式.基于大型通用有限元软件ANSYS,本文提出了一种大跨度桥梁抖振时域实用分析方法.在该方法中,自激力是通过向有限元模型中添加气动刚度和气动阻尼单元的方式实现的.以深圳湾大桥为例,计算了该桥在不同风速和迎角下的抖振位移响应,与传统频域分析结果的比较表明了该方法的可行性.
抖振是大跨度桥梁易发生的一种风振形式.基于大型通用有限元软件ANSYS,本文提出了一种大跨度桥梁抖振时域实用分析方法.在该方法中,自激力是通过向有限元模型中添加气动刚度和气动阻尼单元的方式实现的.以深圳湾大桥为例,计算了该桥在不同风速和迎角下的抖振位移响应,与传统频域分析结果的比较表明了该方法的可行性.
2004, 22(4): 461-465.
摘要：
在超声速边界层的计算中,出口处边界条件的好坏会影响计算结果.Poinsot & S.K.Lele曾建议一种NSCBC方法,被广泛采用.但它是建立在局部一维无粘关系假定上的,而靠近壁面处速度梯度大,粘性作用不能忽略,上述假定不一定成立.此外,NSCBC方法在亚声速区需要下游的某些已知条件,而这往往是无法预先知道的,从而不得不做一些假定.本文将分析这样做可能出现的问题,并对某一类问题,提出一种可以得到更好结果的方法.
在超声速边界层的计算中,出口处边界条件的好坏会影响计算结果.Poinsot & S.K.Lele曾建议一种NSCBC方法,被广泛采用.但它是建立在局部一维无粘关系假定上的,而靠近壁面处速度梯度大,粘性作用不能忽略,上述假定不一定成立.此外,NSCBC方法在亚声速区需要下游的某些已知条件,而这往往是无法预先知道的,从而不得不做一些假定.本文将分析这样做可能出现的问题,并对某一类问题,提出一种可以得到更好结果的方法.
2004, 22(4): 466-469.
摘要：
在国内首座非定常风洞中,研究了不同幅值、不同频率脉动的非定常自由来流对飞机模型气动特性的影响.结果表明,当来流减速时,升力系数增加,当来流加速时,升力系数减小;随着迎角的增大,非定常自由来流对升力系数的影响增大;当来流脉动频率一定时,随脉动幅值的增大,升力系数增大;当脉动幅值一定时,随脉动频率增大,升力系数增大.因此在研究飞机机动飞行的同时,考虑飞行速度变化的影响是十分必要的.
在国内首座非定常风洞中,研究了不同幅值、不同频率脉动的非定常自由来流对飞机模型气动特性的影响.结果表明,当来流减速时,升力系数增加,当来流加速时,升力系数减小;随着迎角的增大,非定常自由来流对升力系数的影响增大;当来流脉动频率一定时,随脉动幅值的增大,升力系数增大;当脉动幅值一定时,随脉动频率增大,升力系数增大.因此在研究飞机机动飞行的同时,考虑飞行速度变化的影响是十分必要的.
2004, 22(4): 470-474.
摘要：
在非结构网格中用LU-SGS隐式算法求解欧拉方程时,两个近似分解因子的项数可能相等,为平衡;也可能不相等,为非平衡.采用对网格重新编号的方法很难达到我们希望的平衡.本文对非平衡性的影响进行了探讨.对二维问题,四边形单元的非结构网格,设计出平衡与非平衡的编号方式.先对标量模型方程分析LU-SGS隐式算法的增长因子,然后通过数值试验来验证这种非平衡性的影响.结果表明,尽管非平衡时也能达到收敛,但平衡却远优于非平衡的情况.
在非结构网格中用LU-SGS隐式算法求解欧拉方程时,两个近似分解因子的项数可能相等,为平衡;也可能不相等,为非平衡.采用对网格重新编号的方法很难达到我们希望的平衡.本文对非平衡性的影响进行了探讨.对二维问题,四边形单元的非结构网格,设计出平衡与非平衡的编号方式.先对标量模型方程分析LU-SGS隐式算法的增长因子,然后通过数值试验来验证这种非平衡性的影响.结果表明,尽管非平衡时也能达到收敛,但平衡却远优于非平衡的情况.
2004, 22(4): 475-480,.
摘要：
本文采用二维数值模拟对脉冲爆轰发动机点火方式中的预爆轰点火进行了分析研究.模拟了预爆轰管在主爆轰管上的两种最基本的布置方式:平行布置方式和正交布置方式.讨论了预爆轰管的布置方式和几何形状对点火的影响.结果表明:正交布置方式优于水平布置方式,因为反射激波在点火过程中起到了非常重要的作用.预爆轰管和主爆轰管的内径比是影响点火的重要参数,而长度的影响则较小.
本文采用二维数值模拟对脉冲爆轰发动机点火方式中的预爆轰点火进行了分析研究.模拟了预爆轰管在主爆轰管上的两种最基本的布置方式:平行布置方式和正交布置方式.讨论了预爆轰管的布置方式和几何形状对点火的影响.结果表明:正交布置方式优于水平布置方式,因为反射激波在点火过程中起到了非常重要的作用.预爆轰管和主爆轰管的内径比是影响点火的重要参数,而长度的影响则较小.
2004, 22(4): 481-485.
摘要：
在大迎角条件下的数值模拟,计算的结果往往是要受多方面因素的综合影响.本文对三角翼大迎角的定常绕流进行一系列数值模拟,详细研究了计算网格等因素对模拟结果的影响.得到了一些有益的结论.
在大迎角条件下的数值模拟,计算的结果往往是要受多方面因素的综合影响.本文对三角翼大迎角的定常绕流进行一系列数值模拟,详细研究了计算网格等因素对模拟结果的影响.得到了一些有益的结论.
2004, 22(4): 486-489,.
摘要：
本文采用数值计算方法,对细长尖拱旋成体大迎角背风侧非对称涡系结构以及与沿轴向交变的侧向力分布的关系进行了研究.求解的是N-S方程,采用Jameson中心格式,湍流采用修正的B-L模型.通过数值模拟揭示了在头部开始产生非对称涡,形成二次涡的机理以及在背风侧形成集中涡的过程.从而说明了采用集中涡来模拟这类流场的正确性,表明了二次涡对多个脱体集中涡形成的贡献.
本文采用数值计算方法,对细长尖拱旋成体大迎角背风侧非对称涡系结构以及与沿轴向交变的侧向力分布的关系进行了研究.求解的是N-S方程,采用Jameson中心格式,湍流采用修正的B-L模型.通过数值模拟揭示了在头部开始产生非对称涡,形成二次涡的机理以及在背风侧形成集中涡的过程.从而说明了采用集中涡来模拟这类流场的正确性,表明了二次涡对多个脱体集中涡形成的贡献.
2004, 22(4): 490-493.
摘要：
本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.
本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.
2004, 22(4): 494-498,.
摘要：
本文在介绍模型表面热流测量方法的基础上,重点介绍了红外热图技术测量的优点,并对红外测温原理、高超声速低密度风洞、红外热像仪、试验模型、数据处理方法等进行了叙述.在M∞＝16,P0＝1500kPa,T0＝923K,α＝0°、30°迎角的试验条件下,获得了半球圆柱模型、升力体模型表面对流换热系数分布,并把半球圆柱模型的红外测热结果与Lees分布结果进行了分析比较,结果表明当40°＜θ＜80°时,半球圆柱模型红外测热结果与Lees分布结果符合得比较好.最后对对流换热系数测量不确定度进行了初步分析.升力体模型对流
本文在介绍模型表面热流测量方法的基础上,重点介绍了红外热图技术测量的优点,并对红外测温原理、高超声速低密度风洞、红外热像仪、试验模型、数据处理方法等进行了叙述.在M∞＝16,P0＝1500kPa,T0＝923K,α＝0°、30°迎角的试验条件下,获得了半球圆柱模型、升力体模型表面对流换热系数分布,并把半球圆柱模型的红外测热结果与Lees分布结果进行了分析比较,结果表明当40°＜θ＜80°时,半球圆柱模型红外测热结果与Lees分布结果符合得比较好.最后对对流换热系数测量不确定度进行了初步分析.升力体模型对流
2004, 22(4): 499-503.
摘要：
本文从湍流模型和数值格式两方面出发,研究了对绕方形柱体的强分离低速湍流的数值模拟.比较了标准k-ε两方程模型和J-B模型,将QUICK格式推广应用于非均匀网格,采用LU分解的强隐式求解技术(LU-SIP)保证了求解的稳定性和效率.计算表明 J-B模型具有明显的优越性,与QUICK格式联合应用可十分显著地改善计算结果.
本文从湍流模型和数值格式两方面出发,研究了对绕方形柱体的强分离低速湍流的数值模拟.比较了标准k-ε两方程模型和J-B模型,将QUICK格式推广应用于非均匀网格,采用LU分解的强隐式求解技术(LU-SIP)保证了求解的稳定性和效率.计算表明 J-B模型具有明显的优越性,与QUICK格式联合应用可十分显著地改善计算结果.
