TransportCoefficients.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#if !defined(LCM_Transport_Coefficients_hpp)
#define LCM_Transport_Coefficients_hpp

#include "Phalanx_ConfigDefs.hpp"
#include "Phalanx_Evaluator_WithBaseImpl.hpp"
#include "Phalanx_Evaluator_Derived.hpp"
#include "Phalanx_MDField.hpp"
#include "Albany_Layouts.hpp"

namespace LCM {
  template<typename EvalT, typename Traits>
  class TransportCoefficients : public PHX::EvaluatorWithBaseImpl<Traits>,
                               public PHX::EvaluatorDerived<EvalT, Traits>  {

  public:

    TransportCoefficients(Teuchos::ParameterList& p,
                         const Teuchos::RCP<Albany::Layouts>& dl);

    void postRegistrationSetup(typename Traits::SetupData d,
                               PHX::FieldManager<Traits>& vm);

    void evaluateFields(typename Traits::EvalData d);

  private:

    typedef typename EvalT::ScalarT ScalarT;
    typedef typename EvalT::MeshScalarT MeshScalarT;

    PHX::MDField<ScalarT,Cell,QuadPoint> c_lattice_;

    PHX::MDField<ScalarT,Cell,QuadPoint> n_trap_;

    PHX::MDField<ScalarT,Cell,QuadPoint> eqps_;

    PHX::MDField<ScalarT,Cell,QuadPoint> temperature_;

    PHX::MDField<ScalarT,Cell,QuadPoint,Dim,Dim> F_;

    PHX::MDField<ScalarT,Cell,QuadPoint> J_;

    PHX::MDField<ScalarT,Cell,QuadPoint> k_eq_;

    PHX::MDField<ScalarT,Cell,QuadPoint> diffusion_coefficient_;

    PHX::MDField<ScalarT,Cell,QuadPoint> convection_coefficient_;

    PHX::MDField<ScalarT,Cell,QuadPoint> c_trapped_;

    PHX::MDField<ScalarT,Cell,QuadPoint> eff_diff_;

    PHX::MDField<ScalarT,Cell,QuadPoint> strain_rate_fac_;

    PHX::MDField<ScalarT,Cell,QuadPoint> total_concentration_;

    PHX::MDField<ScalarT,Cell,QuadPoint,Dim,Dim> F_mech_;

    std::size_t num_pts_;

    std::size_t num_dims_;

    std::size_t  worksetSize;

    bool weighted_average_;

    ScalarT alpha_;

    RealType n_lattice_;

    RealType ideal_gas_constant_;

    RealType trap_binding_energy_;

    RealType a_, b_, c_, avogadros_num_;

    RealType pre_exponential_factor_;

    RealType Q_;

    RealType partial_molar_volume_;

    RealType ref_total_concentration_;

    bool lattice_strain_flag_;

    bool have_eqps_;

  };
}

#endif