Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Softening the Lower Crust: Modes of Syn-Transport Transposition Around and Adjacent to a Deep Crustal Granulite Nappe, Parry Sound Domain, Grenville Province, Ontario, Canada

The Parry Sound domain is a granulite nappe-stack transported cratonward during reactivation of the ductile lower and middle crust in the late convergence of the Mesoproterozoic Grenville orogeny. Field observations suggest the following with respect to the ductile sheath: (1) Formation of a carapace of transposed amphibolite facies gneiss derived from and enveloping the western extremity of the Parry Sound domain and separating it from high-strain gneiss of adjacent allochthons. This ductile sheath formed dynamically around the moving granulite nappe through the development of systems of progressively linked shear zones. (2) Transposition initiated by hydration (amphibolization) of granulite facies gneiss by introduction of fluid along cracks accompanying pegmatite emplacement. Shear zones nucleated along pegmatite margins and subsequently linked and rotated. The source of the pegmatites was most likely subjacent migmatitic and pegmatite-rich units or units over which Parry Sound domain was transported. Comparison of gneisses of the ductile sheath with high-strain layered gneiss of adjacent allochthons show the mode of transposition of penetratively layered gneiss depended on whether or not the gneiss protoliths were amphibolite or granulite facies tectonites before initiation of transposition, resulting in, e.g., folding before shearing, no folding before shearing, respectively. Meter-scale truncation along high-strain gradients at the margins of both types of transposition-related shear zones observed within and marginal to Parry Sound domain mimic features at kilometer scales, implying that apparent truncation by transposition originating in a manner similar to the ductile sheath may be a common feature of deep crustal ductile reworking. Citation: Culshaw, N., C. Gerbi, and J. Marsh (2010), Softening the lower crust: Modes of syn-transport transposition around and adjacent to a deep crustal granulite nappe, Parry Sound domain, Grenville Province, Ontario, Canada, Tectonics, 29, TC5013, doi:10.1029/2009TC002537.

Fractures on comet 67P/Churyumov-Gerasimenko observed by Rosetta/OSIRIS

The Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard the Rosetta spacecraft currently orbiting comet 67P/Churyumov-Gerasimenko has yielded unprecedented views of a comet's nucleus. We present here the first ever observations of meter-scale fractures on the surface of a comet. Some of these fractures form polygonal networks. We present an initial assessment of their morphology, topology, and regional distribution. Fractures are ubiquitous on the surface of the comet's nucleus. Furthermore, they occur in various settings and show different topologies suggesting numerous formation mechanisms, which include thermal insulation weathering, orbital-induced stresses, and possibly seasonal thermal contraction. However, we conclude that thermal insolation weathering is responsible for creating most of the observed fractures based on their morphology and setting in addition to thermal models that indicate diurnal temperature ranges exceeding 200K and thermal gradients of similar to 15K/min at perihelion are possible. Finally, we suggest that fractures could be a facilitator in surface evolution and long-term erosion.

Lacustrine sediments from two setions in northern Greece

We investigated two lignite quarries in northern Greece for orbital and suborbital climate variability. Sections Lava and Vegora are located at the southern and northern boundaries of the Ptolemais Basin, a northwest southeast elongated intramontane basin that contains Upper Miocene to Lower Pliocene lacustrine sediments. Sediments show cyclic alterations of marl-rich (light), and coal-rich or clay-rich (dark) strata on a decimeter to meter scale. First, we established low-resolution ground-truth stratigraphy based on paleomagnetics and biostratigraphy. Accordingly, the lower 67 m and 65 m that were investigated in both sections Vegora and Lava, respectively, belong to the Upper Miocene and cover a time period of 6.85 to 6.57 and 6.46 to 5.98 Ma at sedimentation rates of roughly 14 and 22 cm/ka. In order to obtain a robust and high-resolution chronology, we then tuned carbonate minima (low L* values; high magnetic susceptibility values) to insolation minima. Besides the known dominance of orbital precession and eccentricity, we detected a robust hemi-precessional cycle in most parameters, most likely indicative for monsoonal influence on climate. Moreover, the insolation-forced time series indicate a number of millennial-scale frequencies that are statistically significant with dominant periods of 1.5-8 kyr. Evolutionary spectral analysis indicates that millennial-scale climate variability documented for the Ptolemais Basin resembles the one that is preserved in ice-core records of Greenland. Most cycles show durations of several tens of thousands of years before they diminish or cease. This is surprising because the generally argued cause for Late Quaternary millennial-scale variability is associated with the presence of large ice sheets, which cannot be the case for the Upper Miocene. Possible explanations maybe a direct response to solar forcing, an influence on the formation of North Atlantic Deep Water through the outflow of high-salinity water, or an atmospheric link to the North Atlantic Oscillation.

