Effect of high order interpolation in the stability and efficiency of the time-integration process in vorticity-velocity CFD algorithms

The numerical solution of the incompressible Navier-Stokes Equations offers an effective alternative to the experimental analysis of Fluid-Structure interaction i.e. dynamical coupling between a fluid and a solid which otherwise is very complex, time consuming and very expensive. To have a method which can accurately model these types of mechanical systems by numerical solutions becomes a great option, since these advantages are even more obvious when considering huge structures like bridges, high rise buildings, or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the KLE to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ODE time integration schemes and thus allows us to tackle the various multiphysics problems as separate modules. The current algorithm for KLE employs a structured or unstructured mesh for spatial discretization and it allows the use of a self-adaptive or fixed time step ODE solver while dealing with unsteady problems. This research deals with the analysis of the effects of the Courant-Friedrichs-Lewy (CFL) condition for KLE when applied to unsteady Stoke’s problem. The objective is to conduct a numerical analysis for stability and, hence, for convergence. Our results confirmthat the time step ∆t is constrained by the CFL-like condition ∆t ≤ const. hα, where h denotes the variable that represents spatial discretization.

TIME LAPSE SEISMIC OBSERVATIONS AND EFFECTS OF RESERVOIR COMPRESSIBILITY AT TEAL SOUTH OIL FIELD

One of the original ocean-bottom time-lapse seismic studies was performed at the Teal South oil field in the Gulf of Mexico during the late 1990’s. This work reexamines some aspects of previous work using modern analysis techniques to provide improved quantitative interpretations. Using three-dimensional volume visualization of legacy data and the two phases of post-production time-lapse data, I provide additional insight into the fluid migration pathways and the pressure communication between different reservoirs, separated by faults. This work supports a conclusion from previous studies that production from one reservoir caused regional pressure decline that in turn resulted in liberation of gas from multiple surrounding unproduced reservoirs. I also provide an explanation for unusual time-lapse changes in amplitude-versus-offset (AVO) data related to the compaction of the producing reservoir which, in turn, changed an isotropic medium to an anisotropic medium. In the first part of this work, I examine regional changes in seismic response due to the production of oil and gas from one reservoir. The previous studies primarily used two post-production ocean-bottom surveys (Phase I and Phase II), and not the legacy streamer data, due to the unavailability of legacy prestack data and very different acquisition parameters. In order to incorporate the legacy data in the present study, all three poststack data sets were cross-equalized and examined using instantaneous amplitude and energy volumes. This approach appears quite effective and helps to suppress changes unrelated to production while emphasizing those large-amplitude changes that are related to production in this noisy (by current standards) suite of data. I examine the multiple data sets first by using the instantaneous amplitude and energy attributes, and then also examine specific apparent time-lapse changes through direct comparisons of seismic traces. In so doing, I identify time-delays that, when corrected for, indicate water encroachment at the base of the producing reservoir. I also identify specific sites of leakage from various unproduced reservoirs, the result of regional pressure blowdown as explained in previous studies; those earlier studies, however, were unable to identify direct evidence of fluid movement. Of particular interest is the identification of one site where oil apparently leaked from one reservoir into a “new” reservoir that did not originally contain oil, but was ideally suited as a trap for fluids leaking from the neighboring spill-point. With continued pressure drop, oil in the new reservoir increased as more oil entered into the reservoir and expanded, liberating gas from solution. Because of the limited volume available for oil and gas in that temporary trap, oil and gas also escaped from it into the surrounding formation. I also note that some of the reservoirs demonstrate time-lapse changes only in the “gas cap” and not in the oil zone, even though gas must be coming out of solution everywhere in the reservoir. This is explained by interplay between pore-fluid modulus reduction by gas saturation decrease and dry-frame modulus increase by frame stiffening. In the second part of this work, I examine various rock-physics models in an attempt to quantitatively account for frame-stiffening that results from reduced pore-fluid pressure in the producing reservoir, searching for a model that would predict the unusual AVO features observed in the time-lapse prestack and stacked data at Teal South. While several rock-physics models are successful at predicting the time-lapse response for initial production, most fail to match the observations for continued production between Phase I and Phase II. Because the reservoir was initially overpressured and unconsolidated, reservoir compaction was likely significant, and is probably accomplished largely by uniaxial strain in the vertical direction; this implies that an anisotropic model may be required. Using Walton’s model for anisotropic unconsolidated sand, I successfully model the time-lapse changes for all phases of production. This observation may be of interest for application to other unconsolidated overpressured reservoirs under production.

Omega-3 fatty acids and the neurovascular responses to mental stress in humans

