是的，都是绝对坐标系的

最近在家呆的闲，又没有实验室的资源，，想看看磁流体的流固耦合的代码和书最好中文的，，，求推荐啊～

滑移壁面没有剪切力

@bubbleEL 您好，请问您德高斯分布函数已经实现了吗？最近也在实现这种办法，不知道能否指点一下？另外一直在思考超出网格之后发散的数学原理是什么？不知道您这边是不是已经弄清楚了？

@东岳 是的，还和网格有关，均匀网格和非均匀网格产生的结果也不同。

@程迪 压力都是零梯度边界条件，压力泊松方程的系数矩阵是奇异的，adjustphi让矩阵不奇异了吗？奇异矩阵怎么解的？

可以用吗？如何解决，求分享

@东岳 明白了！谢谢老师

@Exthan 不知道是边界条件设错了还是应该更换求解器

不是的，后面我用komegaSST算了一下这个模型，吧压力求解器改成了PCG之后就收敛了，而且结果比较满意，希望对遇到同样问题的朋友也有参考价值

大家好，我最近在模拟波浪。在现在interfoam底下的case里，对于波的设置条件里，有一项ramptime，我不清楚它的意思，请大家帮忙解答。谢谢🙏

入口边界条件？我有一个问题，最近看到一些模拟射流的把入口延长了一段，加了管内流动，用来避免喷管直接出口处的较大的梯度，东岳老师怎么看@东岳 ，不知道您的入口速度怎么给的

@东岳 我是按这个公式计算的，密度为1，那么u_tau=sqrt(tau_w)，但是我发现，论文里的摩擦速度刚好是我用OF算出来的1.225倍。论文是用fluent算的，fluent里面的空气密度是1.225。我现在想不通 为什么会差这1.225倍。:jingya:

@东岳 感谢🙏🏼 试了一下还是不行，有新的其他函数的错误了😹

回复: 请教个关于编译求解器时报错的问题solidificationFoam.zip :high:

@马乔 在 OF文件读写 中说：

#include "fvCFD.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // // class proxyWriter : public regIOobject { private: scalar value_; proxyWriter(const proxyWriter&) = delete; void operator=(const proxyWriter&) = delete; public: proxyWriter(const IOobject& io, const scalar& value) : regIOobject(io), value_(value) { } virtual bool writeData(Ostream& os) const { os<<value_; return os.good(); } proxyWriter& operator++() { ++value_; return *this; } }; int main(int argc, char *argv[]) { #include "setRootCaseLists.H" #include "createTime.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // proxyWriter writer ( IOobject ( "value", runTime.timeName(), runTime, IOobject::NO_READ, IOobject::AUTO_WRITE ), 0. ); Info<< "\nStarting time loop\n" << endl; while (runTime.loop()) { Info<< "Time = " << runTime.timeName() << nl << endl; ++writer; runTime.write(); runTime.printExecutionTime(Info); } Info<< "End\n" << endl; return 0; } :chigua:

@马乔 @Samuel-Tu 请问这个怎么用？这是重新编译一个lib然后link到controlDict吗？那怎么知道是在这一步的什么时候输出？是计算前还是计算后？

大家好

f29cfccc-509a-4c69-bbd0-3eddb612ac78-图片.png

如图所示的二维空间，蓝色圆形的导热系数与白色区域不同，我可以用blockMesh生成规则网格，利用funkySetFieldsDict在已知圆柱中心、半径情况下定义导热系数。

对于三维问题，我利用3dmax生成了不同直径的堆积球，球心、半径的信息我不清楚，如图所示生成stl文件：

d750e2b8-86fd-41c9-aea4-216a88ba6a85-图片.png

接下来我想计算整体的导热问题，这也是一个各向异性导热问题，我仅需要用blockmesh画网格即可，但如何定义球的物性？

大家好，我在字典文件中看到

xcells #calc "round((($xmax)-($xmin))/$delta)";

想知道#calc 可以使用的标准函数包括哪些呢？哪里能查得到？

最近在研究边界条件的问题，想复制一个边界条件出来修改编译测试。参考了网上的教程，但是还是有些问题，想请教大家。

OpenFOAM里面的边界条件有些只包括一个C文件，一个H文件，这种边界条件按照网上的教程都可以直接编译。比如fixedFluxPressure，我拷贝到新的文件夹test下，创建一个Make文件夹，在其内部创建files和options文件：
files:

fixedFluxPressure/fixedFluxPressureFvPatchScalarField.C LIB = $(FOAM_USER_LIBBIN)/libUsr

options:

EXE_INC = \ -I$(LIB_SRC)/finiteVolume/lnInclude\ LIB_LIBS = \ -lfiniteVolume

编译正常。
但是还有一部分的边界条件包含多个C文件和H文件，这种就编译不成功了。以cyclic为例，拷贝$FOAM_SRC/finiteVolume/fields/fvPatchFields/constraint/cyclic 到testcyclic文件夹，test下面生成Make文件夹，Make里面建立files和options:

files:

cyclic/cyclicFvPatchField.C LIB = $(FOAM_USER_LIBBIN)/libcyclicFvPatchField

options:

EXE_INC = \ -I$(LIB_SRC)/finiteVolume/lnInclude\ LIB_LIBS = \ -lfiniteVolume

执行wmake libso

wmake libso . wmakeLnInclude: linking include files to ./lnInclude Making dependency list for source file cyclicFvPatchField.C g++ -std=c++11 -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -Wno-attributes -O3 -DNoRepository -ftemplate-depth-100 -I/home/xc15/OpenFOAM/OpenFOAM-7/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/xc15/OpenFOAM/OpenFOAM-7/src/OpenFOAM/lnInclude -I/home/xc15/OpenFOAM/OpenFOAM-7/src/OSspecific/POSIX/lnInclude -fPIC -c cyclic/cyclicFvPatchField.C -o Make/linux64GccDPInt32Opt/cyclic/cyclicFvPatchField.o cyclic/cyclicFvPatchField.C:32:1: error: redefinition of ‘Foam::cyclicFvPatchField<Type>::cyclicFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&)’ Foam::cyclicFvPatchField<Type>::cyclicFvPatchField ^ In file included from cyclic/cyclicFvPatchField.H:225:0, from cyclic/cyclicFvPatchField.C:26: cyclic/cyclicFvPatchField.C:32:1: error: ‘Foam::cyclicFvPatchField<Type>::cyclicFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&)’ previously declared here Foam::cyclicFvPatchField<Type>::cyclicFvPatchField ^ cyclic/cyclicFvPatchField.C:44:1: error: redefinition of ‘Foam::cyclicFvPatchField<Type>::cyclicFvPatchField(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const Foam::dictionary&)’ Foam::cyclicFvPatchField<Type>::cyclicFvPatchField ........ make: *** [Make/linux64GccDPInt32Opt/cyclic/cyclicFvPatchField.o] Error 1

error太长，只贴了一部分，但是基本上都是说重复定义和声明

正确的编译方式应该是什么？:mihu:

我想问一下，有没有人在openFoam中植入DBM（Discrete Bubble Model）呢。感谢。