A BLOCK-CENTERED UPWIND APPROXIMATION OF THE SEMICONDUCTOR DEVICE PROBLEM ON A DYNAMICALLY CHANGING MESH

The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.

An Upwind Mixed Volume Element-Fractional Step Method on a Changing Mesh for Compressible Contamination Treatment from Nuclear Waste

In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convection-diffusion system of four partial differential equations:a parabolic equation for the pressure,two convection-diffusion equations for the concentrations of brine and radionuclide and a heat conduction equation for the temperature.The pressure appears within the concentration equations and heat conduction equation,and the Darcy velocity controls the concentrations and the temperature.The pressure is solved by the conservative mixed volume element method,and the order of the accuracy is improved by the Darcy velocity.The concentration of brine and temperature are computed by the upwind mixed volume element method on a changing mesh,where the diffusion is discretized by a mixed volume element and the convection is treated by an upwind scheme.The composite method can solve the convection-dominated diffusion problems well because it eliminates numerical dispersion and nonphysical oscillation and has high order computational accuracy.The mixed volume element has the local conservation of mass and energy,and it can obtain the brine and temperature and their adjoint vector functions simultaneously.The conservation nature plays an important role in numerical simulation of underground fluid.The concentrations of radionuclide factors are solved by the method of upwind fractional step difference and the computational work is decreased by decomposing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.By the theory and technique of a priori estimates of differential equations,we derive an optimal order result in L^(2) norm.Numerical examples are given to show the effectiveness and practicability and the composite method is testified as a powerful tool to solve the well-known actual problem.

AConservative Upwind Approximation on Block-Centered Difference for Chemical Oil Recovery Displacement Problem

A kind of conservative upwind method is discussed for chemical oil recovery displacement in porous media.The mathematical model is formulated by a nonlinear convection-diffusion system dependent on the pressure,Darcy velocity,concentration and saturations.The flow equation is solved by a conservative block-centered method,and the pressure and Darcy velocity are obtained at the same time.The concentration and saturations are determined by convection-dominated diffusion equations,so an upwind approximation is adopted to eliminate numerical dispersion and nonphysical oscillation.Block-centered method is conservative locally.An upwind method with block-centered difference is used for computing the concentration.The saturations of different components are solved by the method of upwind fractional step difference,and the computational work is shortened significantly by dividing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.Using the variation discussion,energy estimates,the method of duality,and the theory of a priori estimates,we complete numerical analysis.Finally,numerical tests are given for showing the computational accuracy,efficiency and practicability of our approach.