Maximum principle preserving high order schemes for convection-dominated diffusion equations

The maximum principle is an important property of solutions to PDE. Correspondingly, it's of great interest for people to design a high order numerical scheme solving PDE with this property maintained. In this thesis, our particular interest is solving convection-dominated diffusion equation. We first review a nonconventional maximum principle preserving(MPP) high order finite volume(FV) WENO scheme, and then propose a new parametrized MPP high order finite difference(FD) WENO framework, which is generalized from the one solving hyperbolic conservation laws. A formal analysis is presented to show that a third order finite difference scheme with this parametrized MPP flux limiters maintains the third order accuracy without extra CFL constraint when the low order monotone flux is chosen appropriately. Numerical tests in both one and two dimensional cases are performed on the simulation of the incompressible Navier-Stokes equations in vorticity stream-function formulation and several other problems to show the effectiveness of the proposed method.

Hamilton decompositions of 6-regular abelian Cayley graphs

In 1969, Lovasz asked whether every connected, vertex-transitive graph has a Hamilton path. This question has generated a considerable amount of interest, yet remains vastly open. To date, there exist no known connected, vertex-transitive graph that does not possess a Hamilton path. For the Cayley graphs, a subclass of vertex-transitive graphs, the following conjecture was made: Weak Lovász Conjecture: Every nontrivial, finite, connected Cayley graph is hamiltonian. The Chen-Quimpo Theorem proves that Cayley graphs on abelian groups flourish with Hamilton cycles, thus prompting Alspach to make the following conjecture: Alspach Conjecture: Every 2k-regular, connected Cayley graph on a finite abelian group has a Hamilton decomposition. Alspach’s conjecture is true for k = 1 and 2, but even the case k = 3 is still open. It is this case that this thesis addresses. Chapters 1–3 give introductory material and past work on the conjecture. Chapter 3 investigates the relationship between 6-regular Cayley graphs and associated quotient graphs. A proof of Alspach’s conjecture is given for the odd order case when k = 3. Chapter 4 provides a proof of the conjecture for even order graphs with 3-element connection sets that have an element generating a subgroup of index 2, and having a linear dependency among the other generators. Chapter 5 shows that if Γ = Cay(A, {s1, s2, s3}) is a connected, 6-regular, abelian Cayley graph of even order, and for some1 ≤ i ≤ 3, Δi = Cay(A/(si), {sj1 , sj2}) is 4-regular, and Δi ≄ Cay(ℤ3, {1, 1}), then Γ has a Hamilton decomposition. Alternatively stated, if Γ = Cay(A, S) is a connected, 6-regular, abelian Cayley graph of even order, then Γ has a Hamilton decomposition if S has no involutions, and for some s ∈ S, Cay(A/(s), S) is 4-regular, and of order at least 4. Finally, the Appendices give computational data resulting from C and MAGMA programs used to generate Hamilton decompositions of certain non-isomorphic Cayley graphs on low order abelian groups.

Statistics of the received power for free space optical channels

