Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

Observations and Modeling of Tropical Planetary Atmospheres

This thesis is a comprised of three different projects within the topic of tropical atmospheric dynamics. First, I analyze observations of thermal radiation from Saturn’s atmosphere and from them, determine the latitudinal distribution of ammonia vapor near the 1.5-bar pressure level. The most prominent feature of the observations is the high brightness temperature of Saturn’s subtropical latitudes on either side of the equator. After comparing the observations to a microwave radiative transfer model, I find that these subtropical bands require very low ammonia relative humidity below the ammonia cloud layer in order to achieve the high brightness temperatures observed. We suggest that these bright subtropical bands represent dry zones created by a meridionally overturning circulation.

Second, I use a dry atmospheric general circulation model to study equatorial superrotation in terrestrial atmospheres. A wide range of atmospheres are simulated by varying three parameters: the pole-equator radiative equilibrium temperature contrast, the convective lapse rate, and the planetary rotation rate. A scaling theory is developed that establishes conditions under which superrotation occurs in terrestrial atmospheres. The scaling arguments show that superrotation is favored when the off-equatorial baroclinicity and planetary rotation rates are low. Similarly, superrotation is favored when the convective heating strengthens, which may account for the superrotation seen in extreme global-warming simulations.

Third, I use a moist slab-ocean general circulation model to study the impact of a zonally-symmetric continent on the distribution of monsoonal precipitation. I show that adding a hemispheric asymmetry in surface heat capacity is sufficient to cause symmetry breaking in both the spatial and temporal distribution of precipitation. This spatial symmetry breaking can be understood from a large-scale energetic perspective, while the temporal symmetry breaking requires consideration of the dynamical response to the heat capacity asymmetry and the seasonal cycle of insolation. Interestingly, the idealized monsoonal precipitation bears resemblance to precipitation in the Indian monsoon sector, suggesting that this work may provide insight into the causes of the temporally asymmetric distribution of precipitation over southeast Asia.

A model for energy and morphology of crystalline grain boundaries with arbitrary geometric character

It has been well-established that interfaces in crystalline materials are key players in the mechanics of a variety of mesoscopic processes such as solidification, recrystallization, grain boundary migration, and severe plastic deformation. In particular, interfaces with complex morphologies have been observed to play a crucial role in many micromechanical phenomena such as grain boundary migration, stability, and twinning. Interfaces are a unique type of material defect in that they demonstrate a breadth of behavior and characteristics eluding simplified descriptions. Indeed, modeling the complex and diverse behavior of interfaces is still an active area of research, and to the author's knowledge there are as yet no predictive models for the energy and morphology of interfaces with arbitrary character. The aim of this thesis is to develop a novel model for interface energy and morphology that i) provides accurate results (especially regarding "energy cusp" locations) for interfaces with arbitrary character, ii) depends on a small set of material parameters, and iii) is fast enough to incorporate into large scale simulations.

In the first half of the work, a model for planar, immiscible grain boundary is formulated. By building on the assumption that anisotropic grain boundary energetics are dominated by geometry and crystallography, a construction on lattice density functions (referred to as "covariance") is introduced that provides a geometric measure of the order of an interface. Covariance forms the basis for a fully general model of the energy of a planar interface, and it is demonstrated by comparison with a wide selection of molecular dynamics energy data for FCC and BCC tilt and twist boundaries that the model accurately reproduces the energy landscape using only three material parameters. It is observed that the planar constraint on the model is, in some cases, over-restrictive; this motivates an extension of the model.

In the second half of the work, the theory of faceting in interfaces is developed and applied to the planar interface model for grain boundaries. Building on previous work in mathematics and materials science, an algorithm is formulated that returns the minimal possible energy attainable by relaxation and the corresponding relaxed morphology for a given planar energy model. It is shown that the relaxation significantly improves the energy results of the planar covariance model for FCC and BCC tilt and twist boundaries. The ability of the model to accurately predict faceting patterns is demonstrated by comparison to molecular dynamics energy data and experimental morphological observation for asymmetric tilt grain boundaries. It is also demonstrated that by varying the temperature in the planar covariance model, it is possible to reproduce a priori the experimentally observed effects of temperature on facet formation.

Finally, the range and scope of the covariance and relaxation models, having been demonstrated by means of extensive MD and experimental comparison, future applications and implementations of the model are explored.