PhononModel< T > Class Template Reference

#include <PhononModel.h>

Inheritance diagram for PhononModel< T >:

Collaboration diagram for PhononModel< T >:

Public Types
typedef NumTypeTraits< T > ::T_Scalar	T_Scalar
typedef Vector< T_Scalar, 3 >	VectorT3
typedef Array< VectorT3 >	VectorT3Array
typedef shared_ptr< VectorT3Array >	T3ptr
typedef Kspace< T >	Tkspace
typedef kvol< T >	Tkvol
typedef pmode< T >	Tmode
typedef Array< T >	Tarray
typedef shared_ptr< Tarray >	Tarrptr
typedef map< int, PhononBC< T > * >	PhononBCMap
typedef Tmode::Mode_ptr	Mode_ptr
typedef Tmode::Reflection	Reflection
typedef Tmode::Reflptr	Reflptr
typedef Tmode::Refl_pair	Refl_pair
typedef Tmode::Refl_Map	Refl_Map
typedef Array< int >	BCcellArray
typedef shared_ptr< BCcellArray >	BCellPtr
typedef vector< BCellPtr >	BCcellList
typedef Array< bool >	BCfaceArray
typedef shared_ptr< BCfaceArray >	BfacePtr
typedef vector< BfacePtr >	BCfaceList

Public Member Functions
	PhononModel (const MeshList &meshes, const GeomFields &geomFields, Tkspace &kspace, PhononMacro &macro)
PhononModelOptions< T > &	getOptions ()
PhononBCMap &	getBCs ()
void	init ()
void	callBoundaryConditions ()
void	updateTL ()
void	COMETupdateTL ()
void	updatee0 ()
void	updateHeatFlux ()
void	initPhononModelLinearization (LinearSystem &ls, Tmode &mode)
void	linearizePhononModel (LinearSystem &ls, Tmode &mode)
bool	advance (const int niter)
void	printTemp ()
T	HeatFluxIntegral (const Mesh &mesh, const int faceGroupId)
int	getIters ()
Public Member Functions inherited from Model
	Model (const MeshList &meshes)
virtual	~Model ()
	DEFINE_TYPENAME ("Model")
virtual map< string, shared_ptr< ArrayBase > > &	getPersistenceData ()
virtual void	restart ()

Private Attributes
const GeomFields &	_geomFields
Tkspace &	_kspace
PhononMacro &	_macro
PhononModelOptions< T >	_options
PhononBCMap	_bcMap
MFRPtr	_initialnorm
int	_niters
BCcellList	_BCells
BCfaceList	_BFaces

Additional Inherited Members
Protected Attributes inherited from Model
const MeshList	_meshes
StorageSiteList	_varSites
StorageSiteList	_fluxSites
map< string, shared_ptr < ArrayBase > >	_persistenceData

Detailed Description

template<class T>
class PhononModel< T >

Definition at line 32 of file PhononModel.h.

Member Typedef Documentation

template<class T >

typedef Array<int> PhononModel< T >::BCcellArray

Definition at line 52 of file PhononModel.h.

template<class T >

typedef vector<BCellPtr> PhononModel< T >::BCcellList

Definition at line 54 of file PhononModel.h.

template<class T >

typedef shared_ptr<BCcellArray> PhononModel< T >::BCellPtr

Definition at line 53 of file PhononModel.h.

template<class T >

typedef Array<bool> PhononModel< T >::BCfaceArray

Definition at line 55 of file PhononModel.h.

template<class T >

typedef vector<BfacePtr> PhononModel< T >::BCfaceList

Definition at line 57 of file PhononModel.h.

template<class T >

typedef shared_ptr<BCfaceArray> PhononModel< T >::BfacePtr

Definition at line 56 of file PhononModel.h.

template<class T >

typedef Tmode::Mode_ptr PhononModel< T >::Mode_ptr

Definition at line 47 of file PhononModel.h.

template<class T >

typedef map<int,PhononBC<T>*> PhononModel< T >::PhononBCMap

Definition at line 46 of file PhononModel.h.

template<class T >

typedef Tmode::Refl_Map PhononModel< T >::Refl_Map

Definition at line 51 of file PhononModel.h.

template<class T >

typedef Tmode::Refl_pair PhononModel< T >::Refl_pair

Definition at line 50 of file PhononModel.h.

template<class T >

typedef Tmode::Reflection PhononModel< T >::Reflection

Definition at line 48 of file PhononModel.h.

template<class T >

typedef Tmode::Reflptr PhononModel< T >::Reflptr

Definition at line 49 of file PhononModel.h.

