CFD Online Logo CFD Online URL
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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 2, 2018, 03:50
Default Velocity output at face centers (instead of cell centers)
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, 18:54
Senior Member
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,112
Rep Power: 60
LuckyTran has a spectacular aura aboutLuckyTran 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 14:46.
LuckyTran is offline   Reply With Quote

Old   February 17, 2020, 10:00
New Member
Join Date: Dec 2018
Posts: 9
Rep Power: 6
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


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 21 August 19, 2021 15:18
[simpleFoam] Floating point exception seuchsy OpenFOAM Running, Solving & CFD 3 November 3, 2016 10:29
How to use "translation" in solidBodyMotionFunction in OpenFOAM rupesh_w OpenFOAM Running, Solving & CFD 5 August 16, 2016 05:27
Output cell face area from all patches niko0807 OpenFOAM Programming & Development 1 October 23, 2015 10:03
[blockMesh] BlockMeshmergePatchPairs hjasak OpenFOAM Meshing & Mesh Conversion 11 August 15, 2008 08:36

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