Merge pull request #37 from aurora-multiphysics/hsaunders1904/heat_conduction_in_infinite_system

Add unsteady heat conduction in infinite system example/test

Author: hsaunders1904

.github/workflows/ci.yml

@@ -40,4 +40,5 @@ jobs:
       shell: bash
       run: |
         . /opt/openfoam/OpenFOAM-12/etc/bashrc || true
+        pip install -r requirements.test.txt
         ./run_tests

README.md

@@ -80,10 +80,10 @@ METHOD=dbg make
 
 Some basic tests can be run with `./run_tests`.
 
-Note that you must have installed the `fluidfoam` Python package:
+Note that you must have the required Python packages installed:
 
 ```bash
-pip install "fluidfoam>=0.2.4"
+pip install -r requirements.test.txt
 ```
 
 ## Quality

include/mesh/Foam2MooseMeshGen.h

@@ -76,8 +76,8 @@ class Foam2MooseMeshAdapter
                         MPI_Comm * comm = nullptr);
   Foam2MooseMeshAdapter() = default;
   ~Foam2MooseMeshAdapter();
-  std::int32_t npoint();
-  std::int32_t nface();
+  int32_t npoint();
+  int32_t nface();
   FoamPoint const & point(uint32_t i);
   FoamFace face(uint32_t i);
   int get_patch_id(std::string const & patch_name);

requirements.test.txt

@@ -0,0 +1,2 @@
+fluidfoam>=0.2.8
+scipy>=1.15.1

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/.gitignore

@@ -0,0 +1,6 @@
+polyMesh/
+0*
+1*
+!0/
+log.*
+processor*/

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/README.md

@@ -0,0 +1,60 @@
+# Heat Conduction in an Infinite System
+
+This simulation consists of a solid and a fluid domain.
+The domains have a shared interface at x = 0.
+The analytical solution for pure diffusion in this configuration
+was derived by Yoshida et al [1].
+
+## Setup
+
+```txt
+    -------------------+-------------------    ↑
+        Solid          |    Fluid              h
+    -------------------+-------------------    ↓
+    ←------- L -------→ ←------- L -------→
+                       |
+                     x = 0
+
+```
+
+|                  | Solid | Fluid |
+|------------------|-------|-------|
+| k $[W/m/K]$      | 1     | 4     |
+| ρ.Cp $[J/m^3/K]$ | 1     | 16    |
+| L $[m]$          | 1     | 1     |
+| h $[m]$          | 0.1   | 0.1   |
+| T $[K]$          | 1     | 0     |
+
+All boundaries (other than the interface) have symmetry boundary conditions.
+The idea here is that `L` is large enough such that the system is
+effectively infinite in the time frame of the simulation.
+
+## Analytical Solution
+
+The analytical solution to the problem is:
+
+$$
+T(x, t) =
+    \begin{cases}
+        T_h
+        - \frac{(T_h - T_c) \sqrt{k_2(\rho C_p)_2}}
+            {\sqrt{k_1(\rho C_p)_1} + \sqrt{k_2 (\rho C_p)_2}}
+        \text{erfc}(\frac{-x}{2\sqrt{k_1 / (\rho C_p)_1 t}}),
+        &x < 0 \text{ (Solid medium)}\\
+
+        T_c
+        + \frac{(T_h - T_c) \sqrt{k_1(\rho C_p)_1}}
+            {\sqrt{k_1(\rho C_p)_1}
+        + \sqrt{k_2 (\rho C_p)_2}}
+        \text{erfc}(\frac{x}{2\sqrt{k_2 / (\rho C_p)_2 t}}),
+        &x > 0 \text{ (Fluid medium)}
+    \end{cases}
+$$
+
+## References
+
+[1] Hiroaki Yoshida, Takayuki Kobayashi, Hidemitsu Hayashi, Tomoyuki Kinjo,
+Hitoshi Washizu, and Kenji Fukuzawa. Boundary condition at a two-phase
+interface in the lattice boltzmann method for the convection diffusion
+equation.
+Physical Review E, 90(1):013303, 2014. doi: https://doi.org/10.1103/PhysRevE.90.013303.

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/fluid-openfoam.i

@@ -0,0 +1,51 @@
+[Mesh]
+    type = FoamMesh
+    case = 'fluid-openfoam'
+    foam_patch = 'left'
+[]
+
+[Variables]
+    [dummy]
+        family = MONOMIAL
+        order = CONSTANT
+        initial_condition = 999
+    []
+[]
+
+[AuxVariables]
+    [fluid_wall_temp]
+        family = MONOMIAL
+        order = CONSTANT
+        initial_condition = 0.075
+    []
+    [solid_heat_flux]
+        family = MONOMIAL
+        order = CONSTANT
+        initial_condition = 0
+    []
+[]
+
+[Problem]
+    type = FoamProblem
+    # Take the heat flux from MOOSE and set it on the OpenFOAM mesh.
+    heat_flux = solid_heat_flux
+    # Take the boundary temperature from OpenFOAM and set it on the MOOSE mesh.
+    foam_temp = fluid_wall_temp
+[]
+
+[Executioner]
+    type = Transient
+    start_time = 0
+    end_time = 0.01
+    dt = 0.0003125
+
+    [TimeSteppers]
+        [foam]
+            type = FoamTimeStepper
+        []
+    []
+[]
+
+[Outputs]
+    exodus = true
+[]

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/fluid-openfoam/0/T

@@ -0,0 +1,51 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  10
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      T;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [0 0 0 1 0 0 0];
+
+internalField   uniform 1e-10;
+
+boundaryField
+{
+    left
+    {
+        type            fixedGradient;
+        gradient        uniform 0.75;
+    }
+    right
+    {
+        type            zeroGradient;
+    }
+    top
+    {
+        type            zeroGradient;
+    }
+    front
+    {
+        type            zeroGradient;
+    }
+    bottom
+    {
+        type            zeroGradient;
+    }
+    back
+    {
+        type            zeroGradient;
+    }
+}
+
+
+// ************************************************************************* //

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/fluid-openfoam/0/U

@@ -0,0 +1,31 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  12
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volVectorField;
+    location    "0";
+    object      U;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 0 1 -1 0 0 0 0 ];
+
+internalField   uniform (0 0 0);
+
+boundaryField
+{
+
+    ".*"
+    {
+        type            fixedValue;
+        value           $internalField;
+    }
+}
+
+// ************************************************************************* //

test/tests/multiapps/unsteady_heat_conduction_in_infinite_system/fluid-openfoam/0/p

@@ -0,0 +1,30 @@
+/*--------------------------------*- C++ -*----------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     | Website:  https://openfoam.org
+    \\  /    A nd           | Version:  12
+     \\/     M anipulation  |
+\*---------------------------------------------------------------------------*/
+FoamFile
+{
+    format      ascii;
+    class       volScalarField;
+    location    "0";
+    object      p;
+}
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+dimensions      [ 1 -1 -2 0 0 0 0 ];
+
+internalField   uniform 0;
+
+boundaryField
+{
+    ".*"
+    {
+        type            calculated;
+        value           $internalField;
+    }
+}
+
+// ************************************************************************* //