[EFC]

Hello everyone.

After a few seconds my simulation crash with this error message.

------------------------------------------------------------------
Solving for fluid region SOLID
DILUPBiCGStab: Solving for Ux, Initial residual = 0.0195625, Final residual = 1.26188e-05, No Iterations 2
DILUPBiCGStab: Solving for Uy, Initial residual = 0.0188434, Final residual = 9.64314e-06, No Iterations 2
DILUPBiCGStab: Solving for Uz, Initial residual = 0.0150732, Final residual = 5.20369e-06, No Iterations 2
wallCondensationSource: applying source 2 to H2O for fieldi: 1
DILUPBiCGStab: Solving for H2O, Initial residual = 0.0108213, Final residual = 3.22714e-09, No Iterations 3
SOLID:AIRE3:T <- region1:AIRE3:T :
heat transfer rate:13694.1
Total mass flux [Kg/s]: 0.0044304
Total mass on the wall [Kg]: 0.0278459
Total heat (>0 leaving the wall to the fluid) [W]: 10339.2
Total heat (>0 leaving the wall to the solid) [W]: 52871.2
Total latent heat released [W]: 10571.3
Total specific heat removed from fluid [W]: 2686.46
walltemperature min:321.627 max:328.747 avg:327.402
SOLID:SSW_EXT:T <- region2:SSW_EXT:T :
heat transfer rate:8115.43
Total mass flux [Kg/s]: 0.0212632
Total mass on the wall [Kg]: 0.0127632
Total heat (>0 leaving the wall to the fluid) [W]: -1404.5
Total heat (>0 leaving the wall to the solid) [W]: 28051.1
Total latent heat released [W]: 50766.5
Total specific heat removed from fluid [W]: 12879.4
walltemperature min:321.12 max:321.26 avg:321.182
SOLID:AIRE2:T <- region1:AIRE2:T :
heat transfer rate:3212.42
Total mass flux [Kg/s]: 0.00211925
Total mass on the wall [Kg]: 0.00183134
Total heat (>0 leaving the wall to the fluid) [W]: 1325.41
Total heat (>0 leaving the wall to the solid) [W]: 2125.84
Total latent heat released [W]: 5049.35
Total specific heat removed from fluid [W]: 1288.37
walltemperature min:322.309 max:358.914 avg:328.309
SOLID:AIRE1:T <- region1:AIRE1:T :
heat transfer rate:179.837
Total mass flux [Kg/s]: 0
Total mass on the wall [Kg]: 0.000128344
Total heat (>0 leaving the wall to the fluid) [W]: 182.747
Total heat (>0 leaving the wall to the solid) [W]: 161.861
Total latent heat released [W]: 0
Total specific heat removed from fluid [W]: 0
walltemperature min:332.974 max:333.719 avg:333.441
wallCondensationSource: applying source 2 to h for fieldi: 0
DILUPBiCGStab: Solving for h, Initial residual = 0.0109704, Final residual = 1.95067e-07, No Iterations 2
[3] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::co nstTransport<Foam::species::thermo<Foam::hConstThe rmo<Foam:erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::calculate() at ??:?
[3] #4 Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::co nstTransport<Foam::species::thermo<Foam::hConstThe rmo<Foam:erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::correct() at ??:?
[3] #5 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
[3] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #7 ?[cuba:03305] *** Process received signal ***
in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/chtMultiRegionFoam"
[cuba:03305] Signal: Floating point exception (8)
[cuba:03305] Signal code: (-6)
[cuba:03305] Failing at address: 0x3e800000ce9
[cuba:03305] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3f040)[0x7f40a9c3f040]
[cuba:03305] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xc7)[0x7f40a9c3efb7]
[cuba:03305] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x3f040)[0x7f40a9c3f040]
[cuba:03305] [ 3] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so(_ZN4Foam11heRho ThermoINS_17rhoReactionThermoENS_13SpecieMixtureIN S_15reactingMixtureINS_14constTransportINS_7specie s6thermoINS_12hConstThermoINS_10perfectGasINS_6spe cieEEEEENS_16sensibleEnthalpyEEEEEEEEEE9calculateE v+0x464)[0x7f40b16d8064]
[cuba:03305] [ 4] /opt/openfoam6/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so(_ZN4Foam11heRho ThermoINS_17rhoReactionThermoENS_13SpecieMixtureIN S_15reactingMixtureINS_14constTransportINS_7specie s6thermoINS_12hConstThermoINS_10perfectGasINS_6spe cieEEEEENS_16sensibleEnthalpyEEEEEEEEEE7correctEv+ 0x1e)[0x7f40b16d825e]
[cuba:03305] [ 5] chtMultiRegionFoam(+0x45eb6)[0x5559bf55deb6]
[cuba:03305] [ 6] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f40a9c21bf7]
[cuba:03305] [ 7] chtMultiRegionFoam(+0x4da2a)[0x5559bf565a2a]
[cuba:03305] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 3 with PID 0 on node cuba exited on signal 8 (Floating point exception).
----------------------------------------------------------------------------------------------------------
My case is a multi-region heat transfer problem with an inlet and outlet patch, it also takes into account condensation on certain walls.The air inside the domain begins to move due to the effect of some heat source, so I would like to know the buoyancy effect and the thermal stratification process throughout the domain.domain.png

After some searching, I believe that this is because some floating point exception, which means a division by zero (I am not quite sure). I checked all the BC and I can't find the error. Another possibility could be the geometry, but I checked it and everything is fine (checkMesh -allTopology -allGeometry). Maybe something is wrong with my thermophysical properties!

It is important to note that I tried simpler geometry and it worked fine.

I will appreciate any help thanks in advance!

Note: if you need a file I will upload it

[EFC]

Hello.

I finally got some advances. So far my case is running... turns out I had apparently bad boundary condition, specifically at the inlet (p_rgh). At the beginning I was using the "fixedFluxPressure" and it seems to get in trouble with the inlet velocity boundary condition "flowRateInletVelocity". Maybe beacause fixedFluxPressure needs flux over the patch (inlet) and it depends on the velocity and at the same time "flowRateInletVelocity" proviedes a velocity BC but is derived from the flux. I am not sure, if you have any comments it would be welcome.

Thanks !