Skip to content

Output heatfluxes (total and maximum) per surface #1534

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 5 commits into from
Feb 4, 2022
Merged

Output heatfluxes (total and maximum) per surface #1534

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
cc9e452
output heatfluxes per surface
pcarruscag Feb 3, 2022
1aae062
fix for MZ outputs per surf
pcarruscag Feb 3, 2022
1372590
you always get more than you bargain for...
pcarruscag Feb 4, 2022
2aff612
Merge remote-tracking branch 'upstream/develop' into output_hf_per_surf
pcarruscag Feb 4, 2022
6325e06
Tobi's suggestion
pcarruscag Feb 4, 2022
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
25 changes: 21 additions & 4 deletions SU2_CFD/include/output/CFlowOutput.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -140,21 +140,38 @@ class CFlowOutput : public CFVMOutput{
* \brief Add aerodynamic coefficients as output fields
* \param[in] config - Definition of the particular problem.
*/
void AddAerodynamicCoefficients(CConfig *config);
void AddAerodynamicCoefficients(const CConfig* config);

/*!
* \brief Set the value of the aerodynamic coefficients
* \param[in] config - Definition of the particular problem.
* \param[in] flow_solver - The container holding all solution data.
*/
void SetAerodynamicCoefficients(CConfig *config, CSolver *flow_solver);
void SetAerodynamicCoefficients(const CConfig* config, const CSolver* flow_solver);

/*!
* \brief Set the value of the rotating frame coefficients (CT, CQ and CMerit).
* \brief Add heat flux coefficients as output fields
* \param[in] config - Definition of the particular problem.
*/
void AddHeatCoefficients(const CConfig* config);

/*!
* \brief Set the value of the heat flux coefficients
* \param[in] config - Definition of the particular problem.
* \param[in] flow_solver - The container holding all solution data.
*/
void SetHeatCoefficients(const CConfig* config, const CSolver* flow_solver);

/*!
* \brief Add rotating frame coefficients as output fields.
*/
void AddRotatingFrameCoefficients();

/*!
* \brief Set the value of the rotating frame coefficients (CT, CQ and CMerit).
* \param[in] flow_solver - The container holding all solution data.
*/
void SetRotatingFrameCoefficients(CConfig *config, CSolver *flow_solver);
void SetRotatingFrameCoefficients(const CSolver* flow_solver);

/*!
* \brief Add CP inverse design output as history fields
Expand Down
2 changes: 1 addition & 1 deletion SU2_CFD/include/output/COutput.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -621,7 +621,7 @@ class COutput {
if (it != historyOutputPerSurface_Map.end()) {
it->second[iMarker].value = value;
} else {
SU2_MPI::Error(string("Cannot find output field with name ") + name, CURRENT_FUNCTION);
SU2_MPI::Error("Cannot find output field with name " + name, CURRENT_FUNCTION);
}
}

Expand Down
2 changes: 1 addition & 1 deletion SU2_CFD/include/solvers/CFVMFlowSolverBase.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2014,7 +2014,7 @@ class CFVMFlowSolverBase : public CSolver {
* \brief Provide the total heat load.
* \return Value of the heat load (viscous contribution).
*/
inline su2double GetTotal_HeatFlux(void) const final { return Total_Heat; }
inline su2double GetTotal_HeatFlux() const final { return Total_Heat; }

/*!
* \brief Provide the total heat load.
Expand Down
42 changes: 8 additions & 34 deletions SU2_CFD/src/output/CFlowCompOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -167,36 +167,14 @@ void CFlowCompOutput::SetHistoryOutputFields(CConfig *config){
AddHistoryOutput("LINSOL_RESIDUAL", "LinSolRes", ScreenOutputFormat::FIXED, "LINSOL", "Residual of the linear solver.");
AddHistoryOutputFields_ScalarLinsol(config);

/// BEGIN_GROUP: ROTATING_FRAME, DESCRIPTION: Coefficients related to a rotating frame of reference.
/// DESCRIPTION: Merit
AddHistoryOutput("FIGURE_OF_MERIT", "CMerit", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "Merit", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CT
AddHistoryOutput("THRUST", "CT", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "CT", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CQ
AddHistoryOutput("TORQUE", "CQ", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "CQ", HistoryFieldType::COEFFICIENT);
/// END_GROUP

/// BEGIN_GROUP: EQUIVALENT_AREA, DESCRIPTION: Equivalent area.
/// DESCRIPTION: Nearfield obj. function
AddHistoryOutput("NEARFIELD_OF", "CNearFieldOF", ScreenOutputFormat::SCIENTIFIC, "EQUIVALENT_AREA", "Nearfield obj. function", HistoryFieldType::COEFFICIENT);
/// END_GROUP

/// /// BEGIN_GROUP: HEAT_COEFF, DESCRIPTION: Heat coefficients on all surfaces set with MARKER_MONITORING.
/// DESCRIPTION: Total heatflux
AddHistoryOutput("TOTAL_HEATFLUX", "HF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total heatflux on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: Maximal heatflux
AddHistoryOutput("MAXIMUM_HEATFLUX", "maxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total maximum heatflux on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// END_GROUP

AddHistoryOutput("MIN_DELTA_TIME", "Min DT", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current minimum local time step");
AddHistoryOutput("MAX_DELTA_TIME", "Max DT", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current maximum local time step");

AddHistoryOutput("MIN_CFL", "Min CFL", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current minimum of the local CFL numbers");
AddHistoryOutput("MAX_CFL", "Max CFL", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current maximum of the local CFL numbers");
AddHistoryOutput("AVG_CFL", "Avg CFL", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current average of the local CFL numbers");

/// /// BEGIN_GROUP: FIXED_CL, DESCRIPTION: Relevant outputs for the Fixed CL mode

/// BEGIN_GROUP: FIXED_CL, DESCRIPTION: Relevant outputs for the Fixed CL mode
if (config->GetFixed_CL_Mode()){
/// DESCRIPTION: Difference between current and target CL
AddHistoryOutput("DELTA_CL", "Delta_CL", ScreenOutputFormat::SCIENTIFIC, "FIXED_CL", "Difference between Target CL and current CL", HistoryFieldType::COEFFICIENT);
Expand All @@ -207,6 +185,7 @@ void CFlowCompOutput::SetHistoryOutputFields(CConfig *config){
/// DESCRIPTION: AOA control command by the CL Driver
AddHistoryOutput("CL_DRIVER_COMMAND", "CL_Driver_Command", ScreenOutputFormat::SCIENTIFIC, "FIXED_CL", "CL Driver's control command", HistoryFieldType::RESIDUAL);
}
/// END_GROUP

if (config->GetDeform_Mesh()){
AddHistoryOutput("DEFORM_MIN_VOLUME", "MinVolume", ScreenOutputFormat::SCIENTIFIC, "DEFORM", "Minimum volume in the mesh");
Expand All @@ -215,23 +194,19 @@ void CFlowCompOutput::SetHistoryOutputFields(CConfig *config){
AddHistoryOutput("DEFORM_RESIDUAL", "DeformRes", ScreenOutputFormat::FIXED, "DEFORM", "Residual of the linear solver for the mesh deformation");
}

/*--- Add analyze surface history fields --- */

AddAnalyzeSurfaceOutput(config);

/*--- Add aerodynamic coefficients fields --- */

AddAerodynamicCoefficients(config);

if (config->GetViscous()) {
AddHistoryOutput("BUFFET", "Buffet", ScreenOutputFormat::SCIENTIFIC, "AERO_COEFF", "Buffet sensor", HistoryFieldType::COEFFICIENT);
}

/*--- Add Cp diff fields ---*/
AddHeatCoefficients(config);

Add_CpInverseDesignOutput();
AddRotatingFrameCoefficients();

/*--- Add nearfield diff fields ---*/
Add_CpInverseDesignOutput();

Add_NearfieldInverseDesignOutput();

Expand Down Expand Up @@ -421,9 +396,6 @@ void CFlowCompOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CSol
}
}

