[mbay101]

Hi everybody,

I m having a hard Time to find the Problem that make OpenFOAM stop simulating my konvektion with chtMultiRegionFoam. It seems that the Problem have something with calcutating h in the Air. That makes my Simulation unstabil.

// * * * * * * * * * * * * *

* * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create fluid mesh for region bottomAir for time = 0
Create solid mesh for region KK for time = 0
*** Reading fluid mesh thermophysical properties for region bottomAir
Adding to thermoFluid
Selecting thermodynamics package 
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
Adding to rhoFluid
Adding to UFluid
Adding to phiFluid
Adding to gFluid
Adding to turbulence
Selecting turbulence model type laminar
Adding to ghFluid
Adding to ghfFluid
Radiation model not active: radiationProperties not found
Selecting radiationModel none
Adding fvOptions
No finite volume options present
*** Reading solid mesh thermophysical properties for region KK
Adding to thermos
Selecting thermodynamics package 
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}
Adding to radiations
Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
Adding fvOptions
No finite volume options present
Time = 1
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0244138, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0235651, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0249624, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.016444, No Iterations 3
Min/max T:292.987 300.021
GAMG: Solving for p_rgh, Initial residual = 0.804056, Final residual = 0.00263329, No Iterations 6
time step continuity errors : sum local = 0.0540164, global = 5.78109e-19, cumulative = 5.78109e-19
Min/max rho:1.15854 1.18635
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0205534, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.988 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 1.69 s ClockTime = 1 s
Time = 2
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.489481, Final residual = 0.0177715, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.215961, Final residual = 0.0114513, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.492604, Final residual = 0.0180205, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.407711, Final residual = 0.0253633, No Iterations 2
Min/max T:293.299 307.626
GAMG: Solving for p_rgh, Initial residual = 0.970476, Final residual = 0.00912021, No Iterations 6
time step continuity errors : sum local = 0.0476966, global = -7.11662e-18, cumulative = -6.53851e-18
Min/max rho:1.12815 1.19433
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.209485, Final residual = 0.00386612, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.952 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 2.25 s ClockTime = 2 s
Time = 3
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.45139, Final residual = 0.0328024, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.55513, Final residual = 0.0539866, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.454886, Final residual = 0.0282488, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.509759, Final residual = 0.0246587, No Iterations 2
Min/max T:31.9756 316.071
GAMG: Solving for p_rgh, Initial residual = 0.723399, Final residual = 0.00572431, No Iterations 5
time step continuity errors : sum local = 0.0547409, global = 3.87611e-18, cumulative = -2.6624e-18
Min/max rho:0.841215 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.109023, Final residual = 0.00190129, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.805 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 2.84 s ClockTime = 3 s
Time = 4
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.415676, Final residual = 0.00778006, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.697505, Final residual = 0.0673045, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.414124, Final residual = 0.0344727, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.628609, Final residual = 0.00417972, No Iterations 4
Min/max T:17.0031 363.37
GAMG: Solving for p_rgh, Initial residual = 0.767669, Final residual = 0.00549893, No Iterations 5
time step continuity errors : sum local = 0.108476, global = 7.559e-18, cumulative = 4.8966e-18
Min/max rho:0.53968 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0721226, Final residual = 0.0012664, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.17 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 3.41 s ClockTime = 3 s
Time = 5
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.34231, Final residual = 0.0171793, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.491816, Final residual = 0.0140244, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.355791, Final residual = 0.0160844, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.463123, Final residual = 0.0248464, No Iterations 2
Min/max T:148.512 347.184
GAMG: Solving for p_rgh, Initial residual = 0.766683, Final residual = 0.00230729, No Iterations 6
time step continuity errors : sum local = 0.0375931, global = 7.12003e-18, cumulative = 1.20166e-17
Min/max rho:0.944051 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0531948, Final residual = 0.0009174, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.29 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 4.05 s ClockTime = 4 s
Time = 6
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.407268, Final residual = 0.0273037, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.43275, Final residual = 0.0111889, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.415631, Final residual = 0.00593571, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.456496, Final residual = 0.0253663, No Iterations 2
Min/max T:212.035 335.633
GAMG: Solving for p_rgh, Initial residual = 0.744447, Final residual = 0.00689067, No Iterations 5
time step continuity errors : sum local = 0.0813514, global = -3.91602e-18, cumulative = 8.10061e-18
Min/max rho:0.983594 1.88713
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0419207, Final residual = 0.000700995, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.455 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 4.65 s ClockTime = 4 s
Time = 7
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.387934, Final residual = 0.0327896, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.490351, Final residual = 0.0114847, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.389361, Final residual = 0.0226725, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.398425, Final residual = 0.00950898, No Iterations 3
Min/max T:119.731 329.884
GAMG: Solving for p_rgh, Initial residual = 0.726224, Final residual = 0.00385633, No Iterations 5
time step continuity errors : sum local = 0.0571994, global = -8.78886e-18, cumulative = -6.88244e-19
Min/max rho:0.587975 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0344731, Final residual = 0.000555654, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.613 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 5.24 s ClockTime = 5 s
Time = 8
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.28233, Final residual = 0.00608968, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.332424, Final residual = 0.0148631, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.293554, Final residual = 0.0232985, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.336824, Final residual = 0.0129785, No Iterations 3
Min/max T:179.143 331.62
GAMG: Solving for p_rgh, Initial residual = 0.747343, Final residual = 0.00702316, No Iterations 5
time step continuity errors : sum local = 0.112819, global = -2.7136e-18, cumulative = -3.40184e-18
Min/max rho:0.913078 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0292409, Final residual = 0.000453907, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.653 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 5.87 s ClockTime = 6 s
Time = 9
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.307796, Final residual = 0.0252584, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.200522, Final residual = 0.00505223, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.324858, Final residual = 0.0291988, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.380318, Final residual = 0.0079069, No Iterations 3
Min/max T:250.067 333.011
GAMG: Solving for p_rgh, Initial residual = 0.704899, Final residual = 0.00422722, No Iterations 5
time step continuity errors : sum local = 0.0450128, global = -9.44612e-18, cumulative = -1.2848e-17
Min/max rho:0.993519 1.40521
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0253593, Final residual = 0.000378493, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.511 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 6.44 s ClockTime = 6 s
Time = 10
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.341598, Final residual = 0.0160249, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.440415, Final residual = 0.0410832, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.348379, Final residual = 0.0247462, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.31586, Final residual = 0.0288227, No Iterations 2
Min/max T:133.182 334.105
GAMG: Solving for p_rgh, Initial residual = 0.678469, Final residual = 0.00227451, No Iterations 5
time step continuity errors : sum local = 0.0222442, global = 4.5616e-19, cumulative = -1.23918e-17
Min/max rho:0.509166 1.66624
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0224042, Final residual = 0.000321794, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 296.18 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 7.02 s ClockTime = 7 s
Time = 11
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.290416, Final residual = 0.0158395, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.394223, Final residual = 0.036334, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.29554, Final residual = 0.0250645, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.27665, Final residual = 0.0269377, No Iterations 2
Min/max T:184.884 334.996
GAMG: Solving for p_rgh, Initial residual = 0.719697, Final residual = 0.00671242, No Iterations 5
time step continuity errors : sum local = 0.0882205, global = 4.84435e-18, cumulative = -7.54745e-18
Min/max rho:0.828506 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0200846, Final residual = 0.000278391, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.961 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 7.59 s ClockTime = 7 s
Time = 12
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.273298, Final residual = 0.0219519, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.221847, Final residual = 0.0197972, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.288302, Final residual = 0.0204731, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.337786, Final residual = 0.0197126, No Iterations 2
Min/max T:245.752 335.745
GAMG: Solving for p_rgh, Initial residual = 0.690738, Final residual = 0.00637376, No Iterations 5
time step continuity errors : sum local = 0.0617498, global = 1.7233e-18, cumulative = -5.82415e-18
Min/max rho:1.03632 1.44389
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0181846, Final residual = 0.000242582, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.826 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 8.16 s ClockTime = 8 s
Time = 13
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.302581, Final residual = 0.00585637, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.386143, Final residual = 0.0377692, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.30196, Final residual = 0.0255285, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.278727, Final residual = 0.025618, No Iterations 2
Min/max T:245.646 1301.59
GAMG: Solving for p_rgh, Initial residual = 0.653636, Final residual = 0.00617912, No Iterations 4
time step continuity errors : sum local = 0.0538323, global = -1.3988e-18, cumulative = -7.22295e-18
Min/max rho:0.262126 1.63589
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0171027, Final residual = 0.000284466, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.769 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 8.75 s ClockTime = 9 s
Time = 14
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.27036, Final residual = 0.0223134, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.321986, Final residual = 0.00629441, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.264422, Final residual = 0.00736481, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.283332, Final residual = 0.0224435, No Iterations 2
Min/max T:225.238 859.95
GAMG: Solving for p_rgh, Initial residual = 0.669059, Final residual = 0.00379205, No Iterations 6
time step continuity errors : sum local = 0.0360521, global = -4.18486e-19, cumulative = -7.64144e-18
Min/max rho:0.414664 1.80028
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.016315, Final residual = 0.000276046, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.743 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 9.35 s ClockTime = 9 s
Time = 15
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.214962, Final residual = 0.0152287, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.122978, Final residual = 0.0105771, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.231511, Final residual = 0.015678, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.239068, Final residual = 0.00535052, No Iterations 3
Min/max T:268.77 679.817
GAMG: Solving for p_rgh, Initial residual = 0.639297, Final residual = 0.00194893, No Iterations 6
time step continuity errors : sum local = 0.020783, global = -2.40931e-18, cumulative = -1.00508e-17
Min/max rho:0.727268 1.84172
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.016017, Final residual = 0.000270371, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.693 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 9.96 s ClockTime = 10 s
Time = 16
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.220373, Final residual = 0.00644444, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.234987, Final residual = 0.0229902, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.231552, Final residual = 0.0109475, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.179052, Final residual = 0.00462697, No Iterations 3
Min/max T:280.286 597.52
GAMG: Solving for p_rgh, Initial residual = 0.610836, Final residual = 0.00289421, No Iterations 5
time step continuity errors : sum local = 0.0171758, global = -5.11318e-18, cumulative = -1.51639e-17
Min/max rho:0.319603 1.48612
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0153671, Final residual = 0.000239808, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.691 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 10.58 s ClockTime = 10 s
Time = 17
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.26343, Final residual = 0.0217329, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.356318, Final residual = 0.0308314, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.280643, Final residual = 0.0263077, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.158732, Final residual = 0.00587434, No Iterations 3
Min/max T:-542.333 553.177
GAMG: Solving for p_rgh, Initial residual = 0.647338, Final residual = 0.00512559, No Iterations 6
time step continuity errors : sum local = 0.0307202, global = -1.32267e-18, cumulative = -1.64866e-17
Min/max rho:0.2 2
Solving for solid region KK
DICPCG: Solving for h, Initial residual = 0.0146789, Final residual = 0.00022046, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 295.707 max(T) [0 0 0 1 0 0 0] 360
ExecutionTime = 11.2 s ClockTime = 11 s
Time = 18
　
Solving for fluid region bottomAir
DILUPBiCG: Solving for Ux, Initial residual = 0.209893, Final residual = 0.00490878, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 0.153147, Final residual = 0.0036681, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.241414, Final residual = 0.00634049, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.220824, Final residual = 0.00727578, No Iterations 3
　
--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded
From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /usr2/sw/OpenFOAM//OpenFOAM-2.2.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy>::*)(double) const) const in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5 
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#6 __libc_start_main in "/lib64/libc.so.6"
#7 
in "/cax/sw/OPENFOAM/LINUX_x86_64/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"

