Implementation of monotonic higher order upwind scheme in KIVA 4

KIVA is a FORTRAN code developed by Los Alamos national lab to simulate complete engine cycle. KIVA is a flow solver code which is used to perform calculation of properties in a fluid flow field. It involves using various numerical schemes and methods to solve the Navier-Stokes equation. This project involves improving the accuracy of one such scheme by upgrading it to a higher order scheme. The numerical scheme to be modified is used in the critical final stage calculation called as rezoning phase. The primitive objective of this project is to implement a higher order numerical scheme, to validate and verify that the new scheme is better than the existing scheme. The latest version of the KIVA family (KIVA 4) is used for implementing the higher order scheme to support handling the unstructured mesh. The code is validated using the traditional shock tube problem and the results are verified to be more accurate than the existing schemes in reference with the analytical result. The convection test is performed to compare the computational accuracy on convective transfer; it is found that the new scheme has less numerical diffusion compared to the existing schemes. A four valve pentroof engine, an example case of KIVA package is used as application to ensure the stability of the scheme in practical application. The results are compared for the temperature profile. In spite of all the positive results, the numerical scheme implemented has a downside of consuming more CPU time for the computational analysis. The detailed comparison is provided. However, in an overview, the implementation of the higher order scheme in the latest code KIVA 4 is verified to be successful and it gives better results than the existing scheme which satisfies the objective of this project.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

EXPERIMENTAL INVESTIGATION OF CERTAIN INTERNAL CONDENSING AND BOILING FLOWS: THEIR SENSITIVITY TO PRESSURE FLUCTUATIONS AND HEAT TRANSFER ENHANCEMENTS

Space-based (satellite, scientific probe, space station, etc.) and millimeter – to – microscale (such as are used in high power electronics cooling, weapons cooling in aircraft, etc.) condensers and boilers are shear/pressure driven. They are of increasing interest to system engineers for thermal management because flow boilers and flow condensers offer both high fluid flow-rate-specific heat transfer capacity and very low thermal resistance between the fluid and the heat exchange surface, so large amounts of heat may be removed using reasonably-sized devices without the need for excessive temperature differences. However, flow stability issues and degradation of performance of shear/pressure driven condensers and boilers due to non-desirable flow morphology over large portions of their lengths have mostly prevented their use in these applications. This research is part of an ongoing investigation seeking to close the gap between science and engineering by analyzing two key innovations which could help address these problems. First, it is recommended that the condenser and boiler be operated in an innovative flow configuration which provides a non-participating core vapor stream to stabilize the annular flow regime throughout the device length, accomplished in an energy-efficient manner by means of ducted vapor re-circulation. This is demonstrated experimentally. Second, suitable pulsations applied to the vapor entering the condenser or boiler (from the re-circulating vapor stream) greatly reduce the thermal resistance of the already effective annular flow regime. For experiments reported here, application of pulsations increased time-averaged heat-flux up to 900 % at a location within the flow condenser and up to 200 % at a location within the flow boiler, measured at the heat-exchange surface. Traditional fully condensing flows, reported here for comparison purposes, show similar heat-flux enhancements due to imposed pulsations over a range of frequencies. Shear/pressure driven condensing and boiling flow experiments are carried out in horizontal mm-scale channels with heat exchange through the bottom surface. The sides and top of the flow channel are insulated. The fluid is FC-72 from 3M Corporation.

Mixing-cell boundary conditions and apparent mass balance errors for advective-dispersive solute transport

