#include <BatteryPCBinaryElectrolyteDiscretization.h>

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

Collaboration diagram for BatteryPCBinaryElectrolyteDiscretization< 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 Gradient< T_Scalar >	TGrad

typedef Array< TGrad >	TGradArray

Public Member Functions
	BatteryPCBinaryElectrolyteDiscretization (const MeshList &meshes, const GeomFields &geomFields, Field &varField, const Field &diffusivityField, const Field &concField, const Field &concGradientField, const bool thermalModel, const int ElectrolyteMeshID)

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

const Field &	_concGradientField

const bool	_thermalModel

const int	_ElectrolyteMeshID

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

Detailed Description

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

Definition at line 26 of file BatteryPCBinaryElectrolyteDiscretization.h.

Member Typedef Documentation

template<class X , class Diag , class OffDiag >

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

Definition at line 34 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 32 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 33 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 37 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 30 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 40 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Gradient<T_Scalar> BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::TGrad

Definition at line 43 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

typedef Array<TGrad> BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::TGradArray

Definition at line 44 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 41 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 42 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 39 of file BatteryPCBinaryElectrolyteDiscretization.h.

template<class X , class Diag , class OffDiag >

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

Definition at line 36 of file BatteryPCBinaryElectrolyteDiscretization.h.

Constructor & Destructor Documentation

template<class X , class Diag , class OffDiag >

BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::BatteryPCBinaryElectrolyteDiscretization	(	const MeshList &	meshes,
		const GeomFields &	geomFields,
		Field &	varField,
		const Field &	diffusivityField,
		const Field &	concField,
		const Field &	concGradientField,
		const bool	thermalModel,
		const int	ElectrolyteMeshID
	)

inline

Definition at line 46 of file BatteryPCBinaryElectrolyteDiscretization.h.

                                                                       :
     Discretization(meshes),
       _geomFields(geomFields),
       _varField(varField),
       _diffusivityField(diffusivityField),
       _concField(concField),
       _concGradientField(concGradientField),
       _thermalModel(thermalModel),
       _ElectrolyteMeshID(ElectrolyteMeshID)
       {}

Member Function Documentation

template<class X , class Diag , class OffDiag >

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

inlinevirtual

Implements Discretization.

Definition at line 64 of file BatteryPCBinaryElectrolyteDiscretization.h.

References BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_concField, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_concGradientField, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_diffusivityField, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_ElectrolyteMeshID, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_geomFields, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_thermalModel, BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_varField, GeomFields::area, GeomFields::areaMag, GeomFields::coordinate, dot(), Mesh::getAllFaceGroups(), Mesh::getCells(), StorageSite::getCount(), Mesh::getFaceCells(), Mesh::getID(), harmonicAverage(), FaceGroup::site, and GeomFields::volume.

     {
       if (mesh.getID() == _ElectrolyteMeshID)
         {
           // electric potential variables
           const StorageSite& cells = mesh.getCells();
    
           const MultiField::ArrayIndex cVarIndex(&_varField,&cells);
 
           //CCMatrix& matrix = dynamic_cast<CCMatrix&>(mfmatrix.getMatrix(cVarIndex,cVarIndex));    
  
           const VectorT3Array& cellCentroid =
             dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[cells]);
           const TArray& cellVolume =
             dynamic_cast<const TArray&>(_geomFields.volume[cells]);
  
           //DiagArray& diag = matrix.getDiag();
 
           const XArray& xCell = dynamic_cast<const XArray&>(xField[cVarIndex]); // need for temp
           XArray& rCell = dynamic_cast<XArray&>(rField[cVarIndex]);
       
           // lnspeciesConcentration variables
           const TArray& lnSpecConcCell = dynamic_cast<const TArray&>(_concField[cells]);
 
           const TGradArray& lnSpecConcGradCell =
             dynamic_cast<const TGradArray&>(_concGradientField[cells]);
 
           const XArray& diffCell =
             dynamic_cast<const XArray&>(_diffusivityField[cells]);
 
           //const IntArray& ibType = dynamic_cast<const IntArray&>(_geomFields.ibType[cells]);
 
 
           foreach(const FaceGroupPtr fgPtr, mesh.getAllFaceGroups())
             {
               const FaceGroup& fg = *fgPtr;
 
               const StorageSite& faces = fg.site;
               const int nFaces = faces.getCount();
               const CRConnectivity& faceCells = mesh.getFaceCells(faces);
               const VectorT3Array& faceArea =
                 dynamic_cast<const VectorT3Array&>(_geomFields.area[faces]);    
               const TArray& faceAreaMag =
                 dynamic_cast<const TArray&>(_geomFields.areaMag[faces]);
               //const VectorT3Array& faceCentroid =dynamic_cast<const VectorT3Array&>(_geomFields.coordinate[faces]);
               //CCAssembler& assembler = matrix.getPairWiseAssembler(faceCells);
               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])[0];
                   else if (vol1 == 0.)
                     faceDiffusivity = (diffCell[c0])[0];
                   else
                     faceDiffusivity = harmonicAverage((diffCell[c0])[0],(diffCell[c1])[0]);
 
                   T_Scalar faceTemp(300.0);
                   if (_thermalModel)
                     {
                       if (vol0 == 0.)
                         faceTemp = (xCell[c1])[2];
                       else if (vol1 == 0.)
                         faceTemp = (xCell[c0])[2];
                       else
                         faceTemp = harmonicAverage((xCell[c0])[2],(xCell[c1])[2]);
                     }
 
 
                   //convert diffusivity for ln term in battery equation
                   const T_Scalar transportNumber = 0.4;
                   const T_Scalar R = 8.314;
                   const T_Scalar F = 96485.0; 
                   faceDiffusivity = faceDiffusivity*(-2.0)*R*faceTemp/F*(1.0-transportNumber);
              
                   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 TGrad gradF = (lnSpecConcGradCell[c0]*vol0+lnSpecConcGradCell[c1]*vol1)/(vol0+vol1);
 
                   const T_Scalar dFluxSecondary = gradF*secondaryCoeff;
         
                   const T_Scalar dFlux = diffCoeff*(lnSpecConcCell[c1]-lnSpecConcCell[c0]) + dFluxSecondary;
 
                   (rCell[c0])[0] += dFlux;
                   (rCell[c1])[0] -= dFlux;
 
                   //assembler.getCoeff01(f) +=diffCoeff;
                   //assembler.getCoeff10(f) +=diffCoeff;
 
                   //diag[c0] -= diffCoeff;
                   //diag[c1] -= diffCoeff;
 
                   //cout <<"dflux" << dFlux << endl;
                   /*
                     cout << " c0 and c1 " << c0 << " " << c1 << endl;
                     cout << " diffcell " << diffCell[c0] << " " << diffCell[c1] << " " << faceDiffusivity << endl;
                     cout << " diffmetric  " << diffMetric << endl;
                     cout << "diffCoeff  " << diffCoeff << endl;
                   */
         
               
 
                 }
             }
         }
     }

Member Data Documentation

template<class X , class Diag , class OffDiag >

const Field& BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_concField

private

Definition at line 203 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

const Field& BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_concGradientField

private

Definition at line 204 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

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

private

Definition at line 202 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

const int BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_ElectrolyteMeshID

private

Definition at line 206 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

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

private

Definition at line 200 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

const bool BatteryPCBinaryElectrolyteDiscretization< X, Diag, OffDiag >::_thermalModel

private

Definition at line 205 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

template<class X , class Diag , class OffDiag >

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

private

Definition at line 201 of file BatteryPCBinaryElectrolyteDiscretization.h.

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

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

BatteryPCBinaryElectrolyteDiscretization.h

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

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

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

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