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

View factor radiation model. Symmetry boundaries verification..

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 1 Post By gschaider
  • 1 Post By cliffoi

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 12, 2012, 11:02
Default View factor radiation model. Symmetry boundaries verification..
  #1
New Member
 
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13
Dimoon is on a distinguished road
Hello everybody!

I simulated radiation transfer in a cylindrical chamber with 4 heaters symmetrically placed inside. For this I used a domain representing 1/4 of a model with symmetry boundary conditions. Cylinder and heater emissivity is 1. Temperatures 300 and 2500 K respectively.

I expected to get something like this sidewall.jpg. It is the result of simulation with commercial software which seems to be reasonable. Color on the side wall represents heat flux.

When I run similar problem with OpenFoam I get this OF_rad_flux_sym.jpg. Most area of cold wall have negative flux, which is unrealistic.

At the same time I performed simulation with walls having zero emissivity instead of symmetry boundaries. Under this conditions the result obtained with OpenFOAM is correct and very similar to the result of another software OF_rad_flux_0.png. This make me think that the problem is connected with treatment of symmetry boundary conditions in view factor radiation model.

I set symmetry condition for radiation by usual way as "type symmetryPlane;". Probably there is a special procedure to set symmetry for this case, does anybody know?
Dimoon is offline   Reply With Quote

Old   July 15, 2012, 17:26
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51
gschaider will become famous soon enoughgschaider will become famous soon enough
Hello Dmitry!
Quote:
Originally Posted by Dimoon View Post
This make me think that the problem is connected with treatment of symmetry boundary conditions in view factor radiation model.

I set symmetry condition for radiation by usual way as "type symmetryPlane;". Probably there is a special procedure to set symmetry for this case, does anybody know?
I'm afraid that the sad truth is that the view factor model does not support symmetry-BCs at all. At least does a quick glance at the code not show any consideration towards symmetry. And I think it should (according to my limited understanding of the model)

But may I suggest that you ask such questions next time in the "Running, Solving & CFD"-forum as it is better suited for questions on case-setup

Bernhard
pgh likes this.
gschaider is offline   Reply With Quote

Old   February 19, 2013, 18:39
Default
  #3
Member
 
Ivor Clifford
Join Date: Mar 2009
Location: Switzerland
Posts: 94
Rep Power: 17
cliffoi is on a distinguished road
I know this discussion is probably dead by now but I thought I'd add my two cents worth since I have been dealing with this exact issue recently. Dmitry, the problem you are probably encountering is that your view factors have not been calculated for a full enclosure, in which case energy is not necessarily conserved. OpenFOAM does not take symmetry boundaries into account when solving for view factor radiation. It simply takes the view factors you provide and calculates the net heat flux based on these. To fix this you need to compute view factors that include the symmetry effects. The approach I would use for your model is:
  • Copy and rotate your mesh (N faces) to obtain a full 3D model with 4N surfaces
  • Calculate the 4Nx4N view factors using this set of surfaces
  • Apply the rules of composition and distribution for view factors to combine the symmetric faces to obtain a NxN view factor matrix that includes the symmetry effects
Depending on what utilities you are using this could take a bit of effort but I guarantee that it does work.

Regards
Ivor
SamZar likes this.
cliffoi is offline   Reply With Quote

Old   March 26, 2013, 10:14
Default
  #4
New Member
 
Join Date: Jan 2013
Posts: 15
Rep Power: 13
eddi0907 is on a distinguished road
Hello,

I want to use the viewfactor radiation with chtmultiregionfoam.

I found very strange results for the model attached (due to size limitation I had to remove the mesh. Is there another possibillity to upload the complete model (1,2 Mb?)):


It's a "simple" test case, a solid sphere within a rectangular air domain surrounding, where the temperature at the outer air walls is fixed to 1000 K while the initial solid sphere and the inner air domain temperature is 300 K.

Running the case the temperature on the solid sphere becomes not physical at all. inhomogeneous distributed with local temperature drops to 226 K!!!!!
(Please see attached picture)

There should never be a temperature below 300K!

I tried OF version 2.1.1 and 2.2.0.

Any help is more than welcome.

Kind Regards.

Edmund
Attached Images
File Type: jpg temperature_after_10_sec.jpg (54.1 KB, 133 views)
Attached Files
File Type: gz radiating-sphere.tar.gz (20.2 KB, 79 views)

Last edited by eddi0907; March 26, 2013 at 10:29.
eddi0907 is offline   Reply With Quote

Old   March 27, 2013, 04:59
Default
  #5
New Member
 
Join Date: Jan 2013
Posts: 15
Rep Power: 13
eddi0907 is on a distinguished road
I changed the geometry.

Using a cube instead of a sphere it seems to work.

Except for low solid thermal conductivity it crashes initially:

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /usr/local/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

I tried changing featureAngle to 0.001 and nFacesInCoarsestLevel to 10000 in viewFactorsDict without success.

How to make a complex surface work in a radiation viewFactor model?

Or is it impossible actually even with OF 2.2.0 because of this bug: ID 656 http://www.openfoam.org/bugs/ ?
Attached Images
File Type: jpg cube_model.jpg (68.6 KB, 72 views)
File Type: jpg temperature_after_10_sec.jpg (57.3 KB, 71 views)