Geochemistry of the lower sheeted dike complex of ODP Holes 137-504B and 140-504B

Sixty-three samples representing 379 m of sheeted dikes from Deep Sea Drilling Project/Ocean Drilling Program Site 504B have been analyzed for major and selected trace elements by X-ray fluorescence. The samples range from microcrystalline aphyric basalts to moderately phyric (2%-10% phenocrysts) diabase that are typically multiply saturated with plagioclase, olivine, and clinopyroxene, in order of relative abundance. All analyzed samples are classified as Group D compositions with moderate to slightly elevated compatible elements (MgÆ-value = 0.65% ± 0.03%; Al2O3 = 15.5% ± 0.8%; CaO = 13.0% ± 0.3%; Ni = 114 ± 29 ppm), and unusually depleted levels of moderate to highly incompatible elements (Nb < 1 ppm; Zr = 44 ± 7 ppm; Rb < 0.5 ppm; Ba ~ 1 ppm; P2O5 = 0.07% ± 0.02%). These compositions are consistent with a multistage melting of a normal ocean ridge basaltic mantle source followed by extensive fractionation of olivine, plagioclase, and clinopyroxene. Leg 140 aphyric to sparsely phyric (0%-2% phenocrysts) basalts and diabases are compositionally indistinguishable from similarly phyric samples at higher levels in the hole. An examination of the entire crustal section, from the overlying volcanics through the sheeted dikes observed in Leg 140, reveals no significant trends indicating the enrichment or depletion of Costa Rica Rift Zone source magmas over time. Similarly, significant trends toward increased or decreased differentiation cannot be identified, although compositional patterns reflecting variable amounts of phenocryst addition are apparent at various depths. Below ? 1700 mbsf to the bottom of the Leg 140 section, there is a broadly systematic pattern of Zn depletion with depth, the result of high-temperature hydrothermal leaching. This zone of depletion is thought to be a significant source of Zn for the hydrothermal fluids depositing metal sulfides at ridge-crest hydrothermal vents and the sulfide-mineralization zone, located in the transition between pillow lavas and sheeted dikes. Localized zones of intense alteration (60%-95% recrystallization) are present on a centimeter to meter scale in many lithologic units. Within these zones, normally immobile elements Ti, Zr, Y, and rare-earth elements are strongly depleted compared with "fresher" samples centimeters away. The extent of compositional variability of these elements tends to obscure primary igneous trends if the highly altered samples are not identified or removed. At levels up to 40% (or possibly 60%) recrystallization, Ti, Zr, and Y retain their primary signatures. Although the mechanisms are unclear, it is possible that these intense alteration zones are a source of Y and rare-earth elements for the typically rare-earth-element-enriched hydrothermal vent fluids of mid-ocean ridges.

Geochemistry and diatom proxy record of ODP Site 127-797

