|
[Sponsors] |
rhoSimpleFoam MRF cyclicACMI centrifugal compressor |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 16, 2017, 10:56 |
rhoSimpleFoam MRF cyclicACMI centrifugal compressor
|
#1 |
New Member
D. Cook
Join Date: Apr 2013
Posts: 4
Rep Power: 13 |
Hi,
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 Setup 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 Symptoms 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? |
|
Tags |
compressor, cyclicacmi, mrf, rhosimplefoam |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Centrifugal Compressor Efficiency | albertotato | CFX | 3 | May 2, 2020 02:54 |
Temperature for Centrifugal compressor | nnkd | FLUENT | 0 | February 11, 2016 01: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 10:51 |