Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volume of fluid”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

New Approaches to Inferences for Steep-Sided Domes on Venus

New mathematical approaches for the relaxation and emplacement of viscous lava domes are presented and applied to steep-sided domes on Venus. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry for two distinct scenarios. In the first scenario, dome relaxation is explored assuming a constant volume of fluid (i.e. lava) has been rapidly emplaced onto the surface. Cooling of lava is represented by a time-variable viscosity and singularities inherent in previous models for dome relaxation have been eliminated. At the onset of relaxation, bulk dynamic viscosities lie in the range between 10(exp 10) - 10(exp 16) Pa s, consistent with basaltic-andesite to rhyolitic compositions. Plausible relaxation times range from 5 to 5000 years, depending on initial lava viscosity. The first scenario, however, is only valid during the final stages of dome relaxation and does not consider the time taken for lava to be extruded onto the surface. In the second scenario, emplacement and growth of a steep-sided dome is considered when the volume of lava on the surface increases over time (i.e. time-variable volume approach). The volumetric flow rate may depend on an arbitrary power of the dome thickness, thus embracing Newtonian as well as other rheologies for describing terrestrial and planetary mass flows. The approach can be used to distinguish between basic flow rate models for fluid emplacement. The formalism results in radial expansion of a dome proportional to t(sup 1/2), consistent with the diffusive nature of the governing equation. The flow at the front is shown to thicken as the front advances for a constant rate of lava supply. Emplacement times are intimately correlated with the bulk rheology. Comparison of the theoretical profiles with the shape of a typical dome on Venus indicates that a Newtonian bulk rheology is most appropriate, consistent with prior studies. However, results here suggest a bulk dynamic viscosity of 10(exp 12) - 10(exp 13) Pa s and emplacement times of approximately 2-16 years. Both scenarios investigated give emplacement times significantly less than prior estimates and compositions consistent with basaltic andesite.

steep-sided domes↗
NASA NTRS2001

CFD Simulation of Liquid Rocket Engine Injectors

Detailed design issues associated with liquid rocket engine injectors and combustion chamber operation require CFD methodology which simulates highly three-dimensional, turbulent, vaporizing, and combusting flows. The primary utility of such simulations involves predicting multi-dimensional effects caused by specific injector configurations. SECA, Inc. and Engineering Sciences, Inc. have been developing appropriate computational methodology for NASA/MSFC for the past decade. CFD tools and computers have improved dramatically during this time period; however, the physical submodels used in these analyses must still remain relatively simple in order to produce useful results. Simulations of clustered coaxial and impinger injector elements for hydrogen and hydrocarbon fuels, which account for real fluid properties, is the immediate goal of this research. The spray combustion codes are based on the FDNS CFD code' and are structured to represent homogeneous and heterogeneous spray combustion. The homogeneous spray model treats the flow as a continuum of multi-phase, multicomponent fluids which move without thermal or velocity lags between the phases. Two heterogeneous models were developed: (1) a volume-of-fluid (VOF) model which represents the liquid core of coaxial or impinger jets and their atomization and vaporization, and (2) a Blob model which represents the injected streams as a cloud of droplets the size of the injector orifice which subsequently exhibit particle interaction, vaporization, and combustion. All of these spray models are computationally intensive, but this is unavoidable to accurately account for the complex physics and combustion which is to be predicted, Work is currently in progress to parallelize these codes to improve their computational efficiency. These spray combustion codes were used to simulate the three test cases which are the subject of the 2nd International Workshop on-Rocket Combustion Modeling. Such test cases are considered by these investigators to be very valuable for code validation because combustion kinetics, turbulence models and atomization models based on low pressure experiments of hydrogen air combustion do not adequately verify analytical or CFD submodels which are necessary to simulate rocket engine combustion. We wish to emphasize that the simulations which we prepared for this meeting are meant to test the accuracy of the approximations used in our general purpose spray combustion models, rather than represent a definitive analysis of each of the experiments which were conducted. Our goal is to accurately predict local temperatures and mixture ratios in rocket engines; hence predicting individual experiments is used only for code validation. To replace the conventional JANNAF standard axisymmetric finite-rate (TDK) computer code 2 for performance prediction with CFD cases, such codes must posses two features. Firstly, they must be as easy to use and of comparable run times for conventional performance predictions. Secondly, they must provide more detailed predictions of the flowfields near the injector face. Specifically, they must accurately predict the convective mixing of injected liquid propellants in terms of the injector element configurations.

Farmer, Richard↗
NASA NTRS1992

A finite-volume numerical method to calculate fluid forces and rotordynamic coefficients in seals

