CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

Velocity output at face centers (instead of cell centers)

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

Like Tree3Likes
  • 3 Post By LuckyTran

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 2, 2018, 02:50
Default Velocity output at face centers (instead of cell centers)
  #1
New Member
 
Join Date: Apr 2018
Posts: 2
Rep Power: 0
myuse2 is on a distinguished road
If I understand correctly with FVM CFD the velocities (fluxes) are computed at face centers while pressures are defined in the cell centers (so as to avoid checkerboard effects).

The output from OpenFOAM however is only cell centered, I wonder if it somehow possible to extract the un-interpolated (raw) face centered velocity vector instead of an interpolated cell centered one?
myuse2 is offline   Reply With Quote

Old   April 4, 2018, 17:54
Default
  #2
Senior Member
 
Lucky Tran
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 4,191
Rep Power: 50
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
The face fluxes are interpolated.

The fluxes are the phi field. In many solvers, if you include a phi file at the start it will write it thereafter (just like rho). If you have an incompressible solver then the phi is the velocity. If compressible you have to calculate the velocity from phi/rho but the rho must be interpolated.

Here is an function object that writes the rho and phi fields. Alternatively you can also go to createFields.H in the solver, change the autoWrite and recompile.

Do you want the fluxes or the face velocities? Because it differs by the face area.
LeeRuns, Marpole and skyoon like this.

Last edited by LuckyTran; April 6, 2018 at 13:46.
LuckyTran is offline   Reply With Quote

Old   February 17, 2020, 09:00
Default
  #3
New Member
 
SOURAV GANGULI
Join Date: Dec 2018
Posts: 9
Rep Power: 4
sourav0921 is on a distinguished road
Hello guys,


I able to print the phi value for eaxh writing time step in pisoFoam. But can you please tell me whether is there any direct way to calculate phi per unit cell area ?
sourav0921 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] Error in mesh writing helios ANSYS Meshing & Geometry 19 January 8, 2019 10:56
[simpleFoam] Floating point exception seuchsy OpenFOAM Running, Solving & CFD 3 November 3, 2016 09:29
How to use "translation" in solidBodyMotionFunction in OpenFOAM rupesh_w OpenFOAM Running, Solving & CFD 5 August 16, 2016 04:27
Output cell face area from all patches niko0807 OpenFOAM Programming & Development 1 October 23, 2015 09:03
[blockMesh] BlockMeshmergePatchPairs hjasak OpenFOAM Meshing & Mesh Conversion 11 August 15, 2008 07:36


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