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

Mixed Convection Turbulent Boundary Layer Problem

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

Like Tree1Likes
  • 1 Post By morard

Reply
 
LinkBack Thread Tools Display Modes
Old   March 4, 2013, 18:52
Default Mixed Convection Turbulent Boundary Layer Problem
  #1
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
Hello Foamers... I am having a problem with my mixed convection turbulent boundary layer case. It seems the free stream velocity that I prescribe at the bottom of my computational domain is unable to "flow through" the domain. As you can see in the attached file, I prescribed a free stream of 0.2m/s at the southern boundary with a 10cm entrance length and the plate is at a higher temperature than the ambient. The more I look into it, I suspect that the BC at the east (opposite to the wall) might be the problem. The BC for U there is pressureInletOutletVelocity and pd is totalPressure. I am employing the pressureInletOutletVelocity to allow for entrainment due to the upflow of mass near the wall. Can anyone give me some feedback on what might be causing this. Two more things, I have expanded the domain and my simulation crashes, and when I prescribe the internal field to 0.2m/s, the BC on the east does not like it and it forces the streamwise velocity to zero.
Attached Images
File Type: jpg MixedConvection_FreeStream_0_2m_s.jpg (38.5 KB, 53 views)
deji is offline   Reply With Quote

Old   March 8, 2013, 15:26
Default
  #2
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
So no one can give me any feedback at all?
deji is offline   Reply With Quote

Old   March 9, 2013, 15:24
Default
  #3
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
It seems like there's nothing in OpenFOAM that caters to combined convection flows. Am I correct? All flows aren't strictly forced or natural convection.
deji is offline   Reply With Quote

Old   March 15, 2013, 10:22
Default
  #4
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
Problem solved!!! I had to make some changes to the modified pressure equation pd. It feels great. If anyone in the future encounters a problem in simulating the classical mixed convection turbulent boundary layer, I might be able to assist.
deji is offline   Reply With Quote

Old   November 27, 2013, 15:26
Default
  #5
Member
 
Jubayer
Join Date: Oct 2009
Location: The University of Western Ontario, London, Ontario
Posts: 42
Blog Entries: 1
Rep Power: 7
cm_jubayer is on a distinguished road
Hi deji,

Can you please let us know the way you solved this mixed convection problem? I am also trying to simulate a low Re flow where the natural convection cannot be neglected compared to the forced convection. Thanks.

Jubayer
cm_jubayer is offline   Reply With Quote

Old   November 27, 2013, 15:32
Default
  #6
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
Hey. It would be profoundly difficult for me to assist you without any specific questions.
deji is offline   Reply With Quote

Old   November 27, 2013, 16:33
Default
  #7
Member
 
Jubayer
Join Date: Oct 2009
Location: The University of Western Ontario, London, Ontario
Posts: 42
Blog Entries: 1
Rep Power: 7
cm_jubayer is on a distinguished road
Hi deji,

You said you made some changes to the pressure. What are the changes? What are your pressure boundary conditions? or, if possible, Can you please upload a test case so that I can look into your boundary conditions? Thanks.

Jubayer
cm_jubayer is offline   Reply With Quote

Old   March 13, 2014, 05:32
Default
  #8
New Member
 
Join Date: Oct 2013
Posts: 3
Rep Power: 3
Thomas William is on a distinguished road
Hi Deji,

I am trying to simulate water flow in a storage tank with forced and natural convection.

transient, compressible: variable density (T), turbulent

Im not having great results with buoyantPimpleFoam? What solver did you use and you mention some pressure modifications?

Thank you
Tom
Thomas William is offline   Reply With Quote

Old   March 13, 2014, 10:56
Default
  #9
Senior Member
 
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 7
deji is on a distinguished road
I started with a solver called fireFoam, and I eventually had to modify the code to enable the computation of mixed convection turbulent wall flows. I called this code lowMachFoam. Regarding the pressure modification, I essentially took out the compressibility for stability. I think you may have success with the utilization of the fireFoam code, give it a try!

Cheers,
Deji
deji is offline   Reply With Quote

Old   April 9, 2014, 03:27
Default
  #10
Member
 
Dejan Morar
Join Date: Nov 2010
Posts: 78
Rep Power: 7
morard is on a distinguished road
Hi,