Late Quaternary sediments of the Japan Sea are characterized by centimeter- to meter-scale alternations of dark and light layers which are synchronous basinwide. High-resolution analyses of the sediments from Ocean Drilling Program site 797 reveal that deposition of the meter-scale alternations reflect variations in paleoceanographic conditions which were closely associated with glacio-eustatic sea level changes through the modulation of the volume and character of the influx to the sea through the Tsushima Strait. The centimeter- to decimeter-scale alternations reflect millennial-scale variations which are possibly associated with Dansgaard-Oeschger (D-O) cycles, with each dark layer appearing to correspond to an interstadial. This variability is attributed to the development of a humid climate in central to eastern Asia and the consequent increase in discharge from the Huanghe and Changjiang Rivers during interstadials. This caused expansion of the East China Sea coastal water (ECSCW), which penetrated more strongly into the Japan Sea. The increased influence of the lower-salinity, nutrient-enriched ECSCW reduced deep water ventilation and enhanced the surface productivity, leading to the development of anoxic bottom waters and deposition of the dark layers. Thus the centimeter- to decimeter-scale alternations of the dark and light layers record wet and dry cycles in central to eastern Asia possibly associated with D-O cycles.

Sediment texture, carbonate content, and carbonate mineralogy of ODP Holes 182-1130A and 182-1132B

At Sites 1130 and 1132 of Ocean Drilling Program Leg 182 in the Great Australian Bight, we recovered an expanded Pleistocene section dominated by packstone and wackestone, deposited at unusually high rates of >20 cm/k.y. Shipboard observations detected an intermittent meter-scale alternation of light gray intervals with olive-gray intervals. Meter-scale samples were collected from the upper 250 m at both sites and decimeter-scale samples from four selected 2.5- to 4.0-m intervals in order to determine the texture and composition of sediments deposited along the upper slope throughout the Quaternary. Detailed textural and compositional data are presented from a total of 540 samples collected from both sites. Results indicate a general coarsening upward at both sites, with an accompanying upcore increase in high-Mg calcite (HMC) and aragonite and a decrease in low-Mg calcite (LMC). Samples collected at decimeter-scale intervals substantiate that the alternating light gray and olive-gray units detected on board ship are lithologically distinct. Light gray units consist of an LMC-rich silt, whereas olive-gray units consist of an aragonite and HMC-rich sand and silt. Sediment sources as well as timing and controls of this cyclic depositional pattern will be the subject of further investigations.

Grain size distribution of Oligocene and Miocene sediments from the Indus Fan

The Owen Ridge south of Oman represents oceanic crust that was uplifted by compressional tectonic forces in the early Miocene. Build-out of the Indus Fan led to deposition of a thick sequence of turbidites over the site of the Ridge during the late Oligocene and early Miocene. Early Miocene uplift of the Ridge led to a pelagic cap of nannofossil chalks. Two short sequences of turbidites from the pre- and syn-uplift phases were chosen for detailed grain size analysis. The upper Oligocene section at Site 731 is composed of thin (centimeter-decimeter scale) graded mud turbidites separated by relatively thick (decimeter-meter scale) intervals of homogeneous, non-bioturbated clayey siltstones. These finer intervals are unusually silt-rich (about 60%) for ungraded material and were probably deposited as undifferentiated muds from a series of turbidity current tails. By contrast, the lower Miocene section at Site 722 is comprised of a sequence of interbedded turbidites and hemipelagic carbonates. Sharp-based silt turbidites are overlain by burrow-mottled marly nannofossil chalks. The Oligocene sequence may have accumulated in an overbank setting on the middle fan - the local topographic position favoring frequent deposition from turbidity current tails and occasional deposition from the body of a turbidity flow. Uplift of the Ridge in the early Miocene led to pelagic carbonate deposition interrupted only by turbidity currents capable of overcoming a topographic barrier. Further uplift eventually led to entirely pelagic carbonate deposition.

Mineraly of oceanic crust from the Kane Transform, Mid-Atlantic Ridge

