共查询到20条相似文献,搜索用时 359 毫秒
1.
为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。 相似文献
2.
3点逐步r型自适应网格算法 总被引:1,自引:0,他引:1
为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率. 相似文献
3.
针对动边界流固耦合的数值模拟问题,基于浸入有限元方法提出了一种耦合流场特征和几何特征的笛卡儿网格局部加密自适应方法,克服了单个自适应指示因子无法精确捕捉固体运动的特征的不足。在耦合自适应策略中,分别以流场涡量和固体位置作为流场和几何信息指示因子来驱动网格自适应。通过方腔顶盖驱动圆盘流动算例,以圆盘体积守恒和特征点的运动轨迹验证耦合自适应方法的优势。计算结果表明:仅基于流动特征的自适应不能很好地保证圆盘的体积守恒;仅基于几何特征的自适应无法有效追踪圆盘的轨迹;而耦合自适应策略能同时较好地保证两项指标的计算精度,在保证总体计算自由度不变的情况下,圆盘区域速度散度2-范数降低了一个数量级,圆盘的轨迹误差2-范数降低了2个数量级。 相似文献
4.
在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波. 相似文献
5.
迎风格式在二维非结构网格中的应用 总被引:2,自引:0,他引:2
将原来在结构网格中提出的AUSM+ (AdvectionUpstreamSplittingMethod)格式 ,推广应用于二维非结构网格中 .利用格林 高斯公式对控制体内的变量进行线性重构 ,获得空间二阶精度 ,并采用Barth型限制器以抑制数值解的振荡 .为提高计算效率 ,采用了自适应多重网格法 .最后给出了NACA0 0 1 2翼型和平板激波反射的几个Euler方程的算例 ,并进行了分析比较 . 相似文献
6.
超声速民机较高的飞行速度导致发动机短舱与机翼机身存在强烈的气动干扰,且高速喷流会改变主流的激波系结构,进而影响声爆强度。以超声速民机为研究对象,对不同短舱布局开展数值模拟研究,揭示发动机短舱位置、数量对近场过压信号特性的影响规律。数值模拟运用有限体积法在直角网格上求解流体控制方程进行,通过伴随自适应网格加密生成高分辨率网格以精确捕捉近场过压信号。结果表明:发动机进气道唇口激波、尾喷口后缘激波及尾喷流与机体机翼引起的激波系存在强烈的相互干扰,一定程度上增加了近场过压幅值,从而增强了声爆强度。短舱沿弦向前移及置于机翼上方可以有效降低声爆强度,沿展向外移通过抑制尾部各激波的合并也可有效降低声爆强度;在相同总推力前提下,相比三发构型,双发构型能有效降低后体激波强度,但较大尺寸的短舱引起较强的进气道唇口激波。综合考虑喷流噪声和气动阻力因素,翼下双发布局对于新一代超声速民机并非最佳选项。 相似文献
7.
将原来在结构网格中提出的AUSM+(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较. 相似文献
8.
采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础. 相似文献
9.
两种流体控制方案矢量喷管内流场计算及分析 总被引:9,自引:2,他引:7
在激波诱导和喉道倾斜2种流体控制方案下,对推力矢量喷管的二维内流场进行了数值模拟.用局部加密的四边形结构网格对流场进行网格划分之后,采用二阶迎风离散格式和二方程湍流模型求解强守恒形式的Navier-Stokes方程,通过对内流场的模拟获得流体控制推力矢量喷管方案的推力矢量特性.计算结果表明:对于激波诱导矢量控制方案,在小落压比大注气量且后缝注气的情况下所产生的矢量角最大;对于喉道倾斜矢量控制方案,在扩张片上以与主流成某一角度的方向注入气流产生的矢量角最大且没有带来较大的推力损失;喉道倾斜矢量控制方案因其推力损失较小且能够降低喷管的重量与造价而更具发展前景. 相似文献
跨声速副翼效率一直是静弹分析领域的热点和难点问题之一。目前,基于计算流体力学(CFD)/计算结构动力学(CSD)耦合的高精度静弹分析方法用于此类问题时还存在网格变形鲁棒性以及分析结果缺乏有效验证等问题。针对上述问题,提出了基于虚拟网格及虚拟位移的网格变形方法,对于迭代中出现的非物理振荡、非一致收敛问题,采用了松弛迭代以及部件载荷综合残差收敛方法。基于上述方法,分析了某型战斗机的跨声速(Ma=0.95)副翼效率,给出了静弹变形对翼面激波位置、激波强度、压力分布的影响以及副翼效率的弹性修正系数。为验证分析结果,开展了静弹试飞辨识,两者吻合良好,表明本文所提方法可以满足复杂构型跨声速副翼效率高精度静弹分析的需求,对于提高静弹工程设计能力具有重要意义。 相似文献
11.
LU-SGS因有较高的鲁棒性和小的内存需求而得到广泛应用,然而用于混合网格计算前需要进行网格排序和分组来实现算法并行;此外,LU-SGS格式收敛效率不高。针对这些缺点,本文提出了一种改进的适用于复杂混合网格的雅可比迭代方法,无需网格排序和分组就可实现算法的并行化,且有较快的收敛速度。该方法编程实现简单,易于采用OpenMP实现并行。算例研究表明,相比于LU-SGS格式,在各来流条件下,本文提出的方法收敛速度更快,鲁棒性好,并行和串行结果一致,且内存需求增加很少。 相似文献
12.