A numerical method to calculate rotordynamic coefficients of seals is presented. The flow in a seal is solved by using a finite-volume formulation of the full Navier-Stokes equations with appropriate turbulence models. The seal rotor is perturbed along a diameter such that the position of the rotor is a sinusoidal function of time. The resulting flow domain changes with time, and the time-dependent flow in the seal is solved using a space conserving moving grid formulation. The time-varying fluid pressure reaction forces are then linked with the rotor center displacement, velocity and acceleration to yield the rotordynamic coefficients. Results for an annular seal are presented, and compared with experimental data and other more simplified numerical methods.

Athavale, M. M.↗
NASA NTRSDate not supplied

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between liquid and vapor phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures, temperatures, and pressurant requirements agree with experimental data, with pressurant mass predictions within 16\% of observed values. This study highlights the importance of selecting appropriate turbulence models to accurately simulate complex flow and heat transfer phenomena during tank pressurization and expulsion. By enhancing the accuracy and reliability of CFD models for liquid hydrogen under cryogenic conditions, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Validation of a Two-Phase CFD Model for Autogenous Pressurization and Expulsion

This paper presents a two-phase computational fluid dynamics (CFD) model for simulating autogenous pressurization and expulsion in a cryogenic liquid hydrogen storage tank, utilizing a volume-of-fluid (VOF) approach combined with the kinetics-based Schrage equation to accurately capture the phase boundary and interfacial heat, mass, and momentum transfer between vapor and liquid phases. The model is validated against experimental data from NASA's K-site tank facility, specifically experimental case 225, which involved autogenous pressurization and controlled expulsion of liquid hydrogen. Various turbulence models are evaluated to assess their influence on model accuracy. The CFD simulations successfully replicate key thermodynamic behaviors observed during the experiments, including pressure evolution, temperature profiles, and phase-change dynamics at the vapor-liquid interface. The predicted tank pressures and temperatures agree well with experimental data, whereas the pressurant mass predictions are within 16% of observed values. This study emphasizes the importance of selecting appropriate turbulence models to accurately simulate the complex flow and heat transfer phenomena during tank autogenous pressurization and expulsion. By improving the accuracy and reliability of CFD models for these processes, this research contributes to developing efficient cryogenic propellant management strategies for future space missions.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗
NASA NTRS2004

Numerical, Analytical, Experimental Study of Fluid Dynamic Forces in Seals: Description of Incompressible Fluid Seal Codes ICYL and IFACE - Volume 4

The objectives of the program were to develop computational fluid dynamics (CFD) codes and simpler industrial codes for analyzing and designing advanced seals for air-breathing and space propulsion engines. The CFD code SCISEAL is capable of producing full three-dimensional flow field information for a variety of cylindrical configurations. An implicit multidomain capability allow the division of complex flow domains to allow optimum use of computational cells. SCISEAL also has the unique capability to produce cross-coupled stiffness and damping coefficients for rotordynamic computations. The industrial codes consist of a series of separate stand-alone modules designed for expeditious parametric analyses and optimization of a wide variety of cylindrical and face seals. Coupled through a Knowledge-Based System (KBS) that provides a user-friendly Graphical User Interface (GUI), the industrial codes are PC based using an OS/2 operating system. These codes were designed to treat film seals where a clearance exists between the rotating and stationary components. Leakage is inhibited by surface roughness, small but stiff clearance films, and viscous pumping devices. The codes have demonstrated to be a valuable resource for seal development of future air-breathing and space propulsion engines.

Liang, Anita D.↗
NASA NTRSDate not supplied

Validation of Ullage Collapse Due to Violent Lateral Slosh

Understanding the coupling between thermodynamics and propellant sloshing is important for cryogenic upper stages, especially during the lift-off phase and in missions including multiple engine restarts, as the fluid dynamic condition of the propellant is mainly influenced by various flight maneuvers. This paper describes the development and validation of a Volume of Fluid (VOF)-based CFD solver, Loci/STREAM-VOF, to the ullage collapse under violent later slosh conditions. Loci/STREAM-VOF is shown to be able to capture detailed slosh dynamics, including slosh frequency, non-linear slosh damping, and violent surface breakup due to forced acceleration. The comparison of surface shape and small drops formed during large amplitude slosh with experimental imagery shows excellent agreement throughout the sloshing period. Comparison of temperature histories in the ullage, at the gas-liquid interface, and inside liquid with the experiment conducted at JAXA also show very good agreement. The simulation points to a stratified temperature field inside the liquid near the interface, critical to the ullage pressure after the collapse. The predicted pressure collapse matches well with the experimental value due to sloshing. Parametric sensitivity studies suggest that high-fidelity Loci/STREM-VOF simulations can play an essential role in predicting the actual thermodynamics of ullage collapse when following the recommended simulation methodology.

