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

Heat source in porous zone

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 5, 2009, 12:28
Default Heat source in porous zone
  #1
Member
 
Thomas Wolfanger
Join Date: Mar 2009
Location: South West Germany
Posts: 62
Rep Power: 17
anger is on a distinguished road
Hello Foamers,

I want to model a flow with a porous zone which works well using rhoPorousSimpleFoam in OF 1.6.x.
In a next step, I want to add a heat sink to the porous medium. Is it a good idea to use (and probably extend) the porousZone machinery to this end or is it a better idea to do something on the solver level?
A possible approach is to modify the hEqn with an appropriate source term and multiply it with a field value which corresponds to 1 in the porous region and 0 elsewhere, but this seems not to be very elegant.
What do you think?

Best regards,
-Thomas
anger is offline   Reply With Quote

Old   June 8, 2011, 09:04
Default
  #2
New Member
 
Dima Risch
Join Date: Jun 2011
Location: Cologne
Posts: 22
Rep Power: 14
dima is on a distinguished road
hallo Thomas,
i want as well to add a heat source to the porous zone like you, and want to ask how you handle this?

kind regards,
Dima
dima is offline   Reply With Quote

Old   June 17, 2011, 04:47
Default
  #3
New Member
 
Dima Risch
Join Date: Jun 2011
Location: Cologne
Posts: 22
Rep Power: 14
dima is on a distinguished road
great job guys
looks like what i am looking for is implemented in OF 2.0.0

best regards, Dima
dima is offline   Reply With Quote

Old   June 17, 2011, 10:43
Default
  #4
New Member
 
Gian Maria Di Stefano
Join Date: Jan 2011
Location: University of Applied Science Manno (CH)
Posts: 10
Rep Power: 15
Gian Maria is on a distinguished road
What theory model are you using?

I have the same problem, but I have resolved it applying a modify to hentalpy equation, changing the laplacian therm
Gian Maria is offline   Reply With Quote

Old   June 17, 2011, 11:15
Default
  #5
New Member
 
Dima Risch
Join Date: Jun 2011
Location: Cologne
Posts: 22
Rep Power: 14
dima is on a distinguished road
at the moment i am using just test models, like pitzDaily or a simple inlet outlet channel (tube) with a pZone in it

i was going to resolve it by applying a modify to hEqn as well, changing the source term -> not finished at the moment, just at the beginning

in the new version of OpenFoam there is a thermalPorousZone model (fixedTemperature)

in my further work i will try to add a new model something like fixedEnthaly
but my knowledge of the OF source structure is low at the moment, so i hope there is not a lot of change to do
dima is offline   Reply With Quote

Old   June 22, 2011, 05:56
Default
  #6
New Member
 
Join Date: Mar 2009
Posts: 29
Rep Power: 17
ayoros is on a distinguished road
Hi Dima,

Did you succeed in using the fixedTemperature model for thermalPorousZone ?

When I try it on the rhoPorousMRFSimpleFoam tutorial, it gives me strange behaviour with T (near 0 after the thermal porousZone) and with rho (limited at 1.5 by fvSolution also after the porousZone).

I get this behaviour with different values for the temperature coefficient (even 0).

Could you tell me if you experienced this model ?

Best regards,
Fabien
ayoros is offline   Reply With Quote

Old   June 22, 2011, 08:01
Default
  #7
New Member
 
Dima Risch
Join Date: Jun 2011
Location: Cologne
Posts: 22
Rep Power: 14
dima is on a distinguished road
Hi Fabien,

i didn't succeed in using the angledDuct case too, the same problems occur

so something might be wrong to the calculation (hDiag and hSource) in fixedTemperature.C
or the (rho) settings in the fvSolution ...

but i am not familiar with that at the moment

regards, Dima
dima is offline   Reply With Quote

Old   June 30, 2011, 09:35
Default
  #8
New Member
 
Dima Risch
Join Date: Jun 2011
Location: Cologne
Posts: 22
Rep Power: 14
dima is on a distinguished road
hello,

i changed the fixedTemperature.C

...
forAll(zones, zoneI)
{
const labelList& cells = mesh.cellZones()[zones[zoneI]];

forAll(cells, i)
{
hDiag[cells[i]] += 1/rate*V[cells[i]]*rho[cells[i]];

hSource[cells[i]] += V[cells[i]]*rho[cells[i]]*T_;

}
}

it looks much better now, but not correct at all
would be nice if someone who knows how to handle this could reply

regards, Dima
dima is offline   Reply With Quote

Old   January 26, 2012, 08:26
Default
  #9
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 16
andrea.pasquali is on a distinguished road
Hi Dima,
I'm workin on thermal porous media too.
I checked the OF 2.1 and the fixedTemperature eqn. source is fixed with Cp:

Quote:
...
tmp<volScalarField> Cp = thermo.Cp();
// TODO: generalize for non-fixedTemperature methods
const scalar rate = 1e6;
forAll(zones, zoneI)
{
const labelList& cells = mesh.cellZones()[zones[zoneI]];
forAll(cells, i)
{
hDiag[cells[i]] += rate*V[cells[i]]*rho[cells[i]];
hSource[cells[i]] += 00098 rate*V[cells[i]]*rho[cells[i]]*Cp()[cells[i]]*T_;
}
}
...
I did not check if is correct yet...
I'm trying also to implement non-fixedT, maybe a powerLaw, but I don't know which kind of law I can use.
Do you know something about it?
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   May 4, 2012, 06:38
Default
  #10
New Member
 
Join Date: Mar 2009
Posts: 29
Rep Power: 17
ayoros is on a distinguished road
Hi Andrea,

Did you manage to implement a non-fixedT law in the thermalPorousZone ?

If yes, would you mind sharing it with me ?

I tried by myself but I'm stuck with hDiag and hSource, I don't unserstand how it works with hEqn.

No matter what I try, I got weird results in the temperature field (mostly very low or very high T field).

Thanks,
Fabien
ayoros is offline   Reply With Quote

Old   May 4, 2012, 16:59
Default
  #11
Senior Member
 
Andrea Pasquali
Join Date: Sep 2009
Location: Germany
Posts: 142
Rep Power: 16
andrea.pasquali is on a distinguished road
Hi Fabien,
no I did not implement it... what law did you use for it?
I'm also interesting to add different thermal conductivity (kappa field) just inside the porous medium.
I'm thinking at: kappa_eff = porosity*kappa_porous_medium + (1-porosity)*kappa_fluid.
Do you know something about this?

Thanks

Andrea
__________________
Andrea Pasquali
andrea.pasquali is offline   Reply With Quote

Old   December 16, 2013, 10:49
Default
  #12
Member
 
nadine moussa
Join Date: Mar 2012
Posts: 30
Rep Power: 14
nadine is on a distinguished road
Hello Andrea,

did you managed to calculate an effective conductivity as function of the porosity?

I am so intrested in this subject!

thanks a lot,
Nadine
nadine is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
species mass source in porous media ? PK FLUENT 0 February 16, 2007 12:12
porous medium as heat source Ralf Schmidt FLUENT 1 November 4, 2005 14:12
UDF Scalar Code: HT 1 Greg Perkins FLUENT 8 October 20, 2000 13:40
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 14, 2000 00:03
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 04:43


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