CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

dnsFoam divergent solution

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 23, 2010, 05:56
Post dnsFoam divergent solution
  #1
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Hello everyone.

I have been working with dnsFoam, trying to simulate a flow between 2 walls. In a second time, I would like to focus on the boundary layer (laminar near the inlet, then transient and finally fully turbulent).
My simulations are three-dimensional, as I read on this forum a 2-D simulation does not fit with the physics.

My solution blows up after a few iterations. Here are some details about the case:

I created the mesh with blockMesh. It is a 64*64*64 mesh with 2 walls (top and bottom). I will refine the mesh but at present I try to run successfully this case.
I am using 4 patches : 1 inlet, 1 outlet and front and back faces.

Inlet velocity is 2m/s and I set the same velocity for the outlet as a test, since the primitive type "calculated" is not accepted.
Pressure is zeroGradient for every patches.

When i run my case this is the problem i am having:

Starting time loop
Time = 0.025

Number of forced K = 12
k(0.025) = 0
epsilon(0.025) = 0.03072
U.f(0.025) = 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 3.6888e-10, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 3.23237e-06, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 3.04559e-06, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 9.60124e-07, No Iterations 105
time step continuity errors : sum local = 3.05364e-07, global = 2.54977e-19, cumulative = 2.54977e-19
ExecutionTime = 35.17 s ClockTime = 37 s

Time = 0.05

Number of forced K = 12
k(0.05) = 7.03496
epsilon(0.05) = 0.10727
U.f(0.05) = 0.00848164
DILUPBiCG: Solving for Ux, Initial residual = 0.942152, Final residual = 6.78512e-06, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.447296, Final residual = 5.41539e-06, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.48381, Final residual = 5.13443e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.512288, Final residual = 7.8354e-07, No Iterations 124
time step continuity errors : sum local = 5.0736e-07, global = 3.11369e-19, cumulative = 5.66345e-19
ExecutionTime = 67.12 s ClockTime = 70 s

Time = 0.075

Number of forced K = 12
k(0.075) = 7.10567
epsilon(0.075) = 0.0846326
U.f(0.075) = 0.0171798
DILUPBiCG: Solving for Ux, Initial residual = 0.66248, Final residual = 1.28027e-07, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.478977, Final residual = 2.44826e-06, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.390294, Final residual = 2.05455e-06, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.494019, Final residual = 9.00221e-07, No Iterations 89
time step continuity errors : sum local = 1.2928e-06, global = -4.11385e-20, cumulative = 5.25207e-19
ExecutionTime = 94.88 s ClockTime = 100 s

Time = 0.1

Number of forced K = 12
k(0.1) = 29.9182
epsilon(0.1) = 0.428527
U.f(0.1) = 0.0275886
DILUPBiCG: Solving for Ux, Initial residual = 0.599916, Final residual = 4.98354e-06, No Iterations 29
DILUPBiCG: Solving for Uy, Initial residual = 0.674577, Final residual = 6.3118e-06, No Iterations 35
DILUPBiCG: Solving for Uz, Initial residual = 0.651333, Final residual = 8.71581e-06, No Iterations 32
DICPCG: Solving for p, Initial residual = 0.415806, Final residual = 5.96243e-07, No Iterations 101
time step continuity errors : sum local = 1.9221e-06, global = -5.15559e-19, cumulative = 9.64744e-21
ExecutionTime = 157.27 s ClockTime = 166 s

Time = 0.125

Number of forced K = 12
k(0.125) = 282.499
epsilon(0.125) = 4.69267
U.f(0.125) = 0.0176418
DILUPBiCG: Solving for Ux, Initial residual = 0.733885, Final residual = 0.347319, No Iterations 1001
DILUPBiCG: Solving for Uy, Initial residual = 0.907019, Final residual = 0.69476, No Iterations 1001
DILUPBiCG: Solving for Uz, Initial residual = 0.905652, Final residual = 1.06002, No Iterations 1001
DICPCG: Solving for p, Initial residual = 0.892429, Final residual = 9.26503e-07, No Iterations 136
time step continuity errors : sum local = 2.51378e-05, global = 1.83659e-17, cumulative = 1.83755e-17
ExecutionTime = 1350.73 s ClockTime = 1571 s

Time = 0.15

Number of forced K = 12
k(0.15) = 121341
epsilon(0.15) = 14327.3
U.f(0.15) = 0.497558

**** I stopped here, clearly diverging ****

As you can see, global and cumulative errors remain acceptable but k and epsilon diverge in all my tests. I have looked through the forum for similar problems but still cannot fix my problem. I would be very greatful for any suggestions to improve my case.

Also, is there any relevant boundary condition for the outlet ?
NorbertB is offline   Reply With Quote

Old   June 23, 2010, 07:38
Default
  #2
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 369
Rep Power: 10
arjun is on a distinguished road
"Pressure is zeroGradient for every patches."


if pressure is zero gradient then pressure equation is singular for incompressible case.

It might be reason for your problems, try fixed pressure at outlet.

(ideally solver should take care of such situation).
arjun is offline   Reply With Quote

Old   June 24, 2010, 04:17
Default
  #3
New Member
 
Join Date: Jun 2010
Posts: 23
Rep Power: 7
NorbertB is on a distinguished road
Hi arjun,

Thank you very much for your reply, I checked it out but it seems dnsFoam does not care about that since a DNS does not make incompressible assumption.

However, I realized that my deltaT was not small enough to keep a courant number below 1. So I changed it and now the solution converges. I should have know it ...

Now I wonder which boundary conditions could fit better to my problem, especially for the inlet/outlet.

If anybody has got an idea.
NorbertB is offline   Reply With Quote

Reply

Tags
boundary condition, divergent, dnsfoam, k-epsilon

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
smoothSolver diverges - solution in using PBiCG solver? makaveli_lcf OpenFOAM Running, Solving & CFD 3 September 11, 2013 12:44
Does FOAM converge to exact solution in Laplace equation? santiagomarquezd OpenFOAM Running, Solving & CFD 0 March 26, 2010 08:52
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16


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