@Voynich 在 fvOptions(Urel)和SRF->Su()具体说明在哪儿？ 中说：

SRF->Su()是和坐标系相关的源项，就是你公式的最后两项，fvOptions(Urel)是自定义的源项。

谢谢，昨天研究了一下已经找到了
4737059c-7918-4a8a-b356-8a5771973def-image.png

我在使用subsetMesh命令时，出现以下提示，温度场无法正常设置，有哪位朋友遇到过吗
Subsetting volScalarField (T--> FOAM Warning :
From Foam::fixedValueFvPatchField<Type>::fixedValueFvPatchField(const Foam::fixedValueFvPatchField<Type> &, const Foam::fvPatch &, const Foam::DimensionedField<Type, Foam::volMesh> &, const Foam::fvPatchFieldMapper &) [with Type = double]
in file fields/fvPatchFields/basic/fixedValue/fixedValueFvPatchField.C at line 81
On field subsetT patch FIBER patchField fixedValue : mapper does not map all values.
To avoid this warning fully specify the mapping in derived patch fields.
alpha.liquid p p_rgh)

@random_ran 老师您好！借楼向您请教一些问题。我目前进行二维圆柱绕流的计算，Re=1e4和5e4。圆柱直径D=0.1m，介质是水(nu=1e-6)，计算域和生成网格（STARCCM生成并导入）如图所示：

图片1.png
图片3.png

湍流模型为 RAS kOmegaSST，边界条件参考Chalmers大学教程设置：

图片2.png

计算结果如下：最大库朗数为0.27，最大壁面y+<0.3，阻力系数误差<1%，St误差约为30%。St计算公式为：fD/U，f由升力时历曲线fft得到。

fvSchemes：