I m thankfull for anything I get
Regards

[wyldckat]

Greetings mbay101,

I noticed that you've asked a related question here: http://www.cfd-online.com/Forums/ope...tml#post443898 post #27

OK, there isn't much information to work with here. The only thing I'm able to see is that your temperature ranges are all over the place - just look at the last few iterations:

Min/max T:179.143 331.62
Min/max T:250.067 333.011
Min/max T:133.182 334.105
Min/max T:184.884 334.996
Min/max T:245.752 335.745
Min/max T:245.646 1301.59
Min/max T:225.238 859.95
Min/max T:268.77 679.817
Min/max T:280.286 597.52
Min/max T:-542.333 553.177

The "rho" values are also all over the place, but a bit more realistic.

Please provide more information about your case, as explained here: http://www.cfd-online.com/Forums/ope...-get-help.html

Best regards,
Bruno

[mbay101]

Hi wyldckat,

thank you for offering your Help. I solve the Problem but im not satisfied with my solution and to be frank with you I don t quite understand it .
After I saw that my solution is very unstabil, I started to change variable in my case and I had succes when i increase the value of mu in thermophsicalProperties. The Problem is that I m trying to simulate the Air and when i use 1,87e-01 insteed of 1,87e-05 that change the velocity of the Air in my Geometrie Re = L.w/nu with nu is kinematic viscosity.

