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

Discrepancy with phi(mass flowrate) calculation at the patch with buoyantSimpleFoam s

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 16, 2021, 03:51
Default Discrepancy with phi(mass flowrate) calculation at the patch with buoyantSimpleFoam s
  #1
New Member
 
venu
Join Date: May 2019
Posts: 8
Rep Power: 7
Venu Angirekula is on a distinguished road
Description:
The problem is related to the way OpenFOAM evaluates the phi(mass flowrate) at
the patch with buoyantSimpleFoam solver. In my system, I implemented a
codedFixedValue BC for inlet velocity, which enters into the space with some
angle from the patch normal. However, the value of phi at the patch produced by
the OF is less than the Original value. Upon observation, it seems that only the
velocity component perpendicular to the patch was considered to generate the
phi. Which disturbs the conservation of heat in the syste

Steps to Reproduce:
I have implemented the velocity profile at the inlet patch in such way that all
velocity vectors are entering into the space radially with 70deg from the patch
normal, as shown in the attached image.

Theoretically mass flowrate calculations:

Volumetric flowrate of patch = 0.02831 m3/sec
Inlet Temperature = 13 deg C
Corresponding Density,rho = 1.238 kg/m3
Supply mass flowrate = 0.03485 kg/sec

But from the OF, with 'flowRatePatch' utility, Mass flowrate = 0.01194 kg/sec.
Which is in accord with direct summation of 'phi' values on the patch.
Nevertheless, there is large difference between the actual supply mass flowrate
and OF generated flowrate.
Upon further investigation, based on "compressibleCreatePhi.H" we found the
calculation of phi
surfaceScalarField phi ( IOobject ( "phi", runTime.timeName(), mesh,
IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), linearInterpolate(rho*U) &
mesh.Sf() );
The calculation of phi on the face :

The surface area vector, Sf ( 0,0,4.032e-05)
The translated velocity vector,U ( 0.93872955,-0.042532896,-0.34202014)
Mass flowrate = (Sf.U)*rho = (0 + 0 + 1.379e-05)*1.238

Please find the case directory here:
https://drive.google.com/file/d/11Fj...ew?usp=sharing.

Observations & Queries:

Why the phi calculations does not taking the contribution of components other
than patch normal velocity.
Based in the above calculations, there is large difference between the actual
supply mass flowrate and OF generated flowrate.
The generated low phi values have also impacted the heat balance of the system.
Venu Angirekula is offline   Reply With Quote

Reply

Tags
massflowrate, openfoam, phi


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
[CGNS] CGNS converters available mbeaudoin OpenFOAM Meshing & Mesh Conversion 137 December 14, 2018 05:20
Near wall treatment in k-omega SST Arnoldinho OpenFOAM Running, Solving & CFD 38 March 8, 2017 14:48
[Other] How to create an MRF zone ? aminem OpenFOAM Meshing & Mesh Conversion 2 December 8, 2014 11:45
Possible Bug in pimpleFoam (or createPatch) (or fluent3DMeshToFoam) cfdonline2mohsen OpenFOAM 3 October 21, 2013 10:28
Warning 097- AB Siemens 6 November 15, 2004 05:41


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