pyvale 2025.5.3__cp311-cp311-manylinux_2_17_x86_64.manylinux2014_x86_64.whl

pyvale/simcases/case15.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 = "case15.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;
+
+

pyvale/simcases/case15.i

@@ -0,0 +1,150 @@
+#-------------------------------------------------------------------------
+# pyvale: gmsh,monoblock,3mat,thermal,transient
+#-------------------------------------------------------------------------
+
+#-------------------------------------------------------------------------
+#_* MOOSEHERDER VARIABLES - START
+
+endTime = 25
+timeStep = 1
+
+coolantTemp=100.0      # degC
+heatTransCoeff=125.0e3 # W.m^-2.K^-1
+surfHeatFlux=10.0e6    # W.m^-2
+timeConst = 1   # s
+
+# 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
+
+# Pure Copper pg.134 at 250degC
+cuDensity = 8829.0  # kg.m^-3
+cuThermCond = 384.0 # W.m^-1.K^-1
+cuSpecHeat = 406.0  # J.kg^-1.K^-1
+
+# Tungsten pg.224 at 600degC
+wDensity = 19150.0  # kg.m^-3
+wThermCond = 127.0 # W.m^-1.K^-1
+wSpecHeat = 147.0  # J.kg^-1.K^-1
+
+# Mesh file string
+mesh_file = 'case15.msh'
+
+#** MOOSEHERDER VARIABLES - END
+#-------------------------------------------------------------------------
+
+[Mesh]
+  type = FileMesh
+  file = ${mesh_file}
+[]
+
+[Variables]
+  [temperature]
+    family = LAGRANGE
+    order = SECOND
+    initial_condition = ${coolantTemp}
+  []
+[]
+
+[Kernels]
+  [heat_conduction]
+    type = HeatConduction
+    variable = temperature
+  []
+  [time_derivative]
+    type = HeatConductionTimeDerivative
+    variable = temperature
+  []
+[]
+
+[Materials]
+  [cucrzr_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${cucrzrThermCond}
+    specific_heat = ${cucrzrSpecHeat}
+    block = 'pipe-cucrzr'
+  []
+
+  [copper_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${cuThermCond}
+    specific_heat = ${cuSpecHeat}
+    block = 'interlayer-cu'
+  []
+
+  [tungsten_thermal]
+    type = HeatConductionMaterial
+    thermal_conductivity = ${wThermCond}
+    specific_heat = ${wSpecHeat}
+    block = 'armour-w'
+  []
+
+  [cucrzr_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${cucrzrDensity}
+    block = 'pipe-cucrzr'
+  []
+
+  [copper_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${cuDensity}
+    block = 'interlayer-cu'
+  []
+
+  [tungsten_density]
+    type = GenericConstantMaterial
+    prop_names = 'density'
+    prop_values = ${wDensity}
+    block = 'armour-w'
+  []
+[]
+
+[BCs]
+  [heat_flux_in]
+    type = FunctionNeumannBC
+    variable = temperature
+    boundary = 'bc-top-heatflux'
+    function = '${fparse surfHeatFlux}*(1-exp(-(1/${timeConst})*t))'
+  []
+  [heat_flux_out]
+    type = ConvectiveHeatFluxBC
+    variable = temperature
+    boundary = 'bc-pipe-heattransf'
+    T_infinity = ${coolantTemp}
+    heat_transfer_coefficient = ${heatTransCoeff}
+  []
+[]
+
+[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]
+  [temp_max]
+      type = NodalExtremeValue
+      variable = temperature
+  []
+  [temp_avg]
+      type = AverageNodalVariableValue
+      variable = temperature
+  []
+[]
+
+[Outputs]
+  exodus = true
+[]