
[Sponsors] 
February 1, 2010, 09:17 
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 
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 kepsilon 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 1e04; 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((1A(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 1e06; relTol 0.1; } pFinal { solver PCG; preconditioner DIC; tolerance 1e06; relTol 0; } U { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } UFinal { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } k { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } epsilon { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } R { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0; } nuTilda { solver PBiCG; preconditioner DILU; tolerance 1e05; 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; } 

February 2, 2010, 06:13 

#2 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 
Any suggestions?


February 2, 2010, 06:32 

#3 
Senior Member
Jiang
Join Date: Oct 2009
Location: Japan
Posts: 186
Rep Power: 17 

February 2, 2010, 07:52 

#4  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 
Quote:
2I'll try to enlarge epsilon values and then I will post if there are some different results Thanks 

February 2, 2010, 08:07 

#5 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 
No changes after enlarging epsilon by a 10^2 factor...


September 4, 2012, 11:58 

#6 
Member
Join Date: Jun 2011
Posts: 80
Rep Power: 15 
vkrastev, did you finde how to solve your problem?? I am stucked with the same thing!!
Thank you in advance! 

Thread Tools  Search this Thread 
Display Modes  


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 ESIOpenCFD  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 ke model implementation  Dean Schrage  Main CFD Forum  3  August 11, 2001 01:02 