CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Internal Flow Instability with layers

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 11, 2019, 03:54
Default Internal Flow Instability with layers
  #1
New Member
 
Francesco
Join Date: Jul 2014
Posts: 26
Rep Power: 10
Zephiro88 is on a distinguished road
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
File Type: zip CFTIII.zip (24.7 KB, 1 views)

Last edited by Zephiro88; June 11, 2019 at 05:47.
Zephiro88 is offline   Reply With Quote

Old   June 20, 2019, 04:20
Default
  #2
New Member
 
Francesco
Join Date: Jul 2014
Posts: 26
Rep Power: 10
Zephiro88 is on a distinguished road
Really? No one that can help me?
:/
Zephiro88 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Issues on the simulation of high-speed compressible flow within turbomachinery dowlee OpenFOAM Running, Solving & CFD 11 August 6, 2021 06:40
Internal flow simulation BCs SteveGoat FLUENT 0 May 10, 2012 06:41
Help !!! : Computation/Calculation of Heat flow through several layers of Glazing glazing Main CFD Forum 5 October 12, 2010 06:28
internal flow and external flow ? Pathway0320 FLUENT 1 November 17, 2006 03:37
Two phase Flow instability shaoping Main CFD Forum 1 September 8, 2003 13:25


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