首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 156 毫秒
1.
采用流体体积分数的混合型多流体数值模型和高精度多介质Piecewise Parabolic Method(PPM)算法模拟一般形式Mie-Grüneisen状态方程的复杂多流体可压缩流动.根据多流体接触界面无振荡原则设计高精度计算格式,采用双波近似和两层迭代算法求解一般状态方程的Riemann问题,利用Strang分裂方法将一维算法直接推广到高维情况,不需要动能校正.二维多流体接触间断平移问题的模拟结果表明本算法在接触间断附近压力和速度没有振荡,界面数值耗散很小,被控制在两到三个网格之内.本文还模拟了二维Riemann问题和激波同多流体交界面的相互作用,算例结果表明本文方法可以有效地处理接触间断、激波等物理问题,且具有耗散小精度高的特点.  相似文献   

2.
一种求解激波问题的中心差分-WENO混合方法研究   总被引:1,自引:1,他引:0       下载免费PDF全文
为提高求解包含激波问题的计算精度和效率,发展了一种低耗散、高效率的中心差分-WENO(Weighted Essentially Non-Oscillatory)混合格式,Navier-Stokes方程的无黏项在光滑流场区域采用六阶中心差分格式离散,而间断附近采用五阶WENO格式求解;基于密度设计了一种新型格式开关实现两种格式在光滑区域和间断之间自动切换,确保数值解在间断附近基本无振荡。针对一维激波熵波作用问题验证了混合格式的低耗散特性;对二维Riemann问题的研究表明发展的基于密度的格式开关更为合理,混合格式的计算效率较WENO格式明显提高;将其应用到激波诱导燃烧问题中,混合格式能很好地捕捉定常流场中激波和化学反应锋面的位置以及非定常流场的振荡频率。  相似文献   

3.
求解多维欧拉方程的二阶旋转输运格式   总被引:1,自引:0,他引:1  
提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.  相似文献   

4.
本文应用最小耗散的混合格式及MacCormack二阶格式加θq╱2耗散修正格式计算了一维激波管Riemann问题、二维激波管中激波与楔的反射流场。这两个格式的特点是对任何初值其耗散项永远不等于零,因而,计算得到的解为物理解,而且除间断外,均具有一致二阶精度,计算出的单波很精确,并且激波分辨度高。本文计算了楔角θ~*=30°、45°,激波马赫数M_s=5.29、10的情况,给出了流场的非定常过程及流场的自模解,一维结果与Riemann问题理论解、二维结果与实验结果的比较,说明计算是成功的。  相似文献   

5.
高分辨率格式在非定常无粘可压缩流动中的应用研究   总被引:1,自引:0,他引:1  
本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。  相似文献   

6.
使用高阶间断Galerkin格式求解守恒律方程组时,激波附近的Gibbs效应容易导致非物理解的产生。为抑制这一现象,必须构造合理的限制器对数值解进行处理。目前间断Galerkin格式中的限制器多源于有限体积法,在非结构网格上只对低阶导数项进行限制,对高阶导数项则很难给出普适判据。文章对间断Galerkin解进行广义Fourier展开,实现不同频域范围内的谱分解;在新的模板坐标系下描述各阶方向导数的变化规律;结合当前单元和相邻单元的信息,分层限制各阶方向导数,实现对非物理解的抑制。通过求解Euler方程,对二维Riemann问题、翼型周围的亚、跨声速流动问题、前台阶问题以及超燃冲压发动机内流场激波反射问题进行数值模拟,检验了新型限制器的可靠性以及向高阶格式推广的可行性。  相似文献   

7.
本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。  相似文献   