ddtSchemes { default Euler; } gradSchemes { default cellLimited Gauss linear 0.5; //grad(p) Gauss linear; grad(U) cellLimited Gauss linear 1; } divSchemes { default none; div(phi,U) Gauss linearUpwindV grad(U); div(phi,k) Gauss linearUpwind default; div(phi,omega) Gauss linearUpwind default; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear limited 1; } interpolationSchemes { default linear; } snGradSchemes { default limited 1; }

fvSolution

solvers { p { solver PCG; preconditioner DIC; tolerance 1e-08; relTol 0; } pFinal { solver PCG; preconditioner DIC; tolerance 1e-08; relTol 0; } U { solver BiCGStab; preconditioner DILU; tolerance 1e-07; relTol 0; } k { solver BiCGStab; preconditioner DILU; tolerance 1e-07; relTol 0; } epsilon { solver BiCGStab; preconditioner DILU; tolerance 1e-07; relTol 0; } omega { solver BiCGStab; preconditioner DILU; tolerance 1e-07; relTol 0; } } PISO { nCorrectors 3; nNonOrthogonalCorrectors 3; pRefCell 0; pRefValue 0; }

跪求老师指点迷津，如何提高St计算精度。

刚接触OpenFoam不就，疑问可能有点多
请问在pimpleFoam对松弛因子和* Final的松弛因子进行设置，如果设置了残差控制，在满足跳出外循环的残差条件之后是不是才会采用* Final设置的松弛因子进行计算？
如果达到设置的外循环次数后都没满足残差要求的话是不是最后一次循环还是用的原来的松弛因子计算？
我在其他地方有看到说最后一次需要采用一次松弛因子为1.0计算才能获得正确的结果，松弛因子设置为1.0是不是就相当于不进行松弛设置？
但是我在我自己写的一个求解器中，不在fvSolution中进行松弛设置和设置松弛因子为1.0得到的结果不一样，这又是为什么？

太酷炫了老铁，各种软件熟记于心

虚拟机上openfoam10和11没有问题

@李东岳 在 非定常CFL数相关问题 中说：

就是简单的一个传输的标量，附加一个s的传输方程

哦哦，谢谢老师！！

@bestucan 谢谢大佬

只有一个图。那就修网格吧。

最近编写了一个拉格朗日DPMFoam下面的颗粒长大的求解器，在OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/intermediate/lnInclude对源文件进行了修改。
源文件，求解器均编译成功，使用新求解器时候不能运行stochasticDispersionRAS的随机游走模型，有如下报错

Starting time loop Courant Number mean: 0 max: 0 Time = 0.01 Evolving kinematicCloud Solving 3-D cloud kinematicCloud #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib64/libc.so.6" #3 std::string::rfind(char, unsigned long) const in "/lib64/libstdc++.so.6" #4 Foam::IOobject::group() const at ??:? #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? #7 ? at ??:? #8 ? at ??:? #9 ? at ??:? #10 __libc_start_main in "/lib64/libc.so.6" #11 ? at ??:? Segmentation fault (core dumped）

修改了一下算法，算简单一点也有报错

Courant Number mean: 0 max: 0 Time = 0.001 Evolving kinematicCloud Solving 3-D cloud kinematicCloud [2] #[3] #0 0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[1] #0 Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:? at ??:? at ??:? at ??:? [3] #1 [1] #[2] #1 Foam::sigSegv::sigHandler(int)1 Foam::sigSegv::sigHandler(int)Foam::sigSegv::sigHandler(int)[0] #1 Foam::sigSegv::sigHandler(int) at ??:? [3] #2 ? at ??:? [2] #2 at ??:? [1] #2 ?? at ??:? [0] #2 ? in "/lib64/libc.so.6" [1] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [3] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [2] #3 std::string::rfind(char, unsigned long) const in "/lib64/libc.so.6" [0] #3 std::string::rfind(char, unsigned long) const in "/lib64/libstdc++.so.6" in "/lib64/libstdc++.so.6" [1] #4 Foam::IOobject::group() const[3] #4 Foam::IOobject::group() const in "/lib64/libstdc++.so.6" [2] #4 Foam::IOobject::group() const in "/lib64/libstdc++.so.6" [0] #4 Foam::IOobject::group() const at ??:? [1] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [3] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [2] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [0] #5 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::kModel() const at ??:? [1] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [3] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [2] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [0] #6 Foam::DispersionRASModel<Foam::KinematicCloud<Foam::Cloud<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > >::cacheFields(bool) at ??:? [1] #7 at ??:? [2] #7 at ??:? [3] #7 at ??:? [0] #7 ???? at ??:? [1] #8 at ??:? [2] #8 at ??:? [3] #8 at ??:? [0] #8 ???? at ??:? [1] #9 at ??:? [2] #9 at ??:? [3] #9 at ??:? [0] #9 ???? at ??:? [1] #10 __libc_start_main at ??:? [3] #10 __libc_start_main at ??:? [2] #10 __libc_start_main at ??:? [0] #10 __libc_start_main in "/lib64/libc.so.6" [1] #11 in "/lib64/libc.so.6" [3] #11 in "/lib64/libc.so.6" [2] #11 in "/lib64/libc.so.6" [0] #11 ???? at ??:? at ??:? at ??:? at ??:?

看起来是和Collisionparcel文件有关系？？？？？为什么我改一下颗粒直接变化还和碰撞模型的文件产生了影响？？我修改完了之后，不用随机游走模型倒是能跑，但是之后我是算例是需要运动随机游走模型的，所以我想搞想清楚原因，是不是我还要对CollisonParcel文件进行修改

@coolhhh 谢谢老师！

下面是我需要计算的方程：

fvScalarMatrix TbEqn ( - fvm::div(-phi, Tb) + fvm::laplacian(DT, Tb) ); TbEqn -= Q; TbEqn.solve();

其中，DT为热传导系数，Q为数值1的标量场，phi通过SIMPLE求解。
设置出口边界条件为：

outlet { type groovyBC; variables "Tk=DT;h=U&normal();Tinf=0;f=1/(1+Tk/(h*mag(delta())));"; valueExpression "Tinf"; gradientExpression "0"; fractionExpression "f"; value uniform 0; }

设置的divSchemes为：

div(-phi,Tb) Gauss linearUpwind grad(Tb);

使用的solvers为：

"(Tb)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-6; relTol 0; }

上述方程是用于拓扑优化求解的，但是优化过程中，经常出现浮点数错误“Signal: Floating point exception (8)”，如何有效解决这个问题呢？
尝试过在h*mag(delta()中添加一个小数，但是也同样出现浮点数报错。

@jlx 请问您在fe40中安装好wave2foam了没

@youhaoyu 我是使用compressibleVoF中的throttle算例更改了模型和边界条件

@sxz0823 您好，可以加个联系方式请教您一下turbinefoam吗，我想向您请教一些问题，我的qq是：935437658 :140:

大佬们，想请教下为什么OpenFOAM自定义并行计算时，初始化（这边初始化是指下述的代码）有时候很快2-3s，有时候需要很久1-2min（测试的时候没有其他程序运行）。
下述是我自定义的并行计算初始化代码：

scalarField MyField(5*N,0); List<scalarField> AllMyField_List(Pstream::nProcs()); AllMyField_List[Pstream::myProcNo()] = MyField; Pstream::gatherList(AllMyField_List); Pstream::scatterList(AllMyField_List); scalarField AllMyField(ListListOps::combine<scalarField>(AllMyField_List, accessOp<scalarField>()));

其中，N为所有网格单元的数量，N=86400。程序是有个变量，我将其分配到各个线程进行计算，在总线程中进行收集。
此外，程序中不止这一个变量，有五六个需要同样操作的变量。之前测试过，如果N很小的时候，初始化过程很快，但是当N增大时，初始化时间就完全不一样了。想请教下如何解决这个bug。

@chenboyao 天哪，膜拜大佬，成功解决了TAT

换课题对于学生来说是非常不友好的。有时候被换了也没招..

接触角不是在三相(固/液/气)的接触点才有吗? 两相流体的界面处应该没有接触角吧?

楼上说的对。版本差的太大，移植相对麻烦。