CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

rhoSimpleFoam MRF cyclicACMI centrifugal compressor

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

LinkBack Thread Tools Display Modes
Old   June 16, 2017, 10:56
Default rhoSimpleFoam MRF cyclicACMI centrifugal compressor
New Member
D. Cook
Join Date: Apr 2013
Posts: 1
Rep Power: 0
dcook is on a distinguished road

This is my first post...I have been experimenting with OpenFOAM for a few years and have been able to find a wealth of information via Google, these forums, research papers, tutorials, etc to solve just about any problem I have run into...until now.

Using OpenFOAM-4.1

Geometry Summary
Centrifugal Fan/Compressor with impeller and spiral regions
Regions are meshed separately and coupled through cyclicACMI

Impeller region is MRF at -309 rad/s, spiral region is stationary
Impeller inlet is totalPressure at 101325
Spiral outlet; have tried several with same results (fixedValue, inletFlowRateVelocity, coded BC to act as restriction, etc).
Turbulence: k-omegaSST
All walls employ wallFunctions
Interface patches are cyclicACMI
Thermophysical Properties hePsiThermo, pureMixture, perfectGas, const transport, hConst, sensibleEnthalpy, standard air (nMoles 1 | molWeight 28.9 | Cp 1007 | Hf 0 | mu 1.8e-05 | Pr 0.7

I have been able to run simpleFoam and rhoSimpleFoam successfully for various boundary conditions but the results are non-physical.
For a typical system (no work), p/rho^gamma is constant for multiple states. I have checked this on the simulation and this seems to be true and this is confirmed by using inlet and outlet patch run-time post processing averages of p,U,T,rho,h and finding the the enthalpy is approximately equal at the inlet and outlet. Problem is that when work is done on the fluid (as would happen in a compressor), the total enthalpy at the outlet should be higher than the total enthalpy at the inlet. In all cases that I have tried, the outlet temperature is slightly lower than the inlet temperature but this is likely due to a slightly higher velocity on the outlet.

I have a higher pressure at the outlet than I do at the inlet (~97 kPa at inlet, ~110 kPa at outlet)...using the relation that T1/T0 = (P1/P0)^((gamma - 1)/gamma), we should expect at least isentropic compression temperature of the gas to about 320K but this is not observed.

I have tested this with OpenFOAM-Extend-4 version using the steadyCompressibleMRFFoam solver and the enthalpy/temperature behave as expected. I have tried porting the rhoCompressibleMRFFoam into an OpenFOAM-4.1 compatible version (changing as little of the extend code as possible) but the issue returns. This makes me think that possibly something in the extend version of MRF or psiThermo is different that is causing this non-physical behavior but I am just not as adept at scouring the code to determine the differences.

Can anyone provide some insight here?
dcook is offline   Reply With Quote


compressor, cyclicacmi, mrf, rhosimplefoam

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
Centrifugal Compressor Efficiency albertotato CFX 0 May 25, 2016 11:44
Temperature for Centrifugal compressor nnkd FLUENT 0 February 11, 2016 02:55
Centrifugal compressor mass flow error Attesz CFX 18 May 27, 2012 10:17
how to solve the diverage of high speed centrifugal compressor, CFD code is STAR CCM layth STAR-CCM+ 3 May 21, 2012 05:48
[ICEM] Meshing a centrifugal compressor impeller Mitpostdoc ANSYS Meshing & Geometry 8 February 25, 2011 11:51

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