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

Duct structure in Wind Tunnel: Continuity problem?

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

Like Tree2Likes
  • 1 Post By JR22
  • 1 Post By JR22

Reply
 
LinkBack Thread Tools Display Modes
Old   March 31, 2013, 16:20
Default Duct structure in Wind Tunnel: Continuity problem?
  #1
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 125
Rep Power: 8
JR22 will become famous soon enough
Hi Foamers,

I am trying to model a duct on a pole within a wind tunnel. I am using helyxOS initially for meshing and problem setup as I develop the skills to correctly use OpenFoam. I am using simpleFoam, SST-kappa-omega turbulence model, and about 8000 iterations. I am trying to measure the flow coming into and out of the duct, and at the center of it.

The duct is composed of a nozzle, a center section, and a diffuser (See top of figure). Using Paraview, I placed 3 planes, one at the front, one at the center, and one at the back. I then use paraview to resample to source (project onto the surfaces, called sources in paraview), then I integrate on them (See bottom of figure).

I took the integrated values for the velocity vector and compared them for several mesh qualities. Towards the higher quality meshing (bottom of the table), results seem to not vary too much. However, the velocity vector in the front is NOT the same as the velocity vector at the back. Does this make sense?

Thanks in advance.



Attached Images
File Type: jpg duct_within_wind_tunnel.jpg (67.5 KB, 62 views)
File Type: png duct_within_wind_tunne_results.png (25.0 KB, 60 views)
immortality likes this.

Last edited by JR22; March 31, 2013 at 19:02. Reason: added more details
JR22 is offline   Reply With Quote

Old   April 1, 2013, 08:09
Default Appropriate mesh size for the SST-kappa-omega
  #2
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 125
Rep Power: 8
JR22 will become famous soon enough
As I increase the mesh size in the inner portion of the duct, the picture starts to look clearer. Here is an updated image of the 3 surfaces in the duct + the data from the last run (ran overnight for 7 hours). The back surface looks like the image the commissioner uses to call batman, but any resemblance to reality is mere coincidence. The updated table here shows that the integration values are changing and starting to look more like reality. I was expecting more acceleration at the center of the nozzle-diffuser.

Any advice on mesh refinement?
Is there something wrong with my model besides mesh size?
Is there a more appropriate way to integrate (putting a patch that allows flow through it? -- is there something like that?)



Attached Images
File Type: jpg duct_within_wind_tunnel_2.jpg (61.1 KB, 55 views)
File Type: png duct_within_wind_tunne_results_2.png (26.8 KB, 55 views)
JR22 is offline   Reply With Quote

Old   April 1, 2013, 09:25
Default Maybe the outlet is not an outlet?
  #3
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 125
Rep Power: 8
JR22 will become famous soon enough
While going back to another thread ( visualise vector field ), I tried to put the money where my mouth was, and found that actually the outlet on the duct is not really an outlet in the strict sense. I mapped the glyphs on the outlet surface of the duct, and there is actually some flow going in. So it is possible that all of the flow that goes through the center doesn't necessarily addup to the outlet of the duct, or should it? Does this make sense? See the image below.

Attached Images
File Type: jpg duct_within_wind_tunnel_3.jpg (87.5 KB, 53 views)
JR22 is offline   Reply With Quote

Old   April 1, 2013, 09:51
Default
  #4
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
hi
an interesting job.how do you resample and integrate on planes as you told?could you explain me a little?
thanks.
immortality is offline   Reply With Quote

Old   April 1, 2013, 10:39
Default Integrate on Mapped Surface using Paraview
  #5
Senior Member
 
JR22's Avatar
 
Jose Rey
Join Date: Oct 2012
Posts: 125
Rep Power: 8
JR22 will become famous soon enough
I posted the steps to get the Glyphs on the mapped surface on the thread in the link (bottom of post). I am replacing the Glyphs with the Integrate here. If you are doing both, they have to both be children of the ResampleWithDataset Object.



1. Create a Plane from the Source menu,
2. Change the X and Y Resolution of the plane to something like 200x100 (x,y). By default, the resolution is 1x1, which is wrong if you are going to do step 3 below.
3. Use the Resample with Dataset filter (the source is the created plane, and the input is your entire dataset).
4. Use Integrate Variables filter on your ResampleWithDataset Object. The integration results will show on a new pane with a spreadsheet kind of look. You can go back to changing the resolution of your source (part of step 2) until your results don't vary much.

Paraphrased from (that one shows how to do the glyphs):
visualise vector field
zaynah04 likes this.

Last edited by JR22; April 1, 2013 at 10:59. Reason: Link to surface mapping for arrow/glyph
JR22 is offline   Reply With Quote

Reply

Tags
continuity, duct, helyx-os, simplefoam, wind tunnel

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
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
Wind Tunnel Blockage question ohma0043 CFX 4 February 13, 2011 10:12
Transient simulation not converging skabilan OpenFOAM Running, Solving & CFD 12 September 17, 2007 17:48
wind engineering problem justin Main CFD Forum 0 February 20, 2006 20:23
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


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