8.
对于包含激波、剪切层等复杂结构的流动问题,为了精确模拟剪切层等精细结构,且保证激波计算的稳定性,必须采用低耗散且强鲁棒的数值通量方法。传统的HLL近似Riemann求解器的耗散性较大,Roe、HLLEM和HLLC等近似Riemann求解器在计算某些含有强激波的物理问题时会出现非物理解,容易导致不稳定。针对这一问题,本文在Riemann求解器中通过合理设计反扩散矩阵,发展了一类具有自适应反扩散的新型Riemann求解器,并将其应用到高阶加权紧致格式,实现了高阶精度求解。通过典型数值算例验证了新型方法的计算精度和稳定性,结果表明本文提出的新型自适应反扩散Riemann求解器克服了传统Riemann求解器的缺陷,既能准确识别剪切层等精细结构,又能保证激波解的稳定性。  相似文献   

9.
闫超  陈靓 《航空学报》1996,17(Z1):68-71
 用不同的数值格式,包括Roe的FDS格式和vanLeer的FVS格式并配备不同限制器,求解二维N-S方程,模拟研究了激波/边界层干扰的流动问题,分析了模拟出的流动现象,仔细研究了其中的格式效应,得出同前人不尽一致的结论:vanLeer格式对限制器的依赖性比Roe格式更强;在本文情况下,minmodlimiter的激波、膨胀波的分辨率比albadalimiter还好  相似文献   

10.
双重加权实质无波动激波捕捉格式的改进和应用   总被引:2,自引:1,他引:2  
本文重新研究了双重加权实质无波动(DWENO)激波捕捉格式的构造过程,对模板光滑系数的定义,以及两重权值的计算方法进行了改进,并通过一维Riemann问题测试了这些改进,确定了参数的选择范围。最后,本文将改进后的DWENO(r=3)格式进一步应用于一维、二维典型问题的数值模拟,并比较了DWENO格式和WENO格式的精度和分辨率。  相似文献   

11.
本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。  相似文献   

12.
针对再入飞行器激波层辐射加热,应用有限体积法计算辐射传输方程,引入了一种边界处理方式——虚网格方法,根据化学非平衡流场数据,计算出飞行器表面辐射加热热流。通过分析发现:采用虚网格方法处理辐射传输边界,能够确保真正的边界条件得到应用,满足真实的物理意义,且辐射传输计算的空间网格、立体角网格、初始温度值的收敛性很好,计算结果精度误差小于2.5%。利用虚网格方法计算飞船再入返回时的辐射加热,结果表明计算所得壁面辐射热流与其他研究者辐射热流计算结果符合良好。   相似文献   

13.
针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环在水下的燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波退化为激波在水中的传播及衰减过程。研究结果表明:燃气泡前期受外界水环境阻滞作用呈现“豌豆状”形态,充分发展阶段气液交界面逐渐失稳,在达到最大尺寸后开始收缩并在中心轴线位置出现凸出的射流;水下爆轰燃气射流发展过程中同时存在脱离燃气泡的水中前导激波和管口燃气泡内的高压区两部分,水中前导激波在传播过程中压力迅速衰减至常压量级,而管口燃气泡内则一直保持较高压力;中心轴线区域气液交界面处反射激波的回传使管口附近出现回击现象,并导致前导激波波阵面上压力峰值逐渐出现在30°方向上。  相似文献   

14.
在气-液界面Richtmyer-Meshkov(R-M)不稳定性研究的基础上,利用不同流体间密度的差异,构造了空气-硅油-水、空气-酒精-硅油两种流体粘度方向迥异的气-液-液三相界面,实验中以氮气作为高压驱动气体,在不同激波马赫数下对这两种界面R-M不稳定性后期尖钉与气泡区发展进行了测量与统计,通过对实验数据的分析,得出了相关规律,并用这些规律与已有的单层气-液界面R-M不稳定性研究成果作比较,得出了异同点。同时,还研究了两种三相界面在R-M不稳定发展中的差异。实验结果表明:当流体的粘度梯度方向(从小到大的方向)与激波方向一致时,界面失稳更加明显,湍流混合更为显著。  相似文献   

