SourceDiscretizationdissi< X, Diag, OffDiag > Class Template Reference

#include <SourceDiscretizationdissi.h>

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

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

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

Public Member Functions
	SourceDiscretizationdissi (const MeshList &meshes, const GeomFields &geomFields, Field &varField, const Field &velocityField, Field &muField, Field &energyField, Field &densityField, Field &sourcedField, Field &sourcecField, Field &sourcepField, const Field &gradientField)
KeModelOptions< T_Scalar > &	getOptions ()
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 MeshList	_meshes
const GeomFields &	_geomFields
Field &	_varField
const Field &	_velocityField
Field &	_muField
Field &	_energyField
Field &	_densityField
Field &	_sourcedField
Field &	_sourcecField
Field &	_sourcepField
const Field &	_gradientField
KeModelOptions< T_Scalar >	_options

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

Detailed Description

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

Definition at line 22 of file SourceDiscretizationdissi.h.

Member Typedef Documentation

template<class X , class Diag , class OffDiag >

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

Definition at line 35 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 33 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 34 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 25 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 31 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 27 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 28 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

typedef Array<Gradient<VectorT3> > SourceDiscretizationdissi< X, Diag, OffDiag >::VGradArray

Definition at line 30 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

typedef Gradient<VectorT3> SourceDiscretizationdissi< X, Diag, OffDiag >::VGradType

Definition at line 29 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 32 of file SourceDiscretizationdissi.h.

Constructor & Destructor Documentation

template<class X , class Diag , class OffDiag >

SourceDiscretizationdissi< X, Diag, OffDiag >::SourceDiscretizationdissi ( const MeshList &  meshes, const GeomFields &  geomFields, Field &  varField, const Field &  velocityField, Field &  muField, Field &  energyField, Field &  densityField, Field &  sourcedField, Field &  sourcecField, Field &  sourcepField, const Field &  gradientField  )

inline

Definition at line 38 of file SourceDiscretizationdissi.h.

                                                        :

      Discretization(meshes),
      _geomFields(geomFields),
      _varField(varField),
      _velocityField(velocityField),
      _muField(muField),
      _energyField(energyField),
      _densityField(densityField),
      _sourcedField(sourcedField),
      _sourcecField(sourcecField),
      _sourcepField(sourcepField),

      _gradientField(gradientField)

   {}

Member Function Documentation

template<class X , class Diag , class OffDiag >

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

inlinevirtual

Implements Discretization.

Definition at line 67 of file SourceDiscretizationdissi.h.

