共查询到20条相似文献,搜索用时 331 毫秒
1.
本文采用气动中心张涵信教授提出的无波动、无自由参数耗散格式(NND)进行了如下的计算:航天飞机头部段对称剖面内的二维无粘绕流;航天飞机头部段三维无粘绕流;轴对称凹陷外形的三维无粘绕流。在物面坐标系下,用代数方程生成计算网格,利用特征相容关系式处理激波和物面边界。为了比较,还用SCM格式进行了计算。结果表明,本文所采用的NND格式,在捕捉流场内激波时不会出现非物理振荡,其收敛速度和所获得的计算结果都是令人满意的。 相似文献
2.
3.
本文给出了用NS方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全NS方程;湍流模型采用双层代数涡粘性模型。使用近似因式分解ADI差分格式离散求解,网格是用保角变换方法生成的相对翼型固定的C型网格。本文给出两类典型非定常绕流数值模拟结果:翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较,表明了本方法准确、高效的特点。 相似文献
4.
二维非定常欧拉方程的有限元解法 总被引:1,自引:0,他引:1
本文通过求解时间相关的欧拉方程,对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法,并用VanLeer逆风通量分裂格式来计算控制体通量。隐式格式的采用,大大提高了计算效率;同时,通过插值使得空间离散达到了二阶精度。另外,本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动,从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了NACA0012翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好, 相似文献
5.
6.
由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。 相似文献
7.
8.
激波-边界层-分离流相互干扰三维湍流的数值模拟 总被引:1,自引:0,他引:1
本文采用数值方法求解时间相关三维可压缩雷诺平均Navier-Stokes方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Badwin-Lomax两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,VanLeer矢通量格式计算无粘通量,中心差分法计算粘性通量,LUSGS时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。 相似文献
9.
风对建筑绕流流动的数值方法 总被引:4,自引:0,他引:4
本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。 相似文献
10.
11.
对一种新式的风机叶片即轴向前倾叶片流场进行数值计算, 采用时间推进法计算无粘流场, 采用非正交曲线坐标系附面层积分方程求解附面层参数, 得到了一些有意义的结果。 相似文献
12.
采用具有7阶精度的weighted essentially non-oscillatory(WENO)差分格式,直接求解可压缩二维非定常N-S方程组,研究了NACA0012翼型平面叶栅低雷诺数流动的特征.直接模拟及与文献对比的结果表明:叶栅尾缘涡脱落的形成过程与圆柱绕流涡脱落的形成过程非常相似.平面叶栅尾迹区的2阶统计量与孤立翼型尾迹区的2阶统计量具有相同的分布特征,但前者的强度显著大于后者.周期性的涡脱落不仅在上下翼面形成非定常分离,也使得尾迹区某点的总压发生准周期性的变化.随着栅距的减小,翼型上的平均分离位置向前缘移动;尾迹区某点的总压变化频率及其幅值均显著地增加;而且栅距越小,速度脉动2阶统计量反而越大. 相似文献
13.
利用自主研发的混合网格Euler/RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。 相似文献
14.
15.
以具有压力分裂形式的简化N S方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制(passivecontrolofshock boundarylayerinteraction)。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。 相似文献
16.
二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解 总被引:1,自引:1,他引:0
为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍. 相似文献
17.
钝锥大攻角超声速分离流场的数值模拟 总被引:2,自引:2,他引:2
最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。 相似文献
18.
研究飞机外挂物在按给定轨迹下沉旋转的离机过渡过程中,流场非结构网格的调整和再生成方法以及过渡过程欧拉方程有限元数值解法。中指出外挂物运动时流场网格调整或再生成仅需在围绕在外挂物的一小区内进行以减少计算量,而且固定网格的欧拉方程的有限元数值解法也可适用于运动网格情况,只要将固定网格的守恒通量项Fk用运动网格的Ek来代替,并建立了Fk和Ek之间的关系。末给出了带外挂物飞机的实例计算结果。 相似文献
19.
20.
具有后挡板的喷流流场的数值研究 总被引:1,自引:0,他引:1
本文对发射台底面与后挡板之间的喷流流场进行了数值研究。采用迎风的Petrov-Galerkin有限元方法求解欧拉方程。描述了喷流出口处及挡板附近的流场细节。计算结果已成功地应用于导弹发射台的设计。 相似文献