template<class T >

typedef shared_ptr<VectorT3Array> PhononModel< T >::T3ptr

Definition at line 40 of file PhononModel.h.

template<class T >

typedef NumTypeTraits<T>::T_Scalar PhononModel< T >::T_Scalar

Definition at line 37 of file PhononModel.h.

template<class T >

typedef Array<T> PhononModel< T >::Tarray

Definition at line 44 of file PhononModel.h.

template<class T >

typedef shared_ptr<Tarray> PhononModel< T >::Tarrptr

Definition at line 45 of file PhononModel.h.

template<class T >

typedef Kspace<T> PhononModel< T >::Tkspace

Definition at line 41 of file PhononModel.h.

template<class T >

typedef kvol<T> PhononModel< T >::Tkvol

Definition at line 42 of file PhononModel.h.

template<class T >

typedef pmode<T> PhononModel< T >::Tmode

Definition at line 43 of file PhononModel.h.

template<class T >

typedef Vector<T_Scalar,3> PhononModel< T >::VectorT3

Definition at line 38 of file PhononModel.h.

template<class T >

typedef Array<VectorT3> PhononModel< T >::VectorT3Array

Definition at line 39 of file PhononModel.h.

Constructor & Destructor Documentation

template<class T >

PhononModel< T >::PhononModel ( const MeshList &  meshes, const GeomFields &  geomFields, Tkspace &  kspace, PhononMacro &  macro  )

inline

Definition at line 59 of file PhononModel.h.

References PhononModel< T >::_BCells, PhononModel< T >::_bcMap, PhononModel< T >::_BFaces, Model::_meshes, PhononBC< T >::bcType, Mesh::getBoundaryFaceGroups(), Mesh::getCells(), StorageSite::getCount(), Mesh::getFaces(), StorageSite::getSelfCount(), FaceGroup::groupType, and FaceGroup::id.

                                                                                                    :
  
  Model(meshes),
    _geomFields(geomFields),
    _kspace(kspace),
    _macro(macro),
    _BCells(),
    _BFaces()
      { 
      //setting boundary conditions
        const int numMeshes = _meshes.size();
  
      for (int n=0; n<numMeshes; n++) //mesh loop beg
        {
          const Mesh& mesh = *_meshes[n];
          const StorageSite& cells=mesh.getCells();
          const StorageSite& faces=mesh.getFaces();
          const int cellCount=cells.getSelfCount();
          const int faceCount=faces.getCount();

          BCellPtr BCptr(new BCcellArray(cellCount));
          _BCells.push_back(BCptr);
          BCptr->zero();

          BfacePtr BFptr(new BCfaceArray(faceCount));
          BFptr->zero();
          _BFaces.push_back(BFptr);
          
          foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())  //facegroup loop beg
            {
              const FaceGroup& fg = *fgPtr;
              if (_bcMap.find(fg.id) == _bcMap.end())
                {
                  PhononBC<T> *bc(new PhononBC<T>());
                  
                  _bcMap[fg.id] = bc;
                  if((fg.groupType == "wall"))
                    {
                      bc->bcType = "reflecting";
                    }
                  else if (fg.groupType == "velocity-inlet")
                    {
                      bc->bcType = "temperature";
                    }
                  else
                    throw CException("PhononModel: unknown face group type "
                                     + fg.groupType);
                }
            }
        }
    
    }

Member Function Documentation

template<class T >

bool PhononModel< T >::advance ( const int  niter )

inline

Definition at line 530 of file PhononModel.h.