References SourceDiscretizationdissi< X, Diag, OffDiag >::_densityField, SourceDiscretizationdissi< X, Diag, OffDiag >::_energyField, SourceDiscretizationdissi< X, Diag, OffDiag >::_geomFields, SourceDiscretizationdissi< X, Diag, OffDiag >::_gradientField, SourceDiscretizationdissi< X, Diag, OffDiag >::_muField, SourceDiscretizationdissi< X, Diag, OffDiag >::_options, SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcecField, SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcedField, SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcepField, SourceDiscretizationdissi< X, Diag, OffDiag >::_varField, KeModelOptions< T >::c1mu, KeModelOptions< T >::c2mu, Mesh::getCells(), StorageSite::getCount(), CRMatrix< T_Diag, T_OffDiag, X >::getDiag(), MultiFieldMatrix::getMatrix(), and GeomFields::volume.

  {


    const StorageSite& cells = mesh.getCells();
    const MultiField::ArrayIndex cVarIndex(&_varField,&cells);

    const VGradArray& vGrad =
      dynamic_cast<const VGradArray&>(_gradientField[cells]);

    const TArray& muCell =
      dynamic_cast<const TArray&>(_muField[cells]);

    const TArray& eCell =
      dynamic_cast<const TArray&>(_varField[cells]);

    const TArray& kCell =
      dynamic_cast<const TArray&>(_energyField[cells]);

    const TArray & rhoCell =
      dynamic_cast<const TArray&>(_densityField[cells]);

    TArray& sourceCell =
      dynamic_cast<TArray&>(_sourcedField[cells]);

    TArray& sourcecCell =
      dynamic_cast<TArray&>(_sourcecField[cells]);

    TArray& sourcepCell =
      dynamic_cast<TArray&>(_sourcepField[cells]);




    TArray& rCell = 
      dynamic_cast<XArray&>(rField[cVarIndex]);


    const int nCells = cells.getCount();
    const TArray& cellVolume = dynamic_cast<const TArray&>(_geomFields.volume[cells]);

    CCMatrix& matrix = dynamic_cast<CCMatrix&>(mfmatrix.getMatrix(cVarIndex,cVarIndex));

    DiagArray& diag = matrix.getDiag();
    T_Scalar two(2.0);
    T_Scalar C1mu = _options.c1mu;
    T_Scalar C2mu = _options.c2mu;
    for(int n=0; n<nCells; n++)
    {
        const VGradType& vg = vGrad[n];
        VGradType vgSquare = vGrad[n];
        T_Scalar sourcecoeff1  = (muCell[n]*eCell[n]*C1mu)/kCell[n];
        T_Scalar sourcecoeff2  = (C2mu*pow(eCell[n],two)*rhoCell[n])/kCell[n];
        T_Scalar dsc1 = (muCell[n]*C1mu)/kCell[n];
        T_Scalar dsc2 = (C2mu*2*eCell[n]*rhoCell[n])/kCell[n];

         T_Scalar sum = 0;
 
        for(int i=0;i<3;i++)
        {
          for(int j=0;j<3;j++)
           {
           T_Scalar x =  vg[i][j]*vg[i][j]+vg[i][j]*vg[j][i] ;
           sum += x;
           //vgSquare[i][j] =  vg[i][j]*vg[i][j]+vg[i][j]*vg[j][i] ;
          // sum += vgSquare[i][j];

           }
        }
        
        sourceCell[n] = sum*sourcecoeff1-sourcecoeff2;
        T_Scalar ds = sum*dsc1-dsc2;

        sourcecCell[n] = sourceCell[n]- ds*eCell[n];
        sourcepCell[n] = ds;
        rCell[n] +=sourcecCell[n]*cellVolume[n];
        diag[n] -=cellVolume[n]*sourcepCell[n];

    }

}

template<class X , class Diag , class OffDiag >

KeModelOptions<T_Scalar>& SourceDiscretizationdissi< X, Diag, OffDiag >::getOptions ( )

inline

Definition at line 64 of file SourceDiscretizationdissi.h.

References SourceDiscretizationdissi< X, Diag, OffDiag >::_options.

{return _options;}

Member Data Documentation

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_densityField

private

Definition at line 157 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_energyField

private

Definition at line 156 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

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

private

Definition at line 152 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

const Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_gradientField

private

Definition at line 161 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

const MeshList SourceDiscretizationdissi< X, Diag, OffDiag >::_meshes

private

Definition at line 151 of file SourceDiscretizationdissi.h.

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_muField

private

Definition at line 155 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

KeModelOptions<T_Scalar> SourceDiscretizationdissi< X, Diag, OffDiag >::_options

private

Definition at line 162 of file SourceDiscretizationdissi.h.

Referenced by SourceDiscretizationdissi< X, Diag, OffDiag >::discretize(), and SourceDiscretizationdissi< X, Diag, OffDiag >::getOptions().

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcecField

private

Definition at line 159 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcedField

private

Definition at line 158 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_sourcepField

private

Definition at line 160 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

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

private

Definition at line 153 of file SourceDiscretizationdissi.h.

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

template<class X , class Diag , class OffDiag >

const Field& SourceDiscretizationdissi< X, Diag, OffDiag >::_velocityField

private

Definition at line 154 of file SourceDiscretizationdissi.h.

---

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

• SourceDiscretizationdissi.h