In many field or laboratory situations, well-mixed reservoirs like, for instance, injection or detection wells and gas distribution or sampling chambers define boundaries of transport domains. Exchange of solutes or gases across such boundaries can occur through advective or diffusive processes. First we analyzed situations, where the inlet region consists of a well-mixed reservoir, in a systematic way by interpreting them in terms of injection type. Second, we discussed the mass balance errors that seem to appear in case of resident injections. Mixing cells (MC) can be coupled mathematically in different ways to a domain where advective-dispersive transport occurs: by assuming a continuous solute flux at the interface (flux injection, MC-FI), or by assuming a continuous resident concentration (resident injection). In the latter case, the flux leaving the mixing cell can be defined in two ways: either as the value when the interface is approached from the mixing-cell side (MC-RT -), or as the value when it is approached from the column side (MC-RT +). Solutions of these injection types with constant or-in one case-distance-dependent transport parameters were compared to each other as well as to a solution of a two-layer system, where the first layer was characterized by a large dispersion coefficient. These solutions differ mainly at small Peclet numbers. For most real situations, the model for resident injection MC-RI + is considered to be relevant. This type of injection was modeled with a constant or with an exponentially varying dispersion coefficient within the porous medium. A constant dispersion coefficient will be appropriate for gases because of the Eulerian nature of the usually dominating gaseous diffusion coefficient, whereas the asymptotically growing dispersion coefficient will be more appropriate for solutes due to the Lagrangian nature of mechanical dispersion, which evolves only with the fluid flow. Assuming a continuous resident concentration at the interface between a mixing cell and a column, as in case of the MC-RI + model, entails a flux discontinuity. This flux discontinuity arises inherently from the definition of a mixing cell: the mixing process is included in the balance equation, but does not appear in the description of the flux through the mixing cell. There, only convection appears because of the homogeneous concentration within the mixing cell. Thus, the solute flux through a mixing cell in close contact with a transport domain is generally underestimated. This leads to (apparent) mass balance errors, which are often reported for similar situations and erroneously used to judge the validity of such models. Finally, the mixing cell model MC-RI + defines a universal basis regarding the type of solute injection at a boundary. Depending on the mixing cell parameters, it represents, in its limits, flux as well as resident injections. (C) 1998 Elsevier Science B.V. All rights reserved.

Intracranial meningiomas associated with cervical syringohydromyelia in a cat

Eine dreizehnjährige, weibliche, nicht kastrierte Hauskatze wurde zur chirurgischen Therapie eines multiplen Meningeoms in der Kleintierklinik des Tierspital Bern vorgestellt. Eine leichtgradige generalisierte Ataxie wurde beobachtet, die Propriozeption war generalisiert herabgesetzt und der Drohreflex war beidseitig reduziert. Es wurden eine prä- und eine postoperative Magnetrezonanztomographie durchgeführt. Drei supratentoriale extra-axiale Raumforderungen wurden diagnostiziert. Eine vierte Masse wurde infratentorial extra-axial über der linken zerebellären Hemisphere lokalisiert. Eine Herniation des kaudoventralen Kleinhirn (Vermis) von ungefähr einem Zentimeter Länge durch das Foramen magnum wurde beobachtet. Eine zervikale Syringohydromyelie wurde als Zufallsbefund diagnostiziert. Die Meningeome wurden durch 3 Kraniotomiestellen entfernt. Zwei Jahre nach der Operation ist die Katze normal. Anhand der vorhandenen Literatur wird die tumor-assoziierte Syringohydromyelie besprochen. Die Therapie der Syringohydromyelie sollte gegen den kausalen pathologischen Prozess (z.B. intrakranieller Tumor) für die Liquorzirkulationsstörung gerichtet sein.

Diversity in T cell costimulation by alpha4beta1 integrin

