Residual stress effects during near-threshold fatigue crack growth

The fatigue crack propagation behaviour of a low alloy, boron-containing steel has been examined after austenitizing at 900°C or 1250°C and tempering at a range of temperatures up to 400°C. Fatigue threshold values were found to vary with austenitizing and tempering treatment in a range between 3.3 to 6 MPa √m when tested at a stress ratio (R) of 0.2. Crack propagation rates in the Paris regime were insensitive to heat treatment variations. The crack propagation path was essentially transgranular in all conditions with small regions of intergranular facets appearing at growth rates around the knee of the da/dN vs ΔK curve. The crack front shape showed marked retardation in the centre of the specimen at low tempering temperatures. Experimental determinations and computer predictions of residual stress levels in the specimens indicated that this was due to a central residual compressive stress resulting from differential cooling rates and the volume change associated with the martensite transformation. The results are discussed in terms of microstructural and residual stress effects on fatigue behaviour. © 1987.

Lichen growth and lichenometry

Lichenometry is one of the most widely used methods available for dating the surface age of various substrata including rock surfaces, boulders, walls, and archaeological remains. It depends on the assumption that if the lag time before colonisation of a substratum by a lichen is known and lichen age can be estimated, then a minimum date can be obtained by measuring the diameter (or another property related to size) of the largest lichen at the site. Lichen age can be determined by variety of methods including calibrating lichen size against surfaces of known age (‘indirect lichenometry’), by constructing a growth rate-size curve from direct measurement of lichen growth (‘direct lichenometry’), using radio-carbon (RC) dating, and from lichen ‘growth rings’. This chapter describes: (1) lichen growth rates and longevity, (2) methods of estimating lichen age, (3) the methodology of lichenometry and (4) applications of lichenometry. Despite its limitations, lichenometry is likely to continue to play an important role in dating a variety of surfaces and also in providing data that contribute to the debate regarding global warming and climate change.

Derivation of Context-free Stochastic L-Grammar Rules for Promoter Sequence Modeling Using Support Vector Machine

Formal grammars can used for describing complex repeatable structures such as DNA sequences. In this paper, we describe the structural composition of DNA sequences using a context-free stochastic L-grammar. L-grammars are a special class of parallel grammars that can model the growth of living organisms, e.g. plant development, and model the morphology of a variety of organisms. We believe that parallel grammars also can be used for modeling genetic mechanisms and sequences such as promoters. Promoters are short regulatory DNA sequences located upstream of a gene. Detection of promoters in DNA sequences is important for successful gene prediction. Promoters can be recognized by certain patterns that are conserved within a species, but there are many exceptions which makes the promoter recognition a complex problem. We replace the problem of promoter recognition by induction of context-free stochastic L-grammar rules, which are later used for the structural analysis of promoter sequences. L-grammar rules are derived automatically from the drosophila and vertebrate promoter datasets using a genetic programming technique and their fitness is evaluated using a Support Vector Machine (SVM) classifier. The artificial promoter sequences generated using the derived L- grammar rules are analyzed and compared with natural promoter sequences.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.

Mechanistic modeling of fracture in asphalt mixtures under compressive loading

When an asphalt mixture is subjected to a destructive compressive load, it experiences a sequence of three deformation stages, as follows: the (1) primary, (2) secondary, and (3) tertiary stages. Most literature research focuses on plastic deformation in the primary and secondary stages, such as prediction of the flow number, which is in fact the initiation of the tertiary stage. However, little research effort has been reported on the mechanistic modeling of the damage that occurs in the tertiary stage. The main objective of this paper is to provide a mechanistic characterizing method for the damage modeling of asphalt mixtures in the tertiary stage. The preliminary study conducted by the writers illustrates that deformation during the tertiary flow of the asphalt mixtures is principally caused by the formation and propagation of cracks, which was signaled by the increase of the phase angle in the tertiary phase. The strain caused by the growth of cracks is the viscofracture strain, which can be obtained by conducting the strain decomposition of the measured total strain in the destructive compressive test. The viscofracture strain is employed in the research reported in this paper to mechanistically characterize the time-dependent fracture (viscofracture) of asphalt mixtures in compression. By using the dissipated pseudostrain energy-balance principle, the damage density and true stress are determined and both are demonstrated to increase with load cycles in the tertiary stage. The increased true stress yields extra viscoplastic strain, which is the reason why the permanent deformation is accelerated by the occurrence of cracks. To characterize the evolution of the viscofracture in the asphalt mixtures in compression, a pseudo J-integral Paris' law in terms of damage density is proposed and the material constants in the Paris' law are determined, which can be employed to predict the fracture of asphalt mixtures in compression. © 2013 American Society of Civil Engineers.

On the Mathematical Modelling of Microbial Growth: Some Computational Aspects

We propose a new approach to the mathematical modelling of microbial growth. Our approach differs from familiar Monod type models by considering two phases in the physiological states of the microorganisms and makes use of basic relations from enzyme kinetics. Such an approach may be useful in the modelling and control of biotechnological processes, where microorganisms are used for various biodegradation purposes and are often put under extreme inhibitory conditions. Some computational experiments are performed in support of our modelling approach.

