Code contains numerical solution of drift diffusion equations system for gunn diode with ballast load and ideal voltage source. In spite of this, the driftdiffusion dd model of electronic device operation is still used in nearly all lineof business device simulations 1. We present the numerical methods and simulations used to solve a charge transport problem in semiconductor physics. Anisotropic effects are incorporated into the basic drift diffusion equations by representing the mobilities with second rank tensors e.
This tutorial example computes the electron number density and mean electron energy in a drift tube. To get the numerical solution, the cranknicolson finite difference method is constructed, which is secondorder accurate in time and space. Wiastesca twodimensional semiconductor analysis package is a program system for the numerical simulation of charge transfer processes in semiconductor structures, especially also in semiconductor lasers. The implemented algorithm is based on the driftdiffusion approach, a set of equations which describes the microscopic operation of the semiconductor devices from a physics point of view.
Combined electromagnetic and drift diffusion models for. This phenomenon can be particularly egregious when the system should not be diffusive at all, for example an ideal fluid acquiring some spurious viscosity in a numerical model. Numerical simulation of groundwater pollution problems based on convection diffusion equation lingyu li, zhe yin college of mathematics and statistics, shandong normal university, jinan, china abstract the analytical solution of the convection diffusion equation is considered by twodimensional fourier transform and the inverse fourier transform. A quick short form for the diffusion equation is ut. Steadystate simulation of semiconductor devices using. The first part examines semiclassical transport methods, including drift diffusion, hydrodynamic, and monte carlo methods for solving the boltzmann transport equation. Simulation of stable domain mode generation in gunn diode.
In this work, we develop numerical simulation of the 2d drift diffusion model ddm. Details regarding numerical implementation and sample codes are provided as templates for sophisticated simulation software. Physics based models that allow for good matching between simulated behavior and measurements of insb material parameters have been developed. The 2d drift diffusion model of the software sibidif, thus, has been modi. Taylorgalerkin bspline finite element method for the onedimensional advection diffusion equation, numerical methods for. To solve the basic equations of the model, we developed the explicit and implicit techniques of driftdiffusion numerical simulation and applied software. Numerical simulation of reactiondiffusion equations on. In this method, the equations are solved sequentially. The problem is described by a wignerpoisson kinetic system we have recently proposed and whose results are in good agreement with known experiments. Numerical method for a 2d drift diffusion model arising in. Comparison and numerical treatment of generalised nernstplanck models, computer physics communications 196 2015, pp. Solution to the driftdiffusion equation maths partner. Abstract we present a new software simulation tool sesame, which solves the driftdiffusionpoisson equations in 1 and 2dimensions.
In particular, we discuss the qualitative properties of. Handbook of optoelectronic device modeling and simulation, j. A method for simulation of insb devices in a commercially available drift diffusion simulator has been demonstrated. Numerical diffusion is a difficulty with computer simulations of continua such as fluids wherein the simulated medium exhibits a higher diffusivity than the true medium. Instances when driftdiffusion equation can represent the trend or predict the mean behavior of the transport properties feature length of the semiconductors smaller than the mean free path of the carriers instances when drift diffusion equations are accurate quasisteady state assumption holds no transient effects. Prime members enjoy free twoday delivery and exclusive access to music, movies, tv shows, original audio series, and kindle books. We present a new approach to simulate the transport of charges across organicorganic layer interfaces in organic semiconductor devices.
Currently, ness contains a drift diffusion dd, kubogreenwood kg, and nonequilibrium greens function negf solvers. In the mathematical modeling and numerical simulation of semiconductor devices, the drift di. In this chapter, we start with a brief introduction to numerical simulation of transport phenomena. As first, we present a onedimensional 1 d pin diode structure simulation achieved by solving the drift diffusion model ddm. Simulation of insb devices using driftdiffusion equations. Introduction simulation tools capable of numerically characterizing semiconductor devices play a vital role in devicesystem design frameworks used by the electronics industry as.
The charge transport models and corresponding equations are implemented within the numerical solvers module which are solved selfconsistently with poisson equation. In this model, we consider doped semiconductor superlattices in which electrons. In this work, we present a numerical model to solve the drift diffusion equations coupled with electromagnetic model, where all simulations codes are implemented using matlab code software. In order to save space, the verification is shown in the reference 12. The convectiondiffusion equation is a combination of the diffusion and convection advection equations, and describes physical phenomena where particles, energy, or other physical quantities are transferred inside a physical system due to two processes. Refined driftdiffusion model for the simulation of charge. The non linear poisson equation has been discretized using the galerkin finite element method following a displacement formulation and the continuity equations have been treated using the eafe scheme. A robust finite element vector formulation of the semiconductor drift diffusion equations incorporat simulation of semiconductor processes and devices. The simplest ntype doping profile for diode has been chosen to make it able to generate. Di proc package provides a simulation of di usion processes and the di erences methods of simulation of solutions for stochastic di erential equations sdes of the itos type, in nancial and actuarial modeling and other areas of applications, for example the stochastic modeling and simulation of pollutant dispersion. The drift diffusion equation and its applications in. Institute of computational physics icp team interdisciplinary team of 8 physicists. Advanced simulation methods for charge transport in oleds.
Numerical simulation of groundwater pollution problems. In the numerical solution, the linearized equations are then decoupled with the gummels method. The drift diffusion equations are the most widely used model to describe semiconductor devices today. These aids are suitable for the stationary and timedomain simulation of injection lasers and photodetectors with various electrophysical, constructive, and technological parameters at different control actions. It is based on the drift diffusion model and considers a multitude of additional physical. Journal of surfactants and detergents jsd, a journal of the american oil chemists society aocs publishes basic and applied scientific research related to the interfacial behavior of petrochemical and oleochemical surfactants and detergent ingredients and how this relates to performance in different applications. It can be derived from the eulerpoisson equations when the relaxation time goes to 0. These simulation tools typically require other software frameworks e. The vast majority of device simulations are normally based on the numerical solution of approximate models which are related to the boltzmann equation, coupled to poissons equation for. Indeed, the continuity equations describe the evolution of carriers in the silicon nanowire film sourcechanneldrain in order to maintain a constant current along the film in which there is no charge accumulation. The charge transport in the nanowire mosfet simulated here. Index termsdiscontinuous galerkin method, drift diffusion equations, multiphysics modeling, poisson equation, semiconductor device modeling.
The bulk of the literature on mathematical models for device simulation is concerned with this nonlinear system of partial differential equations and numerical software for its solution is commonplace at practically every research facility in the field. Electrons are released due to thermionic emission on the left boundary with an. Analytical solutions of transient driftdiffusion in pn. A robust finite element vector formulation of the semiconductor drift diffusion equations incorporating anisotropic transport properties. More precisely in this work the gummel map algorithm is employed to solve the drift diffusion model for semiconductors. Wiastesca modeling and simulation of semiconductor devices product. Where ei is the unit vector along the i axis and e, is defined as en, vivi + pl,p,v,v 2. Solution of drift diffusion equations are conducted with fast implicit finitedifference method euler. The algorithm solves the whole set of equations numerically, determining the electrostatic potential and carrier densities at each point within the onedimensional geometry depicting the simulated structure. The drift diffusion interface solves a pair of reactionadvection diffusion equations, one for the electron density and the other for the mean electron energy.
Driftdiffusion simulation of highspeed optoelectronic. We present the software and method cadiff which allows for numerical simulation of interdiffusion and reactive diffusion processes. The reason is that the dd model has adequately explained or predicted the behavior of commercially important electronic devices through this rapid technology advancement. This system consists of the continuity equations for particle densities and a poissonequation for electrostatic potential.
Sesame is distributed both as an open source python package and as a standalone executable for windows. Numerical solution of advectiondiffusion equation using a. Quantum corrections to classical models bryan biegel nas device modeling workshop, august 78, 1997 projects wigner function and transfermatrix modeling of macroscopic quantum devices in 3d quantum corrections to classical drift diffusion and hydrodynamic models in 3d. In this study, we present a numerical scheme to so lve the driftdiffusion traffic flow model under the steady state. The famous diffusion equation, also known as the heat equation, reads. The computed results showed that the use of the current method in the simulation is very applicable for the solution of the advection diffusion equation. Numerical simulation shows excellent agreement with the analytical solution. Hsunjung cho and shihching lo department of transportationtechnology and management national chiao tung university 1001, ta hsueh rd. Numerical methods for a quantum driftdiffusion equation. It is shown that the method keeps the kinetics from overshooting the stable branches when a large time step is used in the simulation. In this work we consider the driftdiffusion model for which the current density is given by equation 5. We consider mathematical models that express certain conservation principles and consist of convection diffusion reactionequations written in integral, differential, or weak form. A guide to numerical methods for transport equations.
95 458 49 1083 1146 1086 1309 569 234 212 545 885 980 996 90 57 1266 1358 552 843 360 1380 646 998 539 552 352 481 1170 424 1238 1495 303 1293 1127 1226 1451 434 1203 265 1012