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

Appropiate BC's for the inlet/outlet in an air chamber

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 31, 2015, 11:59
Default Appropiate BC's for the inlet/outlet in an air chamber
  #1
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Hello foamers,

I'm working on a case with chtMultiRegionFoam solver where I have some solid regions that are covered with a vertical air chamber in order to keep the solid regions warm and protected from the solar radiation. The air chamber keeps the surface of the solid that it is in contact with warm by natural convection (no air is pumped inside the chamber).

My main doubt is what sould be the proper BC's for velocity, pressure and temperature in the inlet and outlet patches? I haven't still worked with opened domains with OF until now...

Many thanks in advance. Any tip or hint will be much appreciated.

Regards,

Alex
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   January 31, 2015, 14:01
Default
  #2
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 286
Blog Entries: 6
Rep Power: 17
student666 is on a distinguished road
Hi, in these days I'm running a test case where i have open space too. I set the boundary of my domain far away from the zone i want to study. I think that it's good if you set freestream for U, and set the external boundaries at fixed pressure(your inital conditon for p). Anyway I'm still looking for the best ones. Look for C++ documentation at openfoam.org for other BC. Bye
student666 is offline   Reply With Quote

Old   February 4, 2015, 18:11
Default
  #3
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Hello M.C.

Thanks for your insight! Yestrday I did some tests with the inletOutlet BC and even combining inletOutlet with outletInlet at a time for both inlet and outlet. I haven't still found an answer to my problem since I can't clearly understand how theses BC's work. I'm not aware on what exactly mean inletValue and outletValue and what exactly they do... Any explanataion on that? I looked into the code but I found nothing clear to my understanding...

Many thanks!

Alex
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   February 5, 2015, 01:26
Default
  #4
Senior Member
 
M. C.
Join Date: May 2013
Location: Italy
Posts: 286
Blog Entries: 6
Rep Power: 17
student666 is on a distinguished road
Hi,

Inlet value for outletInlet BC set the velocity value at those cells if U vector is entering domain, and viceversa for inletOutlet. This is what i learned searching theads' forum. I suppose, maybe i'm wrong, these BC have been developed to take into account recirculating effects at inlet & outlet.
Have you tryed to set very low values for velocity at inlet & zeroGradient at outlet? i mean 0.001m/s so you're smoothing those effects. Changes something?
student666 is offline   Reply With Quote

Old   February 22, 2015, 13:46
Default
  #5
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Hello,

Thanks for your tips M.C. Finally I think I found a specification of BC's that may be the most suitable one for my case. At the very beggining I was trying to define the inletOutlet BC to U and here came my mistake. In my case I can't define U because this is the variable I don't have, I know that both inlet and outlet are at atmospherical pressure instead. Thus, I defined
  • p_rgh as fixedValue,
  • p as calculated,
  • U as pressureInletOutletVelocity and
  • T as inletOutlet in the outlet and fixedValue at the inlet

I also tried some other options such as defining U at the inlet as pressureInletVelocity. However, I think that the way I do it in the end is more correct since air can get in and out of the domain at the inlet and the outlet. So far, this specification seems to make air behave more or less as expected. The only problem I found so far is that some cells at the inlet have a low value of temperature, something totally absurd since at the patch the temperature is fixed and hotter (more than say 25 K, sometimes even more) than in the cell center, in the attached pict you can notice it. The temperature should raise inside the chamber along the z axis since in the external patch I set up a fixed heat flux (100W/m). In the picture you can see the geometry of my problem (Clarification: the inlet is the lower patch and the outlet is the upper one).


Maybe it has something to do with the fact that I only tested this specification during only 500 seconds, which is, in my opinion, a too small period of time taking into account that the initial conditions are not very physical (fixed value of 311.15 K in the air region and 303.15 K at the adjacent solid region). I will try a longer run of say 50000 seconds to check what happens with the temperature at the cells in the inlet patch at the end of the run and I will report it here.

If anyone thinks that this is not the best specification of BC's for my case and knows any other way that can better suit my needs, please, don't hesitate to let me know! I would be much grateful.

Regards,

Alex
Attached Images
File Type: jpg T.jpg (58.9 KB, 27 views)
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!

Last edited by zfaraday; March 2, 2015 at 17:32. Reason: 1.-Clarification added. 2.-Picture attached.
zfaraday is offline   Reply With Quote

Old   March 3, 2015, 15:03
Default
  #6
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Dear foamers,

Finally I did manage to find a more proper specification for the BC's of my case. I changed the definition for T at the inlet as shown in the following list:
  • p_rgh as fixedValue,
  • p as calculated,
  • U as pressureInletOutletVelocity (pressureNormalInletOutletVelocity was also tested with very similar results) and
  • T as inletOutlet in the outlet and also inletOutlet at the inlet. inletValue was set in the inlet with the same value used before in fixedValue BC, since T at the inlet is known. Some values of inletValue in the outlet have been tested in order to find the most suitable specifiacition.

Changing the inlet (lower patch) specification for T I solve the problem of some cells adjacent to the patch having super low and unrealistic values of T. Now the T distribution at the inlet looks nice and realistic (attached you can see U and T distributions for the case of inletValue at the outlet -upper patch- equal to the inletValue at the inlet at t=501s).







