Skip to content

Add Velocity volume output for compressible flow solvers #2077

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
pcarruscag merged 12 commits into from
Aug 7, 2023
Merged

Add Velocity volume output for compressible flow solvers #2077

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
7dfc9c2
Add Velocity Output To CFlowCompOutput.cpp
MastermindScope Jul 10, 2023
9615189
Add Velocity Volume Output To CFlowCompOutput.cpp
MastermindScope Jul 11, 2023
b5636c4
Added Velocity Output To CNEMOCompOutput.cpp
MastermindScope Jul 11, 2023
cd58b68
Fixed spelling in description
MastermindScope Jul 11, 2023
ced15ae
Added my name to AUTHORS.md
MastermindScope Jul 11, 2023
37d0954
Update SU2_CFD/src/output/CNEMOCompOutput.cpp
pcarruscag Jul 11, 2023
a459e25
Merge branch 'develop' into patch-1
pcarruscag Jul 11, 2023
03623bb
Added SetVolumeOutputValue for Velocity to CFlowCompOutput.cpp
MastermindScope Jul 11, 2023
d0b7922
Add SetVolumeOutputValue for velocity to CNEMOCompOutput.cpp
MastermindScope Jul 11, 2023
a31b800
Update CFlowCompOutput.cpp with suggestions from EvertBunschoten
MastermindScope Jul 17, 2023
6c156e2
Update CNEMOCompOutput.cpp with suggestions from EvertBunschoten and …
MastermindScope Jul 17, 2023
b3e4bdd
Update CFlowCompOutput.cpp to reflect pcarruscag's suggestions
MastermindScope Jul 17, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
1 change: 1 addition & 0 deletions AUTHORS.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -132,6 +132,7 @@ Thomas D. Economon
Tim Albring
TobiKattmann
Trent Lukaczyk
Vinzenz Götz
VivaanKhatri
Wally Maier
Y. Chandukrishna
Expand Down
11 changes: 11 additions & 0 deletions SU2_CFD/src/output/CFlowCompOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -218,6 +218,7 @@ void CFlowCompOutput::SetVolumeOutputFields(CConfig *config){
AddVolumeOutput("DENSITY", "Density", "SOLUTION", "Density");
AddVolumeOutput("MOMENTUM-X", "Momentum_x", "SOLUTION", "x-component of the momentum vector");
AddVolumeOutput("MOMENTUM-Y", "Momentum_y", "SOLUTION", "y-component of the momentum vector");

if (nDim == 3)
AddVolumeOutput("MOMENTUM-Z", "Momentum_z", "SOLUTION", "z-component of the momentum vector");
AddVolumeOutput("ENERGY", "Energy", "SOLUTION", "Energy");
Expand All @@ -237,6 +238,11 @@ void CFlowCompOutput::SetVolumeOutputFields(CConfig *config){
AddVolumeOutput("TEMPERATURE", "Temperature", "PRIMITIVE", "Temperature");
AddVolumeOutput("MACH", "Mach", "PRIMITIVE", "Mach number");
AddVolumeOutput("PRESSURE_COEFF", "Pressure_Coefficient", "PRIMITIVE", "Pressure coefficient");
AddVolumeOutput("VELOCITY-X", "Velocity_x", "PRIMITIVE", "x-component of the velocity vector");
AddVolumeOutput("VELOCITY-Y", "Velocity_y", "PRIMITIVE", "y-component of the velocity vector");

if (nDim == 3)
AddVolumeOutput("VELOCITY-Z", "Velocity_z", "PRIMITIVE", "z-component of the velocity vector");

// Datadriven fluid model
if(config->GetKind_FluidModel() == DATADRIVEN_FLUID){
Expand Down Expand Up @@ -330,6 +336,11 @@ void CFlowCompOutput::LoadVolumeData(CConfig *config, CGeometry *geometry, CSolv

const su2double factor = solver[FLOW_SOL]->GetReferenceDynamicPressure();
SetVolumeOutputValue("PRESSURE_COEFF", iPoint, (Node_Flow->GetPressure(iPoint) - solver[FLOW_SOL]->GetPressure_Inf())/factor);
SetVolumeOutputValue("VELOCITY-X", iPoint, Node_Flow->GetVelocity(iPoint, 0));
SetVolumeOutputValue("VELOCITY-Y", iPoint, Node_Flow->GetVelocity(iPoint, 1));
if (nDim == 3){
SetVolumeOutputValue("VELOCITY-Z", iPoint, Node_Flow->GetVelocity(iPoint, 2));
}

if(config->GetKind_FluidModel() == DATADRIVEN_FLUID){
SetVolumeOutputValue("EXTRAPOLATION", iPoint, Node_Flow->GetDataExtrapolation(iPoint));
Expand Down
7 changes: 7 additions & 0 deletions SU2_CFD/src/output/CNEMOCompOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -242,6 +242,10 @@ void CNEMOCompOutput::SetVolumeOutputFields(CConfig *config){
AddVolumeOutput("PRESSURE", "Pressure", "PRIMITIVE", "Pressure");
AddVolumeOutput("TEMPERATURE_TR", "Temperature_tr", "PRIMITIVE", "Temperature_tr");
AddVolumeOutput("TEMPERATURE_VE", "Temperature_ve", "PRIMITIVE", "Temperature_ve");
AddVolumeOutput("VELOCITY-X", "Velocity_x", "PRIMITIVE", "x-component of the velocity vector");
AddVolumeOutput("VELOCITY-Y", "Velocity_y", "PRIMITIVE", "y-component of the velocity vector");
if (nDim == 3)
AddVolumeOutput("VELOCITY-Z", "Velocity_z", "PRIMITIVE", "z-component of the velocity vector");

AddVolumeOutput("MACH", "Mach", "PRIMITIVE", "Mach number");
AddVolumeOutput("PRESSURE_COEFF", "Pressure_Coefficient", "PRIMITIVE", "Pressure coefficient");
Expand Down Expand Up @@ -313,10 +317,13 @@ void CNEMOCompOutput::LoadVolumeData(CConfig *config, CGeometry *geometry, CSolv

SetVolumeOutputValue("MOMENTUM-X", iPoint, Node_Flow->GetSolution(iPoint, nSpecies));
SetVolumeOutputValue("MOMENTUM-Y", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+1));
SetVolumeOutputValue("VELOCITY-X", iPoint, Node_Flow->GetPrimitive(iPoint, nSpecies));
Comment thread
EvertBunschoten marked this conversation as resolved.
SetVolumeOutputValue("VELOCITY-Y", iPoint, Node_Flow->GetPrimitive(iPoint, nSpecies+1));
if (nDim == 3){
SetVolumeOutputValue("MOMENTUM-Z", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+2));
SetVolumeOutputValue("ENERGY", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+3));
SetVolumeOutputValue("ENERGY_VE", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+4));
SetVolumeOutputValue("VELOCITY-Z", iPoint, Node_Flow->GetPrimitive(iPoint, nSpecies+2));
Comment thread
pcarruscag marked this conversation as resolved.
} else {
SetVolumeOutputValue("ENERGY", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+2));
SetVolumeOutputValue("ENERGY_VE", iPoint, Node_Flow->GetSolution(iPoint, nSpecies+3));
Expand Down