pyvale 2025.5.3__cp311-cp311-win_amd64.whl

pyvale/simcases/case20.i

@@ -0,0 +1,250 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,3Dstcgmsh,1mat,thermomechanical,steady,
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+# Thermal Loads/BCs
+coolantTemp = 150.0      # degC
+heatTransCoeff = 125.0e3 # W.m^-2.K^-1
+
+surfHeatPower = 0.25e3     # W
+blockLeng = 49.5e-3
+blockWidth = 37e-3
+surfArea = ${fparse blockLeng*blockWidth}   # m^2
+surfHeatFlux = ${fparse surfHeatPower/surfArea} # W.m^-2
+
+# Material Properties:
+# Thermal Props: SS316L @ 400degC
+ss316LDensity = 7770.0  # kg.m^-3
+ss316LThermCond = 19.99 # W.m^-1.K^-1
+ss316LSpecHeat = 556.0  # J.kg^-1.K^-1
+
+# Mechanical Props: SS316L @ 400degC
+ss316LEMod = 168e9       # Pa
+ss316LPRatio = 0.3      # -
+
+# Thermo-mechanical coupling
+stressFreeTemp = 150 # degC
+ss316LThermExp = 17.8e-6 # 1/degC
+
+# Element properties
+elem_order = 'SECOND'
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+
+[GlobalParams]
+    displacements = 'disp_x disp_y disp_z'
+[]
+
+[Mesh]
+    type = FileMesh
+    file = 'case20.msh'
+[]
+
+[Variables]
+    [temperature]
+        family = LAGRANGE
+        order = ${elem_order}
+        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]
+    [copper_thermal]
+        type = HeatConductionMaterial
+        thermal_conductivity = ${ss316LThermCond}
+        specific_heat = ${ss316LSpecHeat}
+    []
+    [copper_density]
+        type = GenericConstantMaterial
+        prop_names = 'density'
+        prop_values = ${ss316LDensity}
+    []
+    [copper_elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${ss316LEMod}
+        poissons_ratio = ${ss316LPRatio}
+    []
+    [copper_expansion]
+        type = ComputeThermalExpansionEigenstrain
+        temperature = temperature
+        stress_free_temperature = ${stressFreeTemp}
+        thermal_expansion_coeff = ${ss316LThermExp}
+        eigenstrain_name = thermal_expansion_eigenstrain
+    []
+
+    [stress]
+        type = ComputeFiniteStrainElasticStress # ComputeLinearElasticStress or ComputeFiniteStrainElasticStress
+    []
+[]
+
+[BCs]
+    [heat_flux_out]
+        type = ConvectiveHeatFluxBC
+        variable = temperature
+        boundary = 'bc-pipe-htc'
+        T_infinity = ${coolantTemp}
+        heat_transfer_coefficient = ${heatTransCoeff}
+    []
+    [heat_flux_in]
+        type = NeumannBC
+        variable = temperature
+        boundary = 'bc-top-heatflux'
+        value = ${surfHeatFlux}
+    []
+
+    # Lock disp_y for whole base
+    [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-base-c-loc-xyz'
+        value = 0.0
+    []
+    [mech_bc_c_dispz]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-base-c-loc-xyz'
+        value = 0.0
+    []
+
+    # Lock z dof along x axis
+    [mech_bc_px_dispz]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-base-nx-loc-z'
+        value = 0.0
+    []
+
+    # Lock x dof along z
+    [mech_bc_pz_dispx]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-base-pz-loc-x'
+        value = 0.0
+    []
+    [mech_bc_nz_dispx]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-base-nz-loc-x'
+        value = 0.0
+    []
+[]
+
+[Preconditioning]
+   [smp]
+       type = SMP
+       full = true
+   []
+[]
+
+# LF-PersonalLaptop AMD 8 core / 8 thread
+# Trans, Precon=ON, NEWTON, pctype=lu,  solve time with 7 mpi tasks = 229.18s
+# Trans, Precon=OFF, NEWTON, pctype=lu,  solve time with 7 mpi tasks = 226.52s
+# Steady, Precon=OFF, NEWTON, pctype=lu,  solve time with 7 mpi tasks = 226.52s
+
+# LF-WorkLaptop AMD 8 core/ 16 threads
+# Steady, Precon=OFF, NEWTON, pctype=lu,  solve time with 8 mpi tasks = 275s
+# Steady, Precon=OFF, NEWTON, pctype=lu,  solve time with 6 mpi tasks = ~350s
+# Steady, Precon=OFF, NEWTON, pctype=lu,  solve time with 4 mpi tasks = 458s
+
+[Executioner]
+    type = Steady
+    #---------------------------------------------------------------------------
+    solve_type = PJFNK   # PJNFK or NEWTON
+    l_max_its = 100       # default = 1000
+    nl_max_its = 100
+    l_tol = 1e-5          # default = 1e-5
+    nl_abs_tol = 1e-6     # default = 1e-50, set 1e-6
+    nl_rel_tol = 1e-6     # default = 1e-8, set 1e-6
+
+    line_search = none
+    petsc_options_iname = '-pc_type -pc_hypre_type'
+    petsc_options_value = 'hypre boomeramg'
+
+    #---------------------------------------------------------------------------
+    # solve_type = 'NEWTON' # NEWTON or PJNFK
+    # petsc_options_iname = '-pc_type'
+    # petsc_options_value = 'lu'
+
+    # l_max_its = 100
+    # nl_max_its = 100
+    # nl_rel_tol = 1e-6
+    # nl_abs_tol = 1e-6
+    # l_tol = 1e-5
+
+    #---------------------------------------------------------------------------
+    #end_time= ${endTime}
+    #dt = ${timeStep}
+[]
+
+[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/case21.geo

@@ -0,0 +1,74 @@
+//==============================================================================
+// Gmsh 3D cylinder imaging test case
+// author: Lloyd Fletcher (scepticalrabbit)
+//==============================================================================
+// Always set to OpenCASCADE - circles and boolean opts are much easier!
+SetFactory("OpenCASCADE");
+
+// Allows gmsh to print to terminal in vscode - easier debugging
+General.Terminal = 1;
+
+// View options - not required when
+Geometry.PointLabels = 1;
+Geometry.CurveLabels = 1;
+Geometry.SurfaceLabels = 1;
+Geometry.VolumeLabels = 0;
+
+//-------------------------------------------------------------------------
+//_* MOOSEHERDER VARIABLES - START
+file_name = "case21.msh";
+
+// Geometric variables
+cyl_height = 25e-3;
+cyl_diam = 25e-3;
+
+// Must be an integer
+elem_order = 1;
+mesh_ref = 4;
+mesh_size = 2.5e-3/mesh_ref;
+num_threads = 4;
+//** MOOSEHERDER VARIABLES - END
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
+// Calculated / Fixed Variables
+cyl_rad = cyl_diam/2;
+tol = mesh_size/4; // Used for bounding box selection tolerance
+
+//------------------------------------------------------------------------------
+// Geometry Definition
+v1 = newv;
+Cylinder(v1) = {0.0,0.0,0.0, // center location of first face
+                0.0,cyl_height,0.0, // vector defining direction
+                cyl_rad,2*Pi};
+
+//------------------------------------------------------------------------------
+// Mesh Sizing
+MeshSize{ PointsOf{ Volume{:}; } } = mesh_size;
+
+//------------------------------------------------------------------------------
+// Physical Volumes and Surfaces
+Physical Volume("cyl-vol") = {Volume{:}};
+
+// Physical surface for mechanical BC for disp_y
+Physical Surface("cyl-surf-vis") = {1};
+Physical Surface("bc-top-disp") = {2};
+Physical Surface("bc-base-disp") = {3};
+
+//------------------------------------------------------------------------------
+// Global meshing
+Mesh.Algorithm = 6;
+Mesh.Algorithm3D = 10;
+
+General.NumThreads = num_threads;
+Mesh.MaxNumThreads1D = num_threads;
+Mesh.MaxNumThreads2D = num_threads;
+Mesh.MaxNumThreads3D = num_threads;
+
+Mesh.ElementOrder = elem_order;
+Mesh 3;
+
+//------------------------------------------------------------------------------
+// Save and exit
+Save Str(file_name);
+Exit;

pyvale/simcases/case21.i

@@ -0,0 +1,155 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,mechanical,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+endTime = 10
+timeStep = 1
+
+# Mechanical Loads/BCs
+topDispRate = ${fparse -1.5e-3 / endTime}  # m/s
+
+# Mechanical Props: SS316L @ 20degC
+ss316LEMod = 200e9       # Pa
+ss316LPRatio = 0.3      # -
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[GlobalParams]
+    displacements = 'disp_x disp_y disp_z'
+[]
+
+[Mesh]
+    type = FileMesh
+    file = 'case21.msh'
+[]
+
+[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 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'
+    []
+[]
+
+[BCs]
+    [bottom_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+    [bottom_y]
+        type = DirichletBC
+        variable = disp_y
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+    [bottom_z]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+
+
+    [top_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-top-disp'
+        value = 0.0
+    []
+    [top_y]
+        type = FunctionDirichletBC
+        variable = disp_y
+        boundary = 'bc-top-disp'
+        function = '${topDispRate}*t'
+    []
+    [top_z]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-top-disp'
+        value = 0.0
+    []
+[]
+
+[Materials]
+    [elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${ss316LEMod}
+        poissons_ratio = ${ss316LPRatio}
+    []
+    [stress]
+        type = ComputeFiniteStrainElasticStress
+    []
+[]
+
+[Preconditioning]
+    [SMP]
+        type = SMP
+        full = true
+    []
+[]
+
+[Executioner]
+    type = Transient
+
+    solve_type = 'NEWTON'
+    petsc_options = '-snes_converged_reason'
+    petsc_options_iname = '-pc_type -pc_hypre_type'
+    petsc_options_value = 'hypre boomeramg'
+
+    l_max_its = 500
+    l_tol = 1e-6
+
+    nl_max_its = 100
+    nl_rel_tol = 1e-6
+    nl_abs_tol = 1e-6
+
+    end_time= ${endTime}
+    dt = ${timeStep}
+
+    [Predictor]
+        type = SimplePredictor
+        scale = 1
+    []
+[]
+
+
+[Postprocessors]
+    [react_y_bot]
+        type = SidesetReaction
+        direction = '0 1 0'
+        stress_tensor = stress
+        boundary = 'bc-base-disp'
+    []
+    [react_y_top]
+        type = SidesetReaction
+        direction = '0 1 0'
+        stress_tensor = stress
+        boundary = 'bc-top-disp'
+    []
+
+    [disp_y_max]
+        type = NodalExtremeValue
+        variable = disp_y
+    []
+    [disp_x_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+    [disp_z_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]

pyvale/simcases/case22.geo

@@ -0,0 +1,82 @@
+//==============================================================================
+// Gmsh 3D plate imaging test case
+// author: Lloyd Fletcher (scepticalrabbit)
+//==============================================================================
+// Always set to OpenCASCADE - circles and boolean opts are much easier!
+SetFactory("OpenCASCADE");
+
+// Allows gmsh to print to terminal in vscode - easier debugging
+General.Terminal = 1;
+
+// View options - not required when
+Geometry.PointLabels = 1;
+Geometry.CurveLabels = 1;
+Geometry.SurfaceLabels = 1;
+Geometry.VolumeLabels = 0;
+
+//-------------------------------------------------------------------------
+//_* MOOSEHERDER VARIABLES - START
+file_name = "case22.msh";
+
+// Geometric variables
+plate_thick = 2e-3;
+plate_height = 40e-3;
+plate_width = 30e-3;
+
+plate_thick_layers = 2;
+
+// Must be an integer
+elem_order = 1;
+mesh_ref = 1;
+mesh_size = 1e-3/mesh_ref;
+num_threads = 4;
+//** MOOSEHERDER VARIABLES - END
+//------------------------------------------------------------------------------
+
+//------------------------------------------------------------------------------
+// Calculated / Fixed Variables
+tol = mesh_size/4; // Used for bounding box selection tolerance
+
+//------------------------------------------------------------------------------
+// Geometry Definition
+s1 = news;
+Rectangle(s1) =
+    {-plate_width/2,0.0,0.0,
+      plate_width,plate_height};
+
+//------------------------------------------------------------------------------
+// Mesh Sizing
+MeshSize{ PointsOf{ Surface{:}; } } = mesh_size;
+Transfinite Surface{Surface{:}};
+Recombine Surface{Surface{:}};
+
+Extrude{0.0,0.0,plate_thick}{
+    Surface{:}; Layers{plate_thick_layers}; Recombine;
+}
+
+//------------------------------------------------------------------------------
+// Physical Volumes and Surfaces
+Physical Volume("plate-vol") = {Volume{:}};
+
+Physical Surface("plate-surf-vis-front") = {6};
+Physical Surface("plate-surf-vis-back") = {1};
+Physical Surface("bc-top-disp") = {4};
+Physical Surface("bc-base-disp") = {2};
+
+//------------------------------------------------------------------------------
+// Global meshing
+Mesh.Algorithm = 6;
+Mesh.Algorithm3D = 10;
+
+General.NumThreads = num_threads;
+Mesh.MaxNumThreads1D = num_threads;
+Mesh.MaxNumThreads2D = num_threads;
+Mesh.MaxNumThreads3D = num_threads;
+
+Mesh.ElementOrder = elem_order;
+Mesh 3;
+
+//------------------------------------------------------------------------------
+// Save and exit
+Save Str(file_name);
+Exit;

pyvale/simcases/case22.i

@@ -0,0 +1,140 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,mechanical,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+endTime = 10
+timeStep = 1
+
+# Mechanical Loads/BCs
+topDispRate = ${fparse 1e-3 / endTime}  # m/s
+
+# Mechanical Props: SS316L @ 20degC
+ss316LEMod = 200e9       # Pa
+ss316LPRatio = 0.3      # -
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[GlobalParams]
+    displacements = 'disp_x disp_y disp_z'
+[]
+
+[Mesh]
+    type = FileMesh
+    file = 'case22.msh'
+[]
+
+[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 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'
+    []
+[]
+
+[BCs]
+    [bottom_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+    [bottom_y]
+        type = DirichletBC
+        variable = disp_y
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+    [bottom_z]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-base-disp'
+        value = 0.0
+    []
+
+
+    [top_x]
+        type = DirichletBC
+        variable = disp_x
+        boundary = 'bc-top-disp'
+        value = 0.0
+    []
+    [top_y]
+        type = FunctionDirichletBC
+        variable = disp_y
+        boundary = 'bc-top-disp'
+        function = '${topDispRate}*t'
+    []
+    [top_z]
+        type = DirichletBC
+        variable = disp_z
+        boundary = 'bc-top-disp'
+        value = 0.0
+    []
+[]
+
+[Materials]
+    [elasticity]
+        type = ComputeIsotropicElasticityTensor
+        youngs_modulus = ${ss316LEMod}
+        poissons_ratio = ${ss316LPRatio}
+    []
+    [stress]
+        type = ComputeFiniteStrainElasticStress
+    []
+[]
+
+[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 = 'bc-base-disp'
+    []
+    [react_y_top]
+        type = SidesetReaction
+        direction = '0 1 0'
+        stress_tensor = stress
+        boundary = 'bc-top-disp'
+    []
+
+    [disp_y_max]
+        type = NodalExtremeValue
+        variable = disp_y
+    []
+    [disp_x_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+    [disp_z_max]
+        type = NodalExtremeValue
+        variable = disp_x
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]