T cell activation requires antigen-specific T cell receptor signals that spatially and temporally coincide with a second costimulatory signal. CD28 and α4β1 integrin both function as T cell costimulators, but their individual mechanisms remain elusive. By directly comparing CD3-dependent functions and signaling pathways employed by these two costimulatory receptors, aspects of their individual signaling mechanisms are explored. We determined that CD28 and α4β1 integrins both use Src-family kinase Lck and MAPK Erk, but to different extents and functional ends. After identifying functional differences between CD28 and integrin costimulatory pathways, the focus of the study turned to integrin signaling in naïve and memory T cell subsets. CD45RO T cells are fully co-activated by natural β1 integrin ligands fibronectin (FN) and VCAM-1, β1 monoclonal antibody 33B6, as well as α4β1 monoclonal antibody 19H8 which binds a combinatorial epitope of the α4β1 heterodimer. While CD28 fully costimulates CD45RA T cells, the degree of activation from integrin ligands varies. FN costimulates CD3-dependent proliferation, IL-2 secretion, and early activation markers CD25 and CD69. However, β1 antibody 33B6, which binds to the same T cell integrins (α4β1 and α5β1) as natural ligand FN, failed to costimulate proliferation or IL-2 in the CD45RA subset, but retained the ability to regulate CD25 and CD69. Unique aspects of 19H8 signaling involve early Erk activation and IL-2 independent proliferation. Signaling defects through 33B6 ligation correlates with poor adhesion under fluid flow conditions, suggesting a cytoskeletal basis for signaling. All together, these data provide evidence for a mechanism of α4β1 integrin signaling and describe functional differences between naïve and memory T cells. ^

Geophysical logging on sediment core and hole CRP-3 from the Ross Sea, Antarctica

Cape Roberts drillhole CRP-3 in the northern part of McMurdo Sound (Ross Sea, Antarctica) targeted the western margin of the Victoria Land basin to investigate Neogene to Palaeogene climatic and tectonic history by obtaining continuous core and downhole logs (Cape Roberts Science Team, 2000). The CRP-3 drillhole extended to 939.42 mbsf (meters below seafloor) at a water depth of 297 m. The first downhole measurements after drilling were the temperature and salinity logs. Both were measured at the beginning and at the end of each of the three logging phases. Although an equilibrium temperature state may not have been fully reached after drilling, the temperature and salinity profiles seem to be scarcely disturbed. The average overall temperature gradient calculated from all temperature measurements is 28.5 K/km; remarkably lower than the temperature gradients found in other boreholes in the western Ross See and the Transantarctic Mountains. Anomalies in the salinity profiles at the beginning of each logging phase were no longer present at the end of the corresponding logging phase. This pattern indicates that drilling mud invaded the formation during drilling operations and flowed back into the borehole after drilling ceased. Thus, zones of temperature and salinity anomalies identify permeable zones in the formation and may be pathways for fluid flow. Radiogenic heat production, calculated from the radionuclide contents, is relatively low, with average values between 0.5 and 1.0 pW/m3. The highest values (up to 2 µW/m3) were obtained for the lower part of the Beacon Sandstone below 855 mbsf. The heat flow component due to radiogenic heat production integrated over the entire borehole is 0.7 mW/m2. Thermal conductivities range from 1.3 to 3 W/mK with an average value of 2.1 W/mK over the Tertiary section. Together with the average temperature gradient of 28.5 K/km this yields an average heat flow value of 60 mW/m2.

(Table 1) Summary of sampling points described and main petrographic, morphological and geochemical characteristics