15.
This paper deals with data association using three sets of passive linear array sonars (PLAS) geometrically positioned in a Y-shaped configuration, fixed in an underwater environment. The data association problem is directly transformed into a 3D assignment, which is known to be NP hard. For generic passive sensors, it can be solved using conventional algorithms, while in PLAS, it becomes a formidable task due to the presence of bearing ambiguity. Thus, the central issue of the problem in PLAS is how to eliminate the bearing ambiguity without increasing tracking error. To solve this problem, the 3D assignment algorithm used the likelihood value of only those observed bearing measurements is modified by incorporating frequency information in consecutive time-aligned scans. The region of possible ghost targets is first established by the geometrical relation of PLAS with respect to target. The ghost targets are then confirmed and eliminated by generating multiple observations in consecutive scans. Representative simulations demonstrate the effectiveness of the proposed approach.  相似文献   

16.
《中国航空学报》2020,33(9):2295-2312
In this paper, a hybrid Lattice Boltzmann Flux Solver (LBFS) with an improved switch function is proposed for simulation of integrated hypersonic fluid-thermal-structural problems. In the solver, the macroscopic Navier–Stokes equations and structural heat transfer equation are discretized by the finite volume method, and the numerical fluxes at the cell interface are reconstructed by the local solution of the Boltzmann equation. To compute the numerical fluxes, two equilibrium distribution functions are introduced. One is the D1Q4 discrete velocity model for calculating the inviscid flux across the cell interface of Navier–Stokes equations, and the other is the D2Q4 model for evaluating the flux of structural energy equation. In this work, a new dual thermal resistance model is proposed to calculate the thermal properties at the fluid–solid interface. The accuracy and stability of the present hybrid solver are validated by simulating several numerical examples, including the fluid-thermal-structural problem of cylindrical leading edge. Numerical results show that the present solver can accurately predict the thermal properties of hypersonic fluid-thermal-structural problems and has the great potential for solving fluid-thermal-structural problems of long-endurance high-speed vehicles.  相似文献   

17.
伽辽金有限元法是解跨声速绕流的有效方法,对超临界流动,为了自动捕捉到激波,通常采用了上风技术,但导致了气流通过激波时质量守恒条件得不到满足。本文通过Bateman变分原理和加罚泛函找到了强制满足质量守恒约束条件的附加项,大大改善了用伽辽金有限元法解超临界绕流的计算结果。  相似文献   

18.
《中国航空学报》2021,34(2):279-287
The free piston shock tunnel is a type of shock tunnel with high performance. For this type of tunnel, the influence mechanism of shock wave attenuation on tailored operation is explored by numerical simulation and theoretical analysis. By introducing the normalized velocity, the simple constraint equation for shock wave under the tailored operation is deduced. Moreover, the real gas effect is also taken into account in this equation. Based on the equation, the tailored operation of shock tunnels can be predicted with very few calculations. The present study shows that the change rate of the thermodynamic state of the gas behind the shock wave is inconsistent with the attenuation rate of the shock wave, which is the fundamental reason why the wind tunnel achieves tailored operation at a lower Mach number of shock waves. This lower Mach number of shock waves differs from the corresponding ideal value by a factor, which is about the square root of shock attenuation rate.  相似文献   

19.
低速下水滴撞击固体表面运动模拟   总被引:2,自引:1,他引:1       下载免费PDF全文
以水滴撞击飞机部件造成结冰问题为背景,运用level set界面追踪方法,把不可压黏性方程组中的表面张力用level set函数表示,并利用2阶变密度投影算法求解二维不可压Navier-Stokes方程组,数值模拟了低速下的水气界面追踪问题.为了能够处理密度比为1 000的水气界面问题,在界面处引入了一定宽度的光滑过渡带宽,进而消除了由于界面两边大密度比引起的数值震荡,由此较为成功地模拟了大水滴在不同表面张力时撞击平板后的形态变化及其运动轨迹.数值实验结果表明本算法可以有效地处理界面处密度比和黏性很大的不可压多相流问题.考虑界面处的表面张力,本算法可模拟不同表面张力下的界面运动变化.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号