[Zephiro88]

Hello Foamers,

i write because i have a little trouble simulating an initernal flow in automotive intake. The simulation runs quite good if in do not create any layers and the results are comparable with the real model, but i'm not sure if i can get the correct separation in the divergent pipe. If i add layers, regardless the number and dimension of them, the solution is very bad, even with the turbulence model deactivated. What i saw is that the solver create a lot of sposts with a coupled very high and very low pressure area (about 10^9) while the real field is around 10^3.
It's a high reynolds simulations (U around 35 m/s)
I use porousSimpleFoam, in the attachment the setup


I addition, i noticed that after 35 itaration, even with all residual that seem to converge, the calculation blow up calculating nut.

Time = 26

smoothSolver:  Solving for Ux, Initial residual = 0.000360454, Final residual = 2.68908e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00201428, Final residual = 1.41121e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00174693, Final residual = 1.35745e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00773537, Final residual = 5.47576e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00135782, Final residual = 4.40613e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000275302, Final residual = 1.26415e-07, No Iterations 5
time step continuity errors : sum local = 7.0658e-06, global = 4.88706e-08, cumulative = 8.94817e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000551431, Final residual = 5.31541e-07, No Iterations 4
bounding omega, min: -4321.26 max: 3.64328e+06 average: 12543.6
DILUPBiCG:  Solving for k, Initial residual = 0.00678389, Final residual = 4.88653e-06, No Iterations 4
bounding k, min: -0.140953 max: 84.8419 average: 32.3289
ExecutionTime = 355.22 s  ClockTime = 356 s

Time = 27

smoothSolver:  Solving for Ux, Initial residual = 0.00034189, Final residual = 2.57485e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00192808, Final residual = 1.35612e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00164708, Final residual = 1.28467e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00741083, Final residual = 5.36712e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00130256, Final residual = 4.19825e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000266348, Final residual = 1.17273e-07, No Iterations 5
time step continuity errors : sum local = 6.59847e-06, global = 8.35547e-09, cumulative = 9.03172e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.00053183, Final residual = 3.22368e-07, No Iterations 5
bounding omega, min: -4318.24 max: 3.64328e+06 average: 12457.5
DILUPBiCG:  Solving for k, Initial residual = 0.00660004, Final residual = 5.25832e-06, No Iterations 4
bounding k, min: -0.113085 max: 84.4356 average: 31.7955
ExecutionTime = 368.54 s  ClockTime = 369 s

Time = 28

smoothSolver:  Solving for Ux, Initial residual = 0.000325084, Final residual = 2.46973e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00184644, Final residual = 1.30387e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00155611, Final residual = 1.21564e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00710823, Final residual = 5.05342e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00123285, Final residual = 3.9896e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000250351, Final residual = 1.11734e-07, No Iterations 5
time step continuity errors : sum local = 6.32847e-06, global = -1.28765e-08, cumulative = 8.90296e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000514867, Final residual = 2.06585e-07, No Iterations 5
bounding omega, min: -4315.65 max: 3.64328e+06 average: 12373.8
DILUPBiCG:  Solving for k, Initial residual = 0.00642881, Final residual = 5.95303e-06, No Iterations 4
bounding k, min: -0.173261 max: 84.1734 average: 31.2758
ExecutionTime = 381.89 s  ClockTime = 382 s

Time = 29

smoothSolver:  Solving for Ux, Initial residual = 0.000309925, Final residual = 2.37422e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00177129, Final residual = 1.25653e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.0014724, Final residual = 1.15341e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0068706, Final residual = 4.77188e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00117536, Final residual = 3.81833e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000242323, Final residual = 1.09718e-07, No Iterations 5
time step continuity errors : sum local = 6.25516e-06, global = -5.36581e-09, cumulative = 8.8493e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000498696, Final residual = 1.53754e-07, No Iterations 5
bounding omega, min: -4313.32 max: 3.64327e+06 average: 12292.3
DILUPBiCG:  Solving for k, Initial residual = 0.00626999, Final residual = 4.58916e-07, No Iterations 5
bounding k, min: -0.181461 max: 84.031 average: 30.7691
ExecutionTime = 395.31 s  ClockTime = 396 s

Time = 30

smoothSolver:  Solving for Ux, Initial residual = 0.000295946, Final residual = 2.2838e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00169958, Final residual = 1.20999e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00139532, Final residual = 1.09484e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0066614, Final residual = 4.49864e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.0011361, Final residual = 3.61009e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000240563, Final residual = 1.11431e-07, No Iterations 5
time step continuity errors : sum local = 6.3949e-06, global = -8.9467e-09, cumulative = 8.75983e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000483913, Final residual = 1.2305e-07, No Iterations 5
bounding omega, min: -4311.14 max: 3.64326e+06 average: 12212.8
DILUPBiCG:  Solving for k, Initial residual = 0.00612316, Final residual = 4.30676e-07, No Iterations 5
bounding k, min: -0.153298 max: 83.899 average: 30.2742
ExecutionTime = 408.83 s  ClockTime = 409 s

Time = 31

