Thermally coupled Constitutive Laws

applications/ConvectionDiffusionApplication/python_scripts/coupled_structural_thermal_solver.py

@@ -67,6 +67,7 @@ def __init__(self, model, custom_settings):
 
         from KratosMultiphysics.ConvectionDiffusionApplication import python_solvers_wrapper_convection_diffusion
         self.thermal_solver = python_solvers_wrapper_convection_diffusion.CreateSolverByParameters(self.model,self.settings["thermal_solver_settings"],"OpenMP")
+        self.is_dynamic = self.settings["structural_solver_settings"]["solver_type"].GetString() == "Dynamic"
 
     def AddVariables(self):
         # Import the structural and thermal solver variables. Then merge them to have them in both structural and thermal solvers.
@@ -149,11 +150,23 @@ def Predict(self):
         self.thermal_solver.Predict()
 
     def SolveSolutionStep(self):
-        fluid_is_converged = self.structural_solver.SolveSolutionStep()
+        KratosMultiphysics.Logger.PrintInfo("\n" + "\t" +"Solving Structural part..." + "\n")
+        solid_is_converged = self.structural_solver.SolveSolutionStep()
+
+        self.RemoveConvectiveVelocity()
+
+        KratosMultiphysics.Logger.PrintInfo("\n" + "\t" + "Solving Thermal part..." + "\n")
         thermal_is_converged = self.thermal_solver.SolveSolutionStep()
 
-        return (fluid_is_converged and thermal_is_converged)
+        KratosMultiphysics.Logger.PrintInfo("+--------------------------------------------------------------+")
+
+        return (solid_is_converged and thermal_is_converged)
 
     def FinalizeSolutionStep(self):
         self.structural_solver.FinalizeSolutionStep()
         self.thermal_solver.FinalizeSolutionStep()
+
+    def RemoveConvectiveVelocity(self):
+        if self.is_dynamic:
+            for node in self.thermal_solver.GetComputingModelPart().Nodes:
+                node.SetSolutionStepValue(KratosMultiphysics.MESH_VELOCITY, node.GetSolutionStepValue(KratosMultiphysics.VELOCITY))