/*---------------------------------------------------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | \\ / A nd | www.openfoam.com \\/ M anipulation | ------------------------------------------------------------------------------- Copyright (C) 2011-2014 OpenFOAM Foundation Copyright (C) 2019-2025 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 . \*---------------------------------------------------------------------------*/ #include "processorFvPatch.H" #include "addToRunTimeSelectionTable.H" #include "transformField.H" // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // namespace Foam { defineTypeNameAndDebug(processorFvPatch, 0); addToRunTimeSelectionTable(fvPatch, processorFvPatch, polyPatch); } // * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * // void Foam::processorFvPatch::makeWeights(scalarField& w) const { if (Pstream::parRun()) { // The face normals point in the opposite direction on the other side scalarField neighbFaceCentresCn ( ( procPolyPatch_.neighbFaceAreas() /(mag(procPolyPatch_.neighbFaceAreas()) + VSMALL) ) & ( procPolyPatch_.neighbFaceCentres() - procPolyPatch_.neighbFaceCellCentres()) ); w = neighbFaceCentresCn /((nf()&coupledFvPatch::delta()) + neighbFaceCentresCn); } else { w = 1.0; } } Foam::tmp Foam::processorFvPatch::delta() const { if (Pstream::parRun()) { // To the transformation if necessary if (parallel()) { return coupledFvPatch::delta() - ( procPolyPatch_.neighbFaceCentres() - procPolyPatch_.neighbFaceCellCentres() ); } else { return coupledFvPatch::delta() - transform ( forwardT(), ( procPolyPatch_.neighbFaceCentres() - procPolyPatch_.neighbFaceCellCentres() ) ); } } else { return coupledFvPatch::delta(); } } Foam::tmp Foam::processorFvPatch::interfaceInternalField ( const labelUList& internalData ) const { return patchInternalField(internalData); } Foam::tmp Foam::processorFvPatch::interfaceInternalField ( const labelUList& internalData, const labelUList& faceCells ) const { auto tpfld = tmp::New(this->size()); patchInternalField(internalData, faceCells, tpfld.ref()); return tpfld; } void Foam::processorFvPatch::initInternalFieldTransfer ( const Pstream::commsTypes commsType, const labelUList& iF ) const { send(commsType, patchInternalField(iF)()); } void Foam::processorFvPatch::initInternalFieldTransfer ( const Pstream::commsTypes commsType, const labelUList& iF, const labelUList& faceCells ) const { send(commsType, interfaceInternalField(iF, faceCells)()); } Foam::tmp Foam::processorFvPatch::internalFieldTransfer ( const Pstream::commsTypes commsType, const labelUList& ) const { return receive