Mixed finite-element methods for elliptic convection-dominated problems arising in semiconductor physics

In this work we develop and analyze an adaptive numerical scheme for simulating a class of macroscopic semiconductor models. At first the numerical modelling of semiconductors is reviewed in order to classify the Energy-Transport models for semiconductors that are later simulated in 2D. In this class of models the flow of charged particles, that are negatively charged electrons and so-called holes, which are quasi-particles of positive charge, as well as their energy distributions are described by a coupled system of nonlinear partial differential equations. A considerable difficulty in simulating these convection-dominated equations is posed by the nonlinear coupling as well as due to the fact that the local phenomena such as "hot electron effects" are only partially assessable through the given data. The primary variables that are used in the simulations are the particle density and the particle energy density. The user of these simulations is mostly interested in the current flow through parts of the domain boundary - the contacts. The numerical method considered here utilizes mixed finite-elements as trial functions for the discrete solution. The continuous discretization of the normal fluxes is the most important property of this discretization from the users perspective. It will be proven that under certain assumptions on the triangulation the particle density remains positive in the iterative solution algorithm. Connected to this result an a priori error estimate for the discrete solution of linear convection-diffusion equations is derived. The local charge transport phenomena will be resolved by an adaptive algorithm, which is based on a posteriori error estimators. At that stage a comparison of different estimations is performed. Additionally a method to effectively estimate the error in local quantities derived from the solution, so-called "functional outputs", is developed by transferring the dual weighted residual method to mixed finite elements. For a model problem we present how this method can deliver promising results even when standard error estimator fail completely to reduce the error in an iterative mesh refinement process.

Convergence results for stochastic particle systems with social interaction

We consider stochastic individual-based models for social behaviour of groups of animals. In these models the trajectory of each animal is given by a stochastic differential equation with interaction. The social interaction is contained in the drift term of the SDE. We consider a global aggregation force and a short-range repulsion force. The repulsion range and strength gets rescaled with the number of animals N. We show that for N tending to infinity stochastic fluctuations disappear and a smoothed version of the empirical process converges uniformly towards the solution of a nonlinear, nonlocal partial differential equation of advection-reaction-diffusion type. The rescaling of the repulsion in the individual-based model implies that the corresponding term in the limit equation is local while the aggregation term is non-local. Moreover, we discuss the effect of a predator on the system and derive an analogous convergence result. The predator acts as an repulsive force. Different laws of motion for the predator are considered.

Efficiency limits for silicon/perovskite tandem solar cells: a theoretical model

The primary goal of this work is related to the extension of an analytic electro-optical model. It will be used to describe single-junction crystalline silicon solar cells and a silicon/perovskite tandem solar cell in the presence of light-trapping in order to calculate efficiency limits for such a device. In particular, our tandem system is composed by crystalline silicon and a perovskite structure material: metilammoniumleadtriiodide (MALI). Perovskite are among the most convenient materials for photovoltaics thanks to their reduced cost and increasing efficiencies. Solar cell efficiencies of devices using these materials increased from 3.8% in 2009 to a certified 20.1% in 2014 making this the fastest-advancing solar technology to date. Moreover, texturization increases the amount of light which can be absorbed through an active layer. Using Green’s formalism it is possible to calculate the photogeneration rate of a single-layer structure with Lambertian light trapping analytically. In this work we go further: we study the optical coupling between the two cells in our tandem system in order to calculate the photogeneration rate of the whole structure. We also model the electronic part of such a device by considering the perovskite top cell as an ideal diode and solving the drift-diffusion equation with appropriate boundary conditions for the silicon bottom cell. We have a four terminal structure, so our tandem system is totally unconstrained. Then we calculate the efficiency limits of our tandem including several recombination mechanisms such as Auger, SRH and surface recombination. We focus also on the dependence of the results on the band gap of the perovskite and we calculare an optimal band gap to optimize the tandem efficiency. The whole work has been continuously supported by a numerical validation of out analytic model against Silvaco ATLAS which solves drift-diffusion equations using a finite elements method. Our goal is to develop a simpler and cheaper, but accurate model to study such devices.

