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

MRFSimpleFoam wind turbine case diverges

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 22, 2012, 08:59
Default MRFSimpleFoam wind turbine case diverges
  #1
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
Hi, everyone.

I am new in OF(OF 2.1.0).
I try to run the Small Horizontal Axis Wind Turbine case using MRFSimpleFoam.

First, I simulated that case with the condition that is no rotation and no velocity, and the result is good and there is no problem. The continuity error is very very small that I expected.

However, the problem occured when I simulated that case with no rotation and 7m/s inlet-velocity. After about 30 iterations, the continuity error started to increase more than 1 and finally it became very very big. The results also was very wrong when I view the velocity and pressure fields using paraView due to bad continuity error.
I already tried to change the fvScheme and fvSolver and simulated a lot of cases with the same mesh file. but Always I failed to solve because of the same reason, 'time step continuity error'.

Before simulating this case, I had the AMI problems and resolved it by making a new mesh. So, there is no AMI problems in this case, I already checked it. And my final goal is to use pimpleDyMFoam as soon as possible.

Therefore, I need your help.
Now, I'm typing it with remote control to my computer in lab. so I couldn't show you the detail of my case. But I try to give you the details as mush as possible.

I want you to reply about
1. how I can reduce the continuity error
2. Whether the fvSchemes and fvSolution is proper or not
and 3. better ideas to solve this case.

