mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-12-28 03:37:59 +00:00
8a8b5db977
- manual revert for 44d0fbd594.
Unless required for post-processing the "value" entry should not
written if it is not mandatory for reading.
This is especially crucial with an 'empty' patch field where the
field (size 0) has a different size from the patch.
- the changes made align fvsPatchField behaviour with fvPatchField
329 lines
8.7 KiB
C
329 lines
8.7 KiB
C
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration |
|
|
\\ / A nd | www.openfoam.com
|
|
\\/ M anipulation |
|
|
-------------------------------------------------------------------------------
|
|
Copyright (C) 2011-2017 OpenFOAM Foundation
|
|
Copyright (C) 2019 OpenCFD Ltd.
|
|
-------------------------------------------------------------------------------
|
|
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 <http://www.gnu.org/licenses/>.
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "fvMatrix.H"
|
|
#include "cyclicFvPatchField.H"
|
|
#include "transformField.H"
|
|
#include "volFields.H"
|
|
|
|
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
|
|
|
template<class Type>
|
|
Foam::cyclicFvPatchField<Type>::cyclicFvPatchField
|
|
(
|
|
const fvPatch& p,
|
|
const DimensionedField<Type, volMesh>& iF
|
|
)
|
|
:
|
|
coupledFvPatchField<Type>(p, iF),
|
|
cyclicPatch_(refCast<const cyclicFvPatch>(p))
|
|
{}
|
|
|
|
|
|
template<class Type>
|
|
Foam::cyclicFvPatchField<Type>::cyclicFvPatchField
|
|
(
|
|
const fvPatch& p,
|
|
const DimensionedField<Type, volMesh>& iF,
|
|
const dictionary& dict,
|
|
const bool needValue
|
|
)
|
|
:
|
|
coupledFvPatchField<Type>(p, iF, dict, IOobjectOption::NO_READ),
|
|
cyclicPatch_(refCast<const cyclicFvPatch>(p, dict))
|
|
{
|
|
if (!isA<cyclicFvPatch>(p))
|
|
{
|
|
FatalIOErrorInFunction(dict)
|
|
<< " patch type '" << p.type()
|
|
<< "' not constraint type '" << typeName << "'"
|
|
<< "\n for patch " << p.name()
|
|
<< " of field " << this->internalField().name()
|
|
<< " in file " << this->internalField().objectPath()
|
|
<< exit(FatalIOError);
|
|
}
|
|
|
|
if (needValue)
|
|
{
|
|
this->evaluate(Pstream::commsTypes::blocking);
|
|
}
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
Foam::cyclicFvPatchField<Type>::cyclicFvPatchField
|
|
(
|
|
const cyclicFvPatchField<Type>& ptf,
|
|
const fvPatch& p,
|
|
const DimensionedField<Type, volMesh>& iF,
|
|
const fvPatchFieldMapper& mapper
|
|
)
|
|
:
|
|
coupledFvPatchField<Type>(ptf, p, iF, mapper),
|
|
cyclicPatch_(refCast<const cyclicFvPatch>(p))
|
|
{
|
|
if (!isA<cyclicFvPatch>(this->patch()))
|
|
{
|
|
FatalErrorInFunction
|
|
<< "\n patch type '" << p.type()
|
|
<< "' not constraint type '" << typeName << "'"
|
|
<< "\n for patch " << p.name()
|
|
<< " of field " << this->internalField().name()
|
|
<< " in file " << this->internalField().objectPath()
|
|
<< exit(FatalError);
|
|
}
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
Foam::cyclicFvPatchField<Type>::cyclicFvPatchField
|
|
(
|
|
const cyclicFvPatchField<Type>& ptf
|
|
)
|
|
:
|
|
cyclicLduInterfaceField(),
|
|
coupledFvPatchField<Type>(ptf),
|
|
cyclicPatch_(ptf.cyclicPatch_)
|
|
{}
|
|
|
|
|
|
template<class Type>
|
|
Foam::cyclicFvPatchField<Type>::cyclicFvPatchField
|
|
(
|
|
const cyclicFvPatchField<Type>& ptf,
|
|
const DimensionedField<Type, volMesh>& iF
|
|
)
|
|
:
|
|
coupledFvPatchField<Type>(ptf, iF),
|
|
cyclicPatch_(ptf.cyclicPatch_)
|
|
{}
|
|
|
|
|
|
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
|
|
|
template<class Type>
|
|
Foam::tmp<Foam::Field<Type>>
|
|
Foam::cyclicFvPatchField<Type>::patchNeighbourField() const
|
|
{
|
|
const Field<Type>& iField = this->primitiveField();
|
|
const labelUList& nbrFaceCells =
|
|
cyclicPatch().cyclicPatch().neighbPatch().faceCells();
|
|
|
|
auto tpnf = tmp<Field<Type>>::New(this->size());
|
|
auto& pnf = tpnf.ref();
|
|
|
|
|
|
if (doTransform())
|
|
{
|
|
forAll(pnf, facei)
|
|
{
|
|
pnf[facei] = transform
|
|
(
|
|
forwardT()[0], iField[nbrFaceCells[facei]]
|
|
);
|
|
}
|
|
}
|
|
else
|
|
{
|
|
forAll(pnf, facei)
|
|
{
|
|
pnf[facei] = iField[nbrFaceCells[facei]];
|
|
}
|
|
}
|
|
|
|
return tpnf;
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
const Foam::cyclicFvPatchField<Type>&
|
|
Foam::cyclicFvPatchField<Type>::neighbourPatchField() const
|
|
{
|
|
const GeometricField<Type, fvPatchField, volMesh>& fld =
|
|
static_cast<const GeometricField<Type, fvPatchField, volMesh>&>
|
|
(
|
|
this->primitiveField()
|
|
);
|
|
|
|
return refCast<const cyclicFvPatchField<Type>>
|
|
(
|
|
fld.boundaryField()[this->cyclicPatch().neighbPatchID()]
|
|
);
|
|
}
|
|
|
|
|
|
|
|
template<class Type>
|
|
void Foam::cyclicFvPatchField<Type>::updateInterfaceMatrix
|
|
(
|
|
solveScalarField& result,
|
|
const bool add,
|
|
const lduAddressing& lduAddr,
|
|
const label patchId,
|
|
const solveScalarField& psiInternal,
|
|
const scalarField& coeffs,
|
|
const direction cmpt,
|
|
const Pstream::commsTypes commsType
|
|
) const
|
|
{
|
|
const labelUList& nbrFaceCells =
|
|
lduAddr.patchAddr
|
|
(
|
|
this->cyclicPatch().neighbPatchID()
|
|
);
|
|
|
|
solveScalarField pnf(psiInternal, nbrFaceCells);
|
|
|
|
// Transform according to the transformation tensors
|
|
transformCoupleField(pnf, cmpt);
|
|
|
|
|
|
const labelUList& faceCells = lduAddr.patchAddr(patchId);
|
|
|
|
// Multiply the field by coefficients and add into the result
|
|
this->addToInternalField(result, !add, faceCells, coeffs, pnf);
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
void Foam::cyclicFvPatchField<Type>::updateInterfaceMatrix
|
|
(
|
|
Field<Type>& result,
|
|
const bool add,
|
|
const lduAddressing& lduAddr,
|
|
const label patchId,
|
|
const Field<Type>& psiInternal,
|
|
const scalarField& coeffs,
|
|
const Pstream::commsTypes
|
|
) const
|
|
{
|
|
const labelList& nbrFaceCells =
|
|
lduAddr.patchAddr
|
|
(
|
|
this->cyclicPatch().neighbPatchID()
|
|
);
|
|
|
|
Field<Type> pnf(psiInternal, nbrFaceCells);
|
|
|
|
// Transform according to the transformation tensors
|
|
transformCoupleField(pnf);
|
|
|
|
const labelUList& faceCells = lduAddr.patchAddr(patchId);
|
|
|
|
// Multiply the field by coefficients and add into the result
|
|
this->addToInternalField(result, !add, faceCells, coeffs, pnf);
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
void Foam::cyclicFvPatchField<Type>::write(Ostream& os) const
|
|
{
|
|
fvPatchField<Type>::write(os);
|
|
}
|
|
|
|
|
|
template<class Type>
|
|
void Foam::cyclicFvPatchField<Type>::manipulateMatrix
|
|
(
|
|
fvMatrix<Type>& matrix,
|
|
const label mat,
|
|
const direction cmpt
|
|
)
|
|
{
|
|
if (this->cyclicPatch().owner())
|
|
{
|
|
label index = this->patch().index();
|
|
|
|
const label globalPatchID =
|
|
matrix.lduMeshAssembly().patchLocalToGlobalMap()[mat][index];
|
|
|
|
const Field<scalar> intCoeffsCmpt
|
|
(
|
|
matrix.internalCoeffs()[globalPatchID].component(cmpt)
|
|
);
|
|
|
|
const Field<scalar> boundCoeffsCmpt
|
|
(
|
|
matrix.boundaryCoeffs()[globalPatchID].component(cmpt)
|
|
);
|
|
|
|
const labelUList& u = matrix.lduAddr().upperAddr();
|
|
const labelUList& l = matrix.lduAddr().lowerAddr();
|
|
|
|
const labelList& faceMap =
|
|
matrix.lduMeshAssembly().faceBoundMap()[mat][index];
|
|
|
|
forAll (faceMap, faceI)
|
|
{
|
|
label globalFaceI = faceMap[faceI];
|
|
|
|
const scalar boundCorr = -boundCoeffsCmpt[faceI];
|
|
const scalar intCorr = -intCoeffsCmpt[faceI];
|
|
|
|
matrix.upper()[globalFaceI] += boundCorr;
|
|
matrix.diag()[u[globalFaceI]] -= boundCorr;
|
|
matrix.diag()[l[globalFaceI]] -= intCorr;
|
|
|
|
if (matrix.asymmetric())
|
|
{
|
|
matrix.lower()[globalFaceI] += intCorr;
|
|
}
|
|
}
|
|
|
|
if (matrix.psi(mat).mesh().fluxRequired(this->internalField().name()))
|
|
{
|
|
matrix.internalCoeffs().set
|
|
(
|
|
globalPatchID, intCoeffsCmpt*pTraits<Type>::one
|
|
);
|
|
matrix.boundaryCoeffs().set
|
|
(
|
|
globalPatchID, boundCoeffsCmpt*pTraits<Type>::one
|
|
);
|
|
|
|
const label nbrPathID = this->cyclicPatch().neighbPatchID();
|
|
|
|
const label nbrGlobalPatchID =
|
|
matrix.lduMeshAssembly().patchLocalToGlobalMap()[mat][nbrPathID];
|
|
|
|
matrix.internalCoeffs().set
|
|
(
|
|
nbrGlobalPatchID, intCoeffsCmpt*pTraits<Type>::one
|
|
|
|
);
|
|
matrix.boundaryCoeffs().set
|
|
(
|
|
nbrGlobalPatchID, boundCoeffsCmpt*pTraits<Type>::one
|
|
);
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|