用于瞬态对流扩散PDE的简单而通用的FVM求解器
一个简单的有限体积工具
这段代码是化学/石油工程师开发一种简单的工具来求解对流扩散方程的一般形式的结果:
αgeneral / equationt +∇。(uϕ)+∇。( -D∇ϕ)+ βϕ =γ
在1D,1D轴对称(径向),2D,2D轴对称(圆柱)和3D域上的简单均匀/不均匀网格上。
该代码在整个或部分边界上接受Dirichlet,Neumann和Robin边界条件(可通过更改以下等式中的a,b和c来实现):
a∇ϕ.n + bϕ = c。
它还接受周期性边界条件。
该代码的主要目的是为那些尝试使用数学模型,以1D形式数值求解模型,将其与解析解进行比较,并以最少的数量将其数字代码扩展至2D和3D的人员提供方便的工具。一维代码中的修改。
离散化方案包括
*中心差
*对流项的迎风方案
*具有许多通量限制器的对流项的 TVD方案
首先,请转到“ Test”文件夹并运行测试脚本。
一些微积分函数(散度,梯度等)和平均技术(算术平均值,谐波平均值等)可用,特别有助于解决非线性或耦合方程式或实现显式方案。
我已使用该代码解决了PDE的耦合非线性系统。您可以在Examples / advanced文件夹中找到其中一些。
“