[more informations about the best but failed case that I've solved]
the log file and files in '0' folder is attached.

1. output from the terminal near the problem

PHP Code:
Time 12
smoothSolver
:  Solving for UxInitial residual 0.24178, Final residual 4.66485e-09No Iterations 24
smoothSolver
:  Solving for UyInitial residual 0.236673, Final residual 8.17346e-09No Iterations 22
smoothSolver
:  Solving for UzInitial residual 0.232616, Final residual 5.05721e-09No Iterations 24
GAMG
:  Solving for pInitial residual 0.00439615, Final residual 4.39105e-06No Iterations 936
GAMG
:  Solving for pInitial residual 8.99011e-06, Final residual 9.74321e-07No Iterations 10
GAMG
:  Solving for pInitial residual 2.92193e-06, Final residual 9.18676e-07No Iterations 3
time step continuity errors 
sum local 1.34624, global = 0.000523313cumulative 0.000717847
smoothSolver
:  Solving for omegaInitial residual 0.330051, Final residual 4.44257e-07No Iterations 8
smoothSolver
:  Solving for kInitial residual 0.0234427, Final residual 2.08236e-07No Iterations 10
ExecutionTime 
11239.2 s  ClockTime 11277 s
Time 
13
smoothSolver
:  Solving for UxInitial residual 0.394172, Final residual 9.36361e-09No Iterations 24
smoothSolver
:  Solving for UyInitial residual 0.236768, Final residual 2.40871e-09No Iterations 24
smoothSolver
:  Solving for UzInitial residual 0.355896, Final residual 9.66341e-09No Iterations 24
GAMG
:  Solving for pInitial residual 0.00107708, Final residual 2.64334e-06No Iterations 1000
GAMG
:  Solving for pInitial residual 3.54345e-05, Final residual 1.29183e-06No Iterations 1000
GAMG
:  Solving for pInitial residual 9.14944e-06, Final residual 9.99899e-07No Iterations 39
time step continuity errors 
sum local 0.686033, global = 0.000218238cumulative 0.000936085
smoothSolver
:  Solving for omegaInitial residual 0.425149, Final residual 3.32176e-07No Iterations 12
smoothSolver
:  Solving for kInitial residual 0.0224997, Final residual 4.76061e-07No Iterations 8
ExecutionTime 
12351.3 s  ClockTime 12392 s
Time 
14
smoothSolver
:  Solving for UxInitial residual 0.195003, Final residual 2.3199e-09No Iterations 24
smoothSolver
:  Solving for UyInitial residual 0.173547, Final residual 9.25896e-09No Iterations 22
smoothSolver
:  Solving for UzInitial residual 0.19108, Final residual 2.86689e-09No Iterations 24
GAMG
:  Solving for pInitial residual 0.00154395, Final residual 1.54298e-06No Iterations 907
GAMG
:  Solving for pInitial residual 3.66792e-06, Final residual 8.65631e-07No Iterations 2
GAMG
:  Solving for pInitial residual 1.34099e-06, Final residual 7.70434e-07No Iterations 1
time step continuity errors 
sum local 4.70061, global = -0.00500594cumulative = -0.00406985
smoothSolver
:  Solving for omegaInitial residual 0.333175, Final residual 5.88616e-07No Iterations 6
smoothSolver
:  Solving for kInitial residual 0.0389014, Final residual 2.82821e-07No Iterations 8
ExecutionTime 
12859 s  ClockTime 12902 s
Time 
15
smoothSolver
:  Solving for UxInitial residual 0.178065, Final residual 6.43032e-09No Iterations 22
smoothSolver
:  Solving for UyInitial residual 0.214259, Final residual 2.16551e-09No Iterations 24
smoothSolver
:  Solving for UzInitial residual 0.196981, Final residual 7.3701e-09No Iterations 22
GAMG
:  Solving for pInitial residual 0.000258065, Final residual 9.99113e-07No Iterations 917
GAMG
:  Solving for pInitial residual 3.00182e-07, Final residual 3.00182e-07No Iterations 0
GAMG
:  Solving for pInitial residual 3.00182e-07, Final residual 3.00182e-07No Iterations 0
time step continuity errors 
sum local 53.2202, global = -0.00137399cumulative = -0.00544385
smoothSolver
:  Solving for omegaInitial residual 0.333346, Final residual 3.29517e-07No Iterations 8
smoothSolver
:  Solving for kInitial residual 0.217336, Final residual 2.4148e-07No Iterations 12
ExecutionTime 
13368.8 s  ClockTime 13413 s
Time 
16
smoothSolver
:  Solving for UxInitial residual 0.33311, Final residual 2.71606e-09No Iterations 22
smoothSolver
:  Solving for UyInitial residual 0.278035, Final residual 9.89641e-09No Iterations 22
smoothSolver
:  Solving for UzInitial residual 0.299317, Final residual 3.8482e-09No Iterations 22
GAMG
:  Solving for pInitial residual 7.77667e-06, Final residual 9.36191e-07No Iterations 6
GAMG
:  Solving for pInitial residual 1.15861e-06, Final residual 9.26386e-07No Iterations 1
GAMG
:  Solving for pInitial residual 8.8794e-07, Final residual 8.8794e-07No Iterations 0
time step continuity errors 
sum local 138.064, global = -0.143375cumulative = -0.148819
smoothSolver
:  Solving for omegaInitial residual 0.333087, Final residual 1.17684e-07No Iterations 6
smoothSolver
:  Solving for kInitial residual 0.0501571, Final residual 3.91316e-07No Iterations 12
ExecutionTime 
13391.3 s  ClockTime 13436 s
Time 
17
smoothSolver
:  Solving for UxInitial residual 0.167577, Final residual 2.71625e-09No Iterations 24
smoothSolver
:  Solving for UyInitial residual 0.132945, Final residual 3.52546e-09No Iterations 22
smoothSolver
:  Solving for UzInitial residual 0.154564, Final residual 2.4897e-09No Iterations 24
GAMG
:  Solving for pInitial residual 1.56889e-05, Final residual 8.79822e-07No Iterations 2
GAMG
:  Solving for pInitial residual 6.76478e-09, Final residual 6.76478e-09No Iterations 0
GAMG
:  Solving for pInitial residual 6.76478e-09, Final residual 6.76478e-09No Iterations 0
time step continuity errors 
sum local 427.601, global = -0.153837cumulative = -0.302656
smoothSolver
:  Solving for omegaInitial residual 0.333002, Final residual 6.87238e-07No Iterations 10
smoothSolver
:  Solving for kInitial residual 0.432678, Final residual 3.14185e-07No Iterations 10
ExecutionTime 
13411.8 s  ClockTime 13456 s 
2) structure of my case
I attached the file.