References PhononModel< T >::_initialnorm, PhononModel< T >::_kspace, PhononModel< T >::_niters, PhononModel< T >::_options, PhononModel< T >::callBoundaryConditions(), pmode< T >::getfield(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), LinearSystem::initAssembly(), PhononModel< T >::initPhononModelLinearization(), LinearSystem::initSolve(), PhononModel< T >::linearizePhononModel(), LinearSystem::postSolve(), PhononModel< T >::updatee0(), LinearSystem::updateSolution(), and PhononModel< T >::updateTL().

  {
    bool keepGoing(true);
    for(int n=0; n<niter; n++)  
      {
        const int klength =_kspace.getlength();
        MFRPtr rNorm;
    
        for(int k=0; k<klength;k++)
          {
            
            Tkvol& kv=_kspace.getkvol(k);
            const int mlength=kv.getmodenum();
            
            for(int m=0;m<mlength;m++)
              {
                
                Tmode& mode=kv.getmode(m);
                
                LinearSystem ls;
                initPhononModelLinearization(ls,mode);
                ls.initAssembly();
                linearizePhononModel(ls,mode);
                ls.initSolve();
                
                MFRPtr kNorm(_options.getPhononLinearSolver().solve(ls));
                
                if (!rNorm)
                  rNorm = kNorm;
                else
                  {
                    // find the array for the 0the direction residual and
                    // add the current residual to it
                    Tkvol& kvol0=_kspace.getkvol(0);
                    Tmode& mode0=kvol0.getmode(0);
                    Field& efield0=mode0.getfield();
                    Field& efield=mode.getfield();
                    
                    ArrayBase& rArray0 = (*rNorm)[efield0];
                    ArrayBase& rArrayd = (*kNorm)[efield];
                    rArray0 += rArrayd;
                  }
                
                ls.postSolve();
                ls.updateSolution();
                
                _options.getPhononLinearSolver().cleanup();
                
              }
          }
                
        if (!_initialnorm) _initialnorm = rNorm;

        if (_niters < 5)
           _initialnorm->setMax(*rNorm);
            
 
        MFRPtr normRatio((*rNorm)/(*_initialnorm));     

        if(_niters % _options.showResidual == 0)
          cout << _niters << ": " << *rNorm <<endl;

        _niters++;
        if ((*rNorm < _options.absTolerance)||(*normRatio < _options.relTolerance )) 
          {
            //cout << endl;
            //cout << "Final Residual at "<<_niters<<": "<< *rNorm <<endl;
            keepGoing=false;
            break;}
        
        updateTL();     //update macroparameters
        updatee0();
        callBoundaryConditions();
      }

    return keepGoing;
  }

template<class T >

void PhononModel< T >::callBoundaryConditions ( )

inline

Definition at line 245 of file PhononModel.h.

References PhononModel< T >::_bcMap, PhononModel< T >::_geomFields, PhononModel< T >::_kspace, Model::_meshes, PhononModel< T >::_options, PhononInterface< X >::applyInterfaceCondition(), PhononBoundary< X >::applyReflectingWall(), PhononBoundary< X >::applyTemperatureWall(), PhononBC< T >::bcType, Mesh::getBoundaryFaceGroups(), Mesh::getInterfaceGroups(), FloatVarDict< T >::getVal(), FaceGroup::id, and FaceGroup::site.

Referenced by PhononModel< T >::advance().

  { 
    const int numMeshes = _meshes.size();
    for (int n=0; n<numMeshes; n++)
      {
        const Mesh& mesh = *_meshes[n];
        
        foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
          {
            const FaceGroup& fg = *fgPtr;
            const StorageSite& faces = fg.site;
            const PhononBC<T>& bc = *_bcMap[fg.id];
            
            PhononBoundary<T> pbc(faces, mesh,_geomFields,_kspace,_options,fg.id);
            
            if (bc.bcType == "reflecting")
              { 
     
                FloatValEvaluator<T>
                  bReflection(bc.getVal("specifiedReflection"),faces);

                pbc.applyReflectingWall(bReflection);
              }
            else if(bc.bcType=="temperature")
              {
                
                FloatValEvaluator<T>
                  bTemperature(bc.getVal("specifiedTemperature"),faces);
          
                pbc.applyTemperatureWall(bTemperature);
              } 
            else
              {
                cout<<"Couldn't find boundary condition"<<endl;
                break;
              }
          };
        foreach(const FaceGroupPtr fgPtr, mesh.getInterfaceGroups())
          {
            int otherMeshID = 0;
            T_Scalar btrans = _options["transmissivity0to1"];
            T_Scalar brefl = 1.0 - _options["transmissivity1to0"];
            if (n==0)
              {
                otherMeshID = 1;
                btrans = _options["transmissivity1to0"];
                brefl = 1.0 - _options["transmissivity0to1"];
              }
            const Mesh& otherMesh = *_meshes[otherMeshID];
            const FaceGroup& fg = *fgPtr;
            PhononInterface<T> pInt(fg, mesh, otherMesh,_geomFields,_kspace,_kspace,_options);

            pInt.applyInterfaceCondition(brefl,btrans);
          };
      };
  };

template<class T >

void PhononModel< T >::COMETupdateTL ( )

inline

Definition at line 341 of file PhononModel.h.

