@zhou 谢谢，不过这个方法还是挺难搞得，真难搞这东西

@学流体的小明 在 周期性边界明渠的两点相关性分析Ruu、Rvv、Rww如何计算？ 中说：

取一个或多个平行于壁面的面（Kim的图也是不同高度的，y^+=5.39和149.23），使用surfaces工具输出这个面上的速度，然后每个面计算速度的展向、流向空间相关函数就行了。
可以直接计算二维的相关。
也可以固定z坐标，一条线一条线地计算流向的相关函数，再做个平均。

明白了，谢谢大佬

这个map我没使用过，具体我还真不了解这到底怎么回事

各位大佬好！
我现在做的是柔性体与流体的流固耦合问题，现在涉及到了网格无关性的验证问题
我想咨询下，我的算例安排为与实际算例相同的流固耦合问题好，还是说将柔性体交界面视为壁面，单纯进行流体计算好？
我的算例较为简单，只是流体与简单几何之间的耦合问题，柔性体变形不会特别离谱

/---------------------------------------------------------------------------\

\ / F ield OpenFOAM: The Open Source CFD Toolbox \ / O peration Website: https://openfoam.org \ / A nd Version: 6 \/ M anipulation

*---------------------------------------------------------------------------*/
Build : 6-fa1285188035
Exec : reactingFoam
Date : Dec 24 2024
Time : 10:45:39
Host : "DESKTOP-MGCJK2F"
PID : 3256
I/O : uncollated
Case : /home/mywsl/OpenFOAM/mywsl-6/run/turbulent
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Selecting chemistryReader foamChemistryReader
elements not defined in "/home/mywsl/OpenFOAM/mywsl-6/run/turbulent/constant/chem_Li_2007"
Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating reaction model

Selecting combustion model EDC
Selecting chemistry solver
{
solver ode;
method standard;
}

StandardChemistryModel: Number of species = 9 and reactions = 21
Selecting ODE solver seulex
using integrated reaction rate
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
deltaT = 1.111111111e-05
Time = 1.11111e-05

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 7.237527843e-08, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 6.496027629e-08, No Iterations 1
DILUPBiCGStab: Solving for H, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for H2, Initial residual = 0.9999999548, Final residual = 1.249524636e-08, No Iterations 1
DILUPBiCGStab: Solving for O, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for OH, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for H2O, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for O2, Initial residual = 0.9999997617, Final residual = 1.249524395e-08, No Iterations 1
DILUPBiCGStab: Solving for HO2, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for H2O2, Initial residual = 1, Final residual = 1.249524693e-08, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 0.001698464143, Final residual = 6.399191709e-07, No Iterations 1
min/max(T) = 300, 1800.356388
DICPCG: Solving for p, Initial residual = 0.8154141575, Final residual = 0.05901504271, No Iterations 4
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.1648808522, global = -0.04664958196, cumulative = -0.04664958196
DICPCG: Solving for p, Initial residual = 0.5937630302, Final residual = 9.1448632e-07, No Iterations 49
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.031207163e-06, global = 1.108666077e-07, cumulative = -0.0466494711
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.7898692099, Final residual = 7.356031497e-07, No Iterations 13
bounding epsilon, min: -57.68897069 max: 199.8018274 average: 2.853407683
DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 7.436499697e-07, No Iterations 14
bounding k, min: -53.76693743 max: 137.3801042 average: 1.662759794
ExecutionTime = 7.11 s ClockTime = 7 s

