Fixed unused initial_porosity

idaholab · Dec 1, 2022 · 58fe0c8 · 58fe0c8
1 parent 56a65a9
commit 58fe0c8
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/...iemodel/electrothermomechs/simple_geometries/yttria_electrothermomechanical_singleblock.i b/...iemodel/electrothermomechs/simple_geometries/yttria_electrothermomechanical_singleblock.i
@@ -1,5 +1,6 @@
 ## Units in the input file: m-Pa-s-K
 
+initial_porosity = 0.38
 initial_temperature = 300 #roughly 600C where the pyrometer kicks in
 
 [GlobalParams]
@@ -62,8 +63,7 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
   [specific_heat_capacity_va]
     initial_condition = 842.2 # at 1500K #568.73 at 1000K #447.281 # at 293K
   []
-  [density_va]
-    initial_condition = 3106.2 ##5010.0*(1-${initial_porosity}) #in kg/m^3
+  [density_va] # initial condition set in ICs block below
   []
   [heat_transfer_radiation]
   []
@@ -73,6 +73,14 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
   []
 []
 
+[ICs]
+  [density_va_ic]
+    type = FunctionIC
+    variable = density_va
+    function = '5010.0*(1-${initial_porosity})' # 3106.2 in kg/m^3
+  []
+[]
+
 [Modules]
   [TensorMechanics/Master]
     [graphite]