RobinMatrixCorrection.h

Go to the documentation of this file.

 
#ifndef included_AMP_RobinMatrixCorrection
#define included_AMP_RobinMatrixCorrection
 
#include "AMP/operators/boundary/BoundaryOperator.h"
#include "AMP/operators/boundary/libmesh/NeumannVectorCorrection.h"
#include "AMP/operators/boundary/libmesh/NeumannVectorCorrectionParameters.h"
#include "AMP/operators/boundary/libmesh/RobinMatrixCorrectionParameters.h"
 
#include "AMP/discretization/createLibmeshElements.h"
 
// Libmesh files
#include "libmesh/quadrature.h"
 
#include <string>
 
namespace AMP::Operator {
 
class RobinMatrixCorrection : public BoundaryOperator
{
public:
    explicit RobinMatrixCorrection( std::shared_ptr<const OperatorParameters> params );
 
    void setVariable( const std::shared_ptr<AMP::LinearAlgebra::Variable> &var )
    {
        d_variable = var;
    }
 
    virtual ~RobinMatrixCorrection() {}
 
    std::string type() const override { return "RobinMatrixCorrection"; }
 
    void apply( AMP::LinearAlgebra::Vector::const_shared_ptr,
                AMP::LinearAlgebra::Vector::shared_ptr ) override
    {
        // Do Nothing
    }
 
    void reset( std::shared_ptr<const OperatorParameters> params ) override;
 
    void addRHScorrection( AMP::LinearAlgebra::Vector::shared_ptr rhs ) override;
 
 
protected:
    Discretization::createLibmeshElements libmeshElements;
 
    std::vector<short int> d_boundaryIds;
 
    std::vector<std::vector<unsigned int>> d_dofIds;
 
    std::vector<std::vector<double>> d_robinValues;
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_variable;
 
    double d_hef; // Convective Coefficient
 
    double d_alpha; // pre-factor solid flux
 
    double d_beta;
 
    double d_gamma;
 
    const std::vector<libMesh::Real> *d_JxW;
 
    const std::vector<std::vector<libMesh::Real>> *d_phi;
 
    std::shared_ptr<libMesh::FEType> d_feType;
 
    std::shared_ptr<libMesh::FEBase> d_fe;
 
    std::shared_ptr<libMesh::QBase> d_qrule;
 
    std::string d_qruleOrderName;
 
    libMeshEnums::Order d_feTypeOrder;
    libMeshEnums::FEFamily d_feFamily;
    libMeshEnums::QuadratureType d_qruleType;
    libMeshEnums::Order d_qruleOrder;
 
    AMP::LinearAlgebra::Vector::shared_ptr d_Frozen;
 
    std::shared_ptr<RobinPhysicsModel> d_robinPhysicsModel;
 
    std::shared_ptr<AMP::Discretization::DOFManager> d_dofManager;
 
private:
    std::shared_ptr<NeumannVectorCorrection> d_NeumannCorrection;
    std::shared_ptr<NeumannVectorCorrectionParameters> d_NeumannParams;
};
} // namespace AMP::Operator
 
#endif