The Semiconductor Interface includes quasi-Fermi level, density-gradient, linear, and logarithmic finite element formulations and a finite volume formulation. The formulation used is selected in the
Discretization section since the shape functions that can be used are directly related to the formulation employed.
The finite volume formulation uses constant shape functions, whilst the finite element formulations can use either linear or quadratic shape functions. In the different formulations the carrier concentration dependent variables (by default Ne and
Ph) represent different quantities. In the linear finite element and finite volume formulations
Ne = N and
Ph = P, where
N is the electron concentration and
P is the hole concentration. For the logarithmic finite element formulation
Ne = ln(N) and
Ph = ln(P). For the quasi-Fermi level and density-gradient formulation, the quasi-Fermi levels for the electrons and holes are the dependent variables.
To change the formulation, first expand the Discretization section. Then under
Discretization select a
Formulation (as in
Figure 2-1). Each formulation has advantages and disadvantages.
Any variables that involve expressions directly derived from the variables in Table 2-1 can also be used in expressions, for example, the electric field,
semi.norm
E, semi.EX, semi.EY, semi.EZ, or the total current,
semi.norm
J, semi.JX, semi.JY, semi.JZ.
The finite element formulation typically solves faster than the finite volume formulation. One reason is that, for an identical mesh, the finite element method with linear shape functions typically results in fewer degrees of freedom. In 2D, for triangular mesh elements, the number of degrees of freedom for the finite element method with linear shape functions is approximately half that for a finite volume discretization. Coupling to other physics interfaces is straightforward and variables can be differentiated using the d operator. The finite element method is an energy conserving method and thus current conservation is not implicit in the technique. Current conservation for the linear formulation is poor and this formulation is provided primarily for reasons of backward compatibility. Current conservation in the log formulation is much better but still not as good as the finite volume method. In order to help with numerical stability, a Galerkin least-squares stabilization method is included. This method usually enhances the ability to achieve a converged solution, particularly when using the linear formulation. However, it can be preferable to disable the stabilization, since the additional numerical diffusion that the technique introduces can produce slightly unphysical results. As a result of the reduced gradients in the dependent variables obtained when using the log formulation, stabilization is often not required when using this technique.