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

Solar Load and fvDOM Radiation Modeling

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

Like Tree8Likes
  • 3 Post By Shadab_Ilahi
  • 3 Post By Shadab_Ilahi
  • 1 Post By Shadab_Ilahi
  • 1 Post By YannGuyot

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 22, 2020, 09:17
Smile Solar Load and fvDOM Radiation Modeling
  #1
New Member
 
MD SHADAB
Join Date: Oct 2018
Posts: 7
Rep Power: 6
Shadab_Ilahi is on a distinguished road
Hello Everyone,


Currently, I am simulating a room which has a one glass wall through which solar radiation enters into the room. I am using an OpenFOAMv-1906.


I am able to give the correct input for solar load calculation which I can see its effects through q_primaryRadiation_0 variable and when I run this simulation along with fvDOM radiation model The surfaces on which solar radiation are falling showing the very less temperature compare to other unexposed surfaces, for reference I have attached the comparative T and Solar Load on the surface images.

I am using an externalWallHeatFluxTemperature boundary condition by which I can specify the flux on all Walls, on adiabatic walls q=0 and qr =qr similarly, on some other walls I specified q=x, qr=qr, where x is some value.

When we apply only Solar Load model alone the q_primaryRadiation_0 is equal to qr and it is a positive value and it gives very high temperatures on exposed surfaces. Whereas with Solar Load + fvDOM, the qr is negative and it gives the low temperature on exposed surfaces.

My main question is why the sun directs exposed surfaces showing ver less temperature when I choose Solar Load + Radiation Model?

For Reference purpose, I am attaching the T and radiationProperties and boundaryRadiationProperties files also.

Looking for your help and thanks in advance.

Regards
Shadab
Attached Images
File Type: jpg T-SolarLoad.jpg (54.1 KB, 207 views)
File Type: jpg qr.jpg (38.3 KB, 143 views)
Attached Files
File Type: zip T.zip (683 Bytes, 108 views)
File Type: zip radiationProperties.zip (1.5 KB, 128 views)
File Type: zip boundaryRadiationProperties.zip (765 Bytes, 127 views)
francescomarra, SharpS and Abs11 like this.
Shadab_Ilahi is offline   Reply With Quote

Old   March 31, 2020, 10:13
Default
  #2
Member
 
Yann Guyot
Join Date: Feb 2019
Posts: 37
Rep Power: 5
YannGuyot is on a distinguished road
Hello,
Could you find an answer to your problem?
cheers
YannGuyot is offline   Reply With Quote

Old   April 2, 2020, 09:50
Default
  #3
New Member
 
MD SHADAB
Join Date: Oct 2018
Posts: 7
Rep Power: 6
Shadab_Ilahi is on a distinguished road
Hi,

I have resolved that issue of showing very low temperature on sun direct exposed surfaces.

The flux/power what I am specifying by using externalWallHeatFluxTemperature BC is for convection transfer only, As we don't know on which faces sun radiation will fall and final load will be the sum of convention plus solar loads.

So I slightly changed the externalWallHeatFluxTemperature BC which consider the solar load by looking into solar load file which is generated by running the solar load radiation model. The solar load written on the boundary patches only and for inside the domain, it is zero.

Similarly, I changed the solver also which carry this solar load field with it, actually, once a solar field is generated it is acting as heat source field which is applied only on the boundaries and once the simulation becomes stable run it with radiation model. Because there will be a very high temperature on the boundaries.

Regards,
Shadab
Shadab_Ilahi is offline   Reply With Quote

Old   April 2, 2020, 10:01
Default
  #4
Member
 
Yann Guyot
Join Date: Feb 2019
Posts: 37
Rep Power: 5
YannGuyot is on a distinguished road
Hello,