Do you have any idea why is my case working only with high kinematic viscosity.

Best regards,

[mbay101]

By the way,

sorry if i didn t follow http://www.cfd-online.com/Forums/ope...-get-help.html to explain my problem. This Problem is not case relevant. I always get the same problem with diffrent cases. I made a case to describe the situation. Here I tried to simulate a worm Solid in Air (free Konvektion) laminar flow.

I mesh my geometry with salome and checkMesh look very good and i never got i problem with the Mesh. I use chtMultiRegionSimpleFoam and i took the system Files from the tutorial (fvScheme, fvSolution and controlDic)

I tried to simulate the cases with laminar or turbulent kEpsilon Modell. both ways i needed to increase the viscosity nu to get the case working. how dose that effect my air velocity? or anything else in my case.

thank you for any help

[mbay101]

Sorry I forgot the case!!

OUT is the Air boundary for the limit of the Region. I can use also wall, but it dosen t seems that big of a diffrence between empty and wall with fixedValue (0 0 0) in U.

Thank you again.

[wyldckat]

Hi mbay101,

Well, if you increase the viscosity 1000x, then it acts more like a heavy liquid or like honey or something like that, which turns the flow into something very laminar and with a good+nice heat flow.

Given that your case seems to be only one solid and one fluid region, I suggest that you try to work based on this simplified multi-region case (has two fluid regions and one solid): http://openfoamwiki.net/index.php/Ge..._-_planeWall2D
The suggestion is to:

