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

Problem going to second order

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 7, 2008, 19:36
Default I have a fully converged solut
  #1
Member
 
Andrew Burns
Join Date: Mar 2009
Posts: 36
Rep Power: 8
andrewburns is on a distinguished road
I have a fully converged solution (residuals between 1e-6 and 1e-8) using first order upwind discretisation. I want to switch to second order to get more accurate drag results however every time I try switching to a second order scheme my k and epsilon values go unbounded (negative) within a few iterations, stay that way and my forces go dodgy.

My discretization schemes are currently:

gradSchemes
{
default Gauss linear;
grad(p) Gauss midPoint;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss linear;
div(R) Gauss linear;
div(phi,nuTilda) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian (nuEff,U) Gauss linear corrected;
laplacian ((1|A(U)),p) Gauss linear limited 1;
laplacian (DkEff,k) Gauss linear limited 1;
laplacian (DepsilonEff,epsilon) Gauss linear limited 1;
laplacian (DREff,R) Gauss linear corrected;
laplacian (DnuTildaEff,nuTilda) Gauss linear corrected;
}

Running a simple steady-state solution in simpleFoam. My current relaxation factors are

p 0.3
U 0.5
k 0.3
epsilon 0.2
R 0.7
nuTilda 0.7

I've also tried:

QUICK
linearUpwind (which threw me an error for some reason)
limitedVanLeer (strictly bound to be positive for k and epsilon, which it appeared to ignore)
SFCD

Of all of those the most stable one appeared to be SFCD but even so k and epsilon unbounded fairly rapidly.

Can someone help me out? I know other people must be able to perform successful second order solutions in openfoam but I just can't make it work.
andrewburns is offline   Reply With Quote

Old   February 11, 2008, 18:28
Default Can anyone help me here? It se
  #2
Member
 
Andrew Burns
Join Date: Mar 2009
Posts: 36
Rep Power: 8
andrewburns is on a distinguished road
Can anyone help me here? It seems that if I take any fully converged first order solution and then switch to second order my k and epsilon (and omega when using komegaSST) go unbounded and stay there.
andrewburns is offline   Reply With Quote

Reply

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
Higher order polynomial problem... Markas FLUENT 0 February 20, 2008 17:44
1st, 2nd order or other Ivan FLUENT 3 January 28, 2008 10:45
4th order discretization of 2nd order difference dusky.He Main CFD Forum 3 February 14, 2007 09:00
first order time, second order space sven CFX 1 October 27, 2005 13:05
1st order temporal & 2nd order spatial Prateep Chatterjee FLUENT 0 January 19, 2003 01:31


All times are GMT -4. The time now is 12:27.