共查询到20条相似文献,搜索用时 31 毫秒
1.
提出Lagrange柱坐标高阶中心型守恒格式.基于用对守恒律的单调迎风算法(MUSCL)构造的高阶子网格压力,引入了柱坐标高阶体权子网格力和柱坐标高阶面权子网格力,构造了柱坐标高阶体权中心型守恒格式和柱坐标高阶面权中心型格式.柱坐标高阶体权中心型守恒格式满足动量守恒、能量守恒,但不能确定保持一维球对称性.柱坐标高阶面权中心型格式满足能量守恒,保持一维球对称性.两种格式里,格点速度以与网格面的数值通量相容的方式计算.对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange柱坐标高阶中心型守恒格式的有效性和精确性. 相似文献
2.
提出了Lagrange中心型守恒气体动力学格式.引入了当前时刻子网格密度与当前时刻网格声速产生的网格分片常数压力.初始网格密度乘以初始子网格体积得到子网格质量,这些子网格质量除以当前时刻子网格体积得到当前时刻子网格密度.应用网格分片常数压力,构造了满足动量守恒、总能量守恒的Lagrange中心型守恒气体动力学格式,格点速度以与网格面的数值通量相容的方式计算.对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange中心型守恒气体动力学格式的有效性和精确性. 相似文献
3.
一类交错网格的Gauss型格式 总被引:1,自引:0,他引:1
本文在交错网格的情况下 ,利用 Gauss型求积公式构造了一类不需解 Riemann问题的求解一维单个双曲守恒律的二阶显式 Gauss型差分格式 ,证明了该格式在CFL条件限制下为 TVD格式 ,并证明了这类格式的收敛性 ,然后将格式推广到方程组的情形 .由于在交错网格的情况下构造的这类差分格式 ,不需要求解 Riemann问题 ,因此这类格式与诸如 Harten等的 TVD格式相比具有如下优点 :由于不需要完整的特征向量系 ,因此可用于求解弱双曲方程组 ,计算更快、编程更加简便等 . 相似文献
4.
在交错网格型Lagrange(拉格朗日)流体力学算法中,通常采用人为粘性捕捉激波,人为粘性的好坏对于计算结果至关重要.研究了一种基于子网格边界处近似Riemann解的新型人为粘性.新人为粘性能够满足动量守恒和熵不等式.利用子网格边界速度差中引入的限制器,新人为粘性能够区分激波和等熵压缩,并能满足球对称问题中的波面不变性.新人为粘性在典型模型数值模拟及惯性约束聚变黑腔整体数值模拟中取得了较好的结果. 相似文献
5.
基于WENO(Weighted Essentially Non-Oscillatory)的思想,提出了一种在非结构网格上求解二维Hamilton-Jacobi(简称H-J)方程的数值方法.该方法利用Abgrall提出的数值通量,在每个三角形单元上构造三次加权插值多项式,得到了一个求解H-J方程的高阶精度格式.数值实验结果表明,该方法计算速度较快,具有较高的精度,而且对导数间断有较高的分辨率. 相似文献
6.
本文讨论时间导数项含小参数的抛物方程.我们依Бахволов构造非均匀网格的差分格式,并证明了格式的一阶一致收敛性.给出了数值结果. 相似文献
7.
8.
本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对,格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量,但不能唯一确定方程组的熵守恒通量,却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响. 相似文献
9.
10.
11.
本文在星形多边形网格上, 构造了扩散方程新的单调有限体积格式.该格式与现有的基于非线性两点流的单调格式的主要区别是, 在网格边的法向流离散模板中包含当前边上的点, 在推导离散法向流的表达式时采用了定义于当前边上的辅助未知量, 这样既可适应网格几何大变形, 同时又兼顾了当前网格边上物理量的变化. 在光滑解情形证明了离散法向流的相容性.对于具有强各向异性、非均匀张量扩散系数的扩散方程, 证明了新格式是单调的, 即格式可以保持解析解的正性. 数值结果表明在扭曲网格上, 所构造的格式是局部守恒和保正的, 对光滑解有高于一阶的精度, 并且, 针对非平衡辐射限流扩散问题, 数值结果验证了新格式在计算效率和守恒精度上优于九点格式. 相似文献
12.
13.
田芳 《数学的实践与认识》2015,(4):268-275
基于非均匀网格上函数的泰勒级数展开,结合残参量修正法,推导了非均匀网格上对流扩散方程的高阶指数型紧致差分格式,选取的算例表明,格式兼有高精度和高分辨率的优点,能够很好的适用于大梯度变化,计算区域中含边界层和对流占优区域中的流动问题的求解. 相似文献
14.
高阶精度加权紧致非线性格式(WCNS)越来越广泛地应用于复杂流动数值模拟.WCNS可以与多种无粘通量分裂方法结合起来使用.但是,常见的通量分裂方法都是基于低阶格式发展起来的,目前还不清楚哪些通量分裂方法最适合WCNS,也不知道这些方法与高阶格式结合时将会产生什么效果.表面热流计算是高超声速流动数值模拟的难点之一,为了在热流计算时选择合适的通量,研究了多种通量分裂方法的耗散大小.每种通量都可以表示成中心部分与耗散部分之和.这些通量的中心部分相同且非常简单,但是耗散部分较为复杂,且不同的通量分裂方法可导致不同的耗散表达式.通过对通量耗散进行分析可以发现耗散大小与网格界面两侧的物理量跳跃近似线性正相关.数值计算表明高阶格式得到的网格界面左右两侧的物理量跳跃通常远比低阶格式小,因而带来的通量耗散小.通过3个典型算例考察了通量耗散对热流计算的影响,其中包括高超激波/边界层干扰算例.基于对van Leer通量、Steger-Warming通量、KFVS通量、Roe通量、AUSM类通量和HLL类通量的考察,给出了通量选择建议. 相似文献
15.
16.
17.
激波捕捉差分方法研究 总被引:1,自引:1,他引:0
在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试.该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式.用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试.通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率. 相似文献
18.
提出了一种基于非结构自适应网格的二维Euler方程的数值解法.采用有限体积法进行空间离散,通量计算采用Jamson中心格式,使得它适用于任意多边形计算单元.为了得到定常解,采用一种显式的四步Runge-Kutta迭代方法对时间进行积分.根据流场参数的变化梯度确定加密边,由加密准则进行自适应网格剖分,然后得到分布合理的加密过后的网格.求解二维Euler方程,对NACA0012翼型进行了数值模拟,通过对自适应前后的数值解的对比,说明所建立的方法是正确的. 相似文献
19.
发展了一种基于高阶迎风格式和嵌套网格捕捉直升机悬停旋翼涡尾迹的方法.无粘通量采用Roe Reimann求解器,使用改进的5阶加权基本无振荡(WENO)格式对交界面左右状态进行高阶插值,并与MUSCL插值进行比较.为便于捕捉尾迹和实施周期性边界条件,计算采用结构嵌套网格,其中高质量的旋翼网格完全嵌套于背景网格中.当解达到近似收敛后在桨尖涡分布区域对背景网格进行加密,如此经过3次得到优化的背景网格.考虑到WENO格式插值的特点,提出了搜索3层洞边界和人工外边界的方法以便插值的直接进行.用该方法对一跨音速和一亚音速悬停旋翼粘性流场进行了数值计算.数值结果表明:所发展方法对涡尾迹具有很高的捕捉能力;与MUSCL格式相比,WENO格式具有较低的数值耗散. 相似文献
20.
本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度. 相似文献