#include <SourceDiscretizationene.h>

Inheritance diagram for SourceDiscretizationene< X >:

Collaboration diagram for SourceDiscretizationene< X >:

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

Public Member Functions
	SourceDiscretizationene (const MeshList &meshes, const GeomFields &geomFields, Field &varField, const Field &velocityField, Field &muField, Field &dissipationField, Field &densityField, Field &sourcekField, const Field &gradientField)

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 &	_dissipationField

Field &	_densityField

Field &	_sourcekField

const Field &	_gradientField

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

Detailed Description

template<class X>
class SourceDiscretizationene< X >

Definition at line 22 of file SourceDiscretizationene.h.

Member Typedef Documentation

template<class X >

typedef NumTypeTraits<X>::T_Scalar SourceDiscretizationene< X >::T_Scalar

Definition at line 25 of file SourceDiscretizationene.h.

template<class X >

typedef Array<T_Scalar> SourceDiscretizationene< X >::TArray

Definition at line 31 of file SourceDiscretizationene.h.

template<class X >

typedef Vector<T_Scalar,3> SourceDiscretizationene< X >::VectorT3

Definition at line 27 of file SourceDiscretizationene.h.

template<class X >

typedef Array<VectorT3> SourceDiscretizationene< X >::VectorT3Array

Definition at line 28 of file SourceDiscretizationene.h.

template<class X >

typedef Array<Gradient<VectorT3> > SourceDiscretizationene< X >::VGradArray

Definition at line 30 of file SourceDiscretizationene.h.

template<class X >

typedef Gradient<VectorT3> SourceDiscretizationene< X >::VGradType

Definition at line 29 of file SourceDiscretizationene.h.

template<class X >

typedef Array<X> SourceDiscretizationene< X >::XArray

Definition at line 32 of file SourceDiscretizationene.h.

Constructor & Destructor Documentation

template<class X >

SourceDiscretizationene< X >::SourceDiscretizationene	(	const MeshList &	meshes,
		const GeomFields &	geomFields,
		Field &	varField,
		const Field &	velocityField,
		Field &	muField,
		Field &	dissipationField,
		Field &	densityField,
		Field &	sourcekField,
		const Field &	gradientField
	)

inline

Definition at line 35 of file SourceDiscretizationene.h.

                                                       :
 
       Discretization(meshes),
       _geomFields(geomFields),
       _varField(varField),
       _velocityField(velocityField),
       _muField(muField),
       _dissipationField(dissipationField),
       _densityField(densityField),
       _sourcekField(sourcekField),
       _gradientField(gradientField)
 
    {}

Member Function Documentation

template<class X >

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

inlinevirtual

Implements Discretization.

Definition at line 56 of file SourceDiscretizationene.h.

References SourceDiscretizationene< X >::_densityField, SourceDiscretizationene< X >::_dissipationField, SourceDiscretizationene< X >::_geomFields, SourceDiscretizationene< X >::_gradientField, SourceDiscretizationene< X >::_muField, SourceDiscretizationene< X >::_sourcekField, SourceDiscretizationene< X >::_varField, Mesh::getCells(), StorageSite::getCount(), and GeomFields::volume.

   {
 
 
     const StorageSite& cells = mesh.getCells();
     const MultiField::ArrayIndex cVarIndex(&_varField,&cells);
     const TArray& cellVolume = dynamic_cast<const TArray&>(_geomFields.volume[cells]);
 
     const VGradArray& vGrad =
       dynamic_cast<const VGradArray&>(_gradientField[cells]);
 
     const TArray& muCell =
       dynamic_cast<TArray&>(_muField[cells]);
 
     const TArray& eCell =
       dynamic_cast<TArray&>(_dissipationField[cells]);
 
     const TArray & rhoCell =
       dynamic_cast<const TArray&>(_densityField[cells]);
 
     TArray & sourceCell =
       dynamic_cast<TArray&>(_sourcekField[cells]);
 
   
     TArray& rCell =
       dynamic_cast<XArray&>(rField[cVarIndex]);
    
   
     const int nCells = cells.getCount();  
     
     
     for(int n=0; n<nCells; n++)
     {
         const VGradType& vg = vGrad[n];
         VGradType vgSquare = vGrad[n]; 
         T_Scalar rhoe = eCell[n]*rhoCell[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];
          }
         }
 /*
         cout << "sum" << sum << endl;
         cout << "muCell" << muCell[n] << endl;
         cout << "rhoe" << rhoCell[n] << endl;
         cout << "epsilon" << eCell[n] << endl; 
  */
         
        sourceCell[n] =  sum*muCell[n]-rhoe;
         //cout << "rsource" << sourceCell[n]*cellVolume[n] << endl;
         rCell[n] +=sourceCell[n]*cellVolume[n];
     }
 }