fixes for simple Elasticity + Heat

Author: pcarruscag

SU2_CFD/src/iteration/CDiscAdjFEAIteration.cpp

@@ -296,7 +296,7 @@ void CDiscAdjFEAIteration::SetDependencies(CSolver***** solver, CGeometry**** ge
   auto structural_numerics = numerics[iZone][iInst][MESH_0][FEA_SOL];
 
   /*--- Some numerics are only instanciated under these conditions ---*/
-  bool fsi = config[iZone]->GetFSI_Simulation();
+  bool fsi = config[iZone]->GetFSI_Simulation() || config[iZone]->GetMultizone_Problem();
   bool nonlinear = config[iZone]->GetGeometricConditions() == LARGE_DEFORMATIONS;
   bool de_effects = config[iZone]->GetDE_Effects() && nonlinear;
   bool element_based = dir_solver->IsElementBased() && nonlinear;

SU2_CFD/src/variables/CFEAVariable.cpp

@@ -6,7 +6,7 @@
  *
  * SU2 Project Website: https://su2code.github.io
  *
- * The SU2 Project is maintained by the SU2 Foundation 
+ * The SU2 Project is maintained by the SU2 Foundation
  * (http://su2foundation.org)
  *
  * Copyright 2012-2020, SU2 Contributors (cf. AUTHORS.md)
@@ -39,7 +39,8 @@ CFEAVariable::CFEAVariable(const su2double *val_fea, unsigned long npoint, unsig
   bool discrete_adjoint   = config->GetDiscrete_Adjoint();
   bool refgeom            = config->GetRefGeom(); // Reference geometry needs to be stored
   bool dynamic_analysis   = config->GetTime_Domain();
-  bool fsi_analysis       = config->GetFSI_Simulation();
+  bool multizone          = config->GetMultizone_Problem();
+  bool fsi_analysis       = config->GetFSI_Simulation() || multizone;
 
   VonMises_Stress.resize(nPoint) = su2double(0.0);
 
@@ -85,8 +86,7 @@ CFEAVariable::CFEAVariable(const su2double *val_fea, unsigned long npoint, unsig
 
   if (prestretch_fem) Prestretch.resize(nPoint,nVar);
 
-  if (config->GetMultizone_Problem())
-    Set_BGSSolution_k();
+  if (multizone) Set_BGSSolution_k();
 }
 
 void CFEAVariable::SetSolution_Vel_time_n() { Solution_Vel_time_n = Solution_Vel; }