Existence of solutions and asymtptotic limits of the Euler-Poisson and the quantum drift-diffusion systems

My work concerns two different systems of equations used in the mathematical modeling of semiconductors and plasmas: the Euler-Poisson system and the quantum drift-diffusion system. The first is given by the Euler equations for the conservation of mass and momentum, with a Poisson equation for the electrostatic potential. The second one takes into account the physical effects due to the smallness of the devices (quantum effects). It is a simple extension of the classical drift-diffusion model which consists of two continuity equations for the charge densities, with a Poisson equation for the electrostatic potential. Using an asymptotic expansion method, we study (in the steady-state case for a potential flow) the limit to zero of the three physical parameters which arise in the Euler-Poisson system: the electron mass, the relaxation time and the Debye length. For each limit, we prove the existence and uniqueness of profiles to the asymptotic expansion and some error estimates. For a vanishing electron mass or a vanishing relaxation time, this method gives us a new approach in the convergence of the Euler-Poisson system to the incompressible Euler equations. For a vanishing Debye length (also called quasineutral limit), we obtain a new approach in the existence of solutions when boundary layers can appear (i.e. when no compatibility condition is assumed). Moreover, using an iterative method, and a finite volume scheme or a penalized mixed finite volume scheme, we numerically show the smallness condition on the electron mass needed in the existence of solutions to the system, condition which has already been shown in the literature. In the quantum drift-diffusion model for the transient bipolar case in one-space dimension, we show, by using a time discretization and energy estimates, the existence of solutions (for a general doping profile). We also prove rigorously the quasineutral limit (for a vanishing doping profile). Finally, using a new time discretization and an algorithmic construction of entropies, we prove some regularity properties for the solutions of the equation obtained in the quasineutral limit (for a vanishing pressure). This new regularity permits us to prove the positivity of solutions to this equation for at least times large enough.

External drivers of biogeochemical cycles in a tropical montane forest in Ecuador

Successful conservation of tropical montane forest, one of the most threatened ecosystems on earth, requires detailed knowledge of its biogeochemistry. Of particular interest is the response of the biogeochemical element cycles to external influences such as element deposition or climate change. Therefore the overall objective of my study was to contribute to improved understanding of role and functioning of the Andean tropical montane forest. In detail, my objectives were to determine (1) the role of long-range transported aerosols and their transport mechanisms, and (2) the role of short-term extreme climatic events for the element budget of Andean tropical forest. In a whole-catchment approach including three 8-13 ha microcatchments under tropical montane forest on the east-exposed slope of the eastern cordillera in the south Ecuadorian Andes at 1850-2200 m above sea level I monitored at least in weekly resolution the concentrations and fluxes of Ca, Mg, Na, K, NO3-N, NH4-N, DON, P, S, TOC, Mn, and Al in bulk deposition, throughfall, litter leachate, soil solution at the 0.15 and 0.3 m depths, and runoff between May 1998 and April 2003. I also used meteorological data from my study area collected by cooperating researchers and the Brazilian meteorological service (INPE), as well as remote sensing products of the North American and European space agencies NASA and ESA. My results show that (1) there was a strong interannual variation in deposition of Ca [4.4-29 kg ha-1 a-1], Mg [1.6-12], and K [9.8-30]) between 1998 and 2003. High deposition changed the Ca and Mg budgets of the catchments from loss to retention, suggesting that the additionally available Ca and Mg was used by the ecosystem. Increased base metal deposition was related to dust outbursts of the Sahara and an Amazonian precipitation pattern with trans-regional dry spells allowing for dust transport to the Andes. The increased base metal deposition coincided with a strong La Niña event in 1999/2000. There were also significantly elevated H+, N, and Mn depositions during the annual biomass burning period in the Amazon basin. Elevated H+ deposition during the biomass burning period caused elevated base metal loss from the canopy and the organic horizon and deteriorated already low base metal supply of the vegetation. Nitrogen was only retained during biomass burning but not during non-fire conditions when deposition was much smaller. Therefore biomass burning-related aerosol emissions in Amazonia seem large enough to substantially increase element deposition at the western rim of Amazonia. Particularly the related increase of acid deposition impoverishes already base-metal scarce ecosystems. As biomass burning is most intense during El Niño situations, a shortened ENSO cycle because of global warming likely enhances the acid deposition at my study forest. (2) Storm events causing near-surface water flow through C- and nutrient-rich topsoil during rainstorms were the major export pathway for C, N, Al, and Mn (contributing >50% to the total export of these elements). Near-surface flow also accounted for one third of total base metal export. This demonstrates that storm-event related near-surface flow markedly affects the cycling of many nutrients in steep tropical montane forests. Changes in the rainfall regime possibly associated with global climate change will therefore also change element export from the study forest. Element budgets of Andean tropical montane rain forest proved to be markedly affected by long-range transport of Saharan dust, biomass burning-related aerosols, or strong rainfalls during storm events. Thus, increased acid and nutrient deposition and the global climate change probably drive the tropical montane forest to another state with unknown consequences for its functions and biological diversity.

Microemulsion: prediction of the phase diagram with a modified Helfrich free energy

Microemulsions are thermodynamically stable, macroscopically homogeneous but microscopically heterogeneous, mixtures of water and oil stabilised by surfactant molecules. They have unique properties like ultralow interfacial tension, large interfacial area and the ability to solubilise other immiscible liquids. Depending on the temperature and concentration, non-ionic surfactants self assemble to micelles, flat lamellar, hexagonal and sponge like bicontinuous morphologies. Microemulsions have three different macroscopic phases (a) 1phase- microemulsion (isotropic), (b) 2phase-microemulsion coexisting with either expelled water or oil and (c) 3phase- microemulsion coexisting with expelled water and oil.rnrnOne of the most important fundamental questions in this field is the relation between the properties of the surfactant monolayer at water-oil interface and those of microemulsion. This monolayer forms an extended interface whose local curvature determines the structure of the microemulsion. The main part of my thesis deals with the quantitative measurements of the temperature induced phase transitions of water-oil-nonionic microemulsions and their interpretation using the temperature dependent spontaneous curvature [c0(T)] of the surfactant monolayer. In a 1phase- region, conservation of the components determines the droplet (domain) size (R) whereas in 2phase-region, it is determined by the temperature dependence of c0(T). The Helfrich bending free energy density includes the dependence of the droplet size on c0(T) as

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.