1. Keep the mesh as-is on this example case and adapt the fluid/solid properties from your case.
2. Then adjust the "blockMeshDict" file to make it more similar to your own geometry.
3. Then finally, use your mesh from Salome.

This way it'll be easier to isolate which step is breaking the run.

Best regards,
Bruno

[mbay101]

Hi Bruno,
thank you for replying. I will try your Suggestion. Another question: I always bring my Source Heat in the system by giving a fixedValue Temprature on a surface. I would Like to insert now a Heat Flux. can you show me how the Code looks like in OpenFOAM 2.2.0?

Inlet_HeatFlux
{
type compressible::turbulentHeatFluxTemperature;
heatSource flux;
q uniform 17;
kappa solidThermo;
kappaName none;
value uniform 297;
}
but this way i can t see any Temperature getting in my System what m i missing here?

[wyldckat]

Quick answer: I think what you are looking for is this: http://foam.sourceforge.net/docs/cpp...8.html#details

I found it through the modules section, under the "Wall boundary Conditions" subsection: http://www.openfoam.org/version2.2.0/documentation.php

[Sherlock_1812]

Hi Bruno,

I also have a "Foam::error:rintStack(Foam::Ostream&)" error appearing when I run my case. However I'm not able to run even one time step.

I have a modified heatTransfer solver that uses a user_defined BC for two of the four boundaries in the geometry. blockMesh compiles well. I am running a laminar case and so I've turned turbulence off in the RASproperties. However I get the following when I execute the solver.

Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Reading field T

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model laminar
Reading field alphat

Calculating field g.h

No finite volume options present


SIMPLE: convergence criteria
    field p_rgh     tolerance 0.01
    field U     tolerance 0.0001
    field T     tolerance 0.01
    field "(k|epsilon|omega)"     tolerance 0.001


Starting time loop

