利用OpenFOAM开发一个求解波动方程的的求解器。
波动方程表达式为:
式中,c为波速,h为波高。
1 创建文件
利用快速创建文件结构。
文件结构如下所示。
2 源代码
头文件中实现物理量的读取
源文件
3 测试案例
采用二维模型,几何尺寸为2x2m,网格尺寸0.02 m,计算域中心位置有一个0.2x0.2m的区域,将h的值设置为1,其他区域h=0。
这里挑重要的文件进行说明,细节可直接查看案例文件。
文件中需要添加h的求解方式
字典文件进行初始化
文件中指定波速
文件指定待求量h的初始值与边界值
计算结果如下图所示。
(完毕)
领取专属 10元无门槛券
私享最新 技术干货