Gabbroic cumulates drilled south of the Kane Transform Fault on the slow-spread Mid-Atlantic Ridge preserve up to three discrete magnetization components. Here we use absolute age constraints derived from the paleomagnetic data to develop a model for the magmatic construction of this section of the lower oceanic crust. By comparing the paleomagnetic data with mineral compositions, and based on thermal models of local reheating, we infer that magmas that began crystallizing in the upper mantle intruded into the lower oceanic crust and formed meter-scale sills. Some of these magmas were crystal-laden and the subsequent expulsion of interstitial liquid from them produced '"cumulus" sills. These small-scale magmatic injections took place over at least 210 000 years and at distances of ~3 km from the ridge axis and may have formed much of the lower crust. This model explains many of the complexities described in this area and can be used to help understand the general formation of oceanic crust at slow-spread ridges.

Mineralogy of ODP Site 105-645

Cores from the upper 70 meters below seafloor (mbsf) (upper Pleistocene) at Ocean Drilling Program (ODP) Site 645 in Baffin Bay show dramatic meter-scale changes in color and mineralogy. Below this interval, mineralogical changes are more gradual to the top of the Miocene at about 550 mbsf. The Pliocene-Pleistocene section can be divided into five facies: Facies 1 - massive, poorly sorted, gravel-bearing muds; Facies 2 - gray silty clays and silty muds; Facies 3 - laminated detricarbonate silty muds; Facies 4 - silty sand and sandy silt; and Facies 5 - poorly sorted muddy sands and silty muds. Facies 4 and 5 are restricted to the Pliocene section below depths of about 275 mbsf. The mineralogical/color cycles in the upper 70 mbsf are the result of alternations between Facies 2 and three lithotypes of Facies 1: lithotype A - tan-colored, carbonate-rich, gravel-bearing mud; lithotype B - weak, red-colored, gravel-bearing mud rich in sedimentary rock fragments; and lithotype C - gray, gravel-bearing mud. A fourth lithotype, D, is restricted to depths of 168-275 mbsf and is dark gray, carbonate-poor, gravel-bearing mud. We believe that all lithotypes of Facies 1 and the sand and gravel fractions of Facies 2 and 3 were deposited by ice rafting. Depositional processes for Facies 4 and 5 probably include ice rafting and bottom- and turbidity-current transport. Data from petrographic analyses of light and heavy sand-sized grains and X-ray analyses of silt- and clay-size fractions suggest that tan-colored sediments (lithotype A of Facies 1; Facies 3) were derived mainly from Paleozoic carbonates of Ellesmere, Devon, and northern Baffin islands. Weak red sediments (lithotype B) contain significant red sedimentary clasts, reworked quartzarenite grains and clasts, and rounded colorless garnets, all derived from Proterozoic sequences of the Borden and Thule basins, and from minor Mesozoic red beds. Other sediments in the upper 335 mbsf at Site 645 contain detritus from a heterogeneous mixture of sources, including Precambrian shield terranes around Baffin Bay. Sediments from 335 to 550 mbsf (Facies 5) are rich in friable sedimentary clasts and detrital micas and contain glauconite and, in a few samples, reworked diatoms. These components suggest derivation from poorly consolidated Mesozoic-Tertiary sediments in coastal outcrops and beneath the modern shelves of northeastern Baffin Island and western Greenland. For the upper Pleistocene section (about 0-100 mbsf), marked mineralogical cyclicity is attributed to fluctuating glacial margins, calving rates, and iceberg melting rates, particularly around the northern end of Baffin Bay. Tan-colored, carbonate-rich units were derived at times of maximum advance of glaciers on Ellesmere and Devon islands, during relatively warm intervals induced by incursion of warm Atlantic surface water into the bay. At the beginning of these warmer episodes, most icebergs were contributed by glaciers near sea level around the Arctic channels, which resulted in deposition of weak red, ice-rafted units rich in Proterozoic sedimentary clasts.

Beryllium-10 and aluminium-26 concentration and production rates in glacial deposits from the Quartermain Mountains