Thank you, I was asking you since I want to do something similar.
I am not specifically inserted in solar load, but a have a very hot source in my domain and I want to include mirrors and windows inside the domain (as baffles) to speculary reflect or transmit radiation, and I have seen that you use such things in your boundaryRadiationProperties.
I have been trying to use such file in my simulations but it seems to not be taken into account in the solver.
There is no tutorial in the source about such problem, so I was wondering if you could explain me a bit about that.
Thank you a lot in advance
Yann
YannGuyot is offline   Reply With Quote

Old   April 3, 2020, 05:50
Default
  #5
New Member
 
MD SHADAB
Join Date: Oct 2018
Posts: 7
Rep Power: 6
Shadab_Ilahi is on a distinguished road
Hi Yann,

To consider the radiation effects on your boundary condition use externalWallHeatFluxTemperature BC with (qr qr) instead of (qr none) when radiation is on.

You can't see the radiation effect on constant temperature and zero gradient type of boundary condition. So, if you want to specify adiabatic wall use flux or power mode and give zero value and make (qr qr) in BC. Be make sure heat flux or power what you specified is the steady heat transfer through that patch with a combination of both convection plus radiation.

Currently, OpenFOAM supports grey diffusive/non-grey diffusive walls only, it means you can see the effect of radiation emission and reflection from the walls but not transmit through it, which is available in fluent as the semi-transparent boundary condition. But by using the above boundary condition we can interact with the environment with available qr on the patch and can give heat through radiation mode also.

I hope it solves your problem!

Cheers
Shadab
francescomarra, lumpor and wht like this.
Shadab_Ilahi is offline   Reply With Quote

Old   April 3, 2020, 06:32
Default
  #6
Member
 
Yann Guyot
Join Date: Feb 2019
Posts: 37
Rep Power: 5
YannGuyot is on a distinguished road
Hello Shadab,


I am little bit confused, because in your case you use this entry:
Code:
ref_wall
{
    type   opaqueReflective;

    // Fraction of the reflected is diffussive
   fd          0.0; // 0: all specular 1: all diffusive

    wallAbsorptionEmissionModel
    {
        type   multiBandAbsorption;
        absorptivity   (0.03 0.07);
        emissivity     (0.03 0.07);
    };
}

or



Code:
    type        transparent;
    wallAbsorptionEmissionModel
    {    
        type            constantAbsorption;
        absorptivity    0.0;
        emissivity      1.0;
    };

Which if I understand well, should do speculat reflection or tramsmissivity, no?


I also have another question, I already use externalWallHeatFluxTemperature BC with (qr qr), but I don't really understand how the material properties of the heated wall taken into account?


Cheers
Yann
YannGuyot is offline   Reply With Quote

Old   April 4, 2020, 04:45
Default
  #7
New Member
 
MD SHADAB
Join Date: Oct 2018
Posts: 7
Rep Power: 6
Shadab_Ilahi is on a distinguished road
Hi Yann,

The type of radiation wall properties that I have specified in boundary radiation properties such as opaqueReflective and transparent are specifically used only with solar load model.

opaqueReflective : This boundary condition gives tranparent = 0, absorptivity= x and reflectivity = 1-absorptivty. another factor fd is there which decide how the reflective component distribute as diffusive or a specular.

transparent : complete transparent =1.

As solar load is not an alone radiation model it must be run with another radiation model such as P1, fvDOMor viewFactor. As I am not familiar with other models except fvDOM So, I can't comment much on it.


In fvDOM model when we discretize the domain into multiple solid angles then we solve RTE for multiple radiation intensity rays, for that we provide greyDiffusiveRadiation or nonGreyDiffusiveRadiation BC which is for only Opaque walls which absorbs and reflect only and it read the emissivity property (emissivity = absorptivity) from boundaryRadiationProperties file only. It performs transmissivity calculation only solar load is true otherwise it is a simple opaque wall BC. for more information visit the BC implementation.

So, with the best of my knowledge, you can't apply a transparent/semitransparent type of BC with fvDOM in OpenFOAMv-6/v1912 current versions.

Hope I have clear your doubts.


Cheers,
Shadab
wht likes this.
Shadab_Ilahi is offline   Reply With Quote