References PhononModel< T >::_kspace, PhononModel< T >::_macro, Model::_meshes, PhononMacro::e0, Mesh::getCells(), StorageSite::getCount(), Kspace< T >::NewtonSolve(), and PhononMacro::temperature.

  {
    const int numMeshes = _meshes.size();
    for (int n=0; n<numMeshes; n++)
      {
        const Mesh& mesh = *_meshes[n];
        const StorageSite& cells = mesh.getCells();
        const int numcells = cells.getCount();
        Tarray& TL=dynamic_cast<Tarray&>(_macro.temperature[cells]);
        Tarray& e0Array=dynamic_cast<Tarray&>(_macro.e0[cells]);
        
        for(int c=0;c<numcells;c++)
          _kspace.NewtonSolve(TL[c],e0Array[c]);
      }
  }

template<class T >

PhononBCMap& PhononModel< T >::getBCs ( )

inline

Definition at line 113 of file PhononModel.h.

References PhononModel< T >::_bcMap.

{return _bcMap;}

template<class T >

int PhononModel< T >::getIters ( )

inline

Definition at line 668 of file PhononModel.h.

References PhononModel< T >::_niters.

{return _niters;}

template<class T >

PhononModelOptions<T>& PhononModel< T >::getOptions ( )

inline

Definition at line 112 of file PhononModel.h.

References PhononModel< T >::_options.

{return _options;}

template<class T >

T PhononModel< T >::HeatFluxIntegral ( const Mesh &  mesh, const int  faceGroupId  )

inline

Definition at line 623 of file PhononModel.h.

References PhononModel< T >::_geomFields, PhononModel< T >::_kspace, GeomFields::area, Mesh::getBoundaryFaceGroups(), Mesh::getCells(), StorageSite::getCount(), kvol< T >::getdk3(), Kspace< T >::getDK3(), Mesh::getFaceCells(), pmode< T >::getfield(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), pmode< T >::getv(), FaceGroup::id, and FaceGroup::site.

  {
    T r(0.);
    bool found = false;
    const T DK3=_kspace.getDK3();
    foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
    {
        const FaceGroup& fg = *fgPtr;
        if (fg.id == faceGroupId)
        {
          const StorageSite& faces = fg.site;
          const int nFaces = faces.getCount();
          const StorageSite& cells = mesh.getCells();
          const CRConnectivity& faceCells=mesh.getFaceCells(faces);
          const Field& areaField=_geomFields.area;
          const VectorT3Array& faceArea=dynamic_cast<const VectorT3Array&>(areaField[faces]);
         
            
          for(int k=0;k<_kspace.getlength();k++)
            {
              Tkvol& kv=_kspace.getkvol(k);
              int modenum=kv.getmodenum();
              for(int m=0;m<modenum;m++)
                {
                  VectorT3 vg=kv.getmode(m).getv();
                  T dk3=kv.getdk3();
                  Field& efield=kv.getmode(m).getfield();
                  const Tarray& eval=dynamic_cast<const Tarray&>(efield[cells]);
                  for(int f=0; f<nFaces; f++)
                    {
                      const VectorT3 An=faceArea[f];
                      const int c1=faceCells(f,1);
                      const T vgdotAn=An[0]*vg[0]+An[1]*vg[1]+An[2]*vg[2];
                      r += eval[c1]*vgdotAn*dk3/DK3;
                    }
                  found=true;
                }
            }
        }
    }
    if (!found)
      throw CException("getHeatFluxIntegral: invalid faceGroupID");
    return r*DK3;
  }

template<class T >

void PhononModel< T >::init ( )

inlinevirtual

Implements Model.

Definition at line 115 of file PhononModel.h.

