solidThermophysicalTransportModel: new thermophysical transport model…

… for solids to handle isotropic and anisotropic is a consistent, general and extensible manner, replacing the horrible hacks which were in solidThermo. This is entirely consistent with thermophysicalTransportModel for fluids and provides the q() and divq() for the solid energy conservation equations. The transport model and properties are specified in the optional thermophysicalTransport dictionary, the default model being isotropic if this dictionary file is not present, thus providing complete backward-compatibility for the common isotropic cases. Anisotropic thermal conductivity is now handled in a much more general manner by the anisotropic model: Class Foam::solidThermophysicalTransportModels::anisotropic Description Solid thermophysical transport model for anisotropic thermal conductivity The anisotropic thermal conductivity field is evaluated from the solid material anisotropic kappa specified in the physicalProperties dictionary transformed into the global coordinate system using default coordinate system and optionally additional coordinate systems specified per-zone in the thermophysicalProperties dictionary. Usage Example of the anisotropic thermal conductivity specification in thermophysicalProperties with two zone-based coordinate systems in addition to the default: \verbatim model anisotropic; // Default coordinate system coordinateSystem { type cartesian; origin (0 0 0); coordinateRotation { type cylindrical; e3 (1 0 0); } } // Optional zone coordinate systems zones { coil1 { type cartesian; origin (0.1 0.2 0.7); coordinateRotation { type cylindrical; e3 (0.5 0.866 0); } } coil2 { type cartesian; origin (0.4 0.5 1); coordinateRotation { type cylindrical; e3 (0.866 0.5 0); } } } \endverbatim This development required substantial rationalisation of solidThermo, coordinateSystems and updates to the solid solver module, solidDisplacementFoam, the wallHeatFlux functionObject, thermalBaffle and all coupled thermal boundary conditions.
OpenFOAM · Oct 19, 2022 · cdaaa61 · cdaaa61
1 parent cc868b6
commit cdaaa61
Show file tree

Hide file tree

Showing 80 changed files with 1,576 additions and 612 deletions.
diff --git a/applications/solvers/combustion/XiFoam/Make/options b/applications/solvers/combustion/XiFoam/Make/options
@@ -17,7 +17,7 @@ EXE_LIBS = \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
     -lthermophysicalTransportModels \
-    -lfluidThermophysicalModels \
+    -lcoupledThermophysicalTransportModels \
     -lmulticomponentThermophysicalModels \
     -lspecie \
     -llaminarFlameSpeedModels \

diff --git a/applications/solvers/combustion/XiFoam/PDRFoam/Make/options b/applications/solvers/combustion/XiFoam/PDRFoam/Make/options
@@ -41,7 +41,7 @@ EXE_LIBS = \
     -lspecie \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels \
+    -lcoupledThermophysicalTransportModels \
     -llaminarFlameSpeedModels \
     -lfiniteVolume \
     -lfvModels \

diff --git a/applications/solvers/compressible/rhoCentralFoam/Make/options b/applications/solvers/compressible/rhoCentralFoam/Make/options
@@ -18,5 +18,5 @@ EXE_LIBS = \
     -lrhoCentralFoam \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels \
+    -lcoupledThermophysicalTransportModels \
     -lmeshTools
diff --git a/applications/solvers/compressible/rhoPorousSimpleFoam/Make/options b/applications/solvers/compressible/rhoPorousSimpleFoam/Make/options
@@ -16,7 +16,7 @@ EXE_LIBS = \
     -lspecie \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels \
+    -lcoupledThermophysicalTransportModels \
     -lfiniteVolume \
     -lsampling \
     -lmeshTools \

diff --git a/applications/solvers/heatTransfer/thermoFoam/Make/options b/applications/solvers/heatTransfer/thermoFoam/Make/options
@@ -18,4 +18,4 @@ EXE_LIBS = \
     -lspecie \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels
+    -lcoupledThermophysicalTransportModels
diff --git a/applications/solvers/modules/fluid/fluid/Make/options b/applications/solvers/modules/fluid/fluid/Make/options
@@ -17,7 +17,7 @@ LIB_LIBS = \
     -lfluidThermophysicalModels \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels \
