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

Boundary conditions for external natural convection (chtMultiRegionFoam)

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

Like Tree2Likes
  • 2 Post By Bloerb

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 11, 2016, 08:52
Default Boundary conditions for external natural convection (chtMultiRegionFoam)
  #1
New Member
 
Joris C.
Join Date: Jan 2013
Posts: 29
Rep Power: 13
Coris is on a distinguished road
I am having trouble in defining the correct boundary conditions for an external natural convection problem. In this context, external is that the domain is relatively large and open at several boundaries. At these open boundaries, I don't know which pressure BC I have to apply.

The example case is in the attachment, along with a drawing of the setup. The case can be run using sh buildMesh.sh and chtMultiRegionFoam.

I am looking for the boundary conditions for p_rgh at the Inlet and Outlet. I have tried various combinations of totalPressure, prghPressure and prghTotalPressure, but none seem to work (mostly non-converging cases).

The only solution I see is the new prghTotalHydrostaticPressure BC added to fireFoam: https://github.com/OpenFOAM/OpenFOAM...0c819854ed27aa
but it is not yet available for chtMRF (http://bugs.openfoam.org/view.php?id=2132)

Any thoughts ?
Attached Images
File Type: png naturalConvectionTest.png (13.5 KB, 111 views)
Attached Files
File Type: zip naturalConvectionTest.zip (12.6 KB, 46 views)
Coris is offline   Reply With Quote

Old   August 11, 2016, 09:52
Default
  #2
Senior Member
 
Join Date: Sep 2013
Posts: 353
Rep Power: 21
Bloerb will become famous soon enough
I do not think this is a boundary problem. I did the following:

Running this with a courant number of 5 is crazy for this kind of simulation. A value above 1 should only be chosen if your boss tells you he wants pretty pictures by yesterday. I lowered that to 0.1. Also lowered your initital dT. In your schemes I changed div(phi,K) to upwind for good measures. I changed your pressure outlet to zeroGradient activated your heat source and it seems to work on first glance.

Also:
What are these for? Are you using a modified solver or am I missing something? Those shouldn't do anything with this solver.
Code:
    hydrostaticInitialization yes;
    nHydrostaticCorrectors 5;
Coris and wht like this.
Bloerb is offline   Reply With Quote

Old   August 11, 2016, 10:38
Default
  #3
New Member
 
Joris C.
Join Date: Jan 2013
Posts: 29
Rep Power: 13
Coris is on a distinguished road
Quote:
Originally Posted by Bloerb View Post
I do not think this is a boundary problem. I did the following:

Running this with a courant number of 5 is crazy for this kind of simulation. A value above 1 should only be chosen if your boss tells you he wants pretty pictures by yesterday. I lowered that to 0.1. Also lowered your initital dT.
Thanks a lot! That probably explains why a lot of my simulations did not converge at all.

Quote:
Originally Posted by Bloerb View Post
I changed your pressure outlet to zeroGradient activated your heat source and it seems to work on first glance.
So I have 2 zeroGradient patches, but then I don't have a reference for the pressure?
I tested this, and it diverges..
I am getting really high temperatures and the simulations breaks off.

Quote:
Originally Posted by Bloerb View Post
Also:
What are these for? Are you using a modified solver or am I missing something? Those shouldn't do anything with this solver.
Code:
    hydrostaticInitialization yes;
    nHydrostaticCorrectors 5;
At one point, I made the effort to incorporate the enhancements from fireFoam to chtMRF. These are the settings needed for that.
Coris is offline   Reply With Quote

Old   August 11, 2016, 11:22
Default
  #4
Senior Member
 
Join Date: Sep 2013
Posts: 353
Rep Power: 21
Bloerb will become famous soon enough
I also lowered the value for your heat source. 100W on such a small area seemed a bit high. I put that down to 1W. Most likely not necessary. I have run it for 10 seconds and the results were plausible. Velocity rising up on the heated wall, as well as rising temperature. You could also set a fixedValue for p_rgh on your inlet. It seems far enough away.
Bloerb is offline   Reply With Quote

Old   August 11, 2016, 11:33
Default
  #5
New Member
 
Joris C.
Join Date: Jan 2013
Posts: 29
Rep Power: 13
Coris is on a distinguished road
I did exactly the same and I am having these issues at ~ 1.05 s..

Could you post the case with your changes ?
Coris is offline   Reply With Quote

Old   May 27, 2021, 19:57
Default Pressure outlet boundary conditions for natural convection
  #6
Member
 
Mondal131211's Avatar
 
Mondal
Join Date: Sep 2018
Location: Canberra ACT
Posts: 68
Rep Power: 8
Mondal131211 is on a distinguished road
Quote:
Originally Posted by Bloerb View Post
I do not think this is a boundary problem. I did the following:

Running this with a courant number of 5 is crazy for this kind of simulation. A value above 1 should only be chosen if your boss tells you he wants pretty pictures by yesterday. I lowered that to 0.1. Also lowered your initital dT. In your schemes I changed div(phi,K) to upwind for good measures. I changed your pressure outlet to zeroGradient activated your heat source and it seems to work on first glance.

Also:
What are these for? Are you using a modified solver or am I missing something? Those shouldn't do anything with this solver.
Code:
    hydrostaticInitialization yes;
    nHydrostaticCorrectors 5;
Hi Bloerb,

Hope you are doing fine. Though it is a relatively old post, currently I am trying to solve a problem pretty much similar to that discussed in this thread. I found you suggested using zero gradient P_rgh boundary condition for the outlet boundary. After having a lot of hectic tries, I followed your steps and found a pretty good solution with this BC. I validated my model with the experimental result as well. However, the reason for writing this message to you is, is it physical to use zero gradient pressure at the outlet boundary? It might be possible, but it beyond my knowledge. I am very interested to know the reason of using a zero gradient pressure at the outlet as I want to stick with this pressure BC combination. This is what I was looking for. Is it possible to give me a brief explanation of using zero gradient pressure at the outlet? I would really very grateful to you.

Or if anybody in this group can answer my question, would appreciate it.

Regards,
Razon
Mondal131211 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
Radiation in semi-transparent media with surface-to-surface model? mpeppels CFX 11 August 22, 2019 07:30
Difficulty In Setting Boundary Conditions Moinul Haque CFX 4 November 25, 2014 17:30
Boundary conditions for natural convection in vertical channel Kader OpenFOAM Running, Solving & CFD 0 May 26, 2014 14:50
boundary conditions for Natural Convection problem Nav CFX 13 June 6, 2011 07:37
Natural convection - Inlet boundary condition max91 CFX 1 July 29, 2008 20:28


All times are GMT -4. The time now is 17:43.