Calculation of adjoint based sensitivities at wall faces. More...

Inheritance diagram for sensitivitySurface:

Collaboration diagram for sensitivitySurface:

[legend]

Public Member Functions
	TypeName ("surface")
	Runtime type information. More...

	sensitivitySurface (const fvMesh &mesh, const dictionary &dict, incompressibleAdjointSolver &adjointSolver)
	Construct from components. More...

virtual	~sensitivitySurface ()=default
	Destructor. More...

void	read ()
	Read controls and update solver pointers if necessary. More...

virtual bool	readDict (const dictionary &dict)
	Read dict if changed. More...

void	computeDerivativesSize ()
	Compute the number of faces on sensitivityPatchIDs_. More...

virtual void	accumulateIntegrand (const scalar dt)
	Accumulate sensitivity integrands. More...

virtual void	assembleSensitivities ()
	Assemble sensitivities. More...

virtual void	clearSensitivities ()
	Zero sensitivity fields and their constituents. More...

autoPtr< adjointEikonalSolver > &	getAdjointEikonalSolver ()
	Get adjoint eikonal solver. More...

virtual void	write (const word &baseName=word::null)
	Write sensitivity maps. More...

bool	getIncludeObjective () const
	Get access to the includeObjective bool. More...

bool	getIncludeSurfaceArea () const
	Get access to the includeSurfaceArea bool. More...

void	setIncludeObjective (const bool includeObjective)
	Set includeObjective bool. More...

void	setIncludeSurfaceArea (const bool includeSurfaceArea)
	Set includeSurfaceArea bool. More...

Public Member Functions inherited from adjointSensitivity
	TypeName ("adjointSensitivity")
	Runtime type information. More...

	declareRunTimeSelectionTable (autoPtr, adjointSensitivity, dictionary,(const fvMesh &mesh, const dictionary &dict, incompressibleAdjointSolver &adjointSolver),(mesh, dict, adjointSolver))

	adjointSensitivity (const fvMesh &mesh, const dictionary &dict, incompressibleAdjointSolver &adjointSolver)
	Construct from components. More...

virtual	~adjointSensitivity ()=default
	Destructor. More...

const incompressibleVars &	primalVars () const
	Get primal variables. More...

const incompressibleAdjointVars &	adjointVars () const
	Get adjoint variables. More...

const incompressibleAdjointSolver &	adjointSolver () const
	Get adjoint solver. More...

virtual const scalarField &	calculateSensitivities ()
	Calculates and returns sensitivity fields. More...

const scalarField &	getSensitivities () const
	Returns the sensitivity fields. More...

tmp< volTensorField >	computeGradDxDbMultiplier ()
	Compute the volTensorField multiplying grad(dxdb) for the volume-based approach to compute shape sensitivity derivatives. More...

tmp< volVectorField >	adjointMeshMovementSource ()
	Compute source term for adjoint mesh movement equation. More...

Public Member Functions inherited from sensitivity
	TypeName ("sensitivity")
	Runtime type information. More...

	sensitivity (const fvMesh &mesh, const dictionary &dict)
	Construct from components. More...

virtual	~sensitivity ()=default
	Destructor. More...

const dictionary &	dict () const
	Return the construction dictionary. More...

Protected Member Functions
void	addGeometricSens ()
	Add sensitivities from dSd/db and dnf/db computed at points and mapped to faces. More...

void	setSuffixName ()
	Set suffix name for sensitivity fields. More...

void	smoothSensitivities ()
	Smooth sensitivity derivatives based on the computation of the 'Sobolev gradient'. More...

scalar	computeRadius (const faMesh &aMesh)
	Compute the physical smoothing radius based on the average boundary face 'length'. More...

Protected Attributes
bool	includeSurfaceArea_
	Include surface area in sens computation. More...

bool	includePressureTerm_
	Include the adjoint pressure term in sens computation. More...

bool	includeGradStressTerm_
	Include the term containing the grad of the stress at the boundary. More...