+    -lcoupledThermophysicalTransportModels \
     -lfiniteVolume \
     -lmeshTools \
     -lsampling \

diff --git a/applications/solvers/modules/fluid/multicomponentFluid/Make/options b/applications/solvers/modules/fluid/multicomponentFluid/Make/options
@@ -28,7 +28,7 @@ LIB_LIBS = \
     -lmulticomponentThermophysicalModels \
     -lmomentumTransportModels \
     -lcompressibleMomentumTransportModels \
-    -lthermophysicalTransportModels \
+    -lcoupledThermophysicalTransportModels \
     -lfluidMulticomponentThermophysicalTransportModels \
     -lfiniteVolume \
     -lfvModels \

diff --git a/applications/solvers/modules/solid/solid/Make/options b/applications/solvers/modules/solid/solid/Make/options
@@ -2,13 +2,14 @@ EXE_INC = \
     -I$(LIB_SRC)/physicalProperties/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \
-    -I$(LIB_SRC)/finiteVolume/cfdTools \
+    -I$(LIB_SRC)/ThermophysicalTransportModels/solidThermophysicalTransportModels/lnInclude \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude \
     -I$(LIB_SRC)/sampling/lnInclude
 
 LIB_LIBS = \
     -lsolidThermo \
+    -lsolidThermophysicalTransportModels \
     -lspecie \
     -lfiniteVolume \
     -lmeshTools \

diff --git a/applications/solvers/modules/solid/solid/solid.C b/applications/solvers/modules/solid/solid/solid.C
@@ -76,9 +76,6 @@ void Foam::solvers::solid::correctDiNum()
 }
 
 
-// * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
-
-
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
 Foam::solvers::solid::solid(fvMesh& mesh)
@@ -90,6 +87,8 @@ Foam::solvers::solid::solid(fvMesh& mesh)
 
     T(thermo.T()),
 
