diff --git a/src/regionModels/thermalBaffle/Make/files b/src/regionModels/thermalBaffle/Make/files deleted file mode 100644 index 08ca06e864..0000000000 --- a/src/regionModels/thermalBaffle/Make/files +++ /dev/null @@ -1,5 +0,0 @@ -thermalBaffle/thermalBaffle.C - -derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.C - -LIB = $(FOAM_LIBBIN)/libthermalBaffle diff --git a/src/regionModels/thermalBaffle/Make/options b/src/regionModels/thermalBaffle/Make/options deleted file mode 100644 index b4e70826a2..0000000000 --- a/src/regionModels/thermalBaffle/Make/options +++ /dev/null @@ -1,26 +0,0 @@ -EXE_INC = \ - -I$(LIB_SRC)/regionModels/regionModel/lnInclude \ - -I$(LIB_SRC)/physicalProperties/lnInclude \ - -I$(LIB_SRC)/thermophysicalModels/solidThermo/lnInclude \ - -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \ - -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \ - -I$(LIB_SRC)/radiationModels/lnInclude \ - -I$(LIB_SRC)/ThermophysicalTransportModels/thermophysicalTransportModel/lnInclude \ - -I$(LIB_SRC)/ThermophysicalTransportModels/fluid/lnInclude \ - -I$(LIB_SRC)/ThermophysicalTransportModels/solid/lnInclude \ - -I$(LIB_SRC)/ThermophysicalTransportModels/coupledThermophysicalTransportModels/lnInclude \ - -I$(LIB_SRC)/dynamicMesh/lnInclude \ - -I$(LIB_SRC)/mesh/extrudeModel/lnInclude \ - -I$(LIB_SRC)/finiteVolume/lnInclude \ - -I$(LIB_SRC)/meshTools/lnInclude - -LIB_LIBS = \ - -lregionModels \ - -lsolidThermo \ - -lmomentumTransportModels \ - -lsolidThermophysicalTransportModels \ - -lcoupledThermophysicalTransportModels \ - -lradiationModels \ - -ldynamicMesh \ - -lfiniteVolume \ - -lmeshTools diff --git a/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.C b/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.C deleted file mode 100644 index 9937ffd1a5..0000000000 --- a/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.C +++ /dev/null @@ -1,410 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see . - -\*---------------------------------------------------------------------------*/ - -#include "thermalBaffleFvPatchScalarField.H" -#include "mappedExtrudedWallPolyPatch.H" -#include "symmetryPolyPatch.H" -#include "addToRunTimeSelectionTable.H" - -// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * / - -bool Foam::thermalBaffleFvPatchScalarField::primary() const -{ - return patch().boundaryMesh().mesh().name() == polyMesh::defaultRegion; -} - - -bool Foam::thermalBaffleFvPatchScalarField::owner() const -{ - return - primary() - && patch().index() < patch().boundaryMesh()[nbrPatch_].index(); -} - - -void Foam::thermalBaffleFvPatchScalarField::checkPatches() const -{ - if (!primary()) return; - - const polyPatch& pp = patch().patch(); - const polyPatch& nbrPp = patch().patch().boundaryMesh()[nbrPatch_]; - - // The patches should be of mapped type - auto checkPatchIsMapped = [&](const polyPatch& pp) - { - if (!isA(pp)) - { - FatalErrorInFunction - << "Patch field of type \"" << typeName - << "\" specified for patch \"" << pp.name() << "\" of field \"" - << internalField().name() << "\", but this patch is not of " - << "type \"" << mappedPatchBase::typeName << "\"" - << exit(FatalError); - } - }; - checkPatchIsMapped(pp); - checkPatchIsMapped(nbrPp); - - const mappedPatchBase& mpb = refCast(pp); - const mappedPatchBase& nbrMpb = refCast(nbrPp); - - // The patches should neighbour a different region - auto checkPatchMapsDifferentRegion = [&](const mappedPatchBase& mpb) - { - if (mpb.sameRegion()) - { - FatalErrorInFunction - << "Patch field of type \"" << typeName - << "\" specified for patch \"" << pp.name() << "\" of field \"" - << internalField().name() << "\", but this patch maps to " - << "another patch in the same region. It should map to a " - << "different region; i.e., that of the thermal baffle model." - << exit(FatalError); - } - }; - checkPatchMapsDifferentRegion(mpb); - checkPatchMapsDifferentRegion(nbrMpb); - - // The neighbour region of this patch and it's neighbour should be the - // same, i.e., that of the thermal baffle model - if (mpb.nbrRegionName() != nbrMpb.nbrRegionName()) - { - FatalErrorInFunction - << "Patch fields of type \"" << typeName - << "\" specified for patches \"" << pp.name() << "\" and \"" - << nbrPp.name() << "\" of field \"" << internalField().name() - << "\", but these patches map to different regions \"" - << mpb.nbrRegionName() << "\" and \"" - << nbrMpb.nbrRegionName() << ". They should map to the same " - << "region; i.e., that of the thermal baffle model." - << exit(FatalError); - } - - // The neighbour patch of this patch and it's neighbour should be - // different, i.e., they should map opposite ends of the thermal baffle - // mesh - if (mpb.nbrPatchName() == nbrMpb.nbrPatchName()) - { - FatalErrorInFunction - << "Patch fields of type \"" << typeName - << "\" specified for patches \"" << pp.name() << "\" and \"" - << nbrPp.name() << "\" of field \"" << internalField().name() - << "\", but these patches map to the same patch; \"" - << mpb.nbrPatchName() << "\" of region \"" - << mpb.nbrRegionName() << ". They should map to different " - << "patches, as these will become the patches at opposite ends of " - << "the extruded baffle mesh." << exit(FatalError); - } -} - - -void Foam::thermalBaffleFvPatchScalarField::checkPatchFields() const -{ - if (!primary()) return; - - const fvPatch& fvp = patch(); - const fvPatch& nbrFvp = patch().boundaryMesh()[nbrPatch_]; - - const fvPatchScalarField& nbrTp = - nbrFvp.lookupPatchField(internalField().name()); - - // The neighbour patch field should be of the same type - if (!isA(nbrTp)) - { - FatalErrorInFunction - << "Patch field of type \"" << typeName - << "\" specified for patch \"" << fvp.name() << "\" of field \"" - << internalField().name() << "\" but the field on the " - << "neighbouring patch \"" << nbrFvp.name() - << "\" is of a different type. Both should be of type \"" - << typeName << "\"." - << exit(FatalError); - } - - // The neighbour patch field's neighbour should be this patch - const thermalBaffleFvPatchScalarField& nbrTBp = - refCast(nbrTp); - if (nbrTBp.nbrPatch_ != patch().name()) - { - FatalErrorInFunction - << "Patch field of type \"" << typeName - << "\" on patch \"" << fvp.name() << "\" of field \"" - << internalField().name() << "\" is specified to neighbour " - << "patch \"" << nbrPatch_ << "\", but this patch does not " - << "reciprocally neighbour patch \"" << fvp.name() << "\"" - << exit(FatalError); - } -} - - -Foam::autoPtr -Foam::thermalBaffleFvPatchScalarField::initBaffleMesh() const -{ - if (!owner()) - { - FatalErrorInFunction - << "Baffle mesh is only available to the owner patch in the " - << "primary region" << exit(FatalError); - } - - checkPatches(); - - const fvMesh& mesh = patch().boundaryMesh().mesh(); - - const mappedPatchBase& mpb = - refCast(patch().patch()); - - const mappedPatchBase nbrMpb = - refCast - (patch().patch().boundaryMesh()[nbrPatch_]); - - const List patchNames - ({ - mpb.nbrPatchName(), - nbrMpb.nbrPatchName(), - "sides" - }); - - const List patchTypes - ({ - mappedWallPolyPatch::typeName, - mappedExtrudedWallPolyPatch::typeName, - symmetryPolyPatch::typeName - }); - - List patchDicts(3); - forAll(patchDicts, patchi) - { - patchDicts[patchi].set("nFaces", 0); - patchDicts[patchi].set("startFace", 0); - } - patchDicts[0].add("neighbourRegion", mesh.name()); - patchDicts[0].add("neighbourPatch", patch().name()); - patchDicts[1].add("neighbourRegion", mesh.name()); - patchDicts[1].add("neighbourPatch", nbrPatch_); - patchDicts[1].add("bottomPatch", patchNames[0]); - - List patchPtrs(3); - forAll(patchPtrs, patchi) - { - patchPtrs[patchi] = polyPatch::New - ( - patchTypes[patchi], - patchNames[patchi], - patchDicts[patchi], - patchi, - mesh.boundaryMesh() - ).ptr(); - } - - dictionary dict(dict_); - dict.add("columnCells", false); - - return - autoPtr - ( - new extrudePatchMesh - ( - mesh, - patch(), - dict, - mpb.nbrRegionName(), - patchPtrs - ) - ); -} - - -Foam::autoPtr -Foam::thermalBaffleFvPatchScalarField::initBaffle() const -{ - if (!owner()) - { - FatalErrorInFunction - << "Baffle model is only available to the owner patch in the " - << "primary region" << exit(FatalError); - } - - checkPatches(); - - const fvMesh& mesh = patch().boundaryMesh().mesh(); - - const mappedPatchBase& mpb = - refCast(patch().patch()); - - dictionary dict(dict_); - dict.add("regionName", mpb.nbrRegionName()); - - return autoPtr - ( - new regionModels::thermalBaffle("thermalBaffle", mesh, dict) - ); -} - - -// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // - -Foam::thermalBaffleFvPatchScalarField::thermalBaffleFvPatchScalarField -( - const fvPatch& p, - const DimensionedField& iF -) -: - coupledTemperatureFvPatchScalarField(p, iF), - dict_(dictionary::null), - nbrPatch_(word::null), - baffleMeshPtr_(), - bafflePtr_() -{} - - -Foam::thermalBaffleFvPatchScalarField::thermalBaffleFvPatchScalarField -( - const fvPatch& p, - const DimensionedField& iF, - const dictionary& dict -) -: - coupledTemperatureFvPatchScalarField(p, iF, dict), - dict_(dict), - nbrPatch_(primary() ? dict.lookup("neighbourPatch") : word::null), - baffleMeshPtr_(owner() ? initBaffleMesh().ptr() : nullptr), - bafflePtr_(owner() ? initBaffle().ptr() : nullptr) -{} - - -Foam::thermalBaffleFvPatchScalarField::thermalBaffleFvPatchScalarField -( - const thermalBaffleFvPatchScalarField& ptf, - const fvPatch& p, - const DimensionedField& iF, - const fvPatchFieldMapper& mapper -) -: - coupledTemperatureFvPatchScalarField - ( - ptf, - p, - iF, - mapper - ), - dict_(ptf.dict_), - baffleMeshPtr_(), - bafflePtr_() -{} - - -Foam::thermalBaffleFvPatchScalarField::thermalBaffleFvPatchScalarField -( - const thermalBaffleFvPatchScalarField& ptf, - const DimensionedField& iF -) -: - coupledTemperatureFvPatchScalarField(ptf, iF), - dict_(ptf.dict_), - baffleMeshPtr_(), - bafflePtr_() -{} - - -// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // - -void Foam::thermalBaffleFvPatchScalarField::autoMap -( - const fvPatchFieldMapper& m -) -{ - coupledTemperatureFvPatchScalarField::autoMap(m); -} - - -void Foam::thermalBaffleFvPatchScalarField::rmap -( - const fvPatchScalarField& ptf, - const labelList& addr -) -{ - coupledTemperatureFvPatchScalarField::rmap(ptf, addr); -} - - -void Foam::thermalBaffleFvPatchScalarField::reset -( - const fvPatchScalarField& ptf -) -{ - coupledTemperatureFvPatchScalarField::reset(ptf); -} - - -void Foam::thermalBaffleFvPatchScalarField::updateCoeffs() -{ - if (this->updated()) - { - return; - } - - checkPatchFields(); - - if (owner()) - { - bafflePtr_->evolve(); - } - - coupledTemperatureFvPatchScalarField::updateCoeffs(); -} - - -void Foam::thermalBaffleFvPatchScalarField::write(Ostream& os) const -{ - coupledTemperatureFvPatchScalarField::write(os); - - if (owner()) - { - forAllConstIter(dictionary, dict_, iter) - { - os << *iter; - } - } - else if (primary()) - { - writeEntry(os, "neighbourPatch", nbrPatch_); - } -} - - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -namespace Foam -{ - makePatchTypeField - ( - fvPatchScalarField, - thermalBaffleFvPatchScalarField - ); -} - - -// ************************************************************************* // diff --git a/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.H b/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.H deleted file mode 100644 index af696162d4..0000000000 --- a/src/regionModels/thermalBaffle/derivedFvPatchFields/thermalBaffle/thermalBaffleFvPatchScalarField.H +++ /dev/null @@ -1,287 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see . - -Class - Foam::compressible::thermalBaffleFvPatchScalarField - -Description - This boundary condition provides a coupled temperature condition between - multiple mesh regions. - - The regions are generally referred to as the: - - primary region, - - and baffle region. - - The primary region creates the baffle region and evolves its energy - equation either: - - 1-D, normal to each patch face - - 2-D, normal and tangential components - - The thermodynamic properties of the baffle material are specified via - dictionary entries on the master patch. - -Usage - Example of the boundary condition specification: - \verbatim - - { - type compressible::thermalBaffle; - - // Underlying coupled boundary condition - Tnbr T; - qrNbr qr;//or none.Name of qr field on neighbourregion - qr none;// or none.Name of qr field on localregion - value uniform 300; - - // Baffle region name - regionName baffleRegion; - - // Solid thermo in solid region - thermoType - { - type heSolidThermo; - mixture pureMixture; - transport constIso; - thermo eConst; - equationOfState rhoConst; - specie specie; - energy sensibleInternalEnergy; - } - - mixture - { - specie - { - molWeight 20; - } - transport - { - kappa 0.01; - } - thermodynamics - { - Hf 0; - Cv 15; - } - density - { - rho 80; - } - } - - radiation - { - radiationModel opaqueSolid; - absorptionEmissionModel none; - scatterModel none; - } - - // Extrude model for new region - - extrudeModel linearNormal; - nLayers 50; - expansionRatio 1; - columnCells false; //3D or 1D - linearNormalCoeffs - { - thickness 0.02; - } - - } - \endverbatim - - Slave patch on primary region: - \verbatim - - { - type compressible::thermalBaffle; - value uniform 300; - \endverbatim - - Patches on baffle region: - \verbatim - bottom - { - type compressible::thermalBaffle; - value uniform 300; - } - - top - { - type compressible::thermalBaffle; - value uniform 300; - } - \endverbatim - -See also - Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField - Foam::regionModels::thermalBaffle - -SourceFiles - thermalBaffleFvPatchScalarField.C - -\*---------------------------------------------------------------------------*/ - -#ifndef thermalBaffleFvPatchScalarField_H -#define thermalBaffleFvPatchScalarField_H - -#include "thermalBaffle.H" -#include "extrudePatchMesh.H" -#include "coupledTemperatureFvPatchScalarField.H" - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -namespace Foam -{ - -/*---------------------------------------------------------------------------*\ - Class thermalBaffleFvPatchScalarField Declaration -\*---------------------------------------------------------------------------*/ - -class thermalBaffleFvPatchScalarField -: - public coupledTemperatureFvPatchScalarField -{ - // Private Data - - //- Dictionary - dictionary dict_; - - //- The neighbouring baffle patch name - const word nbrPatch_; - - //- Thermal baffle mesh extruded from this patch - autoPtr baffleMeshPtr_; - - //- Thermal baffle model - autoPtr bafflePtr_; - - - // Private Member Functions - - //- Is this patch in the primary region? - bool primary() const; - - //- Is this patch the owner? - bool owner() const; - - //- Check consistency of the patch settings - void checkPatches() const; - - //- Check consistency of the patch field settings - void checkPatchFields() const; - - //- Construct and return the thermal baffle mesh - autoPtr initBaffleMesh() const; - - //- Construct and return the thermal baffle model - autoPtr initBaffle() const; - - -public: - - //- Runtime type information - TypeName("compressible::thermalBaffle"); - - - // Constructors - - //- Construct from patch and internal field - thermalBaffleFvPatchScalarField - ( - const fvPatch&, - const DimensionedField& - ); - - //- Construct from patch, internal field and dictionary - thermalBaffleFvPatchScalarField - ( - const fvPatch&, - const DimensionedField&, - const dictionary& - ); - - //- Construct by mapping given - // thermalBaffleFvPatchScalarField onto a new patch - thermalBaffleFvPatchScalarField - ( - const thermalBaffleFvPatchScalarField&, - const fvPatch&, - const DimensionedField&, - const fvPatchFieldMapper& - ); - - //- Construct as copy setting internal field reference - thermalBaffleFvPatchScalarField - ( - const thermalBaffleFvPatchScalarField&, - const DimensionedField& - ); - - //- Construct and return a clone setting internal field reference - virtual tmp clone - ( - const DimensionedField& iF - ) const - { - return tmp - ( - new thermalBaffleFvPatchScalarField(*this, iF) - ); - } - - - // Member Functions - - // Mapping functions - - //- Map (and resize as needed) from self given a mapping object - // Used to update fields following mesh topology change - virtual void autoMap(const fvPatchFieldMapper&); - - //- Reverse map the given fvPatchField onto this fvPatchField - // Used to reconstruct fields - virtual void rmap(const fvPatchScalarField&, const labelList&); - - //- Reset the fvPatchField to the given fvPatchField - // Used for mesh to mesh mapping - virtual void reset(const fvPatchScalarField&); - - - //- Update the coefficients associated with the patch field - virtual void updateCoeffs(); - - //- Write - virtual void write(Ostream&) const; -}; - - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -} // End namespace Foam - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -#endif - -// ************************************************************************* // diff --git a/src/regionModels/thermalBaffle/thermalBaffle/thermalBaffle.C b/src/regionModels/thermalBaffle/thermalBaffle/thermalBaffle.C deleted file mode 100644 index 60cd9c7b21..0000000000 --- a/src/regionModels/thermalBaffle/thermalBaffle/thermalBaffle.C +++ /dev/null @@ -1,318 +0,0 @@ -/*---------------------------------------------------------------------------*\ - ========= | - \\ / F ield | OpenFOAM: The Open Source CFD Toolbox - \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2011-2022 OpenFOAM Foundation - \\/ M anipulation | -------------------------------------------------------------------------------- -License - This file is part of OpenFOAM. - - OpenFOAM is free software: you can redistribute it and/or modify it - under the terms of the GNU General Public License as published by - the Free Software Foundation, either version 3 of the License, or - (at your option) any later version. - - OpenFOAM is distributed in the hope that it will be useful, but WITHOUT - ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or - FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License - for more details. - - You should have received a copy of the GNU General Public License - along with OpenFOAM. If not, see . - -\*---------------------------------------------------------------------------*/ - -#include "thermalBaffle.H" -#include "fvm.H" -#include "fvcDiv.H" -#include "absorptionEmissionModel.H" -#include "zeroGradientFvPatchFields.H" -#include "wedgePolyPatch.H" -#include "emptyPolyPatch.H" -#include "FaceCellWave.H" -#include "LayerInfoData.H" - -// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // - -namespace Foam -{ -namespace regionModels -{ - -// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // - -defineTypeNameAndDebug(thermalBaffle, 0); - - -// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * // - -tmp thermalBaffle::calcThickness() const -{ - if (intCoupledPatchIDs_.size() != 2) - { - FatalErrorInFunction - << "Mesh \"" << regionMesh().name() - << "\" does not have exactly two coupled patches" - << exit(FatalError); - } - - // Initial faces from which to wave - DynamicList