Memosa-FVM  0.2
Public Types | Public Member Functions | Private Attributes | List of all members
TrapBandTunnelingDiscretization< X, Diag, OffDiag > Class Template Reference

#include <TrapBandTunnelingDiscretization.h>

Inheritance diagram for TrapBandTunnelingDiscretization< X, Diag, OffDiag >:
Inheritance graph
[legend]
Collaboration diagram for TrapBandTunnelingDiscretization< X, Diag, OffDiag >:
Collaboration graph
[legend]

Public Types

typedef NumTypeTraits< X >
::T_Scalar 		T_Scalar
 
typedef Array< T_Scalar > TArray
 
typedef CRMatrix< Diag,
OffDiag, X > 		CCMatrix
 
typedef CCMatrix::DiagArray DiagArray
 
typedef CCMatrix::OffDiagArray OffDiagArray
 
typedef Array< X > XArray
 
typedef Array< Vector
< T_Scalar, 3 > > 		VectorT3Array
 

Public Member Functions

 TrapBandTunnelingDiscretization (const MeshList &meshes, const GeomFields &geomFields, const Field &varField, const Field &electricField, const Field &conductionbandField, const ElectricModelConstants< T_Scalar > &constants)
 
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
 
const Field & _varField
 
const Field & _electricField
 
const Field & _conductionbandField
 
const ElectricModelConstants
< T_Scalar > & 		_constants
 

Additional Inherited Members

- Protected Attributes inherited from Discretization
const MeshList & _meshes
 

Detailed Description

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

Definition at line 21 of file TrapBandTunnelingDiscretization.h.

Member Typedef Documentation

template<class X , class Diag , class OffDiag >
typedef CRMatrix<Diag,OffDiag,X> TrapBandTunnelingDiscretization< X, Diag, OffDiag >::CCMatrix

Definition at line 26 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef CCMatrix::DiagArray TrapBandTunnelingDiscretization< X, Diag, OffDiag >::DiagArray

Definition at line 27 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef CCMatrix::OffDiagArray TrapBandTunnelingDiscretization< X, Diag, OffDiag >::OffDiagArray

Definition at line 28 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef NumTypeTraits<X>::T_Scalar TrapBandTunnelingDiscretization< X, Diag, OffDiag >::T_Scalar

Definition at line 24 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef Array<T_Scalar> TrapBandTunnelingDiscretization< X, Diag, OffDiag >::TArray

Definition at line 25 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef Array<Vector<T_Scalar, 3> > TrapBandTunnelingDiscretization< X, Diag, OffDiag >::VectorT3Array

Definition at line 30 of file TrapBandTunnelingDiscretization.h.

template<class X , class Diag , class OffDiag >
typedef Array<X> TrapBandTunnelingDiscretization< X, Diag, OffDiag >::XArray

Definition at line 29 of file TrapBandTunnelingDiscretization.h.

Constructor & Destructor Documentation

template<class X , class Diag , class OffDiag >
TrapBandTunnelingDiscretization< X, Diag, OffDiag >::TrapBandTunnelingDiscretization ( const MeshList &  meshes,
const GeomFields &  geomFields,
const Field &  varField,
const Field &  electricField,
const Field &  conductionbandField,
const ElectricModelConstants< T_Scalar > &  constants 
)
inline

Definition at line 32 of file TrapBandTunnelingDiscretization.h.

37  :
38  Discretization(meshes),
39  _geomFields(geomFields),
40  _varField(varField),
41  _electricField(electricField),
42  _conductionbandField(conductionbandField),
43  _constants(constants)
44  {}
Discretization::Discretization
Discretization(const MeshList &meshes)
Definition: Discretization.cpp:8
TrapBandTunnelingDiscretization::_constants
const ElectricModelConstants< T_Scalar > & _constants
Definition: TrapBandTunnelingDiscretization.h:247

Member Function Documentation

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

Implements Discretization.

Definition at line 46 of file TrapBandTunnelingDiscretization.h.

References TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_conductionbandField, TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_constants, TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_electricField, TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_geomFields, TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_varField, GeomFields::coordinate, ElectricModelConstants< T >::electron_trapdepth, fabs(), Mesh::getCellCells(), Mesh::getCells(), StorageSite::getCount(), CRConnectivity::getCount(), CRMatrix< T_Diag, T_OffDiag, X >::getDiag(), MultiFieldMatrix::getMatrix(), StorageSite::getSelfCount(), HBAR_SI, mag(), ME, PI, PositiveValueOf(), QE, sqrt(), and GeomFields::volume.