At the outlet, though, both distributions are not as nice and smooth as they are in the inlet. My guess is that the problem comes from the need to define a value of T for the incoming flux. I tried to solve this problem by increasing the value of the inletValue at the outlet since T at the outlet is higher than it is at the inlet and I thought that it made sense to define a higher value of T for the incoming flux at the outlet. However, this approach does not solve my problem but it gets perhaps even worst. Another approach I tried in order to solve this issue at the outlet was to define the upper patch as zeroGradient. Using this approach I don't need to define any value of T for the incoming flux so when I first thought about this approach I thought that this was going to be the best approach to solve my case. It was not like that, I got the worst and more unphysical results instead. In the following posts I will post the results for both approaches so that you can notice what I talk about.
Attached Images
File Type: jpg T.jpg (22.5 KB, 280 views)
File Type: jpg U.jpg (24.5 KB, 278 views)
File Type: jpg U_vec.jpg (54.3 KB, 278 views)
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!

Last edited by zfaraday; March 3, 2015 at 15:17. Reason: time at which results are taken added
zfaraday is offline   Reply With Quote

Old   March 3, 2015, 15:07
Default
  #7
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Results obtained with an inletValue at the outlet patch 5C higher than in the inlet patch (t=501).





Attached Images
File Type: jpg T.jpg (21.4 KB, 280 views)
File Type: jpg U.jpg (23.9 KB, 280 views)
File Type: jpg U_vec.jpg (53.4 KB, 280 views)
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   March 3, 2015, 15:12
Default
  #8
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Results obtained with a zeroGradient condition for T at the outlet (t=51s).






Note: No 0 is missing in the time, I had to manually stop the simulation after a lot of iterations because it was taking too long because of the U distribution that forced timeStep to be too small in order to keep Courant Number below 1.
Attached Images
File Type: jpg T.jpg (21.8 KB, 280 views)
File Type: jpg U.jpg (22.2 KB, 279 views)
File Type: jpg U_vec.jpg (33.5 KB, 280 views)
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!

Last edited by zfaraday; March 3, 2015 at 15:21. Reason: Note added
zfaraday is offline   Reply With Quote

Old   March 3, 2015, 15:33
Default
  #9
Senior Member
 
Alex
Join Date: Oct 2013
Posts: 337
Rep Power: 21
zfaraday will become famous soon enough
Now that the results have been shown I'd like to get some advice on what the best approach would be to solve my case. I don't know whether the first one I presented is the best approach or if there is another approach I haven't been able to find out that would give better results for my case. Any hint or advice will be very welcome and appreciated.

In case you need more info about my case to evaluate if I did something wrong, don't hesitate to ask! I only attached info about the inlet and the outlet because I think that the problem comes from their specification, but I may have specified wrong another patch or another thing that I may miss.

Many thanks in advance.

Best regards,

Alex
__________________
Web site where I present my Master's Thesis: foamingtime.wordpress.com

The case I talk about in this site was solved with chtMultiRegionSimpleFoam solver and involves radiation. Some basic tutorials are also resolved step by step in the web. If you are interested in these matters, you are invited to come in!
zfaraday is offline   Reply With Quote

Old   March 15, 2018, 08:17
Default
  #10
Member
 
Marco Bernardes
Join Date: May 2009
Posts: 57
Rep Power: 17
masb is on a distinguished road
Quote:
Originally Posted by zfaraday View Post
Dear foamers,

Finally I did manage to find a more proper specification for the BC's of my case. I changed the definition for T at the inlet as shown in the following list:
  • p_rgh as fixedValue,
  • p as calculated,
  • U as pressureInletOutletVelocity (pressureNormalInletOutletVelocity was also tested with very similar results) and
  • T as inletOutlet in the outlet and also inletOutlet at the inlet. inletValue was set in the inlet with the same value used before in fixedValue BC, since T at the inlet is known. Some values of inletValue in the outlet have been tested in order to find the most suitable specifiacition.

Changing the inlet (lower patch) specification for T I solve the problem of some cells adjacent to the patch having super low and unrealistic values of T. Now the T distribution at the inlet looks nice and realistic (attached you can see U and T distributions for the case of inletValue at the outlet -upper patch- equal to the inletValue at the inlet at t=501s).







At the outlet, though, both distributions are not as nice and smooth as they are in the inlet. My guess is that the problem comes from the need to define a value of T for the incoming flux. I tried to solve this problem by increasing the value of the inletValue at the outlet since T at the outlet is higher than it is at the inlet and I thought that it made sense to define a higher value of T for the incoming flux at the outlet. However, this approach does not solve my problem but it gets perhaps even worst. Another approach I tried in order to solve this issue at the outlet was to define the upper patch as zeroGradient. Using this approach I don't need to define any value of T for the incoming flux so when I first thought about this approach I thought that this was going to be the best approach to solve my case. It was not like that, I got the worst and more unphysical results instead. In the following posts I will post the results for both approaches so that you can notice what I talk about.
Hi zfaraday,

I'm trying to set up the inletOutlet BC for the temperature. I did as follow:

type inletOutlet;
value $internalField;
inletValue $internalField;

or
type inletOutlet;
value 273;
inletValue 573;

but I get the folloeing error message:

--> FOAM FATAL ERROR:

request for surfaceScalarField phi from objectRegistry region0 failed
available objects of type surfaceScalarField are
0()

From function const Type& Foam:bjectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>]
in file /home/ubuntu/OpenFOAM/OpenFOAM-5.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.

FOAM aborting

Could you kindly hellp me with this issue?

thanks!
masb is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, inlet boundary condition, outlet boundary condition

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
Air Lift Pump BCs ramesh.gannamani FLUENT 0 July 20, 2009 17:10
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59
Hydrogen Air combustion in a combustion chamber popi CFX 7 July 11, 2007 18:40
moving water from a chamber using air Matt Mason CFX 0 December 2, 2003 19:46
Help!!!problems in air jet into water bath kim FLUENT 4 June 9, 2003 07:04


All times are GMT -4. The time now is 12:08.