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

Channel with Symmetry planes

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 22, 2011, 10:40
Default Channel with Symmetry planes
  #1
New Member
 
sona cs
Join Date: Aug 2010
Posts: 8
Rep Power: 6
sona is on a distinguished road
Hey all,

I'm trying to run a steady state simpleFoam case on a square channel (only quarter part) with symmetric planes. The mesh is generated in blockMeshDict.The model used is kEpsilon with wall function. The fluid used is water.The Reynolds no is 20000.The axial velocity profile in the radial direction is plotted in terms of y+ Vs U+(Axial velocity normalized with friction velocity) and compared with universal velocity profile. But the prediction in the buffer region (y+= 8-30) is not matching with universal velocity profile. Hope you can help me.

Thanks and Regards

Sona

Here I attached the blockMeshDict, U,p,k,epsilon fields, ControlDict,fvSolution and fvSchemes files used.

blockMeshDict

convertToMeters 1.0;

vertices
(
(0 0 0)
(0.005 0 0)
(0.005 0.005 0)
(0 0.005 0)
(0 0 0.5)
(0.005 0 0.5)
(0.005 0.005 0.5)
(0 0.005 0.5)

);



blocks
(
hex (0 1 2 3 4 5 6 7) (50 50 100) simpleGrading (1 1 1)

);

patches
(
patch inlet
(
(0 1 2 3)
)
patch outlet
(
(4 5 6 7)
)
wall walls
(
(1 5 6 2)
(2 6 7 3)
)
symmetryPlane sides
(
(0 1 5 4)
(0 3 7 4)
)

);

mergePatchPairs
(
);

U


internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (0 0 2);
}

outlet
{
type zeroGradient;
}

walls
{
type fixedValue;
value uniform (0 0 0);
}

sides
{
type symmetryPlane;
}
}

p file


internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 0;
}

walls
{
type zeroGradient;
}

sides
{
type symmetryPlane;
}
}

k
dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.01286;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.01286;
}
outlet
{
type zeroGradient;
}
walls
{
type kqRWallFunction;
value uniform 0.01286;
}

sides
{
type symmetryPlane;
}
}

epsilon

dimensions [0 2 -3 0 0 0 0];

internalField uniform 0.3399;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.3399;
}
outlet
{
type zeroGradient;
}
walls
{
type epsilonWallFunction;
value uniform 0.3399;
}

sides
{
type symmetryPlane;
}
}

sona is offline   Reply With Quote

Old   November 22, 2011, 10:41
Default
  #2
New Member
 
sona cs
Join Date: Aug 2010
Posts: 8
Rep Power: 6
sona is on a distinguished road
Further Details:

fvSolution
solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0;
smoother GaussSeidel;
nPreSweeps 0;
nPostSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 100;
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-08;
relTol 0;
}

epsilon
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-08;
relTol 0;
}

}

SIMPLE
{
nNonOrthogonalCorrectors 0;
convergence 1e-06;
}

relaxationFactors
{
p 0.3;
U 0.7;
k 0.7;
epsilon 0.7;

}

fvSchemes

divSchemes
{
default none;
div(phi,U) Gauss linearUpwindV Gauss linear;
div(phi,k) Gauss linearUpwind Gauss;
div(phi,epsilon) Gauss linearUpwind Gauss;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}
sona is offline   Reply With Quote

Old   November 22, 2011, 10:43
Default
  #3
New Member
 
sona cs
Join Date: Aug 2010
Posts: 8
Rep Power: 6
sona is on a distinguished road
The graph (comparison of axial velocity with universal velocity profile)
sona is offline   Reply With Quote

Old   November 22, 2011, 10:44
Default
  #4
New Member
 
sona cs
Join Date: Aug 2010
Posts: 8
Rep Power: 6
sona is on a distinguished road
sona is offline   Reply With Quote

Old   November 22, 2011, 10:50
Default
  #5
New Member
 
sona cs
Join Date: Aug 2010
Posts: 8
Rep Power: 6
sona is on a distinguished road
Graph of universal velocity profile from simulation and from correlation is attached
Attached Images
File Type: jpg U.jpg (37.3 KB, 34 views)
sona is offline   Reply With Quote

Old   February 7, 2013, 19:33
Default Y plus values
  #6
Member
 
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 6
sdharmar is on a distinguished road
Hi,
How did you get this yplus values
BR,
Suranga.
sdharmar 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
Layer addition problem near symmetry planes bigbang OpenFOAM Native Meshers: snappyHexMesh and Others 6 December 8, 2011 07:35
Boundary Conditions while using Symmetry planes garychem CFX 2 March 18, 2011 20:22
Simple Question Regarding Symmetry Planes Atella CFX 3 April 11, 2010 06:44
Simple Question Regarding Symmetry Planes Atella Main CFD Forum 0 April 9, 2010 10:58
Symmetry plane for LES in channel Flo Main CFD Forum 3 March 12, 2009 12:40


All times are GMT -4. The time now is 23:09.