我想要对热管中的一个单元进行仿真，上壁面给定稳定的热源，下方为多孔介质区域，并假设初始状态液相含量为1。中间两侧为气槽，蒸发出来的气体通过气槽流入冷凝段。
现在我对固体壁面和气槽处的交界面（还有下面多孔区域和气槽的交界面）的温度边界条件不是很清楚应该怎么设置。稳态状态下气槽的温度近似3K的过热度，但是我这个是瞬态仿真，气槽的温度也是逐渐升高的。请问我应该怎么动态设置这个温度边界~有做过相关算例的前辈也可以分享一下您的边界设置~谢谢
1628062229(1).png

@李东岳 感谢李老师的快速回复。
不过我还是不太确定，尝试了以下语句：

label patchID = this->mesh_.boundaryMesh().findPatchID("A_to_B"); fvPatch& patch = this->mesh_.boundaryMesh()[patchID];

编译报错如下：

error: type 'Foam::fvPatch' does not provide a call operator refCast<const mappedPatchBase>(patch().patch());

该如何解决呢？网上查了半天没有类似的问题。

@冠竹 可以在ICEM导出.msh的时候，scale一下；也可以在转化网格的时候，用fluent3DMeshToFoam XXX.msh -scale 0.001缩放一下。

@tidedrinker 感谢分享！！！:146:

有没有做过水锤的老师，自己在流体力学书上找到了连续性方程和运动方程，在openfoam中也找到了可压缩求解器，但是看不懂求解器啥意思

@zx 我想计算管道内的水锤，因为管道有调压井，所以又有自由液面

看到tutorials里面只有层流的例子，这个不能计算紊流吗？

那三个文件中
phases (water air);

pMin 10000;

sigma 0.07;
这两个变量不知道什么意思
还有水的文件中
equationOfState
{
C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
}

这个也不知道是啥，我这个算例是不是不用变

@cresendo 太感谢了。我赶紧试试

@wying 我仔细看了一下，不是一个时间步长内，结构运动响应和流体域都求解多次，应该是一个时间步内，结构运动响应和流体域都求解一次。目前openfoam大部分求解器都是一个时间步内求解一次，只有少部分分为自时间步求解。

@hurricane007 谢谢老铁

@李东岳
https://www.cfd-online.com/Forums/openfoam-solving/175539-cyclic-patchtype-inletoutlet-boundary-condition.html
我是这个帖子看到有人提到他用inletoulet和cyclic 一起用的，不过他说并行有问题，
所以我试了一下串行的时候确实可以， 并行我倒是没试过。
我感觉两个物理量，一个循环，另一个就普通的进出口，在理论上也没啥问题吧？只是openfoam给了patchtype以后，有的边界条件就能在cyclic patch type上用了。

@李东岳 谢谢老师

@李东岳 网格步长0.4米，时间步长只有0.001，科朗数只有0.01，这是网格质量的原因么

谢谢老师。
我现在将k和omega的初始场的值设为0.1和10，计算就没有发散了，虽然计算的收敛速度较慢，但把松弛因子调小后也还是在不断减小。
fluent里面初始化k和omega场时也是按入口计算的值，为什么到OpenFoam里面这样初始化会导致收敛性变差呢？

@我是河滩 您好，个人刚接触openfoam不久，也需要设置一个随时间变化的边界条件，苦寻无果。如果您方便的话，能否给与一下详细的代码设置或者指点一番，不胜感激！！！

@李东岳 谢谢东岳教授

@李梦瑶 可以用cyclicAMI边界试试，或者用blockMeshDict画网格，pointwise画出的网格有时候用cyclic也行，icem的网格很难成功

@李东岳 好的谢谢老师

@supersoldier 你好，这个2.1是2.1.1版本吗？或者说2.1.1可以直接拿来用吗？

@tens 如果是在初始条件里写的：1.假设在alpha.water文件里写的，不确定是否可同时对saltwater和air初始化。2.writeEntry这句没写，不确定是否有影响

@oitocfd 好的，谢谢您的解答！