View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0004242 | OpenFOAM | Bug | public | 2025-06-10 10:10 | 2025-06-10 15:20 |
| Reporter | Giorgio Ceresa | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | closed | Resolution | suspended | ||
| Platform | linux64 | OS | RockyLinux | OS Version | 8.10 |
| Product Version | 12 | ||||
| Summary | 0004242: Unphysical temperature in the outlet using MRF and compressible fluid | ||||
| Description | The evaluated test case, NASA Rotor 37, is an axial flow compressor rotor that has been extensively studied for aerodynamic research and validation by computational fluid dynamics (CFD). The analysis of turbulent and subsonic compressible conditions with air as an ideal gas shows incorrect results for the temperature field: - Lower static temperature at the outlet - Total temperature remains unchanged between inlet and outlet. The compressor blade provides work to the fluid; therefore, the expected result for the energy equation is an increase in static temperature and total temperature at the outlet. It seems that the solver does not take the work done by the blades into account in the energy equation using steady-state MRF modeling. | ||||
| Steps To Reproduce | Run Allrun in case directory and visualize temperature and totalTemperature results in Paraview. | ||||
| Additional Information | Refer to attached presentation and case. #Function for creation of totalTemperaure field dynamicCode1 { libs ("libutilityFunctionObjects.so, libthermophysicalProperties.so"); type coded; enabled yes; writeControl writeTime; writeInterval 1; codeWrite #{ // Lookup const volScalarField& T = mesh().lookupObject<volScalarField>("T"); const volVectorField& U = mesh().lookupObject<volVectorField>("U"); const volScalarField& Cp = mesh().lookupObject<volScalarField>("Cp"); volScalarField totalTemperature ( IOobject ( "totalTemperature", mesh().time().name(), mesh(), IOobject::READ_IF_PRESENT, IOobject::AUTO_WRITE ), T + pow(mag(U), 2) / (2 * Cp) ); // Write totalTemperature.write(); #}; } | ||||
| Tags | No tags attached. | ||||
|
|
|
|
|
Will you be contributing to OpenFOAM maintenance funding to cover the work you would like us to do on this for you? https://openfoam.org/maintenance https://openfoam.org/news/funding-2025/ |
|
|
Request for free development, no funding is offered. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2025-06-10 10:10 | Giorgio Ceresa | New Issue | |
| 2025-06-10 10:10 | Giorgio Ceresa | File Added: Rotor37_Compressor_clean.tar.xz | |
| 2025-06-10 10:10 | Giorgio Ceresa | File Added: Bug report.pdf | |
| 2025-06-10 10:28 | henry | Note Added: 0013574 | |
| 2025-06-10 10:29 | henry | Additional Information Updated | |
| 2025-06-10 14:37 | henry | Assigned To | => henry |
| 2025-06-10 14:37 | henry | Status | new => closed |
| 2025-06-10 14:37 | henry | Resolution | open => suspended |
| 2025-06-10 14:37 | henry | Note Added: 0013577 |