bool	includeTransposeStresses_
	Include the transpose part of the adjoint stresses. More...

bool	useSnGradInTranposeStresses_
	Use snGrad in the transpose part of the adjoint stresses. More...

bool	includeDivTerm_
	Include the term from the deviatoric part of the stresses. More...

bool	includeDistance_
	Include distance variation in sens computation. More...

bool	includeMeshMovement_
	Include mesh movement variation in sens computation. More...

bool	includeObjective_
	Include terms directly emerging from the objective function. More...

bool	writeGeometricInfo_
	Write geometric info for use by external programs. More...

bool	smoothSensitivities_
	Smooth sensitivity derivatives based on the computation of the 'Sobolev gradient'. More...

autoPtr< adjointEikonalSolver >	eikonalSolver_

autoPtr< adjointMeshMovementSolver >	meshMovementSolver_

autoPtr< volVectorField >	nfOnPatchPtr_

autoPtr< volVectorField >	SfOnPatchPtr_

autoPtr< volVectorField >	CfOnPatchPtr_

Protected Attributes inherited from adjointSensitivity
scalarField	derivatives_

incompressibleAdjointSolver &	adjointSolver_

const incompressibleVars &	primalVars_

incompressibleAdjointVars &	adjointVars_

objectiveManager &	objectiveManager_

Protected Attributes inherited from sensitivity
const fvMesh &	mesh_

dictionary	dict_

autoPtr< volScalarField >	fieldSensPtr_

Additional Inherited Members
Static Public Member Functions inherited from adjointSensitivity
static autoPtr< adjointSensitivity >	New (const fvMesh &mesh, const dictionary &dict, incompressibleAdjointSolver &adjointSolver)
	Return a reference to the selected turbulence model. More...

Detailed Description

Calculation of adjoint based sensitivities at wall faces.

Reference:

    The computation of the sensitivity derivatives  follows the (E-)SI
    formulation of

        Kavvadias, I. S., Papoutsis-Kiachagias, E. M., 
        Giannakoglou, K. C. (2015).
        On the proper treatment of grid sensitivities in continuous adjoint
        methods for shape optimization.
        Journal of computational physics, 301, 1-18.
        https://doi.org/10.1016/j.jcp.2015.08.012

    whereas their smoothing based on the computation of the 'Sobolev
    gradient' is derived from

        Vassberg J. C., Jameson A. (2006).
        Aerodynamic Shape Optimization Part I: Theoretical Background.
        VKI Lecture Series, Introduction to Optimization and
        Multidisciplinary Design, Brussels, Belgium, 8 March, 2006.

Source files

Definition at line 77 of file sensitivitySurfaceIncompressible.H.

Constructor & Destructor Documentation

◆ sensitivitySurface()

sensitivitySurface	(	const fvMesh &	mesh,
		const dictionary &	dict,
		incompressibleAdjointSolver &	adjointSolver
	)

Construct from components.

Definition at line 399 of file sensitivitySurfaceIncompressible.C.

References IOobjectOption::AUTO_WRITE, mesh, IOobjectOption::NO_READ, and Foam::read().

Here is the call graph for this function:

◆ ~sensitivitySurface()

virtual ~sensitivitySurface ( )

virtualdefault

Destructor.

Member Function Documentation

◆ addGeometricSens()

void addGeometricSens ( )

protected

Add sensitivities from dSd/db and dnf/db computed at points and mapped to faces.

Definition at line 55 of file sensitivitySurfaceIncompressible.C.

References fvMesh::boundary(), polyMesh::boundaryMesh(), polyMesh::faces(), forAll, objectiveManager::getObjectiveFunctions(), Foam::I, sensitivitySurface::includeObjective_, deltaBoundary::makeFaceCentresAndAreas_d(), sensitivity::mesh_, primitiveMesh::nPoints(), adjointSensitivity::objectiveManager_, p, Foam::foamVersion::patch, face::points(), polyMesh::points(), UList< T >::size(), syncTools::syncPointList(), and Foam::Zero.

