SubchannelFourEqNonlinearOperator.h

Go to the documentation of this file.

 
#ifndef included_AMP_SubchannelFourEqNonlinearOperator
#define included_AMP_SubchannelFourEqNonlinearOperator
 
#include "AMP/operators/Operator.h"
#include "AMP/operators/subchannel/SubchannelOperatorParameters.h"
 
 
namespace AMP::Operator {
 
class SubchannelFourEqNonlinearOperator : public Operator
{
public:
    typedef std::unique_ptr<AMP::Mesh::MeshElement> ElementPtr;
 
    explicit SubchannelFourEqNonlinearOperator( std::shared_ptr<const OperatorParameters> params );
 
    virtual ~SubchannelFourEqNonlinearOperator() {}
 
    std::string type() const override { return "SubchannelFourEqNonlinearOperator"; }
 
    void apply( AMP::LinearAlgebra::Vector::const_shared_ptr u,
                AMP::LinearAlgebra::Vector::shared_ptr f ) override;
 
    void reset( std::shared_ptr<const OperatorParameters> params ) override;
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getInputVariable() const override
    {
        return d_inpVariable;
    }
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> getOutputVariable() const override
    {
        return d_outVariable;
    }
 
    std::shared_ptr<AMP::LinearAlgebra::VectorSelector> selectOutputVector() const override;
 
    std::shared_ptr<AMP::LinearAlgebra::VectorSelector> selectInputVector() const override;
 
    void setVector( AMP::LinearAlgebra::Vector::shared_ptr frozenVec )
    {
        d_cladTemperature = frozenVec;
    }
 
    std::shared_ptr<SubchannelPhysicsModel> getSubchannelPhysicsModel()
    {
        return d_subchannelPhysicsModel;
    }
 
    double getInletTemperature() { return d_Tin; }
 
    double getOutletPressure() { return d_Pout; }
 
    auto getParams() { return d_params; }
 
    void getLateralFaces( std::shared_ptr<AMP::Mesh::Mesh>,
                          std::map<AMP::Mesh::Point, ElementPtr> &,
                          std::map<AMP::Mesh::Point, ElementPtr> & );
 
    std::map<AMP::Mesh::Point, double> getGapWidths( std::shared_ptr<AMP::Mesh::Mesh>,
                                                     const std::vector<double> &,
                                                     const std::vector<double> &,
                                                     const std::vector<double> & );
 
    void getAxialFaces( const AMP::Mesh::MeshElement &, ElementPtr &, ElementPtr & );
 
    void fillSubchannelGrid( std::shared_ptr<AMP::Mesh::Mesh> ); // function to fill the subchannel
                                                                 // data for all processors
 
    int getSubchannelIndex( double x,
                            double y ); // function to give unique index for each subchannel
 
protected:
    std::shared_ptr<OperatorParameters>
    getJacobianParameters( AMP::LinearAlgebra::Vector::const_shared_ptr u ) override;
 
    std::shared_ptr<const SubchannelOperatorParameters> d_params;
 
    std::shared_ptr<SubchannelPhysicsModel> d_subchannelPhysicsModel;
 
private:
    bool d_initialized;
 
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_inpVariable;
    std::shared_ptr<AMP::LinearAlgebra::Variable> d_outVariable;
    std::shared_ptr<AMP::LinearAlgebra::Vector> d_cladTemperature;
 
    bool d_forceNoConduction; // option to force conduction terms to zero; used for testing
    bool d_forceNoTurbulence; // option to force turbulence terms to zero; used for testing
    bool d_forceNoHeatSource; // option to force heat source terms to zero; used for testing
    bool
        d_forceNoFriction; // option to force friction and form loss terms to zero; used for testing
 
    double d_Pout;           // exit pressure [Pa]
    double d_Tin;            // inlet temperature [K]
    double d_mass;           // inlet global mass flow rate [kg/s]
    double d_win;            // inlet lateral mass flow rate [kg/s]
    double d_gamma;          // fission heating coefficient
    double d_theta;          // channel angle [rad]
    double d_turbulenceCoef; // proportionality constant relating turbulent momentum to turbulent
                             // energy transport
    double d_reynolds;       // reynolds number
    double d_prandtl;        // prandtl number
    double d_KG;             // lateral form loss coefficient
 
    std::string d_frictionModel; // friction model
    double d_friction;           // friction factor
    double d_roughness;          // surface roughness [m]
 
    std::vector<double> d_channelDiam; // Channel hydraulic diameter using the wetted perimeter
    std::vector<double> d_channelArea; // Channel flow area
    std::vector<double> d_rodDiameter; // Average rod diameter for each subchannel
    std::vector<double> d_rodFraction; // Fraction of a rod in each subchannel
    std::vector<double> d_channelMass; // Mass flow rate for each subchannel [kg/s]
 
    size_t d_NGrid;                 // number of grid spacers
    std::vector<double> d_zMinGrid; // z min positions of each grid spacer
    std::vector<double> d_zMaxGrid; // z max positions of each grid spacer
    std::vector<double> d_lossGrid; // loss coefficients for each grid spacer
 
    std::string d_source;            // heat source type
    std::string d_heatShape;         // heat shape used if heat source type is "totalHeatGeneration"
    double d_Q;                      // (sum of rod powers)/4 for each subchannel
    std::vector<double> d_QFraction; // fraction of max rod power in each subchannel
 
    std::vector<double> d_x, d_y, d_z;
    std::vector<bool> d_ownSubChannel; // does this processor own this subchannel (multiple
                                       // processors may own a subchannel)?
    std::vector<std::vector<ElementPtr>> d_subchannelElem; // List of elements in each subchannel
    std::vector<std::vector<ElementPtr>>
        d_subchannelFace;    // List of z-face elements in each subchannel
    size_t d_numSubchannels; // number of subchannels
 
    double Volume( double, double );              // evaluates specific volume
    double Temperature( double, double );         // evaluates temperature
    double ThermalConductivity( double, double ); // evaluates thermal conductivity
    double DynamicViscosity( double, double );    // evaluates dynamic viscosity
    double Enthalpy( double, double );            // evaluates specific enthalpy
 
    ElementPtr getAxiallyAdjacentLateralFace( const AMP::Mesh::MeshElement *,
                                              const AMP::Mesh::MeshElement &,
                                              const std::map<AMP::Mesh::Point, ElementPtr> & );
};
 
 
} // namespace AMP::Operator
 
#endif