|
[Sponsors] |
Adding Radiations in chtMultiRegionSimpleFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 12, 2019, 11:53 |
|
#21 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Well, that will take a while
Are the dimenions of the case such small? You have 0.0005 X 0.0005 X 0.0005 m³ Each side has just 0.5 mm length! Is it realy such small? Are you using OF7.0? Regards Peter |
|
August 13, 2019, 03:51 |
|
#22 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Well, that will take a while
Quote:
I am sorry I was re-scaling it in Allrun, that's why it became so small. Now I removed that line from Allrun script. I am using OpenFOAM 4.1. |
||
August 16, 2019, 22:08 |
|
#23 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Here is your case in OF7.0
Regards Peter Last edited by peterhess; September 8, 2019 at 08:53. Reason: Case updated! |
|
August 21, 2019, 09:48 |
|
#24 |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Thank you so much for your help. I have one general question here: If I have two regions in my geometry, just like I have Now, (one is solid and one is air). And I am putting radiations on Solid region using viewFactor radiation model. and those radations will radiate from solid to air. Do I need to put radiation model also on air region to accept the radiations coming from the solid region? Thank you |
|
August 21, 2019, 11:10 |
|
#25 | ||
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Quote:
Activating the radiation in Solid calculates the radiation in solid. i.e. inside the solid! Anyway you can try it in the last model I uploaded by activating the radiation in Solid and deactivating it in Fluid ans see what happens... Quote:
There is a radiation coming from the Fluid-Side of the interface between the Solid and Fluid. If you delete the solid region and defined the same temperature on the Fluid-Side of the interface, you will get the identical results! Activating the radiation in a region lets Radiation in this region calculated. Viewfactor method does not calculates the interaction between two (or more) regions. The only method for my poor knowledge that is able doing that is FVDOM! Example: The radiation from a LED in a car head lamp through the polycarbonate front cover. Anyway, I dont know if OF is able to do that for "Semi transparent medium". Regards Peter |
|||
August 21, 2019, 11:22 |
|
#26 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Quote:
I think I am completely out of track. what I am doing is the following. I put viewFactor radiation model in the Air region. and OpaqueSolid model in the solid region. And I am putting the temperature of 500K at the solid region. Now as we know the temperature on the solid is 500K, and the temperature of the air is 300K, and I put the velocity (U) of the air equals to zero, so I am expecting to see the heat transfer through radiation from solid to air. I think this is wrong approach? And what happened is that I saw the heat transfer in paraview, but when I removed the radiations from the simulations, the results remained the same. It means that the radiations are not working. My only focus is, to simulate the radiations coming out of solid(because it is HOT), into the air. Just like we have heaters in our home. And I am unable to decide, which radiation model will work in this case. Thank you |
||
August 21, 2019, 13:33 |
|
#27 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
For sure the simulations with and without radiation are not identical!
They seams to be identical, but the wall heat flux shows the difference... Because you have constant walls temperatures: - The temperature fields will be identical - The very weak velocity field will be also identical because the very weak free convection fields is getting its buoyancy from identical temperature fields - The density is also identical, cause everywhere an identical temperature fields govers - The pressure is also identical because the temperature and density are identical (perfect gas low) Calculate the heat flux from the solid to the fluid (or at the walls of the fluid) with and without the radiation and you will find the difference. Alternativly you could change the constant temperature of the solid to constant heat generation. In this case the temperature of the solid/fluid interface will be much higher without radiation in comparing with radiation! Regards Peter |
|
August 21, 2019, 18:01 |
|
#28 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Quote:
From your this post, I understood the following: 1. If I am putting radiation (for example: viewFactor model) in one region, then it just calculates radiation inside this region, not outside this. 2. If I want to see radiations from one region to the other, then I should use fvDOM radiation model. Please let me know, If I understood wrong. Thank you |
||
August 22, 2019, 12:21 |
|
#29 | ||
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Quote:
Every region has its own viewFactorDict that just calculates the viewFactor inside it. There is no viewFactorDict for regions exchange. If Exchange between the regions to be recognized, then fvDoM must be used. By the way, you can not calculate the radiation in one region and turn it off in the other! The "other" region, in our case Solid MUST also turn the radiation on, but it as opaqueRegion to be defined! That is why opaqueRadtiation is at all exists as a radiation model... Like that the region with viewFactor calculate the radiation in it and the other region with opaque namely does not.. Quote:
Regards Peter Last edited by peterhess; August 22, 2019 at 16:41. |
|||
August 22, 2019, 17:45 |
|
#30 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Quote:
One more question please, can we also use P1 radiation model in our case? because I checked one paper and in that, the geometry was closely related to our geometry, and he is using P1 model. In case if you want to have a look, the link is given below: http://www.tfd.chalmers.se/~hani/kur...Foam_final.pdf Thank you |
||
August 22, 2019, 18:00 |
|
#31 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
||
August 25, 2019, 13:21 |
|
#32 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Quote:
Hi. I have one general question, how can we see that our radiation is working? because I have run the same case, with radiation and without radiation, but it looks exactly the same when I see temperature in Paraview. I don't know how I can check that my radiation is working or NOT? And one more question, Is it necessary to use topoSet to make cellSets when using view Factor radiation model? Thank you |
||
August 25, 2019, 14:06 |
|
#33 | ||
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Quote:
To see the difference with and without radiation, you need to calculate the heat flux leaving the solid and resieved by the outer walls of fluid. I told this way before... To calculate the wallHeatFlux leaving the interface (Fluid/Solid) and recieved by the outer boundaries of the fluid you need to use the utility: WallHeatFlux https://www.openfoam.com/documentati...lHeatFlux.html Quote:
topoSet just tells which cells belongs to which region and which boundaries belongs to which region. Regards Peter |
|||
August 25, 2019, 14:19 |
|
#34 | |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Quote:
Thank you so much for your reply. I am sorry I am asking you stupid question, but I didn't get how to use wallHeatFlux in my simulation? Where do I need to mention it? Thank you |
||
August 25, 2019, 17:35 |
|
#35 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
||
August 28, 2019, 06:30 |
|
#36 | |
New Member
Matteo Quirino
Join Date: Feb 2019
Posts: 6
Rep Power: 7 |
Quote:
Thanks a lot this is very useful to know! Can you please tell me where did you find those information? I searched for keywords like "angle" and "feat" in the shootRays.H, viewFacotrsGen.C and faceAgglomerate.C files but I couldn't find anything relevant. My C++ skills are not the best so I might have missed something. Furthermore, when I change the featAngle the rays seems that they do not change, I attached two figures, the first one with featAngle 1 and the second with 180, they are the same, I don't understand why. The explanation you gave about the faces is confirmed by the last picture, there is a cube with 36 faces, I used nFacesInCoarsestLevel 20, so 36/20=1.8 and indeed the faces are gatherd 2 by 2. Still I cannot visualize in paraview what concerns the rays. Regards Matteo Last edited by MatteoQ; August 28, 2019 at 10:10. |
||
August 30, 2019, 16:46 |
|
#37 | ||
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Hello MatteoQ!
And sorry for the late answer. Quote:
Well the source of my informations is actualy the help of Ansys-Fluent... The S2S (surface to surface) Radiation method is very good described there. See Clustering: http://www.afs.enea.it/project/neptu...th/node116.htm Quote:
If you compared the size of the F file in the region, where the view factor is calculated, between different values for the angle, you will find that the size is identical. That is because the number of patches is identical. Anyway, if you compared the values inside the F file for different angles, you will notice that they are slightly different! I can not comment that because I did not made any visualisation for the rays bevore Last edited by peterhess; September 2, 2019 at 01:45. |
|||
September 3, 2019, 03:59 |
|
#38 |
New Member
Matteo Quirino
Join Date: Feb 2019
Posts: 6
Rep Power: 7 |
Thanks again for you answer ,
I will have a look inside the F files to observe the changes you mentioned and I will investigate deeper in the openfoam codes to see if I find something. Your answers have been very usefull to me. Regards MatteoQ |
|
September 9, 2019, 09:26 |
|
#39 |
Senior Member
Raza Javed
Join Date: Apr 2019
Location: Germany
Posts: 183
Rep Power: 7 |
Here are the images... |
|
September 9, 2019, 10:25 |
|
#40 |
Senior Member
Peter Hess
Join Date: Apr 2011
Location: Austria
Posts: 250
Rep Power: 17 |
Well, you have vacuum!!!
The values on the scale are absolute values. Take a look to rho and you will see that rho in paraFoam is very low! rho has an avarage value of 6e-5. Regards Peter |
|
Tags |
chtmulitregionfoam, chtmultiregionsimpefoam, openfoam, radiation |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Incorporating pyrolysis Model in chtMultiRegionFoam | Archana V | OpenFOAM Programming & Development | 31 | April 2, 2022 02:11 |
multiRegionHeater error | ordinary | OpenFOAM Running, Solving & CFD | 2 | June 9, 2020 17:43 |
chtMultiRegionSimpleFoam: crash on parallel run | student666 | OpenFOAM Running, Solving & CFD | 3 | April 20, 2017 11:05 |
conjugateHeatFoam + interFoam | farhagim | OpenFOAM Programming & Development | 15 | July 19, 2016 07:55 |
chtMultiRegionSimpleFoam | samiam1000 | OpenFOAM Running, Solving & CFD | 39 | March 31, 2016 08:43 |