DiffusionTransportModel.h

Go to the documentation of this file.

#ifndef included_AMP_DiffusionTransportModel
#define included_AMP_DiffusionTransportModel
 
#include "AMP/materials/Material.h"
#include "AMP/materials/Property.h"
#include "AMP/operators/ElementPhysicsModel.h"
 
#include <map>
#include <memory>
#include <string>
#include <vector>
 
// Libmesh headers
DISABLE_WARNINGS
#include "libmesh/point.h"
ENABLE_WARNINGS
 
 
namespace AMP::Operator {
 
typedef ElementPhysicsModelParameters DiffusionTransportModelParameters;
 
class DiffusionTransportModel : public ElementPhysicsModel
{
public:
    explicit DiffusionTransportModel(
        std::shared_ptr<const DiffusionTransportModelParameters> params );
 
    virtual ~DiffusionTransportModel() {}
 
    // For LinearOperator's Reset/Init
 
    virtual void preLinearAssembly() {}
 
    virtual void postLinearAssembly() {}
 
    virtual void preLinearElementOperation() {}
 
    virtual void postLinearElementOperation() {}
 
    virtual void preLinearGaussPointOperation() {}
 
    virtual void postLinearGaussPointOperation() {}
 
    // For NonlinearOperator's Init
 
    virtual void preNonlinearInitElementOperation() {}
 
    virtual void postNonlinearInitElementOperation() {}
 
    virtual void nonlinearInitGaussPointOperation() {}
 
    // For NonlinearOperator's Apply
 
    virtual void preNonlinearAssembly() {}
 
    virtual void postNonlinearAssembly() {}
 
    virtual void preNonlinearAssemblyElementOperation() {}
 
    virtual void postNonlinearAssemblyElementOperation() {}
 
    virtual void preNonlinearAssemblyGaussPointOperation() {}
 
    virtual void postNonlinearAssemblyGaussPointOperation() {}
 
    // For NonlinearOperator's Reset
 
    virtual void globalReset() {}
 
    virtual void preNonlinearReset() {}
 
    virtual void postNonlinearReset() {}
 
    virtual void preNonlinearResetElementOperation() {}
 
    virtual void postNonlinearResetElementOperation() {}
 
    virtual void preNonlinearResetGaussPointOperation() {}
 
    virtual void postNonlinearResetGaussPointOperation() {}
 
    virtual void nonlinearResetGaussPointOperation( const double * ) {}
 
    // For NonlinearOperator's GetJacobianParameters
 
    virtual void preNonlinearJacobian() {}
 
    virtual void postNonlinearJacobian() {}
 
    virtual void preNonlinearJacobianElementOperation() {}
 
    virtual void postNonlinearJacobianElementOperation() {}
 
    virtual void preNonlinearJacobianGaussPointOperation() {}
 
    virtual void postNonlinearJacobianGaussPointOperation() {}
 
    virtual void nonlinearJacobianGaussPointOperation( const double * ) {}
 
    static std::shared_ptr<std::vector<double>>
    bilogTransform( const std::vector<double> &u, const double a, const double b );
 
    static void bilogScale( std::vector<double> &u, const double a, const double b );
 
    virtual void getTransport( std::vector<double> &result,
                               std::map<std::string, std::shared_ptr<std::vector<double>>> &args,
                               const std::vector<libMesh::Point> &Coordinates = d_DummyCoords );
 
    std::shared_ptr<AMP::Materials::Material> getMaterial() { return d_material; }
    std::shared_ptr<AMP::Materials::Property> getProperty() { return d_property; }
 
    bool isaTensor() { return d_IsTensor; }
 
protected:
    std::shared_ptr<AMP::Materials::Material> d_material;
 
    std::shared_ptr<AMP::Materials::Property> d_property;
 
    bool d_UseBilogScaling;
 
    std::string d_BilogVariable;
 
    bool d_BilogScaleCoefficient;
 
protected: // used to be private
    std::array<double, 2> d_BilogRange;
    double d_BilogEpsilonRangeLimit;
 
    static const std::vector<libMesh::Point> d_DummyCoords;
 
    bool d_IsTensor;
};
} // namespace AMP::Operator
 
#endif