Merge pull request #20 from LIBRA-project/simplify-neutronics

removed unneeded stuff

Author: RemDelaporteMathurin

.github/workflows/ci.yml

@@ -25,7 +25,9 @@ jobs:
       - name: Run neutronics
         shell: bash -l {0}
         working-directory: analysis/neutron
-        run: python openmc_model.py
+        run: |
+          python openmc_model.py
+          jupyter-nbconvert --to notebook postprocessing.ipynb --execute
 
       - name: Run tritium model
         shell: bash -l {0}

analysis/neutron/.gitkeep

@@ -26,34 +26,3 @@ def get_exp_cllif_density(temp, LiCl_frac=0.695):
     rho = a + b * C + c * temp + d * C**2 + e * C**3 + f * temp * C**2 + g * C * temp**2
 
     return rho
-
-
-def calculate_cylinder_volume(radius, height):
-    volume = math.pi * radius**2 * height
-    return volume
-
-
-def translate_surface(
-    surface: (
-        openmc.XPlane | openmc.YPlane | openmc.ZPlane | openmc.Plane | openmc.Sphere
-    ),
-    dx: float,
-    dy: float,
-    dz: float,
-) -> openmc.Surface:
-    """Translate an OpenMC surface by dx, dy, dz."""
-    if isinstance(surface, openmc.XPlane):
-        surface.x0 += dx
-    elif isinstance(surface, openmc.YPlane):
-        surface.y0 += dy
-    elif isinstance(surface, openmc.ZPlane):
-        surface.z0 += dz
-    elif isinstance(surface, openmc.Plane):
-        surface.d += surface.a * dx + surface.b * dy + surface.c * dz
-    elif isinstance(surface, openmc.Sphere):
-        surface.x0 += dx
-        surface.y0 += dy
-        surface.z0 += dz
-    else:
-        raise TypeError(f"Unsupported surface type: {type(surface)}")
-    return surface

analysis/neutron/helpers.py

@@ -1,6 +1,6 @@
 import openmc
-import vault
 from libra_toolbox.neutronics.neutron_source import A325_generator_diamond
+from libra_toolbox.neutronics import vault
 import helpers
 
 
@@ -519,4 +519,4 @@ def baby_model():
     with open(processed_data_file, "w") as f:
         json.dump(existing_data, f, indent=4)
 
-    print(f"Processed data stored in {processed_data_file}")
+    print(f"Processed data stored in {processed_data_file}")

analysis/neutron/openmc_model.py

@@ -195,7 +195,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
    "outputs": [
     {
@@ -212,7 +212,7 @@
    "source": [
     "# Visualization\n",
     "\n",
-    "from vault import Air\n",
+    "from libra_toolbox.neutronics.vault import Air\n",
     "from openmc_model import air\n",
     "\n",
     "x_c = 587\n",