pyvale 2025.5.3__cp311-cp311-manylinux_2_17_i686.manylinux2014_i686.whl

pyvale/simcases/case12.i

@@ -0,0 +1,329 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,monoblock,3mat,therm--mechanical,steady,
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+coolantTemp=100.0      # degC
+heatTransCoeff=125.0e3 # W.m^-2.K^-1
+surfHeatFlux=4.67e6    # W.m^-2, Taken from Adel's first paper
+
+# Material Properties:
+# Thermal Props: Copper-Chromium-Zirconium pg.148 at 200degC
+cucrzrDensity = 8816.0  # kg.m^-3
+cucrzrThermCond = 343.0 # W.m^-1.K^-1
+cucrzrSpecHeat = 407.0  # J.kg^-1.K^-1
+
+# Thermal Props: Pure (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
+
+# Thermal Props: Tungsten at 600degC
+wDensity = 19150.0  # kg.m^-3
+wThermCond = 127.0 # W.m^-1.K^-1
+wSpecHeat = 147.0  # J.kg^-1.K^-1
+
+# Mechanical Props: Copper-Chromium-Zirconium at 200degC
+cucrzrEMod = 123e9       # Pa
+cucrzrPRatio = 0.33      # -
+
+# Mechanical Props: OFHC Copper at 250degC
+cuEMod = 108e9       # Pa
+cuPRatio = 0.33      # -
+
+# Mechanical Props: Tungsten at 250degC
+wEMod = 387e9       # Pa
+wPRatio = 0.29      # -
+
+# Thermo-mechanical coupling
+stressFreeTemp = 20 # degC
+cucrzrThermExp =  17.7e-6 # 1/degC
+cuThermExp = 17.8e-6 # 1/degC
+wThermExp = 4.72e-6 # 1/degC
+
+# Mesh file string
+mesh_file = 'case12.msh'
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[GlobalParams]
+  displacements = 'disp_x disp_y disp_z'
+[]
+
+[Mesh]
+  type = FileMesh
+  file = ${mesh_file}
+[]
+
+[Variables]
+  [temperature]
+    family = LAGRANGE
+    order = SECOND
+    initial_condition = ${coolantTemp}
+  []
+[]
+
+[Kernels]
+  [heat_conduction]
+    type = HeatConduction
+    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 = 'vonmises_stress strain_xx strain_xy strain_xz strain_yx strain_yy strain_yz strain_zx strain_zy strain_zz stress_xx stress_xy stress_xz stress_yx stress_yy stress_yz stress_zx stress_zy stress_zz max_principal_strain mid_principal_strain min_principal_strain'
+  []
+[]
+
+[Materials]
+  [cucrzr_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${cucrzrThermCond}
+    specific_heat = ${cucrzrSpecHeat}
+    block = 'pipe-cucrzr'
+  []
+  [cucrzr_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${cucrzrDensity}
+    block = 'pipe-cucrzr'
+  []
+  [cucrzr_elasticity]
+    type = ComputeIsotropicElasticityTensor
+    youngs_modulus = ${cucrzrEMod}
+    poissons_ratio = ${cucrzrPRatio}
+    block = 'pipe-cucrzr'
+  []
+  [cucrzr_expansion]
+    type = ComputeThermalExpansionEigenstrain
+    temperature = temperature
+    stress_free_temperature = ${stressFreeTemp}
+    thermal_expansion_coeff = ${cucrzrThermExp}
+    eigenstrain_name = thermal_expansion_eigenstrain
+    block = 'pipe-cucrzr'
+  []
+
+
+  [copper_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${cuThermCond}
+    specific_heat = ${cuSpecHeat}
+    block = 'interlayer-cu'
+  []
+  [copper_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${cuDensity}
+    block = 'interlayer-cu'
+  []
+  [copper_elasticity]
+    type = ComputeIsotropicElasticityTensor
+    youngs_modulus = ${cuEMod}
+    poissons_ratio = ${cuPRatio}
+    block = 'interlayer-cu'
+  []
+  [copper_expansion]
+    type = ComputeThermalExpansionEigenstrain
+    temperature = temperature
+    stress_free_temperature = ${stressFreeTemp}
+    thermal_expansion_coeff = ${cuThermExp}
+    eigenstrain_name = thermal_expansion_eigenstrain
+    block = 'interlayer-cu'
+  []
+
+
+  [tungsten_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${wThermCond}
+    specific_heat = ${wSpecHeat}
+    block = 'armour-w'
+  []
+  [tungsten_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${wDensity}
+    block = 'armour-w'
+  []
+  [tungsten_elasticity]
+    type = ComputeIsotropicElasticityTensor
+    youngs_modulus = ${wEMod}
+    poissons_ratio = ${wPRatio}
+    block = 'armour-w'
+  []
+  [tungsten_expansion]
+    type = ComputeThermalExpansionEigenstrain
+    temperature = temperature
+    stress_free_temperature = ${stressFreeTemp}
+    thermal_expansion_coeff = ${wThermExp}
+    eigenstrain_name = thermal_expansion_eigenstrain
+    block = 'armour-w'
+  []
+
+  [stress]
+    type = ComputeFiniteStrainElasticStress # ComputeLinearElasticStress or ComputeFiniteStrainElasticStress
+  []
+[]
+
+[BCs]
+  [heat_flux_in]
+    type = NeumannBC
+    variable = temperature
+    boundary = 'bc-top-heatflux'
+    value = ${surfHeatFlux}
+  []
+  [heat_flux_out]
+    type = ConvectiveHeatFluxBC
+    variable = temperature
+    boundary = 'bc-pipe-heattransf'
+    T_infinity = ${coolantTemp}
+    heat_transfer_coefficient = ${heatTransCoeff}
+  []
+
+  # Lock disp_y for base
+    # NOTE: if locking y on base need to comment all disp_y conditions below
+    [mech_bc_c_dispy]
+      type = DirichletBC
+     variable = disp_y
+      boundary = 'bc-base-disp'
+      value = 0.0
+  []
+
+  # Lock all disp DOFs at the center of the block
+  [mech_bc_c_dispx]
+      type = DirichletBC
+      variable = disp_x
+      boundary = 'bc-c-point-xyz-mech'
+      value = 0.0
+  []
+  #[mech_bc_c_dispy]
+  #    type = DirichletBC
+  #    variable = disp_y
+  #    boundary = 'bc-c-point-xyz-mech'
+  #    value = 0.0
+  #[]
+  [mech_bc_c_dispz]
+      type = DirichletBC
+      variable = disp_z
+      boundary = 'bc-c-point-xyz-mech'
+      value = 0.0
+  []
+
+  # Lock disp yz along the x (left-right) axis
+  #[mech_bc_l_dispy]
+  #    type = DirichletBC
+  #    variable = disp_y
+  #    boundary = 'bc-l-point-yz-mech'
+  #    value = 0.0
+  #[]
+  [mech_bc_l_dispz]
+      type = DirichletBC
+      variable = disp_z
+      boundary = 'bc-l-point-yz-mech'
+      value = 0.0
+  []
+  #[mech_bc_r_dispy]
+  #    type = DirichletBC
+  #    variable = disp_y
+  #    boundary = 'bc-r-point-yz-mech'
+  #    value = 0.0
+  #[]
+  [mech_bc_r_dispz]
+      type = DirichletBC
+      variable = disp_z
+      boundary = 'bc-r-point-yz-mech'
+      value = 0.0
+  []
+
+  # Lock disp xy along the z (front-back) axis
+  [mech_bc_f_dispx]
+      type = DirichletBC
+      variable = disp_x
+      boundary = 'bc-f-point-xy-mech'
+      value = 0.0
+  []
+  #[mech_bc_f_dispy]
+  #    type = DirichletBC
+  #    variable = disp_y
+  #    boundary = 'bc-f-point-xy-mech'
+  #    value = 0.0
+  #[]
+  [mech_bc_b_dispx]
+      type = DirichletBC
+      variable = disp_x
+      boundary = 'bc-b-point-xy-mech'
+      value = 0.0
+  []
+  #[mech_bc_b_dispy]
+  #    type = DirichletBC
+  #    variable = disp_y
+  #    boundary = 'bc-b-point-xy-mech'
+  #    value = 0.0
+  #[]
+[]
+
+[Preconditioning]
+  [smp]
+    type = SMP
+    full = true
+  []
+[]
+
+[Executioner]
+  type = Steady
+  solve_type = 'PJFNK'
+  petsc_options_iname = '-pc_type -pc_hypre_type'
+  petsc_options_value = 'hypre    boomeramg'
+[]
+
+
+[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
+  []
+  [disp_z_max]
+    type = NodalExtremeValue
+    variable = disp_z
+  []
+
+  [strain_xx_max]
+      type = ElementExtremeValue
+      variable = strain_xx
+  []
+  [strain_yy_max]
+      type = ElementExtremeValue
+      variable = strain_yy
+  []
+  [strain_zz_max]
+    type = ElementExtremeValue
+    variable = strain_zz
+  []
+[]
+
+
+[Outputs]
+  exodus = true
+[]

pyvale/simcases/case13.i

@@ -0,0 +1,140 @@
+#-------------------------------------------------------------------------
+# pyvale: simple,2Dplate,1mat,thermal,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+endTime = 200
+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   # s
+
+# Material Properties: Pure (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
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[Mesh]
+    [generated]
+        type = GeneratedMeshGenerator
+        dim = 2
+        nx = ${nElemX}
+        ny = ${nElemY}
+        xmax = ${lengX}
+        ymax = ${lengY}
+        elem_type = ${eType}
+    []
+[]
+
+[Variables]
+    [temperature]
+        initial_condition = ${coolantTemp}
+    []
+[]
+
+[Kernels]
+    [heat_conduction]
+        type = HeatConduction
+        variable = temperature
+    []
+    [time_derivative]
+        type = HeatConductionTimeDerivative
+        variable = temperature
+    []
+[]
+
+[Materials]
+    [copper_thermal]
+        type = HeatConductionMaterial
+        thermal_conductivity = ${cuThermCond}
+        specific_heat = ${cuSpecHeat}
+    []
+    [copper_density]
+        type = GenericConstantMaterial
+        prop_names = 'density'
+        prop_values = ${cuDensity}
+    []
+[]
+
+[BCs]
+    [heat_flux_out_left]
+        type = ConvectiveHeatFluxBC
+        variable = temperature
+        boundary = 'left'
+        T_infinity = ${coolantTemp}
+        heat_transfer_coefficient = ${heatTransCoeff}
+    []
+    [heat_flux_in_right]
+        type = FunctionNeumannBC
+        variable = temperature
+        boundary = 'right'
+        function = '${fparse surfHeatFlux}*(1-exp(-(1/${timeConst})*t))'
+    []
+    [heat_flux_out_bot]
+        type = ConvectiveHeatFluxBC
+        variable = temperature
+        boundary = 'bottom'
+        T_infinity = ${coolantTemp}
+        heat_transfer_coefficient = ${heatTransCoeff}
+    []
+    [heat_flux_in_top]
+        type = FunctionNeumannBC
+        variable = temperature
+        boundary = 'top'
+        function = '${fparse surfHeatFlux}*(1-exp(-(1/${timeConst})*t))'
+    []
+[]
+
+[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 # TODO: check this helps
+    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]
+    [max_temp]
+        type = NodalExtremeValue
+        variable = temperature
+    []
+    [avg_temp]
+        type = AverageNodalVariableValue
+        variable = temperature
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]

pyvale/simcases/case14.i

@@ -0,0 +1,159 @@
+#-------------------------------------------------------------------------
+# pyvale: simple,2Dplate,mechanical,steady,
+#-------------------------------------------------------------------------
+# NOTE: default 2D MOOSE solid mechanics is plane strain
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+# NOTE: only used for transient solves
+endTime = 60    # s
+timeStep = 5    # s
+
+# Geometric Properties
+lengX = 50e-3  # m
+lengY = 100e-3   # m
+
+# Mesh Properties
+nElemX = 10
+nElemY = 20
+eType = QUAD8 # QUAD4 for 1st order, QUAD8 for 2nd order
+
+# Mechanical Loads/BCs
+topDispRate = ${fparse 1e-3 / 60}  # m/s
+
+# Material Properties: OFHC Copper 250degC
+cuEMod= 108e9   # Pa
+cuPRatio = 0.33     # -
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[GlobalParams]
+    displacements = 'disp_x disp_y'
+[]
+
+[Mesh]
+    [generated]
+        type = GeneratedMeshGenerator
+        dim = 2
+        nx = ${nElemX}
+        ny = ${nElemY}
+        xmax = ${lengX}
+        ymax = ${lengY}
+        elem_type = ${eType}
+    []
+[]
+
+[Modules/TensorMechanics/Master]
+    [all]
+        strain = SMALL
+        incremental = true
+        add_variables = true
+        material_output_family = MONOMIAL   # MONOMIAL, LAGRANGE
+        material_output_order = FIRST       # CONSTANT, FIRST, SECOND,
+        generate_output = 'vonmises_stress stress_xx stress_yy stress_xy strain_xx strain_yy strain_xy max_principal_strain mid_principal_strain min_principal_strain'
+    []
+[]
+
+[BCs]
+    [bottom_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bottom'
+        value = 0.0
+    []
+    [bottom_y]
+        type = DirichletBC
+        variable = disp_y
+        boundary = 'bottom'
+        value = 0.0
+    []
+    [top_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'top'
+        value = 0.0
+    []
+    [top_y]
+        type = FunctionDirichletBC
+        variable = disp_y
+        boundary = 'top'
+        function = '${topDispRate}*t'
+    []
+[]
+
+[Materials]
+    [elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${cuEMod}
+        poissons_ratio = ${cuPRatio}
+    []
+    [stress]
+        type = ComputeFiniteStrainElasticStress # ComputeLinearElasticStress
+    []
+[]
+
+[Preconditioning]
+    [SMP]
+        type = SMP
+        full = true
+    []
+[]
+
+[Executioner]
+    type = Transient
+    solve_type = 'PJFNK'
+    petsc_options_iname = '-pc_type -pc_hypre_type'
+    petsc_options_value = 'hypre boomeramg'
+    end_time= ${endTime}
+    dt = ${timeStep}
+[]
+
+
+[Postprocessors]
+    [react_y_bot]
+        type = SidesetReaction
+        direction = '0 1 0'
+        stress_tensor = stress
+        boundary = 'bottom'
+    []
+    [react_y_top]
+        type = SidesetReaction
+        direction = '0 1 0'
+        stress_tensor = stress
+        boundary = 'top'
+    []
+
+    [disp_y_max]
+        type = NodalExtremeValue
+        variable = disp_y
+    []
+    [disp_x_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+
+    [max_yy_stress]
+        type = ElementExtremeValue
+        variable = stress_yy
+    []
+
+    [strain_yy_avg]
+        type = ElementAverageValue
+        variable = strain_yy
+    []
+    [strain_xx_avg]
+        type = ElementAverageValue
+        variable = strain_xx
+    []
+
+    [stress_vm_max]
+        type = ElementExtremeValue
+        variable = vonmises_stress
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]