Skip to content

Commit fa88da8

Browse files
thomasdickthomasdick
authored
Merge branch 'develop' into feature_Sobolev_smoothing_solver
2 parents 5598ad2 + 648cf4c commit fa88da8

10 files changed

Lines changed: 56 additions & 89 deletions

File tree

‎Common/src/CConfig.cpp‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -5856,14 +5856,14 @@ void CConfig::SetOutput(SU2_COMPONENT val_software, unsigned short val_izone) {
58565856
break;
58575857
case MAIN_SOLVER::NEMO_EULER:
58585858
if (Kind_Regime == ENUM_REGIME::COMPRESSIBLE) cout << "Compressible two-temperature thermochemical non-equilibrium Euler equations." << endl;
5859-
if(Kind_FluidModel == SU2_NONEQ){
5859+
if (Kind_FluidModel == SU2_NONEQ){
58605860
if ((GasModel != "N2") && (GasModel != "AIR-5") && (GasModel != "ARGON"))
58615861
SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5, ARGON.", CURRENT_FUNCTION);
58625862
}
58635863
break;
58645864
case MAIN_SOLVER::NEMO_NAVIER_STOKES:
58655865
if (Kind_Regime == ENUM_REGIME::COMPRESSIBLE) cout << "Compressible two-temperature thermochemical non-equilibrium Navier-Stokes equations." << endl;
5866-
if(Kind_FluidModel == SU2_NONEQ){
5866+
if (Kind_FluidModel == SU2_NONEQ){
58675867
if ((GasModel != "N2") && (GasModel != "AIR-5") && (GasModel != "ARGON"))
58685868
SU2_MPI::Error("The GAS_MODEL given as input is not valid. Choose one of the options: N2, AIR-5, ARGON.", CURRENT_FUNCTION);
58695869
}

‎SU2_CFD/include/output/CNEMOCompOutput.hpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -45,7 +45,7 @@ class CNEMOCompOutput final: public CFlowOutput {
4545
* \brief Constructor of the class
4646
* \param[in] config - Definition of the particular problem.
4747
*/
48-
CNEMOCompOutput(CConfig *config, unsigned short nDim);
48+
CNEMOCompOutput(const CConfig *config, unsigned short nDim);
4949

5050
/*!
5151
* \brief Load the history output field values

‎SU2_CFD/include/solvers/CSolver.hpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1253,7 +1253,7 @@ class CSolver {
12531253
* \param[in] val_marker - Surface marker where the boundary condition is applied.
12541254
*/
12551255
inline virtual void BC_Smoluchowski_Maxwell(CGeometry *geometry,
1256-
CSolver **solution_container,
1256+
CSolver **solver_container,
12571257
CNumerics *conv_numerics,
12581258
CNumerics *visc_numerics,
12591259
CConfig *config,

‎SU2_CFD/src/drivers/CDriver.cpp‎

Lines changed: 0 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -3123,7 +3123,6 @@ bool CFluidDriver::Monitor(unsigned long ExtIter) {
31233123

31243124
switch (config_container[ZONE_0]->GetKind_Solver()) {
31253125
case MAIN_SOLVER::EULER: case MAIN_SOLVER::NAVIER_STOKES: case MAIN_SOLVER::RANS:
3126-
StopCalc = integration_container[ZONE_0][INST_0][FLOW_SOL]->GetConvergence(); break;
31273126
case MAIN_SOLVER::NEMO_EULER: case MAIN_SOLVER::NEMO_NAVIER_STOKES:
31283127
StopCalc = integration_container[ZONE_0][INST_0][FLOW_SOL]->GetConvergence(); break;
31293128
case MAIN_SOLVER::HEAT_EQUATION:

‎SU2_CFD/src/fluid/CSU2TCLib.cpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -781,7 +781,7 @@ vector<su2double>& CSU2TCLib::ComputeSpeciesEve(su2double val_T, bool vibe_only)
781781
}
782782
Eel = Ru/MolarMass[iSpecies] * (num/denom);
783783
}
784-
if(vibe_only == true) {eves[iSpecies] = Ev;}
784+
if (vibe_only == true) {eves[iSpecies] = Ev;}
785785
else {eves[iSpecies] = Ev + Eel;}
786786
}
787787

‎SU2_CFD/src/output/CNEMOCompOutput.cpp‎

Lines changed: 9 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -4,22 +4,12 @@
44
* \author W. Maier, R. Sanchez
55
* \version 7.2.1 "Blackbird"
66
*
7-
* The current SU2 release has been coordinated by the
8-
* SU2 International Developers Society <www.su2devsociety.org>
9-
* with selected contributions from the open-source community.
7+
* SU2 Project Website: https://su2code.github.io
108
*
11-
* The main research teams contributing to the current release are:
12-
* - Prof. Juan J. Alonso's group at Stanford University.
13-
* - Prof. Piero Colonna's group at Delft University of Technology.
14-
* - Prof. Nicolas R. Gauger's group at Kaiserslautern University of Technology.
15-
* - Prof. Alberto Guardone's group at Polytechnic University of Milan.
16-
* - Prof. Rafael Palacios' group at Imperial College London.
17-
* - Prof. Vincent Terrapon's group at the University of Liege.
18-
* - Prof. Edwin van der Weide's group at the University of Twente.
19-
* - Lab. of New Concepts in Aeronautics at Tech. Institute of Aeronautics.
9+
* The SU2 Project is maintained by the SU2 Foundation
10+
* (http://su2foundation.org)
2011
*
21-
* Copyright 2012-2018, Francisco D. Palacios, Thomas D. Economon,
22-
* Tim Albring, and the SU2 contributors.
12+
* Copyright 2012-2022, SU2 Contributors (cf. AUTHORS.md)
2313
*
2414
* SU2 is free software; you can redistribute it and/or
2515
* modify it under the terms of the GNU Lesser General Public
@@ -40,9 +30,9 @@
4030
#include "../../../Common/include/geometry/CGeometry.hpp"
4131
#include "../../include/solvers/CSolver.hpp"
4232

43-
CNEMOCompOutput::CNEMOCompOutput(CConfig *config, unsigned short nDim) : CFlowOutput(config, nDim, false) {
33+
CNEMOCompOutput::CNEMOCompOutput(const CConfig *config, unsigned short nDim) : CFlowOutput(config, nDim, false) {
4434

45-
turb_model = config->GetKind_Turb_Model();
35+
turb_model = config->GetKind_Turb_Model();
4636
nSpecies = config->GetnSpecies();
4737

4838
/*--- Set the default history fields if nothing is set in the config file ---*/
@@ -118,6 +108,7 @@ CNEMOCompOutput::CNEMOCompOutput(CConfig *config, unsigned short nDim) : CFlowOu
118108
void CNEMOCompOutput::SetHistoryOutputFields(CConfig *config){
119109

120110
/// BEGIN_GROUP: RMS_RES, DESCRIPTION: The root-mean-square residuals of the SOLUTION variables.
111+
/// DESCRIPTION: Root-mean square residual of the species densities.
121112
for(iSpecies = 0; iSpecies < nSpecies; iSpecies++)
122113
AddHistoryOutput("RMS_DENSITY_" + std::to_string(iSpecies), "rms[Rho_" + std::to_string(iSpecies) + "]", ScreenOutputFormat::FIXED, "RMS_RES", "Root-mean square residual of the species density " + std::to_string(iSpecies) + ".", HistoryFieldType::RESIDUAL);
123114
/// DESCRIPTION: Root-mean square residual of the momentum x-component.
@@ -279,7 +270,7 @@ void CNEMOCompOutput::SetVolumeOutputFields(CConfig *config){
279270
AddVolumeOutput("MACH", "Mach", "PRIMITIVE", "Mach number");
280271
AddVolumeOutput("PRESSURE_COEFF", "Pressure_Coefficient", "PRIMITIVE", "Pressure coefficient");
281272

282-
if (config->GetKind_Solver() == MAIN_SOLVER::NEMO_NAVIER_STOKES){
273+
if (config->GetViscous()) {
283274
AddVolumeOutput("LAMINAR_VISCOSITY", "Laminar_Viscosity", "PRIMITIVE", "Laminar viscosity");
284275

285276
AddVolumeOutput("SKIN_FRICTION-X", "Skin_Friction_Coefficient_x", "PRIMITIVE", "x-component of the skin friction vector");
@@ -367,7 +358,7 @@ void CNEMOCompOutput::LoadVolumeData(CConfig *config, CGeometry *geometry, CSolv
367358
su2double factor = 1.0/(0.5*solver[FLOW_SOL]->GetDensity_Inf()*VelMag);
368359
SetVolumeOutputValue("PRESSURE_COEFF", iPoint, (Node_Flow->GetPressure(iPoint) - solver[FLOW_SOL]->GetPressure_Inf())*factor);
369360

370-
if (config->GetKind_Solver() == MAIN_SOLVER::NEMO_NAVIER_STOKES){
361+
if (config->GetViscous()){
371362
SetVolumeOutputValue("LAMINAR_VISCOSITY", iPoint, Node_Flow->GetLaminarViscosity(iPoint));
372363
}
373364

‎SU2_CFD/src/solvers/CEulerSolver.cpp‎

Lines changed: 0 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -8434,8 +8434,6 @@ void CEulerSolver::SetFreeStream_TurboSolution(CConfig *config) {
84348434

84358435
}
84368436

8437-
8438-
84398437
void CEulerSolver::PreprocessAverage(CSolver **solver, CGeometry *geometry, CConfig *config, unsigned short marker_flag) {
84408438

84418439
unsigned long iVertex, iPoint;
@@ -8507,7 +8505,6 @@ void CEulerSolver::PreprocessAverage(CSolver **solver, CGeometry *geometry, CCon
85078505
}
85088506
}
85098507

8510-
85118508
#ifdef HAVE_MPI
85128509

85138510
/*--- Add information using all the nodes ---*/

‎SU2_CFD/src/solvers/CNEMOEulerSolver.cpp‎

Lines changed: 31 additions & 51 deletions
Original file line numberDiff line numberDiff line change
@@ -1026,7 +1026,6 @@ void CNEMOEulerSolver::SetNondimensionalization(CConfig *config, unsigned short
10261026

10271027
bool unsteady = (config->GetTime_Marching() != TIME_MARCHING::STEADY);
10281028
bool viscous = config->GetViscous();
1029-
bool dynamic_grid = config->GetGrid_Movement();
10301029
bool gravity = config->GetGravityForce();
10311030
bool turbulent = false;
10321031
bool tkeNeeded = ((turbulent) && (config->GetKind_Turb_Model() == TURB_MODEL::SST));
@@ -1428,27 +1427,11 @@ void CNEMOEulerSolver::BC_Sym_Plane(CGeometry *geometry, CSolver **solver_contai
14281427
unsigned short iDim, jDim, iSpecies, iVar, jVar;
14291428
unsigned long iPoint, iVertex;
14301429

1431-
su2double *Normal = nullptr, Area, UnitNormal[3], *NormalArea,
1432-
**Jacobian_b, **DubDu,
1433-
rho, cs, P, rhoE, rhoEve, conc, *u, *dPdU;
1430+
/*--- Allocate the necessary vector structures ---*/
1431+
su2double Normal[MAXNDIM] = {0.0}, UnitNormal[MAXNDIM] = {0.0};
14341432

14351433
bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
14361434

1437-
/*--- Allocate arrays ---*/
1438-
Normal = new su2double[nDim];
1439-
NormalArea = new su2double[nDim];
1440-
Jacobian_b = new su2double*[nVar];
1441-
DubDu = new su2double*[nVar];
1442-
u = new su2double[nDim];
1443-
1444-
for (iVar = 0; iVar < nVar; iVar++) {
1445-
Jacobian_b[iVar] = new su2double[nVar];
1446-
DubDu[iVar] = new su2double[nVar];
1447-
}
1448-
1449-
/*--- Get species molar mass ---*/
1450-
auto& Ms = FluidModel->GetSpeciesMolarMass();
1451-
14521435
/*--- Loop over all the vertices on this boundary (val_marker) ---*/
14531436
for (iVertex = 0; iVertex < geometry->nVertex[val_marker]; iVertex++) {
14541437
iPoint = geometry->vertex[val_marker][iVertex]->GetNode();
@@ -1460,15 +1443,14 @@ void CNEMOEulerSolver::BC_Sym_Plane(CGeometry *geometry, CSolver **solver_contai
14601443
geometry->vertex[val_marker][iVertex]->GetNormal(Normal);
14611444

14621445
/*--- Calculate parameters from the geometry ---*/
1463-
Area = GeometryToolbox::Norm(nDim, Normal);
1446+
su2double Area = GeometryToolbox::Norm(nDim, Normal);
14641447

14651448
for (iDim = 0; iDim < nDim; iDim++){
1466-
NormalArea[iDim] = -Normal[iDim];
14671449
UnitNormal[iDim] = -Normal[iDim]/Area;
14681450
}
14691451

14701452
/*--- Retrieve the pressure on the vertex ---*/
1471-
P = nodes->GetPressure(iPoint);
1453+
su2double P = nodes->GetPressure(iPoint);
14721454

14731455
/*--- Apply the flow-tangency b.c. to the convective flux ---*/
14741456
for (iSpecies = 0; iSpecies < nSpecies; iSpecies++)
@@ -1485,25 +1467,30 @@ void CNEMOEulerSolver::BC_Sym_Plane(CGeometry *geometry, CSolver **solver_contai
14851467
/*--- If using implicit time-stepping, calculate b.c. contribution to Jacobian ---*/
14861468
if (implicit) {
14871469

1470+
/*--- Allocate arrays needed for implicit solver ---*/
1471+
su2double Velocity[MAXNDIM] = {0.0};
1472+
1473+
/*--- Get species molar mass ---*/
1474+
auto& Ms = FluidModel->GetSpeciesMolarMass();
1475+
14881476
/*--- Initialize Jacobian ---*/
14891477
for (iVar = 0; iVar < nVar; iVar++)
14901478
for (jVar = 0; jVar < nVar; jVar++)
14911479
Jacobian_i[iVar][jVar] = 0.0;
14921480

14931481
/*--- Calculate state i ---*/
1494-
rho = nodes->GetDensity(iPoint);
1495-
rhoE = nodes->GetSolution(iPoint,nSpecies+nDim);
1496-
rhoEve = nodes->GetSolution(iPoint,nSpecies+nDim+1);
1497-
dPdU = nodes->GetdPdU(iPoint);
1482+
su2double rho = nodes->GetDensity(iPoint);
1483+
su2double rhoE = nodes->GetSolution(iPoint,nSpecies+nDim);
1484+
su2double rhoEve = nodes->GetSolution(iPoint,nSpecies+nDim+1);
1485+
auto dPdU = nodes->GetdPdU(iPoint);
14981486
for (iDim = 0; iDim < nDim; iDim++)
1499-
u[iDim] = nodes->GetVelocity(iPoint,iDim);
1487+
Velocity[iDim] = nodes->GetVelocity(iPoint,iDim);
15001488

1501-
conc = 0.0;
1489+
su2double conc = 0.0;
15021490
for (iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
1503-
cs = nodes->GetMassFraction(iPoint,iSpecies);
1491+
su2double cs = nodes->GetMassFraction(iPoint,iSpecies);
15041492
conc += cs * rho/Ms[iSpecies];
15051493

1506-
/////// NEW //////
15071494
for (iDim = 0; iDim < nDim; iDim++) {
15081495
Jacobian_i[nSpecies+iDim][iSpecies] = dPdU[iSpecies] * UnitNormal[iDim];
15091496
Jacobian_i[iSpecies][nSpecies+iDim] = cs * UnitNormal[iDim];
@@ -1512,7 +1499,7 @@ void CNEMOEulerSolver::BC_Sym_Plane(CGeometry *geometry, CSolver **solver_contai
15121499

15131500
for (iDim = 0; iDim < nDim; iDim++) {
15141501
for (jDim = 0; jDim < nDim; jDim++) {
1515-
Jacobian_i[nSpecies+iDim][nSpecies+jDim] = u[iDim]*UnitNormal[jDim]
1502+
Jacobian_i[nSpecies+iDim][nSpecies+jDim] = Velocity[iDim]*UnitNormal[jDim]
15161503
+ dPdU[nSpecies+jDim]*UnitNormal[iDim];
15171504
}
15181505
Jacobian_i[nSpecies+iDim][nSpecies+nDim] = dPdU[nSpecies+nDim] *UnitNormal[iDim];
@@ -1533,17 +1520,6 @@ void CNEMOEulerSolver::BC_Sym_Plane(CGeometry *geometry, CSolver **solver_contai
15331520
}
15341521
}
15351522
}
1536-
delete [] Normal;
1537-
delete [] NormalArea;
1538-
delete [] u;
1539-
1540-
for (iVar = 0; iVar < nVar; iVar++) {
1541-
delete [] Jacobian_b[iVar];
1542-
delete [] DubDu[iVar];
1543-
}
1544-
1545-
delete [] Jacobian_b;
1546-
delete [] DubDu;
15471523
}
15481524

15491525
void CNEMOEulerSolver::BC_Far_Field(CGeometry *geometry, CSolver **solver_container,
@@ -1553,7 +1529,7 @@ void CNEMOEulerSolver::BC_Far_Field(CGeometry *geometry, CSolver **solver_contai
15531529
unsigned short iDim;
15541530
unsigned long iVertex, iPoint, Point_Normal;
15551531

1556-
su2double *V_infty, *V_domain, *U_domain,*U_infty;
1532+
su2double *V_infty, *V_domain, *U_infty, *U_domain;
15571533

15581534
/*--- Set booleans from configuration parameters ---*/
15591535
bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
@@ -1566,6 +1542,9 @@ void CNEMOEulerSolver::BC_Far_Field(CGeometry *geometry, CSolver **solver_contai
15661542
for (iVertex = 0; iVertex < geometry->nVertex[val_marker]; iVertex++) {
15671543
iPoint = geometry->vertex[val_marker][iVertex]->GetNode();
15681544

1545+
/*--- Allocate the value at the infinity ---*/
1546+
V_infty = GetCharacPrimVar(val_marker, iVertex);
1547+
15691548
/*--- Check if the node belongs to the domain (i.e, not a halo node) ---*/
15701549
if (geometry->nodes->GetDomain(iPoint)) {
15711550

@@ -1954,8 +1933,9 @@ void CNEMOEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container
19541933
bool gravity = config->GetGravityForce();
19551934
bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
19561935

1957-
su2double *U_domain = new su2double[nVar]; su2double *U_outlet = new su2double[nVar];
1958-
su2double *V_domain = new su2double[nPrimVar]; su2double *V_outlet = new su2double[nPrimVar];
1936+
su2double *U_domain;
1937+
su2double *U_outlet = new su2double[nVar];
1938+
su2double *V_domain, *V_outlet;
19591939
su2double *Normal = new su2double[nDim];
19601940
su2double *Ys = new su2double[nSpecies];
19611941

@@ -1973,6 +1953,9 @@ void CNEMOEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container
19731953
for (iVertex = 0; iVertex < geometry->nVertex[val_marker]; iVertex++) {
19741954
iPoint = geometry->vertex[val_marker][iVertex]->GetNode();
19751955

1956+
/*--- Allocate the value at the outlet ---*/
1957+
V_outlet = GetCharacPrimVar(val_marker, iVertex);
1958+
19761959
/*--- Check if the node belongs to the domain (i.e., not a halo node) ---*/
19771960
if (geometry->nodes->GetDomain(iPoint)) {
19781961

@@ -1987,8 +1970,8 @@ void CNEMOEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container
19871970
UnitNormal[iDim] = Normal[iDim]/Area;
19881971

19891972
/*--- Current solution at this boundary node ---*/
1990-
for (iVar = 0; iVar < nVar; iVar++) U_domain[iVar] = nodes->GetSolution(iPoint, iVar);
1991-
for (iVar = 0; iVar < nPrimVar; iVar++) V_domain[iVar] = nodes->GetPrimitive(iPoint, iVar);
1973+
U_domain = nodes->GetSolution(iPoint);
1974+
V_domain = nodes->GetPrimitive(iPoint);
19921975

19931976
/*--- Initialize solution at outlet ---*/
19941977
for (iVar = 0; iVar < nVar; iVar++) U_outlet[iVar] = 0.0;
@@ -2178,10 +2161,7 @@ void CNEMOEulerSolver::BC_Outlet(CGeometry *geometry, CSolver **solver_container
21782161
}
21792162

21802163
/*--- Free locally allocated memory ---*/
2181-
delete [] U_domain;
21822164
delete [] U_outlet;
2183-
delete [] V_domain;
2184-
delete [] V_outlet;
21852165
delete [] Normal;
21862166
delete [] Ys;
21872167

@@ -2406,7 +2386,6 @@ void CNEMOEulerSolver::BC_Supersonic_Outlet(CGeometry *geometry, CSolver **solve
24062386
su2double *U_outlet, *U_domain;
24072387

24082388
bool implicit = (config->GetKind_TimeIntScheme() == EULER_IMPLICIT);
2409-
bool dynamic_grid = config->GetGrid_Movement();
24102389
string Marker_Tag = config->GetMarker_All_TagBound(val_marker);
24112390

24122391
su2double *Normal = new su2double[nDim];
@@ -2427,6 +2406,7 @@ void CNEMOEulerSolver::BC_Supersonic_Outlet(CGeometry *geometry, CSolver **solve
24272406
U_domain = nodes->GetSolution(iPoint);
24282407

24292408
/*--- Allocate the value at the outlet ---*/
2409+
V_outlet = GetCharacPrimVar(val_marker, iVertex);
24302410
V_outlet = V_domain;
24312411
U_outlet = U_domain;
24322412

‎SU2_CFD/src/solvers/CNEMONSSolver.cpp‎

Lines changed: 8 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -96,8 +96,12 @@ void CNEMONSSolver::Preprocessing(CGeometry *geometry, CSolver **solver_containe
9696
SetPrimitive_Limiter(geometry, config);
9797
}
9898

99+
/*--- Compute vorticity and strain mag. ---*/
100+
99101
ComputeVorticityAndStrainMag(*config, iMesh);
100102

103+
/*--- Compute the TauWall from the wall functions ---*/
104+
101105
if (wall_functions) {
102106
SetTau_Wall_WF(geometry, solver_container, config);
103107
}
@@ -319,10 +323,9 @@ void CNEMONSSolver::BC_HeatFluxNonCatalytic_Wall(CGeometry *geometry,
319323
su2double zero[MAXNDIM] = {0.0};
320324
nodes->SetVelocity_Old(iPoint, zero);
321325

322-
for (auto iDim = 0u; iDim < nDim; iDim++){
326+
for (auto iDim = 0u; iDim < nDim; iDim++)
323327
LinSysRes(iPoint, nSpecies+iDim) = 0.0;
324-
nodes->SetVal_ResTruncError_Zero(iPoint,nSpecies+iDim);
325-
}
328+
nodes->SetVel_ResTruncError_Zero(iPoint);
326329

327330
/*--- Apply viscous residual to the linear system ---*/
328331
LinSysRes.SubtractBlock(iPoint, Res_Visc);
@@ -635,10 +638,9 @@ void CNEMONSSolver::BC_IsothermalNonCatalytic_Wall(CGeometry *geometry,
635638
/*--- Initialize viscous residual to zero ---*/
636639
for (auto iVar = 0u; iVar < nVar; iVar ++) {Res_Visc[iVar] = 0.0;}
637640

638-
for (auto iDim = 0u; iDim < nDim; iDim++){
641+
for (auto iDim = 0u; iDim < nDim; iDim++)
639642
LinSysRes(iPoint, nSpecies+iDim) = 0.0;
640-
nodes->SetVal_ResTruncError_Zero(iPoint,nSpecies+iDim);
641-
}
643+
nodes->SetVel_ResTruncError_Zero(iPoint);
642644

643645
/*--- Calculate the gradient of temperature ---*/
644646
su2double Ti = nodes->GetTemperature(iPoint);

0 commit comments

Comments
 (0)