Free space optical (FSO) communication links can experience extreme signal degradation due to atmospheric turbulence induced spatial and temporal irradiance fuctuations (scintillation) in the laser wavefront. In addition, turbulence can cause the laser beam centroid to wander resulting in power fading, and sometimes complete loss of the signal. Spreading of the laser beam and jitter are also artifacts of atmospheric turbulence. To accurately predict the signal fading that occurs in a laser communication system and to get a true picture of how this affects crucial performance parameters like bit error rate (BER) it is important to analyze the probability density function (PDF) of the integrated irradiance fuctuations at the receiver. In addition, it is desirable to find a theoretical distribution that accurately models these ?uctuations under all propagation conditions. The PDF of integrated irradiance fuctuations is calculated from numerical wave-optic simulations of a laser after propagating through atmospheric turbulence to investigate the evolution of the distribution as the aperture diameter is increased. The simulation data distribution is compared to theoretical gamma-gamma and lognormal PDF models under a variety of scintillation regimes from weak to very strong. Our results show that the gamma-gamma PDF provides a good fit to the simulated data distribution for all aperture sizes studied from weak through moderate scintillation. In strong scintillation, the gamma-gamma PDF is a better fit to the distribution for point-like apertures and the lognormal PDF is a better fit for apertures the size of the atmospheric spatial coherence radius ρ0 or larger. In addition, the PDF of received power from a Gaussian laser beam, which has been adaptively compensated at the transmitter before propagation to the receiver of a FSO link in the moderate scintillation regime is investigated. The complexity of the adaptive optics (AO) system is increased in order to investigate the changes in the distribution of the received power and how this affects the BER. For the 10 km link, due to the non-reciprocal nature of the propagation path the optimal beam to transmit is unknown. These results show that a low-order level of complexity in the AO provides a better estimate for the optimal beam to transmit than a higher order for non-reciprocal paths. For the 20 km link distance it was found that, although minimal, all AO complexity levels provided an equivalent improvement in BER and that no AO complexity provided the correction needed for the optimal beam to transmit. Finally, the temporal power spectral density of received power from a FSO communication link is investigated. Simulated and experimental results for the coherence time calculated from the temporal correlation function are presented. Results for both simulation and experimental data show that the coherence time increases as the receiving aperture diameter increases. For finite apertures the coherence time increases as the communication link distance is increased. We conjecture that this is due to the increasing speckle size within the pupil plane of the receiving aperture for an increasing link distance.

Will Whole-Tree Harvest of Jack Pine (Pinus banksiana) deplete soil nutrients in low-productivity sand soils?

In 2009 and 2010 a study was conducted on the Hiawatha National Forest (HNF) to determine if whole-tree harvest (WTH) of jack pine would deplete the soil nutrients in the very coarse-textured Rubicon soil. WTH is restricted on Rubicon sand in order to preserve the soil fertility, but the increasing construction of biomass-fueled power plants is expected to increase the demand for forest biomass. The specific objectives of this study were to estimate biomass and nutrient content of above- and below-ground tree components in mature jack pine (Pinus banksiana) stands growing on a coarse-textured, low-productivity soil, determine pools of total C and N and exchangeable soil cations in Rubicon sand, and to compare the possible impacts of conventional stem-only harvest (CH) and WTH on soil nutrient pools and the implications for productivity of subsequent rotations. Four even-aged jack pine stands on Rubicon soil were studied. Allometric equations were used to estimate above-ground biomass and nutrients, and soil samples from each stand were taken for physical and chemical analysis. Results indicate that WTH will result in cation deficits in all stands, with exceptionally large Ca deficits occurring in two stands. Where a deficit does not occur, the cation surplus is small and, chemical weathering and atmospheric deposition is not anticipated to replace the removed cations. CH will result in a surplus of cations, and will likely not result in productivity declines during the first rotation. However even under CH, the surplus is small, and chemical weathering and atmospheric deposition will not supply enough cations for the second rotation.

Water transport in complex, non-wetting porous layers with applications to water management in low temperature fuel cells

An experimental setup was designed to visualize water percolation inside the porous transport layer, PTL, of proton exchange membrane, PEM, fuel cells and identify the relevant characterization parameters. In parallel with the observation of the water movement, the injection pressure (pressure required to transport water through the PTL) was measured. A new scaling for the drainage in porous media has been proposed based on the ratio between the input and the dissipated energies during percolation. A proportional dependency was obtained between the energy ratio and a non-dimensional time and this relationship is not dependent on the flow regime; stable displacement or capillary fingering. Experimental results show that for different PTL samples (from different manufacturers) the proportionality is different. The identification of this proportionality allows a unique characterization of PTLs with respect to water transport. This scaling has relevance in porous media flows ranging far beyond fuel cells. In parallel with the experimental analysis, a two-dimensional numerical model was developed in order to simulate the phenomena observed in the experiments. The stochastic nature of the pore size distribution, the role of the PTL wettability and morphology properties on the water transport were analyzed. The effect of a second porous layer placed between the porous transport layer and the catalyst layer called microporous layer, MPL, was also studied. It was found that the presence of the MPL significantly reduced the water content on the PTL by enhancing fingering formation. Moreover, the presence of small defects (cracks) within the MPL was shown to enhance water management. Finally, a corroboration of the numerical simulation was carried out. A threedimensional version of the network model was developed mimicking the experimental conditions. The morphology and wettability of the PTL are tuned to the experiment data by using the new energy scaling of drainage in porous media. Once the fit between numerical and experimental data is obtained, the computational PTL structure can be used in different types of simulations where the conditions are representative of the fuel cell operating conditions.

