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).

What is the shape of a uniformly massive object that generates a gravitational potential equivalent to that of two equal point-masses? If the weight of each point-mass is sufficiently small compared to the distance between the points then the answer is a pair of balls of equal radius, one centered at each of the two points, but otherwise it is a certain domain of revolution about the axis passing through the two points. The existence and uniqueness of such a domain is known, but an explicit parameterization is known only in the plane where the region is referred to as a Neumann oval. We construct a four-dimensional “Neumann ovaloid”, solving explicitly this inverse potential problem.

In this paper we investigate staggered discontinuous Galerkin method for the Helmholtz equation with large wave number on general quadrilateral and polygonal meshes. The method is highly flexible by allowing rough grids such as the trapezoidal grids and highly distorted grids, and at the same time, is numerical flux free. Furthermore, it allows hanging nodes, which can be simply treated as additional vertices. By exploiting a modified duality argument, the stability and convergence can be proved under the condition that \kappa h is sufficiently small, where \kappa h is the wave number and \kappa h is the mesh size. Error estimates for both the scalar and vector variables in \kappa h norm are established. Several numerical experiments are tested to verify our theoretical results and to present the capability of our method for capturing singular solutions.

In this paper, we investigate the superconvergence property of the local discontinuous Galerkin (LDG) methods for solving one-dimensional linear time dependent fourth order problems. We prove that the error between the LDG solution and a particular projection of the exact solution, $\bar{e}_u$, achieves $(k+\frac32)$-th order superconvergence when polynomials of degree $k$ ($k\ge 1$) are used. Numerical experiments of $P^k$ polynomials, with $1\le k \le3$, are displayed to demonstrate the theoretical results, which show that the error $\bar{e}_u$ actually achieves $(k+2)$-th order superconvergence, indicating that the error bound for $\bar{e}_u$ obtained in this paper is sub-optimal. Initial-boundary value problems, nonlinear equations and solutions having singularities are numerically investigated to verify that the conclusions hold true for very general cases.

We present a new formulaiton of the incompressible Navier-Stokes equation in terms of an auxiliary field that differs from the velocity by a gauge transformation. The gauge freedom allows us to assign simple and specific boundary conditions for both the auxiliary field and the gauge field, thus eliminating the issue of pressure boundary condition in the usual primitive variable formulation. The resulting dynamic and kinematic equations can then be solved by standard methods for heat and Poisson equations. A normal mode analysis suggests that in contrast to the classical projection method, the gauge method does not suffer from the problem of numerical boundary layers. Thus the subtleties on the spatial discretization for the projection methods are removed. Consequently, the projection step in the gauge method can be accomplished by standard Poissin solves. We demonstrate the efficiency and accuracy of the gauge method by several numerical examples, including the flow past cylinder.