Referenced by sensitivitySurface::accumulateIntegrand().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setSuffixName()

void setSuffixName ( )

protected

Set suffix name for sensitivity fields.

Definition at line 193 of file sensitivitySurfaceIncompressible.C.

References adjointSensitivity::adjointVars_, sensitivity::dict(), sensitivitySurface::includeMeshMovement_, word::null, shapeSensitivitiesBase::setSuffix(), and variablesSet::solverName().

Referenced by sensitivitySurface::write().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ smoothSensitivities()

void smoothSensitivities ( )

protected

Smooth sensitivity derivatives based on the computation of the 'Sobolev gradient'.

Definition at line 214 of file sensitivitySurfaceIncompressible.C.

References dictionary::add(), aMesh(), fvMesh::boundary(), TimePaths::caseSystem(), sensitivitySurface::computeRadius(), polyMesh::dbDir(), DebugInfo, sensitivity::dict(), Foam::dimArea, Foam::dimless, Foam::constant::electromagnetic::e, Foam::endl(), Time::findInstance(), Foam::gMax(), Foam::Info, iters(), Foam::fam::laplacian(), Foam::mag(), sensitivity::mesh_, faMesh::meshSubDir, IOobjectOption::MUST_READ, Foam::nl, IOobjectOption::NO_READ, IOobjectOption::NO_WRITE, patchID, GeometricField< Type, PatchField, GeoMesh >::primitiveFieldRef(), Time::printExecutionTime(), IOobjectOption::READ_IF_PRESENT, autoPtr< T >::ref(), autoPtr< T >::reset(), Foam::fam::Sp(), Foam::sqr(), Foam::sqrt(), dictionary::subDictOrAdd(), fvMesh::time(), Time::timeName(), IOobject::typeHeaderOk(), vsm(), and Foam::Zero.

Referenced by sensitivitySurface::assembleSensitivities().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeRadius()

scalar computeRadius ( const faMesh & aMesh )

protected

Compute the physical smoothing radius based on the average boundary face 'length'.

Definition at line 378 of file sensitivitySurfaceIncompressible.C.

References aMesh(), polyMesh::bounds(), sensitivity::dict(), forAll, Foam::gAverage(), polyMesh::geometricD(), dictionary::getOrDefault(), sensitivity::mesh_, polyMesh::nGeometricD(), Foam::pow(), and boundBox::span().

Referenced by sensitivitySurface::smoothSensitivities().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ TypeName()

TypeName ( "surface" )

Runtime type information.

◆ read()

void read ( )

Read controls and update solver pointers if necessary.

Definition at line 495 of file sensitivitySurfaceIncompressible.C.

Referenced by SIBase::read().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ readDict()

bool readDict ( const dictionary & dict )

virtual

Read dict if changed.

Reimplemented from sensitivity.

Definition at line 555 of file sensitivitySurfaceIncompressible.C.

References sensitivity::dict(), sensitivitySurface::eikonalSolver_, sensitivitySurface::meshMovementSolver_, and sensitivity::readDict().

Referenced by SIBase::readDict().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeDerivativesSize()

void computeDerivativesSize ( )

virtual

Compute the number of faces on sensitivityPatchIDs_.

Reimplemented from sensitivity.

Definition at line 576 of file sensitivitySurfaceIncompressible.C.

References fvMesh::boundary(), adjointSensitivity::derivatives_, sensitivity::mesh_, List< T >::setSize(), and UPtrList< T >::size().

Here is the call graph for this function:

◆ accumulateIntegrand()

void accumulateIntegrand ( const scalar dt )

virtual

Accumulate sensitivity integrands.

Implements adjointSensitivity.

Definition at line 587 of file sensitivitySurfaceIncompressible.C.

Referenced by SIBase::accumulateIntegrand().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ assembleSensitivities()

void assembleSensitivities ( )

