DamageConstitutiveRelation.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

#ifndef _DAMAGECONSTITUTIVERELATION_H_
#define _DAMAGECONSTITUTIVERELATION_H_

#include "ConstitutiveRelation.H"
#include "AMRDamage.H"
#include "NamespaceHeader.H"

class DamageConstitutiveRelation : public ConstitutiveRelation
{
  friend class AMRDamage;

  ConstitutiveRelation* m_undamagedConstitutiveRelation;
  AMRDamage* m_damageModel;

public:
  
  DamageConstitutiveRelation(ConstitutiveRelation* a_ptr, AMRDamage* a_damage)
  {
    m_undamagedConstitutiveRelation = a_ptr->getNewConstitutiveRelation();
    m_damageModel = a_damage;
  }

  virtual ~DamageConstitutiveRelation();
  

  virtual void computeMu(LevelData<FArrayBox>& a_mu,
                         const LevelData<FArrayBox>& a_vel, const Real& a_scale,
                         const LevelData<FArrayBox>* a_crseVel,
                         int a_nRefCrse,
                         const LevelData<FArrayBox>& a_A,
                         const LevelSigmaCS& a_coordSys,
                         const ProblemDomain& a_domain,
                         const IntVect& a_ghostVect = IntVect::Zero) const;
  

  virtual void computeDissipation(LevelData<FArrayBox>& a_dissipation,
                                  const LevelData<FArrayBox>& a_vel, 
                                  const LevelData<FArrayBox>* a_crseVel,
                                  int a_nRefCrse,
                                  const LevelData<FArrayBox>& a_A,
                                  const LevelSigmaCS& a_coordSys,
                                  const ProblemDomain& a_domain,
                                  const IntVect& a_ghostVect = IntVect::Zero) const;

  // computes face-centered \f$\mu_{AS}\f$ based on cell-centered velocity
  virtual void computeFaceMu(LevelData<FluxBox>& a_mu,
                             LevelData<FArrayBox>& a_vel, const Real& a_scale,
                             const LevelData<FArrayBox>* a_crseVel,
                             int a_nRefCrse,
                             const LevelData<FluxBox>& a_A, 
                             const LevelSigmaCS& a_coordSys,
                             const ProblemDomain& a_domain,
                             const IntVect& a_ghostVect = IntVect::Zero) const;


 


  virtual ConstitutiveRelation* getNewConstitutiveRelation() const;

private:

  DamageConstitutiveRelation() : m_undamagedConstitutiveRelation(NULL)
                                      
  {;}



  int level(const ProblemDomain&  a_domain) const
  {
    CH_assert(m_damageModel != NULL); // for now, we will say there *has* to be a damage model
    //FIXME this is a bit klugey... we need to work out access various levels from
    // m_damageModel but we typically don't know which level that we are (no level argument), so for
    //now (maybe we will change the ConstitutiveRelation interface later) 
    //work it out by comparing domain boxes
    int lev = 0;
    while ( a_domain.domainBox() != 
            m_damageModel->damage(lev)->disjointBoxLayout().physDomain().domainBox())
      {
        lev++;
        CH_assert( lev <= m_damageModel->finestLevel());
      }
    return lev;
  }

  static void  modifyMu
  (FArrayBox& a_damage,
   FArrayBox& a_mu,
   const FArrayBox& a_transportedDamage, 
   const FArrayBox& a_thck,
   const Box& a_box) ;

  static void  computeLocalDamage 
  (FArrayBox& a_damage,
   const FArrayBox& a_mu,
   const FArrayBox& a_gradU, 
   const FArrayBox& a_thck,
   const FArrayBox& a_topg,
   const FArrayBox& a_water,
   const Real& a_rhoi,
   const Real& a_rhow,
   const Real& a_gravity,
   const Real& a_sealevel,
   const Box& a_box);

  static void  computeLocalDamageVT
  (FArrayBox& a_damage,
   const FArrayBox& a_vt,
   const FArrayBox& a_thck,
   const FArrayBox& a_topg,
   const FArrayBox& a_water,
   const Real& a_rhoi,
   const Real& a_rhow,
   const Real& a_gravity,
   const Real& a_sealevel,
   const Box& a_box);

  
  // computes face-centered effective viscoity  and a damage coeffcient based on cell-centered velocity
  virtual void computeFaceMuDamage(LevelData<FluxBox>& a_mu,
                                   LevelData<FluxBox>& a_damage,
                                   LevelData<FArrayBox>& a_vel, const Real& a_scale,
                                   const LevelData<FArrayBox>* a_crseVel,
                                   int a_nRefCrse,
                                   const LevelData<FluxBox>& a_A, 
                                   const LevelSigmaCS& a_coordSys,
                                   const ProblemDomain& a_domain,
                                   const IntVect& a_ghostVect = IntVect::Zero) const;


};

#include "NamespaceFooter.H"
#endif