zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge branch 'issue-2005-objectiveManager-masked-function' into 'develop'

Browse Source 
COMP:  objective manager masks regIOobject::write function (#2005)

See merge request Development/openfoam!427
This commit is contained in:
Mark Olesen
2021-03-05 12:36:06 +00:00
parent 175d35de22 50d21b32f2
commit 22832b87cb
4 changed files with 7 additions and 228 deletions
Download Patch File Download Diff File Expand all files Collapse all files

167
src/optimisation/adjointOptimisation/adjoint/adjointBoundaryConditions/adjointBoundaryCondition/adjointBoundaryConditionTemplates.C
View File

@ -1,167 +0,0 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2007-2019 PCOpt/NTUA
Copyright (C) 2013-2019 FOSS GP
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 "emptyFvPatch.H"
#include "adjointBoundaryCondition.H"
#include "adjointSolverManager.H"
#include "HashTable.H"
#include "surfaceInterpolationScheme.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
template<class Type>
tmp
<
Field<typename Foam::outerProduct<Foam::vector, Type>::type>
>
adjointBoundaryCondition::computePatchGrad(word name)
{
// Return field
typedef typename outerProduct<vector, Type>::type GradType;
auto tresGrad = tmp<Field<GradType>>::New(patch_.size(), Zero);
auto& resGrad = tresGrad.ref();
const labelList& faceCells = patch_.faceCells();
const fvMesh& mesh = patch_.boundaryMesh().mesh();
const cellList& cells = mesh.cells();
// Go through the surfaceInterpolation scheme defined in gradSchemes for
// consistency
const GeometricField<Type, fvPatchField, volMesh>& field =
mesh.lookupObject<volVectorField>(name);
// Gives problems when grad(AdjointVar) is computed using a limited scheme,
// since it is not possible to know a priori how many words to expect in the
// stream.
// Interpolation scheme is now read through interpolation schemes.
/*
word gradSchemeName ("grad(" + name + ')');
Istream& is = mesh.gradScheme(gradSchemeName);
word schemeData(is);
*/
tmp<surfaceInterpolationScheme<Type>> tinterpScheme
(
surfaceInterpolationScheme<Type>::New
(
mesh,
mesh.interpolationScheme("interpolate(" + name + ")")
)
);
GeometricField<Type, fvsPatchField, surfaceMesh> surfField
(
tinterpScheme().interpolate(field)
);
// Auxiliary fields
const surfaceVectorField& Sf = mesh.Sf();
tmp<vectorField> tnf = patch_.nf();
const vectorField& nf = tnf();
const scalarField& V = mesh.V();
const labelUList& owner = mesh.owner();
// Compute grad value of cell adjacent to the boundary
forAll(faceCells, fI)
{
const label cI = faceCells[fI];
const cell& cellI = cells[cI];
for (const label faceI : cellI) // global face numbering
{
label patchID = mesh.boundaryMesh().whichPatch(faceI);
if (patchID == -1) //face is internal
{
const label own = owner[faceI];
tensor flux = Sf[faceI]*surfField[faceI];
if (cI == own)
{
resGrad[fI] += flux;
}
else
{
resGrad[fI] -= flux;
}
}
else // Face is boundary. Covers coupled patches as well
{
if (!isA<emptyFvPatch>(mesh.boundary()[patchID]))
{
const fvPatch& patchForFlux = mesh.boundary()[patchID];
const label boundaryFaceI = faceI - patchForFlux.start();
const vectorField& Sfb = Sf.boundaryField()[patchID];
resGrad[fI] +=
Sfb[boundaryFaceI]
*surfField.boundaryField()[patchID][boundaryFaceI];
}
}
}
resGrad[fI] /= V[cI];
}
// This has concluded the computation of the grad at the cell next to the
// boundary. We now need to compute the grad at the boundary face
const fvPatchField<Type>& bField = field.boundaryField()[patch_.index()];
resGrad = nf*bField.snGrad() + (resGrad - nf*(nf & resGrad));
return tresGrad;
}
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
template<class Type>
adjointBoundaryCondition::adjointBoundaryCondition
(
const fvPatch& p,
const DimensionedField<Type, volMesh>& iF,
const word& solverName
)
:
patch_(p),
managerName_("objectiveManager" + solverName),
adjointSolverName_(solverName),
simulationType_("incompressible"),
boundaryContrPtr_(nullptr),
addATCUaGradUTerm_(nullptr)
{
// Set the boundaryContribution pointer
setBoundaryContributionPtr();
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// ************************************************************************* //

54
src/optimisation/adjointOptimisation/adjoint/adjointBoundaryConditions/adjointZeroInlet/adjointZeroInletFvPatchFieldsFwd.C
View File

@ -1,54 +0,0 @@
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2007-2019 PCOpt/NTUA
Copyright (C) 2013-2019 FOSS GP
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/>.
\*---------------------------------------------------------------------------*/
#ifndef adjointZeroInletFvPatchFieldsFwd_H
#define adjointZeroInletFvPatchFieldsFwd_H
#include "fieldTypes.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
template<class Type> class adjointZeroInletFvPatchField;
makePatchTypeFieldTypedefs(adjointZeroInlet);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //

8
src/optimisation/adjointOptimisation/adjoint/objectiveManager/objectiveManager/objectiveManager.C
View File

@ -5,8 +5,8 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2007-2020 PCOpt/NTUA
Copyright (C) 2013-2020 FOSS GP
Copyright (C) 2007-2021 PCOpt/NTUA
Copyright (C) 2013-2021 FOSS GP
Copyright (C) 2019 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
@ -252,7 +252,7 @@ scalar objectiveManager::print()
}
bool objectiveManager::write
bool objectiveManager::writeObjectives
(
const scalar weightedObjective,
const bool valid
@ -289,7 +289,7 @@ void objectiveManager::updateAndWrite()
updateNormalizationFactor();
update();
scalar weightedObjective = print();
write(weightedObjective);
writeObjectives(weightedObjective);
}

6
src/optimisation/adjointOptimisation/adjoint/objectiveManager/objectiveManager/objectiveManager.H
View File

@ -5,8 +5,8 @@
\\ / A nd | www.openfoam.com
\\/ M anipulation |
-------------------------------------------------------------------------------
Copyright (C) 2007-2020 PCOpt/NTUA
Copyright (C) 2013-2020 FOSS GP
Copyright (C) 2007-2021 PCOpt/NTUA
Copyright (C) 2013-2021 FOSS GP
Copyright (C) 2019 OpenCFD Ltd.
-------------------------------------------------------------------------------
License
@ -146,7 +146,7 @@ public:
scalar print();
//- Write objective function history
virtual bool write
virtual bool writeObjectives
(
const scalar weightedObjective,
const bool valid = true