pyvale 2025.5.3__cp311-cp311-macosx_13_0_x86_64.whl

pyvale/simcases/case10.i

@@ -0,0 +1,257 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,3Dstcgmsh,1mat,thermomechanical,steady,
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+# Thermal Loads/BCs
+coolantTemp = 100.0      # degC
+heatTransCoeff = 125.0e3 # W.m^-2.K^-1
+surfHeatFlux = 5.0e6    # W.m^-2
+
+# Material Properties:
+# Thermal Props:OFHC) Copper at 250degC
+ss316LDensity = 8829.0  # kg.m^-3
+ss316LThermCond = 384.0 # W.m^-1.K^-1
+ss316LSpecHeat = 406.0  # J.kg^-1.K^-1
+
+# Mechanical Props: OFHC Copper 250degC
+ss316LEMod = 108e9       # Pa
+ss316LPRatio = 0.33      # -
+
+# Thermo-mechanical coupling
+stressFreeTemp = 150 # degC
+cuThermExp = 17.8e-6 # 1/degC
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+
+[GlobalParams]
+    displacements = 'disp_x disp_y disp_z'
+[]
+
+[Mesh]
+    type = FileMesh
+    file = 'case10.msh'
+[]
+
+[Variables]
+    [temperature]
+        family = LAGRANGE
+        order = FIRST
+        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 = ${cuThermExp}
+        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 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'
+    #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
+    []
+
+    #[strain_xx_avg]
+    #    type = ElementAverageValue
+    #    variable = strain_xx
+    #[]
+    #[strain_yy_avg]
+    #    type = ElementAverageValue
+    #    variable = strain_yy
+    #[]
+[]
+
+[Outputs]
+    exodus = true
+[]

pyvale/simcases/case11.geo

