Srlife moose interface

.github/workflows/formatting.yml

@@ -1,5 +1,5 @@
 name: Check code formatting
-on: [push, pull_request]
+on: [push, pull_request, workflow_dispatch]
 jobs:
   check-formatting:
     runs-on: ubuntu-latest

.github/workflows/test.yml

@@ -1,5 +1,5 @@
 name: tests
-on: [push, pull_request]
+on: [push, pull_request, workflow_dispatch]
 jobs:
   test-linux:
     runs-on: ubuntu-latest

.gitmodules

@@ -0,0 +1,6 @@
+[submodule "moose"]
+	path = moose
+	url = https://github.com/idaholab/moose.git
+[submodule "nemlapp"]
+	path = nemlapp
+	url = https://github.com/sagarbhatt0904/nemlapp.git

install.sh

@@ -0,0 +1,81 @@
+#!/bin/bash
+set -e
+
+# Parse arguments
+BUILD_MOOSE=false
+if [[ "$1" == "--with-moose" ]]; then
+    BUILD_MOOSE=true
+    echo "Installing srlife with MOOSE..."
+    # Initialize MOOSE submodule
+    git submodule update --init moose
+    cd moose
+    git checkout master
+    cd ../
+else
+    echo "Installing srlife only..."
+fi
+
+# Get repo root
+SRLIFE_DIR="$(pwd)"
+ENV_NAME="srlifeMoose"
+MOOSE_JOBS=32
+
+# Create conda environment
+echo "Creating conda environment..."
+# ensure conda shell functions are available in this script
+eval "$(conda shell.bash hook)"
+
+# Create the environment only if it doesn't already exist
+if conda env list | awk '!/^#/ && NF>0 {print $1}' | grep -Fxq "$ENV_NAME"; then
+    echo "Conda environment '$ENV_NAME' already exists; skipping creation."
+else
+    conda create -n "$ENV_NAME" moose-dev=2025.09.18=mpich
+    #TODO:install seacas --- maybe add this to requirements.txt later
+    conda install moose-seacas
+fi
+
+
+# Activate environment
+conda activate $ENV_NAME
+
+# Install Python dependencies
+echo "Installing Python dependencies..."
+pip3 install --user wheel
+pip3 install -r requirements.txt
+
+
+
+
+# Build MOOSE (if specified)
+if [[ "$BUILD_MOOSE" == "true" ]]; then
+    echo "Building MOOSE..."
+    cd $SRLIFE_DIR/moose/test
+    make -j$MOOSE_JOBS
+
+    # Set environment variables
+    echo "Installing Thermohydraylics module..."
+    cd $SRLIFE_DIR/moose/modules/thermal_hydraulics/
+    make -j$MOOSE_JOBS
+
+    echo "Installing nemlapp for solid mechanics..."
+    cd $SRLIFE_DIR/nemlapp
+    make -j$MOOSE_JOBS
+fi
+
+echo "Installation complete!"
+echo ""
+if [[ "$BUILD_MOOSE" == "true" ]]; then
+    echo "To use add the following to your rc file:"
+    echo "  conda activate $ENV_NAME"
+    echo "  export PYTHONPATH=$SRLIFE_DIR/"
+    echo "  export NEML_DIR=<PATH_TO_NEML>/neml"
+    echo "  export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:<PATH_TO_NEML>/neml/lib"
+    echo "  export MOOSE_DIR=$SRLIFE_DIR/moose/"
+    echo "  export MOOSE_THM=$SRLIFE_DIR/moose/modules/thermal_hydraulics/thermal_hydraulics-opt"
+    echo "  export NEMLAPP=$SRLIFE_DIR/nemlapp/nemlapp-opt"
+    echo "  export MOOSE_MPI=$(which mpirun)"
+    echo "  export MOOSE_NPROCS=<SET NUM OF MOOSE JOBS>"
+else
+    echo "To use:"
+    echo "  conda activate $ENV_NAME"
+fi

requirements.txt

@@ -1,5 +1,5 @@
 numpy
-scipy==1.9.2
+scipy>=1.9.2
 h5py
 vtk
 pytest

srlife/damage.py