References PhononModel< T >::_BCells, PhononModel< T >::_bcMap, PhononModel< T >::_BFaces, PhononModel< T >::_geomFields, PhononModel< T >::_initialnorm, PhononModel< T >::_kspace, PhononModel< T >::_macro, Model::_meshes, PhononModel< T >::_niters, PhononModel< T >::_options, Field::addArray(), GeomFields::area, GeomFields::areaMag, pmode< T >::calce0(), PhononMacro::e0, Kspace< T >::findSpecs(), Mesh::getBoundaryFaceGroups(), Mesh::getCells(), StorageSite::getCount(), kvol< T >::getdk3(), Kspace< T >::getDK3(), pmode< T >::gete0field(), Mesh::getFaceCells(), pmode< T >::getfield(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), StorageSite::getOffset(), pmode< T >::getreflmap(), pmode< T >::getresid(), pmode< T >::getv(), FaceGroup::id, FaceGroup::site, sqrt(), and PhononMacro::temperature.

  {
    
    _initialnorm = MFRPtr();
    _niters=0;
    const int numMeshes=_meshes.size();
    
    for (int n=0;n<numMeshes;n++)  //mesh loop beg
      {
        const Mesh& mesh=*_meshes[n];
        const int numK=_kspace.getlength();
        const StorageSite& cells=mesh.getCells();
        const int numcells=cells.getCount();
        
        //initialize lattice temperature
        shared_ptr<Tarray> TLcell(new Tarray(numcells));
        const T Tinit=_options["initialTemperature"];
        *TLcell=Tinit;
        _macro.temperature.addArray(cells,TLcell);

        T e0sum=0.;
        
        for (int k=0;k<numK;k++)  //kspace loop beg
          {
            Tkvol& kv=_kspace.getkvol(k);
            const int numM=kv.getmodenum();
            const T dk3=kv.getdk3();
            
            for (int m=0;m<numM;m++) //mode loop beg
              {
                //initialize each phonon mode
                Tmode& mode=kv.getmode(m);
                Field& efield=mode.getfield();
                Field& e0field=mode.gete0field();
                Field& resfield=mode.getresid();
                Tarrptr e0var=shared_ptr<Tarray>(new Tarray(numcells));
                Tarrptr evar=shared_ptr<Tarray>(new Tarray(numcells));
                Tarrptr resid=shared_ptr<Tarray>(new Tarray(numcells));
                const T einit=mode.calce0(Tinit);
                e0sum+=einit*dk3;
                *evar=einit;
                *e0var=einit;
                *resid=0.;
                efield.addArray(cells,evar);
                e0field.addArray(cells,e0var);
                resfield.addArray(cells,resid);
              }; //mode loop end
          }; //kspace loop end

        e0sum=e0sum/_kspace.getDK3();
        shared_ptr<Tarray> e0cell(new Tarray(numcells));
        *e0cell=e0sum;
        _macro.e0.addArray(cells,e0cell);
        shared_ptr<Tarray> e0ResidCell(new Tarray(numcells));
        *e0ResidCell=0.;
        //      _macro.e0Residual.addArray(cells,e0ResidCell);

        // Compute reflections--later, should be moved inside above loop
        
        for (int k=0;k<numK;k++)  //kspace loop beg
          {
            Tkvol& kv=_kspace.getkvol(k);
            const int numM=kv.getmodenum();
            //const T dk3=kv.getdk3();
            
            for (int m=0;m<numM;m++) //mode loop beg
              { 
                Tmode& mode=kv.getmode(m);
                Refl_Map& rmap=mode.getreflmap();
                VectorT3 vg=mode.getv();
                T vmag=sqrt(pow(vg[0],2)+pow(vg[1],2)+pow(vg[2],2));
                VectorT3 si=vg/vmag;
                VectorT3 so;

                foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
                  {
                    const FaceGroup& fg = *fgPtr;
                    const StorageSite& faces = fg.site;
                    const Field& AreaMagField=_geomFields.areaMag;
                    const Tarray& AreaMag=dynamic_cast<const Tarray&>(AreaMagField[faces]);
                    const Field& AreaDirField=_geomFields.area;
                    const VectorT3Array& AreaDir=
                      dynamic_cast<const VectorT3Array&>(AreaDirField[faces]);
                              
                    const VectorT3 n=AreaDir[0]/AreaMag[0];
                    const T sidotn=si[0]*n[0]+si[1]*n[1]+si[2]*n[2];

                    if (sidotn > T_Scalar(0.0))
                      {
                        so=si-2.*(si[0]*n[0]+si[1]*n[1]+si[2]*n[2])*n;
                        T soMag=sqrt(pow(so[0],2)+pow(so[1],2)+pow(so[2],2));
                        so/=soMag;
                        so*=vmag;
                        Refl_pair refls;
                        Refl_pair reflsFrom;
                        _kspace.findSpecs(n,m,k,refls);
                        rmap[fg.id]=refls;
                        const int k1=refls.first.second;
                        Tmode& mode2=_kspace.getkvol(k1).getmode(m);
                        Refl_Map& rmap2=mode2.getreflmap();
                        reflsFrom.first.second=-1;
                        reflsFrom.second.second=k;
                        rmap2[fg.id]=reflsFrom;
                      }
                  }
              }
          }

        BCcellArray& BCArray=*(_BCells[n]);
        BCfaceArray& BCfArray=*(_BFaces[n]);
        foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
          {
            const FaceGroup& fg = *fgPtr;
            if((_bcMap[fg.id]->bcType == "reflecting"))
              {
                const CRConnectivity& BfaceCells=mesh.getFaceCells(fg.site);
                const int faceCount=fg.site.getCount();
                const int offSet=fg.site.getOffset();
                for(int i=0;i<faceCount;i++)
                  {
                    int cell1=BfaceCells(i,0);
                    BCArray[cell1]=1;
                  }
                for(int i=offSet;i<offSet+faceCount;i++)
                  BCfArray[i]=true;
              } 
          }     
      }
  }

