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

Question about the Energy Equation (Turbomachinery)

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

LinkBack Thread Tools Search this Thread Display Modes
Old   May 11, 2024, 05:47
Question Question about the Energy Equation (Turbomachinery)
New Member
Join Date: Jun 2023
Location: Austria
Posts: 9
Rep Power: 2
Lor_enz is on a distinguished road
Hello Foamers,

Maybe someone can clarify or have a suggestion what is wrong with my temperature field. I am solving a compressible flow through a turbomachinery (radial turbine). The model includes a MRF zone, and between the zones a mixing plane treatment is included.
The mesh was generated using Ansys TurboGrid and was transferred to OpenFOAM. The same mesh was used in a CFX simulation to compare results with literature (APU Test Case).

Case Setup:
  • Inlet: total pressure|total temperature with pressureVelocityInletOutlet
  • Outlet: total pressure|inletOutlet T with pressureVelocityInletOutlet
  • Periodics: cyclicAMI
  • Interfaces: stator<>rotor, rotor<>diffuser -> mixing plane
  • Walls: fixedValue for U, else zeroGradient and wall functions
  • komega SST turbulence model, y+ values in a good range
  • BC were ramped up
  • MRF rotational speed 7508.41 rad/s
  • Thermo: psiThermo, sutherland, sensibleEnthalpy, fluid: air

Primarily I am using rhoSimpleFoam to validate the setup and compare the results. Pressure and velocity fields tend towards the same results obtained using CFX. However, the temperature field is fully off. I do not see any cooling inside the turbine, although the pressure and velocity reduce along the flow path.

Energy equation in rhoSimpleFoam
 fvScalarMatrix EEqn
        fvm::div(phi, he)
      + (
   == "e"
          ? fvc::div(phi, volScalarField("Ekp", 0.5*magSqr(U) + p/rho))
          : fvc::div(phi, volScalarField("K", 0.5*magSqr(U)))
      - fvm::laplacian(turbulence->alphaEff(), he)
        fvOptions(rho, he)
For my understanding, the kinetic part of the energy equation should reduce along the blade.
In the UEqn the MRF zone is considered adding the Coriolis and centrifugal force to the system using the term

OpenFOAM solves in absolute velocity formulation. Therefore, quite high velocities are reached inside the rotational domain. Could this affect the kinetic part of the simulation?

Am I missing some term in the energy equation or the UEqn that reduces the energy in the system due to the MRF zone? Pressure gradients at the blade surfaces do look good, and I see a moment acting on those surfaces (post processed using the force library).

Attached the pressure, velocity, tempreaure and total tempreature field along the flow path of the turbine.

looking forward to the discussion,
Attached Images
File Type: png uField.png (126.1 KB, 8 views)
File Type: png pField.png (105.5 KB, 6 views)
File Type: png tField.png (146.6 KB, 5 views)
File Type: png tTotalField.png (118.7 KB, 5 views)
Lorenz H.

PhD Student
Montanuniversity Leoben
Lor_enz is offline   Reply With Quote


energy equation, mrf, rhosimplefoam, turbomachery

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
Energy equation - where is the turbulent kinetic energy? usv001 OpenFOAM Programming & Development 1 January 25, 2022 15:04
Energy equation for compressible flow majkl OpenFOAM Running, Solving & CFD 5 February 10, 2021 15:17
Fluent do not use my velocity field(by UDF) to solve energy equation tangleiplus Fluent UDF and Scheme Programming 6 January 21, 2019 21:28
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Need help:about energy equation in CFX Stein CFX 4 July 2, 2009 22:31

All times are GMT -4. The time now is 16:27.