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

Dealing with BC's in OF 1.6

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 1, 2010, 09:17
Default Dealing with BC's in OF 1.6
  #1
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Hello everybody,
I am an almost brand new OF user and I think I have some problems to set appropriate boundary conditons for a 2d incompressible aerodinamic case. The case in question is made up by a rectangular ambient simulating the wind tunnel, with inlet and outlet sections at the two sides. Immersed into the "tunnel" there is a 2d solid profile (initially I started simulating the flow around a very simple geometry, for instance a rectangular section). The inlet value of the velocity field has to be of about 4.28 m/s, and the Reynolds number referred to profile's lenght is of about 10^5. The solver I wish to use is the pisoFoam standard one for incompressible and unsteady flows, initially coupled with a k-epsilon standard RAS model. I tried to set the BC's and the other dictionaries (contrloDict, fvSchemes and fvSolution) having a look to the tutorials, but the results I have obtained are not good, because the velocity field seems more like a potential case one rather than a turbulent one, and also because the solution reaches convergence almost instantly. Down below I post the bc's, an image of the velocity field and the other dictionaries contents, so I will be quite happy if someone of us could tell me what I'm doing wrong...
Thank you in advance


controlDict dictionary:

application pisoFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 0.1;

deltaT 1e-04;

writeControl runTime;

writeInterval 0.01;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;


fvSchemes dictionary:

ddtSchemes
{
default Euler;
}

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

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

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

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

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}


fvSolution dictionary:

solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}

pFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}

U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

UFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

}

PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
}


BC's for the velocity field:

internalField uniform (0 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform (4.28 0 0);
}

OUTLET
{
type zeroGradient;
}

SOLID WALLS
{
type fixedValue;
value uniform (0 0 0);
}


BC's for the kinematic pressure field:

internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}

OUTLET
{
type fixedValue;
value uniform 0;
}

SOLID WALLS
{
type zeroGradient;
}


BC's for the TKE (k):

internalField uniform 0.0687;

boundaryField
{
INLET
{
type fixedValue;
value uniform 0.0687;
}
OUTLET
{
type zeroGradient;
}
SOLID WALLS
{
type kqRWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0687;
}


BC's for the TKE dissipation rate (epsilon):

internalField uniform 0.109;

boundaryField
{
INLET
{
type fixedValue;
value uniform 0.109;
}
OUTLET
{
type zeroGradient;
}
SOLID WALLS
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.109;
}

BC's for the Reynolds stresses (R):

internalField uniform (0 0 0 0 0 0);

boundaryField
{
INLET
{
type fixedValue;
value uniform (0 0 0 0 0 0);
}

OUTLET
{
type zeroGradient;
}

SOLID WALLS
{
type kqRWallFunction;

}

BC's for the turbulent viscosity (nut):

internalField uniform 0;

boundaryField
{
INLET
{
type calculated;
value uniform 0;
}
OUTLET
{
type calculated;
value uniform 0;
}
SOLID WALLS
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}
Attached Images
File Type: jpg RettPisoVel.jpg (17.3 KB, 37 views)
vkrastev is offline   Reply With Quote

Old   February 2, 2010, 06:13
Default
  #2
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Any suggestions?
vkrastev is offline   Reply With Quote

Old   February 2, 2010, 06:32
Default
  #3
Senior Member
 
Jiang
Join Date: Oct 2009
Location: Japan
Posts: 186
Rep Power: 17
panda60 is on a distinguished road
Quote:
Originally Posted by vkrastev View Post
Any suggestions?

Do you have initilized your velocity field with potentialFoam ?
and the initial value of epsion could be larger , for example 100 times than the real value.
Good luck.
panda60 is offline   Reply With Quote

Old   February 2, 2010, 07:52
Default
  #4
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
Quote:
Originally Posted by panda60 View Post
Do you have initilized your velocity field with potentialFoam ?
and the initial value of epsion could be larger , for example 100 times than the real value.
Good luck.
1-No, I don't made any initialization of the velocity field. Is that always required for a transient simulation or it is simply a "good practice"?

2-I'll try to enlarge epsilon values and then I will post if there are some different results

Thanks
vkrastev is offline   Reply With Quote

Old   February 2, 2010, 08:07
Default
  #5
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20
vkrastev is on a distinguished road
No changes after enlarging epsilon by a 10^2 factor...
vkrastev is offline   Reply With Quote

Old   September 4, 2012, 11:58
Default
  #6
Member
 
Join Date: Jun 2011
Posts: 80
Rep Power: 15
maalan is on a distinguished road
vkrastev, did you finde how to solve your problem?? I am stucked with the same thing!!
Thank you in advance!
maalan 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
Porous media and adjacent BCs Carno FLUENT 1 January 14, 2010 00:54
BCs in an FSI simulation of a flexble pipe vivekcfd OpenFOAM 10 January 8, 2010 04:14
OpenFOAM Version 1.6 Released opencfd OpenFOAM Announcements from ESI-OpenCFD 0 July 27, 2009 17:55
BCs for Pressure Correction Equation (SIMPLE) Bharath Somayaji Main CFD Forum 1 March 1, 2006 06:12
Need Assistance: BCs for k-e model implementation Dean Schrage Main CFD Forum 3 August 11, 2001 01:02


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