I'm trying to simulate a turbulent mixed convection case so that I have decided to post my question here instead of creating a new thread.

My case is turbulent mixed convection to water in a square duct. Heating is asymmetric (just one wall is heated while the other three are not). For
temperature I use fixedGradient at the heated wall while the other walls are adiabatic. fixedValue at the inlet and advective at the outlet are
use for both temperature and velocity. For p_rgh I use buoyantPressure at the walls, zeroGradient at the inlet and fixedValue at the oulet
(I have also tried totalPressure at the outlet without success). Solver is buoyantPimpleFoam (term dpdt is excluded from the energy equation).
SGS model is homogeneousDynOneEqEddy. OpenFOAM version is 2.1.x.

The problem is: as soon as I introduce the heating at the wall, it comes to the unphysical increase of the velocity at the outlet and to the back-flow
from the outlet boundary.

Unfortunately, I cannot upload an image at the moment. The geometry is alike the one form Deji with the difference that I have four walls in my case.

Does anybody know what might cause the problem and how to solve it?

Regards,
Dejan

Last edited by morard; April 15, 2014 at 05:10.
morard is offline   Reply With Quote

Old   April 15, 2014, 05:24
Default
  #11
Member
 
Dejan Morar
Join Date: Nov 2010
Posts: 78
Rep Power: 7
morard is on a distinguished road
Here are the images at the outlet boundary.
Attached Images
File Type: png U_no_heating.png (43.5 KB, 16 views)
File Type: png U_with_heating.png (33.4 KB, 17 views)
morard is offline   Reply With Quote

Old   May 26, 2014, 15:20
Default
  #12
New Member
 
korichi Abdelkader
Join Date: May 2014
Posts: 4
Rep Power: 3
Kader is on a distinguished road
Dear Morad,
I have similar difficulty, I think the problem in u, p, p_rgh boundary conditions at the inlet and outlet. I have tried differents combinaisons, but can't find the adequates boundary.
best regards
Kader is offline   Reply With Quote

Old   May 27, 2014, 03:25
Default
  #13
Member
 
Dejan Morar
Join Date: Nov 2010
Posts: 78
Rep Power: 7
morard is on a distinguished road
Dear Kader,

The only combination of boundary conditions that works for my case is:

p_rgh

inlet
{
type totalPressure;
U U;
phi phi;
rho rho;
psi none;
gamma 1.4;
p0 uniform 101325;
value uniform 101325;
}

outlet
{
type buoyantPressure;
gradient uniform 0;
value uniform 101325;
}

U

inlet
{
type pressureInletOutletVelocity;
phi phi;
value uniform (0 0 0.044);
}

outlet
{
type pressureInletOutletVelocity;
phi phi;
value uniform (0 0 0);
}

Four outlet velocity you can also apply inletOutlet or advective. But, I do not know how to specify inlet velocity at some specific value.

Regards,
Dejan
mgg likes this.
morard is offline   Reply With Quote

Old   May 27, 2014, 06:24
Default
  #14
New Member
 
korichi Abdelkader
Join Date: May 2014
Posts: 4
Rep Power: 3
Kader is on a distinguished road
Thank you Morad, I will try it
Kader is offline   Reply With Quote

Old   May 29, 2014, 09:51
Default
  #15
New Member
 
korichi Abdelkader
Join Date: May 2014
Posts: 4
Rep Power: 3
Kader is on a distinguished road
Dear Dejan,
I have used the recomanded boundary condition, I have the flowing error message:
request for volScalarField rho from objectRegistry rgion0 failed availabe objects of type volScalarField are
13
(rhok
nut ...
I think the error in p_rgh boundary condition
have you any idea
thanks
Kader
Kader 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
problem with boundary layer usefkeyghobadi OpenFOAM Running, Solving & CFD 4 April 23, 2012 10:25
Boundary Layer roughness/ low reynolds wall treatment Luigi_ STAR-CCM+ 1 March 14, 2012 09:40
Set inlet boundary conditions in turbulent boundary layer v2f run kostas Main CFD Forum 0 March 9, 2012 11:24
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
[GAMBIT] 3D Boundary Layer Laminar and Turbulent meshing Harald D ANSYS Meshing & Geometry 1 July 7, 2009 06:20


All times are GMT -4. The time now is 15:29.