@@ -1033,7 +1033,7 @@ def calculate_volume_flaw_flattened_eq_stress(
         else:
             sigma_e[sigma_e < 0] = 0
             g = (
-                np.trapz(
+                np.trapezoid(
                     (sigma_e / (sigma_e_max + self.tolerance)) ** Navg[..., None, None],
                     time,
                     axis=0,
@@ -1196,7 +1196,7 @@ def calculate_surface_flaw_flattened_eq_stress(
         else:
             sigma_e[sigma_e < 0] = 0
             g = (
-                np.trapz(
+                np.trapezoid(
                     (sigma_e / (sigma_e_max + self.tolerance))
                     ** Navg[..., None, None, None],
                     time,
@@ -1474,7 +1474,7 @@ def calculate_volume_flaw_element_log_reliability(
         # Calculating ratio of cyclic stress to max cyclic stress (in one cycle)
         # for time-independent and time-dependent cases
         g = (
-            np.trapz(
+            np.trapezoid(
                 (pstress / (pstress_max + 1.0e-14)) ** Navg[..., None],
                 time,
                 axis=0,
@@ -1556,7 +1556,7 @@ def calculate_surface_time_dep_pstress(
         # Calculating ratio of cyclic stress to max cyclic stress (in one cycle)
         # for time-independent and time-dependent cases
         g = (
-            np.trapz(
+            np.trapezoid(
                 (surf_pstress / (surf_pstress_max + 1.0e-14)) ** Navg[..., None, None],
                 time,
                 axis=0,
@@ -1743,7 +1743,7 @@ def calculate_volume_flaw_avg_normal_stress(
             sigma_n_0 = sigma_n
         else:
             g = (
-                np.trapz(
+                np.trapezoid(
                     (sigma_n / (sigma_n_max + self.tolerance)) ** Navg[..., None, None],
                     time,
                     axis=0,
@@ -1890,7 +1890,7 @@ def calculate_surface_flaw_time_dep_normal_stress(
             sigma_n_0 = sigma_n
         else:
             g = (
-                np.trapz(
+                np.trapezoid(
                     (sigma_n / (sigma_n_max + self.tolerance))
                     ** Navg[..., None, None, None],
                     time,

srlife/data/damage/SiC.xml

@@ -46,12 +46,12 @@
   </base>
   <cares type="StandardModel">
     <threshold_vol>
-      <temperatures>298.15 1073.15 1273.15 1473.15 1673.15 1773.15</temperatures>
-      <values>0.0 0.0 0.0 0.0 0.0 0.0</values>
+      <temperatures>298.15 811.15 1089.15 1366.15 1566.15</temperatures>
+      <values>0.0 0.0 0.0 0.0 0.0</values>
     </threshold_vol>
     <threshold_surf>
-      <temperatures>298.15 1073.15 1273.15 1473.15 1673.15 1773.15</temperatures>
-      <values>0.0 0.0 0.0 0.0 0.0 0.0</values>
+      <temperatures>298.15 811.15 1089.15 1366.15 1566.15</temperatures>
+      <values>0.0 0.0 0.0 0.0 0.0</values>
     </threshold_surf>
     <strength_vol>
       <temperatures>298.15 811.15 1089.15 1366.15 1566.15</temperatures>
@@ -88,4 +88,100 @@
       <values>520 520 750 210 210</values>
     </fatigue_Bs>
   </cares>
+  <tested_2P type="StandardModel">
+    <threshold_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>0.0 0.0 0.0 0.0</values>
+    </threshold_vol>
+    <threshold_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>0.0 0.0 0.0 0.0</values>
+    </threshold_surf>
+    <strength_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>515.3 515.3 515.3 515.3</values>
+    </strength_vol>
+    <strength_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>711.8 711.8 711.8 711.8</values>
+    </strength_surf>
+    <modulus_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>8.29 8.29 8.29 8.29</values>
+    </modulus_vol>
+    <modulus_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>6.75 6.75 6.75 6.75</values>
+    </modulus_surf>
+    <c_bar>0.82</c_bar>
+    <nu>0.16</nu>
+    <fatigue_Nv>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>63.65 63.65 63.65 63.65</values>
+      <!--values>60.93 60.93 60.93 60.93</values-->
+    </fatigue_Nv>
+    <fatigue_Ns>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>63.65 63.65 63.65 63.65</values>
+      <!--values>60.93 60.93 60.93 60.93</values-->
+    </fatigue_Ns>
+    <fatigue_Bv>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>14.29 2.84 2.84 2.84</values>
+      <!--values>11.38 1.21 8.93 8.93</values-->
+    </fatigue_Bv>
+    <fatigue_Bs>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>14.29 2.84 2.84 2.84</values>
+      <!--values>11.38 1.21 8.93 8.93</values-->
+    </fatigue_Bs>
+  </tested_2P>
+  <tested_3P type="StandardModel">
+    <threshold_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>252.2 252.2 252.2 252.2</values>
+    </threshold_vol>
+    <threshold_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>191.9 191.9 191.9 191.9</values>
+    </threshold_surf>
+    <strength_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>1690.2 1690.2 1690.2 1690.2</values>
+    </strength_vol>
+    <strength_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>1291.6 1291.6 1291.6 1291.6</values>
+    </strength_surf>
+    <modulus_vol>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>1.11 1.11 1.11 1.11</values>
+    </modulus_vol>
+    <modulus_surf>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>2.02 2.02 2.02 2.02</values>
+    </modulus_surf>
+    <c_bar>0.82</c_bar>
+    <nu>0.16</nu>
+    <fatigue_Nv>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>63.65 63.65 63.65 63.65</values>
+      <!--values>60.93 60.93 60.93 60.93</values-->
+    </fatigue_Nv>
+    <fatigue_Ns>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>63.65 63.65 63.65 63.65</values>
+      <!--values>60.93 60.93 60.93 60.93</values-->
+    </fatigue_Ns>
+    <fatigue_Bv>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>14.29 2.84 2.84 2.84</values>
+      <!--values>11.38 1.21 8.93 8.93</values-->
+    </fatigue_Bv>
+    <fatigue_Bs>
+      <temperatures>298.15 873.15 1073.15 1273.15</temperatures>
+      <values>14.29 2.84 2.84 2.84</values>
+      <!--values>11.38 1.21 8.93 8.93</values-->
+    </fatigue_Bs>
+  </tested_3P>
 </models>