Targeted sampling on the Dolgovskoy Mound (northern Shatsky Ridge) revealed the presence of spectacular laterally extensive and differently shaped authigenic carbonates. The sampling stations were selected based on sidescan sonar and profiler images that show patchy backscatter and irregular and discontinuous reflections in the near subsurface. The interpretation of acoustic data from the top part of the mound supports the seafloor observations and the sampling that revealed the presence of a complex subsurface plumbing system characterized by carbonates and gas. The crusts sampled consist of carbonate cemented layered hemipelagic sedimentary Unit 1 associated with several centimetres thick microbial mats. Three different carbonate morphologies were observed: (a) tabular slabs, (b) subsurface cavernous carbonates consisting of void chambers up to 20 cm**3 in size and (c) chimney and tubular conduits vertically oriented or forming a subhorizontal network in the subsurface. The methanogenic origin of the carbonates is established based on visual observations of fluids seepage structures, 13C depletion of the carbonates (d13C varying between -36.7 per mil and -27.4 per mil), and by thin carbonate layers present within the thick microbial mats. Laboratory experiments with a Hele-Shaw cell were conducted in order to simulate the gas seepage through contrasting grain size media present on the seafloor. Combined petrography, visual observations and sandbox simulations allowed a characterization of the dynamics and the structures of the plumbing system in the near subsurface. Based on sample observations and the experiments, three observed morphologies of authigenic carbonates are interpreted, respectively, as (a) Darcian porous flow through the finely laminated clayey/coccolith-rich layers, (b) gas accumulation chambers at sites where significant fluid escape was impeded by thicker clayey layers forming the laminated Unit1 and (c) focussed vertical fluid venting and subhorizontal migration of overpressured fluids released from (b). The Hele-Shaw cell experiments represent a promising tool for investigating shallow fluid flow pathways in marine systems.

Tab. 1: Trace-element concentrations in ore samples collected from mercury mines and deposits in southwestern Alaska

A belt of small but numerous mercury deposits extends for about 500 km in the Kuskokwim River region of southwestern Alaska. The southwestern Alaska mercury belt is part of widespread mercury deposits of the circum Pacific region that are similar to other mercury deposits throughout the world because they are epithermal with formation temperatures of about 200 °C, the ore is dominantly cinnabar with Hg-Sb-As±Au geochemistry, and mineralized forms include vein, vein breccias, stockworks, replacements, and disseminations. The southwestern Alaska mercury belt has produced about 1400 t of mercury, which is small on an international scale. However, additional mercury deposits are likely to be discovered because the terrain is topographically low with significant vegetation cover. Anomalous concentrations of gold in cinnabar ore suggest that gold deposits are possible in higher temperature environments below some of the Alaska mercury deposits. We correlate mineralization of the southwestern Alaska mercury deposits with Late Cretaceous and early Tertiary igneous activity. Our 40Ar/39Ar ages of 70 ±3 Ma from hydrothermal sericites in the mercury deposits indicate a temporal association of igneous activity and mineralization. Furthermore, we suggest that our geological ancl geochemical data from the mercury deposits indicate that ore fluids were generated primarily in surrounding sedimentary wall rocks when they were cut by Late Cretaceous and early Tertiary intrusions. In our ore genesis model, igneous activity provided the heat to initiate dehydration reactions and expel fluids from hydrous minerals and formational waters in the surrounding sedimentary wall rocks, causing thermal convection and hydrothermal fluid flow through permeable rocks and along fractures and faults. Our isotopic data from sulfide and alteration minerals of the mercury deposits indicate that ore fluids were derived from multiple sources, with most ore fluids originating from the sedimentary wall rocks.

Strike-slip faults mediate the rise of crustal-derived fluids and mud volcanism in the deep sea

We report on newly discovered mud volcanoes located at about 4500 m water depth 90 km west of the deformation front of the accretionary wedge of the Gulf of Cadiz, and thus outside of their typical geotectonic environment. Seismic data suggest that fluid flow is mediated by a >400-km-long strike-slip fault marking the transcurrent plate boundary between Africa and Eurasia. Geochemical data (Cl, B, Sr, 87Sr/86Sr, Delta18O, DeltaD) reveal that fluids originate in oceanic crust older than 140 Ma. On their rise to the surface, these fluids receive strong geochemical signals from recrystallization of Upper Jurassic carbonates and clay-mineral dehydration in younger terrigeneous units. At present, reports of mud volcanoes in similar deep-sea settings are rare, but given that the large area of transform-type plate boundaries has been barely investigated, such pathways of fluid discharge may provide an important, yet unappreciated link between the deeply buried oceanic crust and the deep ocean.

Granulometry, mineralogy and geochemistry of sediments from ODP Hole 169-1037B and Site 169-1038