Hypertension is the most prevalent form of cardiovascular disease (CVD) in the world, and is known to increase the risk for developing other diseases. Recently, the American Heart Association introduced a new classification of blood pressure, prehypertension (PHT). The criteria for PHT include a systolic of 120-139 mmHg and/or a diastolic blood pressure of 80-89 mmHg. It has been observed that individuals with PHT have a higher risk of developing hypertension later in life. Therefore, it is important to understand the mechanisms contributing to PHT in order to possibly prevent hypertension. Omega-3 fatty acids found in fish oils have been suggested as a means of lowering blood pressure. However, little is known on the effects of fish oil in PHT humans. Therefore we conducted two studies. In Study 1 we investigated PHT and normotensive (NT) individuals during a mental stress task. Mental stress is known to contribute to the development of hypertension. In Study 2 PHT and NT subjects were placed in an eight week double-blind placebo controlled study in which subjects consumed 9g/day of either fish oil or placebo (olive oil) in addition to their regular diets. Subjects were tested during a resting baseline (seated and supine), 5 minutes of a mental stress task, and 5 minutes of recovery both pre and post supplementation. We measured arterial pressure (AP), heart rate (HR), muscle sympathetic nerve activity (MSNA), and forearm and calf vascular responses. In Study 1 PHT demonstrated augmented AP and blunted vasodilation during mental stress, but MSNA did not change. In Study 2, fish oil did not directly influence blood pressure, MSNA or vascular responses to mental stress. However, it became clear that fish oil had an effect on some but not all subjects (both PHT and NT). Specifically, subjects who experienced a reduced blood pressure response to fish oil also demonstrated a decrease in MSNA and HR during mental stress. Collectively, the investigations in this dissertation had several novel findings. First, PHT individuals demonstrate an augmented pressor and blunted vascular response to mental stress, a response that may be contributing to the development of hypertension. Second, fish oil does not consistently lower resting blood pressure, but the interindividual responses may be related to MSNA. Third, fish oil attenuated the heart rate and MSNA responses and to mental stress in both PHT and NT. In conclusion, we found that there are both similarities and differences in the way PHT and NT individuals respond to mental stress and fish oil.

Planar magneto-photonic and gradient-photonic structures : crystals and metamaterials.

In the field of photonics, two new types of material structures, photonic crystals and metamaterials, are presently of great interest. Both are studied in the present work, which focus on planar magnetic materials in the former and planar gradient metamaterials in the latter. These planar periodic structures are easy to handle and integrate into optical systems. The applications are promising field for future optical telecommunication systems and give rise to new optical, microwave and radio technologies. The photonic crystal part emphasizes the utilization of magnetic material based photonic crystals due to its remarkable magneto-optical characteristics. Bandgaps tuning by magnetic field in bismuth-gadolinium-substituted lutetium iron garnet (Bi0.8 Gd0.2 Lu2.0 Fe5 O12) based one- dimensional photonic crystals are investigated and demonstrated in this work. Magnetic optical switches are fabricated and tested. Waveguide formulation for band structure in magneto photonic crystals is developed. We also for the first time demonstrate and test two- dimensional magneto photonic crystals optical. We observe multi-stopbands in two- dimensional photonic waveguide system and study the origin of multi-stopbands. The second part focus on studying photonic metamaterials and planar gradient photonic metamaterial design. We systematically study the effects of varying the geometry of the fishnet unit cell on the refractive index in optical frequency. It is the first time to design and demonstrate the planar gradient structure in the high optical frequency. Optical beam bending using planar gradient photonic metamaterials is observed. The technologies needed for the fabrication of the planar gradient photonic metamaterials are investigated. Beam steering devices, shifter, gradient optical lenses and etc. can be derived from this design.

Analysis of atmospheric transport and the effects of seasonal and interannual transport variability on measurements made at the Pico Mountain observatory

A non-hierarchical K-means algorithm is used to cluster 47 years (1960–2006) of 10-day HYSPLIT backward trajectories to the Pico Mountain (PM) observatory on a seasonal basis. The resulting cluster centers identify the major transport pathways and collectively comprise a long-term climatology of transport to the observatory. The transport climatology improves our ability to interpret the observations made there and our understanding of pollution source regions to the station and the central North Atlantic region. I determine which pathways dominate transport to the observatory and examine the impacts of these transport patterns on the O3, NOy, NOx, and CO measurements made there during 2001–2006. Transport from the U.S., Canada, and the Atlantic most frequently reaches the station, but Europe, east Africa, and the Pacific can also contribute significantly depending on the season. Transport from Canada was correlated with the North Atlantic Oscillation (NAO) in spring and winter, and transport from the Pacific was uncorrelated with the NAO. The highest CO and O3 are observed during spring. Summer is also characterized by high CO and O3 and the highest NOy and NOx of any season. Previous studies at the station attributed the summer time high CO and O3 to transport of boreal wildfire emissions (for 2002–2004), and boreal fires continued to affect the station during 2005 and 2006. The particle dispersion model FLEXPART was used to calculate anthropogenic and biomass-burning CO tracer values at the station in an attempt to identify the regions responsible for the high CO and O3 observations during spring and biomass-burning impacts in summer.

Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