smoothSolver:  Solving for Ux, Initial residual = 0.000283184, Final residual = 2.19967e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.001632, Final residual = 1.16703e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00132533, Final residual = 1.04066e-06, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00650281, Final residual = 4.33612e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00110712, Final residual = 3.57645e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000236736, Final residual = 1.02741e-07, No Iterations 5
time step continuity errors : sum local = 5.93511e-06, global = -1.41155e-08, cumulative = 8.61868e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000469047, Final residual = 1.07789e-07, No Iterations 5
bounding omega, min: -4309.05 max: 3.64324e+06 average: 12135.1
DILUPBiCG:  Solving for k, Initial residual = 0.00598726, Final residual = 4.05422e-07, No Iterations 5
bounding k, min: -0.282153 max: 83.7111 average: 29.7903
ExecutionTime = 422.19 s  ClockTime = 423 s

Time = 32

smoothSolver:  Solving for Ux, Initial residual = 0.000271328, Final residual = 2.11995e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00156938, Final residual = 1.12591e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.0012604, Final residual = 9.90803e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00632632, Final residual = 4.14839e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00106764, Final residual = 3.47664e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000231964, Final residual = 1.05552e-07, No Iterations 5
time step continuity errors : sum local = 6.13756e-06, global = -9.06377e-09, cumulative = 8.52804e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000455417, Final residual = 9.65435e-08, No Iterations 5
bounding omega, min: -4307.06 max: 3.64322e+06 average: 12059.2
DILUPBiCG:  Solving for k, Initial residual = 0.00586181, Final residual = 3.83679e-07, No Iterations 5
bounding k, min: -0.237226 max: 83.4489 average: 29.3167
ExecutionTime = 435.67 s  ClockTime = 436 s

Time = 33

smoothSolver:  Solving for Ux, Initial residual = 0.00026037, Final residual = 2.04723e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00151031, Final residual = 1.08595e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00120007, Final residual = 9.4425e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.0061501, Final residual = 3.91647e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00102908, Final residual = 3.24892e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000218712, Final residual = 9.25201e-08, No Iterations 5
time step continuity errors : sum local = 5.41543e-06, global = 8.21374e-09, cumulative = 8.61018e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000441576, Final residual = 9.47487e-08, No Iterations 5
bounding omega, min: -4305.22 max: 3.6432e+06 average: 11985
DILUPBiCG:  Solving for k, Initial residual = 0.00574568, Final residual = 3.63054e-07, No Iterations 5
bounding k, min: -0.200491 max: 83.0856 average: 28.8529
ExecutionTime = 449.09 s  ClockTime = 450 s

Time = 34

smoothSolver:  Solving for Ux, Initial residual = 0.000250105, Final residual = 1.97783e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00145509, Final residual = 1.0485e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00114521, Final residual = 9.01007e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00599303, Final residual = 3.83651e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00100317, Final residual = 3.10672e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000220824, Final residual = 1.01027e-07, No Iterations 5
time step continuity errors : sum local = 5.95364e-06, global = 4.70079e-09, cumulative = 8.65719e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000428687, Final residual = 9.15733e-08, No Iterations 5
bounding omega, min: -4303.47 max: 3.64318e+06 average: 11912.5
DILUPBiCG:  Solving for k, Initial residual = 0.00563812, Final residual = 3.44069e-07, No Iterations 5
bounding k, min: -0.334767 max: 82.635 average: 28.3982
ExecutionTime = 462.61 s  ClockTime = 463 s

Time = 35

smoothSolver:  Solving for Ux, Initial residual = 0.000240681, Final residual = 1.91502e-07, No Iterations 10
smoothSolver:  Solving for Uy, Initial residual = 0.00140329, Final residual = 1.01314e-06, No Iterations 10
smoothSolver:  Solving for Uz, Initial residual = 0.00109507, Final residual = 8.61904e-07, No Iterations 10
GAMG:  Solving for p, Initial residual = 0.00586675, Final residual = 3.70285e-06, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.000988506, Final residual = 3.06432e-07, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.000213891, Final residual = 9.0713e-08, No Iterations 5
time step continuity errors : sum local = 5.38421e-06, global = 2.52467e-09, cumulative = 8.68243e-07
DILUPBiCG:  Solving for omega, Initial residual = 0.000415573, Final residual = 9.30143e-08, No Iterations 5
bounding omega, min: -4301.89 max: 3.64316e+06 average: 11841.5
DILUPBiCG:  Solving for k, Initial residual = 0.00553833, Final residual = 3.25078e-07, No Iterations 5
bounding k, min: -0.495215 max: 82.1776 average: 27.9522
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > > >::F2() const at ??:?
#6  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > > >::F23() const at ??:?
#7  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > > >::correctNut(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8  Foam::kOmegaSSTBase<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > > >::correct() at ??:?
#9  ? in "/home/bmc/OpenFOAM/OpenFOAM-v1712/platforms/linux64GccDPInt32Opt/bin/porousSimpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? in "/home/bmc/OpenFOAM/OpenFOAM-v1712/platforms/linux64GccDPInt32Opt/bin/porousSimpleFoam"
Eccezione in virgola mobile (core dump creato)
bmc@openfoam:~/OpenFOAM/bmc-v1712/run/Lambonew2$

[Zephiro88]

Really? No one that can help me?
:/