MixtureThermalExpansion_Def.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  //
//*****************************************************************//

#include "Teuchos_TestForException.hpp"
#include "Phalanx_DataLayout.hpp"

#include "Intrepid_FunctionSpaceTools.hpp"

namespace LCM {

//**********************************************************************
template<typename EvalT, typename Traits>
MixtureThermalExpansion<EvalT, Traits>::
MixtureThermalExpansion(const Teuchos::ParameterList& p) :
  biotCoefficient       (p.get<std::string>     ("Biot Coefficient Name"),
         p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") ),
  porosity       (p.get<std::string>                   ("Porosity Name"),
         p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") ),
  alphaSkeleton  (p.get<std::string>                   ("Skeleton Thermal Expansion Name"),
         p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") ),
  alphaPoreFluid       (p.get<std::string>      ("Pore-Fluid Thermal Expansion Name"),
                 p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") ),
  J           (p.get<std::string>                   ("DetDefGrad Name"),
                 p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") ),
  mixtureThermalExpansion      (p.get<std::string>    ("Mixture Thermal Expansion Name"),
         p.get<Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout") )
{
  this->addDependentField(biotCoefficient);
  this->addDependentField(porosity);
  this->addDependentField(J);
  this->addDependentField(alphaPoreFluid);
  this->addDependentField(alphaSkeleton);

  this->addEvaluatedField(mixtureThermalExpansion);

  this->setName("Mixture Thermal Expansion"+PHX::TypeString<EvalT>::value);

  Teuchos::RCP<PHX::DataLayout> scalar_dl =
    p.get< Teuchos::RCP<PHX::DataLayout> >("QP Scalar Data Layout");
  std::vector<PHX::DataLayout::size_type> dims;
  scalar_dl->dimensions(dims);
  numQPs  = dims[1];
}

//**********************************************************************
template<typename EvalT, typename Traits>
void MixtureThermalExpansion<EvalT, Traits>::
postRegistrationSetup(typename Traits::SetupData d,
                      PHX::FieldManager<Traits>& fm)
{
  this->utils.setFieldData(mixtureThermalExpansion,fm);
  this->utils.setFieldData(biotCoefficient,fm);
  this->utils.setFieldData(porosity,fm);
  this->utils.setFieldData(J,fm);
  this->utils.setFieldData(alphaSkeleton,fm);
  this->utils.setFieldData(alphaPoreFluid,fm);
}

//**********************************************************************
template<typename EvalT, typename Traits>
void MixtureThermalExpansion<EvalT, Traits>::
evaluateFields(typename Traits::EvalData workset)
{
  // Compute Strain tensor from displacement gradient
  for (std::size_t cell=0; cell < workset.numCells; ++cell) {
    for (std::size_t qp=0; qp < numQPs; ++qp) {

        mixtureThermalExpansion(cell,qp) = (biotCoefficient(cell,qp)*J(cell,qp)-porosity(cell,qp))
                                       *alphaSkeleton(cell,qp) +
                                       porosity(cell,qp)*alphaPoreFluid(cell,qp);

    }
  }

}

//**********************************************************************
}