Non-invasive method to predict viscosity and size using total internal reflection fluorescence microscopy (TIRFM) system

This study develops an automated analysis tool by combining total internal reflection fluorescence microscopy (TIRFM), an evanescent wave microscopic imaging technique to capture time-sequential images and the corresponding image processing Matlab code to identify movements of single individual particles. The developed code will enable us to examine two dimensional hindered tangential Brownian motion of nanoparticles with a sub-pixel resolution (nanoscale). The measured mean square displacements of nanoparticles are compared with theoretical predictions to estimate particle diameters and fluid viscosity using a nonlinear regression technique. These estimated values will be confirmed by the diameters and viscosities given by manufacturers to validate this analysis tool. Nano-particles used in these experiments are yellow-green polystyrene fluorescent nanospheres (200 nm, 500 nm and 1000 nm in diameter (nominal); 505 nm excitation and 515 nm emission wavelengths). Solutions used in this experiment are de-ionized (DI) water, 10% d-glucose and 10% glycerol. Mean square displacements obtained near the surface shows significant deviation from theoretical predictions which are attributed to DLVO forces in the region but it conforms to theoretical predictions after ~125 nm onwards. The proposed automation analysis tool will be powerfully employed in the bio-application fields needed for examination of single protein (DNA and/or vesicle) tracking, drug delivery, and cyto-toxicity unlike the traditional measurement techniques that require fixing the cells. Furthermore, this tool can be also usefully applied for the microfluidic areas of non-invasive thermometry, particle tracking velocimetry (PTV), and non-invasive viscometry.

Cloning and characterization of the genes involved in cambial growth in response to elevated [CO2]

Aspen (Populus tremuloides) trees growing under elevated [CO2] at a free-air CO2 enrichment (FACE) site have produced significantly more biomass compared to control trees. The molecular mechanisms underlying the observed increase in biomass productivity was investigated by producing transcriptomic profiles of the vascular cambium zone (VCZ) and leaves, followed by a comparative study to identify genes and pathways that showed significant changes following long-term exposure to elevated [CO2]. This study is mainly to verify if genetic modification of a few selected candidate genes including CAP1, CKX6, and ASML2 that are expressed in vascular cambium in response to elevated [CO2] can cause the changes in plant growth and development. To this end, these three genes were cloned into both sense and antisense constructs. Then antisense and sense transgenic lines of above-mentioned genes were developed. 15 events were generated for 5 constructs, which were confirmed with regular PCR and RT-PCR. Confirmed plants were planted in greenhouse for growth and phenotypic characterization. The expression of CAP1, CKX6 and ASML2 in antisense plants was measured by real-time RT-PCR, and the changes caused by gene interference in cambial growth were studies by analyzing the microscopic sections made from the antisense transgenic plants. It has been found that 1) CAP1 is mainly expressed in xylem and root. 2) RNAi suppression of CAP1 significantly affected height and diameter. 3) CAP1, ASML2 and CKX6 affected xylem and phloem cell proliferation and elongation. Due to the delay in regenerating sense transgenic plants, the characterization of sense transgenic plants is limited to growth only.

Understanding volcanic processes using UV camera measurements of sulfur dioxide and coincident infrasound and seismicity

The exsolution of volatiles from magma maintains an important control on volcanic eruption styles. The nucleation, growth, and connectivity of bubbles during magma ascent provide the driving force behind eruptions, and the rate, volume, and ease of gas exsolution can affect eruptive activity. Volcanic plumes are the observable consequence of this magmatic degassing, and remote sensing techniques allow us to quantify changes in gas exsolution. However, until recently the methods used to measure volcanic plumes did not have the capability of detecting rapid changes in degassing on the scale of standard geophysical observations. The advent of the UV camera now makes high sample rate gas measurements possible. This type of dataset can then be compared to other volcanic observations to provide an in depth picture of degassing mechanisms in the shallow conduit. The goals of this research are to develop a robust methodology for UV camera field measurements of volcanic plumes, and utilize this data in conjunction with seismoacoustic records to illuminate degassing processes. Field and laboratory experiments were conducted to determine the effects of imaging conditions, vignetting, exposure time, calibration technique, and filter usage on the UV camera sulfur dioxide measurements. Using the best practices determined from these studies, a field campaign was undertaken at Volcán de Pacaya, Guatemala. Coincident plume sulfur dioxide measurements, acoustic recordings, and seismic observations were collected and analyzed jointly. The results provide insight into the small explosive features, variations in degassing rate, and plumbing system of this complex volcanic system. This research provides useful information for determining volcanic hazard at Pacaya, and demonstrates the potential of the UV camera in multiparameter studies.