共查询到20条相似文献,搜索用时 31 毫秒
1.
本文工作是在 Mac Cormack通量分裂格式 [2 ] 基础上发展的一种有限面积通量分裂隐式格式 ,其特点是在求解隐式离散化方程时 ,采用往返扫描一次的 Gauss- Seidel线松驰迭代方法 ,避免了对时间步长的任何限制。为提高定常解精度 ,格式的显式右端项采用二阶精度的离散。在数值求解跨音速涡轮平面叶栅问题中 ,对壁面边界作了较仔细的隐式处理。数值计算表明本方法保持了 Mac Cormack格式具有的高收敛速率 ,(约 30个时间步即可达到定常解 )而每一时间步计算量约减少一半。数值结果与实验结果符合得很好。 相似文献
2.
提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。 相似文献
3.
一、引 言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法… 相似文献
4.
5.
本文发展了二维非定常不可压N-S方程的局部网格加密方法(LMR):需加密区域预先指定;仅对空间方向加密,分粗细二层网格,细网格覆盖在粗网格上;粗和细网格上分别用显式和隐式差分格式。结合压力修正法类的余量型差分格式,恰当地处理了粗细网格之间的信息传递,使得粗细网格交接面上质量守恒,即满足连续性条件。粗网格通过插值给细网格边界值影响细网格;细网格通过粗细网格压力Poisson方程耦合求解影响粗网格,并且由于压力Poisson方程从动量方程形成,在细网格覆盖下的粗网格上,压力Poisson方程是细网格方程的组合,从而耦合求解时粗细网格压力Poisson方程不需迭代。 本文中计算了二维方腔内的自然对流模型问题,Rayleigh数是10~6。粗网格用显式ULWC格式,细网格用隐式余量型近似因式分解格式。Poisson方程用快速直接算法PO1STG(在FISHPACK中)。粗网格距△=1/16,细网格距△=1/64,加密1/4区域,计算到定常局部网格加密方法所需CPU时间比同等均匀网格(△=1/64)节约一倍以上。 相似文献
6.
7.
小轿车绕流场数值计算和分析 总被引:8,自引:0,他引:8
本文求解了雷诺平均Navier-Stokes方程,湍流的模拟采用了标准的两方程模型,修正的两方程模型,亚格子模型的大涡模拟三种要用有限体积方法离散求解Navier-Stokes方程。对流项的计算采用三阶迎风格式。为了提高求解过程的稳定性,采用显隐式混合格式来获得主对角元素区域占优。 相似文献
8.
用CSCM格式求解激波反射问题 总被引:1,自引:0,他引:1
本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。 相似文献
9.
结构动力学方程的显式与隐式数值计算 总被引:2,自引:0,他引:2
从理论上详细推导了动力学方程的显式隐式时间积分数值计算方法,针对经典中心差分格式在求解阻尼结构计算效率降低的问题,研究提出了一种新的显式积分格式,理论证明该数值积分格式能够处理各种复杂阻尼状态而不降低计算效率。对于不同的时间步长通过设定调节加速度参数,能得到比经典显式计算格式和NEWMARK逐步时间积分更高的精度,最后通过数值算例对三种计算格式的理论分析进行了验证。 相似文献
10.
11.
12.
本文介绍一种求解涡度流函数形式的定常不可压缩Navier-Stokes方程的半隐式指数型差分格式。涡度方程用上述格式求解,而流函数方程用多层网格法求解。这种组合求解方式具有很好的稳定性及较快的收敛速度。本文对Re=100,400,1000,3000,5000的驱动方腔及Re=400,1000的驱动长方腔进行了数值模拟,结果与现有的计算相吻合。 相似文献
13.
用三维N—S方程计算钝体绕流 总被引:8,自引:5,他引:3
本文采用有限体积方法求解三维低速N-S方程,研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组,采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法,采用面松驰迭代来实现方程的求解,湍流的计算采用了标准的双方程模型,本文给出了钝体和类小汽车外形的计算结果,钝体的计算结果和实验结果进行了对比,与其它方法的计算结果也进行了对比,它们之 相似文献
14.
首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法(FRLBM)。采用两种求解方法:一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程(直接法);另一种是先执行碰撞步,再求解纯对流方程(分步法)。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。 相似文献
15.
16.
[1]和[2]曾用几种差分格式求解一维N-S方程和欧拉方程,计算表明,Hopscotch格式(跳点格式)计算速度最快。为探讨这种格式在多维流场特别是跨音速流场计算中应用的效果,本文将此差分格式应用于有限体上,成为跳点有限体法,并用该方法求解二维欧拉方程组。对若干个喷管和叶栅所做的计算得到了满意的结果,计算时间较目前常用的差分格式少得多。 相似文献
17.
本文采用算子附加修正方法,提出了一个气动计算中简单而有效的差分格式。此差分格式的特点是保持原显式格式的精度,计算过程有显式格式的简单,又有隐式格式的特点,即对于任意步长此格式都是稳定的。本文利用此格式求解了N-S方程,给出了二维激波与边界层干扰的分离流计算结果,并与MacCormack的隐式格式进行了比较。计算结果表明,本文提出的格式与MacCormack隐式格式相比,大大简化了计算过程,节省了计算机时。 相似文献
18.
19.
二维粘性不可压缩流场的分块耦合求解方法 总被引:1,自引:1,他引:1
本文用分块耦合的方法,并采用S.Abdallah^[1,2]和M.L.Mansour^[3]的在非交错网络上隐式求解原始变量的Navier-Stokes方程的提法和离散格式,成功地解决了二维复杂流场的数值求解问题,获得了空腔流、T型叉管流、NACA0020翼型绕流和方块钝体绕流的稳态解,并进行了讨论。此方法还可推广到三维复杂流场的计算问题。 相似文献
20.
DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。 相似文献