Invariant barriers to reactive front propagation in fluid flows

We present theory and experiments on the dynamics of reaction fronts in two-dimensional, vortex-dominated flows, for both time-independent and periodically driven cases. We find that the front propagation process is controlled by one-sided barriers that are either fixed in the laboratory frame (time-independent flows) or oscillate periodically (periodically driven flows). We call these barriers burning invariant manifolds (BIMs), since their role in front propagation is analogous to that of invariant manifolds in the transport and mixing of passive impurities under advection. Theoretically, the BIMs emerge from a dynamical systems approach when the advection-reaction-diffusion dynamics is recast as an ODE for front element dynamics. Experimentally, we measure the location of BIMs for several laboratory flows and confirm their role as barriers to front propagation.

Calcium waves driven by "sensitization" wave-fronts

OBJECTIVE: Cellular Ca(2+) waves are understood as reaction-diffusion systems sustained by Ca(2+)-induced Ca(2+) release (CICR) from Ca(2+) stores. Given the recently discovered sensitization of Ca(2+) release channels (ryanodine receptors; RyRs) of the sarcoplasmic reticulum (SR) by luminal SR Ca(2+), waves could also be driven by RyR sensitization, mediated by SR overloading via Ca(2+) pump (SERCA), acting in tandem with CICR. METHODS: Confocal imaging of the Ca(2+) indicator fluo-3 was combined with UV-flash photolysis of caged compounds and the whole-cell configuration of the patch clamp technique to carry out these experiments in isolated guinea pig ventricular cardiomyocytes. RESULTS: Upon sudden slowing of the SERCA in cardiomyocytes with a photoreleased inhibitor, waves indeed decelerated immediately. No secondary changes of Ca(2+) signaling or SR Ca(2+) content due to SERCA inhibition were observed in the short time-frame of these experiments. CONCLUSIONS: Our findings are consistent with Ca(2+) loading resulting in a zone of RyR 'sensitization' traveling within the SR, but inconsistent with CICR as the predominant mechanism driving the Ca(2+) waves. This alternative mode of RyR activation is essential to fully conceptualize cardiac arrhythmias triggered by spontaneous Ca(2+) release.

Predictive oncology: a review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth.

Empirical evidence and theoretical studies suggest that the phenotype, i.e., cellular- and molecular-scale dynamics, including proliferation rate and adhesiveness due to microenvironmental factors and gene expression that govern tumor growth and invasiveness, also determine gross tumor-scale morphology. It has been difficult to quantify the relative effect of these links on disease progression and prognosis using conventional clinical and experimental methods and observables. As a result, successful individualized treatment of highly malignant and invasive cancers, such as glioblastoma, via surgical resection and chemotherapy cannot be offered and outcomes are generally poor. What is needed is a deterministic, quantifiable method to enable understanding of the connections between phenotype and tumor morphology. Here, we critically assess advantages and disadvantages of recent computational modeling efforts (e.g., continuum, discrete, and cellular automata models) that have pursued this understanding. Based on this assessment, we review a multiscale, i.e., from the molecular to the gross tumor scale, mathematical and computational "first-principle" approach based on mass conservation and other physical laws, such as employed in reaction-diffusion systems. Model variables describe known characteristics of tumor behavior, and parameters and functional relationships across scales are informed from in vitro, in vivo and ex vivo biology. We review the feasibility of this methodology that, once coupled to tumor imaging and tumor biopsy or cell culture data, should enable prediction of tumor growth and therapy outcome through quantification of the relation between the underlying dynamics and morphological characteristics. In particular, morphologic stability analysis of this mathematical model reveals that tumor cell patterning at the tumor-host interface is regulated by cell proliferation, adhesion and other phenotypic characteristics: histopathology information of tumor boundary can be inputted to the mathematical model and used as a phenotype-diagnostic tool to predict collective and individual tumor cell invasion of surrounding tissue. This approach further provides a means to deterministically test effects of novel and hypothetical therapy strategies on tumor behavior.

