CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

problem with total enthalpy solved by rhoporousmrfsimplefoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 26, 2012, 08:59
Default problem with total enthalpy solved by rhoporousmrfsimplefoam
  #1
New Member
 
Join Date: Oct 2011
Posts: 4
Rep Power: 5
mump0 is on a distinguished road
Hi all!

I'm solving a compressor with rhoporousmrfsimplefoam and my problem is,that my real total enthalpy is below the isentrope one at the outlet. So i have an efficieny factor above 1(what a nice factor ), but it should be near by 0.7(CFX). At the outlet my pressure is near by the value of CFX(a little bit to high), total temperature is to low and my density is to high. To find my mistake much faster i created a case with a simple pipe and the same settings and solver. I changed some BC's and other settings, but there is still the same problem. When i solve it laminar with low velocity the result looks fine. Is there a problem with my komegasst-model? Or maybe my BC's are not correct?
These are the files of the pipe-case:
Quote:
p:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1.0133e5;
boundaryField
{
OUTLET
{
type zeroGradient;
}
WALL
{
type zeroGradient;
}
INLET
{
type totalPressure;
U U;
psi psi;
phi phi;
gamma 1.4;
p0 uniform 1.0133e5;
value uniform 1.0133e5;
}
}
// ************************************************** *********************** //

U:

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
OUTLET
{
type flowRateInletVelocity;
flowRate -0.2;
value uniform (0 0 0);
}
WALL
{
type fixedValue;
value uniform (0 0 0);
}
INLET
{
type zeroGradient;
}
}
// ************************************************** *********************** //

T:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 298.15;
boundaryField
{
OUTLET
{
type zeroGradient;
}
INLET
{
type totalTemperature;
U U;
psi psi;
phi phi;
gamma 1.4;
T0 uniform 298.15;
value uniform 298.15;
}
WALL
{
type zeroGradient;
}
}
// ************************************************** *********************** //

k:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1;
boundaryField
{
OUTLET
{
type zeroGradient;
}
INLET
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.05;
value uniform 1;
}
WALL
{
type compressible::kqRWallFunction;
value uniform 1;
}
}

// ************************************************** *********************** //

omega:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 500;
boundaryField
{
OUTLET
{
type zeroGradient;
}
INLET
{
type fixedValue;
value uniform 500;
}
WALL
{
type compressible:megaWallFunction;
value uniform 500;
}
}

// ************************************************** *********************** //

alphat and mut are default.
Does someone has any idea?

Regards, Eric
mump0 is offline   Reply With Quote

Old   February 14, 2012, 05:21
Default
  #2
New Member
 
Join Date: Oct 2011
Posts: 4
Rep Power: 5
mump0 is on a distinguished road
edit: In the case with the simple pipe the total temperature at the outlet is lower than the total temperature at the inlet, but it should be the same, because its adiabatic. I think my BC's are correct, aren't they? So is there a problem with the h-equation or something else?

Thank you in advance!

Regards, Eric
mump0 is offline   Reply With Quote

Old   February 16, 2012, 06:31
Default
  #3
dre
New Member
 
Martin
Join Date: Feb 2012
Location: Germany
Posts: 13
Rep Power: 0
dre is on a distinguished road
I encountered a similar issue. I was simulating a compressor wheel with rhoPorousMRFSimpleFoam using OpenFOAM 2.1.x. Despite compression, the temperature decreased by about 30K. Previous simulations with openFOAM 2.0.x resulted in an increase of about 30K.

So I took a closer look at the hEqn.H of the solver:

OpenFOAM 2.1.x
Code:
    fvScalarMatrix hEqn
    (
        fvm::div(phi, h)
      - fvm::Sp(fvc::div(phi), h)
      - fvm::laplacian(turbulence->alphaEff(), h)
     ==
      - fvc::div(phi, 0.5*magSqr(U), "div(phi,K)")
    );
OpenFOAM 2.0.x
Code:
    fvScalarMatrix hEqn
    (
        fvm::div(phi, h)
      - fvm::Sp(fvc::div(phi), h)
      - fvm::laplacian(turbulence->alphaEff(), h)
     ==
        fvc::div(phi/fvc::interpolate(rho), rho/psi, "div(U,p)")
      - (rho/psi)*fvc::div(phi/fvc::interpolate(rho))
    );
So there has been a change in the enthalpy equation. However, I don't know the reason for that. Can someone help to clarify this?
dre is offline   Reply With Quote

Old   January 28, 2013, 03:06
Default
  #4
Member
 
prasant
Join Date: Jan 2013
Posts: 29
Rep Power: 4
prasant is on a distinguished road
Hello All,

Did you get any information regarding correction in hEqn If so, please share the information with me. I tried to modified the hEqn as per Eric comments. But still outlet temperature is decreasing.

Regards
Prasant.
prasant is offline   Reply With Quote

Reply

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
Gambit - meshing over airfoil wrapping (?) problem JFDC FLUENT 1 July 11, 2011 05:59
BuoyantFoam and total pressure bc problem is a bug ariorus OpenFOAM Bugs 3 January 28, 2008 05:20
meaning static enthalpy and total enthalpy mspark CD-adapco 1 April 1, 2004 06:20
looking for a problem could be solved by CFX5.5.1 Hason Main CFD Forum 1 October 16, 2002 16:01
total pressure change in a pipe flow Im FLUENT 4 February 6, 2002 20:09


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