CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Radiative Heat Transfer

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 17, 2003, 09:44
Default Radiative Heat Transfer
  #1
Manosh
Guest
 
Posts: n/a
Having problem with Discrete Ordinates Mathod (DOM). I want to discuss with someone who are doing the same. Can you please reply me then I will write details.

Regards Manosh
  Reply With Quote

Old   February 17, 2003, 11:05
Default Re: Radiative Heat Transfer
  #2
Rami
Guest
 
Posts: n/a
Hi Manosh,

You are invited to share your problems with me, either here or by e-mail.

Rami
  Reply With Quote

Old   February 18, 2003, 10:21
Default Re: Radiative Heat Transfer
  #3
nuray kayakol
Guest
 
Posts: n/a
What's the problem?

Nuray Kayakol
  Reply With Quote

Old   February 20, 2003, 13:07
Default Re: Radiative Heat Transfer
  #4
Manosh
Guest
 
Posts: n/a
Dear Rami/Nuray

The problem is having 'Negative Intencities'! How do you reduce negative intencities in (DOM)?

Manosh
  Reply With Quote

Old   February 20, 2003, 13:47
Default Re: Radiative Heat Transfer
  #5
nuray kayakol
Guest
 
Posts: n/a
Hi,

This is very well-known problem of DOM when cell-centered differencing scheme is used. It appears when the walls are very hot compared to gas for furnace applications. The common practice is to use a negative intensity fix-up procedure (use of up-wind differencing scheme )for the cells having negative intensity. In your coding you need to change the value of alpha (averaging coeff. which is 0.5 for cell-centered diff. scheme) from 0.5 to 1.0 when negative intensity is encountered.

Here is a simple example; Consider a control volume next to a hot wall. (Intensity values are completely unrealistic for the sake of simple calculation)

I_bw(wall intensity)=11 I_p(cell-centered intensity)=5. I_bw is known from wall Temp. I_p is calculated from recurrence of DOM equation. I_cf(cell-faced intensity leaving the control volume)=?

I_p=(I_bw+I_cf)/2

This equation tells you that the intensity at the centre is the arithmatic average of cell-faced intencities. =(11-1)/2=5 Therofore I_cf=-1 which is unrealistic.

Then accept I_cf=I_p

Read following paper for further information

N. Selšuk and N. Kayakol, Evaluation of Angular Quadrature and Spatial Differencing Schemes for Discrete Ordinates Method in Rectangular Furnaces, Proceeding of 31st National Heat Transfer Conference, ASME HDT-Vol.325, 3, 151-158, Houston, Texas (1996).

Best Regards N. Kayakol
  Reply With Quote

Old   February 23, 2003, 06:54
Default Re: Radiative Heat Transfer
  #6
Rami
Guest
 
Posts: n/a
Dear Manosh,

In the method we developed (Fiterman, A., Ben-Zvi, R., and Kribus, A., DOTS: Pseudo-Time-Stepping solution of the Discrete-Ordinate Equations, Numerical Heat Transfer B 25, 163-183, 1999) you are guaranteed that negative intensities are impossible if non-negative intensities are used as the initial guess. For details and proof please refer to the Appendix of the above reference.

Besides that, this method is very efficient due to its explicitness and use of mult-grid. Recently we demonstrated (in a paper accepted by Numerical Heat Transfer B, but not published yet) it also has high parallel efficiency.

Rami
  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
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Radiative Heat Transfer kamakura117 STAR-CCM+ 0 July 19, 2010 19:51
Convective / Conductive Heat Transfer in Hypersonic flows enigma Main CFD Forum 2 November 1, 2009 23:53
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55
radiative heat transfer Tim FLUENT 0 February 20, 2004 14:28


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