Time = 1e-05

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigSegv::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::fv::gaussGrad<Foam::Vector<double> >::correctBoundaryConditions(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#4  Foam::fv::gaussGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#5  Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#6  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#7  Foam::incompressible::RASModels::laminar::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"
#8  
 in "/home/srivathsan/OpenFOAM/srivathsan-2.2.2/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoamTemp"
#9  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#10  
 in "/home/srivathsan/OpenFOAM/srivathsan-2.2.2/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoamTemp"

My doubts: 1. Is this a solver specific error or generic? This is because I get a similar error when I run other modified solvers.
2. How do i rectify it?
3. What does the line #0 mean in the error?

[Sherlock_1812]

I am able to see that the the solver is not able to access libfiniteVolume.so and libincompressibleRASModels.so for some reason there by not being able to use fvSchemes properly.

Otherwise, I'm not able to make any headway with the above error. Any idea anyone?

[wyldckat]

Hi Srivaths,

The stack trace is that last part of the output you're seeing, which starts with the #0, down to whichever is the last # number. It gives you a trace of the subroutine call history, which lead to the current crash. The history is in reverse, namely, the first event was #10 and crashed somewhere near #0.

So, if we look from #0 to #10, here's what they mean:

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#0 - this is the method that handles the printing to the output screen/log file of this stack trace. In other words, you're seeing this stack trace, thanks to this method

#1  Foam::sigSegv::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"

#1 - SIGSEGV: http://en.wikipedia.org/wiki/Segmentation_fault
This is the reason why printStack was called.

#2  Uninterpreted:

#2 - Some uninterpreted machine code was found. Yes, that's exactly what it means: it's not interpreted for human comprehension

#3  Foam::fv::gaussGrad<Foam::Vector<double> >::correctBoundaryConditions(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"

#3 - Remember when you configure in "system/fvSchemes" something about a "gauss grad"? This is a method related to that. It is usually used for solving the equations the solver has defined.
This one "corrects boundary conditions", which is why it's called correctBoundaryConditions

#4  Foam::fv::gaussGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"

#4 - Also part of "gauss grad", but this one does the grad calculation... which at some point call correctBoundaryConditions.

#5  Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"

#5 - This what called the "gauss grad", since this gradScheme::grad is the more generic method that handles the "grad schemes".

#6  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"

#6 - Pretty confusing line, isn't it? This requires a person to be very well trained in coding C++, in order to figure out where the method's name is located! If you can't find it, the answer is this: the method name here is "Foam::fvc::grad"

#7  Foam::incompressible::RASModels::laminar::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so"

#7 - Nice... this one pretty much says it all:

1. divDevReff - the method.
2. laminar - the class for the turbulence model in question.
3. RASModels - the name-space for the RAS Models
4. incompressible - the name-space for the incompressible models.
5. Foam - the main name-space Foam

#8  
 in "/home/srivathsan/OpenFOAM/srivathsan-2.2.2/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoamTemp"

#8 - buoyantBoussinesqSimpleFoamTemp - this is your modified solver's name

#9  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"

#9 - Linux/library C magic Sorry, what meant to say is that this is one of the main libraries that enables us to run C/C++ code.

#10  
 in "/home/srivathsan/OpenFOAM/srivathsan-2.2.2/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoamTemp"

#10 - where all of it started... from the solver binary. Yes, it needs to start with some basic code that says something like: "I'm an executable and will libc.so be so kind as to execute the C code part of my binary form?"

As for the solution? I don't know. You haven't provided the specific boundary conditions (and initial field values) you're using, nor did you indicate the configurations you have in "fvSchemes" and "fvSolution". Nor did you give the equation you've added to the solver (and where exactly you added it), so I have no clue if there is something wrong in it.
Therefore, all I can do is guess... and my guess is that there is a value or field being initiated with 0... although this would lead to a SIGFPE, and not a SIGSEGV ...

Best regards,
Bruno

[Sherlock_1812]

Hi Bruno,

Thank you so much for that very detailed post . There were few lines which I couldn't understand in the error but you've explained them all.

I have only renamed the solver that way with my own boundary condition for a particular patch. However, I have my RASproperties file reading laminar and I've just kept the default setting for the fields k, epsilon etc. Let me take a while to go back to my case and see if I've missed anything.

Thanks a ton, again!

[Sherlock_1812]

Hi,

I've had a look at my case to find out the source of the error and correct it, but I'm not able to. The following are my fvSchemes and fvSolutions files.

About the case: Its a modified buoyant solver with my boundary condition for a free surface in the geometry. I've kept the RASproperties file set to laminar and have retained the default settings in the fvSchemes and fvSolutions folder.

I'm sure I'm missing something really simple, but what is it?

fvSchemes

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      bounded Gauss upwind;
    div(phi,T)      bounded Gauss upwind;
    div(phi,k)      bounded Gauss upwind;
    div(phi,epsilon) bounded Gauss upwind;
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         none;
    laplacian(nuEff,U) Gauss linear corrected;
    laplacian(Dp,p_rgh) Gauss linear corrected;
    laplacian(alphaEff,T) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p_rgh           ;
}

fvSolutions:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-08;
        relTol          0.01;
    }

    "(U|T|k|epsilon|R)"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-05;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;

    residualControl
    {
        p_rgh           1e-2;
        U               1e-4;
        T               1e-2;

        // possibly check turbulence fields
        "(k|epsilon|omega)" 1e-3;
    }
}

relaxationFactors
{
    fields
    {
        p_rgh           0.7;
    }
    equations
    {
        U               0.3;
        T               0.5;
        "(k|epsilon|R)" 0.7;
    }
}

[wyldckat]

Hi Srivaths,

OK, I've given a quick read to your description and it's still not enough information

I'm guessing here, but I think the problem is related to your custom boundary condition. Have a look at the following page: http://openfoamwiki.net/index.php/HowTo_debugging

Best regards,
Bruno