pyvale 2025.5.3__cp311-cp311-win32.whl

pyvale/simcases/case18.i

@@ -0,0 +1,254 @@
+#-------------------------------------------------------------------------
+# pyvale: simple,2Dplate,1mat,thermomechanical,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+endTime = 300
+timeStep = 5
+
+# Geometric Properties
+lengX = 100e-3  # m
+lengY = 50e-3   # m
+
+# Mesh Properties
+nElemX = 20
+nElemY = 10
+eType = QUAD4 # QUAD4 for 1st order, QUAD8 for 2nd order
+
+# Thermal Loads/BCs
+coolantTemp = 20.0      # degC
+heatTransCoeff = 125.0e3 # W.m^-2.K^-1
+surfHeatFlux = 500.0e3    # W.m^-2
+timeConst = 1
+
+# Material Properties:
+# Thermal Props:OFHC) Copper at 250degC
+cuDensity = 8829.0  # kg.m^-3
+cuThermCond = 384.0 # W.m^-1.K^-1
+cuSpecHeat = 406.0  # J.kg^-1.K^-1
+
+# Mechanical Props: OFHC Copper 250degC
+cuEMod = 108e9       # Pa
+cuPRatio = 0.33      # -
+
+# Thermo-mechanical coupling
+stressFreeTemp = 20 # degC
+cuThermExp = 17.8e-6 # 1/degC
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+# Selection tolerance for bounding box node selection
+sTol = ${fparse lengX/(nElemX*4)}
+
+[GlobalParams]
+    displacements = 'disp_x disp_y'
+[]
+
+[Mesh]
+    [generated_mesh]
+        type = GeneratedMeshGenerator
+        dim = 2
+        nx = ${nElemX}
+        ny = ${nElemY}
+        xmax = ${lengX}
+        ymax = ${lengY}
+        elem_type = ${eType}
+    []
+
+    [node_bottom_left]
+        type = BoundingBoxNodeSetGenerator
+        input = generated_mesh
+        bottom_left = '${fparse 0-sTol}
+                       ${fparse 0-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse 0+sTol}
+                     ${fparse 0+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = bottom_left_node
+    []
+    [node_top_left]
+        type = BoundingBoxNodeSetGenerator
+        input = node_bottom_left
+        bottom_left = '${fparse 0-sTol}
+                       ${fparse lengY-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse 0+sTol}
+                     ${fparse lengY+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = top_left_node
+    []
+    [node_bottom_right]
+        type = BoundingBoxNodeSetGenerator
+        input = node_top_left
+        bottom_left = '${fparse lengX-sTol}
+                       ${fparse 0-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse lengX+sTol}
+                     ${fparse 0+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = bottom_right_node
+    []
+[]
+
+[Variables]
+    [temperature]
+        family = LAGRANGE
+        order = FIRST
+        initial_condition = ${coolantTemp}
+    []
+[]
+
+[Kernels]
+    [heat_conduction]
+        type = HeatConduction
+        variable = temperature
+    []
+    [time_derivative]
+        type = HeatConductionTimeDerivative
+        variable = temperature
+    []
+[]
+
+[Modules/TensorMechanics/Master]
+    [all]
+        strain = SMALL                      # SMALL or FINITE
+        incremental = true
+        add_variables = true
+        material_output_family = MONOMIAL   # MONOMIAL, LAGRANGE
+        material_output_order = FIRST       # CONSTANT, FIRST, SECOND,
+        automatic_eigenstrain_names = true
+        generate_output = 'strain_xx strain_xy strain_xz strain_yx strain_yy strain_yz strain_zx strain_zy strain_zz'
+    []
+[]
+
+[Materials]
+    [copper_thermal]
+        type = HeatConductionMaterial
+        thermal_conductivity = ${cuThermCond}
+        specific_heat = ${cuSpecHeat}
+    []
+    [copper_density]
+        type = GenericConstantMaterial
+        prop_names = 'density'
+        prop_values = ${cuDensity}
+    []
+    [copper_elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${cuEMod}
+        poissons_ratio = ${cuPRatio}
+    []
+    [copper_expansion]
+        type = ComputeThermalExpansionEigenstrain
+        temperature = temperature
+        stress_free_temperature = ${stressFreeTemp}
+        thermal_expansion_coeff = ${cuThermExp}
+        eigenstrain_name = thermal_expansion_eigenstrain
+    []
+    [stress]
+        type = ComputeFiniteStrainElasticStress # ComputeLinearElasticStress or ComputeFiniteStrainElasticStress
+    []
+[]
+
+[BCs]
+    [heat_flux_out]
+        type = ConvectiveHeatFluxBC
+        variable = temperature
+        boundary = 'left'
+        T_infinity = ${coolantTemp}
+        heat_transfer_coefficient = ${heatTransCoeff}
+    []
+    [heat_flux_in]
+        type = FunctionNeumannBC
+        variable = temperature
+        boundary = 'right'
+        function = '${fparse surfHeatFlux}*(1-exp(-(1/${timeConst})*t))'
+    []
+
+    [left_disp_y]
+        type = DirichletBC
+        variable = disp_y
+        boundary = 'left'
+        value = 0.0
+    []
+    [left_disp_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'left'
+        value = 0.0
+    []
+[]
+
+[Preconditioning]
+    [smp]
+        type = SMP
+        full = true
+    []
+[]
+
+[Executioner]
+    type = Transient
+
+    solve_type = PJFNK   # PJNFK or NEWTON
+    l_max_its = 100       # default = 1000
+    l_tol = 1e-6          # default = 1e-5
+    nl_abs_tol = 1e-6     # default = 1e-50
+    nl_rel_tol = 1e-6     # default = 1e-8
+
+    line_search = none
+    petsc_options_iname = '-pc_type -pc_hypre_type'
+    petsc_options_value = 'hypre boomeramg'
+
+    start_time = 0.0
+    end_time = ${endTime}
+    dt = ${timeStep}
+
+    [Predictor]
+      type = SimplePredictor
+      scale = 1
+    []
+[]
+
+
+[Postprocessors]
+    [temp_max]
+        type = NodalExtremeValue
+        variable = temperature
+    []
+    [temp_avg]
+        type = AverageNodalVariableValue
+        variable = temperature
+    []
+
+    [disp_x_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+    [disp_y_max]
+        type = NodalExtremeValue
+        variable = disp_y
+    []
+
+    [strain_xx_max]
+        type = ElementExtremeValue
+        variable = strain_xx
+    []
+    [strain_yy_max]
+        type = ElementExtremeValue
+        variable = strain_yy
+    []
+
+    [strain_xx_avg]
+        type = ElementAverageValue
+        variable = strain_xx
+    []
+    [strain_yy_avg]
+        type = ElementAverageValue
+        variable = strain_yy
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]

pyvale/simcases/case18_1.i

@@ -0,0 +1,254 @@
+#-------------------------------------------------------------------------
+# pyvale: simple,2Dplate,1mat,thermomechanical,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+endTime = 300
+timeStep = 5
+
+# Geometric Properties
+lengX = 100e-3  # m
+lengY = 50e-3   # m
+
+# Mesh Properties
+nElemX = 20
+nElemY = 10
+eType = QUAD4 # QUAD4 for 1st order, QUAD8 for 2nd order
+
+# Thermal Loads/BCs
+coolantTemp = 20.0      # degC
+heatTransCoeff = 125.0e3 # W.m^-2.K^-1
+surfHeatFlux = 500.0e3    # W.m^-2
+timeConst = 1
+
+# Material Properties:
+# Thermal Props:OFHC) Copper at 250degC
+cuDensity = 8829.0  # kg.m^-3
+cuThermCond = 384.0 # W.m^-1.K^-1
+cuSpecHeat = 406.0  # J.kg^-1.K^-1
+
+# Mechanical Props: OFHC Copper 250degC
+cuEMod = 108e9       # Pa
+cuPRatio = 0.33      # -
+
+# Thermo-mechanical coupling
+stressFreeTemp = 20 # degC
+cuThermExp = 17.8e-6 # 1/degC
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+# Selection tolerance for bounding box node selection
+sTol = ${fparse lengX/(nElemX*4)}
+
+[GlobalParams]
+    displacements = 'disp_x disp_y'
+[]
+
+[Mesh]
+    [generated_mesh]
+        type = GeneratedMeshGenerator
+        dim = 2
+        nx = ${nElemX}
+        ny = ${nElemY}
+        xmax = ${lengX}
+        ymax = ${lengY}
+        elem_type = ${eType}
+    []
+
+    [node_bottom_left]
+        type = BoundingBoxNodeSetGenerator
+        input = generated_mesh
+        bottom_left = '${fparse 0-sTol}
+                       ${fparse 0-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse 0+sTol}
+                     ${fparse 0+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = bottom_left_node
+    []
+    [node_top_left]
+        type = BoundingBoxNodeSetGenerator
+        input = node_bottom_left
+        bottom_left = '${fparse 0-sTol}
+                       ${fparse lengY-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse 0+sTol}
+                     ${fparse lengY+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = top_left_node
+    []
+    [node_bottom_right]
+        type = BoundingBoxNodeSetGenerator
+        input = node_top_left
+        bottom_left = '${fparse lengX-sTol}
+                       ${fparse 0-sTol}
+                       ${fparse 0-sTol}'
+        top_right = '${fparse lengX+sTol}
+                     ${fparse 0+sTol}
+                     ${fparse 0+sTol}'
+        new_boundary = bottom_right_node
+    []
+[]
+
+[Variables]
+    [temperature]
+        family = LAGRANGE
+        order = FIRST
+        initial_condition = ${coolantTemp}
+    []
+[]
+
+[Kernels]
+    [heat_conduction]
+        type = HeatConduction
+        variable = temperature
+    []
+    [time_derivative]
+        type = HeatConductionTimeDerivative
+        variable = temperature
+    []
+[]
+
+[Modules/TensorMechanics/Master]
+    [all]
+        strain = SMALL                      # SMALL or FINITE
+        incremental = true
+        add_variables = true
+        material_output_family = MONOMIAL   # MONOMIAL, LAGRANGE
+        material_output_order = FIRST       # CONSTANT, FIRST, SECOND,
+        automatic_eigenstrain_names = true
+        generate_output = 'strain_xx strain_xy strain_xz strain_yx strain_yy strain_yz strain_zx strain_zy strain_zz'
+    []
+[]
+
+[Materials]
+    [copper_thermal]
+        type = HeatConductionMaterial
+        thermal_conductivity = ${cuThermCond}
+        specific_heat = ${cuSpecHeat}
+    []
+    [copper_density]
+        type = GenericConstantMaterial
+        prop_names = 'density'
+        prop_values = ${cuDensity}
+    []
+    [copper_elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${cuEMod}
+        poissons_ratio = ${cuPRatio}
+    []
+    [copper_expansion]
+        type = ComputeThermalExpansionEigenstrain
+        temperature = temperature
+        stress_free_temperature = ${stressFreeTemp}
+        thermal_expansion_coeff = ${cuThermExp}
+        eigenstrain_name = thermal_expansion_eigenstrain
+    []
+    [stress]
+        type = ComputeFiniteStrainElasticStress # ComputeLinearElasticStress or ComputeFiniteStrainElasticStress
+    []
+[]
+
+[BCs]
+    [heat_flux_out]
+        type = ConvectiveHeatFluxBC
+        variable = temperature
+        boundary = 'left'
+        T_infinity = ${coolantTemp}
+        heat_transfer_coefficient = ${heatTransCoeff}
+    []
+    [heat_flux_in]
+        type = FunctionNeumannBC
+        variable = temperature
+        boundary = 'right'
+        function = '${fparse surfHeatFlux}*(1-exp(-(1/${timeConst})*t))'
+    []
+
+    [left_disp_y]
+        type = DirichletBC
+        variable = disp_y
+        boundary = 'left'
+        value = 0.0
+    []
+    [left_disp_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'left'
+        value = 0.0
+    []
+[]
+
+[Preconditioning]
+    [smp]
+        type = SMP
+        full = true
+    []
+[]
+
+[Executioner]
+    type = Transient
+
+    solve_type = PJFNK   # PJNFK or NEWTON
+    l_max_its = 100       # default = 1000
+    l_tol = 1e-6          # default = 1e-5
+    nl_abs_tol = 1e-6     # default = 1e-50
+    nl_rel_tol = 1e-6     # default = 1e-8
+
+    line_search = none
+    petsc_options_iname = '-pc_type -pc_hypre_type'
+    petsc_options_value = 'hypre boomeramg'
+
+    start_time = 0.0
+    end_time = ${endTime}
+    dt = ${timeStep}
+
+    [Predictor]
+      type = SimplePredictor
+      scale = 1
+    []
+[]
+
+
+[Postprocessors]
+    [temp_max]
+        type = NodalExtremeValue
+        variable = temperature
+    []
+    [temp_avg]
+        type = AverageNodalVariableValue
+        variable = temperature
+    []
+
+    [disp_x_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+    [disp_y_max]
+        type = NodalExtremeValue
+        variable = disp_y
+    []
+
+    [strain_xx_max]
+        type = ElementExtremeValue
+        variable = strain_xx
+    []
+    [strain_yy_max]
+        type = ElementExtremeValue
+        variable = strain_yy
+    []
+
+    [strain_xx_avg]
+        type = ElementAverageValue
+        variable = strain_xx
+    []
+    [strain_yy_avg]
+        type = ElementAverageValue
+        variable = strain_yy
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]