SetHistoryOutputValue("TOTAL_HEATFLUX", flow_solver->GetTotal_HeatFlux());
SetHistoryOutputValue("MAXIMUM_HEATFLUX", flow_solver->GetTotal_MaxHeatFlux());

SetHistoryOutputValue("MIN_DELTA_TIME", flow_solver->GetMin_Delta_Time());
SetHistoryOutputValue("MAX_DELTA_TIME", flow_solver->GetMax_Delta_Time());

Expand Down Expand Up @@ -462,9 +434,11 @@ void CFlowCompOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CSol
SetHistoryOutputValue("BUFFET", flow_solver->GetTotal_Buffet_Metric());
}

SetHeatCoefficients(config, flow_solver);

/*--- Set rotating frame coefficients --- */

SetRotatingFrameCoefficients(config, flow_solver);
SetRotatingFrameCoefficients(flow_solver);

/*--- Set Cp diff fields ---*/

Expand Down
60 changes: 19 additions & 41 deletions SU2_CFD/src/output/CFlowIncOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -153,25 +153,6 @@ void CFlowIncOutput::SetHistoryOutputFields(CConfig *config){
if (config->AddRadiation()) AddHistoryOutput("BGS_RAD_ENERGY", "bgs[E_Rad]", ScreenOutputFormat::FIXED, "BGS_RES", "BGS residual of the radiative energy.", HistoryFieldType::RESIDUAL);
/// END_GROUP

/// BEGIN_GROUP: ROTATING_FRAME, DESCRIPTION: Coefficients related to a rotating frame of reference.
/// DESCRIPTION: Merit
AddHistoryOutput("FIGURE_OF_MERIT", "CMerit", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "Merit", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CT
AddHistoryOutput("THRUST", "CT", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "CT", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CQ
AddHistoryOutput("TORQUE", "CQ", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "CQ", HistoryFieldType::COEFFICIENT);
/// END_GROUP

/// BEGIN_GROUP: HEAT_COEFF, DESCRIPTION: Heat coefficients on all surfaces set with MARKER_MONITORING.
/// DESCRIPTION: Total heatflux
AddHistoryOutput("TOTAL_HEATFLUX", "HF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total heatflux on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: Maximal heatflux
AddHistoryOutput("MAXIMUM_HEATFLUX", "maxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total maximum heatflux on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: Temperature
if (heat || weakly_coupled_heat)
AddHistoryOutput("AVG_TEMPERATURE", "Temp", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total avg. temperature on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// END_GROUP

/// DESCRIPTION: Angle of attack
AddHistoryOutput("AOA", "AoA", ScreenOutputFormat::SCIENTIFIC,"AOA", "Angle of attack");
/// DESCRIPTION: Linear solver iterations
Expand Down Expand Up @@ -199,14 +180,17 @@ void CFlowIncOutput::SetHistoryOutputFields(CConfig *config){
AddHistoryOutput("STREAMWISE_HEAT", "SWHeat", ScreenOutputFormat::FIXED, "STREAMWISE_PERIODIC", "Integrated heat for streamwise periodic flow");
}

/*--- Add analyze surface history fields --- */

AddAnalyzeSurfaceOutput(config);

/*--- Add aerodynamic coefficients fields --- */

AddAerodynamicCoefficients(config);

if (heat || weakly_coupled_heat) {
AddHeatCoefficients(config);
AddHistoryOutput("AVG_TEMPERATURE", "Temp", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Average temperature on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
Copy link
Copy Markdown
Contributor

@TobiKattmann TobiKattmann Feb 4, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just as a side note: For an AvgTemp OF for the flow solver one needs to use SURFACE_STATIC_TEMPERATURE (which works with MARKER_ANALYZE) and not AVG_TEMPERATURE. For solid domains it is AVG_TEMPERATURE though.

Copy link
Copy Markdown
Member Author

@pcarruscag pcarruscag Feb 4, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

What about weakly_coupled_heat?

Copy link
Copy Markdown
Contributor

@TobiKattmann TobiKattmann Feb 4, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yep, you're right. There one would use AVG_TEMPERATURE

Copy link
Copy Markdown
Member Author

@pcarruscag pcarruscag Feb 4, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

🤦 I won't stick my hand in that xD

}

AddRotatingFrameCoefficients();

}

void CFlowIncOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CSolver **solver) {
Expand Down Expand Up @@ -239,25 +223,19 @@ void CFlowIncOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CSolv
}

if (weakly_coupled_heat){
SetHistoryOutputValue("TOTAL_HEATFLUX", heat_solver->GetTotal_HeatFlux());
SetHistoryOutputValue("MAXIMUM_HEATFLUX", heat_solver->GetTotal_MaxHeatFlux());
SetHistoryOutputValue("AVG_TEMPERATURE", heat_solver->GetTotal_AvgTemperature());
SetHistoryOutputValue("RMS_TEMPERATURE", log10(heat_solver->GetRes_RMS(0)));
SetHistoryOutputValue("MAX_TEMPERATURE", log10(heat_solver->GetRes_Max(0)));
SetHeatCoefficients(config, heat_solver);
SetHistoryOutputValue("AVG_TEMPERATURE", heat_solver->GetTotal_AvgTemperature());
SetHistoryOutputValue("RMS_TEMPERATURE", log10(heat_solver->GetRes_RMS(0)));
SetHistoryOutputValue("MAX_TEMPERATURE", log10(heat_solver->GetRes_Max(0)));
if (multiZone) SetHistoryOutputValue("BGS_TEMPERATURE", log10(heat_solver->GetRes_BGS(0)));
}
if (heat){
SetHistoryOutputValue("TOTAL_HEATFLUX", flow_solver->GetTotal_HeatFlux());
SetHistoryOutputValue("MAXIMUM_HEATFLUX", flow_solver->GetTotal_MaxHeatFlux());
SetHistoryOutputValue("AVG_TEMPERATURE", flow_solver->GetTotal_AvgTemperature());
if (nDim == 3) SetHistoryOutputValue("RMS_TEMPERATURE", log10(flow_solver->GetRes_RMS(4)));
else SetHistoryOutputValue("RMS_TEMPERATURE", log10(flow_solver->GetRes_RMS(3)));

if (nDim == 3) SetHistoryOutputValue("MAX_TEMPERATURE", log10(flow_solver->GetRes_Max(4)));
else SetHistoryOutputValue("MAX_TEMPERATURE", log10(flow_solver->GetRes_Max(3)));
if (multiZone){
if (nDim == 3) SetHistoryOutputValue("BGS_TEMPERATURE", log10(flow_solver->GetRes_BGS(4)));
else SetHistoryOutputValue("BGS_TEMPERATURE", log10(flow_solver->GetRes_BGS(3)));
if (heat) {
SetHeatCoefficients(config, flow_solver);
SetHistoryOutputValue("AVG_TEMPERATURE", flow_solver->GetTotal_AvgTemperature());
SetHistoryOutputValue("RMS_TEMPERATURE", log10(flow_solver->GetRes_RMS(nDim + 1)));
SetHistoryOutputValue("MAX_TEMPERATURE", log10(flow_solver->GetRes_Max(nDim + 1)));
if (multiZone) {
SetHistoryOutputValue("BGS_TEMPERATURE", log10(flow_solver->GetRes_BGS(nDim + 1)));
}
}

Expand Down Expand Up @@ -296,7 +274,7 @@ void CFlowIncOutput::LoadHistoryData(CConfig *config, CGeometry *geometry, CSolv

/*--- Set rotating frame coefficients --- */

SetRotatingFrameCoefficients(config, flow_solver);
SetRotatingFrameCoefficients(flow_solver);

/*--- Keep this as last, since it uses the history values that were set. ---*/

Expand Down
54 changes: 49 additions & 5 deletions SU2_CFD/src/output/CFlowOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
/*!
* \file CFlowOutput.cpp
* \brief Main subroutines for compressible flow output
* \brief Common functions for flow output.
* \author R. Sanchez
* \version 7.3.0 "Blackbird"
*
Expand Down Expand Up @@ -1083,7 +1083,7 @@ void CFlowOutput::LoadSurfaceData(CConfig *config, CGeometry *geometry, CSolver
SetVolumeOutputValue("Y_PLUS", iPoint, solver[FLOW_SOL]->GetYPlus(iMarker, iVertex));
}

void CFlowOutput::AddAerodynamicCoefficients(CConfig *config){
void CFlowOutput::AddAerodynamicCoefficients(const CConfig* config) {

/// BEGIN_GROUP: AERO_COEFF, DESCRIPTION: Sum of the aerodynamic coefficients and forces on all surfaces (markers) set with MARKER_MONITORING.
/// DESCRIPTION: Drag coefficient
Expand Down Expand Up @@ -1141,7 +1141,7 @@ void CFlowOutput::AddAerodynamicCoefficients(CConfig *config){
AddHistoryOutput("COMBO", "ComboObj", ScreenOutputFormat::SCIENTIFIC, "COMBO", "Combined obj. function value.", HistoryFieldType::COEFFICIENT);
}

void CFlowOutput::SetAerodynamicCoefficients(CConfig *config, CSolver *flow_solver){
void CFlowOutput::SetAerodynamicCoefficients(const CConfig* config, const CSolver* flow_solver){

SetHistoryOutputValue("DRAG", flow_solver->GetTotal_CD());
SetHistoryOutputValue("LIFT", flow_solver->GetTotal_CL());
Expand Down Expand Up @@ -1183,14 +1183,58 @@ void CFlowOutput::SetAerodynamicCoefficients(CConfig *config, CSolver *flow_solv
SetHistoryOutputValue("AOA", config->GetAoA());
}

void CFlowOutput::SetRotatingFrameCoefficients(CConfig *config, CSolver *flow_solver) {
void CFlowOutput::AddHeatCoefficients(const CConfig* config) {

if (!config->GetViscous()) return;

/// BEGIN_GROUP: HEAT, DESCRIPTION: Heat coefficients on all surfaces set with MARKER_MONITORING.
/// DESCRIPTION: Total heatflux
AddHistoryOutput("TOTAL_HEATFLUX", "HF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total heatflux on all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: Maximal heatflux
AddHistoryOutput("MAXIMUM_HEATFLUX", "maxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Maximum heatflux across all surfaces set with MARKER_MONITORING.", HistoryFieldType::COEFFICIENT);

vector<string> Marker_Monitoring;
for (auto iMarker = 0u; iMarker < config->GetnMarker_Monitoring(); iMarker++) {
Marker_Monitoring.push_back(config->GetMarker_Monitoring_TagBound(iMarker));
}
/// DESCRIPTION: Total heatflux
AddHistoryOutputPerSurface("TOTAL_HEATFLUX_ON_SURFACE", "HF", ScreenOutputFormat::SCIENTIFIC, "HEAT_SURF", Marker_Monitoring, HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: Total heatflux
AddHistoryOutputPerSurface("MAXIMUM_HEATFLUX_ON_SURFACE", "maxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT_SURF", Marker_Monitoring, HistoryFieldType::COEFFICIENT);
/// END_GROUP
}

void CFlowOutput::SetHeatCoefficients(const CConfig* config, const CSolver* flow_solver) {

if (!config->GetViscous()) return;

SetHistoryOutputValue("TOTAL_HEATFLUX", flow_solver->GetTotal_HeatFlux());
SetHistoryOutputValue("MAXIMUM_HEATFLUX", flow_solver->GetTotal_MaxHeatFlux());

for (auto iMarker = 0u; iMarker < config->GetnMarker_Monitoring(); iMarker++) {
SetHistoryOutputPerSurfaceValue("TOTAL_HEATFLUX_ON_SURFACE", flow_solver->GetSurface_HF_Visc(iMarker), iMarker);
SetHistoryOutputPerSurfaceValue("MAXIMUM_HEATFLUX_ON_SURFACE", flow_solver->GetSurface_MaxHF_Visc(iMarker), iMarker);
}
}

void CFlowOutput::AddRotatingFrameCoefficients() {
/// BEGIN_GROUP: ROTATING_FRAME, DESCRIPTION: Coefficients related to a rotating frame of reference.
/// DESCRIPTION: Merit
AddHistoryOutput("FIGURE_OF_MERIT", "CMerit", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "Thrust over torque", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CT
AddHistoryOutput("THRUST", "CT", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "Thrust coefficient", HistoryFieldType::COEFFICIENT);
/// DESCRIPTION: CQ
AddHistoryOutput("TORQUE", "CQ", ScreenOutputFormat::SCIENTIFIC, "ROTATING_FRAME", "Torque coefficient", HistoryFieldType::COEFFICIENT);
/// END_GROUP
}

void CFlowOutput::SetRotatingFrameCoefficients(const CSolver* flow_solver) {

SetHistoryOutputValue("THRUST", flow_solver->GetTotal_CT());
SetHistoryOutputValue("TORQUE", flow_solver->GetTotal_CQ());
SetHistoryOutputValue("FIGURE_OF_MERIT", flow_solver->GetTotal_CMerit());
}


void CFlowOutput::Add_CpInverseDesignOutput(){

AddHistoryOutput("INVERSE_DESIGN_PRESSURE", "Cp_Diff", ScreenOutputFormat::FIXED, "CP_DIFF", "Cp difference for inverse design");
Expand Down
4 changes: 2 additions & 2 deletions SU2_CFD/src/output/CHeatOutput.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -104,8 +104,8 @@ void CHeatOutput::SetHistoryOutputFields(CConfig *config){
AddHistoryOutput("BGS_TEMPERATURE", "bgs[T]", ScreenOutputFormat::FIXED, "BGS_RES", "Block-Gauss-Seidel residual of the temperature", HistoryFieldType::RESIDUAL);

AddHistoryOutput("TOTAL_HEATFLUX", "HF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total heatflux on all surfaces defined in MARKER_MONITORING", HistoryFieldType::COEFFICIENT);
AddHistoryOutput("MAXIMUM_HEATFLUX", "MaxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total maximal heatflux on all surfaces defined in MARKER_MONITORING", HistoryFieldType::COEFFICIENT);
AddHistoryOutput("AVG_TEMPERATURE", "AvgTemp", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Total average temperature on all surfaces defined in MARKER_MONITORING", HistoryFieldType::COEFFICIENT);
AddHistoryOutput("MAXIMUM_HEATFLUX", "MaxHF", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Maximum heatflux on all surfaces defined in MARKER_MONITORING", HistoryFieldType::COEFFICIENT);
AddHistoryOutput("AVG_TEMPERATURE", "AvgTemp", ScreenOutputFormat::SCIENTIFIC, "HEAT", "Average temperature on all surfaces defined in MARKER_MONITORING", HistoryFieldType::COEFFICIENT);
AddHistoryOutput("CFL_NUMBER", "CFL number", ScreenOutputFormat::SCIENTIFIC, "CFL_NUMBER", "Current value of the CFL number");

AddHistoryOutput("COMBO", "ComboObj", ScreenOutputFormat::SCIENTIFIC, "COMBO", "Combined obj. function value.", HistoryFieldType::COEFFICIENT);
Expand Down
Loading