Old   April 7, 2020, 15:21
Default
  #8
Member
 
Yann Guyot
Join Date: Feb 2019
Posts: 37
Rep Power: 5
YannGuyot is on a distinguished road
Hello,


Thank you very much for your complete answer, I still don't understand how such BC can only be used for solar load if it uses fvdom as a radiation model anyway, but there is no big deal and I will consider that I can't use reflective walls in my simulations.
Regarding my other question, I still don't really understand how externalwallheatflux BC works, I put flux equal to zero as you suggested but the temperature goes really high (with qr=300W/m2), I looked deeper into the code and found something about relaxation and qrRelaxation but I don't really get their meaning and how it affects results.
Thank you again in advance
Yann
tian likes this.
YannGuyot is offline   Reply With Quote

Old   July 27, 2020, 23:02
Default
  #9
Member
 
chengan.wang
Join Date: Jan 2016
Location: china
Posts: 47
Rep Power: 9
wangchengan2003 is on a distinguished road
Send a message via Skype™ to wangchengan2003
Dear Shadab,

I want to know that which solver you used for your problem?

Best regards

Chengan
wangchengan2003 is offline   Reply With Quote

Old   August 11, 2020, 09:56
Default
  #10
New Member
 
Sammy Breen
Join Date: Jul 2020
Posts: 4
Rep Power: 4
SBreen is on a distinguished road
Hello,

i am just starting out with the solar load feature of OpenFOAM and i was wondering if you could share your case files, because I am working with a simialar setup. Thank you very much.
SBreen is offline   Reply With Quote

Old   November 24, 2021, 00:56
Default
  #11
Member
 
Bushra Rasheed
Join Date: Dec 2020
Posts: 97
Rep Power: 4
B_R_Khan is on a distinguished road
Quote:
Originally Posted by Shadab_Ilahi View Post
Hi Yann,

The type of radiation wall properties that I have specified in boundary radiation properties such as opaqueReflective and transparent are specifically used only with solar load model.

opaqueReflective : This boundary condition gives tranparent = 0, absorptivity= x and reflectivity = 1-absorptivty. another factor fd is there which decide how the reflective component distribute as diffusive or a specular.

transparent : complete transparent =1.

As solar load is not an alone radiation model it must be run with another radiation model such as P1, fvDOMor viewFactor. As I am not familiar with other models except fvDOM So, I can't comment much on it.


In fvDOM model when we discretize the domain into multiple solid angles then we solve RTE for multiple radiation intensity rays, for that we provide greyDiffusiveRadiation or nonGreyDiffusiveRadiation BC which is for only Opaque walls which absorbs and reflect only and it read the emissivity property (emissivity = absorptivity) from boundaryRadiationProperties file only. It performs transmissivity calculation only solar load is true otherwise it is a simple opaque wall BC. for more information visit the BC implementation.

So, with the best of my knowledge, you can't apply a transparent/semitransparent type of BC with fvDOM in OpenFOAMv-6/v1912 current versions.

Hope I have clear your doubts.


Cheers,
Shadab
Hi!

Thanks for explaining solar load radiation boundaries. Can you tell what does opaqueDiffuisve boundary condition do? Isn't it the same as opaque reflective since opaque reflective too has diffusive fd option?

THanks
B_R_Khan is offline   Reply With Quote

Old   June 10, 2022, 13:36
Default
  #12
New Member
 
Abbes
Join Date: Apr 2022
Posts: 5
Rep Power: 2
Abs11 is on a distinguished road
Quote:
Originally Posted by Shadab_Ilahi View Post
Hello Everyone,


Currently, I am simulating a room which has a one glass wall through which solar radiation enters into the room. I am using an OpenFOAMv-1906.


I am able to give the correct input for solar load calculation which I can see its effects through q_primaryRadiation_0 variable and when I run this simulation along with fvDOM radiation model The surfaces on which solar radiation are falling showing the very less temperature compare to other unexposed surfaces, for reference I have attached the comparative T and Solar Load on the surface images.

