A contribution to the simulation of Vlasov-based models

Abstract

Esta tesis está dedicada al desarrollo, aplicación y test de métodos para la simulación numérica de problemas procedentes de la física y de la ingeniería electrónica. La principal herramienta aplicada a lo largo de todo el trabajo es la ecuación de Vlasov (transporte) en la forma de la Boltzmann Transport Equation (BTE) para la descripción del transporte de partículas cargadas en plasmas y dispositivos electrónicos: las cargas se mueven bajo el efecto de un campo de fuerza y sufren scattering debido a otras cargas o fonones (pseudo-partículas que describen de manera efectiva las vibraciones de los iones del retículo cristalino).<br/>La BTE ha de ser acoplada con una ecuación o sistema de ecuaciones para calcular el campo de fuerza: para estructuras simples se usa la ecuación de Poisson; para plasmas, donde los efectos magnéticos no se pueden despreciar debido a las altas velocidades de las partículas, se usa la fuerza de Lorentz, por lo cual se han de resolver las ecuaciones de Maxwell; en nanoestructuras, por ejemplo transistores con dimensiones confinadas, la ecuación de Poisson necesita ser acoplada con la ecuación de Schrödinger para la descripción de las dimensiones cuánticas y para la descomposición en sub-bandas, o niveles de energía.<br/>Las colisiones son el scattering que las cargas padecen debido a las interacciones con otras cargas o con el retículo cristalino fijo, representado en forma de fonones. En la tesis se emplean diversos operadores de scattering: los más simples son operadores lineales de relajación; se estudia un modelo para la simulación de semiconductores donde se tienen en cuenta colisiones con fonones acústicos, en aproximación elástica, y fonones ópticos.<br/>Tras la introducción, en el primer capítulo se desarrollan los métodos numéricos más importantes: primero un método de interpolación no oscilante (PWENO), necesario para evitar las oscilaciones producidas por la reconstrucción por polinomios de Lagrange, que incrementa la variación total cuando aparecen choques: las oscilaciones en el espacio de fases son características del problema, pero si el método añade oscilaciones espúreas (es decir, debidas al método en sí), entonces el resultado numérico no tiene sentido, o simplemente explota. El segundo método numérico fundamental es la técnica de splitting: cuando se resuelve un problema complicado, si se puede dividir en sub-problemas y resolverlos por separado, entonces se puede reconstruir una aproximación para el problema completo; esta técnica se usa para el time splitting (separación de la parte de transporte y de colisión) y el splitting dimensional (dividir el espacio de fases en posición y velocidad). La tercera herramienta fundamental es un sólver para advección lineal: se usan dos métodos, uno basado en trazar hacia atrás las características a nivel puntual y otro basado en reconstruir valores integrales en segmentos en lugar de puntos; el primero controla mejor las oscilaciones, el segundo fuerza la conservación de masa.<br/>En el capítulo 2 estos métodos se aplican a algunos tests conocidos para averiguar su solidez.<br/>En el capítulo 3 estos métodos se aplican a la simulación de un diodo, y los resultados se comparan con resultados anteriores obtenidos por esquemas Runge-Kutta basados en diferencias finitas para aproximar las derivadas parciales.<br/>El capítulo 4 está dedicado a la construcción y simulación de modelos intermedios entre una ecuación cinética, con operador de colisión de tipo relajación, y su aproximación más grosera, ésta última siendo la ecuación del calor. Para obtener modelos intermedios, se busca un cierre de las ecuaciones de los momentos de orden cero y uno. Se proponen esquemas "asymptotic-preserving" para la ecuación cinética, que evitan la stiffness de la parte de advección a través de una descomposición de la función de distribución en su media más fluctuaciones. En cuanto a las clausuras de las ecuaciones de los momentos, se proponen esquemas de relajación para aislar las no-linealidades. Estos métodos son aplicados a un test conocido, el Su-Olson test.<br/>El último capítulo está dedicado a la simulación de un MOSFET (Metal Oxide Semiconductor Field Effect Transistor) 2D de dimensión nanométrica en el que los electrones se comportan como partículas en una dimensión y como ondas en las dimensiones confinadas. La descomposición en sub-bandas se realiza a través de una ecuación de Schrödinger 1D en estado estacionario. Las dimensiones, así como las sub-bandas, están acopladas por la ecuación de Poisson en la expresión de la densidad, y por el operador de colisión. Se propone un sólver microscópico para estados transitorios, basado en técnicas de splitting para las BTEs (una para cada nivel de energía), métodos de características para el transporte y una iteración de tipo Newton para resolver el problema acoplado Schrödinger-Poisson para el cálculo del campo de fuerza.