Courant Number mean: 0.3263677773 max: 38.86161238
deltaT = 5.71632726e-08
Time = 1.11683e-05

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.01088173554, Final residual = 9.291750018e-07, No Iterations 2
DILUPBiCGStab: Solving for Uy, Initial residual = 0.126974428, Final residual = 5.413457541e-07, No Iterations 3
DILUPBiCGStab: Solving for H, Initial residual = 0.05031914993, Final residual = 2.526806224e-07, No Iterations 3
DILUPBiCGStab: Solving for H2, Initial residual = 0.04988384318, Final residual = 2.536567037e-07, No Iterations 3
DILUPBiCGStab: Solving for O, Initial residual = 0.05023601744, Final residual = 2.528851275e-07, No Iterations 3
DILUPBiCGStab: Solving for OH, Initial residual = 0.0498739039, Final residual = 2.53709804e-07, No Iterations 3
DILUPBiCGStab: Solving for H2O, Initial residual = 0.04983585131, Final residual = 2.537599816e-07, No Iterations 3
DILUPBiCGStab: Solving for O2, Initial residual = 0.04987406878, Final residual = 2.536836372e-07, No Iterations 3
DILUPBiCGStab: Solving for HO2, Initial residual = 0.04887920295, Final residual = 2.557918423e-07, No Iterations 3
DILUPBiCGStab: Solving for H2O2, Initial residual = 0.04849527326, Final residual = 2.565684649e-07, No Iterations 3
DILUPBiCGStab: Solving for h, Initial residual = 0.05497148043, Final residual = 4.9437153e-07, No Iterations 3
min/max(T) = 299.9248087, 1796.462561
DICPCG: Solving for p, Initial residual = 0.01518365283, Final residual = 3.016271454e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.870692989e-09, global = -5.333936149e-09, cumulative = -0.04664947643
DICPCG: Solving for p, Initial residual = 0.0003488105302, Final residual = 1.108179983e-09, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.557155586e-10, global = 8.320508309e-11, cumulative = -0.04664947635
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.1451083224, Final residual = 2.148726213e-07, No Iterations 2
DILUPBiCGStab: Solving for k, Initial residual = 0.8290663092, Final residual = 5.98405794e-08, No Iterations 4
ExecutionTime = 8.97 s ClockTime = 9 s

Courant Number mean: 0.0013035774 max: 0.1228067562
deltaT = 6.644108124e-08
Time = 1.12347e-05

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.1832008417, Final residual = 3.515359245e-07, No Iterations 8
DILUPBiCGStab: Solving for Uy, Initial residual = 0.8416763107, Final residual = 7.756203776e-07, No Iterations 10
DILUPBiCGStab: Solving for H, Initial residual = 0.4428609204, Final residual = 9.062138129e-07, No Iterations 9
DILUPBiCGStab: Solving for H2, Initial residual = 0.4410214426, Final residual = 9.736737868e-07, No Iterations 9
DILUPBiCGStab: Solving for O, Initial residual = 0.4424908421, Final residual = 9.30398611e-07, No Iterations 9
DILUPBiCGStab: Solving for OH, Initial residual = 0.4409439717, Final residual = 9.745009796e-07, No Iterations 9
DILUPBiCGStab: Solving for H2O, Initial residual = 0.4408222366, Final residual = 9.738422502e-07, No Iterations 9
DILUPBiCGStab: Solving for O2, Initial residual = 0.4409741964, Final residual = 9.739327596e-07, No Iterations 9
DILUPBiCGStab: Solving for HO2, Initial residual = 0.4368143603, Final residual = 7.794393661e-07, No Iterations 9
DILUPBiCGStab: Solving for H2O2, Initial residual = 0.4352225269, Final residual = 9.850169941e-07, No Iterations 9
DILUPBiCGStab: Solving for h, Initial residual = 0.4603404275, Final residual = 8.710686777e-07, No Iterations 10
min/max(T) = 299.9999982, 1791.283236
DICPCG: Solving for p, Initial residual = 0.04989128258, Final residual = 1.728078903e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.9720677e-09, global = -2.805246459e-09, cumulative = -0.04664947915
DICPCG: Solving for p, Initial residual = 0.00040507004, Final residual = 2.635572617e-10, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 6.291142477e-11, global = 5.449968679e-11, cumulative = -0.0466494791
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.2709179929, Final residual = 2.790707667e-07, No Iterations 9
DILUPBiCGStab: Solving for k, Initial residual = 0.6253639963, Final residual = 3.507848237e-07, No Iterations 10
ExecutionTime = 10.81 s ClockTime = 11 s

