研制了用5阶FD-WENO格式(WEN05)及2阶Godunov格式(MUSCL)求解舣曲守恒律组的应用软件。通过求懈若干 Riemann问题及较复杂的一维激波相互碰撞问题对这些软什进行测试和定量比较.发现对于Sod Riemann问题,两种格式都易于算出具有较高精度和较高分辨率的数值结粜。
2023-04-10 12:03:31 3.38MB 自然科学 论文
1
威诺 WENO(加权基本非振荡)方案
2023-03-22 15:10:17 61KB MATLAB
1
朱莉娅-韦诺 任意阶的 WENO 重建算法的优化 Julia 实现。基于 Dumbser、Hidalgo 和 Zanotti 在非保守双曲系统的高阶时空自适应 WENO 有限体积方案(DOI 10.1016/j.cma.2013.09.022) 中的工作。 用法 u您的输入数据必须采用 shape数组的形式(nx,ny,nz,nvar)。 nvar是每个单元格中包含的变量数,并且nx,ny,nz分别是x、y、z 轴中的单元格数。如果您的网格是 1D 或 2D,请根据需要设置ny=1和/或nz=1。选择整数N,其中N+1是所需的精度顺序。WENO 重建的系数通过调用获得: julia> weno(u,N); 该调用每次都会构造 WENO 系数矩阵M1,M2,M3,M4和振荡指标矩阵。Σ如果需要许多 WENO 重建,则预先计算这些实体的计算效率会更高: julia> M1,M2,M3,M4 = coefficient_matrices(N); julia> Σ = oscillation_indicator(N); julia> chΣT = chol(Σ)'; julia> wen
2022-06-10 09:06:27 5KB julia 算法
双曲性守恒方程组采用高阶、高分辨率的WENO格式时有两类分裂方法,即逐点分裂和特征分裂。本文基于后者,对特征分裂重构时强间断和接触间断位置出现的振荡情况进行研究,对重构变量加以改进,发现改进后的WENO格式克服了间断处的振荡,然后以LU-SGS为子迭代的双时间步法求解Euler方程,选用一维Sod、二维前台阶和双马赫反射算例,并与Roe格式计算结果进行对比,发现WENO格式分辨率更高,耗散更小。
2022-05-17 16:56:24 998KB 自然科学 论文
1
一维和二维域中线性对流方程的 WENO(加权基本非振荡)方案。
2022-05-16 22:00:30 46KB matlab
1
有限差分方法,利用WENO重构求解二维Euler方程,包含WENOWENO-Z、WENO-ZN等格式 大学生课程设计 课程设计 自己大二写的课程设计
2022-04-22 18:05:10 19KB 重构
weno代码
2022-04-17 16:02:52 287KB weno
1
由 介绍的五阶 WENO 方案的一维实现 [1] 舒志旺。 “本质上是非振荡的,加权本质上是非振荡的双曲线守恒定律。” 非线性双曲方程的高级数值逼近。 斯普林格,柏林,海德堡,1998 年。325-432。 和 [2]江,广山和吴成钦。 “理想磁流体动力学方程的高阶 WENO 有限差分格式。” 计算物理学杂志 150.2 (1999): 561-594。 本代码旨在作为该方法实现的指南。 它举例说明了有限差分 (FD) 和有限体积 (FV) 方法的按分量重建的实现。 在这次更新中,我还包括 FV 方法中的特征重建。 与往常一样,此代码背后的理念是可读而不是高效。 在这里,我将这个例子献给所有开始数值方法之路的 CFD 学生。 曼努埃尔·A·迪亚兹(2018 年 6 月)
2022-04-07 20:37:53 103KB matlab
1
WENO,TVD等差分格式解欧拉方程验证黎曼问题和双马赫反射问题的求解程序 Matlab。 算法种类多,移植方便,直接运行
2022-03-30 21:46:52 24.6MB CFD weno
1
此程序用于计算二维黎曼问题,采用五阶WENO格式空间离散和三阶TVD Runge-Kutta时间离散,同时采用LF分裂进行求解,生成的数据可以直接导入tecplot中进行绘图