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

simpleFoam bounding and time step continuity errors

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

Like Tree1Likes
  • 1 Post By alberto

Reply
 
LinkBack Thread Tools Display Modes
Old   November 21, 2011, 16:48
Default simpleFoam bounding and time step continuity errors
  #1
plm
New Member
 
RW
Join Date: Nov 2011
Posts: 22
Rep Power: 5
plm is on a distinguished road
Hey all,
Quite new to OF so hope I can find some help here........!

I'm trying to run a simpleFoam case on an aerofoil which I am meshing progressively finer and finer with gmsh. I am having trouble, however, seemingly when the mesh gets to a particular level of fineness.....

I am using the spalmart-allmaras model.

The timestep continuity errors shoot up massively as does the bounding of nuTilda (negative value). This is causing an exception error...

I'm also getting a similar problem in another case using k-epsilon modelling, where the bounding of epsilon and k is similarly negative and very large (again occuring when the mesh becomes fine enough).

Here is my fvSchemes file

Quote:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default steadyState;
}

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

divSchemes
{
default none;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,nuTilda) Gauss linearUpwind grad(nuTilda);
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
laplacian(1,p) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


// ************************************************** *********************** //


Can anyone give me any hints/comments?
plm is offline   Reply With Quote

Old   November 22, 2011, 03:50
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I would apply a limiter to the gradients: cellLimited Gauss linear 1; On unstructured grids, use least-squares.

Also, you might want to check your under-relaxation factors for the variables that become unbounded.

Best,
hua1015 likes this.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 22, 2011, 07:59
Default
  #3
plm
New Member
 
RW
Join Date: Nov 2011
Posts: 22
Rep Power: 5
plm is on a distinguished road
Hi alberto,
I have tried changing my gradSchemes to
Quote:
gradSchemes
{
default cellLimited Gauss linear 1;
grad(p) cellLimited Gauss linear 1;
grad(U) cellLimited Gauss linear 1;
}
And have also tried using leastSquares (I am indeed using an unstructured grid).

I've also tried reducing the relaxation factor for nuTilda (I think this is the correct thing to do but would welcome comments on why). I have

Quote:
relaxationFactors
{
default 0;
p 0.3;
U 0.7;
nuTilda 0.3;
}


I am still getting the same error though..... here is an excerpt

Quote:
Time = 148

smoothSolver: Solving for Ux, Initial residual = 0.610081558956, Final residual = 0.0400084362249, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.731335804101, Final residual = 0.0677856085541, No Iterations 4
GAMG: Solving for p, Initial residual = 0.842172295544, Final residual = 0.00537577427988, No Iterations 1
time step continuity errors : sum local = 7.39747532602e+98, global = 2.29685796544e+84, cumulative = 2.27139561979e+84
smoothSolver: Solving for nuTilda, Initial residual = 0.609600451293, Final residual = 0.00103376546081, No Iterations 8
bounding nuTilda, min: -3.43444039237e+96 max: 1.75239938872e+94 average: -7.57227581501e+91
ExecutionTime = 135.71 s ClockTime = 141 s

Time = 149

smoothSolver: Solving for Ux, Initial residual = 0.602726092692, Final residual = 0.000172607699897, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.620480876287, Final residual = 0.000181276687112, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0112844600912, Final residual = 0.000615512502725, No Iterations 3
time step continuity errors : sum local = 3.28149820498e+97, global = -1.10712455583e+82, cumulative = 2.26032437423e+84
smoothSolver: Solving for nuTilda, Initial residual = 0.182007863134, Final residual = 0.000705812932937, No Iterations 2
bounding nuTilda, min: -1.33871035339e+93 max: 1.22671944841e+94 average: 2.75116459072e+89
ExecutionTime = 136.62 s ClockTime = 142 s

Time = 150

smoothSolver: Solving for Ux, Initial residual = 0.614010972757, Final residual = 0.0373659328091, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.712775524613, Final residual = 0.0644015111041, No Iterations 4

And then it crashes with the exception error......

Hope you can help

Regards,
plm
plm is offline   Reply With Quote

Old   November 22, 2011, 16:06
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Does the code run without problems if you turn off the turbulence model? It seems none of the equations is converging. I would start checking the setup of the boundary conditions, the mesh quality (checkMesh), ...

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   November 22, 2011, 16:50
Default
  #5
plm
New Member
 
RW
Join Date: Nov 2011
Posts: 22
Rep Power: 5
plm is on a distinguished road
Hi alberto, thanks once again for the help!

I've tried running without the turbulence model and it appears I am having problems with my mesh... checkMesh turns up this error

Quote:
Checking geometry...
Overall domain bounding box (-15 -15 0) (15 15 1)
Mesh (non-empty, non-wedge) directions (1 1 0)
Mesh (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
***Boundary openness (2.14168942912e-06 3.6649966099e-06 7.81860757915e-20) possible hole in boundary description.
***Open cells found, max cell openness: 1, number of open cells 11
<<Writing 11 non closed cells to set nonClosedCells
Minumum face area = 8.07879992847e-12. Maximum face area = 1.21577118801. Face area magnitudes OK.
Min volume = 8.07879991934e-12. Max volume = 0.549816126703. Total volume = 706.22812922. Cell volumes OK.
Mesh non-orthogonality Max: 178.844177523 average: 2.93472826671
***Number of non-orthogonality errors: 11.
<<Writing 11 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 11 faces are incorrectly oriented.
<<Writing 11 faces with incorrect orientation to set wrongOrientedFaces
Max skewness = 1.00000101083 OK.

Failed 4 mesh checks.

I will continue to investigate but would welcome any comments
plm is offline   Reply With Quote

Old   November 22, 2011, 17:12
Default
  #6
plm
New Member
 
RW
Join Date: Nov 2011
Posts: 22
Rep Power: 5
plm is on a distinguished road
alberto,
I seem to be getting a problem with undefined faces in OF when using gmshToFoam which I think is causing problems later on....

Would it be possible for you to take a look at my .geo file and see what you think - I'm not sure if you're familiar with gmsh but I can't spot any problems...

Regards,
plm
Attached Files
File Type: txt naca.txt (7.9 KB, 37 views)
plm is offline   Reply With Quote

Old   November 22, 2011, 22:20
Default
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I am not very familiar with gmsh, sorry.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Reply

Tags
aerofoil, bounding, gmsh, simplefoam, time step continuity

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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
Multiple floating objects CKH OpenFOAM 10 September 21, 2011 23:13
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
directMapped problem panda60 OpenFOAM Bugs 4 July 8, 2010 10:23


All times are GMT -4. The time now is 17:26.