template<class T >

void PhononModel< T >::initPhononModelLinearization ( LinearSystem &  ls, Tmode &  mode  )

inline

Definition at line 422 of file PhononModel.h.

References Model::_meshes, MultiField::addArray(), MultiFieldMatrix::addMatrix(), Field::getArrayPtr(), Mesh::getCellCells(), Mesh::getCells(), pmode< T >::getfield(), LinearSystem::getMatrix(), and LinearSystem::getX().

Referenced by PhononModel< T >::advance().

 {
   const int numMeshes = _meshes.size();
   for (int n=0; n<numMeshes; n++)
     {
       const Mesh& mesh = *_meshes[n];
       const StorageSite& cells = mesh.getCells();
       
       Field& fnd = mode.getfield();
       
       MultiField::ArrayIndex vIndex(&fnd,&cells); 
       
       ls.getX().addArray(vIndex,fnd.getArrayPtr(cells));
       
       const CRConnectivity& cellCells = mesh.getCellCells();
       
       shared_ptr<Matrix> m(new CRMatrix<T,T,T>(cellCells));
       
       ls.getMatrix().addMatrix(vIndex,vIndex,m);
     }
   
  }; 

template<class T >

void PhononModel< T >::linearizePhononModel ( LinearSystem &  ls, Tmode &  mode  )

inline

Definition at line 445 of file PhononModel.h.

References PhononModel< T >::_bcMap, PhononModel< T >::_geomFields, Model::_meshes, GeomFields::area, GeomFields::areaMag, PhononBC< T >::bcType, LinearSystem::getB(), Mesh::getBoundaryFaceGroups(), Mesh::getCells(), StorageSite::getCount(), pmode< T >::gete0field(), Mesh::getFaceCells(), pmode< T >::getfield(), LinearSystem::getMatrix(), pmode< T >::gettau(), pmode< T >::getv(), LinearSystem::getX(), FaceGroup::id, Linearizer::linearize(), and FaceGroup::site.

Referenced by PhononModel< T >::advance().

   {
     // _velocityGradientModel.compute();
     DiscrList discretizations;
    
     Field& fnd = mode.getfield();  //field for e"
     Field& e0fld=mode.gete0field(); //field for e0
     T tau=mode.gettau();
     VectorT3 vg=mode.getv();
    
     shared_ptr<Discretization>
       colldisc(new PhononCollisionDiscretization<T,T,T>(_meshes, _geomFields,fnd,e0fld,tau)); 
     
     discretizations.push_back(colldisc);
     
     shared_ptr<Discretization> 
       convdisc(new PhononConvectionDiscretization<T,T,T> (_meshes,_geomFields,fnd,vg));

     discretizations.push_back(convdisc);
     
     Linearizer linearizer;
     
     linearizer.linearize(discretizations,_meshes,ls.getMatrix(),
                          ls.getX(), ls.getB());
     
     

     // boundary conditions
     const int numMeshes = _meshes.size();
     for (int n=0; n<numMeshes; n++)
       {
         const Mesh& mesh = *_meshes[n];
         const StorageSite& cells = mesh.getCells();
         foreach(const FaceGroupPtr fgPtr, mesh.getBoundaryFaceGroups())
           {
             const FaceGroup& fg = *fgPtr;
             const StorageSite& faces = fg.site;
             const int nFaces = faces.getCount();
             const PhononBC<T>& bc = *_bcMap[fg.id];
             Tarray& dsf = dynamic_cast< Tarray&>(fnd[cells]);
             BaseGenericPhononBCS<T,T,T> gpbc(faces, mesh, _geomFields,
                                               fnd,
                                               ls.getMatrix(),
                                               ls.getX(),
                                               ls.getB());
             
             
             const CRConnectivity& faceCells = mesh.getFaceCells(faces);                         
             const Field& areaMagField = _geomFields.areaMag;
             const Tarray& faceAreaMag = dynamic_cast<const Tarray &>(areaMagField[faces]);
             const Field& areaField = _geomFields.area;
             const VectorT3Array& faceArea=dynamic_cast<const VectorT3Array&>(areaField[faces]); 
             
             if (( bc.bcType == "CopyBC")) 
               {
                 for(int f=0; f< nFaces; f++)
                   {
                     gpbc.applyExtrapolationBC(f);
                   }
               }     
             else{
               for(int f=0; f< nFaces; f++)
                 {
                   const VectorT3 en = faceArea[f]/faceAreaMag[f];
                   const T v_dot_en = vg[0]*en[0]+vg[1]*en[1]+vg[2]*en[2];
                   
                   if(v_dot_en < T_Scalar(0.0))
                     {               
                       //incoming direction - dirchlet bc
                       const int c1= faceCells(f,1);    // boundary cell
                       T bvalue =dsf[c1];
                       gpbc.applyDirichletBC(f,bvalue);
                     }
                   else
                     {
                     //outgoing direction - extrapolation bc
                     gpbc.applyExtrapolationBC(f);
                     }    
                 }
             }   
           }
       }

   }