Last edited by eddi0907; March 27, 2013 at 08:27.
eddi0907 is offline   Reply With Quote

Old   March 26, 2014, 09:36
Default
  #6
New Member
 
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12
MangoMango is on a distinguished road
@eddi0907

..can you post your casefile in this thread? I would be really glad, because I can't find a proper viewFactor Tutorial case (the rectangle one).

Thanks in advance

Best regards,

Alex
MangoMango is offline   Reply With Quote

Old   March 31, 2014, 01:48
Default
  #7
New Member
 
Join Date: Jan 2013
Posts: 15
Rep Power: 13
eddi0907 is on a distinguished road
Hi Alex,

you can try viewfacctorgen_sphere_bug.tar.gz

on http://www.openfoam.org/mantisbt/view.php?id=656

Regards

Edmund
eddi0907 is offline   Reply With Quote

Old   March 31, 2014, 11:41
Default
  #8
New Member
 
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12
MangoMango is on a distinguished road
Hi Edmund,

thank you very much. I will give it a shot

regards

Alex
MangoMango is offline   Reply With Quote

Old   March 31, 2014, 12:50
Default
  #9
New Member
 
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12
MangoMango is on a distinguished road
Hi Edmund,

not sure and correct me if i'm wrong, but shouldn't be the 5th of the Allrun-skript line "faceAgglomerate -dict constant/viewFactorsDict" instead of "faceAgglomerate -dict constant/viewFactorsDict"

The viewFactorDict file is inconsistent too. { <- too much.

If I correct both - there is still something odd. I have too run Allrun twice and then there is only F00: 1.02709 shown... (using of-2.3.0)

Maybe my corrections are wrong. Because..

This errormsg pops up too:
Determining initial surface intersections
-----------------------------------------



--> FOAM FATAL ERROR:
Number of cells in mesh:8000 does not equal size of cellLevel:15844
This might be because of a restart with inconsistent cellLevel.

From function hexRef8::getLevel0EdgeLength() const
in file polyTopoChange/polyTopoChange/hexRef8.C at line 358.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hexRef8::getLevel0EdgeLength() const at ??:?
#3 Foam::hexRef8::hexRef8(Foam:olyMesh const&, bool) at ??:?
#4 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double, bool, Foam::refinementSurfaces const&, Foam::refinementFeatures const&, Foam::shellSurfaces const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Aborted (core dumped)


Please don't get me wrong. I don't want to play the smartypants...I only got stuck and need help. As a newbee its tough to see through this type of thermophysicalmodel

It would be really glad of you if you could help me on this

Best regards
Alex

p.s.: I ran this case on two machines to crosscheck if there are any compiler/build errors (but both gcc and icc versions are showing the same behaviour)
Attached Files
File Type: gz view_factor_testcase_sphere_bug-2.tar.gz (51.7 KB, 27 views)
MangoMango is offline   Reply With Quote

Old   April 29, 2014, 02:12
Default
  #10
New Member
 
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13
Dimoon is on a distinguished road
Quote:
Originally Posted by MangoMango View Post
@eddi0907

..can you post your casefile in this thread? I would be really glad, because I can't find a proper viewFactor Tutorial case (the rectangle one).

Thanks in advance

Best regards,

Alex
Hi Alex! I can not find it. But you can use standard tutorial: /tutorials/heatTransfer/chtMultiRegionSimpleFoam/multiRegionHeaterRadiation

I used it to prepare that case.

best regards,
Dmitry
Dimoon is offline   Reply With Quote

Old   March 23, 2021, 00:34
Default
  #11
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 107
Rep Power: 5
dasith0001 is on a distinguished road
Hi all,

just want to confirm following on the tread, the ViewFactor radiation model is still not directly applicable to a domain with a symmetryPlane ?

if not applicable, could you please direct me to a case on how to overcome the issue?

Thank you in advance
Dasith
dasith0001 is offline   Reply With Quote

Old   April 6, 2021, 03:31
Default
  #12
New Member
 
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13
Dimoon is on a distinguished road
Quote:
Originally Posted by dasith0001 View Post
Hi all,

just want to confirm following on the tread, the ViewFactor radiation model is still not directly applicable to a domain with a symmetryPlane ?

if not applicable, could you please direct me to a case on how to overcome the issue?

Thank you in advance
Dasith
I have not checked the recent versions, but couple years ago there was no changes which resolve the issue.

You can overcome it if you simulate the whole domain and do not use symmetry BC.
Dimoon is offline   Reply With Quote

Old   April 6, 2021, 03:38
Default
  #13
Senior Member
 
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 107
Rep Power: 5
dasith0001 is on a distinguished road
Hi Dmitry,

Thank you for your reply. I moved to FvDom model even it takes too much machine time. It is still better than running the whole domain.

Cheers,
Dasith
dasith0001 is offline   Reply With Quote

Reply


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
LES and combustion model Margherita Cadorin CFX 0 October 29, 2008 05:24
Low-Reynolds viscous model verification Carlos Xisto FLUENT 0 October 28, 2008 17:34
Kato-Launder model sam Main CFD Forum 13 September 21, 2006 10:15
Can We Model Spherical Symmetry in CFX ashishraval CFX 2 May 1, 2004 14:14
Advanced Turbulence Modeling in Fluent, Realizable k-epsilon Model Jonas Larsson FLUENT 5 March 13, 2000 03:27


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