@李东岳

李老师您好。

目前出口是采用如下方程

这里Uc尝试过：

1. 使用上下层速度的平均值1.5。
2. 在X方向上采集最靠近边界的三个点的速度值进行平均，作为Uc。

如果采用2的话，Uc是一个随着Y坐标以及时间变换的函数Uc(y,t)

Phi是压力泊松方程构造的中间量，这个查了好多文章都说的很模糊，我也就直接用了零梯度。

各位老师好。目前我所模拟的是两层不同速度流体的混合问题。
N-S方程的求解采用的是Fractional Step Method，时间步进为三步龙格库塔法。
上下速度边界采用的是零梯度，左右采用周期性边界条件，进口指定一个速度Profile，出口采用的是对流流出的边界条件。

在计算过程中我遇到了一个问题，当压力结构到达出口边界时，会发生严重的扭曲。

这里主要有几个疑问：

1. 差分格式，库朗数，时间积分方式会是造成这个现象的原因吗？
2. 出口的速度边界采用了对流流出边界条件，出口的压力（以及构造压力泊松方程时出现的标量：速度散度）应该采用什么样的边界条件？
3. 我试着统计了一下压力泊松方程的残差，发现即使迭代100步，残差值 r = Ax-b 依然在0.1这个量级，这正常吗？

自己在读本科的时候，大概是对科研一点想法也没有的吧。从小到大被教育好好学习，上了大学就轻松了，于是我就真的大学轻松了四年，之后在找个工作，拿一份能够养活自己的薪水，大概这就是自己的人生了。然而事实却告诉我，并不能这样，2014年恰好遭遇行业严冬，因为学校牌子多少还算是个985，所以虽然自身条件不怎么样，我也顺利进入了一所国企，但是实际的薪资却是差强人意，面对高额的房价再估算一下自己工资未来的涨幅，无奈之下只能工作之余复习考研。

当时的自己特别想出国读书，看看外面的世界，然而遭到了家里的强烈反对，家长纷纷表示按我本科的表现想要出国读书实在是好高骛远，后来我又表示最少换个专业换个更好的学校，至少避免重蹈本科工作的覆辙，又被家里拒绝了。我个人性格又软（可能这就是我一直做不好的事情的原因），于是就听家里的话老实报回了本校本专业。当时或许也是觉得已经没有退路，反而激发了我的斗志，最后考出了400多分的成绩，让家人后悔不已，没有听我当初的话让我拼一把，不过当时我是觉得本校也有本校的好处，毕竟朋友都在这里生活上也比较开心。于是开开心心地区复试了。

复试前便是选导师，由于本科没有实验室经历，也就凭本科上课的感觉选了现在的导师。导师年纪不大，但是却非常有为，每年论文稳定产出，学术称号也是一堆，让我觉得我的选择是正确的。当时刚进实验室时，心里是想着虽然硕士没能出去，至少博士能够出去读。然而很可惜，事实却很骨感，导师初始对我很热情，然而得知我出国读博而不是继续跟他的想法之后，我立即变成了“别人的学生”，进入不管不问偶尔打杂的时间。无奈下我开始自己找课题自己弄，也就是这个时期了解到了OpenFOAM。当然最初只是单纯觉得这个做出来图“很好看”，外加许多编程的地方比较多，自己也多少有点编程的天赋，于是就入了这个坑，顺便还能再结合本专业的内容进行研究。因为周围也没有搞这个的，也就自己摸索，外加的确Openfoam在外网开源的资料比较多，也就对照的看一下。由于之前只写过一些面向过程的程序，记得刚看Openfoam这种大量运用面向对象特性的程序感到无从下手，不过毕竟只要愿意下功夫，还是会有收获的。自己也算磕磕绊绊地上手了，虽然自己没有涉及过CFD和FVM底层的东西，但最后也能自己根据自己需要写一些功能，修改一下求解器，也算是能够上手了。这个过程或许是非常痛苦的，然而当有所进步时，也同时是最开心的，自己也是在这个过程中第一次有了做科研很有趣的想法，后来自费参加了东岳老师的暑期培训，也逐渐思考着未来自己能不能向纯CFD算法靠拢。然而个人的能力毕竟是有限的，尤其是对于我们这种工科非理科，除非是算法上很大创新，否则光有数值模拟没有试验很难发出什么高质量的文章。结果满打满算，自己即将毕业时不过手握两篇EI水会，虽然凭借自己努力考出了比较高的语言成绩，然而这是博士并非硕士，语言成绩已经不是什么重要因素，自己的论文又实在不够看，在申请季收到的回复寥寥无几，要么是“实验室已满”，要么是“我没有钱如果你自费我们再考虑你”，于是自己现在还是前路未卜。