virtual

Assemble sensitivities.

Implements adjointSensitivity.

Definition at line 803 of file sensitivitySurfaceIncompressible.C.

Referenced by sensitivityBezier::assembleSensitivities(), and sensitivityVolBSplines::assembleSensitivities().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ clearSensitivities()

void clearSensitivities ( )

virtual

Zero sensitivity fields and their constituents.

Reimplemented from adjointSensitivity.

Definition at line 896 of file sensitivitySurfaceIncompressible.C.

References shapeSensitivitiesBase::clearSensitivities(), adjointSensitivity::clearSensitivities(), sensitivitySurface::eikonalSolver_, sensitivitySurface::includeDistance_, sensitivitySurface::includeMeshMovement_, and sensitivitySurface::meshMovementSolver_.

Referenced by SIBase::clearSensitivities().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getAdjointEikonalSolver()

autoPtr< adjointEikonalSolver > & getAdjointEikonalSolver ( )

Get adjoint eikonal solver.

Definition at line 913 of file sensitivitySurfaceIncompressible.C.

References sensitivitySurface::eikonalSolver_.

◆ write()

void write ( const word & baseName = word::null )

virtual

Write sensitivity maps.

Reimplemented from adjointSensitivity.

Definition at line 919 of file sensitivitySurfaceIncompressible.C.

References sensitivitySurface::CfOnPatchPtr_, sensitivitySurface::nfOnPatchPtr_, sensitivitySurface::setSuffixName(), sensitivitySurface::SfOnPatchPtr_, shapeSensitivitiesBase::write(), adjointSensitivity::write(), and sensitivitySurface::writeGeometricInfo_.

Referenced by SIBase::write().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ getIncludeObjective()

bool getIncludeObjective ( ) const

inline

Get access to the includeObjective bool.

Definition at line 31 of file sensitivitySurfaceIncompressibleI.H.

References sensitivitySurface::includeObjective_.

Referenced by SIBase::read().

Here is the caller graph for this function:

◆ getIncludeSurfaceArea()

bool getIncludeSurfaceArea ( ) const

inline

Get access to the includeSurfaceArea bool.

Definition at line 37 of file sensitivitySurfaceIncompressibleI.H.

References sensitivitySurface::includeSurfaceArea_.

◆ setIncludeObjective()

void setIncludeObjective ( const bool includeObjective )

inline

Set includeObjective bool.

Definition at line 44 of file sensitivitySurfaceIncompressibleI.H.

Referenced by SIBase::read().

Here is the caller graph for this function:

◆ setIncludeSurfaceArea()

void setIncludeSurfaceArea ( const bool includeSurfaceArea )

inline

Set includeSurfaceArea bool.

Definition at line 53 of file sensitivitySurfaceIncompressibleI.H.

Referenced by SIBase::read().

Here is the caller graph for this function:

Member Data Documentation

◆ includeSurfaceArea_

bool includeSurfaceArea_

protected

Include surface area in sens computation.

Definition at line 90 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), sensitivitySurface::getIncludeSurfaceArea(), and sensitivitySurface::read().

◆ includePressureTerm_

bool includePressureTerm_

protected

Include the adjoint pressure term in sens computation.

Definition at line 95 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), and sensitivitySurface::read().

◆ includeGradStressTerm_

bool includeGradStressTerm_

protected

Include the term containing the grad of the stress at the boundary.

Definition at line 100 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), and sensitivitySurface::read().

◆ includeTransposeStresses_

bool includeTransposeStresses_

protected

Include the transpose part of the adjoint stresses.

Definition at line 105 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), and sensitivitySurface::read().

◆ useSnGradInTranposeStresses_

bool useSnGradInTranposeStresses_

protected

Use snGrad in the transpose part of the adjoint stresses.

Definition at line 110 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), and sensitivitySurface::read().

◆ includeDivTerm_

bool includeDivTerm_

protected

Include the term from the deviatoric part of the stresses.