3) boundary
BLADEFIRST
{
type wall;
nFaces 28088;
startFace 32433504;
}
BLADESECOND
{
type wall;
nFaces 28088;
startFace 32461592;
}
HUB
{
type wall;
nFaces 8924;
startFace 32489680;
}
AMI_DISK_FRONT
{
type cyclicAMI;
nFaces 64812;
startFace 32498604;
matchTolerance 0.001;
neighbourPatch AMI_DOMAIN_FRONT;
transform noOrdering;
}
AMI_DISK_BACK
{
type cyclicAMI;
nFaces 64812;
startFace 32563416;
matchTolerance 0.001;
neighbourPatch AMI_DOMAIN_BACK;
transform noOrdering;
}
AMI_DISK_SIDE
{
type cyclicAMI;
nFaces 40468;
startFace 32628228;
matchTolerance 0.001;
neighbourPatch AMI_DOMAIN_SIDE;
transform noOrdering;
}
AMI_DOMAIN_FRONT
{
type cyclicAMI;
nFaces 29440;
startFace 32668696;
matchTolerance 0.001;
neighbourPatch AMI_DISK_FRONT;
transform noOrdering;
}
AMI_DOMAIN_BACK
{
type cyclicAMI;
nFaces 19968;
startFace 32698136;
matchTolerance 0.001;
neighbourPatch AMI_DISK_BACK;
transform noOrdering;
}
AMI_DOMAIN_SIDE
{
type cyclicAMI;
nFaces 11776;
startFace 32718104;
matchTolerance 0.001;
neighbourPatch AMI_DISK_SIDE;
transform noOrdering;
}
INLET_VELOCITY
{
type patch;
nFaces 38656;
startFace 32729880;
}
OUTLE_PRESSURE
{
type patch;
nFaces 29184;
startFace 32768536;
}
SYMMETRY
{
type wall;
nFaces 41984;
startFace 32797720;
}

4) MRFZones
FLUID
{
// Fixed patches (by default they 'move' with the MRF zone)
nonRotatingPatches (AMI_DISK_FRONT AMI_DISK_BACK AMI_DISK_SIDE AMI_DOMAIN_FRONT AMI_DOMAIN_BACK AMI_DOMAIN_SIDE);
origin origin [0 1 0 0 0 0 0] (0 0 0);
axis axis [0 0 0 0 0 0 0] (0 0 1);
omega omega [0 0 -1 0 0 0 0] 0;
}

5) fvSchemes
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
// default none;
// div(phi,U) Gauss limitedLinearV 1;
// div(phi,k) Gauss limitedLinear 1;
// div(phi,omega) Gauss limitedLinear 1;
// div((nuEff*dev(T(grad(U))))) Gauss linear;
default none;
div(phi,U) Gauss upwind;
div(phi,U) Gauss linearUpwind grad(U);
div(phi,k) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p ;
}

6) fvSolution
solvers
{
p
{
solver GAMG;
tolerance 1e-06;
relTol 0.001;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 20;
agglomerator faceAreaPair;
mergeLevels 1;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0;
}
k
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-06;
relTol 0;
}
omega
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-06;
relTol 0;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 2;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
U 0.5;
k 0.5;
omega 0.5;
}
}
Attached Images
File Type: jpg 1.jpg (34.8 KB, 57 views)
Attached Files
File Type: gz info.tar.gz (14.8 KB, 32 views)

Last edited by ysh1227; May 22, 2012 at 09:19.
ysh1227 is offline   Reply With Quote

Old   May 6, 2015, 09:51
Default
  #2
Senior Member
 
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 112
Rep Power: 8
be_inspired is on a distinguished road
Did you find the solution?
be_inspired is offline   Reply With Quote

Old   May 7, 2015, 10:13
Default
  #3
New Member
 
Seon Hye
Join Date: Feb 2012
Location: South Korea
Posts: 18
Rep Power: 5
ysh1227 is on a distinguished road
Unfortunately, I do not use openfoam anymore, so I have no idea about that question. Sorry.
ysh1227 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
Can I use Turbogrid to mesh a wind turbine blade? Henry Liu CFX 10 November 10, 2011 06:41
Pressure Drop/Porous Media within Shrouded Wind Turbine Simulation joshdsouza FLUENT 0 March 3, 2010 19:02
Usage of Star-CD in analysis of Wind Turbine Manabendra M. De CD-adapco 5 January 7, 2007 00:42
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
Wind turbine simulation Saturn Main CFD Forum 1 June 12, 2006 03:57


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