I am using an externalWallHeatFluxTemperature boundary condition by which I can specify the flux on all Walls, on adiabatic walls q=0 and qr =qr similarly, on some other walls I specified q=x, qr=qr, where x is some value.

When we apply only Solar Load model alone the q_primaryRadiation_0 is equal to qr and it is a positive value and it gives very high temperatures on exposed surfaces. Whereas with Solar Load + fvDOM, the qr is negative and it gives the low temperature on exposed surfaces.

My main question is why the sun directs exposed surfaces showing ver less temperature when I choose Solar Load + Radiation Model?

For Reference purpose, I am attaching the T and radiationProperties and boundaryRadiationProperties files also.

Looking for your help and thanks in advance.

Regards
Shadab
Dear Shadab,
i'm inerested in your model, what solver did you use, is was BuoyantSimpleFoam?
I wanted to simulate a similar case in steady state with solar load only and I got very high temperatures (about 600 K), I thought it was because of the absence of the radiation model but 600 K is still a lot .
can you please give more information about your model, numerical methods and schemes.
could you please share your case file, it will be really helpful for me to do internal simulation of thermal comfort.
regards
Abbes
Abs11 is offline   Reply With Quote

Old   June 23, 2022, 05:38
Default integration of an evaporative air cooler with a refrigeration unit to enhance the ind
  #13
New Member
 
haidermumtaz
Join Date: Mar 2017
Location: najaf
Posts: 5
Rep Power: 7
haideralshami is on a distinguished road
integration of an evaporative air cooler with a refrigeration unit to enhance the indoor thermal environment and production of clean water in iraq
room 3m x4m and height 3m ,walls from hollow block and roof from curpy steel , with door and window in east wall, there is duct inside room two supply one to input cooled air from evaporative air cooler and other from refrigeration unit with return to refrigeration unit, i want simulate to this system with effect sun solar radiation with outdoor conditions on walls and roof daily and 3monthly (in summer) by any program and how to get on results from as indoor and outdoor conditions as temp. and relative humidity
haideralshami is offline   Reply With Quote

Old   September 10, 2022, 06:03
Post Solar Load
  #14
New Member
 
Aswin Chempanchery
Join Date: Jul 2022
Posts: 2
Rep Power: 0
AshC is on a distinguished road
Quote:
Originally Posted by Shadab_Ilahi View Post
Hi,

I have resolved that issue of showing very low temperature on sun direct exposed surfaces.

The flux/power what I am specifying by using externalWallHeatFluxTemperature BC is for convection transfer only, As we don't know on which faces sun radiation will fall and final load will be the sum of convention plus solar loads.

So I slightly changed the externalWallHeatFluxTemperature BC which consider the solar load by looking into solar load file which is generated by running the solar load radiation model. The solar load written on the boundary patches only and for inside the domain, it is zero.

Similarly, I changed the solver also which carry this solar load field with it, actually, once a solar field is generated it is acting as heat source field which is applied only on the boundaries and once the simulation becomes stable run it with radiation model. Because there will be a very high temperature on the boundaries.

Regards,
Shadab





Hello
I am also trying to simulate the same conditions as you.
I have a glass wall(defined as a solid surface) through which solar radiation enters the system and the system heats up which contains air ,2 solid and water flowing inside one of the solid.
I want to know about the boundary conditions you gave on your glass wall with complete explanation.
I am using OpenFoam 10 in which solar load isn't working i guess.
so I have coupled solar load with other radiation models which i read about in one of the forums.
But to my bad luck i am not getting any heat input.(solar load is not working at all).

If you still have your case with you ,sharing the case files can be helpful too.
I can learn it by my own and you wont be needing to waste your time explaining things to me.

Please please help me on this.
AshC is offline   Reply With Quote

Reply

Tags
fvdom, solar load

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



All times are GMT -4. The time now is 21:02.