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.

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.