The effects of climate on the growth and physiology of tropical rainforest trees

Tropical rainforests account for more than a third of global net primary production and contain more than half of the global forest carbon. Though these forests are a disproportionately important component of the global carbon cycle, the relationship between rainforest productivity and climate remains poorly understood. Understanding the link between current climate and rainforest tree stem diameter increment, a major constituent of forest productivity, will be crucial to efforts at modeling future climate and rainforest response to climate change. This work reports the physiological and stem growth responses to micrometeorological and phenological states of ten species of canopy trees in a Costa Rican wet tropical forest at sub-annual time intervals. I measured tree growth using band dendrometers and estimated leaf and reproductive phenological states monthly. Electronic data loggers recorded xylem sap flow (an indicator of photosynthetic rate) and weather at half-hour intervals. An analysis of xylem sap flow showed that physiological responses were independent of species, which allowed me to construct a general model of weather driven sap flow rates. This model predicted more than eighty percent of climate driven sap flow variation. Leaf phenology influenced growth in three of the ten species, with two of these species showing a link between leaf phenology and weather. A combination of rainfall, air temperature, and irradiance likely provided the cues that triggered leaf drop in Dipteryx panamensis and Lecythis ampla. Combining the results of the sap flow model, growth, and the climate measures showed tree growth was correlated to climate, though the majority of growth variation remained unexplained. Low variance in the environmental variables and growth rates likely contributed to the large amount of unexplained variation. A simple model that included previous growth increment and three meteorological variables explained from four to nearly fifty percent of the growth variation. Significant growth carryover existed in six of the ten species, and rainfall was positively correlated to growth in eight of the ten species. Minimum nighttime temperature was also correlated to higher growth rates in five of the species and irradiance in two species. These results indicate that tropical rainforest tree trunks could act as carbon sinks if future climate becomes wetter and slightly warmer. ^

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

The inclusion of differentially assessed lands in urban growth model calibration: a comparison of two approaches using SLEUTH

This research sought to understand the role that differentially assessed lands (lands in the United States given tax breaks in return for their guarantee to remain in agriculture) play in influencing urban growth. Our method was to calibrate the SLEUTH urban growth model under two different conditions. The first used an excluded layer that ignored such lands, effectively rendering them available for development. The second treated those lands as totally excluded from development. Our hypothesis was that excluding those lands would yield better metrics of fit with past data. Our results validate our hypothesis since two different metrics that evaluate goodness of fit both yielded higher values when differentially assessed lands are treated as excluded. This suggests that, at least in our study area, differential assessment, which protects farm and ranch lands for tenuous periods of time, has indeed allowed farmland to resist urban development. Including differentially assessed lands also yielded very different calibrated coefficients of growth as the model tried to account for the same growth patterns over two very different excluded areas. Excluded layer design can greatly affect model behavior. Since differentially assessed lands are quite common through the United States and are often ignored in urban growth modeling, the findings of this research can assist other urban growth modelers in designing excluded layers that result in more accurate model calibration and thus forecasting.

Modeling Nymphoides architecture: A morphological analysis of Nymphoides aquatica (Menyanthaceae)

• Premise of the study: Species in the aquatic genus Nymphoides have inflorescences that appear to arise from the petioles of floating leaves. The inflorescence-floating leaf complex can produce vegetative propagules and/or additional inflorescences and leaves. We analyzed the morphology of N. aquatica to determine how this complex relates to whole plant architecture and whether whole plant growth is sympodial or monopodial. • Methods: We used dissections, measurements, and microscopic observations of field-collected plants and plants cultivated for 2 years in outdoor tanks in south Florida, USA. • Key results: Nymphoides aquatica had a submerged plagiotropic rhizome that produced floating leaves in an alternate/spiral phyllotaxy. Rhizomes were composed of successive sympodial units that varied in the number of leaves produced before the apex terminated. The basic sympodial unit had a prophyll that subtended a renewal-shoot bud, a short-petioled leaf (SPL) with floating lamina, and an inflorescence; the SPL axillary bud expanded as a vegetative propagule. Plants produced either successive basic sympodial units or expanded sympodia that intercalated long-petioled leaves between the prophyll and the SPL. • Conclusions: Nymphoides aquatica grows sympodially, forming a rhizome composed of successive basic sympodia and expanded sympodial units. Variations on these types of sympodial growth help explain the branching patterns and leaf morphologies described for other Nymphoides species. Monitoring how these two sympodial phases are affected by water depth provides an ecologically meaningful way to assess N. aquatica’s responses to altered hydrology.

Essays in economic growth, political economy and institutions

