Large time behavior framework for the time-increasing weak solutions of bipolar hydrodynamic model of semiconductors

Abstract: In this paper, we consider an isentropic Euler-Poisson equations for the bipolar hydrodynamic model of semiconductor devices, which has a non-flat doping profile and insulating boundary conditions. Using a technical energy method and an entropy dissipation estimate, we present a framework for the large time behavior of time-increasing weak entropy solutions. It is shown that the weak solutions converge to the stationary solutions in L2 norm with exponential decay rate. No regularity and smallness conditions are assumed.


Introduction
In this paper, isentropic Euler-Poisson equations for the bipolar hydrodynamic model of semiconductor devices are considered.This model is as follows here n 1 , n 2 , J 1 , J 2 and E are the unknown functions of the space variable x ∈ [0, 1] and time variable t ≥ 0, representing the electron density, the hole density, the electron current density, the hole current density and the electric field respectively.The function D(x) > 0, usually called the doping profile, stands for the density of impurities in semiconductor devices.In this paper, we assume the doping profile D(x) satisfies (1.2) For the system (1.1), the initial-boundary conditions are described by n i (x, 0) = n i0 (x) ≥ 0, J i (x, 0) = J i0 (x), with the compatibility condition When n 2 = J 2 = 0 in (1.1), the bipolar model turns into the unipolar one, that is (1.5) Recently, many efforts are made on the systems (1.1) and (1.5) to considering the large time behavior of their weak entropy solutions.With the smallness assumption on the amplitude of background electron current, [1] first proved the uniformly bounded density weak entropy solutions of the unipolar hydrodynamic model (1.5), decay exponentially to the stationary solutions.[3] considered a similar problem on the bipolar model with a non-flat doping profile.However, the uniform bounded condition 0 ≤ n i (x, t) ≤ C 0 (1.6) in [1] ( [3])is stiff and still be open although it seems natural from physical point of view.For example, the L ∞ bounds obtained in [2,4,7] grow with time.In this paper, instead of proving the hard bone (1.6), we will give a large time behavior framework for density time-increasing entropy solutions to the bipolar hydrodynamic model (1.1) − (1.3).The related work on unipolar model, we can see the reference [6].We make some preparation work before to introduce the primary result.The vector function (n 1 , n 2 , J 1 , J 2 , E) is a weak solution of problem (1.1) − (1.4), if it satisfies the equation (1.1) in the distributional sense, verifies the restriction (1.3) and (1.4).Furthermore, a weak solution of system (1.1) − (1.4) is called an entropy solution if it satisfies the entropy inequality in the sense of distribution.And (η e , q e ) are mechanical entropy-entropy flux pair satisfying (1.8)

AIMS Mathematics
Volume 2, Issue 1, 102-110 104 The corresponding stationary system of problem (1.1) with the boundary condition (1.10) In reference [5], the author give the following existence and uniqueness Theorem, that is theorem A Problem (1.9) − (1.10) has an unique stationary solution

results
This following Theorem is main result of this paper.Theorem 1 (Large time behavior framework).
(N 1 , N 2 , E)(x, t) be the unique stationary smooth solution.If ) then there exist positive constants T (α), C, and C such that for any t > T (α).
With less regularity of the L ∞ entropy solutions, we can only obtain zero-order estimates.To get the exponential time decay estimate between the entropy solution and the corresponding stationary solution, we need explore the entropy dissipation.

Large time behavior framework of time-increasing entropy solutions
In this part, we will prove the large time behavior framework for the L ∞ entropy solutions, in which the bounds of densities may increase with time, that is Theorem 1. Specifically speaking, for any global entropy solutions of (1.1) − (1.4) with the densities satisfy (2.1), we get an exponential decay rate for the electric field and the relative entropy between the entropy solution and the stationary solution.To this purpose, we introduce new variables Naturally, y i (i = 1, 2) is absolutely continuous in x for a.e t > 0.Moreover, we have From (1.1) and the corresponding stationary equation, we get y i (i = 1, 2) admits the equations Multiplying y i with (3.3) integrating over the spatial domain (0,1) and then adding the results together for i = 1, 2, we have Noticing then we have 3) satisfies, we have for some small positive constant δ and δ = δ 1 − δ < 1.To see this, let ε = 2δ > 0, we have (3.10) (3.11) Now we explore the entropy dissipation estimate.To this end, we introduce the relative entropyentropy flux pair to make full use of the entropy inequality.
The relative entropy-entropy flux are: where , η e and q e are the entropy-entropy flux pair defined in (1.8).
Using the entropy inequality, we have the following estimates on the relative entropy-entropy flux pair (η * , q * ): that is d dt where , 1}. Since then there exists positive constant C 2 such that (3.18) turns into for some positive constant C 3 > 0. * since we consider the large time behavior, without loss of generality, we always assume t > t * .
On the other hand, noticing for some constant C 5 > 0. Thus, we prove Theorem 1.

Remark on the assumption (2.3)
The assumption (2.It is worthy to point out that (4.1) indicates that n i is bounded with respect to x and t.While the assumption (2.3) permits the bounds of n i (i = 1, 2) grow with time.