|
[Sponsors] |
View factor radiation model. Symmetry boundaries verification.. |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 12, 2012, 11:02 |
View factor radiation model. Symmetry boundaries verification..
|
#1 |
New Member
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13 |
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? |
|
July 15, 2012, 17:26 |
|
#2 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Hello Dmitry!
Quote:
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 |
||
February 19, 2013, 18:39 |
|
#3 |
Member
Ivor Clifford
Join Date: Mar 2009
Location: Switzerland
Posts: 94
Rep Power: 17 |
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:
Regards Ivor |
|
March 26, 2013, 10:14 |
|
#4 |
New Member
Join Date: Jan 2013
Posts: 15
Rep Power: 13 |
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 Last edited by eddi0907; March 26, 2013 at 10:29. |
|
March 27, 2013, 04:59 |
|
#5 |
New Member
Join Date: Jan 2013
Posts: 15
Rep Power: 13 |
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/ ? Last edited by eddi0907; March 27, 2013 at 08:27. |
|
March 26, 2014, 09:36 |
|
#6 |
New Member
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12 |
@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 |
|
March 31, 2014, 01:48 |
|
#7 |
New Member
Join Date: Jan 2013
Posts: 15
Rep Power: 13 |
Hi Alex,
you can try viewfacctorgen_sphere_bug.tar.gz on http://www.openfoam.org/mantisbt/view.php?id=656 Regards Edmund |
|
March 31, 2014, 11:41 |
|
#8 |
New Member
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12 |
Hi Edmund,
thank you very much. I will give it a shot regards Alex |
|
March 31, 2014, 12:50 |
|
#9 |
New Member
Join Date: Sep 2013
Location: Germany
Posts: 18
Rep Power: 12 |
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) |
|
April 29, 2014, 02:12 |
|
#10 | |
New Member
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13 |
Quote:
I used it to prepare that case. best regards, Dmitry |
||
March 23, 2021, 00:34 |
|
#11 |
Senior Member
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 107
Rep Power: 5 |
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 |
|
April 6, 2021, 03:31 |
|
#12 | |
New Member
Dmitry Brovin
Join Date: Jul 2012
Posts: 4
Rep Power: 13 |
Quote:
You can overcome it if you simulate the whole domain and do not use symmetry BC. |
||
April 6, 2021, 03:38 |
|
#13 |
Senior Member
Desh
Join Date: Mar 2021
Location: Sydney
Posts: 107
Rep Power: 5 |
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 |
|
|
|
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 |