H Q Yang↗
NASA NTRS2009

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗
NASA NTRS1998

Stress, Aging and Thirst

After growth during adolesence, total body water decreases progressively with aging from 65% of body weight to about 53% of body weight in the 70th decade; a majority of the loss occurs from the extracellular volume, from 42% to about 25%, respectively. Cellular volume also reaches equilibrium in the 70th decade at about 25% of body weight. Various stresses such as exercise, heat and attitude exposure, ad prior dehydration attenuate voluntary fluid intake (involuntary dehydration). Voluntary fluid intake appears to decrease with aging (involuntary dehydration in this sense aging can be considered as a stress. Kidney function and muscle mass (80% water) decrease somewhat with aging, and voluntary fluid intake (thirst) is also attenuated. Thirst is stimulated by increasing osmolality (hypernatremia) of the extracellular fluid and by decreased extracellular volume (mainly plasma volume) which act to increase intracellular fluid volume osmolality to activiate drinking. The latter decreases fluid compartment osmolality which ' It terminates drinking. However, this drinking mechanism seems to be attenuated with aging such that increasing plasma osmolality no longer stimulates fluid intake appropriately. Hypernatremia in the elderly has been associated all too frequently with greater incidence of bacterial infection and increased mortality. Involuntary dehydration can be overcome in young men by acclimation to an intermittent exercise-in-heat training program. Perhaps exercise training in the elderly would also increase voluntary fluid intake and increase muscle mass to enhance retention of water.

Greenleaf, John E.↗
NASA NTRS1982

Computer simulations of postural change, water immersion and bedrest - An integrative approach for understanding the spaceflight response

Mathematical modeling techniques were used to simulate the fluid electrolyte (F-E) responses during gravity unloading. It is shown that the response to weightlessness can best be understood by separately examining the acute (hours to days) and chronic (days to weeks) phases, and assuming the presence of normal, although complex, feedback regulatory processes. Headward shifts of fluid are shown to be primarily responsible for acute body losses of extracellular F-E. Losses of body water are closely related to the volume of fluid shifts from the legs. A diuresis is predicted within the first several hours of hypogravity, and this may be obscured by a reduced F-E intake; on Skylab, early F-E losses occurred primarily by deficit intake.

Leonard, J. I.↗
NASA NTRSDate not supplied

CFD Modeling of Cryogenic Chilldown in a Complex Channel under Normal and Low Gravity Conditions

Future NASA architectures have baselined cryogenic propulsion systems as well as cryogenic fluid management to support lunar missions and ultimately to support future missions to Mars. These missions will require chilling hardware down prior to engine restart as well as chilling lines and tanks prior to transferring and refueling these propulsion elements in orbit. In lieu of expensive tests conducted on-orbit, accurate predictive computational models of these chilldown processes can be used to reduce system and propellant mass as well as mission risk. To gain confidence in these computational models, appropriate anchoring and validation to experimental data in a relevant environment needs to be performed. Recent ground and sub-orbital flight experiments conducted by the Japan Aerospace Exploration Agency (JAXA) investigated chilldown of a complex channel resembling a turbopump bearing cavity at low flow rates. This work presents Computational Fluid Dynamics (CFD) model development of the chilldown experiment employing two-phase flow boiling models available in commercial CFD software STAR-CCM+ using Volume of Fluid (VOF) and the traditional Euler-Euler multiphase flow solvers. Comparisons of the numerical and experimental results under normal and low-gravity conditions are presented. An assessment of solid wall temperatures and phase distribution yielded important insights into multiphase solver choice, dependence on gravity environment, and challenges associated with cryogenic flow boiling prediction and validation.

CFD↗
NASA NTRS1989

Zero-Gravity Vortex Vent and PVT Gaging System

Space Station and satellite reservicing will require the ability to vent gas on orbit from liquid supply or storage tanks and to gage liquid quantity under microgravity conditions. In zero gravity, (zero-g) the vortex vent is capable of venting gas from a tank of liquid containing gas randomly distributed as bubbles. The concept uses a spinning impeller to create centrifugal force inside a vortex tube within a tank. This creates a gas pocket and forces the liquid through a venturi and back into the tank. Gas is then vented from the gas pocket through a liquid detector and then out through an exhaust port. If the liquid detector senses liquid in the vent line, the fluid is directed to the low-pressure port on the venturi and is returned to the tank. The advantages of this system is that it has no rotating seals and is compatible with most corrosive and cryogenic fluids. A prototype was designed and built at the NASA Johnson Space Center and flown on the KC-135 zero-g aircraft. During these test flights, where microgravity conditions are obtained for up to 30 sec, the prototype demonstrated that less than 0.10 percent of the volume of fluid vented was liquid when the tank was half full of liquid. The pressure volume temperature (PVT) gaging system is used in conjunction with the vortex vent to calculate the amount of liquid remaining in a tank under microgravity conditions. The PVT gaging system is used in conjunction with the vortex vent to gage liquid quantity in zero or low gravity. The system consists of a gas compressor, accumulator, and temperature and pressure instrumentation. To measure the liquid in a tank a small amount of gas is vented from the tank to the compressor and compressed into the accumulator. Pressure and temperature in the tank and accumulator are measured before and after the gas transfer occurs. Knowing the total volume of the tank, the volume of the accumulator, the volume of the intermediate lines, and initial and final pressures and temperatures, the mass of the gas leaving the tank is equated to the mass of the gas entering the accumulator. The volume of liquid remaining in the tank is calculated using the ideal gas law.

Downey, M. G.↗
NASA NTRSDate not supplied

Modeling of Particle Transport and Impact in Transient, Multiphase Oxidizer-Rich Environments Using Computational Fluid Dynamics and Lagrangian Particle Tracking

Particle impacts from foreign object debris (FOD) in an oxidizer-rich environment have been shown to be a credible and significant risk for ignition and failure of components and vehicles. The risk of particle ignition is dependent on many factors, including particle impact velocity/energy and fluid composition (i.e. fluid phase, species presence, and state). Prediction of particle transport and particle impacts in relevant environments, including transient, multiphase flows, is important for assessing the risk and avoiding repeated expensive system level tests. The capability to model both the fluid dynamics and particle transport in complex three-dimensional geometries and multiphase environments has been developed at Marshall Space Flight Center (MSFC) in the propulsion fluid dynamic branch (ER42) by utilizing computational fluid dynamics (CFD) simulations and Lagrangian particle tracking. Using CFD tools, including density-based Loci/CHEM with a cavitation model and pressure-based Loci/STREAM with a volume of fluids model, MSFC/ER42 has provided particle impact predictions to multiple projects for oxidizer compatibility assessments. Simulations of transient, multiphase environments such as liquid slugs flowing through tubing and cavitating flow in valves have been completed alongside one-way coupled Lagrangian particle tracking to predict particle impacts.

M.V. Fischels↗
NASA NTRS2016

Modeling Ullage Dynamics of Tank Pressure Control Experiment during Jet Mixing in Microgravity

A CFD model for simulating the fluid dynamics of the jet induced mixing process is utilized in this paper to model the pressure control portion of the Tank Pressure Control Experiment (TPCE) in microgravity1. The Volume of Fluid (VOF) method is used for modeling the dynamics of the interface during mixing. The simulations were performed at a range of jet Weber numbers from non-penetrating to fully penetrating. Two different initial ullage positions were considered. The computational results for the jet-ullage interaction are compared with still images from the video of the experiment. A qualitative comparison shows that the CFD model was able to capture the main features of the interfacial dynamics, as well as the jet penetration of the ullage.

active control↗
DOE OSTI2021

A Numerical Analysis of The Pool Boiling Process in Flooded Evaporators Using Different Refrigerants

The pool boiling process can be observed in several energy conversion processes including commercial and industrial refrigeration, industrial air-cooling operations, and power generation. The process becomes more involved when pool boiling in a tube bundle is considered. In the current study, a numerical model is developed to predict the key performance parameters of a flooded evaporator while considering a range of working fluids. A kettle reboiler configuration was considered, and a performance model was developed to account for boiling on individual tubes, merging of vapor bubbles, and movement under gravity. A volume of fluid (VOF) model was used to deal with the different phases in the simulation. Various fluids were considered in this study including HFE-7000, HFE-7300, and water. The trajectories of the bubbles were tracked, and the resulting information has been summarized in physical measurable quantities.

Zhang, Mingkan↗
DOE OSTI2022

Hardness Prediction by Incorporating Heat Transfer and Molten Pool Fluid Flow in a Multi-pass, Multi-layer Weld for Onsite Repair of Grade 91 steel

Introduction: In the current fleet of fossil-fired power plants, creep strength enhanced ferritic steels (CSEF) are used to sustain the harsh service conditions. Enhanced properties of Grade 91 steel result from tempered martensite with a fine distribution of MX and M23C6 carbides. Grade 91 steel is subjected to onsite welding repair to remedy their degradation due to extreme service condition. Knowledge of weld repairability of these steels, such as as-welded hardness distribution, is essential to establishing sound repair procedures. Experimental trial and error tests can consume a lot of time as many welding variables need to be studied. For numerical modelling, most of the multi-pass multi-layer models are based on finite element method, which are limited to solve the heat conduction equation and ignore convective heat transfer due to melt flow. Moreover, the mesh has to be pre-built based on a known or assumed weld cross-section geometry. These finite element based models thus have limited predictive capability as defects are not considered and nugget size are pre-assumed. This research aims at developing a thermal and microstructure evolution model incorporating molten pool dynamics in a multi-pass multi-layer material deposition to predict the as-welded hardness distribution. Technical Approach: All the thermal, physical, and metallurgical properties of Grade 91 as a function of temperature are collected from the literature and inputted into the thermo-fluid model based on Flow-3D, a computational fluid dynamics software. A multi-pass, multi-layer material deposition is simulated where the melting of filler wire into the molten pool is directly considered based on the volume of fluid (VOF) method. The flow behaviour of the molten pool is used to understand the formation of deposition geometry and defects. The temperature profiles during the multi-pass, multi-layer welding are calculated. The results computed using the new model are compared against the experimental data of fusion zone geometry and thermal cycles. Hardness prediction in the heat-affected zone (HAZ) are made using Johnson-Mehl-Avrami (JMA) equation for solid-state phase transformation kinetics. The JMA parameters are extracted from the experimental data available in the literature. For comparison, a standard finite element heat conduction model is also developed to predict the thermal cycles and hardness distribution in the multi-pass, multi-layer weld. Expected Result: Results obtained using the molten pool dynamic simulation versus the finite element heat conduction model are compared. Specifically, the effects of convective heat transfer on the accuracy of the calculated thermal history, bead shape and size, and HAZ hardness distribution are examined.

Modeling, Grade 91 steel, Multi-pass multi-layer d↗
NASA NTRS2014

Search for Fluid Inclusions in a Carbonaceous Chondrite Using a New X-Ray Micro-Tomography Technique Combined with FIB Sampling

Early solar system aqueous fluids are preserved in some H chondrites as aqueous fluid inclusions in halite (e.g., [1]). Although potential fluid inclusions are also expected in carbonaceous chondrites [2], they have not been surely confirmed. In order to search for these fluid inclusions, we have developped a new X-ray micro-tomography technique combined with FIB sampling and applied this techniqu to a carbanaceous chondrite. Experimental: A polished thin section of Sutter's Mill meteorite (CM) was observed with an optical microscope and FE-SEM (JEOL 7001F) for chosing mineral grains of carbonates (mainly calcite) and sulfides (FeS and ZnS) 20-50 microns in typical size, which may have aqueous fluid inclusions. Then, a "house" similar to a cube with a roof (20-30 microns in size) is sampled from the mineral grain by using FIB (FEI Quanta 200 3DS). Then, the house was atached to a thin W-needle by FIB and imaged by a SR-based imaging microtomography system with a Fresnel zone plate at beamline BL47XU, SPring-8, Japan. One sample was imaged at two X-ray energies, 7 and 8 keV, to identify mineral phases (dual-enegy microtomography: [3]). The size of voxel (pixel in 3D) was 50-80 nm, which gave the effective spatial resolution of approx. 200 nm. A terrestrial quartz sample with an aqueous fluid inclusion with a bubble was also examined as a test sample by the same method. Results and discussion: A fluid inclusion of 5-8 microns in quartz was clearly identified in a CT image. A bubble of approx. 4 microns was also identified as refraction contrast although the X-ray absorption difference between fluid and bubble is small. Volumes of the fluid and bubble were obtained from the 3D CT images. Fourteen grains of calcite, two grains of iron sulfide and one grain of (Zn,Fe)S were examined. Ten calcite, one iron sulfide and one (Zn,Fe)S grains have inclusions >1 micron in size (the maximum: approx. 5 microns). The shapes are spherical or irregular. Tiny inclusions (<1 micron) are also present in all the grains examined. These results show that mineral grains have more inclusions than expected from 2D observations. The X-ray absorption of the inclusions shows that they are not solid inclusions. No bubbles were observed inside, indicating that we cannot determine whether they are really aqueous fluids or merely voids. One calcite grain has an inclusion approx. 2 microns in size, which seems to have a bubble and a tiny solid daughter crystal inside (three-phase inclusion). As we know the exact 3D position of the inclusion, we will anlyze the inclusion by SIMS after freezing the sample as has been done for a halite sample [3]. The present technique is useful for finding small inclusions not only in carbonaceous chondrites but also for terrestrial materials.

Tsuchiyama, A.↗