ElecDiffusionDiscretization< X, Diag, OffDiag > Class Template Reference

#include <ElecDiffusionDiscretization.h>

Inheritance diagram for ElecDiffusionDiscretization< X, Diag, OffDiag >:

Collaboration diagram for ElecDiffusionDiscretization< X, Diag, OffDiag >:

Public Types
typedef NumTypeTraits< X > ::T_Scalar	T_Scalar
typedef CRMatrix< Diag, OffDiag, X >	CCMatrix
typedef CCMatrix::DiagArray	DiagArray
typedef CCMatrix::PairWiseAssembler	CCAssembler
typedef Gradient< X >	XGrad
typedef Array< int >	IntArray
typedef Array< X >	XArray
typedef Array< T_Scalar >	TArray
typedef Vector< T_Scalar, 3 >	VectorT3
typedef Array< VectorT3 >	VectorT3Array
typedef Array< XGrad >	GradArray

Public Member Functions
	ElecDiffusionDiscretization (const MeshList &meshes, const GeomFields &geomFields, Field &varField, const Field &diffusivityField, const Field &varGradientField)
void	discretize (const Mesh &mesh, MultiFieldMatrix &mfmatrix, MultiField &xField, MultiField &rField)
Public Member Functions inherited from Discretization
	Discretization (const MeshList &meshes)
virtual	~Discretization ()
	DEFINE_TYPENAME ("Discretization")

Private Attributes
const GeomFields &	_geomFields
Field &	_varField
const Field &	_diffusivityField
const Field &	_varGradientField

Additional Inherited Members
Protected Attributes inherited from Discretization
const MeshList &	_meshes

Detailed Description

template<class X, class Diag, class OffDiag>
class ElecDiffusionDiscretization< X, Diag, OffDiag >

Definition at line 31 of file ElecDiffusionDiscretization.h.

Member Typedef Documentation

template<class X , class Diag , class OffDiag >

typedef CCMatrix::PairWiseAssembler ElecDiffusionDiscretization< X, Diag, OffDiag >::CCAssembler

Definition at line 39 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef CRMatrix<Diag,OffDiag,X> ElecDiffusionDiscretization< X, Diag, OffDiag >::CCMatrix

Definition at line 37 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef CCMatrix::DiagArray ElecDiffusionDiscretization< X, Diag, OffDiag >::DiagArray

Definition at line 38 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<XGrad> ElecDiffusionDiscretization< X, Diag, OffDiag >::GradArray

Definition at line 49 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<int> ElecDiffusionDiscretization< X, Diag, OffDiag >::IntArray

Definition at line 43 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef NumTypeTraits<X>::T_Scalar ElecDiffusionDiscretization< X, Diag, OffDiag >::T_Scalar

Definition at line 35 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<T_Scalar> ElecDiffusionDiscretization< X, Diag, OffDiag >::TArray

Definition at line 45 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Vector<T_Scalar,3> ElecDiffusionDiscretization< X, Diag, OffDiag >::VectorT3

Definition at line 46 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<VectorT3> ElecDiffusionDiscretization< X, Diag, OffDiag >::VectorT3Array

Definition at line 47 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<X> ElecDiffusionDiscretization< X, Diag, OffDiag >::XArray

Definition at line 44 of file ElecDiffusionDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Gradient<X> ElecDiffusionDiscretization< X, Diag, OffDiag >::XGrad

Definition at line 41 of file ElecDiffusionDiscretization.h.

Constructor & Destructor Documentation

template<class X , class Diag , class OffDiag >

ElecDiffusionDiscretization< X, Diag, OffDiag >::ElecDiffusionDiscretization ( const MeshList &  meshes, const GeomFields &  geomFields, Field &  varField, const Field &  diffusivityField, const Field &  varGradientField  )

inline

Definition at line 51 of file ElecDiffusionDiscretization.h.

                                                         :
    Discretization(meshes),
    _geomFields(geomFields),
    _varField(varField),
    _diffusivityField(diffusivityField),
    _varGradientField(varGradientField)
  {}

Member Function Documentation

template<class X , class Diag , class OffDiag >

void ElecDiffusionDiscretization< X, Diag, OffDiag >::discretize ( const Mesh &  mesh, MultiFieldMatrix &  mfmatrix, MultiField &  xField, MultiField &  rField  )

inlinevirtual

Implements Discretization.

Definition at line 63 of file ElecDiffusionDiscretization.h.