Design of a Semi-analytical Method to Describe Plasma-Wall Interactions

Electric probes are objects immersed in the plasma with sharp boundaries which collect of emit charged particles. Consequently, the nearby plasma evolves under abrupt imposed and/or naturally emerging conditions. There could be localized currents, different time scales for plasma species evolution, charge separation and absorbing-emitting walls. The traditional numerical schemes based on differences often transform these disparate boundary conditions into computational singularities. This is the case of models using advection-diffusion differential equations with source-sink terms (also called Fokker-Planck equations). These equations are used in both, fluid and kinetic descriptions, to obtain the distribution functions or the density for each plasma species close to the boundaries. We present a resolution method grounded on an integral advancing scheme by using approximate Green's functions, also called short-time propagators. All the integrals, as a path integration process, are numerically calculated, what states a robust grid-free computational integral method, which is unconditionally stable for any time step. Hence, the sharp boundary conditions, as the current emission from a wall, can be treated during the short-time regime providing solutions that works as if they were known for each time step analytically. The form of the propagator (typically a multivariate Gaussian) is not unique and it can be adjusted during the advancing scheme to preserve the conserved quantities of the problem. The effects of the electric or magnetic fields can be incorporated into the iterative algorithm. The method allows smooth transitions of the evolving solutions even when abrupt discontinuities are present. In this work it is proposed a procedure to incorporate, for the very first time, the boundary conditions in the numerical integral scheme. This numerical scheme is applied to model the plasma bulk interaction with a charge-emitting electrode, dealing with fluid diffusion equations combined with Poisson equation self-consistently. It has been checked the stability of this computational method under any number of iterations, even for advancing in time electrons and ions having different time scales. This work establishes the basis to deal in future work with problems related to plasma thrusters or emissive probes in electromagnetic fields.

Twist grain boundaries in three-dimensional lamellar Turing structures

Steady spatial self-organization of three-dimensional chemical reaction-diffusion systems is discussed with the emphasis put on the possible defects that may alter the Turing patterns. It is shown that one of the stable defects of a three-dimensional lamellar Turing structure is a twist grain boundary embedding a Scherk minimal surface.

Microtubule self-organization is gravity-dependent

Although weightlessness is known to affect living cells, the manner by which this occurs is unknown. Some reaction-diffusion processes have been theoretically predicted as being gravity-dependent. Microtubules, a major constituent of the cellular cytoskeleton, self-organize in vitro by way of reaction-diffusion processes. To investigate how self-organization depends on gravity, microtubules were assembled under low gravity conditions produced during space flight. Contrary to the samples formed on an in-flight 1 × g centrifuge, the samples prepared in microgravity showed almost no self-organization and were locally disordered.

Self-organized criticality: sandpiles, singularities, and scaling.

We present an overview of the statistical mechanics of self-organized criticality. We focus on the successes and failures of hydrodynamic description of transport, which consists of singular diffusion equations. When this description applies, it can predict the scaling features associated with these systems. We also identify a hard driving regime where singular diffusion hydrodynamics fails due to fluctuations and give an explicit criterion for when this failure occurs.

Ecological chaos in the wake of invasion.

Irregularities in observed population densities have traditionally been attributed to discretization of the underlying dynamics. We propose an alternative explanation by demonstrating the evolution of spatiotemporal chaos in reaction-diffusion models for predator-prey interactions. The chaos is generated naturally in the wake of invasive waves of predators. We discuss in detail the mechanism by which the chaos is generated. By considering a mathematical caricature of the predator-prey models, we go on to explain the dynamical origin of the irregular behavior and to justify our assertion that the behavior we present is a genuine example of spatiotemporal chaos.

Comparison of Lagrangian Coherent Structures and Relative Dispersion in a Mixing Layer