现在离毕业只有3个多月，我却手里还没有收到什么像样的offer，加上春招由于疫情缘故也转为网络形式，比起往年形式不容乐观。自己又是进又不得，退又不得，我现在26岁，家里人又一直催我找工作，早点安定下来。我却每天生活在痛苦当中,深感这三年的付出如同泡沫一般付诸东流，自己无论怎么折腾都只是在瞎蹦跶。我又该何去何从呢....真的就要这么放手吗....

可能我这个问题比较奇葩，目标是要在multiphaseInterFoam里引入造波功能，造波功能主要基于waves2Foam。我的思路是这样的：

首先我观察了一下intefoam跟multiphaseInterFoam对于mixture的求解区别。发现主要在于，interfoam是把两相分为主相以及背景相，求解时只求解主相，这样背景相可以简单的进行反向选取,即alpha2=1-alpha1就能获得了。

但是相数超过2，就不存在这个背景相的问题了，在multiphaseMixture里，每一相都单独计算，最后累加出一个总相alphas出来。

对于解速度压力那一套，两者我感觉没有太大区别。

我打算先从两相做起，在multiphaseInterFoam里先实现一般的水-气两相造波，再考虑第三相的问题。也就是打算用multiphaseInterFoam复现interfoam的功能。

我按照以下步骤进行：

1.修改multiphasemixture的代码，提供一个获取各个相的相分数的接口。
2.修改multiphaseInterFoam的代码，将这个接口提供给原先造波消波的松弛区，每一步都对波浪进行修正。
2.修改waves2Foam的代码，对空气设定一个相分数的入射边界（实际上就是波浪相分数边界的反向选取）。
3.在求解器中，每次循环都同时对alpha.water以及alpha.air在边界上进行更新。

比如对于水，边界上入射的相分数可能是(假设2D算例，竖直方向共5个网格)：
(0 0 0.6 1 1)
那么对于空气，就是
(1 1 0.4 0 0)
这样进行累加计算alphas的时候保证边界气液相分数累加为1，我在每一步的文件夹下观察结果，也的确没有发生越界。

但是实际运算的时候，出现了气液相分数越界，总相分数不守恒的现象，最后造出来的波是这种形状的

因为multiphaseInterFoam我想用在化工方面比较多一点，所以我想请问一下我是遗漏了哪些部分吗？

也请@东岳 老师指点。

第二个问题就是，我这个问题其实也类似一个分层流动，从上到下，第一相，第二相，第三相之间，第一相跟第三相是没有任何接触的。有没有划分成两个计算域，在第二相上交换信息的办法？

本人在学习OF之前，使用的商业软件基本都是依赖势流理论，引入格林函数基于边界元理论求解Laplace方程，最后用伯努利方程估算压力。

今天发现OF中也有一个potentialFoam是基于势流理论的，但是似乎跟我了解的势流理论有很大区别，只能计算稳态问题，但是控制方程的确也是Laplace方程。看起来两种求解器不单纯是FVM与边界元这种离散方式的区别，而是还有其他的不同，请问有了解这一块的大佬能讲讲这其中的问题所在吗？

@wwzhao 看起来显示的只是个并行计算的信息？
总而言之谢谢老师的指点了

@wwzhao 老师我想请教一个问题。我看到一些论文把重叠网格分为插值边界单元，贡献单元以及孤点单元，我想问着是否跟OpenFOAM的重叠单元的interpolated from local , mixed local/remote 以及 remoted相对应呢？

我计算的时候遇到了remoted单元不断上升的情况（初始0，计算100s后可能上升到几千），但是我的计算确实也没有发散，能够进行下去，如果remoted对应孤点单元，是不是代表我这个计算设定是不准确的呢？

@chaoscfd

想请问一下我是这是什么回事呢。。。
好像snap阶段结束后..就直接真的没了...
@东岳 顺便请教一下东岳老师

@chaoscfd
这里想请问您几个问题：

1. 建立refineRegions，我这里是用distance的办法定义的refineRegion，离模型xx米内用xx等级这样，我还需要另外用box对其局部定义吗？ 还是说我只需要提高distance的等级就可以了？
   2.第二就是includeAngle，我这里取得是90°，我之前尝试过取180°，的确新选取的边很多都是我想要的，但是也出现了好多我不想要的edge（比如明明就是个连续的面，结果突然出现一个edge斜着贯穿了这个面，类似一个正方型强行被分成2个三角形）。有没有什么办法，只捕捉到我想优化的edge，而抛弃一些莫名其妙捕捉到的edge呢

@东岳 好的 明白了 不过咱的求解器是基于InterFoam自己改的...我试试算的效果如何吧
谢谢老师