On four numerical schemes for a unipolar degenerate drift-diffusion model

International audienceWe consider a unipolar degenerate drift-diffusion system where the relation between the concentration of the charged species c and the chemical potential h is $h(c) = log(c/(1−c))$. For four different finite volume schemes based on four different formulations of the fluxes of the problem, we discuss stability and existence results. For two of them, we report a convergence proof. Numerical experiments illustrate the behaviour of the different schemes

Uncertainty quantification for kinetic models in socio-economic and life sciences

Kinetic equations play a major rule in modeling large systems of interacting particles. Recently the legacy of classical kinetic theory found novel applications in socio-economic and life sciences, where processes characterized by large groups of agents exhibit spontaneous emergence of social structures. Well-known examples are the formation of clusters in opinion dynamics, the appearance of inequalities in wealth distributions, flocking and milling behaviors in swarming models, synchronization phenomena in biological systems and lane formation in pedestrian traffic. The construction of kinetic models describing the above processes, however, has to face the difficulty of the lack of fundamental principles since physical forces are replaced by empirical social forces. These empirical forces are typically constructed with the aim to reproduce qualitatively the observed system behaviors, like the emergence of social structures, and are at best known in terms of statistical information of the modeling parameters. For this reason the presence of random inputs characterizing the parameters uncertainty should be considered as an essential feature in the modeling process. In this survey we introduce several examples of such kinetic models, that are mathematically described by nonlinear Vlasov and Fokker--Planck equations, and present different numerical approaches for uncertainty quantification which preserve the main features of the kinetic solution.Comment: To appear in "Uncertainty Quantification for Hyperbolic and Kinetic Equations

A Lagrangian scheme for the solution of nonlinear diffusion equations using moving simplex meshes

A Lagrangian numerical scheme for solving nonlinear degenerate Fokker{Planck equations in space dimensions d>2 is presented. It applies to a large class of nonlinear diffusion equations, whose dynamics are driven by internal energies and given external potentials, e.g. the porous medium equation and the fast diffusion equation. The key ingredient in our approach is the gradient ow structure of the dynamics. For discretization of the Lagrangian map, we use a finite subspace of linear maps in space and a variational form of the implicit Euler method in time. Thanks to that time discretisation, the fully discrete solution inherits energy estimates from the original gradient ow, and these lead to weak compactness of the trajectories in the continuous limit. Consistency is analyzed in the planar situation, d = 2. A variety of numerical experiments for the porous medium equation indicates that the scheme is well-adapted to track the growth of the solution's support

Multi-dimensional modeling and simulation of semiconductor nanophotonic devices

Self-consistent modeling and multi-dimensional simulation of semiconductor nanophotonic devices is an important tool in the development of future integrated light sources and quantum devices. Simulations can guide important technological decisions by revealing performance bottlenecks in new device concepts, contribute to their understanding and help to theoretically explore their optimization potential. The efficient implementation of multi-dimensional numerical simulations for computer-aided design tasks requires sophisticated numerical methods and modeling techniques. We review recent advances in device-scale modeling of quantum dot based single-photon sources and laser diodes by self-consistently coupling the optical Maxwell equations with semiclassical carrier transport models using semi-classical and fully quantum mechanical descriptions of the optically active region, respectively. For the simulation of realistic devices with complex, multi-dimensional geometries, we have developed a novel hp-adaptive finite element approach for the optical Maxwell equations, using mixed meshes adapted to the multi-scale properties of the photonic structures. For electrically driven devices, we introduced novel discretization and parameter-embedding techniques to solve the drift-diffusion system for strongly degenerate semiconductors at cryogenic temperature. Our methodical advances are demonstrated on various applications, including vertical-cavity surface-emitting lasers, grating couplers and single-photon sources

On the L2L^2 Stability of Finite Volumes for Stationary First Order Systems

International audienceThe aim of this paper is two-folds. Firstly we study first order stationary systems of PDEs of the form $\sum_k A_k\partial_k U + KU = 0$ with $K_\eta tr0$ on $R^d$ . We prove that the classical assumption $K > 0$ is not necessary for the well-posedness of the system and is violated in the particular case of the first order Poisson problem. Secondly we prove the $L^2$ stability of the finite volume discretisations provided the term KU is appropriately discretised on faces. Our result relies on a discrete Gagliardo-Nirenberg-Sobolev inequality to be submitted