An Ecological and Chemical Assessment of Copper Contaminated Streams in the Keweenaw Peninsula

The copper mining boom in Michigan's Upper Peninsula ended in the mid-1960s, but the historical mining still affects the region to this day. Earlier studies conducted in the Keweenaw have shown that trace metals in the sediments negatively affect benthic macroinvertebrate populations. However, because the concentrations of trace metals that are observed to be toxic often differ significantly between the laboratory and the environment, a better method for determining toxic levels of trace metals in the natural environment is desirable in order to establish surface water quality guidelines that effectively protect aquatic life. There were four research objectives for this research project. First, to determine if trace-level concentrations of copper can result in detectable ecological impacts even in the presence of high dissolved organic carbon (DOC). Second, to determine if there is a "safe" concentration of total dissolved copper below which there is little to no ecological impairment. Third, to establish which streams in the Keweenaw Peninsula have been most impacted by elevated levels of total dissolved copper. Fourth, to use this information to evaluate revisions to the water quality criterion for copper that were recently proposed by the Michigan Department of Environmental Quality (MDEQ). In order to collect water quality and macroinvertebrate data, two sampling surveys of approximately 50 streams were completed in the spring and summer of 2012. Our findings demonstrate that negative ecological impacts can be detected even in the presence of high concentrations of DOC. The majority of surveyed streams showed evidence of total dissolved copper concentrations that were elevated above background levels. Our findings suggest that there are detectable negative impacts below the current water quality standard for copper in many Keweenaw streams. The diversity of benthic macroinvertebrates and the number of species present has been reduced as a result of exposure to copper. Additionally, the multimetric approach used by MDEQ is unable to detect copper impairment in local streams due to the use of several insensitive metrics. The proposed changes to the copper criterion would increase the amount of total dissolved copper allowable despite the fact that approximately 25% of streams sampled have aquatic chemistries that would leave them vulnerable to high levels of copper ions.

ELECTRON TRANSPORT IN LOW-DIMENSIONAL NANOSTRUCTURES - THEORETICAL STUDY WITH APPLICATION

Graphene as a carbon monolayer has attracted extensive research interest in recent years. My research work within the frame of density functional theory has suggested that positioning graphene in proximity to h-BN may induce a finite energy gap in graphene, which is important for device applications. For an AB-stacked graphene/BN bilayer, a finite gap is induced at the equilibrium configuration. This induced gap shows a linear relationship with the applied strain. For a graphene/BN/graphene trilayer, a negligible gap is predicted in the ground state due to the overall symmetry of the system. When an electric field is applied, a tunable gap can be obtained for both AAA and ABA stackings. Enhanced tunneling current in the AA-stacked bilayer nanoribbons is predicted compared to either single-layer or AB-stacked bilayer nanoribbons. Interlayer separation between the nanoribbons is shown to have a profound impact on the conducting features. The effect of boron or nitrogen doping on the electronic transport properties of C60 fullerene is studied. The BC59 fullerene exhibits a considerably higher current than the pristine or nitrogen doped fullerenes beyond the applied bias of 1 V, suggesting it can be an effective semiconductor in p-type devices. The interaction between nucleic acid bases - adenine (A), guanine (G), cytosine (C), thymine (T) and uracil (U) - and a hydrogen-passivated silicon nanowire (SiNW) is investigated. The binding energy of the bases with the SiNW shows the order: G > A~C~T~U. This suggests that the interaction strength of a hydrogen passivated SiNW with the nucleic acid bases is nearly the same-G being an exception. The nature of the interaction is suggested to be electrostatic.