CFD中文网

@李东岳 岳哥，别笑话我，求教我咋搞他

各位大佬，现想利用fe中的block couple方法耦合求解，需要耦合的项是laplacian项，但是现在fe5.0中没有该方法，于是我自己添加。一顿操作后，编译遇到如下问题，始终没能定位问题出在哪，求大佬指导一二。

In file included from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/PtrList.H:315,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/List.C:30,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/List.H:291,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/labelList.H:36,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/Pstream.H:43,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/OPstream.H:35,
from /home/dong/foam/foam-extend-5.0/src/foam/lnInclude/parRun.H:35,
from /home/dong/foam/foam-extend-5.0/src/finiteVolume/lnInclude/fvCFD.H:4,
from pbFourEqnsTwoPhaseFoam.C:43:
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/PtrList.C: In instantiation of ‘Foam::PtrList<T>::PtrList(const Foam::PtrList<T>&) [with T = Foam::CoeffField<double>]’:
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/FieldField.C28: required from ‘Foam::FieldField<Field, Type>::FieldField(const Foam::FieldField<Field, Type>&) [with Field = Foam::CoeffField; Type = double]’
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/BlockLduMatrix.C:84:28: required from ‘Foam::BlockLduMatrix<Type>::BlockLduMatrix(const Foam::BlockLduMatrix<Type>&) [with Type = double]’
/home/dong/foam/foam-extend-5.0/src/finiteVolume/lnInclude/BlockLduSystem.C:99:24: required from ‘Foam::BlockLduSystem<blockType, sourceType>::BlockLduSystem(const Foam::BlockLduSystem<blockType, sourceType>&) [with blockType = double; sourceType = double]’
pTEqn.H:73:70: required from here
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/PtrList.C:54:18: error: invalid conversion from ‘Foam::Field<double>’ to ‘Foam::CoeffField<double>’ [-fpermissive]
54 | ptrs_[i] = (a[i]).clone().ptr();

@pengdi 多相密度基方法可参考Saurel的文章：doi:10.1017/S0022112008002061，基于迎风格式应对多可压缩流动问题处理效果较好，求解器是自己开发的。

动量插值里引入的压力梯度三阶导项前面的系数乘了deltaT，如果deltaT太小压力光顺效果会变差。

用CompressibleIsoInterFoam算超声速射流有两个本质上无法解决的问题：

1. 压力基Rhie-Chow插值属于中心插值，可压缩出现shock时会振荡，如果开启了choi修正加入了耗散会把激波抹平（当然看你关注什么，如果只是关注推力之类的统计量可能没啥影响，但是研究流场结构就够呛了）
2. 如果你认真推一遍公式，会发现它的能量（温度）方程压根就是错的，算出负温度不足为奇。另外，我记得thermo里的mixture函数在每个时间步计算时都会initial一下，温度会重新算，可以按“Negative initial”关键字搜索定位一下代码。
   总之吧，这个求解器勉强能用，如果你换套网格、换个格式（比如snGrad用limited）或许到这个时刻就能算通了，然而隐患一直都在，不知啥时候又会出现。
   如果真想把浸没射流算好，建议用多相的density-based solver，我们算的效果长这样：

@yhdthu 大佬别闹，求赐教

@heike256 东岳哥好，我是把两相密度基算法(p-relax, 5方程)用压力基的求解套路实现一下，受制于Rhie-Chow插值，处理压缩性较强的问题不是很给力（会振荡），所以想找个相对简单点的算例测试一下，暂时没找到

如题，最近写了个新的求解器，想验证一下算法，有哪位大佬指教一下

@李东岳 哈哈哈，物理隔绝了

@yfleung jhd上的代码有bug，且效率低，我自己重写可用。
https://cloud.tsinghua.edu.cn/d/05b0149e37904e7ca61c/

@izumi 我理解是这样：SRF给边界是相对速度，所以跟着坐标系一起转。MRF是绝对速度，所以固定不动。

@calrel https://www.cfd-online.com/Forums/openfoam-pre-processing/62140-mapfields-only-internal-fields.html

我看到DFSEM入口有bug？

讨论地址：https://www.cfd-online.com/Forums/openfoam-solving/177711-turbulentdfseminlet.html

有问题的是11式，本来是个无量纲参数，根号出来其实有量纲，不知有人测试过没？

模型还是很简单的，就是一个柱体+方形，但是我用icem画必须要用两个O grid，网格过渡很不好（至少我看起来很难受），不知有什么其他思路？

由于对网格质量和控制的要求，我把流体域分成了两部分进行网格划分，显然二者有个重合的界面。
我先用blockMesh生成了核心区的网格。之后我想用snappyhexmesh生成非核心区网格。
我想请问，是否可以将已经有了的核心区的界面导入，相当于固定了面网格，去生成非核心区网格呢？

我目前没查到相关方法，我试了mergeMesh，但效果很不好，不知各位有何高见？

所以最终结论是PCG+DIC效率最高是么？
PCG+GAMG是不是稳定性好一些？

如题。

@xpqiu 好的，谢谢邱老师，我再试试😄

@xpqiu 好的，我的网格在观测点处太粗了，时间步长选择的和探针的采样频率一致，是不是需要比这个值还小呢？

@xpqiu 是的，定性的可以对上。计算中波动被磨平了很多，您认为是怎么回事呢？