This dissertation analyzes the effects of political and economic institutions on economic development and growth.^ The first essay develops an overlapping-generations political economy model to analyze the incentives of various social groups to finance human capital accumulation through public education expenditures. The contribution of this study to the literature is that it helps explain the observed differences in the economic growth performance of natural resource-abundant countries. The results suggest that the preferred tax rates of the manufacturers on one hand and the political coalition of manufacturers and landowners, on the other hand, are equal to the socially optimal tax rate. However, we show that owners of natural resources prefer an excessively high tax rate, which suppresses aggregate output to a suboptimal level.^ The second essay examines the relationship between the political influence of different social classes and public education spending in panel data estimation. The novel contribution of this paper to the literature is that I proxy the political power and influence of the natural resource owners, manufacturers, and landowners with macroeconomic indicators. The motivation behind this modeling choice is to substantiate the definition of the political power of social classes with economic fundamentals. I use different governance indicators in the estimations to find out how different institutions mediate the overall impact of the political influence of various social classes on public education spending. The results suggest that political stability and absence of violence and rule of law are the important governance indicators.^ The third essay develops a counter argument to Acemoglu et al. (2010) where the thesis is that French institutions and economic reforms fostered economic progress in those German regions invaded by the Napoleonic armies. By providing historical data on urbanization rates used as proxies for economic growth, I demonstrate that similar different rates of economic growth were observed in the regions of France in the post-Napoleonic period as well. The existence of different economic growth rates makes it hard to argue that the differences in economic performance in the German regions that were invaded by the French and those that were spared a similar fate follow from regional differences in economic institutions.^

Site specific growth of metal catalyzed silica nanowires for biological and chemical sensing

In this research the integration of nanostructures and micro-scale devices was investigated using silica nanowires to develop a simple yet robust nanomanufacturing technique for improving the detection parameters of chemical and biological sensors. This has been achieved with the use of a dielectric barrier layer, to restrict nanowire growth to site-specific locations which has removed the need for post growth processing, by making it possible to place nanostructures on pre-pattern substrates. Nanowires were synthesized using the Vapor-Liquid-Solid growth method. Process parameters (temperature and time) and manufacturing aspects (structural integrity and biocompatibility) were investigated. Silica nanowires were observed experimentally to determine how their physical and chemical properties could be tuned for integration into existing sensing structures. Growth kinetic experiments performed using gold and palladium catalysts at 1050°C for 60 minutes in an open-tube furnace yielded dense and consistent silica nanowire growth. This consistent growth led to the development of growth model fitting, through use of the Maximum Likelihood Estimation (MLE) and Bayesian hierarchical modeling. Transmission electron microscopy studies revealed the nanowires to be amorphous and X-ray diffraction confirmed the composition to be SiO2 . Silica nanowires were monitored in epithelial breast cancer media using Impedance spectroscopy, to test biocompatibility, due to potential in vivo use as a diagnostic aid. It was found that palladium catalyzed silica nanowires were toxic to breast cancer cells, however, nanowires were inert at 1μg/mL concentrations. Additionally a method for direct nanowire integration was developed that allowed for silica nanowires to be grown directly into interdigitated sensing structures. This technique eliminates the need for physical nanowire transfer thus preserving nanowire structure and performance integrity and further reduces fabrication cost. Successful nanowire integration was physically verified using Scanning electron microscopy and confirmed electrically using Electrochemical Impedance Spectroscopy of immobilized Prostate Specific Antigens (PSA). The experiments performed above serve as a guideline to addressing the metallurgic challenges in nanoscale integration of materials with varying composition and to understanding the effects of nanomaterials on biological structures that come in contact with the human body.

Simulating everglades national park hydrology and phosphorus transport under existing and future scenarios using numerical modeling

The Florida Everglades has a long history of anthropogenic changes which have impacted the quantity and quality of water entering the system. Since the construction of Tamiami Trail in the 1920's, overland flow to the Florida Everglades has decreased significantly, impacting ecosystems from the wetlands to the estuary. The MIKE Marsh Model of Everglades National Park (M3ENP) is a numerical model, which simulates Everglades National Park (ENP) hydrology using MIKE SHE/MIKE 11software. This model has been developed to determine the parameters that effect Everglades hydrology and understand the impact of specific flow changes on the hydrology of the system. ^ As part of the effort to return flows to the historical levels, several changes to the existing water management infrastructure have been implemented or are in the design phase. Bridge construction scenarios were programed into the M3ENP model to review the effect of these structural changes and evaluate the potential impacts on water levels and hydroperiods in the receiving Northeast Shark Slough ecosystem. These scenarios have shown critical water level increases in an area which has been in decline due to low water levels. Results from this work may help guide future decisions for restoration designs. ^ Excess phosphorus entering Everglades National Park in South Florida may promote the growth of more phosphorus-opportunistic species and alter the food chain from the bottom up. Two phosphorus transport methods were developed into the M3ENP hydrodynamic model to determine the factors affecting phosphorus transport and the impact of bridge construction on water quality. Results showed that while phosphorus concentrations in surface waters decreased overall, some areas within ENP interior may experience an increase in phosphorus loading which the addition of bridges to Tamiami Trail. Finally, phosphorus data and modeled water level data was used to evaluate the spectral response of Everglades vegetation to increasing phosphorus availability using Landsat imagery.^

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.