Courant Number mean: 0.0008942027266 max: 0.05529325678
deltaT = 7.968158396e-08
Time = 1.13144e-05

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.06593954665, Final residual = 3.351747253e-07, No Iterations 8
DILUPBiCGStab: Solving for Uy, Initial residual = 0.4102152305, Final residual = 9.318726704e-07, No Iterations 10
DILUPBiCGStab: Solving for H, Initial residual = 0.2892697203, Final residual = 5.28621965e-07, No Iterations 9
DILUPBiCGStab: Solving for H2, Initial residual = 0.2624968155, Final residual = 4.376897291e-07, No Iterations 9
DILUPBiCGStab: Solving for O, Initial residual = 0.2808888368, Final residual = 5.574822147e-07, No Iterations 9
DILUPBiCGStab: Solving for OH, Initial residual = 0.2591248329, Final residual = 4.295608139e-07, No Iterations 9
DILUPBiCGStab: Solving for H2O, Initial residual = 0.2597666461, Final residual = 4.522225416e-07, No Iterations 9
DILUPBiCGStab: Solving for O2, Initial residual = 0.2612620687, Final residual = 4.395858975e-07, No Iterations 9
DILUPBiCGStab: Solving for HO2, Initial residual = 0.2220485989, Final residual = 5.99379148e-07, No Iterations 9
DILUPBiCGStab: Solving for H2O2, Initial residual = 0.1820866327, Final residual = 6.019580413e-07, No Iterations 10
DILUPBiCGStab: Solving for h, Initial residual = 0.09990726406, Final residual = 5.829150064e-07, No Iterations 10
min/max(T) = 299.9999961, 1784.763549
DICPCG: Solving for p, Initial residual = 0.02262674874, Final residual = 3.398546862e-08, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.084807953e-09, global = -4.286910994e-09, cumulative = -0.04664948339
DICPCG: Solving for p, Initial residual = 0.0001704309813, Final residual = 3.621876029e-10, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.790648795e-11, global = 1.74305778e-11, cumulative = -0.04664948337
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.0897510291, Final residual = 6.434387951e-07, No Iterations 3
DILUPBiCGStab: Solving for k, Initial residual = 0.2857450205, Final residual = 8.582666035e-07, No Iterations 9
ExecutionTime = 12.45 s ClockTime = 12 s

Courant Number mean: 0.001010967894 max: 0.05839075857
deltaT = 9.556638249e-08
Time = 1.141e-05

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 0.01353903567, Final residual = 4.237016717e-07, No Iterations 6
DILUPBiCGStab: Solving for Uy, Initial residual = 0.1965833104, Final residual = 4.895900842e-07, No Iterations 7

统一回复一下：

感谢给的热心解答。

用线提取密度的垂向变化(内波)，表面波的probe暂时不适用这个问题。

有效解决方案：
方案A：用toposet选取一定区域的网格，输出rho，这就有了xyzrho的数据。然后根据需要的xz坐标，做插值即可。
方案B：选取全局的网格，输出cuttingPlane上的数据。
方案C：输出特定等直面的数据。

方案A灵活，但操作复杂；方案B数据最全，但存储空间比较大；方案C直接，对三维工况不友好。

这都是后处理的技巧，能获得需要的数据。

