For a bipolar hydrodynamic model of semiconductors in the form of EulerPoisson equations with Dirichlet or Neumann boundary conditions, in this paper we first heuristically analyze the most probable asymptotic profile (the so-called diffusion waves) and then prove this long-time behavior rigorously. For this, we construct correction functions to show the convergence of the original solution to the diffusion wave with optimal convergence rates by the energy method. Moreover, in the case with Dirichlet boundary condition, when the initial perturbation is in some weighted L^1 space, a faster and optimal convergence rate is also given.

In this paper, high order central finite difference schemes in a finite interval are analyzed for the diffusion equation. Boundary conditions of the initial-boundary value problem (IBVP) are treated by the simplified inverse Lax-Wendroff (SILW) procedure. For the fully discrete case, a third order explicit Runge-Kutta method is used as an example for the analysis. Stability is analyzed by both the GKS (Gustafsson, Kreiss and Sundstr\"om) theory and the eigenvalue visualization method on both semi-discrete and fully discrete schemes. The two different analysis techniques yield consistent results. Numerical tests are performed to demonstrate and validate the analysis results.

In earlier work, Maire developed a class of cell-centered Lagrangian schemes for solving Euler equations of compressible gas dynamics in cylindrical coordinates. These schemes use a node-based discretization of the numerical fluxes. The control volume version has several distinguished properties, including the conservation of mass, momentum and total energy and compatibility with the geometric conservation law (GCL). However it also has a limitation in that it cannot preserve spherical symmetry for one-dimensional spherical flow. An alternative is also given to use the first order area-weighted approach which can ensure spherical symmetry, at the price of sacrificing conservation of momentum. In this paper, we apply the methodology proposed in our recent work to the first order control volume scheme of Maire to obtain the spherical symmetry property. The modified scheme can preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid, and meanwhile it maintains its original good properties such as conservation and GCL. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme in terms of symmetry, non-oscillation and robustness properties.

No abstract uploaded!

No abstract uploaded!