Many glacial deposits in the Quartermain Mountains, Antarctica present two apparent contradictions regarding the degradation of unconsolidated deposits. The glacial deposits are up to millions of years old, yet they have maintained their meter-scale morphology despite the fact that bedrock and regolith erosion rates in the Quartermain Mountains have been measured at 0.1-4.0 m/Ma. Additionally, ground ice persists in some Miocene-aged soils in the Quartermain Mountains even though modeled and measured sublimation rates of ice in Antarctic soils suggest that without any recharge mechanisms ground ice should sublimate in the upper few meters of soil on the order of 10**3 to 10**5 years. This paper presents results from using the concentration of cosmogenic nuclides beryllium-10 (10Be) and aluminum-26 (26Al) in bulk sediment samples from depth profiles of three glacial deposits in the Quartermain Mountains. The measured nuclide concentrations are lower than expected for the known ages of the deposits, erosion alone does not always explain these concentrations, and deflation of the tills by the sublimation of ice coupled with erosion of the overlying till can explain some of the nuclide concentration profiles. The degradation rates that best match the data range 0.7-12 m/Ma for sublimation of ice with initial debris concentrations ranging 12-45% and erosion of the overlying till at rates of 0.4-1.2 m/Ma. Overturning of the tills by cryoturbation, vertical mixing, or soil creep is not indicated by the cosmogenic nuclide profiles, and degradation appears to be limited to within a few centimeters of the surface. Erosion of these tills without vertical mixing may partially explain how some glacial deposits in the Quartermain Mountains maintain their morphology and contain ground ice close to the surface for millions of years.

Chemistry of secondary minerals from the deep sheeted dike complex of ODP Holes 137-504B and 140-504B

Dolerites sampled from the lower sheeted dikes from Hole 504B during Ocean Drilling Program Legs 137 and 140, between 1562.4 and 2000.4 mbsf, were examined to document the mineralogy, petrography, and mineral parageneses associated with secondary alteration, to constrain the thermal history and composition of hydrothermal fluids. The main methods used were mineral chemical analyses by electron microprobe, X-ray diffraction, and cathodoluminescence microscopy. Temperatures of alteration were estimated on the basis of single and/or coexisting mineral chemistry. Permeability is important in controlling the type and extent of alteration in the studied dike section. At the meter-scale, intervals of weakly altered dolerites containing fresh olivine are interpreted as having experienced restricted exposure to hydrothermal fluids. At the centimeter- or millimeter-scale, alteration patches and extensively altered halos adjacent to veins reflect the permeability related to intergranular primary porosity and cracks. Most of the sheeted dike alteration in this case resulted from non-focused, pervasive fluid-rock interaction. This study confirms and extends the previous model for hydrothermal alteration at Hole 504B: hydrothermal alteration at the ridge axis followed by seawater recharge and off-axis alteration. The major new discoveries, all related to higher temperatures of alteration, are: (1) the presence of hydrothermal plagioclase (An80-95), (2) the presence of deuteric and/or hydrothermal diopside, and (3) the general increasing proportion of amphiboles, and particularly magnesio-hornblende with depth. We propose that the dolerites at Hole 504B were altered in five stages. Stage 1 occurred at high temperatures (less than 500° to 700°C) and involved late-magmatic formation of Na- and Ti-rich diopside, the hydrothermal formation of Na, Ti-poor diopside and the hydrothermal formation of an assemblage of An-rich plagioclase + hornblende. Stage 2 occurred at lower temperatures (250°-320°C) and is characterized by the appearance of actinolite, chlorite, chlorite-smectite, and/or talc (in low permeability zones) and albite. During Stage 3, quartz and epidote precipitated from evolved hydrothermal fluids at temperatures between 310° and 320°C. Anhydrite appeared during Stage 4 and likely precipitated directly from heated seawater. Stage 5 occurred off-axis at low temperatures (250°C) with laumontite and prehnite from evolved fluids.