# Internal Flow Instability with layers

 Register Blogs Members List Search Today's Posts Mark Forums Read

June 11, 2019, 03:54
Internal Flow Instability with layers
#1
New Member

Francesco
Join Date: Jul 2014
Posts: 26
Rep Power: 11
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.
Code:
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\$
Attached Files
 CFTIII.zip (24.7 KB, 1 views)

Last edited by Zephiro88; June 11, 2019 at 05:47.

 June 20, 2019, 04:20 #2 New Member   Francesco Join Date: Jul 2014 Posts: 26 Rep Power: 11 Really? No one that can help me? :/

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post dowlee OpenFOAM Running, Solving & CFD 11 August 6, 2021 06:40 SteveGoat FLUENT 0 May 10, 2012 06:41 glazing Main CFD Forum 5 October 12, 2010 06:28 Pathway0320 FLUENT 1 November 17, 2006 03:37 shaoping Main CFD Forum 1 September 8, 2003 13:25

All times are GMT -4. The time now is 02:37.