请问一下关于利用OpenFOAM做海洋内波模拟
-
目前我的做法就是，基于waves2Foam的框架，对其进行扩展，利用kdV方程的解模拟海洋内孤立波。这个跟普通的孤立波相当近似，但是上层的气相变成了密度稍低的水，并且上层的低密度水跟下层的高密度水互不相容（可以把上层的水当做油来看）。
公式参考的是：这两位的论文
个人主要就是改了一下几点：
TransportProperties 里 气相相关属性替换为液相waves2Foam波浪模块，原先只操作液相网格不操作气相，现在两层流体需要对上层的网格也进行操作，对其速度进行赋值。
解决完以上问题后，进行setFields，得到结果如上图所示：
对比论文，这个水质点速度分布我认为是比较合理的，然而.....如果用waveFoam(本质上还是Interfoam，不过多了个速度入口的波浪边界条件)开始计算，结果发现速度会发生突变，导致波面很快变形消失。
t=0.001时刻的速度分布。
、这里我尝试了将波浪方程中的时变项全部去掉，或者把时间步长取得再小一点，都遇到了这种突变。
个人猜测是interFoam进行多相流计算的时候还有什么关键点我没考虑到，或者是interFoam只可以算气液多相流？
请各位赐教。
