只改了一个松弛因子,收敛解就变了 — Rhie–Chow 面速度记住欠松弛的地方
·4 min read
在 SIMPLE 中只改 α_u,收敛后的速度场就会改变。原因在于 Rhie–Chow 面速度混入前一次迭代值的方式,Majumdar 修正可以消除它。
给曲线单元填入真实面积后,∂x/∂x 不再等于 1 —— 对角线叉积面积为何与 Green–Gauss 相容
·4 min read
在曲线网格 FVM 里,把单元中心雅可比当作体积会让 ∂x/∂x 偏离 1。只有对角线叉积面积与直边面向量围成的多边形一致。
网格歪斜时梯度也跟着偏 —— Green–Gauss 与加权最小二乘重构
·4 min read
靠近壁面的扁平单元里,Green–Gauss面插值把梯度拧向直角,距离加权n=1把曲面误差从0.0293降到0.0222。
非结构网格上二阶重构制造的虚假极值 — Barth–Jespersen 与 Venkatakrishnan 限制器
·5 min read
非结构单元的二阶面值一旦越出邻居范围,Barth–Jespersen会削掉斜率,Venkatakrishnan则用K把折点磨圆。
不分裂方向求解对流 — Bell-Colella-Glaz 非分裂迎风格式
·5 min read
方向分裂会把对角对流拧歪。BCG预测子用(1-CFL)斜率权重,再加一项横向修正,把旋转圆盘的峰值留在网格上,而不是靠加密。
任意方向上的欧拉特征结构 — 三族波、两重重根、一处奇点
·4 min read
此时当 n_x 为 0.01 时,左特征向量的分母正好是 100。30° 的面用 R-1,法向躺在 y 上时改选 R-2。
在非结构网格上拯救梯度的权重 — 移动最小二乘(MLS)近似
·5 min read
默认 h=0.07、25 个散点、二阶这一档里,sin 导数 RMSE 为 0.21,大约只占振幅 6.28 的 3%。
歪斜单元让梯度走偏 — Green–Gauss 与 Least-Squares 重构对比
·5 min read
skew为0.9时,打印出的Green–Gauss误差是31%,最小二乘是9%。从对称环开始,比值就停在大约3.4附近。
FDM vs FEM vs FVM:三种离散化技术的本质区别
·6 min read
同一平流-扩散方程下,FDM离散微分算子,FEM把解投到基函数,FVM平衡界面通量,守恒、精度阶和非结构网格因此分道扬镳。