CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

radiation modeling in vacuum

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2011, 10:22
Default radiation modeling in vacuum
  #1
Senior Member
 
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 8
mirko is on a distinguished road
As famously said, in space no-one can hear your scream. But seriously ...

The short question is:

how are field values on patches updated? For example, in laplacianFoam, it is obvious when the fields are calculated. But it is not clear when the patch boundary conditions are applied.

I need to apply patch boundary conditions without calculating the fields.

The long of the story follows:

I need to do a simulation of radiative heat transfer between solids in vacuum.

I have hacked the existing chtMultiRegion solver to create a much simpler solver that involves radiation and conduction only. However, just as vacuum does not carry sound, it does not conduct heat.

I can set the head conductivity to a very low number, but the solver will still solve for the temperature field. And for very large meshes, the CPU cost will add up. I would like to avoid that.

So, I need to model a region with radiation only.

I looked at regionProperties.H/C, and it seems straightforward to extend that code to include a region property of type vacuum.

The next item on the agenda is to introduce boundary conditions, steady state for now: solidToVacuum and vacuumToSolid. I think I can manage that.

My main question is the calculation of radiation flux in the vacuum region. From studying existing solvers (see http://foam.sourceforge.net/docs/cpp/a02314_source.html), the procedure is to
  1. calculate the temperature field in the fluid region.
  2. That updates the temperature on the patches (this happens behind the scenes)
  3. The patch temperature is then used as a boundary condition on the radiation flux.

For the vacuum solver, I would need to skip step 1. I need to update the patch temperature and solve for radiation.

But I don't know how to trigger the calls to the boundary condition updates.

Thanks,

Mirko
mirko is offline   Reply With Quote

Old   September 9, 2011, 12:41
Default
  #2
Senior Member
 
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 397
Rep Power: 12
marupio is on a distinguished road
I'm not entirely sure what exactly would suit your needs... so here's what comes to mind:

fvPatchField::evaluate() looks like it updates the boundary field values.
fvPatchField::initEvaluate() looks related, but I think it has to do with parallel.
fvPatchField::updateCoeffs() looks like it updates matrix coefficients, intended to be used after evaluate(). Do you use an fvMatrix in your solution?

If you are using evaluate(), I'm not clear on how it will work... boundary values typically depend on the neighbouring cell in the internal field. Not sure how that works with your vacuum... but I don't deal with vacuums, so I defer to you.

If you are using an fvMatrix, and need the matrix coefficients to "absorb" the boundary conditions, that's a little more complicated. Boundary conditions affecting the diagonal of the matrix are updated with addBoundaryDiag - see fvMatrixSolve.C. Boundary conditions that require out-of-core multiplication, such as cyclic or processor, these require different treatment by different matrix solvers, and therefore, rather than trying to generically absorb them into the matrix, we have to build a set of boundary and internal coeffs, one for each boundary cell... and then let the solver handle them. Again, see fvMatrixSolve.C.

Hope that helps.

-Dave

Follow me on Twitter @DavidGaden
marupio is offline   Reply With Quote

Old   September 12, 2011, 17:34
Default
  #3
Senior Member
 
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 8
mirko is on a distinguished road
Thanks,

I'll look into it

Mirko
mirko is offline   Reply With Quote

Old   November 16, 2011, 09:32
Default laplacianFoam with radiation ?
  #4
New Member
 
Mitja Majerle
Join Date: Nov 2011
Posts: 2
Rep Power: 0
ssdmmaje is on a distinguished road
I need something similar, a simple laplacian solver (I added extra heat deposition) with radiation capabilities.
I am simultaneously working on chtMultiRegionFoam (it is one week since I am trying to understand how it works - is there any manual on this example ?), and I do not know what will be simpler, to get to understand chtMultiRegionFoam and add extra heat deposition or to somehow include radiation in laplacianFoam (including extraheat was really easy, but radiation is quite a problem). Any recommendations ?

Thank you,
M
ssdmmaje is offline   Reply With Quote

Old   July 2, 2015, 10:20
Default
  #5
New Member
 
Join Date: Mar 2015
Posts: 11
Rep Power: 2
coccone is on a distinguished road
Did you succeed in modifying chtMultiRegion?

I have a similar problem and am very interested in the subject, but I am not so skilled to get hands in the solver.

Any suggestions?
__________________
Use the log-force Luke!
coccone is offline   Reply With Quote

Old   July 2, 2015, 11:21
Default
  #6
New Member
 
Mitja Majerle
Join Date: Nov 2011
Posts: 2
Rep Power: 0
ssdmmaje is on a distinguished road
No, I implemented radiation into laplacianFoam.
Would this help you ?
ssdmmaje is offline   Reply With Quote

Old   July 3, 2015, 05:42
Default
  #7
New Member
 
Join Date: Mar 2015
Posts: 11
Rep Power: 2
coccone is on a distinguished road
Quote:
Originally Posted by ssdmmaje View Post
No, I implemented radiation into laplacianFoam.
Would this help you ?
Unfortunately not, I am bound to cht for different region properties... but thank you though
__________________
Use the log-force Luke!
coccone is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
vacuum for radiation purpose tianchen FLUENT 0 May 24, 2011 16:56
Radiation modeling cloudnqh FLUENT 6 February 15, 2010 10:06
Radiation in vacuum nitin CFX 0 July 6, 2009 01:49
Modeling Thermal Radiation Kaushik Saha FLUENT 0 September 17, 2008 19:12
Question about Radiation Modeling Zhengcai Ye FLUENT 0 October 12, 2004 18:44


All times are GMT -4. The time now is 17:56.