template<class T >

void PhononModel< T >::printTemp ( )

inline

Definition at line 608 of file PhononModel.h.

References PhononModel< T >::_macro, Model::_meshes, Mesh::getCells(), StorageSite::getCount(), and PhononMacro::temperature.

  {
     const int numMeshes = _meshes.size();
     for (int n=0; n<numMeshes; n++)
       {
         const Mesh& mesh = *_meshes[n];
         const StorageSite& cells = mesh.getCells();
         const int numcells=cells.getCount();
         Tarray& TL=dynamic_cast<Tarray&>(_macro.temperature[cells]);
         for(int i=0;i<numcells;i++)
           cout<<TL[i]<<" "<<i<<endl;
       }

  }

template<class T >

void PhononModel< T >::updatee0 ( )

inline

Definition at line 357 of file PhononModel.h.

References PhononModel< T >::_kspace, PhononModel< T >::_macro, Model::_meshes, pmode< T >::calce0(), Mesh::getCells(), StorageSite::getCount(), pmode< T >::gete0field(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), and PhononMacro::temperature.

Referenced by PhononModel< T >::advance().

  {
    
    const int numMeshes = _meshes.size();
    for (int n=0; n<numMeshes; n++)
      {
        const Mesh& mesh = *_meshes[n];
        const StorageSite& cells = mesh.getCells();
        const int numcells = cells.getCount();
        const int numK=_kspace.getlength();
        Tarray& TL=dynamic_cast<Tarray&>(_macro.temperature[cells]);
        
        for(int k=0;k<numK;k++)
          {
            Tkvol& kv=_kspace.getkvol(k);
            const int modenum=kv.getmodenum();
            
            for(int m=0;m<modenum;m++)
              {         
                Tmode& mode=kv.getmode(m);
                Field& e0field=mode.gete0field();
                Tarray& e0_val=dynamic_cast<Tarray&>(e0field[cells]);
                for(int c=0;c<numcells;c++)
                  e0_val[c]=mode.calce0(TL[c]);
              }
          }
        
      }
  }

template<class T >

void PhononModel< T >::updateHeatFlux ( )

inline

Definition at line 387 of file PhononModel.h.

References PhononModel< T >::_kspace, PhononModel< T >::_macro, Model::_meshes, Field::addArray(), Mesh::getCells(), StorageSite::getCount(), kvol< T >::getdk3(), pmode< T >::getfield(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), pmode< T >::getv(), and PhononMacro::heatFlux.

 {
   const int numMeshes = _meshes.size();
   for (int n=0; n<numMeshes; n++)
     {
       VectorT3 zero_vec;
       const Mesh& mesh = *_meshes[n];
       const StorageSite& cells = mesh.getCells();
       const int numcells = cells.getCount();
       T3ptr heatFluxptr=T3ptr(new VectorT3Array(numcells));
       VectorT3Array& heatFlux=*heatFluxptr;
       zero_vec[0]=0.;zero_vec[1]=0.;zero_vec[2]=0.;
       heatFlux=zero_vec;
       
       const int numK=_kspace.getlength();
       for(int k=0;k<numK;k++)
         {
           Tkvol& kv=_kspace.getkvol(k);
           T dk3=kv.getdk3();
           const int modenum=kv.getmodenum();
           for(int m=0;m<modenum;m++)
             {
               Tmode& mode=kv.getmode(m);
               VectorT3 vg=mode.getv();
               Field& efield=mode.getfield();
               const Tarray& eval=dynamic_cast<const Tarray&>(efield[cells]);

               for(int c=0;c<numcells;c++)
                 heatFlux[c]+=eval[c]*vg*dk3;
             }
         }
       _macro.heatFlux.addArray(cells,heatFluxptr);  
     }
 }

template<class T >

void PhononModel< T >::updateTL ( )

