Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BINARY 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Transient Effects in Planar Solidification of Dilute Binary Alloys

The initial transient during planar solidification of dilute binary alloys is studied in the framework of the boundary integral method that leads to the non-linear Volterra integral governing equation. An analytical solution of this equation is obtained for the case of a constant growth rate which constitutes the well-known Tiller's formula for the solute transient. The more physically relevant, constant ramping down temperature case has been studied both numerically and analytically. In particular, an asymptotic analytical solution is obtained for the initial transient behavior. A numerical technique to solve the non-linear Volterra equation is developed and the solution is obtained for a family of the governing parameters. For the rapid solidification condition, growth rate spikes have been observed even for the infinite kinetics model. When recirculating fluid flow is included into the analysis, the spike feature is dramatically diminished. Finally, we have investigated planar solidification with a fluctuating temperature field as a possible mechanism for frequently observed solute trapping bands.

Mazuruk, Konstantin↗
DOE OSTI2015

RVA: A Plugin for ParaView 3.14

RVA is a plugin developed for the 64-bit Windows version of the ParaView 3.14 visualization package. RVA is designed to provide support in the visualization and analysis of complex reservoirs being managed using multi-fluid EOR techniques. RVA, for Reservoir Visualization and Analysis, was developed at the University of Illinois at Urbana-Champaign, with contributions from the Illinois State Geological Survey, Department of Computer Science and National Center for Supercomputing Applications. RVA was designed to utilize and enhance the state-of-the-art visualization capabilities within ParaView, readily allowing joint visualization of geologic framework and reservoir fluid simulation model results. Particular emphasis was placed on enabling visualization and analysis of simulation results highlighting multiple fluid phases, multiple properties for each fluid phase (including flow lines), multiple geologic models and multiple time steps. Additional advanced functionality was provided through the development of custom code to implement data mining capabilities. The built-in functionality of ParaView provides the capacity to process and visualize data sets ranging from small models on local desktop systems to extremely large models created and stored on remote supercomputers. The RVA plugin that we developed and the associated User Manual provide improved functionality through new software tools, and instruction in the use of ParaView-RVA, targeted to petroleum engineers and geologists in industry and research. The RVA web site (http://rva.cs.illinois.edu) provides an overview of functions, and the development web site (https://github.com/shaffer1/RVA) provides ready access to the source code, compiled binaries, user manual, and a suite of demonstration data sets. Key functionality has been included to support a range of reservoirs visualization and analysis needs, including: sophisticated connectivity analysis, cross sections through simulation results between selected wells, simplified volumetric calculations, global vertical exaggeration adjustments, ingestion of UTChem simulation results, ingestion of Isatis geostatistical framework models, interrogation of joint geologic and reservoir modeling results, joint visualization and analysis of well history files, location-targeted visualization, advanced correlation analysis, visualization of flow paths, and creation of static images and animations highlighting targeted reservoir features.

Source record↗
DOE OSTI2017

A collision model for grain-resolving simulations of flows over dense, mobile, polydisperse granular sediment beds

We present a collision model for phase-resolved Direct Numerical Simulations of sediment transport that couple the fluid and particles by the Immersed Boundary Method. Typically, a contact model for these types of simulations comprises a lubrication force for particles in close proximity to another solid object, a normal contact force to prevent particles from overlapping, and a tangential contact force to account for friction. Our model extends the work of previous authors to improve upon the time integration scheme to obtain consistent results for particle-wall collisions. Furthermore, we account for polydisperse spherical particles and introduce new criteria to account for enduring contact, which occurs in many sediment transport situations. This is done without using arbitrary values for physically-defined parameters and by maintaining the full momentum balance of a particle in enduring contact. Lastly, we validate our model against several test cases for binary particle-wall collisions as well as the collective motion of a sediment bed sheared by a viscous flow, yielding satisfactory agreement with experimental data by various authors.

42 ENGINEERING↗
DOE OSTI2023

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗
DOE OSTI2024

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Probing the Incompressibility of Nuclear Matter at Ultrahigh Density through the Prompt Collapse of Asymmetric Neutron Star Binaries

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density K max to within tens of percent. This otherwise inaccessible measure of K max can potentially reveal the presence of hyperons or quarks inside neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Pore-Scale Simulation of Intragranular Diffusion: Effects of Incomplete Mixing on Macroscopic Manifestations

Diffusive mass transfer into and out of intragranular micropores ("intragranular diffusion") plays an important role in the large-scale transport of some groundwater contaminants. We are interested in understanding the combined effect of pore-scale advection and intragranular diffusion on solute transport at the effective porous medium scale. To study this problem, we have developed a three-dimensional pore-scale model of fluid flow and solute transport that incorporates diffusion into and out of intragranular pore spaces. Our model is based on the Smoothed Particle Hydrodynamics (SPH) simulation method, which represents fluid and solid phases by a mesh-free particle discretization. In the pore spaces, fluid flow is simulated by discretizing the Navier-Stokes equations using the SPH approach. Solute transport is represented by advection, diffusion within the fluid phase, and diffusion between the fluid and solid phases. Our model is implemented on large-scale parallel computing hardware, allowing us to simulate millions of computational particles and represent fully three-dimensional systems of pores and grains with arbitrarily complex physical geometry. We have used this model system to perform numerical experiments using various model porous media systems, which allows us to draw comparisons between macroscopic measures computed from the pore-scale simulations (such as breakthrough curves) and those predicted by macroscopic formulations that assume complete mixing over the representative volume. In this paper we present results of 3D simulations of pore-scale flow and transport, including cases with and without intragranular diffusion, in two model porous media, one with randomly-packed uniform spherical grains and a second with randomly-packed spheres drawn from a binary grain size distribution. Breakthrough curves were computed from the 3D simulations at various transport distances. Comparable breakthrough curves were computed using 1D macroscopic models with parameters determined independently to the degree possible. Based on comparisons of the pore-scale and macroscopic model results, we draw two primary conclusions. First, non-Fickian behavior is persistent and ubiquitous at the scales considered, and most cases are better represented by a multi-rate mass transfer model even when there is no distinct secondary porosity (i.e., no intragranular diffusion). This suggests that diffusive mass transfer processes between preferential flow paths and relatively immobile zones within the primary porosity may have significant impact on transport, particular in low-concentration tails. Second, the application of mass transfer rate parameters based on an assumption of well-mixed concentrations at the pore-scale tend to overestimate the amount of mass transfer that occurs in heterogeneous pore geometries in which preferential flow leads to incomplete pore-scale lateral mixing.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2025

The Multiple Paths to Merger of Unequal-mass Black Hole Binaries in the Disks of Active Galactic Nuclei

The accretion disks that power active galactic nuclei (AGN) are thought to house populations of stars and compact objects; after forming binaries these compact objects may merge, begetting gravitational waves (GWs) such as those detected by LIGO and VIRGO. We present a comprehensive study of the early evolution of binaries within AGN disks as their orbits are influenced by the surrounding gas, focusing on eccentric and unequal-mass binaries. Nearly equal-mass binaries behave similarly to their equal-mass counterparts: Prograde binaries inspiral, albeit somewhat slowly, and have their eccentricities damped; retrograde binaries inspiral ∼2–3 times faster than their prograde counterparts, and those with near-equal masses are driven quickly toward near-unity eccentricities. However, the primaries in retrograde binaries with mass ratios of m 2 /m 1 ≲ 0.4 experience significantly weaker headwinds and retain substantial accretion disks that help damp binary eccentricities, slowing binary inspirals. Additionally, we find that while accretion drives prograde binaries toward equal masses thanks to the exchange of material between the primary and secondary accretion disks, retrograde binaries are driven slowly toward more extreme mass ratios. Prograde binaries, and generally those with low mass ratios, likely accrete for multiple e-folding timescales before merger. On the other hand, high-mass-ratio retrograde binaries may merge before accreting substantially, potentially approaching merger with detectable eccentricity. Future ground-based GW observatories, with their broader frequency coverage, should be particularly useful for studying these populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Comprehensive experimental investigation on biomass-glass beads binary fluidization: A data set for CFD model validation

Fluidization behavior of biomass and glass beads binary mixtures in a bubbling fluidized bed was experimentally investigated. Mixtures containing different mass fraction of Loblolly Pine white wood and glass beads were fluidized at different fluidization velocities. The particle properties were characterized in a QICPIC that uses a dynamic image processing method to measure both particle size and sphericity. The minimum fluidization velocity was determined using the pressure drop method. An image processing method was developed to capture the dynamic expanded bed height at a very high frequency. The effect of biomass mass fraction and inlet gas velocity on mixing and segregation behavior was studied and analyzed through pressure drop measurements. Pressure drop fluctuations and expanded bed height fluctuations via fast Fourier transform were analyzed and compared. Here, the complete and accurate experimental data reported in this study could provide a benchmark data set for various computational fluid dynamics models validation, calibration, and identification.

42 ENGINEERING↗
DOE OSTI2022

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗
DOE OSTI2025

pH-SWITCHABLE DYNAMIC BINARY COMPLEXES AS VISCOSITY MODIFYING AGENTS

The study of supramolecular structures has been gaining a lot of momentum in recent years due to their potential application in several interesting areas such as medicine, catalysis, energy, oil and gas, etc. Hydraulic fracturing is one area where the ongoing demand for novel viscosity modifiers has prompted researchers to explore various facets of chemistry. In this study, a dynamic binary complex (DBC) solution is presented, with the potential to be useful in the hydraulic fracturing domain. It has a supramolecular structure that is formed by the self-assembly of oleic acid and diethylenetriamine into an elongated network under alkaline conditions. With less than 2 wt% of constituents in aqueous solution, a viscous gel that exhibits high viscosities even under shear was formed. It was sensitive to pH and salinity, and its zero-shear viscosity could be tuned by a factor of ~280 by changing the pH. Furthermore, its viscous properties were pronounced in the presence of salt. Sand settling tests revealed its potential to hold up sand particles for extended periods of time, longer than 24 hours. In conclusion, this DBC solution system has potential to be utilized as a smart pH-switchable hydraulic fracturing fluid that can be prepared using seawater.

02 PETROLEUM↗
NASA NTRS1998

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗
NASA NTRS2002

PCS: The First Fluid Physics Payload on ISS

The Physics of Colloids in Space (PCS) experiment was accommodated within International Space Station (ISS) EXpedite the PRocessing of Experiments to Space Station (EXPRESS) Rack 2 and was remotely operated from early June 2001 until February 2002 from NASA Glenn Research Center's Telescience Support Center in Cleveland, Ohio and from a remote site at Harvard University in Cambridge, Massachusetts. PCS is an experiment conceived by Professor David A. Weitz of Harvard University (the Principal Investigator), focusing on the behavior of three different classes of colloid mixtures. The sophisticated light scattering instrumentation comprising PCS is capable of color imaging, and dynamic and static light scattering from 11 to 169 degrees, Bragg scattering over the range from 10 to 60 degrees, and laser light scattering at low angles from 0.3 to 6.0 degrees. The PCS instrumentation performed remarkably well, demonstrating a flexibility that enabled experiments to be performed that had not been envisioned prior to launch. While on-orbit, PCS accomplished 2400 hours of science operations, and was declared a resounding success. Each of the eight sample cells worked well and produced interesting and important results. Crystal nucleation and growth and the resulting structures of two binary colloidal crystal alloys were studied, with the long duration microgravity environment of the ISS facilitating extended studies on the growth and coarsening characteristics of the crystals. In another experiment run, the de-mixing of the colloid-polymer critical-point sample was studied as it phase-separates into two phases, one that resembles a gas and one that resembles a liquid. This process was studied over four decades of length scale, from 1 micron to 1 centimeter, behavior that cannot be observed in this sample on Earth because sedimentation would cause the colloids to fall to the bottom of the cell faster than the de-mixing process could occur. Similarly, the study of gelation and aging of another colloid-polymer sample, the colloid-polymer gel, also provided valuable information on gelation mechanisms, as did investigations on the extremely the low concentration silica and polystyrene fractal gel samples.

Doherty, M.↗
NASA NTRS2012

NASA Tech Briefs, November 2012

The topics include: Visual System for Browsing, Analysis, and Retrieval of Data (ViSBARD); Time-Domain Terahertz Computed Axial Tomography NDE System; Adaptive Sampling of Time Series During Remote Exploration; A Tracking Sun Photometer Without Moving Parts; Surface Temperature Data Analysis; Modular, Autonomous Command and Data Handling Software with Built-In Simulation and Test; In-Situ Wire Damage Detection System; Amplifier Module for 260-GHz Band Using Quartz Waveguide Transitions; Wideband Agile Digital Microwave Radiometer; Buckyball Nucleation of HiPco Tubes; FACT, Mega-ROSA, SOLAROSA; An Integrated, Layered-Spinel Composite Cathode for Energy Storage Applications; Engineered Multifunctional Surfaces for Fluid Handling; Polyolefin-Based Aerogels; Adjusting Permittivity by Blending Varying Ratios of SWNTs; Gravity-Assist Mechanical Simulator for Outreach; Concept for Hydrogen-Impregnated Nanofiber/Photovoltaic Cargo Stowage System; DROP: Durable Reconnaissance and Observation Platform; Developing Physiologic Models for Emergency Medical Procedures Under Microgravity; Spectroscopic Chemical Analysis Methods and Apparatus; Low Average Sidelobe Slot Array Antennas for Radiometer Applications; Motion-Corrected 3D Sonic Anemometer for Tethersondes and Other Moving Platforms; Water Treatment Systems for Long Spaceflights; Microchip Non-Aqueous Capillary Electrophoresis (MicronNACE) Method to Analyze Long-Chain Primary Amines; Low-Cost Phased Array Antenna for Sounding Rockets, Missiles, and Expendable Launch Vehicles; Mars Science Laboratory Engineering Cameras; Seismic Imager Space Telescope; Estimating Sea Surface Salinity and Wind Using Combined Passive and Active L-Band Microwave Observations; A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing; Scalable SCPPM Decoder; QuakeSim 2.0; HURON (HUman and Robotic Optimization Network) Multi-Agent Temporal Activity Planner/Scheduler; MPST Software: MoonKommand

Source record↗
DOE OSTI2019

A Fokker–Planck approach to a moment closure for mixing in variable-density turbulence

We develop a theory for the cascade mixing terms in a moment closure approach to binary active scalar mixing in variable-density turbulence. To address the variable-density complications we apply, as a principle and constraint, the conservation of the probability density function (PDF) through a Fokker–Planck equation with bounded sample space whose attractor is the beta PDF with skewness. Mixing is related to a single-point PDF as a realisability principle to provide mathematically rigorous expressions for the small scale statistics in terms of largescale moments. The problem of the unknown small-scale mixing is replaced with the determination of the drift and diffusion terms of a Fokker–Planck equation in a beta-PDF-convergent stochastic process. We find that realisability of a beta-convergent process requires the mixing time-scale ratio, taken as a constant in passive scalar mixing, to be a function of the mean mass fraction, mean fluid density, the Atwood number, the density-volume correlation and moments of the density field. We develop and compare the new model with direct numerical simulations data of non-stationary homogeneous variable-density turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Impact of induced seismic events on seal integrity, Texas Gulf Coast

Recent publications have suggested that large-scale CO 2 injection could trigger earthquakes and that even small- to moderate-sized earthquakes may threaten the seal integrity of the injection zone, and potentially damage buildings and other surface structures. In this study, we compared seal thickness to estimated fault displacement due to a single hypothetical seismic event in a selected area of the Texas Gulf Coast comprising an offshore strip of state waters along two Texas counties. To evaluate the slip generated by a single seismic event, we compiled well log information on shale/sand sequences and seismic information on fault geometric characteristics of a section of Lower Miocene age. The section is thousands of feet thick and is overlain and underlain by marine shales (Amph. B and Anahuac, respectively) that are relatively easy to correlate between wells. The Amph. B. shale is the secondary and ultimate seal for all injection intervals in the Lower Miocene. Given its thickness, no realistic seismic event or small series of seismic events will offset it significantly. However, this may not be true of smaller local primary seals. An analysis of geophysical logs of a total of 71 wells yielded a total of 2,871 sand / shale binary intervals. An analysis of the dedicated 3D seismic survey counted 723 fault traces at five roughly horizontal horizons within the Lower Miocene Fault displacement estimated using the product of the fault length times an uncertain multiplier coefficient assumed to follow a triangular distribution with a 10 -3 to 10 -5 range and a mode of 8 × 10 -5 . We then compared estimated single-event fault displacements to seal thicknesses by means of a Monte-Carlo analysis. Only 1.8% of thickness/displacement pairs display a displacement greater than 20% of the seal thickness. Only 0.26% of the pairs result in a displacement of half the seal thickness and only 0.05% of thickness/displacement pairs result in a clear seal rupture. The next step was to compare the magnitude of the event generated by such a displacement to documented magnitudes of “large” earthquakes generated by waterflooding and fluid disposal. Based on this analysis, we conclude that seismicity that may arise from CO 2 injection appears not to be a serious complication for CO 2 storage integrity, at least in the Gulf Coast area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Oxygen transport properties estimation by classical trajectory–direct simulation Monte Carlo

Coupling direct simulation Monte Carlo (DSMC) simulations with classical trajectory calculations is a powerful tool to improve predictive capabilities of computational dilute gas dynamics. The considerable increase in computational effort outlined in early applications of the method can be compensated by running simulations on massively parallel computers. In particular, Graphics Processing Unit acceleration has been found quite effective in reducing computing time of classical trajectory (CT)-DSMC simulations. The aim of the present work is to study dilute molecular oxygen flows by modeling binary collisions, in the rigid rotor approximation, through an accurate Potential Energy Surface (PES), obtained by molecular beams scattering. The PES accuracy is assessed by calculating molecular oxygen transport properties by different equilibrium and non-equilibrium CT-DSMC based simulations that provide close values of the transport properties. Comparisons with available experimental data are presented and discussed in the temperature range 300–900 K, where vibrational degrees of freedom are expected to play a limited (but not always negligible) role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