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

unphysical behaviour in buoyantSimpleFoam at inclined surfaces

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 12, 2013, 12:09
Default unphysical behaviour in buoyantSimpleFoam at inclined surfaces
  #1
Member
 
hannes
Join Date: Mar 2013
Posts: 41
Rep Power: 5
hanness is on a distinguished road
Hi there,
I wanted to run a case with the buoyantSimpleFoam solver but got some strange behaviour of the flow - at an inclined wall the hot air moves in the direction of gravity. For comparison I ran the same case with a vertical wall and got the expected result. Moreover, running the case with the inclined wall with the buoyantBoussinesqSimpleFoam solver also gives the expected results. - see attached images
inclinedWall_Uz2.jpg

straightWall_Uz.jpg

inclinedWallBoussinesq_Uz2.jpg

The boundary conditions were the following:
p_rgh
internalField uniform 101325;
boundaryField
{
hotWall
{
type buoyantPressure;
value $internalField;
}
fixedWalls
{
type buoyantPressure;
value $internalField;
}
}

U
internalField uniform (0 0 0);
boundaryField
{
hotWall
{
type fixedValue;
value uniform (0 0 0);
}
fixedWalls
{
type fixedValue;
value uniform (0 0 0);
}
}

k
internalField uniform 0.1;
boundaryField
{
hotWall
{
type compressible::kqRWallFunction;
value uniform 0.1;
}
fixedWalls
{
type compressible::kqRWallFunction;
value uniform 0.1;
}
}

epsilon
internalField uniform 0.01;
boundaryField
{
hotWall
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}
fixedWalls
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}
}

mut
internalField uniform 0.01;
boundaryField
{
hotWall
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}
fixedWalls
{
type compressible::epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.01;
}
}

Can anybody tell me whether this is a bug or my boundary conditions are not suitable?

Thanks a lot,
Hannes

Last edited by hanness; March 12, 2013 at 12:19. Reason: error in bc for k
hanness is offline   Reply With Quote

Old   March 14, 2013, 06:52
Default
  #2
RMF
New Member
 
Join Date: Aug 2012
Posts: 11
Rep Power: 5
RMF is on a distinguished road
Hi Hannes,

this looks really strange.
Which schemes did you use?

The fluid is an ideal gas?

Bye RMF
RMF is offline   Reply With Quote

Old   March 14, 2013, 10:49
Default
  #3
Member
 
hannes
Join Date: Mar 2013
Posts: 41
Rep Power: 5
hanness is on a distinguished road
Dear RMF,

yes, the fluid is calculated as ideal gas. The used schemes are as follows:

ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,h) Gauss upwind;
div(phi,K) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div((muEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(muEff,U) Gauss linear uncorrected;
laplacian((rho*(1|A(U))),p_rgh) Gauss linear uncorrected;
laplacian(alphaEff,h) Gauss linear uncorrected;
laplacian(DkEff,k) Gauss linear uncorrected;
laplacian(DepsilonEff,epsilon) Gauss linear uncorrected;
laplacian(DREff,R) Gauss linear uncorrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p_rgh;
}

anybody any idea?
Thanks a lot
Hannes
hanness is offline   Reply With Quote

Old   March 15, 2013, 10:17
Default
  #4
Member
 
hannes
Join Date: Mar 2013
Posts: 41
Rep Power: 5
hanness is on a distinguished road
thanks to Henry I got the solution:
"Inclining the wall makes the mesh non-orthogonal but you choose an uncorrected laplacian but a corrected snGrad. If you correct for non-orthogonality the flow goes in the expected direction."
hanness is offline   Reply With Quote

Old   April 2, 2013, 06:11
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18
immortality is on a distinguished road
hi Hanness.do we use unCorrected in orthogonal meshes only and for non normal faces(non-orthogonal) use corrected?should both laplacian and snGradSchemes be the same?am i correct?when we use limit?
immortality 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
twoPhaseEulerFoam unphysical behaviour GerhardHolzinger OpenFOAM Running, Solving & CFD 13 April 14, 2014 21:02
twoPhaseEulerFoam - unphysical behaviour Lydia OpenFOAM Running, Solving & CFD 14 June 14, 2012 03:21
Modeling both radiation and convection on surfaces - Ansys Transient Thermal R13 s.mishra ANSYS 0 March 31, 2012 04:12
Unphysical behaviour in dynamic mesh/vof model xaero FLUENT 0 April 24, 2011 07:50
Unstable behaviour after long period of stablility plunge11 FLUENT 1 April 6, 2011 09:15


All times are GMT -4. The time now is 16:06.