next up previous
Next: 3.2 Material Models Up: 3.1 The Device Simulator Previous: 3.1.3 Hydrodynamic Model Current

3.1.4 Lattice Heat Flow Equation

Additionally the heat flux equation is solved self consistently with the DD or HD transport equations. The heat flux equation reads:

    div$\displaystyle  {\bf S}_L = -\rho_L \cdot c_L \cdot \frac{\partial {\it T}_\mathrm{L}}{\partial t} + H$ (3.15)

and is solved in combination with Fourier's law:
    $\displaystyle {\bf {S}}_L = -\kappa_L \cdot$   grad$\displaystyle {\it T}_\mathrm{L}$ (3.16)

$ S_L$ denotes the lattice heat flux density, $ \rho_L$ the mass density, $ c_L$ the heat capacity, and $ \kappa_L$ the thermal conductivity. $ {\it T}_\mathrm{L}$ represents the lattice temperature. The implementation of the heat generation term $ H$ depends on the transport model. For the DD case the carriers are assumed to energetically relax instantly which is equivalent to an infinitely small relaxation time. Thus, they deposit the energy accepted by the field instantly to the lattice. For the recombination heat it is assumed that every carrier provides the energy equivalent to the band gap energy [147]. Thus, the generation term $ H$ reads:


    $\displaystyle H = $grad$\displaystyle \bigg(\frac{E_C}{q}-\psi\bigg)\cdot {\bf {J}}_n +$   grad$\displaystyle \bigg(\frac{E_V}{q}-\psi \bigg)\cdot {\bf {J}}_p + R\cdot (E_C - E_V)$ (3.17)

For the hydrodynamic case for $ H$ the following term is used:
    $\displaystyle H = \frac{3\cdot {\it k}_{\mathrm{B}}}{2} \bigg(n\cdot\frac{T_n-{...
...}}+ p\cdot\frac{T_p-{\it T}_\mathrm{L}}{\tau_{w,p}}\bigg) + R \cdot (E_C - E_V)$ (3.18)

(3.19) is obtained from the stationary energy flux balance derived from the second moment solving the Boltzmann equation and the carrier continuity equations (3.2) and (3.3) for the hydrodynamic model used [49,147]. The first two terms of (3.19) represent energy relaxation approaches approximating the scattering terms of the Boltzmann equation [147]. The third term represents an effective recombination term. For the recombination a similar approximation is taken as in the DD case. This is sound and introduces a small error only, as the recombination energies generally are smaller than the Joule like contribution in (3.18) and (3.19). This is different for intensive light irradiation where additional contributions have to be added in order to model the effective recombination term.


next up previous
Next: 3.2 Material Models Up: 3.1 The Device Simulator Previous: 3.1.3 Hydrodynamic Model Current
Quay
2001-12-21