48  {
49  const StorageSite& cells = mesh.getCells();
50  const int nCells = cells.getSelfCount();
51 
52  const TArray& conduction_band = dynamic_cast<const TArray&> (_conductionbandField[cells]);
53 
54  const VectorT3Array& electric_field = dynamic_cast<const VectorT3Array&> (_electricField[cells]);
55 
56  const VectorT3Array& cellCentroid =
57  dynamic_cast<const VectorT3Array& > (_geomFields.coordinate[cells]);
58 
59  const MultiField::ArrayIndex cVarIndex(&_varField,&cells);
60 
61  CCMatrix& matrix = dynamic_cast<CCMatrix&>(mfmatrix.getMatrix(cVarIndex,cVarIndex));
62 
63  const XArray& xCell = dynamic_cast<const XArray&>(_varField[cells]);
64 
65  XArray& rCell = dynamic_cast<XArray&>(rField[cVarIndex]);
66 
67  DiagArray& diag = matrix.getDiag();
68 
69  const TArray& cellVolume =
70  dynamic_cast<const TArray&>(_geomFields.volume[cells]);
71 
72  const CRConnectivity& cellCells = mesh.getCellCells();
73 
74  TArray* ts = (new TArray(cells.getCount()));
75  *ts = 0;
76  TArray& transmission = *ts;
77 
78  const T_Scalar electron_capture_cross = _constants["electron_capture_cross"];
79  const T_Scalar electron_effmass = _constants["electron_effmass"];
80  const int normal = _constants["normal_direction"];
81  const int nTrap = _constants["nTrap"];
82  vector<T_Scalar> electron_trapdepth = _constants.electron_trapdepth;
83  //const T_Scalar epsilon = 1e-18;
84 
85  const int nMax = 200;
86 
87  T_Scalar transmissionHigh(0), transmissionLow(0);
88 
89  int idHigh, idLow, high, low, me, count;
90 
91  bool foundHigh, foundLow, flag;
92 
93  idHigh = idLow = 0;
94 
95  foundHigh = foundLow = false;
96 
97  for(int c=0; c<cells.getCount(); c++){
98  transmission[c] = 0.0;
99  }
100 
101  for(int c=0; c<nCells; c++){
102 
103  for(int i=0; i<nTrap; i++){
104 
105  T_Scalar en = conduction_band[c] - electron_trapdepth[i];
106 
107  transmission[c] = 1.0;
108 
109  //-------------------------------------------------------------------------//
110 
111  high = low = me = c;
112 
113  flag = false; //hit the boundary or not?
114 
115  count = 0;
116 
117  while(flag == false && count < nMax ){
118 
119  const int nbc = cellCells.getCount(me);
120 
121  T_Scalar drmin = 0.0;
122 
123  int neighborUp = 0;
124 
125  for(int nc = 0; nc < nbc; nc++){
126 
127  const int neighbor = cellCells(me, nc);
128  const T_Scalar dr = cellCentroid[me][normal] - cellCentroid[neighbor][normal];
129  if (dr < drmin){
130  drmin = dr;
131  neighborUp = neighbor;
132  }
133  }
134  if (neighborUp < nCells) {
135  high = neighborUp;
136  low = me;
137  me = high;
138  }
139  else flag = true;
140 
141  T_Scalar dX = cellCentroid[me][normal] - cellCentroid[low][normal];
142  T_Scalar factor = -2.0/HBAR_SI * sqrt(2.0*electron_effmass*ME*QE);
143  T_Scalar valueMe = PositiveValueOf( conduction_band[me] - en);
144  T_Scalar valueLow = PositiveValueOf( conduction_band[low] - en);
145  T_Scalar avg = (valueMe + valueLow) / 2.0;
146  T_Scalar exponent = factor * sqrt(avg) * fabs(dX);
147 
148  transmission[me] = transmission[low] * exp(exponent);
149 
150  if (en - conduction_band[me] >0 ){
151 
152  foundHigh = true;
153  idHigh = me;
154  transmissionHigh = transmission[me];
155  // cout << "found high " << c <<" " << me << " " << transmissionHigh << endl;
156  break;
157 
158  }
159  count ++;
160  }
161 
162 #if 0
163  //-------------------------------------------------------------------------//
164 
165  high = low = me = c;
166 
167  flag = false;
168 
169  count = 0;
170 
171  while(flag == false && count < nMax){
172 
173  const int nbc = cellCells.getCount(me);
174 
175  T_Scalar drmin = 0.0;
176 
177  int neighborUp = 0;
178 
179  for(int nc = 0; nc < nbc; nc++){
180  const int neighbor = cellCells(me, nc);
181  const T_Scalar dr = cellCentroid[me][normal] - cellCentroid[neighbor][normal];
182  if (dr > drmin){
183  drmin = dr;
184  neighborUp = neighbor;
185  }
186  }
187 
188  if (neighborUp < nCells) {
189  high = neighborUp;
190  low = me;
191  me = high;
192  }
193  else flag = true;
194 
195  T_Scalar dX = cellCentroid[me][normal] - cellCentroid[low][normal];
196  T_Scalar factor = -2.0/HBAR_SI * sqrt(2.0*electron_effmass*ME*QE);
197  T_Scalar valueMe = PositiveValueOf( conduction_band[me] - en);
198  T_Scalar valueLow = PositiveValueOf( conduction_band[low] - en);
199  T_Scalar avg = (valueMe + valueLow) / 2.0;
200  T_Scalar exponent = factor * sqrt(avg) * fabs(dX);
201 
202  transmission[me] = transmission[low] * exp(exponent);
203 
204  if (en - conduction_band[me] >0 ){
205 
206  foundLow = true;
207  idLow = me;
208  transmissionLow = transmission[me];
209  //cout << "found low " << c <<" " << me << endl;
210  break;
211  }
212  // cout << c << " low " << low << endl;
213  count ++;
214 
215  }
216 #endif
217  //-------------------------------------------------------------------------//
218 
219  const T_Scalar ef = mag(electric_field[c]);
220 
221  const T_Scalar alpha = cellVolume[c] * QE * ef * ef * electron_capture_cross /
222  (16 * PI*PI * HBAR_SI * electron_effmass * electron_trapdepth[i]);
223 
224  if (foundLow == true){
225  rCell[c][i] -= alpha * transmissionLow * xCell[c][i];
226  diag[c](i,i) -= alpha * transmissionLow;
227  //diag[c][i] -= alpha * transmissionLow;
228  rCell[idLow][nTrap] += alpha * transmissionLow * xCell[c][i];
229  }
230 
231  if (foundHigh == true){
232  rCell[c][i] -= alpha * transmissionHigh * xCell[c][i];
233  diag[c](i,i) -= alpha * transmissionHigh;
234  //diag[c][i] -= alpha * transmissionHigh;
235  rCell[idHigh][nTrap] += alpha * transmissionHigh * xCell[c][i];
236  }
237  }
238  }
239  }
CRConnectivity::getCount
int getCount(const int i) const
Definition: CRConnectivity.h:165
MultiFieldMatrix::getMatrix
Matrix & getMatrix(const Index &rowIndex, const Index &colIndex)
Definition: MultiFieldMatrix.cpp:42
PI
#define PI
Definition: PhysicsConstant.h:10
StorageSite::getSelfCount
int getSelfCount() const
Definition: StorageSite.h:40
GeomFields::coordinate
Field coordinate
Definition: GeomFields.h:19
mag
T mag(const Vector< T, 3 > &a)
Definition: Vector.h:260
ME
#define ME
Definition: PhysicsConstant.h:16
sqrt
Tangent sqrt(const Tangent &a)
Definition: Tangent.h:317
HBAR_SI
#define HBAR_SI
Definition: PhysicsConstant.h:22
MultiField::ArrayIndex
pair< const Field *, const StorageSite * > ArrayIndex
Definition: MultiField.h:21
QE
#define QE
Definition: PhysicsConstant.h:18
Mesh::getCellCells
const CRConnectivity & getCellCells() const
Definition: Mesh.cpp:480
Mesh::getCells
const StorageSite & getCells() const
Definition: Mesh.h:109
GeomFields::volume
Field volume
Definition: GeomFields.h:26
PositiveValueOf
const T PositiveValueOf(T input)
Definition: ElectricUtilityFunctions.h:50
fabs
Tangent fabs(const Tangent &a)
Definition: Tangent.h:312
StorageSite::getCount
int getCount() const
Definition: StorageSite.h:39
ElectricModelConstants::electron_trapdepth
vector< T > electron_trapdepth
Definition: ElectricBC.h:77
TrapBandTunnelingDiscretization::_constants
const ElectricModelConstants< T_Scalar > & _constants
Definition: TrapBandTunnelingDiscretization.h:247

Member Data Documentation

template<class X , class Diag , class OffDiag >
const Field& TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_conductionbandField
private

Definition at line 246 of file TrapBandTunnelingDiscretization.h.

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

template<class X , class Diag , class OffDiag >
const ElectricModelConstants<T_Scalar>& TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_constants
private

Definition at line 247 of file TrapBandTunnelingDiscretization.h.

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

template<class X , class Diag , class OffDiag >
const Field& TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_electricField
private

Definition at line 245 of file TrapBandTunnelingDiscretization.h.

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

template<class X , class Diag , class OffDiag >
const GeomFields& TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_geomFields
private

Definition at line 243 of file TrapBandTunnelingDiscretization.h.

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

template<class X , class Diag , class OffDiag >
const Field& TrapBandTunnelingDiscretization< X, Diag, OffDiag >::_varField
private

Definition at line 244 of file TrapBandTunnelingDiscretization.h.

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

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