Central Hill is in the northern part of the Escanaba Trough, which is a sediment-filled rift of southern Gorda Ridge. Central Hill is oriented north-south and is associated with extensive sulfide deposits. Hydrothermal alteration of sediment from Site 1038 was studied through analyses of mineralogy and the chemistry and oxygen isotopic compositions of one nearly pure clay sample. In addition, Site 1037 was drilled to establish the character of the unaltered sedimentary sequence away from the hydrothermal centers of the Northern Escanaba Trough Study Area (NESCA). Mineralogy of the clay-size fraction of turbiditic and hemipelagic sediments of Hole 1037B are predominantly quartz, feldspar, pyroxene, illite, chlorite, and smectite, representing continental-derived material. Cores from Hole 1038I, located within the area of Central Hill but away from known active vent areas, recovered minor amounts of chlorite/smectite mixed-layer clay in the fine fraction, indicating a low-temperature hydrothermal alteration. The 137.4-m-thick sediment section of Hole 1038G is located in an area of low-temperature venting. The uppermost sample is classified as chlorite/smectite mixed layer, which is underlain by chlorite as the dominant mineral. The lowermost deposits of Hole 1038G are also characterized by chlorite/smectite mixed-layer clay. In comparison to Hole 1038I, the mineralogic sequence of Hole 1038G reflects increased chloritization. Intensely altered sediment is almost completely replaced by hydrothermal chlorite in subsurface sediments of Hole 1038H. Alteration to chlorite is characterized by depletion in Na, K, Ti, Ca, Sr, Cs, and Tl and enrichment in Ba. Further, Eu depletion reflects a high-temperature plagioclase alteration. A chlorite 18O value of 2.6 indicates formation at a temperature of ~190°C. It is concluded that the authigenic chlorite in Hole 1038H formed by an active high-temperature fluid flow in the shallow subsurface.

U/Th systematics of authigenic carbonates and d18O records from Hydrate Ridge, off the coast of Oregon, USA

Uranium (U) concentrations and activity ratios (d234U) of authigenic carbonates are sensitive recorders of different fluid compositions at submarine seeps of hydrocarbon-rich fluids ("cold seeps") at Hydrate Ridge, off the coast of Oregon, USA. The low U concentrations (mean: 1.3 ± 0.4 µg/g) and high 234U values (165-317 per mil) of gas hydrate carbonates reflect the influence of sedimentary pore water indicating that these carbonates were formed under reducing conditions below or at the seafloor. Their 230Th/234U ages span a time interval from 0.8 to 6.4 ka and cluster around 1.2 and 4.7 ka. In contrast, chemoherm carbonates precipitate from marine bottom water marked by relatively high U concentrations (mean: 5.2 ± 0.8 µg/g) and a mean d234U ratio of 166 ± 3 per mil. Their U isotopes reflect the d234U ratios of the bottom water being enriched in 234U relative to normal seawater. Simple mass balance calculations based on U concentrations and their corresponding d234U ratios reveal a contribution of about 11% of sedimentary pore water to the bottom water. From the U pore water flux and the reconstructed U pore water concentration a mean flow rate of about 147 ± 68 cm/a can be estimated. 230Th/234U ages of chemoherm carbonates range from 7.3 to 267.6 ka. 230Th/234U ages of two chemoherms (Alvin and SE-Knoll chemoherm) correspond to time intervals of low sealevel stands in marine isotope stages (MIS) 2, 4, 5, 6, 7 and 8. This observation indicates that fluid flow at cold seep sites sensitively reflects pressure changes of the hydraulic head in the sediments. The d18OPDB ratios of the chemoherm carbonates support the hypothesis of precipitation during glacial times. Deviations of the chemoherm d18O values from the marine d18O record can be interpreted as to reflect temporally and spatially varying bottom water and/or vent fluid temperatures during carbonate precipitation between 2.6 and 8.6°C.