inline

Definition at line 302 of file PhononModel.h.

References PhononModel< T >::_kspace, PhononModel< T >::_macro, Model::_meshes, Mesh::getCells(), StorageSite::getCount(), kvol< T >::getdk3(), pmode< T >::getfield(), Kspace< T >::getkvol(), Kspace< T >::getlength(), kvol< T >::getmode(), kvol< T >::getmodenum(), pmode< T >::gettau(), Kspace< T >::NewtonSolve(), and PhononMacro::temperature.

Referenced by PhononModel< T >::advance().

  {

    const int numMeshes = _meshes.size();
    for (int n=0; n<numMeshes; n++)
      {
        const Mesh& mesh = *_meshes[n];
        const StorageSite& cells = mesh.getCells();
        const int numcells = cells.getCount();
        const int numK=_kspace.getlength();
        Tarray& TL=dynamic_cast<Tarray&>(_macro.temperature[cells]);
        Tarrptr e_sumptr=shared_ptr<Tarray>(new Tarray(numcells));
        Tarray& e_sum=*(e_sumptr);
        e_sum=0.;
        
        for(int k=0;k<numK;k++)
          {
            Tkvol& kv=_kspace.getkvol(k);
            const int modenum=kv.getmodenum();
            T dk3=kv.getdk3();
            
            for(int m=0;m<modenum;m++)
              {
                Tmode& mode=kv.getmode(m);
                T tau=mode.gettau();
                Field& efield=mode.getfield();
                Tarray& e_val=dynamic_cast<Tarray&>(efield[cells]);
                
                for(int c=0;c<numcells;c++)
                    e_sum[c]+=e_val[c]*dk3/tau;

              }
          }
        
        for(int c=0;c<numcells;c++)
          _kspace.NewtonSolve(TL[c],e_sum[c]);
      }
  }

Member Data Documentation

template<class T >

BCcellList PhononModel< T >::_BCells

private

Definition at line 679 of file PhononModel.h.

Referenced by PhononModel< T >::init(), and PhononModel< T >::PhononModel().

template<class T >

PhononBCMap PhononModel< T >::_bcMap

private

Definition at line 676 of file PhononModel.h.

Referenced by PhononModel< T >::callBoundaryConditions(), PhononModel< T >::getBCs(), PhononModel< T >::init(), PhononModel< T >::linearizePhononModel(), and PhononModel< T >::PhononModel().

template<class T >

BCfaceList PhononModel< T >::_BFaces

private

Definition at line 680 of file PhononModel.h.

Referenced by PhononModel< T >::init(), and PhononModel< T >::PhononModel().

template<class T >

const GeomFields& PhononModel< T >::_geomFields

private

Definition at line 672 of file PhononModel.h.

Referenced by PhononModel< T >::callBoundaryConditions(), PhononModel< T >::HeatFluxIntegral(), PhononModel< T >::init(), and PhononModel< T >::linearizePhononModel().

template<class T >

MFRPtr PhononModel< T >::_initialnorm

private

Definition at line 677 of file PhononModel.h.

Referenced by PhononModel< T >::advance(), and PhononModel< T >::init().

template<class T >

Tkspace& PhononModel< T >::_kspace

private

Definition at line 673 of file PhononModel.h.

Referenced by PhononModel< T >::advance(), PhononModel< T >::callBoundaryConditions(), PhononModel< T >::COMETupdateTL(), PhononModel< T >::HeatFluxIntegral(), PhononModel< T >::init(), PhononModel< T >::updatee0(), PhononModel< T >::updateHeatFlux(), and PhononModel< T >::updateTL().

template<class T >

PhononMacro& PhononModel< T >::_macro

private

Definition at line 674 of file PhononModel.h.

Referenced by PhononModel< T >::COMETupdateTL(), PhononModel< T >::init(), PhononModel< T >::printTemp(), PhononModel< T >::updatee0(), PhononModel< T >::updateHeatFlux(), and PhononModel< T >::updateTL().

template<class T >

int PhononModel< T >::_niters

private

Definition at line 678 of file PhononModel.h.

Referenced by PhononModel< T >::advance(), PhononModel< T >::getIters(), and PhononModel< T >::init().

template<class T >

PhononModelOptions<T> PhononModel< T >::_options

private

Definition at line 675 of file PhononModel.h.

Referenced by PhononModel< T >::advance(), PhononModel< T >::callBoundaryConditions(), PhononModel< T >::getOptions(), and PhononModel< T >::init().

---

The documentation for this class was generated from the following file:

• PhononModel.h