@@ -0,0 +1,337 @@
+//==============================================================================
+// Gmsh 3D Monoblock
+// author: Lloyd Fletcher (scepticalrabbit)
+//==============================================================================
+SetFactory("OpenCASCADE");
+General.Terminal = 0;
+
+// View options: turns on/off which pieces of geometry we show
+Geometry.Points = 0;
+Geometry.Curves = 0;
+Geometry.Surfaces = 0;
+Geometry.Volumes = 0;
+
+// View options: turns labels on and off for different types of geometry
+Geometry.PointLabels = 0;
+Geometry.CurveLabels = 0;
+Geometry.SurfaceLabels = 0;
+Geometry.VolumeLabels = 0;
+
+//------------------------------------------------------------------------------
+// Variable Definitions
+
+//------------------------------------------------------------------------------
+//_* MOOSEHERDER VARIABLES - START
+file_name = "case11.msh";
+num_threads = 7;
+
+// Specified Geometry variables
+pipe_rad_int = 6e-3;
+pipe_thick = 1.5e-3;
+
+interlayer_thick = 2e-3;
+
+monoblock_depth = 12e-3;
+monoblock_side = 3e-3;
+monoblock_arm_height = 8e-3;
+
+// Specified Mesh variables
+base_divs = 1;
+mesh_ref = 1; //  Must be an integer greater than 0
+
+//_* MOOSEHERDER VARIABLES - END
+//------------------------------------------------------------------------------
+
+// Calculated Geometry Variables
+pipe_rad_ext = pipe_rad_int+pipe_thick;
+
+interlayer_rad_int = pipe_rad_ext;
+interlayer_rad_ext = interlayer_rad_int+interlayer_thick;
+
+monoblock_width = 2*interlayer_rad_ext + 2*monoblock_side;
+monoblock_height = monoblock_width + monoblock_arm_height;
+
+pipe_cent_x = 0.0;
+pipe_cent_y = interlayer_rad_ext + monoblock_side;
+
+// Calculated Mesh Variables
+pipe_sect_nodes = Round(mesh_ref*5); // Must be odd
+pipe_rad_nodes = Round(mesh_ref*5);
+interlayer_rad_nodes = Round(mesh_ref*5);
+monoblock_side_nodes = Round(mesh_ref*5);
+monoblock_arm_nodes = Round(mesh_ref*5);
+monoblock_depth_nodes = Round(mesh_ref*2);
+monoblock_width_nodes = Floor((pipe_sect_nodes-1)/2)+1;
+/*
+// This is a more reasonable mesh refinement for the monoblock but solve time
+// is much longer
+pipe_sect_nodes = Round(mesh_ref*11); // Must be odd
+pipe_rad_nodes = Round(mesh_ref*7);
+interlayer_rad_nodes = Round(mesh_ref*7);
+monoblock_side_nodes = Round(mesh_ref*9);
+monoblock_arm_nodes = Round(mesh_ref*11);
+monoblock_depth_nodes = Round(mesh_ref*5);
+monoblock_width_nodes = Floor((pipe_sect_nodes-1)/2)+1;
+*/
+
+// Calculate approx element size by dividing the circumference
+elem_size = 2*Pi*pipe_rad_int/(4*(pipe_sect_nodes-1));
+tol = elem_size/4; // Used for selection tolerance of bounding boxes
+
+//------------------------------------------------------------------------------
+// Geometry Definition
+s1 = news;
+Rectangle(s1) =
+    {-monoblock_width/2,0.0,0.0,
+    monoblock_width/2,monoblock_width/2};
+
+s2 = news;
+Rectangle(s2) =
+    {-monoblock_width/2,monoblock_width/2,0.0,
+    monoblock_width/2,monoblock_width/2};
+
+s3 = news;
+Rectangle(s3) =
+    {0.0,monoblock_width/2,0.0,
+    monoblock_width/2,monoblock_width/2};
+
+s4 = news;
+Rectangle(s4) =
+    {0.0,0.0,0.0,
+    monoblock_width/2,monoblock_width/2};
+
+sa1 = news;
+Rectangle(sa1) =
+    {-monoblock_width/2,monoblock_width,0.0,
+    monoblock_width/2,monoblock_arm_height};
+
+sa2 = news;
+Rectangle(sa2) =
+    {0.0,monoblock_width,0.0,
+    monoblock_width/2,monoblock_arm_height};
+
+BooleanFragments{ Surface{s1}; Delete; }{ Surface{s2,s3,s4,sa1,sa2}; Delete; }
+
+pipe_block_surfs() = Surface In BoundingBox{
+        -monoblock_width/2-tol,-tol,-tol,
+        monoblock_width+tol,monoblock_width+tol,tol};
+
+c2 = newc; Circle(c2) = {pipe_cent_x,pipe_cent_y,0.0,interlayer_rad_ext};
+cl2 = newcl; Curve Loop(cl2) = {c2};
+s6 = news; Plane Surface(s6) = {cl2};
+
+c3 = newc; Circle(c3) = {pipe_cent_x,pipe_cent_y,0.0,pipe_rad_ext};
+cl3 = newcl; Curve Loop(cl3) = {c3};
+s7 = news; Plane Surface(s7) = {cl3};
+
+c4 = newc; Circle(c4) = {pipe_cent_x,pipe_cent_y,0.0,pipe_rad_int};
+cl4 = newcl; Curve Loop(cl4) = {c4};
+s8 = news; Plane Surface(s8) = {cl4};
+
+BooleanDifference{Surface{pipe_block_surfs(),s6,s7}; Delete;}
+                 { Surface{s8}; Delete;}
+
+//------------------------------------------------------------------------------
+// Transfinite Line/Curve Meshing
+
+cm1() = Curve In BoundingBox{
+    pipe_cent_x-interlayer_rad_ext-tol,pipe_cent_y-interlayer_rad_ext-tol,-tol,
+    pipe_cent_x+interlayer_rad_ext+tol,pipe_cent_y+interlayer_rad_ext+tol,+tol};
+
+For ss In {0:#cm1()-1}
+    Transfinite Curve(cm1(ss)) = pipe_sect_nodes;
+EndFor
+
+cm1() = Curve In BoundingBox{
+    pipe_cent_x+pipe_rad_int-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x+pipe_rad_ext+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm1(0)) = pipe_rad_nodes;
+
+cm2() = Curve In BoundingBox{
+    pipe_cent_x-pipe_rad_ext-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x-pipe_rad_int+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm2(0)) = pipe_rad_nodes;
+
+cm3() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y-pipe_rad_ext-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y-pipe_rad_int+tol,+tol};
+Transfinite Curve(cm3(0)) = pipe_rad_nodes;
+
+cm4() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y+pipe_rad_int-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y+pipe_rad_ext+tol,+tol};
+Transfinite Curve(cm4(0)) = pipe_rad_nodes;
+
+
+cm1() = Curve In BoundingBox{
+    pipe_cent_x+interlayer_rad_int-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x+interlayer_rad_ext+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm1(0)) = interlayer_rad_nodes;
+
+cm2() = Curve In BoundingBox{
+    pipe_cent_x-interlayer_rad_ext-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x-interlayer_rad_int+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm2(0)) = interlayer_rad_nodes;
+
+cm3() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y-interlayer_rad_ext-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y-interlayer_rad_int+tol,+tol};
+Transfinite Curve(cm3(0)) = interlayer_rad_nodes;
+
+cm4() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y+interlayer_rad_int-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y+interlayer_rad_ext+tol,+tol};
+Transfinite Curve(cm4(0)) = interlayer_rad_nodes;
+
+
+cm1() = Curve In BoundingBox{
+    pipe_cent_x+interlayer_rad_ext-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x+monoblock_width/2+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm1(0)) = monoblock_side_nodes;
+
+cm2() = Curve In BoundingBox{
+    pipe_cent_x-monoblock_width/2-tol,pipe_cent_y-tol,-tol,
+    pipe_cent_x-interlayer_rad_ext+tol,pipe_cent_y+tol,+tol};
+Transfinite Curve(cm2(0)) = monoblock_side_nodes;
+
+cm3() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y-monoblock_width/2-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y-interlayer_rad_ext+tol,+tol};
+Transfinite Curve(cm3(0)) = monoblock_side_nodes;
+
+cm4() = Curve In BoundingBox{
+    pipe_cent_x-tol,pipe_cent_y+interlayer_rad_ext-tol,-tol,
+    pipe_cent_x+tol,pipe_cent_y+monoblock_width/2+tol,+tol};
+Transfinite Curve(cm4(0)) = monoblock_side_nodes;
+
+
+cm1() = Curve In BoundingBox{
+    -monoblock_width/2-tol,monoblock_width-tol,-tol,
+    monoblock_width/2+tol,monoblock_width+monoblock_arm_height+tol,+tol};
+
+For ss In {0:#cm1()-1}
+
+    Transfinite Curve(cm1(ss)) = monoblock_arm_nodes;
+EndFor
+
+
+
+cm1() = Curve In BoundingBox{
+    -monoblock_width/2-tol,-tol,-tol,
+    -monoblock_width/2+tol,monoblock_width+tol,+tol};
+
+For ss In {0:#cm1()-1}
+    Transfinite Curve(cm1(ss)) = monoblock_width_nodes;
+EndFor
+
+
+cm2() = Curve In BoundingBox{
+    monoblock_width/2-tol,-tol,-tol,
+    monoblock_width/2+tol,monoblock_width+tol,+tol};
+
+For ss In {0:#cm2()-1}
+    Transfinite Curve(cm2(ss)) = monoblock_width_nodes;
+EndFor
+
+
+// Mesh top, bottom and armour horizontal lines
+cm3() = Curve In BoundingBox{
+    -monoblock_width/2-tol,-tol,-tol,
+    monoblock_width/2+tol,+tol,+tol};
+
+For ss In {0:#cm3()-1}
+    Transfinite Curve(cm3(ss)) = monoblock_width_nodes;
+EndFor
+
+
+cm3() = Curve In BoundingBox{
+    -monoblock_width/2-tol,monoblock_width-tol,-tol,
+    monoblock_width/2+tol,monoblock_width+tol,+tol};
+
+For ss In {0:#cm3()-1}
+    Transfinite Curve(cm3(ss)) = monoblock_width_nodes;
+EndFor
+
+
+cm4() = Curve In BoundingBox{
+    -monoblock_width/2-tol,monoblock_width+monoblock_arm_height+-tol,-tol,
+    monoblock_width/2+tol,monoblock_width+monoblock_arm_height+tol,+tol};
+
+For ss In {0:#cm4()-1}
+    Transfinite Curve(cm4(ss)) = monoblock_width_nodes;
+EndFor
+
+
+//------------------------------------------------------------------------------
+// Transfinite Surface Meshing
+
+// Mesh the pipe and interlayer
+sm1() = Surface In BoundingBox{
+    pipe_cent_x-interlayer_rad_ext-tol,pipe_cent_y-interlayer_rad_ext-tol,-tol,
+    pipe_cent_x+interlayer_rad_ext+tol,pipe_cent_y+interlayer_rad_ext+tol,+tol};
+
+For ss In {0:#sm1()-1}
+    Transfinite Surface{sm1(ss)};
+    Recombine Surface{sm1(ss)};
+EndFor
+
+// Mesh the armour on top of the block
+sm1() = Surface In BoundingBox{
+    -monoblock_width/2-tol,monoblock_width-tol,-tol,
+    monoblock_width/2+tol,monoblock_width+monoblock_arm_height+tol,+tol};
+
+For ss In {0:#sm1()-1}
+    Transfinite Surface{sm1(ss)};
+    Recombine Surface{sm1(ss)};
+EndFor
+
+// Mesh the block around the interlayer
+Transfinite Surface{22} = {15,16,17,13};
+Recombine Surface{22};
+
+Transfinite Surface{27} = {15,24,5,13};
+Recombine Surface{27};
+
+Transfinite Surface{30} = {24,5,28,27};
+Recombine Surface{30};
+
+Transfinite Surface{31} = {27,28,17,16};
+Recombine Surface{31};
+
+//------------------------------------------------------------------------------
+// Extrude the surface mesh to 3D
+Extrude{0.0,0.0,monoblock_depth}{
+    Surface{:}; Layers{monoblock_depth_nodes}; Recombine;
+}
+
+//------------------------------------------------------------------------------
+// Physical Surfaces for Loads and Boundary Condition
+Physical Surface("bc-top-heatflux") = {36,41};
+Physical Surface("bc-pipe-heattransf") = {54,59,83,67};
+
+//------------------------------------------------------------------------------
+// Physical Volumes for Material Defs
+Physical Volume("pipe-cucrzr") = {6,9,5,14};
+Physical Volume("interlayer-cu") = {4,7,10,13};
+Physical Volume("armour-w") = {1,2,8,11,3,12};
+
+//------------------------------------------------------------------------------
+// Global Mesh controls
+Mesh.Algorithm = 6;
+Mesh.Algorithm3D = 10;
+
+General.NumThreads = num_threads;
+Mesh.MaxNumThreads1D = num_threads;
+Mesh.MaxNumThreads2D = num_threads;
+Mesh.MaxNumThreads3D = num_threads;
+
+Mesh.ElementOrder = 2;
+
+Mesh 3;
+
+//------------------------------------------------------------------------------
+// Save and exit
+Save Str(file_name);
+Exit;
+
+