Definition at line 115 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), and sensitivitySurface::read().

◆ includeDistance_

bool includeDistance_

protected

Include distance variation in sens computation.

Definition at line 120 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), sensitivitySurface::assembleSensitivities(), sensitivitySurface::clearSensitivities(), and sensitivitySurface::read().

◆ includeMeshMovement_

bool includeMeshMovement_

protected

Include mesh movement variation in sens computation.

Definition at line 125 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), sensitivitySurface::assembleSensitivities(), sensitivitySurface::clearSensitivities(), sensitivitySurface::read(), and sensitivitySurface::setSuffixName().

◆ includeObjective_

bool includeObjective_

protected

Include terms directly emerging from the objective function.

Definition at line 130 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), sensitivitySurface::addGeometricSens(), sensitivitySurface::getIncludeObjective(), and sensitivitySurface::read().

◆ writeGeometricInfo_

bool writeGeometricInfo_

protected

Write geometric info for use by external programs.

Definition at line 135 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), sensitivitySurface::read(), and sensitivitySurface::write().

◆ smoothSensitivities_

bool smoothSensitivities_

protected

Smooth sensitivity derivatives based on the computation of the 'Sobolev gradient'.

Definition at line 141 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), and sensitivitySurface::read().

◆ eikonalSolver_

autoPtr<adjointEikonalSolver> eikonalSolver_

protected

Definition at line 143 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), sensitivitySurface::assembleSensitivities(), sensitivitySurface::clearSensitivities(), sensitivitySurface::getAdjointEikonalSolver(), sensitivitySurface::read(), and sensitivitySurface::readDict().

◆ meshMovementSolver_

autoPtr<adjointMeshMovementSolver> meshMovementSolver_

protected

Definition at line 145 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::accumulateIntegrand(), sensitivitySurface::assembleSensitivities(), sensitivitySurface::clearSensitivities(), sensitivitySurface::read(), and sensitivitySurface::readDict().

◆ nfOnPatchPtr_

autoPtr<volVectorField> nfOnPatchPtr_

protected

Definition at line 148 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), and sensitivitySurface::write().

◆ SfOnPatchPtr_

autoPtr<volVectorField> SfOnPatchPtr_

protected

Definition at line 149 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), and sensitivitySurface::write().

◆ CfOnPatchPtr_

autoPtr<volVectorField> CfOnPatchPtr_

protected

Definition at line 150 of file sensitivitySurfaceIncompressible.H.

Referenced by sensitivitySurface::assembleSensitivities(), and sensitivitySurface::write().

The documentation for this class was generated from the following files:

src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurface/sensitivitySurfaceIncompressible.H
src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurface/sensitivitySurfaceIncompressible.C
src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurface/sensitivitySurfaceIncompressibleI.H

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ sensitivitySurface()

◆ ~sensitivitySurface()

Member Function Documentation

◆ addGeometricSens()

◆ setSuffixName()

◆ smoothSensitivities()

◆ computeRadius()

◆ TypeName()

◆ read()

◆ readDict()

◆ computeDerivativesSize()

◆ accumulateIntegrand()

◆ assembleSensitivities()

◆ clearSensitivities()

◆ getAdjointEikonalSolver()

◆ write()

◆ getIncludeObjective()

◆ getIncludeSurfaceArea()

◆ setIncludeObjective()

◆ setIncludeSurfaceArea()

Member Data Documentation

◆ includeSurfaceArea_

◆ includePressureTerm_

◆ includeGradStressTerm_

◆ includeTransposeStresses_

◆ useSnGradInTranposeStresses_

◆ includeDivTerm_

◆ includeDistance_

◆ includeMeshMovement_

◆ includeObjective_

◆ writeGeometricInfo_

◆ smoothSensitivities_

◆ eikonalSolver_

◆ meshMovementSolver_

◆ nfOnPatchPtr_

◆ SfOnPatchPtr_

◆ CfOnPatchPtr_