CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Angular momentum flux through faceZone

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

LinkBack Thread Tools Display Modes
Old   January 4, 2010, 08:23
Default Angular momentum flux through faceZone
New Member
Join Date: Apr 2009
Posts: 26
Rep Power: 8
askjak is on a distinguished road
I have a cylindrical geometry (see graphics) and I want to evaluate the angular momentum flux into the cylinder through the 30 ports (among other quantities) using a cylindrical coordinate system. The mesh is conforming between the two zones (red and blue).

My first attempt was to create a cyclic patch in each of the 30 ports and write a postprocessing tool based on patchIntegrate to find the angular momentum flux. This worked well formally, but it is very tedious to create the meshes in a commercial mesher in such a way that I can create the cyclic patches with createPatch (and we need a LOT of meshes). Other issues with this method is handling parallel postprocessing (see Cyclic patches and parallel postprocessing problems), and fundamentally I would prefer to get the angular momentum flux out during the simulation and not “postmortem”.

I have now created a faceZone for each of the 30 ports (see pink on graphics) according to the advice given in Cyclic patches and parallel postprocessing problems. I can now use OpenFOAM-1.6.x/src/postProcessing/functionObjects/field to get information out on a faceZone. But as I understand my problem is slightly more complex than the functionality in fieldValues since I need to know the position of the individual face and multiply phi with U before integration etc.

My question now is the following: What would you recommend? In some way extending the functionObjects?

Attached Images
File Type: png zones4.PNG (36.1 KB, 51 views)
askjak is offline   Reply With Quote

Old   March 5, 2010, 15:42
Default FaceZones
Join Date: Sep 2009
Location: Buenos Aires
Posts: 35
Rep Power: 8
lfbarcelo is on a distinguished road
How did you create the faceZones?? I have made a mesh using snappyHexMesh. I meshed an stl corresponding to a surface with no width, inside a wind Tunnel. The problem is that paraView, in an attempt to show surface pressure, shows this surface as if it was "stained" because the pressure is different on each side of this surface. In order to "create" a certain width, ando solve this problem, I am trying to use "create Baffles", but I can't create the faceZone out of the patch, a face zone is necessary as an input agument for "createBaffles".
lfbarcelo is offline   Reply With Quote


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
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Derivation of Momentum Equation in Integral Form Demonwolf Main CFD Forum 2 October 29, 2009 20:53
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 21:09
Angular momentum conservation? Barry Main CFD Forum 19 April 11, 2002 10:37

All times are GMT -4. The time now is 18:42.