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

CFD Analysis with DO Radiation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 29, 2013, 14:15
Default CFD Analysis with DO Radiation
  #1
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Hi All,

I'm a new member of the forum and an engineering student working with CFD for the first time. I'm new to this so please be patient with me. I work with ANSYS 14.0.

My task is to model a lamp in a duct/box and establish a radiation field with some flow going through this contained space. I've already setup my geometry with inlets and outlets with the lamp in the middle and solved the velocity profiles. My problem is then working with Discrete Ordinate Radiation (DO) emanating from the lamp.

I've tried looking at some manual for basic DO radiation application and followed it quite closely. I found though that my final radiation field looks like a flower (bare with me..). This is very strange as the lamp shape I used is circular so I expected the radiation field to be circular as well. I should probably post some pics.. ** I will do that soon..**

Is there something wrong that I did? Is this a common problem with DO model?

Thanks!
ChaChaLaLa is offline   Reply With Quote

Old   April 29, 2013, 14:38
Default
  #2
Senior Member
 
Join Date: Dec 2011
Location: Madrid, Spain
Posts: 134
Rep Power: 15
michujo is on a distinguished road
Hi. I think this phenomenon is usually referred to in the literature as "ray effect". Try increasing the order of the angular discretization. Could you post a picture of the radiation field and the mesh?

Cheers,
Michujo.
michujo is offline   Reply With Quote

Old   April 29, 2013, 16:18
Default Pictures of Situation
  #3
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Hi. I should have posted pictures with my post. Appologies. Here are some screenshots. Thanks for any help!

Edit: I will look up "ray effect" and try adjusting the angles. Thanks.
Attached Images
File Type: jpg Velocity Profile Vector.jpg (24.7 KB, 38 views)
File Type: jpg Boundary Condition.jpg (14.4 KB, 38 views)
File Type: jpg Boundary Condition 2.jpg (13.9 KB, 36 views)
File Type: jpg Radiation Contour.jpg (13.6 KB, 37 views)
File Type: jpg Radiation Contour 2.jpg (15.5 KB, 32 views)
ChaChaLaLa is offline   Reply With Quote

Old   April 29, 2013, 16:29
Default mesh and geometry
  #4
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Hi I've also got screenshot of my mesh and geometry. Spent a long time on this. I think its good .

Basically, I have a face mapped projected from one side to another with some inflation near edges where things might get tricky. Let me know if you think it can be improved!

Thanks
Attached Images
File Type: jpg Geometry.jpg (14.9 KB, 26 views)
File Type: jpg Mesh1.jpg (35.3 KB, 34 views)
File Type: jpg Mesh2.jpg (43.2 KB, 31 views)
ChaChaLaLa is offline   Reply With Quote

Old   April 29, 2013, 17:36
Default Increased Angle of Discretization
  #5
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Hi, Sorry for the many posts. Last post of the day I promise.

So, I increased the angle of discretization to 8 by 8 (Theta and Phi) and the problem still persists. **defaults were 2 by 2** I tried increasing the pixel as well to 3x3 ** defaults were 1x1** and it made the solution look a bit better but still really off.. (as in its not circular) I am running out of ideas of what to try next

I also can't really find any good information on "ray effect". Could anybody kindly point me in the right direction? Thanks!
ChaChaLaLa is offline   Reply With Quote

Old   April 30, 2013, 10:01
Default
  #6
Senior Member
 
Join Date: Dec 2011
Location: Madrid, Spain
Posts: 134
Rep Power: 15
michujo is on a distinguished road
Hi, I am not an expert on radiation modelling but I think that the ray effect can be mitigated if you come up with a mesh that is somehow aligned with the radiation field. From a radiating cylinder I would expect radiation to propagate along lines of theta=constant.

Could you try building an axisymmetric mesh around your cylinder? Your geometry is not cylindrical so you might not be able to extend it all the way to the walls of your rectangle, but you could do the following:

1) Cut a cylinder off your cube, concentric with your cylindrical lamp. Try getting the cylinder boundary to be as close as possible to the rectangle walls.

2) Put the cylinder and the rest of the cube in the same part in Design Modeler so that you get a conformal mesh.

3) Build an axisymmetric mesh on the cylinder. you can define "Size edges" at the inner and outer boundaries of the cylinder with the same number of elements so that you get mesh lines with theta=constant.

4) Mesh the rest of the rectangle. Hopefully the mesher will build a smooth transition from the cylinder to the rectangle walls, if not try reducing the size of the cylinder.

5) Sweep the face mesh to get the 3D mesh.

I drew a picture of the mesh I'm talking about (it's awful I know).


Also here are some references you might want to have a look at about the ray effect:

http://www.tandfonline.com/doi/abs/1...9#.UX-7kaJFDdU

http://sfera.sollab.eu/downloads/Sch...OM_FVM_MCM.pdf

http://uwspace.uwaterloo.ca/bitstrea...assanzadeh.pdf


Cheers,
Michujo.

P.S: If there's a radiation expert out there please help.
Attached Images
File Type: jpg cylmesh.jpg (39.3 KB, 25 views)
michujo is offline   Reply With Quote

Old   April 30, 2013, 12:44
Default
  #7
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Thanks Michujo. I really appreciate all this help.

I'll go through these papers thoroughly. I had a sinking suspicious feeling inside that the mesh was the problem. I will confirm if that is the case today. Will post update of new mesh if problem was solved. Thanks!

Update:
So, I've taken your advise and created an asymmetric mesh around the cylinder, sweep it, everything's looking good. Solved for velocities no problem. Tried the radiation and this time I got a different shaped "flower" if you will.. **I've worked on many iterations of modifications to try and fix it to no avail** .. ..

I also have some trouble understanding the behavior I'm observing when I change the radiation intensities.. Read the help files and still confused :S

Sometimes FLUENT refuses to solve for incident radiation.. I checked all the absorption coeff and its 0. so.. It isn't being absorbed by anything.. where did all that radiation go??

If there are experts out there, please help! Thanks!
Attached Images
File Type: jpg Square Radiation Contour.jpg (22.7 KB, 24 views)

Last edited by ChaChaLaLa; April 30, 2013 at 17:02.
ChaChaLaLa is offline   Reply With Quote

Old   May 1, 2013, 13:07
Default Success
  #8
New Member
 
Join Date: Apr 2013
Posts: 8
Rep Power: 13
ChaChaLaLa is on a distinguished road
Playing around with the mesh a little more. I managed to get just the profile I was looking for. Hooray! Here is a pic.

Thanks for all the input.
Attached Images
File Type: jpg Success.jpg (40.2 KB, 32 views)
ChaChaLaLa is offline   Reply With Quote

Reply

Tags
ansys 14.0, discrete ordinate, radiation


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
should i use fea or cfd for thermal analysis? mr. bill FLUENT 7 March 3, 2016 15:31
CFD Design...The CFD Future John C. Chien Main CFD Forum 20 November 20, 2015 00:40
Radiation / Convektion / CFD in ANSYS-CFX / Workbench ehrenwirth CFX 8 October 22, 2013 06:12
CFD Analysis of pump remith CFX 7 October 6, 2008 08:57
CFD Symposium (Call for Papers) Chris R. Kleijn Main CFD Forum 0 October 5, 1998 11:25


All times are GMT -4. The time now is 06:53.