adjoint_far_field_lift_response_function.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//
//  Main authors:    Marc Nunez, Inigo Lopez
//
 
#ifndef ADJOINT_LIFT_FAR_FIELD_RESPONSE_FUNCTION_H
#define ADJOINT_LIFT_FAR_FIELD_RESPONSE_FUNCTION_H
 
 
// System includes
 
// External includes
 
// Project includes
#include "adjoint_potential_response_function.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
class KRATOS_API(COMPRESSIBLE_POTENTIAL_FLOW_APPLICATION) AdjointLiftFarFieldResponseFunction : public AdjointPotentialResponseFunction
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(AdjointLiftFarFieldResponseFunction);
 
 
    AdjointLiftFarFieldResponseFunction(ModelPart& rModelPart, Parameters ResponseSettings);
 
    ~AdjointLiftFarFieldResponseFunction() override;
 
 
    void InitializeSolutionStep() override;
 
    double CalculateValue(ModelPart& rModelPart) override;
 
    void CalculateGradient(const Element& rAdjointElement,
                                   const Matrix& rResidualGradient,
                                   Vector& rResponseGradient,
                                   const ProcessInfo& rProcessInfo) override;
 
    void CalculateGradient(const Condition& rAdjointCondition,
                                   const Matrix& rResidualGradient,
                                   Vector& rResponseGradient,
                                   const ProcessInfo& rProcessInfo) override;
 
    void CalculatePartialSensitivity(Element& rAdjointElement,
                                             const Variable<double>& rVariable,
                                             const Matrix& rSensitivityMatrix,
                                             Vector& rSensitivityGradient,
                                             const ProcessInfo& rProcessInfo) override;
 
    void CalculatePartialSensitivity(Condition& rAdjointCondition,
                                             const Variable<double>& rVariable,
                                             const Matrix& rSensitivityMatrix,
                                             Vector& rSensitivityGradient,
                                             const ProcessInfo& rProcessInfo) override;
 
    void CalculatePartialSensitivity(Element& rAdjointElement,
                                             const Variable<array_1d<double, 3>>& rVariable,
                                             const Matrix& rSensitivityMatrix,
                                             Vector& rSensitivityGradient,
                                             const ProcessInfo& rProcessInfo) override;
 
    void CalculatePartialSensitivity(Condition& rAdjointCondition,
                                             const Variable<array_1d<double, 3>>& rVariable,
                                             const Matrix& rSensitivityMatrix,
                                             Vector& rSensitivityGradient,
                                             const ProcessInfo& rProcessInfo) override;
 
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
private:
 
 
 
    double mReferenceChord;
    double mDynamicPressure;
    double mStepSize;
    std::string mFarFieldModelPartName;
    array_1d<double, 3> mFreeStreamVelocity;
    array_1d<double, 3> mWakeNormal;
 
 
 
 
    double ComputeLiftContribution(Element& rElement, const array_1d<double, 3> rNormal, const ProcessInfo& rProcessInfo);
 
 
 
 
 
}; // Class AdjointLiftFarFieldResponseFunction
 
 
 
 
 
} // namespace Kratos.
 
#endif // ADJOINT_LIFT_FAR_FIELD_RESPONSE_FUNCTION_H_INCLUDED