CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLUENT (https://www.cfd-online.com/Forums/fluent/)
-   -   Fluent Radiation/porous media (https://www.cfd-online.com/Forums/fluent/48421-fluent-radiation-porous-media.html)

Schmitt pierre-Louis June 5, 2008 10:13

Fluent Radiation/porous media
 
Hi,

I'm now facing a big problem with fluent. I have to model a high temperature device wich comprises a porous media and a high temperature fluid zone. The problem is that by default a porous zone is considered as a fluid zone for the radiation and the radiation can pass trough the porous medium. It's not possible to consider that one side of an internal face participate to radiation and the other doesn't since it makes a bug during the computation and I can't define this interface as a wall since it will block the flow. There is also a problem using UDF since a UDF for special radaition can't be attached to any boudary condition exept a wall. Do you see a way to overpass this problem, or is it impossible. Thank you Pierre-Louis Schmitt


Rami June 10, 2008 04:26

Re: Fluent Radiation/porous media
 
Hello Pierre-Louis Schmitt,

I work in research and development of volumetric solar receivers. You are facing a problem which is also highly relevant for me. In fact, I tried to get answers from Fluent in these regards for a very long time (few years), but they never answered it satisfactorily.

In the past we have used another commercial CFD package and augmented it with a lot of user coding to treat radiation in a participating medium, inter-surface radiation, porous medium treatment for HT and resistance to flow, chemical kinetics and more. In these models, it was assumed that in the volume occupied by the porous medium two distinct temperatures prevail - one for the fluid and the other for the solid. It was also assumed that the fluid is non-participating in the radiation (which is usually a good approximation), and that the porous medium is represented by its equivalent effective radiative properties. The relevant temperature for this medium emission, according the the above assumptions, is the solid temperature of each cell. This model worked fine as long as the other CFD package (no names mentioned...) did not fail. However it failed whenever the geometry was not trivial and in many other cases related to the problem physics, so we had to abandon it.

As far as I understand, there is no straightforward treatment of something similar in Fluent, as it assumes a SINGLE temperature in a cell, which is the FLUID temperature in the case of a porous medium. As this choice is limited and inappropriate to our needs (and probably to any body's dealing with radiation in a porous medium) it seems to call for a massive UDF coding. Being relatively novice in Fluent, I have not started it yet, and it may take a long time to get there.

As I started, I hoped this will be solved by Fluent, but they don't seem to tend to do it and/or did not understand the need for it.

Sorry I am not helping you here, just sharing my point of view. Maybe someone else (a more experienced user, and hopefully some expert from Fluent support or development teams) will have a constructive solution (or say what I claim is wrong, and this is already solved...).

Rami


Rami June 17, 2008 03:41

Re: Fluent Radiation/porous media
 
Unfortunately, it seems nobody responds. It is regretful that Fluent employees are not responding either... Let's hope someone will soon.

Rami June 24, 2008 01:48

Re: Fluent Radiation/porous media
 
Pierre-Louis Schmitt,

Did you get any response/help on this issue? I would be very grateful if you can share with me and this forum.

Thanks, Rami

Schmitt pierre-Louis June 24, 2008 02:34

Re: Fluent Radiation/porous media
 
Unfortunately I didn't get any response on this problem, I found a kind of "technic" to overpass the problem for simple cases. on a first step you compute a solution considering the porous zone as a solid zone with a thermal conduction equivalent to that of the porous medium, fluid zone is kept normally. Then you print the profile of temperature, reload the case with the porous zone and then impose the profile of temperature on fluid zones.

I think this hypothesis is good when the velocity of the flow in the porous media is low. It's the only way I found to get a simulation of my device.

Pierre-Louis Schmitt

Rami June 25, 2008 02:05

Re: Fluent Radiation/porous media
 
Thank you. This workaround may be useful in your case, but not in more general like mine. Hopefully someone (FLUENT stuff, please respond!) will suggest a comprehensive solution, since the need to consider porous media participating in radiation is quite common in many applications.

Rami March 19, 2009 04:46

Re: Fluent Radiation/porous media
 
Quote:

Originally Posted by Schmitt pierre-Louis
;151747
Unfortunately I didn't get any response on this problem, I found a kind of "technic" to overpass the problem for simple cases. on a first step you compute a solution considering the porous zone as a solid zone with a thermal conduction equivalent to that of the porous medium, fluid zone is kept normally. Then you print the profile of temperature, reload the case with the porous zone and then impose the profile of temperature on fluid zones.

I think this hypothesis is good when the velocity of the flow in the porous media is low. It's the only way I found to get a simulation of my device.

Pierre-Louis Schmitt

I am still trying to have FLUENT development/support constructive response to this issue. As of now, they seem to treat this as a low priority issue in the development, since they consider this as something of little interest to other users and/or merely academic problem.

Are there any others (except me and Pierre-Louis Schmitt) who need this combination of radiative heat transfer in porous media and suffer from this apparent FLUENT deficiency to its treatment?

Rami March 23, 2009 02:54

I have just noted this issue was already discussed for a long time. Here is one thread, dated 2001:
http://www.cfd-online.com/Forums/flu...ous-media.html

So it seems that there is interest in radiative HT in porous media. Any others interested, so we may convince FLUENT to finally do it right?

Rami May 6, 2009 04:56

No one? Am I the only one in need for such features?

waltastisch January 28, 2011 10:43

not alone
 
No you're not alone with this problem. I face the same problem in a high temperature furnace in which a porous body is sintered. It still seems impossible to have radiative heat transfer to the porous zone as well as fluid flow between the porous zone and its surroundings in FLUENT :mad:

Rami January 30, 2011 05:31

It's funny. I have already gave up using Fluent for these kind of problems. Recently, a visiting student at our group has told me he had received a UDF for this kind of problems from Fluent (after long discussions with them). It is called "dual cell method", used in a tutorial called "Light-Off simulation for Catalytic Converter" dated April 1, 2005 by Fluent Inc., and in principle uses a (very inefficient and not elegant) method, in which the mesh is doubled, part for the solid and part for the fluid, and the UDF couples both to enable treating radiation (if I correctly understood). I did not check it myself - since we no longer hold a Fluent license, as it seemed inadequate for our needs - but was told it does the trick (though inefficiently). I just wonder why the Fluent guys did not suggest it to me when it was rellevant...

Maybe you can ask them about it, having all the details above.

waltastisch January 31, 2011 02:50

I was offered this UDF as well, but since we have rather large and complex models, we decided that this is not the right choice.

To our luck, our problem is thermally uncoupled. The gas flow and the reactions within the porous medium are locally chemically important, but of neglible impact to the global temperature or flow field. Therefore we calculate the problem twice with the same model: first time the porous zone is just a solid, thermal radiation is calculated OK this way. The second time we calculate the same problem again, but the temperatures are taken as fixed values from the first run. Therefore radiation does not need to be solved and the porous medium formulation of FLUENT can be used.

This seems all quite OK for our special problem setup, but it is still strange that FLUENT does not provide a standard method to handle this...

Rami January 31, 2011 05:49

Hi Walter,

Out of curiosity: where are you working and what are you specializing in (if you wish to tell, of course)?

The reason I am asking is also due to the fact that we tried very hard to get a solution for our problem (a volumetric solar receiver) with Fluent. This was after we had developed all the required capabilities (radiation and chemical kinetics in a general porous medium) for another package (PHOENICS), but had a lot of problems with PHOENICS itself (mainly lack of generality in the geometry and grid and too many bugs). It took ages to get replies from Fluent, and they never mentioned this dual cell UDF, and in general - seemed like doing us a favor answering at all. I suspect it was since we are in academia, holding a single seat academic license. Your answer may help support my suspicion.

At the moment we are looking again for a suitable package, and I think we will give Fluent another shot, but will also try others, with the hope to get more collaborative and constructive attitude.

Rami

waltastisch January 31, 2011 07:20

Hi Rani

We're commercial users (only single license) and we use our models to calculate global models of high temperature furnaces. Within the oven we have porous solids which show some slow chemical reactions with the oven gas - we didnt' care about these reactions up to now. We want to extend our oven models to see the evolution of these chemical reactions as well. At the current stage we're just planning these steps. It'll take some months to do the programming and see how it works. I'm rather confident we will be able to do it in the way mentioned. To me, the Fluent support seems sometimes a bit slow and sometimes you have to remind them to take care of the question - but its generally OK.

pilou January 31, 2011 07:49

Hi, I finally found a way to overpass this problem

I made a UDF that I copy past here :

DEFINE_PROPERTY(abs_coeff,c,t)
{
#if RP_HOST
real abs;

{
if (THREAD_ID(t) == valeur1)
abs = absp;
else
abs = 0.5;
}

return abs;

#endif

}

This UDF basically attach a different coefficient on the different zones depending on the zone ID. This way, the radiation in the porous medium is not perfectly equal to 0 but is very low and it gives a good approximation of the porous medium behaviour regarding the radiation.

This UDF is ready to be used in parallel processing.
Hope it can help

Pierre-Louis Schmitt (Pilou)

oky February 1, 2011 20:54

Hi everyone,

I'm trying simulation porous media in rectangular channel, but the result isn't suitable with any research.

So, would you help me. I wish someone can check my simulation and give some reports if there is something wrong.

Thank you for your help.
Please send your e-mail, than i will send you my works to to your email.
my email: oky.andytya.net@gmail.com

Regrads,
OKY Andytya P

note:
I use ANSYS Fluent 6.3 [CFD]

pilou February 2, 2011 03:58

Can you give us a few more information about your problem, type of porous medium, are you using a radiation model, what are the values of viscous and inertial resistances you are using...
Maybe we can help without the .cas file.

Regards

mystic_cfd February 3, 2011 02:42

Fluent Radiation/porous media
 
i am also one of the disappointed users of Fluent with regards to porous zone modeling.

if you still have to decide on a code, have look at CFX. In the latest release (13), they claim the following:
"In ANSYS CFX technology, porous CHT objects can now be modeled with separate fluid and solid temperatures. A user can specify the interfacial area density between solid and fluid together with a heat transfer coefficient. Energy is then conducted through the solid based on the solid properties and exchanged with the fluid."

Since CFX and Fluent are set to merge, we only have to wait another 3 or 4 years to have it in Fluent as well... :p

oky February 17, 2011 21:42

For anyone, especially Mr. Pilou

Thank you for helping, i'm a student. I want get the flow and heat transfer characteristic in porous medium.

I have a channel with two baffles porous medium.
From the bottom of channel, i give heater.
The porosity is 15%.
No Radiation.
I'm trying to simulation in Laminar and Turbulent condition.

Hěr0 July 13, 2011 10:35

Quote:

Originally Posted by pilou (Post 292946)
Hi, I finally found a way to overpass this problem

I made a UDF that I copy past here :

DEFINE_PROPERTY(abs_coeff,c,t)
{
#if RP_HOST
real abs;

{
if (THREAD_ID(t) == valeur1)
abs = absp;
else
abs = 0.5;
}

return abs;

#endif

}

This UDF basically attach a different coefficient on the different zones depending on the zone ID. This way, the radiation in the porous medium is not perfectly equal to 0 but is very low and it gives a good approximation of the porous medium behaviour regarding the radiation.

This UDF is ready to be used in parallel processing.
Hope it can help

Pierre-Louis Schmitt (Pilou)

Hi Pilou, can you clarify why you want low absorption coefficient in the porous medium, please? I have the same problem simulating a flow passes through the stainless steel. By default partecipates in radiation is turned on for the porous zone and I know it is considered as a fluid zone. If I want the porous medium radiates as well as the wall of the furnace I think I have to increase the absorption of the porous zone to simulate the increase in heat transfert between fluid and solid due to the radiative contribution. I set high absorption in the SS material properties but the contour of absorption shows only the fluid absorption. Is there a way to pass to fluent the abs_eff as done for k_eff. The same for emissivity of the porous solid.

Another question....how can you set the emissivity of the fluid? I can find any option.

Thank you so much! Let me know, please!

Luigi

pprtyy November 26, 2011 18:02

Hej,

your comments are good. But CFX might be a solution, since it has the function to analysis conjunction of fluid and solid in porous material

fredom May 15, 2013 05:20

Hello everybody,

I am interested for this topic, because I would like to simulate heat transfer between hot wall 1 (473K) to wall 2 which is perforated sheet. The distance between them is 200mm.

The purpose is to define the temperature of the perforated sheet. As the velocity inlet is low 2m/s and it pass paralle to the perforated sheet, should I define the perforated sheet as a porous jump ?

Thanks very much for your comments..

Fredom

Fluent_user_2015 December 15, 2015 13:28

Implementation of radiation in porouse media in fluent
 
Hello everybody,

i am facing the same problem as described above. in want to simulate combustion iside a porouse zone in fluent. Radiation is also very important in my case. Even in the newes version of fluent (16) it seems not to be possible to consider radiation from porous media in porous zones. :confused:

Has anybody a solution for this problem?

Thanks very much for your help!

Rami December 16, 2015 02:57

Dear Fluent_user_2015,



I recently joined another group (hydrology research) and do not use Fluent any more, nor radiation and thermal issues in general. However, I' try to help.


Please find below the (long) correspondence I had with Fluent support that addresses using radiation solvers with thermal non-equilibrium in porous medium. I hope it is helpful.


The trick is described in the 29/01/15 message from SKS:
(rpsetvar 'rad/enable-netm? #t)
and further discussed in following messages. You have plenty of information and I trust it should suffice.

Best luck,
Rami

From: Rami
Sent: Thursday, February 06, 2014 14:42
To: 'SKS'
Subject: RE: SR Owner Notification SR#: ---has been updated


Hi SKS,
This is what I suspected…
Thank you for your efforts,
Rami

From: SKS
Sent: Thursday, February 06, 2014 14:19
To: Rami
Subject: Re: SR Owner Notification SR#: ---has been updated


Hi Rami,

I can reproduce the problem. I had a discussion with Development and found that currently P1 radiation with NETM will not account for solid!!. So, I would recommend, if you want to use P1 radiation, please use the UDF you developed. You may use in-built DO radiation with NETM, as programming DO radiation is not an easy task. Since development is creating compatibility with one radiation model after another, it will take some time for them to make it fully available. Please bare with us.

Thanks
SKS

On Thu, Feb 6, 2014 at 4:25 PM, Rami wrote:
Hi SKS,
The files of my Feb 4 2014 message are attached again. It would be nice (and less error prone!) if it can run using the same UDF, input files (cas and inp.txt) and settings (except P1 or DOM models).
I am grateful for your dedication!
Rami
From: SKS
Sent: Thursday, February 06, 2014 12:47

To: Rami
Subject: Re: SR Owner Notification SR#: ---has been updated
Hi Rami,
Sorry that initially I understood your question very differently. That is the reason my answer did not make sense now.
Now, I tried to run the case file you sent me on Feb 4 2014, I am not able to initialize the case. It seems that the case file is corrupted. Could you please check and re-upload the files?.
By the way, you did not do anything wrong. This kind of problem occurs when it is in the development stage. As FLUENT is getting bigger and bigger, we need to make sure that any piece of new code must work with all combinations and we must validate it. Once I get the correct case file, I will investigate and let you know the workaround if there is any
Thanks
SKS

On Thu, Feb 6, 2014 at 2:48 PM, Rami wrote:
Hi SKS,
Thank you for answering, but this time I do not understand your reply.
As you may see in the files I had sent you on Feb 4, 2014,
·I had enabled NETM + radiation before setting up the case.
·The NETM option was used and created the solid zone.
·I had enabled radiation in the solid zone (and also in the fluid, with fully transparent properties).
·I initially chose DOM for radiation.
·I prescribed BCs (specifically, collimated radiation using the Qirr DEFINE_PROFILE macro from the on the front surface, inlet:012)
·It all ran fine.
·I then only changed the radiation model from DOM to P1 (without any other change in the case file and in the udf), and then the radiation properties disappeared from the solid material (foam), and the radiation BC settings disappeared from all solid surfaces.
Is there anything I did wrong?
Thanks,
Rami
From: SKS
Sent: Thursday, February 06, 2014 04:40

To: Rami
Subject: Re: SR Owner Notification SR#: ---has been updated
Hi Rami,
Did you enable Participation of radiation for solid in cell-zone condition panel?
Since NETM + Radiation is in the process of development, I would recommend you to use the UDF which is specifically written for your purpose, only if enabling "Participation of Radiation" did not help. NETM + Radiation feature will be available only after the compatibility check with all other models and the accuracy of the code is tested extensively. Hopefully, it will be available in R16.
Thanks
SKS
On Tue, Feb 4, 2014 at 7:31 PM, Rami wrote:
Dear SKS,
Thank you so much for your wonderful and prompt help! Everything worked as you suggested.
Now, I still wish to use the built-in P1 model with the scheme command trick to allow using it with a non-equilibrium thermal model (NETM) in order to compare it to experiment described in Sec. 4.2 and Figs. 10-11 of the original reference paper (attached again for your convenience). I imagined that with the new UDF, cas and input files (attached) I need to only switch between the DOM and P1 and that’s it. However, it seems that the solid has neither access to the absorption and scattering coefficients nor to the BCs (the Marshak BCs at walls and possibly the collimated irradiation in the front).
I tried also (in the attached cas file) to put the absorption and scattering coefficients of the solid in the fluid material slots, but the temperature is almost uniform, since the irradiation is absent.
Is there a way to do it, or should I stick with the former udf, modified from the UDF Guide example?
BTW, I have found that I may read the input file in a DEFINE_INIT rather than DEFINE_ON_DEMAND, which is more convenient.
Thanks a lot for all you did thus far and also in advance!
Rami
From: SKS
Sent: Tuesday, February 04, 2014 11:49

To: Rami
Subject: Re: SR Owner Notification SR#: ---has been updated
Hi Rami,
Glad to know that you are happy with the scheme command though it has some issues.
Yes. This kind of problem is expected, as this feature is not fully developed. That’s why I mentioned in my last mail that, this has not been tested extensively. This is the same reason why this feature has been hidden. We provide the scheme command as a workaround for the users who explicitly ask for the combining of NETM and radiation. Developers are working to make this feature available to all users (without the scheme command). It may be available at R16. In R15 we do not have any option. So, please follow the same order to setup the problem.
1. As a workaround for the error, when you try to replace the mesh switch to equilibrium thermal model -> replace mesh -> enable NETM. This will avoid the error. But, this is just a work around. When this model is available in later version, you will be able to use it in the same way you use for the other features. When you are using workbench you should change to equilibrium thermal model before you save the project. Please let me know, if this works
2. I looked at the UDF. It is better to use DEFINE_ON_DEMAND, as you are going to read the file only once in the whole simulation. The problem is that you are defined the variables are real. So, please convert the data types from real to float. It will read the data properly. I have tested it. It does read properly.
Thanks
SKS
On Mon, Feb 3, 2014 at 1:54 PM, Rami wrote:
Hi SKS,
I managed to use your suggestion and it seems to work fine with the DOM, but it is indeed tricky to use, and I have tweaked a lot to find the right order to set up a case using standalone FLUENT 15.0:
1.Start FLUENT.
2.Read the mesh file created by WB.
3.Issue the TUI command (rpsetvar 'rad/enable-netm? #t) .
4.Setup the case, including solving energy and DOM, then define the porous domain as non-equilibrium.
5.Add the UDF, build and load.
Doing it in a different order caused errors of various types, e.g.,
Error: received a fatal signal (Segmentation fault).
Error Object: #f
The case and the UDF files I used are attached.
Now I have two related questions:
1.Suppose I wish to replace the mesh with another mesh (e.g., to assess the grid independence), and the new mesh is, say, the attached FFF.msh . Just reading and replacing the existing case grid does not work, since the dual cell method was applied to enable the non-equilibrium porous zone. The new mesh zones and BCs are therefore yielding errors. How do I do it without having to go through the entire setting process again (with a good chance of missing part of the definitions)?
2.I would like to add reading input from a text file rather than putting them “hard-wired” in the UDF, but failed– all the parameters remained 0 (see the commented lines in the DEFINE_INIT macro within the UDF and the attached inp.txt as an example). The advice from my local agent (using a DEFINE_ON_DEMAND macro also failed). What have I done wrong?
Thank you in advance,
Rami
From: Rami
Sent: Wednesday, January 29, 2014 16:51
To: 'SKS'
Subject: RE: SR Owner Notification SR#: ---has been updated

Hi SKS,
I am still checking and have not fully set the case, but – it seems to work (with DOM)! I will let you know when I reach some conclusions.
Thank you so much!
Rami
From: SKS
Sent: Wednesday, January 29, 2014 13:07
To: Rami
Subject: Re: SR Owner Notification SR#: ---has been updated

Hi Rami,
I spoke internally as well as Radiation development manager and radiation team members. There is a rpvar to enable radiation model with Non-equilibrium thermal model. Before setup the case, type
(rpsetvar 'rad/enable-netm? #t)
in TUI. Then you will be able to enable Non-equilibrium thermal model with Radiation.
NOTE:
  • It is available only in R15.0
  • This has not been tested extensively. That is the reason, it is hidden.
  • IMPORTANT: When you enable Non-equilibrium thermal model and close the cell-zone panel and reopen the same panel, you will notice that, FLUENT switches the Thermal Model to Equilibrium from Non-equilibrium. This is an artifact. I have reported this to development and they will fix this in R16. But, when you enable Non-equilibrium model and just click OK, solver consider Non-equilibrium thermal effect with Radiation. So, make sure that, each time you close cell-zone panel Non-equilibrium thermal effect is enabled.
Please let me know, how it went
Thanks
SKS
On Tue, Jan 28, 2014 at 3:09 PM, Rami wrote:
Hi SKS,
Thank you for your prompt reply! As I have already explained, the P1 is a severe compromise in the problems we wish to cope with, so I will prefer to use it only as the last line.
In the following I will add my response to each point you made.
If you would like to use FLUENT inbuilt Non-equilibrium thermal model and implement DO radiation through UDF, it is tricky. I am checking internally, if it has been done by anybody in ANSYS. I will get back to you, once I get reply.
Please check and let me know.
We can also try implementing energy equation for solid through UDS and use inbuilt radiation model. I haven't tested it or I cannot find any work done by anybody. But, by this way, you will be able to use DO radiation model without any problem.
This is unclear to me: I am having flow in a porous medium (foam), so I need to solve a fluid in a porous medium (mass equation, resistance terms and convection to the solid – all influenced by the porosity) + a solid energy equation (including radiation and convection terms). Please clarify how your suggestion copes with these while using a DOM for the radiation.
To answer your question regarding example in UDF manual, it should be Thread *t0 = THREAD_T0(thread);. Because, t0 is used as cell thread in later lines.
Being novice in C and writing UDFs, I tried your correction rather than mine in a case I ran, and got the same results. Thank you!
Best regards,
Rami
From: SKS
Sent: Monday, January 27, 2014 13:50
To: Rami
Subject: Fwd: SR Owner Notification SR#: ---has been updated

Hi Rami,
Great to know that you have implemented P1 radiation using UDF with the help of example in UDF manual!!.
If you would like to use FLUENT inbuilt Non-equilibrium thermal model and implement DO radiation through UDF, it is tricky. I am checking internally, if it has been done by anybody in ANSYS. I will get back to you, once I get reply.
We can also try implementing energy equation for solid through UDS and use inbuilt radiation model. I haven't tested it or I cannot find any work done by anybody. But, by this way, you will be able to use DO radiation model without any problem.
To answer your question regarding example in UDF manual, it should be Thread *t0 = THREAD_T0(thread);. Because, t0 is used as cell thread in later lines.
Thanks
SKS
/************************************************** *****************/

Technical Details:

Activity Description: I am actually looking for a better radiation solver (such as the DOM) rather than P1

Activity Detail: Dear SKS,

Thank you for your response. I am already aware of all you wrote, but this was not my query.

I need a better radiation solver to cope with the specular reflections and spectral effects of the quartz window. I was therefore asking if the P1 implementation example in the UDF guide (1) has a better option, i.e., is the limitation removed so that all the radiation solvers (especially the DOM) are already available in combination with the non-equilibrium porous media. I understand it is still not so (what a disappointment...).

I am already using a modified P1 UDF including all the details in the papers attached by Ed Moses on my behalf, but this is a very crude model, which does not allow for the spectral effect of the quartz (which is very strong, as it is nearly fully transparent below ~3.5um, and nearly black above it) and its specular reflections. The P1 model will also limit me in prescribing the incident solar irradiation to a more realistic case in which it is neither diffuse nor collimated. It is also known to be a crude and non-conservative approximation (also observed by Wu in the paper sent to you, both in the need to modify the BCs and in the comparison to the experiment).

If you can suggest a better approach than P1 (and not writing a full DOM UDF (2)... ) that overcomes these limitations, please let me know. I would be grateful if you could send it directly to Rami
Regards,
Rami

(1) PS: I suspect the P1 implementation example in the UDF guide has an error in DEFINE_PROFILE(p1_bc, thread, position), and the line "Thread *t0=thread-t0;" should be "Thread *t0=thread;" - is that correct?

(2) BTW, we had written a standalone DOM Fortran solver and implemented it in another CFD package long ago.

Thank you.

User12 May 25, 2016 07:03

Hello,

I am also trying to model radiation and porous media in fluent. Reading this thread helped me a lot to understand the limitiations of fluent regarding this matter.
My model is a sort of catalytic converter where the catalyst zone (porous media) is a metallic honeycomb at elevated temperature. Upstream this zone, the metal casing forms a really sharp cone, so this wall is very close and exposed to the front face of the porous catalyst and could get hot due to radiation from it. My objective is to model how much this casing is going to heat up. I am far from being an expert in radiation but, in my opinion, the front face of the porous media will behave as a solid surface in terms of radiation. Since fluent treats porous media as a fluid, I cannot define an emissivity for this face. Then, I tried to set a non-equilibrium thermal model in the porous media by following Rami's solution (previous post). Thus, a solid zone is created and the face emissivity can be defined. However, this wall appears as an external wall, or at least the boundary conditions that can be set are those corresponding to an external wall, i.e, heat flux with external emissivity, convection/radiation to a free media, and so on. I performed many tests and my conclusion is that the solid body is not interacting with its surroundings. Irrespective of the boundary conditions I try, the solution does not change. My impression is that this fictitious solid which is created is only connected to the corresponding fluid zone of the porous media in order to calculate the different temperature profiles (non-equilibrium in the porous media). Therefore, for my case, Rami's solution is not working. I spoke with fluent support but they cannot give me te solution for this problem.

Searching for an alternative, I built a very simple model in which I use DO and only solve for radiation equations. I found out that when I set a very high absorption coefficient in a gas zone, this behaves in a very similar way as a solid emitting at the same temperature, so maybe if I set a high absorption coefficient in my porous media, I could capture radiation of the front face as if it was solid. I was wondering if this trick could work out for my case, but I am not sure if I will negatively affect the model in other aspects. Can anybody advise on this? Does this make physical sense?

Thanks

mmunige July 23, 2016 05:12

Modelling heat flow in porous media
 
Dear all,

I am modelling heat transfer through a porous medium in fluent (from hot air to porous medium), have defined porous medium in fluid cell zone, resistance coefficients,
porosity, heat transfer coefficients since i am using non-equilibrium thermal model. this is transient 2D simulation. (a vertical column having porous material inside and top inlet for hot air and bottom outlet).

in start temperature distribution along the domain is fine and close to experimental values, but after 300 or 400 seconds temperature does not change or changes very little at all sections of domain and further it goes below the experimental values. I m using both fluid and solid properties as temperature dependent.
please guide me where i am doing wrong so I am not getting the uniform trend of temperature variation along the domain in whole simulation period.
please

Regards

iam September 1, 2016 10:29

Question Modelling heat flow in porous media
 
Dear Ramni,

Thanks a lot for posting your correspondence with Ansys support, it was really helpful to activate the non-thermal equilibrium option.

I am working at the moment in a surface burner which is conformed of a ceramic fibrous porous mat. How could I estimate the Interfacial area density, and the absorpion coefficient for the porous material?

Thanks a lot for your help!!


All times are GMT -4. The time now is 14:54.