The purpose of this work is to incorporate numerically, in a discontinuous Galerkin (DG) solver of a Boltzmann-Poisson model for hot electron transport, an electronic conduction band whose values are obtained by the spherical averaging of the full band structure given by a local empirical pseudopotential method (EPM) around a local minimum of the conduction band for silicon, as a midpoint between a radial band model and an anisotropic full band, in order to provide a more accurate physical description of the electron group velocity and conduction energy band structure in a semiconductor. This gives a better quantitative description of the transport and collision phenomena that fundamentally define the behaviour of the Boltzmann - Poisson model for electron transport used in this work. The numerical values of the derivatives of this conduction energy band, needed for the description of the electron group velocity, are obtained by means of a cubic spline interpolation. The EPM-Boltzmann-Poisson transport with this spherically averaged EPM calculated energy surface is numerically simulated and compared to the output of traditional analytic band models such as the parabolic and Kane bands, numerically implemented too, for the case of 1D $n^+-n-n^+$ silicon diodes with $400nm$ and $50nm$ channels. Quantitative differences are observed in the kinetic moments related to the conduction energy band used, such as mean velocity, average energy, and electric current (momentum).

In this paper, we develop and analyze the Runge-Kutta discontinuous Galerkin (RKDG) method to solve weakly coupled hyperbolic multi-domain problems. Such problems involve transfer type boundary conditions with discontinuous fluxes between different domains, calling for special techniques to prove stability of the RKDG methods. We prove both stability and error estimates for our RKDG methods on simple models, and then apply them to a biological cell proliferation model \cite{EMSC}. Numerical results are provided to illustrate the good behavior of our RKDG methods.

The main purpose of this paper is to give stability analysis and error estimates of the local discontinuous Galerkin (LDG) methods coupled with three specific implicit-explicit (IMEX) Runge-Kutta time discretization methods up to third order accuracy, for solving one-dimensional time-dependent linear fourth order partial differential equations. In the time discretization, all the lower order derivative terms are treated explicitly and the fourth order derivative term is treated implicitly. By the aid of energy analysis, we show that the IMEX-LDG schemes are unconditionally energy stable, in the sense that the time step $\dt$ is only required to be upper-bounded by a constant which is independent of the mesh size $h$. The optimal error estimate is also derived by the aid of the elliptic projection and the adjoint argument. Numerical experiments are given to verify that the corresponding IMEX-LDG schemes can achieve optimal error accuracy.

In this paper, we propose a simple bound-preserving sweeping procedure for conservative numerical approximations. Conservative schemes are of importance in many applications, yet for high order methods, the numerical solutions do not necessarily satisfy maximum principle. This paper constructs a simple sweeping algorithm to enforce the bound of the solutions. It has a very general framework acting as a postprocessing step accommodating many point-based or cell average-based discretizations. The method is proven to preserve the bounds of the numerical solution while conserving a prescribed quantity designated as a weighted average of values from all points. The technique is demonstrated to work well with a spectral method, high order finite difference and finite volume methods for scalar conservation laws and incompressible flows. Extensive numerical tests in 1D and 2D are provided to verify the accuracy of the sweeping procedure.

Simulation of electromagnetic wave propagation in metamaterials leads to more complicated time domain Maxwell's equations than the standard Maxwell's equations in free space. In this paper, we develop and analyze a non-dissipative discontinuous Galerkin (DG) method for solving the Maxwell's equations in Drude metamaterials. Previous discontinuous Galerkin methods in the literature for electromagnetic wave propagation in metamaterials were either non-dissipative but sub-optimal, or dissipative and optimal. Our method uses a different and simple choice of numerical fluxes, achieving provable non-dissipative stability and optimal error estimates simultaneously. We prove the stability and optimal error estimates for both semi- and fully discrete DG schemes, with the leap-frog time discretization for the fully discrete case. Numerical results are given to demonstrate that the DG method can solve metamaterial Maxwell's equations effectively.