+    thermophysicalTransport(solidThermophysicalTransportModel::New(thermo)),
+
     DiNum(0)
 {
     // Read the controls
@@ -174,7 +173,7 @@ void Foam::solvers::solid::thermophysicalPredictor()
         fvScalarMatrix eEqn
         (
             fvm::ddt(rho, e)
-          + thermo.divq(e)
+          + thermophysicalTransport->divq(e)
           ==
             fvModels().source(rho, e)
         );

diff --git a/applications/solvers/modules/solid/solid/solid.H b/applications/solvers/modules/solid/solid/solid.H
@@ -38,7 +38,7 @@ SourceFiles
 #define solid_H
 
 #include "solver.H"
-#include "solidThermo.H"
+#include "solidThermophysicalTransportModel.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -74,6 +74,12 @@ protected:
         volScalarField& T;
 
 
+    // Thermophysical transport
+
+        //- Pointer to the solid thermophysical transport model
+        autoPtr<solidThermophysicalTransportModel> thermophysicalTransport;
+
+
     // Time-step controls
 
         scalar DiNum;

diff --git a/...ns/solvers/multiphase/multiphaseEulerFoam/phaseSystems/phaseModel/phaseModel/phaseModel.H b/...ns/solvers/multiphase/multiphaseEulerFoam/phaseSystems/phaseModel/phaseModel/phaseModel.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2015-2021 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2015-2022 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -349,7 +349,7 @@ public:
 
         // Transport
 
-            //- Effective thermal turbulent diffusivity for temperature
+            //- Effective thermal turbulent conductivity
             //  of mixture for patch [W/m/K]
             virtual tmp<scalarField> kappaEff(const label patchi) const = 0;
 

diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/Make/options b/applications/solvers/stressAnalysis/solidDisplacementFoam/Make/options
@@ -6,13 +6,15 @@ EXE_INC = \
     -I$(LIB_SRC)/physicalProperties/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \
+    -I$(LIB_SRC)/ThermophysicalTransportModels/solidThermophysicalTransportModels/lnInclude \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude
 
 EXE_LIBS = \
-    -lfiniteVolume \
     -lsolidThermo \
     -lsolidDisplacementThermo \
+    -lsolidThermophysicalTransportModels \
+    -lfiniteVolume \
+    -lmeshTools \
     -lfvModels \
-    -lfvConstraints \
-    -lmeshTools
+    -lfvConstraints
diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermophysicalProperties.H b/applications/solvers/stressAnalysis/solidDisplacementFoam/readThermophysicalProperties.H
@@ -2,6 +2,13 @@ Info<< "Reading thermophysical properties\n" << endl;
 
 solidDisplacementThermo thermo(mesh);
 
+autoPtr<solidThermophysicalTransportModel> thermophysicalTransport;
+
+if (thermo.thermalStress())
+{
+    thermophysicalTransport = solidThermophysicalTransportModel::New(thermo);
+}
+
 const volScalarField& E(thermo.E());
 const volScalarField& nu(thermo.nu());
 

diff --git a/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C b/applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2011-2021 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2022 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -37,6 +37,7 @@ Description
 
 #include "fvCFD.H"
 #include "solidDisplacementThermo.H"
+#include "solidThermophysicalTransportModel.H"
 #include "fvModels.H"
 #include "fvConstraints.H"
 
@@ -76,7 +77,7 @@ int main(int argc, char *argv[])
                 fvScalarMatrix TEqn
                 (
                     fvm::ddt(rho, Cp, T)
-                  + thermo.divq(T)
+                  + thermophysicalTransport->divq(T)
                  ==
                     fvModels.source(rho*Cp, T)
                 );

diff --git a/...rs/stressAnalysis/solidDisplacementFoam/solidDisplacementThermo/solidDisplacementThermo.H b/...rs/stressAnalysis/solidDisplacementFoam/solidDisplacementThermo/solidDisplacementThermo.H
@@ -61,7 +61,7 @@ class solidDisplacementThermo
         //- Heat capacity at constant pressure [J/kg/K]
         volScalarField Cp_;
 
-        //- Thermal diffusivity for temperature [W/m/K]
+        //- Thermal conductivity [W/m/K]
         volScalarField kappa_;
 
         //- Youngs modulus [Pa]
@@ -176,7 +176,7 @@ public:
 
         // Access to transport state variables
 
-            //- Thermal diffusivity for temperature of mixture [W/m/K]
+            //- Thermal conductivity of mixture [W/m/K]
             virtual const volScalarField& kappa() const;
 
 
@@ -314,7 +314,7 @@ public:
             //- Thermal diffusivity for energy of mixture for patch [kg/m/s]
             virtual tmp<scalarField> alphahe(const label patchi) const;
 
-            //- Effective thermal turbulent diffusivity for temperature
+            //- Effective thermal turbulent conductivity
             //  of mixture [W/m/K]
             virtual tmp<volScalarField> kappaEff
             (
@@ -327,7 +327,7 @@ public:
                 const volScalarField& alphat
             ) const;
 
-            //- Effective thermal turbulent diffusivity for temperature
+            //- Effective thermal turbulent conductivity
             //  of mixture for patch [W/m/K]
             virtual tmp<scalarField> kappaEff
             (

diff --git a/...olvers/stressAnalysis/solidDisplacementFoam/solidEquilibriumDisplacementFoam/Make/options b/...olvers/stressAnalysis/solidDisplacementFoam/solidEquilibriumDisplacementFoam/Make/options
@@ -6,11 +6,13 @@ EXE_INC = \
     -I$(LIB_SRC)/physicalProperties/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \
+    -I$(LIB_SRC)/ThermophysicalTransportModels/solidThermophysicalTransportModels/lnInclude \
     -I$(LIB_SRC)/finiteVolume/lnInclude \
     -I$(LIB_SRC)/meshTools/lnInclude
 
 EXE_LIBS = \
     -lfiniteVolume \
     -lsolidThermo \
     -lsolidDisplacementThermo \
+    -lsolidThermophysicalTransportModels \
     -lmeshTools
diff --git a/...solidDisplacementFoam/solidEquilibriumDisplacementFoam/solidEquilibriumDisplacementFoam.C b/...solidDisplacementFoam/solidEquilibriumDisplacementFoam/solidEquilibriumDisplacementFoam.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2011-2021 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2022 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -26,8 +26,7 @@ Application
 
 Description
     Steady-state segregated finite-volume solver of linear-elastic,
-    small-strain deformation of a solid body, with optional thermal
-    diffusion and thermal stresses.
+    small-strain deformation of a solid body.
 
     Simple linear elasticity structural analysis code.
     Solves for the displacement vector field D, also generating the
@@ -37,6 +36,7 @@ Description
 
 #include "fvCFD.H"
 #include "solidDisplacementThermo.H"
+#include "solidThermophysicalTransportModel.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 

diff --git a/src/OpenFOAM/meshes/polyMesh/polyMesh.H b/src/OpenFOAM/meshes/polyMesh/polyMesh.H
@@ -307,7 +307,7 @@ protected:
             //- Is the mesh moving
             bool moving_;
 
-            //- Is the mesh topology changing
+            //- Has the mesh topology changed
             bool topoChanged_;
 
 
@@ -528,7 +528,7 @@ public:
                 return moving_;
             }
 
-            //- Is mesh topology changing
+            //- Has the mesh topology changed
             bool topoChanged() const
             {
                 return topoChanged_;

diff --git a/src/ThermophysicalTransportModels/Allwmake b/src/ThermophysicalTransportModels/Allwmake
@@ -8,5 +8,7 @@ wmake $targetType
 wmake $targetType fluidMulticomponentThermo
 wmake $targetType phaseFluidThermo
 wmake $targetType phaseFluidMulticomponentThermo
+wmake $targetType solidThermophysicalTransportModels
+wmake $targetType coupledThermophysicalTransportModels
 
 #------------------------------------------------------------------------------
diff --git a/src/ThermophysicalTransportModels/Make/files b/src/ThermophysicalTransportModels/Make/files
@@ -1,9 +1,6 @@
 thermophysicalTransportModel/thermophysicalTransportModel.C
 
 BCs = derivedFvPatchFields
-$(BCs)/externalTemperature/externalTemperatureFvPatchScalarField.C
-$(BCs)/temperatureCoupledBase/temperatureCoupledBase.C
-$(BCs)/coupledTemperature/coupledTemperatureFvPatchScalarField.C
 $(BCs)/thermalBaffle1D/thermalBaffle1DFvPatchScalarFields.C
 $(BCs)/totalFlowRateAdvectiveDiffusive/totalFlowRateAdvectiveDiffusiveFvPatchScalarField.C
 $(BCs)/externalCoupledTemperatureMixed/externalCoupledTemperatureMixedFvPatchScalarField.C

diff --git a/src/ThermophysicalTransportModels/coupledThermophysicalTransportModels/Make/files b/src/ThermophysicalTransportModels/coupledThermophysicalTransportModels/Make/files
@@ -0,0 +1,5 @@
+externalTemperature/externalTemperatureFvPatchScalarField.C
+patchKappa/patchKappa.C
+coupledTemperature/coupledTemperatureFvPatchScalarField.C
+
+LIB = $(FOAM_LIBBIN)/libcoupledThermophysicalTransportModels
diff --git a/src/ThermophysicalTransportModels/coupledThermophysicalTransportModels/Make/options b/src/ThermophysicalTransportModels/coupledThermophysicalTransportModels/Make/options
@@ -0,0 +1,19 @@
+EXE_INC = \
+    -I$(LIB_SRC)/ThermophysicalTransportModels/lnInclude \
+    -I$(LIB_SRC)/ThermophysicalTransportModels/solidThermophysicalTransportModels/lnInclude \
+    -I$(LIB_SRC)/MomentumTransportModels/momentumTransportModels/lnInclude \
+    -I$(LIB_SRC)/MomentumTransportModels/compressible/lnInclude \
+    -I$(LIB_SRC)/physicalProperties/lnInclude \
+    -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
+    -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \
+    -I$(LIB_SRC)/finiteVolume/lnInclude \
+    -I$(LIB_SRC)/meshTools/lnInclude \
+
+LIB_LIBS = \
+    -lfluidThermophysicalModels \
+    -lsolidThermo \
+    -lthermophysicalTransportModels \
+    -lmomentumTransportModels \
+    -lsolidThermophysicalTransportModels \
+    -lfiniteVolume \
+    -lmeshTools