SubchannelHelpers.h

Go to the documentation of this file.

// This file containts some helper functions for subchannel operators
#ifndef included_AMP_SubchannelHelpers
#define included_AMP_SubchannelHelpers
 
#include "AMP/mesh/Mesh.h"
#include "AMP/mesh/MeshElement.h"
#include "AMP/operators/subchannel/SubchannelOperatorParameters.h"
#include "AMP/operators/subchannel/SubchannelPhysicsModel.h"
#include "AMP/vectors/Vector.h"
 
#include <string>
#include <vector>
 
 
namespace AMP::Operator::Subchannel {
 
 
// Convert the input parameters
std::shared_ptr<const SubchannelOperatorParameters>
    convert( std::shared_ptr<const OperatorParameters> );
 
 
size_t getNumberOfSubchannels( std::shared_ptr<AMP::Mesh::Mesh> subchannel );
 
 
std::shared_ptr<AMP::Mesh::Mesh>
subsetForSubchannel( std::shared_ptr<AMP::Mesh::Mesh> subchannel, size_t i, size_t j );
 
 
void getSubchannelProperties( std::shared_ptr<AMP::Mesh::Mesh> subchannel,
                              const std::vector<double> &clad_x,
                              const std::vector<double> &clad_y,
                              const std::vector<double> &clad_d,
                              std::vector<double> &x,
                              std::vector<double> &y,
                              std::vector<double> &area,
                              std::vector<double> &diam,
                              std::vector<double> &rod_diameter,
                              std::vector<double> &channel_fraction );
 
 
void getCladProperties( AMP::AMP_MPI comm,
                        std::shared_ptr<AMP::Mesh::Mesh> clad,
                        std::vector<double> &x,
                        std::vector<double> &y,
                        std::vector<double> &diam );
 
 
std::vector<double>
getHeatFluxGeneration( std::string shape, std::vector<double> z, double diam, double Q_tot );
 
 
std::vector<double> getHeatFluxGenerationWithDiscretizationError( std::string shape,
                                                                  std::vector<double> z,
                                                                  double diam,
                                                                  double Q_tot );
 
 
std::vector<double> getHeatFluxClad( std::vector<double> z,
                                     std::vector<AMP::Mesh::MeshElementID> face_ids,
                                     double channelDiam,
                                     double reynolds,
                                     double prandtl,
                                     double fraction,
                                     std::shared_ptr<SubchannelPhysicsModel> subchannelPhysicsModel,
                                     AMP::LinearAlgebra::Vector::const_shared_ptr flow,
                                     AMP::LinearAlgebra::Vector::const_shared_ptr clad_temp );
 
 
AMP::LinearAlgebra::Vector::shared_ptr
getCladHydraulicDiameter( std::shared_ptr<AMP::Mesh::Mesh> clad,
                          std::shared_ptr<AMP::Mesh::Mesh> subchannel,
                          AMP::AMP_MPI comm );
 
 
} // namespace AMP::Operator::Subchannel
 
#endif