仿照https://www.cfd-online.com/Forums/openfoam-solving/60471-call-fortran-program-openfoam.html进行Fortran程序的调用，编译过程中报错： from chemistryModel/rho2ChemistryModel/rho2ChemistryModels.C:37:
/home/sy/OpenFOAM/sy-v1706/src/thermophysicalModels/strath/strathSpecie/lnInclude/CEATransportI.H:152:7: error: expected constructor, destructor, or type conversion before ‘(’ token
phypro(&TG, &TE, &PT, &PE, &CON[0], &CP_FINAL[0][0], &H_FINAL[0][0], &MIU_FINAL, &K_FINAL[0]);
但仿照https://www.cnblogs.com/snake553/p/6962386.html编写C++程序可以实现Fortran程序的调用

只修改了粘度出现一下报错想问问该如何解决，粘度为2e-6m2/s
[8] #0 Foam::error::printStack(Foam::Ostream&)[11] #0 [14] #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) at ??:?
[11] #1 Foam::sigFpe::sigHandler(int) at ??:?
[8] #1 Foam::sigFpe::sigHandler(int) at ??:?
[14] #1 Foam::sigFpe::sigHandler(int) at ??:?
[11] #2 ? at ??:?
[8] #2 ? at ??:?
[14] #2 ? in "/lib64/libpthread.so.0"
[8] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libpthread.so.0"
[11] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libpthread.so.0"
[14] #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
at ??:?
[8] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&)[11] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[14] #4 void Foam::multiply<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[8] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[11] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[14] #5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[8] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[11] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[14] #6 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() at ??:?
[8] #7 at ??:?
[11] #7 at ??:?
[14] #7 ??? at ??:?
[8] #8 __libc_start_main at ??:?
[11] #8 __libc_start_main at ??:?
[14] #8 __libc_start_main in "/lib64/libc.so.6"
[8] #9 in "/lib64/libc.so.6"
[11] #9 in "/lib64/libc.so.6"
[14] #9 ??? at ??:?
[ibnode77:65243:0:65243] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
at ??:?
[ibnode77:65246:0:65246] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
at ??:?
[ibnode77:65249:0:65249] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
==== backtrace (tid: 65243) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0 [ibnode77:65243] *** Process received signal ***
[ibnode77:65243] Signal: Floating point exception (8)
[ibnode77:65243] Signal code: (-6)
[ibnode77:65243] Failing at address: 0x4080000fedb
[ibnode77:65243] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b424c49c630]
[ibnode77:65243] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b424c49c4fb]
[ibnode77:65243] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b424c49c630]
[ibnode77:65243] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b424b51044d]
[ibnode77:65243] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b4247910418]
[ibnode77:65243] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b424799b8b6]
[ibnode77:65243] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b42479e7676]
[ibnode77:65243] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65243] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b424c6cb555]
[ibnode77:65243] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65243] *** End of error message ***
==== backtrace (tid: 65246) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0 [ibnode77:65246] *** Process received signal ***
[ibnode77:65246] Signal: Floating point exception (8)
[ibnode77:65246] Signal code: (-6)
[ibnode77:65246] Failing at address: 0x4080000fede
[ibnode77:65246] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b3d893ab630]
[ibnode77:65246] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b3d893ab4fb]
[ibnode77:65246] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b3d893ab630]
[ibnode77:65246] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b3d8841f44d]
[ibnode77:65246] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b3d8481f418]
[ibnode77:65246] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b3d848aa8b6]
[ibnode77:65246] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b3d848f6676]
[ibnode77:65246] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65246] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b3d895da555]
[ibnode77:65246] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65246] *** End of error message ***
==== backtrace (tid: 65249) ====
0 0x000000000000f4fb raise() ???:0
1 0x000000000069a44d Foam::multiply() ???:0
2 0x00000000000d4418 Foam::multiply<Foam::fvPatchField, Foam::volMesh>() ???:0
3 0x000000000015f8b6 Foam::operator*<Foam::fvPatchField, Foam::volMesh>() ???:0
4 0x00000000001ab676 Foam::RASModels::RNGkEpsilon<Foam::IncompressibleTurbulenceModelFoam::transportModel >::correct() ???:0
5 0x000000000042a05e main() ???:0
6 0x0000000000022555 __libc_start_main() ???:0
7 0x000000000042e172 _start() ???:0 [ibnode77:65249] *** Process received signal ***
[ibnode77:65249] Signal: Floating point exception (8)
[ibnode77:65249] Signal code: (-6)
[ibnode77:65249] Failing at address: 0x4080000fee1
[ibnode77:65249] [ 0] /lib64/libpthread.so.0(+0xf630)[0x2b14b9db5630]
[ibnode77:65249] [ 1] /lib64/libpthread.so.0(raise+0x2b)[0x2b14b9db54fb]
[ibnode77:65249] [ 2] /lib64/libpthread.so.0(+0xf630)[0x2b14b9db5630]
[ibnode77:65249] [ 3] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(ZN4Foam8multiplyERNS_5FieldIdEERKNS_5UListIdEES6+0x11d)[0x2b14b8e2944d]
[ibnode77:65249] [ 4] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4Foam8multiplyINS_12fvPatchFieldENS_7volMeshEEEvRNS_14GeometricFieldIdT_T0_EERKS6_S9+0xa8)[0x2b14b5229418]
[ibnode77:65249] [ 5] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(ZN4FoammlINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS7_RKS8+0x156)[0x2b14b52b48b6]
[ibnode77:65249] [ 6] /share/home/wangshuai/OpenFOAM/OpenFOAM-5.x/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels11RNGkEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x386)[0x2b14b5300676]
[ibnode77:65249] [ 7] cfdemSolverPiso[0x42a05e]
[ibnode77:65249] [ 8] /lib64/libc.so.6(__libc_start_main+0xf5)[0x2b14b9fe4555]
[ibnode77:65249] [ 9] cfdemSolverPiso[0x42e172]
[ibnode77:65249] *** End of error message *** Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted. mpirun noticed that process rank 14 with PID 0 on node ibnode77 exited on signal 8 (Floating point exception).