Thesis (Master's)--University of Washington, 2016-06

Migration of additives from thermoplastic polymers

A homologous series of ultra-violet stabilisers containing 2-hydroxybenzophenone (HBP) moiety as a uv absorbing chromophore with varying alkyl chain lengths and sizes were prepared by known chemical synthesis. The strong absorbance of the HBP chromophore was utilized to evaluate the concentration of these stabilisers in low density polyethylene films and concentration of these stabilisers in low density polyethylene films and in relevant solvents by ultra-violet/visible spectroscopy. Intrinsic diffusion coefficients, equilibrium solubilities, volatilities from LDPE films and volatility of pure stabilisers were studied over a temperature range of 5-100oC. The effects of structure, molecular weight and temperature on the above parameters were investigated and the results were analysed on the basis of theoretical models published in the literature. It has been found that an increase in alkyl chain lengths does not change the diffusion coefficients to a significant level, while attachment of polar or branched alkyl groups change their value considerably. An Arrhenius type of relationship for the temperature dependence of diffusion coefficients seems to be valid only for a narrow temperature range, and therefore extrapolation of data from one temperature to another leads to a considerable error. The evidence showed that increase in additive solubility in the polymer is favoured by lower heat of fusions and melting points of additives. This implies the validity of simple regular solution theory to provide an adequate basis for understanding the solubility of additives in polymers The volubility of stabilisers from low density polyethylene films showed that of an additive from a polymer can be expressed in terms of a first-order kinetic equation. In addition the rate of loss of stabilisers was discussed in relation to its diffusion, solubility and volatility and found that all these factors may contribute to the additive loss, although one may be a rate determining factor. Stabiliser migration from LDPE into various solvents and food simulants was studied at temperatures 5, 23, 40 and 70oC; from the plots of rate of migration versus square root time, characteristic diffusion coefficients were obtained by using the solution of Fick's diffusion equations. It was shown that the rate of migration depends primarily on partition coefficients between solvent and the polymer of the additive and also on the swelling action of the contracting media. Characteristic diffusion coefficients were found to approach to intrinsic values in non swelling solvents, whereas in the case of highly swollen polymer samples, the former may be orders of magnitude greater than the latter.

Special Functions and Pathways for Problems in Astrophysics: An Essay in Honor of A.M. Mathai

The paper provides a review of A.M. Mathai's applications of the theory of special functions, particularly generalized hypergeometric functions, to problems in stellar physics and formation of structure in the Universe and to questions related to reaction, diffusion, and reaction-diffusion models. The essay also highlights Mathai's recent work on entropic, distributional, and differential pathways to basic concepts in statistical mechanics, making use of his earlier research results in information and statistical distribution theory. The results presented in the essay cover a period of time in Mathai's research from 1982 to 2008 and are all related to the thematic area of the gravitationally stabilized solar fusion reactor and fractional reaction-diffusion, taking into account concepts of non-extensive statistical mechanics. The time period referred to above coincides also with Mathai's exceptional contributions to the establishment and operation of the Centre for Mathematical Sciences, India, as well as the holding of the United Nations (UN)/European Space Agency (ESA)/National Aeronautics and Space Administration (NASA) of the United States/ Japanese Aerospace Exploration Agency (JAXA) Workshops on basic space science and the International Heliophysical Year 2007, around the world. Professor Mathai's contributions to the latter, since 1991, are a testimony for his social con-science applied to international scientific activity.

Noise-induced standing waves in oscillatory systems with time-delayed feedback

In oscillatory reaction-diffusion systems, time-delay feedback can lead to the instability of uniform oscillations with respect to formation of standing waves. Here, we investigate how the presence of additive, Gaussian white noise can induce the appearance of standing waves. Combining analytical solutions of the model with spatio-temporal simulations, we find that noise can promote standing waves in regimes where the deterministic uniform oscillatory modes are stabilized. As the deterministic phase boundary is approached, the spatio-temporal correlations become stronger, such that even small noise can induce standing waves in this parameter regime. With larger noise strengths, standing waves could be induced at finite distances from the (deterministic) phase boundary. The overall dynamics is defined through the interplay of noisy forcing with the inherent reaction-diffusion dynamics.