This thesis is dedicated to the development, application and test of numerical methods for the numerical simulation of problems arising from physics and electronic engineering. The main tool which is used all along the work is the Vlasov (transport) equation in the form of the Boltzmann Transport Equation (BTE) for the description of the transport and collisions of charged particles in plasmas and electronic devices: charge carriers are driven by a force field and scattered by other carriers or phonons (pseudo-particles giving an effective representation of the oscillating field produced by the vibrating ions).<br/>The BTE must be coupled to an equation or a system of equations for the computation of the force field: for simple structures the Poisson equation is used; for plasmas, where the magnetic phenomena cannot be neglected due to the high velocities of the particles, the Lorentz force is used, so the Maxwell equations have to be solved; for nanostructures, e.g. transistors with confined dimensions, the Poisson equation needs coupling with Schrödinger equation for the description of the quantum dimensions and the decomposition into subbands, or energy levels.<br/>Collisions mean the scattering the carriers suffer due to the interactions with other carriers or the fixed lattice, in form of phonons. All along the thesis several scattering operator are used: the simplest ones are linear relaxation-time operators; a model for the simulation of a semiconductor is studied in which collisions are taken into account with acoustic phonons, in the elastic approximation, and optical phonons.<br/>After the introduction, in the first chapter the most important numerical methods are developed: first of all a pointwise non-oscillatory interpolation method (PWENO) needed to avoid the simple Lagrange polynomial reconstruction, which increases the total variation when shocks appear: oscillations are part of the physics of the problem, but if the method adds spurious, non-physical oscillations, then the numerical result is meaningless, or it simply blows up. The second fundamental numerical method is the splitting technique: when solving a complicated problem, if we are able to subdivide it into sub-problem and solve them for separate, then we can reconstruct an approximation for the complete problem; this technique is used for both time splitting (separate transport from collisions) and dimensional splitting (split the phase space into either dimensions). The third fundamental instrument is the solver for linear advections: two methods are used, one based on pointwise following backwards the characteristics and another one based on reconstructing integral values along segments instead of point values; the first one controls better oscillations, the second one forces mass conservation.<br/>These methods are applied in chapter 2 to some well-known benchmark tests to control their robustness.<br/>In chapter 3 these methods are applied to the simulation of a diode, and the results compared to previous results obtained by Runge-Kutta schemes based on finite differences schemes for the approximation of the partial derivatives.<br/>Chapter 4 is dedicated to the construction and simulation of intermediate models between a kinetic equation, with relaxation-time collision operator, and its coarsest approximation, this one being the heat equations. In order to obtain intermediate models, the moment equations are closed at zeroth and first order. Asymptotic-preserving schemes are proposed for the kinetic equation, which avoid the stiffness of the advection part by decomposing the distribution function into its average plus fluctuations. As for the moment closures, relaxation schemes are proposed in order to confine the non-linearities in the right hand side. These methods are then applied to a known benchmark, the Su-Olson test.<br/>The last chapter is dedicated to the simulation of a nanoscaled 2D MOSFET (Metal Oxide Field Effect Transistor) in which electrons behave as particles in one dimension and as waves in the confined dimensions. The subband decomposition is realized through a stationary-state 1D Schrödinger equation. The dimensions as well as the subbands are coupled by the Poisson equation in the expression of the density and by the collision operator. A transient-state microscopic solver is proposed, based on splitting techniques for the BTE's (one for each energy level), characteristics methods for the transport and a Newton iteration for the solution of the coupled Schrödinger-Poisson system for computing the force field.