已解决，解决方案是scalarDict.set("vv", "vv [0 0 0 0 0 0 0] -0.1");

@Shihang-Chen
你好，我在做类似的代码。想请教伪时间步植入后，计算结果中出现松弛要怎么处理（计算的周期和实际周期不一致）。我是在动量方程中加入的隐式伪时间步。
512ceafd-f193-4d59-a045-c2859969922c-9e6a4889e6591786ea13032244e9111.png

@学流体的小明 后来我看了一下运行日志发现量纲不匹配并不是动量方程的量纲不匹配，是我在定义阻尼区起点的时候没有定义量纲所以会一直报错，后来把阻尼区起点定义量纲之后就可以啦，整个代码是这样的。```

dimensionedScalar r0 ( "r0", dimensionSet(0,0,-1,0,0,0,0), 0 ); dimensionedScalar r1 ( "r1", dimensionSet(0,0,-1,0,0,0,0), 5 ); dimensionedScalar x1 = max(mesh.Cf().component(vector::X))-min(mesh.Cf().component(vector::X)); dimensionedScalar x0 ( "x0", dimensionSet(0,1,0,0,0,0,0), 12 ); volScalarField R ( "R", max(r1*(mesh.C().component(vector::X)-x0)/(x1-x0),r0) ); fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U) - fvm::laplacian(muf,U) - (fvc::grad(U)&fvc::grad(muf)) + rho*U*R );

@李东岳 @欧阳 我计算ODW，一百万网格，氢氨机理s59r356，cpuLoad（redistributionInterval 20; maxImbalance 0.1;）和DLBFoam（只使用负载均衡，不使用pyJac）效果基本相当。

各位大佬，我不小心在Ubunto家目录下执行了foamCleanTutorials，命令行报

mxy@LAPTOP-CODMI2H5:~$ foamCleanTutorials Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/DNS/dnsFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/acoustic/acousticFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/basic/laplacianFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/basic/potentialFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/basic/scalarTransportFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/PDRFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/XiFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/chemFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/coldEngineFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/fireFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/combustion/reactingFoam wclean libso BCs Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/compressible/rhoPimpleAdiabaticFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/compressible/rhoPimpleFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/compressible/rhoSimpleFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/compressible/sonicFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/discreteMethods/dsmc/dsmcFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/discreteMethods/molecularDynamics/mdEquilibrationFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/discreteMethods/molecularDynamics/mdFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/electromagnetics/electrostaticFoam Cleaning application /home/mxy/OpenFOAM/OpenFOAM-v2312/applications/solvers/electromagnetics/magneticFoam

删了我好多，这种情况需要重新wmake求解器吗？？

963df654ce25372fc5621513b9fe905.png 如图是blockMesh中创建的oversetPatch，用于触发在任何其他边界条件之前进行重叠（overset）插值。ccmTo.png 这张图的polyMesh是利用ccmToFoam，将starccm网格导入到OpenFoam的，请问如何在starccm里面加上图1中的oversetPatch呢。@李东岳

是不是边界的名重复了呢

算出来的最低温度一直在降

数值扩散掉了

是的，你试试

@1064168551 在 求助：在codedict中读取求解器中定义的标量 中说：

@李东岳 a6bccb55-b62e-4688-bc61-dc2f9e6911e5-image.png 东岳老师，抱歉再打扰您，我现在改成这个样子，还是会报错，dimensionedScalar好像还是无法lookup

报错在这里，就还是这行const dimensionedScalar vv = mesh.lookupObject<dimensionedScalar>("vv");有问题
56be09a2-6412-49a3-9c59-598dfb8ef235-image.png

遇到了类似的问题，就是水平集模拟两相流时，怎么样利用水平集计算的相界面把一个对流扩散方程隔离开，只在一个相内求解对流扩散方程。

发散的话就是发散了，问题就很多，需要你详细的debug