From dfa7766b5b7075296d4d43eaea66eef84da796a2 Mon Sep 17 00:00:00 2001 From: Uwe Date: Mon, 27 Feb 2023 14:54:58 +0100 Subject: [PATCH] [FEM] [skip ci] fixes for flow example - the heat solver settings were wrong so that the heat was actually not calculated for convection - also set more realistic temperatures --- src/Mod/Fem/femexamples/equation_flow_elmer_2D.py | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/src/Mod/Fem/femexamples/equation_flow_elmer_2D.py b/src/Mod/Fem/femexamples/equation_flow_elmer_2D.py index 9879b79b09..32d09e317d 100644 --- a/src/Mod/Fem/femexamples/equation_flow_elmer_2D.py +++ b/src/Mod/Fem/femexamples/equation_flow_elmer_2D.py @@ -133,10 +133,12 @@ def setup(doc=None, solvertype="elmer"): equation_flow.IdrsParameter = 3 equation_flow.LinearIterativeMethod = "Idrs" equation_flow.LinearPreconditioning = "ILU1" + equation_heat.Convection = "Computed" equation_heat.IdrsParameter = 3 equation_heat.LinearIterativeMethod = "Idrs" equation_heat.LinearPreconditioning = "ILU1" - equation_heat.RelaxationFactor = 0.9 + equation_heat.Priority = 5 + equation_heat.Stabilize = True # material @@ -220,7 +222,7 @@ def setup(doc=None, solvertype="elmer"): # constraint wall temperature Temperature_Wall = ObjectsFem.makeConstraintTemperature(doc, "Temperature_Wall") - Temperature_Wall.Temperature = 293.0 + Temperature_Wall.Temperature = 300.0 Temperature_Wall.NormalDirection = Vector(0, 0, -1) Temperature_Wall.References = [ (BooleanFragments, "Edge2"), @@ -238,7 +240,7 @@ def setup(doc=None, solvertype="elmer"): # constraint heating rod temperature Temperature_HeatingRod = ObjectsFem.makeConstraintTemperature(doc, "Temperature_HeatingRod") - Temperature_HeatingRod.Temperature = 600.0 + Temperature_HeatingRod.Temperature = 373.0 Temperature_HeatingRod.NormalDirection = Vector(0, -1, 0) Temperature_HeatingRod.References = [(BooleanFragments, "Edge1")] analysis.addObject(Temperature_HeatingRod)