References ElecDiffusionDiscretization< X, Diag, OffDiag >::_diffusivityField, ElecDiffusionDiscretization< X, Diag, OffDiag >::_geomFields, ElecDiffusionDiscretization< X, Diag, OffDiag >::_varField, ElecDiffusionDiscretization< X, Diag, OffDiag >::_varGradientField, GeomFields::area, GeomFields::areaMag, GeomFields::coordinate, dot(), Mesh::getAllFaceCells(), Mesh::getCells(), CRMatrix< T_Diag, T_OffDiag, X >::PairWiseAssembler::getCoeff01(), CRMatrix< T_Diag, T_OffDiag, X >::PairWiseAssembler::getCoeff10(), StorageSite::getCount(), CRMatrix< T_Diag, T_OffDiag, X >::getDiag(), Mesh::getFaces(), MultiFieldMatrix::getMatrix(), CRMatrix< T_Diag, T_OffDiag, X >::getPairWiseAssembler(), harmonicAverage(), GeomFields::ibType, Mesh::IBTYPE_BOUNDARY, Mesh::IBTYPE_FLUID, and GeomFields::volume.

  {
    const StorageSite& cells = mesh.getCells();
    const StorageSite& faces = mesh.getFaces();
    const MultiField::ArrayIndex cVarIndex(&_varField,&cells);
    CCMatrix& matrix = dynamic_cast<CCMatrix&>(mfmatrix.getMatrix(cVarIndex,
                                                             cVarIndex));

    const VectorT3Array& faceArea =
      dynamic_cast<const VectorT3Array&>(_geomFields.area[faces]);
    
    const TArray& faceAreaMag =
      dynamic_cast<const TArray&>(_geomFields.areaMag[faces]);

    const VectorT3Array& cellCentroid =
      dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[cells]);

    const VectorT3Array& faceCentroid =
      dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[faces]);

    const TArray& cellVolume =
      dynamic_cast<const TArray&>(_geomFields.volume[cells]);
    
    const CRConnectivity& faceCells = mesh.getAllFaceCells();

    CCAssembler& assembler = matrix.getPairWiseAssembler(faceCells);
    DiagArray& diag = matrix.getDiag();

    const XArray& xCell = dynamic_cast<const XArray&>(xField[cVarIndex]);
    XArray& rCell = dynamic_cast<XArray&>(rField[cVarIndex]);

    const GradArray& xGradCell =
      dynamic_cast<const GradArray&>(_varGradientField[cells]);

    const TArray& diffCell =
      dynamic_cast<const TArray&>(_diffusivityField[cells]);

    const IntArray& ibType = dynamic_cast<const IntArray&>(_geomFields.ibType[cells]);

    const int nFaces = faces.getCount();
    for(int f=0; f<nFaces; f++)
    {
        const int c0 = faceCells(f,0);
        const int c1 = faceCells(f,1);

        T_Scalar vol0 = cellVolume[c0];
        T_Scalar vol1 = cellVolume[c1];
        
        VectorT3 ds=cellCentroid[c1]-cellCentroid[c0];

        // for ib faces ignore the solid cell and use the face centroid for diff metric
        if (((ibType[c0] == Mesh::IBTYPE_FLUID)
             && (ibType[c1] == Mesh::IBTYPE_BOUNDARY)) ||
            ((ibType[c1] == Mesh::IBTYPE_FLUID)
             && (ibType[c0] == Mesh::IBTYPE_BOUNDARY)))
        {
            if (ibType[c0] == Mesh::IBTYPE_FLUID)
            {
                vol1 = 0.;
                ds = faceCentroid[f]-cellCentroid[c0];
            }
            else
            {
                vol0 = 0.;
                ds = cellCentroid[c1]-faceCentroid[f];
            }
        }
        
        T_Scalar faceDiffusivity(1.0);
        if (vol0 == 0.)
          faceDiffusivity = diffCell[c1];
        else if (vol1 == 0.)
          faceDiffusivity = diffCell[c0];
        else
          faceDiffusivity = harmonicAverage(diffCell[c0],diffCell[c1]);
        
        const T_Scalar diffMetric = faceAreaMag[f]*faceAreaMag[f]/dot(faceArea[f],ds);
        const T_Scalar diffCoeff = faceDiffusivity*diffMetric;
        const VectorT3 secondaryCoeff = faceDiffusivity*(faceArea[f]-ds*diffMetric);
        
        const XGrad gradF = (xGradCell[c0]*vol0+xGradCell[c1]*vol1)/(vol0+vol1);

        const X dFluxSecondary = gradF*secondaryCoeff;
        
        const X dFlux = diffCoeff*(xCell[c1]-xCell[c0]) + dFluxSecondary;

        rCell[c0][1] += dFlux[1];
        rCell[c1][1] -= dFlux[1];
        
        assembler.getCoeff01(f)[3] +=diffCoeff;
        assembler.getCoeff10(f)[3] +=diffCoeff;

        diag[c0][3] -= diffCoeff;
        diag[c1][3] -= diffCoeff;
    }
  }

Member Data Documentation

template<class X , class Diag , class OffDiag >

const Field& ElecDiffusionDiscretization< X, Diag, OffDiag >::_diffusivityField

private

Definition at line 163 of file ElecDiffusionDiscretization.h.

Referenced by ElecDiffusionDiscretization< X, Diag, OffDiag >::discretize().

template<class X , class Diag , class OffDiag >

const GeomFields& ElecDiffusionDiscretization< X, Diag, OffDiag >::_geomFields

private

Definition at line 161 of file ElecDiffusionDiscretization.h.

Referenced by ElecDiffusionDiscretization< X, Diag, OffDiag >::discretize().

template<class X , class Diag , class OffDiag >

Field& ElecDiffusionDiscretization< X, Diag, OffDiag >::_varField

private

Definition at line 162 of file ElecDiffusionDiscretization.h.

Referenced by ElecDiffusionDiscretization< X, Diag, OffDiag >::discretize().

template<class X , class Diag , class OffDiag >

const Field& ElecDiffusionDiscretization< X, Diag, OffDiag >::_varGradientField

private

Definition at line 164 of file ElecDiffusionDiscretization.h.

Referenced by ElecDiffusionDiscretization< X, Diag, OffDiag >::discretize().

---

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

• ElecDiffusionDiscretization.h