Search NASA⌕ Search

SEARCH · Search NASA

Results for “Air simulations”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Illinois Compressed Air Energy Storage

Compressed Air Storage Energy (CAES) is one of the few mid- technology readiness level (TRL) energy storage technologies that can address the long-duration infrastructure needed for dealing with variable electric output from renewable energy sources and be reliable backup source for replacing natural gas during supply interruptions. In CAES the goal is to capture and store compressed air in subsurface sedimentary strata when off-peak power is available, or there is a need for grid balancing. The stored high-pressure air is returned to the surface and used to power turbines during reductions in either renewable energy or supply issues with fossil fuels. The Illinois CAES project evaluates the feasibility of capturing surplus electrical energy from renewable sources and off-peak energy at a fossil fuel power plant at the University of Illinois Urbana - Champaign (UIUC) campus. The UIUC Abbott Power Plant uses natural gas and coal to generate electricity (capacity: 35 MWe by coal and 49 MWe by NG). UIUC receives additional electricity from on campus solar farm, and off-campus wind farm. Also, UIUC offsets electricity usage by integrating geothermal energy systems into building heating Also, UIUC offsets steam, hot and chilled water usage by integrating geothermal energy systems into building heating and cooling systems. Furthermore, the two UIUC solar farms (Solar Farm 1 is 21 acres and Solar Farm 2 is 54 acres) to generate 4.68 megawatts (MW) and 12.1 MW, respectively. Campus receives 8.6% of the wind-generated electricity from the Rail Splitter Wind Farm. The project objectives were to design an integrated system to 1) capture surplus electrical energy from renewable sources and the Abbott Power Plant using a CAES system, 2) store both the compressed air and the thermal heat generated by compression in the subsurface as part of an adiabatic system, 3) simulate the movement of the air and heat in the subsurface, 4) recover the compressed air and stored thermal heat to rotate turbine generators during sustained interruption due to weather events or fossil fuel disruptions.

03 NATURAL GAS↗
DOE OSTI2023

A framework for calibrating and validating an air loop dynamic model in an HVAC system in Modelica

The use of Modelica for simulating the dynamic behaviors of building heating, ventilation, and air conditioning (HVAC) systems has gained popularity. Calibration of a model that represents large and complex HVAC systems in Modelica involves the determination of hundreds of parameters using real-world operation data. Considering the coupling effects among various components, such a calibration process is complex and time-consuming. In this study, we propose a systematic framework to efficiently calibrate and validate a complex HVAC system model in Modelica. The framework includes strategies to solicit real-system operation data and to decouple the model. The goal of this calibration framework is to accurately and efficiently determine a set of Modelica model parameters that provide a good match between the simulated and real system behaviors. To demonstrate the validity of the framework, it was applied to calibrate the parameters of a Modelica model of the air loop subsystem of a real AHU-VAV system. The results show that the calibrated model can generate simulated results, including VAV air flow rate, outdoor air flow, and fan power consumption, that match well with the operational data of the real system. The coefficient of the variation of the root mean square error, CV(RMSE), of the air flow rate and power consumption are 14.6 % and 11.5 %, respectively. The results prove that the framework is valid and effective, and it can be used to calibrate other complex HVAC system Modelica models in the future.

Li, Yicheng↗
DOE OSTI2020

Peak load reduction and load shaping in HVAC and refrigeration systems in commercial buildings by using a novel lightweight dynamic priority-based control strategy

Reducing peak power demand in a building can reduce electricity expenses for the building owner and contribute to the efficiency and reliability of the electrical power grid. For the building owner, reduced expenses come from the reduction or elimination of peak power charges on electricity bills. For the power system operator, reducing peak power demand leads to a more predictable load profile and reduces stress on the electric grid system. Herein we present a computationally inexpensive, dynamic, and retrofit-deployable control strategy to effect peak load reduction and load shaping. The effectiveness of the control strategy is examined in a simulation with 80 air-conditioning units and 40 refrigeration units. The results show that a peak demand reduction of 60 kW can be achieved relative to peak demand in a typical set point–based approach. The proposed strategy was deployed in a gymnasium building with four rooftop HVAC units, where it showed over 15% peak demand (kW) reduction savings while maintaining or lowering energy consumption (in kilowatt-hours) relative to the set point–based thermostat controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Heat release surrogates for NH 3 /H 2 /N 2 –air premixed flames

The present study investigates the performance of NH, NH 2 , O-atom, and H-atom as heat release rate (HRR) surrogates for NH 3 /H 2 /N 2 –air premixed flame through simulations. The simulations are conducted across different pressures, reactant inlet temperatures, fuel blend compositions, and equivalence ratios. Cantera is used to simulate one-dimensional (1D) freely propagating flames to investigate the spatial correlations of the species with the HRR. PeleLMeX, a low-Mach direct numerical simulation (DNS) code with Adaptive Mesh Refinement (AMR), is used to simulate two-dimensional (2D) flame-vortex interactions to investigate the temporal correlations including stretch effects. Three different mechanisms (Jiang et al. 2020; Glarborg et al. 2018; Okafor et al. 2018) were considered in the 1D flame simulations, whereas only the Jiang mechanism was considered in the flame-vortex simulations. The HRR surrogate performance for the 2D flames is evaluated at two different locations: (1) the centerline and (2) the flame cusp. The cusp is defined as the region in the flame front with the greatest curvature and the centerline encounters the highest tangential strain rate. The 1D flame results suggest that, although there is not uniformly good spatial correlation for HRR across all flame conditions, NH is the best overall as a HRR surrogate for laminar flames. The 2D flame results, however, suggest that O-atom and H-atom have satisfactory temporal correlations at different conditions—the former for rich flames, the latter for high-pressure flames. Furthermore, these simulations provide guidance to experimental measurements of surrogate HRR markers in unsteady multi-dimensional flames using laser diagnostics to detect species such as NH, O-atom, and H-atom.

Ammonia↗
DOE OSTI2023

Cloud condensation nuclei activity of internally mixed particle populations at a remote marine free troposphere site in the North Atlantic Ocean

Here, this study reports results from research conducted at the Observatory of Mount Pico (OMP), 2225 m above mean sea level on Pico Island in the Azores archipelago in June and July 2017. We investigated the chemical composition, mixing state, and cloud condensation nuclei (CCN) activities of long-range transported free tropospheric (FT) particles. FLEXible PARTicle Lagrangian particle dispersion model (FLEXPART) simulations reveal that most air masses that arrived at the OMP during the sampling period originated in North America and were highly aged (average plume age > 10 days). We probed size-resolved chemical composition, mixing state, and hygroscopicity parameter ($\mathcal{κ}$) of individual particles using computer-controlled scanning electron microscopy with an energy-dispersive X-ray spectrometer (CCSEM-EDX). Based on the estimated individual particle mass from elemental composition, we calculated the mixing state index, $χ$. During our study, FT particle populations were internally mixed ($χ$ of samples are between 53% and 87%), owing to the long atmospheric aging time. We used data from a miniature Cloud Condensation Nucleus Counter (miniCCNC) to derive the hygroscopicity parameter, $\mathcal{κ}$ CCNC . Combining $\mathcal{κ}$ CCNC and FLEXPART, we found that air masses recirculated above the North Atlantic Ocean with lower mean altitude had higher $\mathcal{κ}$ CCNC due to the higher contribution of sea salt particles. We used CCSEM-EDX and phase state measurements to predict single-particle $\mathcal{κ}$ ($\mathcal{κ}$ CCSEM-EDX ) values, which overlap with the lower range of $\mathcal{κ}$ CCNC measured below 0.15 % SS. Therefore, CCSEM-EDX measurements can be useful in predicting the lower bound of $\mathcal{κ}$, which can be used in climate models to predict CCN activities, especially in remote locations where online CCN measurements are unavailable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Enhanced cooling tower technology for power plant efficiency increase and operating flexibility

The Final Report describes the results of the work completed during the “Enhanced Cooling Tower Technology for Power Plant Efficiency Increase and Operating Flexibility” project. The objectives of work were to develop and demonstrate the prototype (up to 100 kWth) of an economically viable “all weather” Sub-Dew Point Cooling Tower (SDPCT) with inlet air cooling and dehumidification including testing in laboratory-controlled environment for simulating various ambient conditions. The inlet air is cooled and dehumidified by a Pressure Dehumidifying System (PDHS) installed at the cooling tower inlet.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Godiva-IV Dosimetry Exercise 2022 (IER-538 CED4A Report)

This report presents the final results of IER 538, The International Nuclear Accident Dosimeter (NAD) Intercomparison Exercise with Godiva-IV. The purpose of the exercise, held in August of 2022, was to test participants’ NADs and dosimetry personnel to the dose estimation requirements established by ANSI/HPS-N13.3 Dosimetry for Criticality Accidents and DOE-STD-1098-2017, Radiological Control. Two prompt critical bursts of the Godiva critical assembly were used to simulate criticality accidents, and NADs from participating laboratories were placed at known locations around Godiva, mounted on either BOttle Manikin ABsorptions (BOMABs) phantoms (to simulate doses received by people) or plates (to simulate doses in free air). Similar exercises have been held in 2016 and 2018 using National Criticality Experiments Research Center’s (NCERC’s) Flattop and Godiva assemblies.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

IER-484 CED3b Report: AFRRI ER1 Dosimetry Characterization

This report presents the preliminary work for IER 484, the characterization of the leakage radia tion field from the TRIGA reactor at the Armed Forces Radiobiology Research Institute (AFRRI). The purpose of the dose characterization is to establish the reference dose needed to execute a NAD inter-comparison exercise. During the week of August 20th, a team from Lawrence Livermore National Laboratory (LLNL), Sandia National Laboratory (SNL), and Institut de Radioprotection et de Sûreté Nucléaire (IRSN) performed various measurements in AFRRI’s TRIGA ER1 (Exposure Room 1) with Passive Neutron Spectrometer (PNS) spheres and nuclear accident dosimeters (NADs) mounted on Bottle Manikin Absorption (BOMAB) phantoms (to simulate doses received by people) and additional NADs mounted on aluminum plates (to simulate doses in free air). Without Atomic Weapons Establishment (AWE), LLNL became responsible for measuring the neutron spectrum using PNS spheres. SNL is responsible for providing photon dose data. LLNL, SNL, and IRSN each deployed their own NADs, which will be compared to the reference dose measured. Final characterization of the dose field based on analyzed data will be presented in the CED4a report.

42 ENGINEERING↗
DOE OSTI2025

CaptureFOAM

CaptureFOAM provides utilities for fast and accurate direct air capture (DAC) contactor simulations, based on the OpenFOAM framework. At the core, CaptureFOAM is a suite of mass transfer and chemistry models built upon the OpenFOAM thin film model, which allows coarser mesh resolution and less constrictive time step requirements when compared to a traditional volume of fluid approach. This allows CaptureFOAM to conduct device-scale simulations for complex geometries including commercial packing materials and novel contactor shapes, matching experimental results in a wide variety of applications.

Kincaid, Kellis [Oak Ridge National Laboratory (OR↗
DOE OSTI2021

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗
DOE OSTI2021

Investigation of Thermal-Hydraulic Effects of Dry Storage Canister Helium Backfill Loss Using the Horizontal Dry Cask Simulator

A previous investigation produced data sets that can be used to benchmark the codes and best practices presently used to determine cladding temperatures and induced cooling air flows in modern horizontal dry storage systems. The horizontal dry cask simulator (HDCS) was designed to generate this benchmark data and add to the existing knowledge base. The objective of the previous HDCS investigation was to capture the dominant physics of a commercial dry storage system in a well-characterized test apparatus for a wide range of operational parameters. The close coupling between the thermal response of the canister system and the resulting induced cooling air flow rate was of particular importance. The previous investigation explored these parameters using helium backfill at 100 kPa and 800 kPa pressure as well as air backfill with a series of simulated decay heats. The helium tests simulated a horizontal dry cask storage system at normal storage conditions with either atmospheric or elevated backfill pressure, while the air tests simulated horizontal storage canisters following a complete loss of helium backfill, in which case the helium would be replaced by air. The present HDCS investigation adds to the previous investigation by exploring steady-state conditions at various stages of the loss of helium backfill from a horizontal dry cask storage system. This is achieved by using helium/air blends as a backfill in the HDCS and running a series of tests using various simulated decay heats to explore the effects of relative helium/air molar concentration on the thermal response of a simulated horizontal dry cask storage system. A total of twenty tests were conducted where the HDCS achieved steady state for various assembly powers, representative of decay heat. The power levels tested were 0.50, 1.00, 2.50, and 5.00 kW. All tests were run at 100 kPa vessel pressure. The backfill gases used in these tests are given in this report as a function of mole fraction of helium (He), balanced by air: 1.0, 0.9, 0.5, 0.1, and 0.0 He. Steady-state conditions (where the steady-state start condition is defined as where the change in temperature with respect to time for the majority of HDCS components is less than or equal to 0.3 K/h) were achieved for all test cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Comparing the cooling effectiveness of operationalisable urban surface combination scenarios for summer heat mitigation

Extreme summer heat in cities exacerbates the vulnerability of urban communities to heatwaves. Vegetative and reflective urban surfaces can help reduce urban heat. This study investigated the impacts of urban trees, green roofs and cool roofs on heat mitigation during average and extreme summer conditions in temperate oceanic Melbourne, Australia. We simulated the city climate using ‘The Air Pollution Model’ (TAPM) at a 1 km spatial resolution over 10 years, which according to our review of the literature, was the most prolonged period for simulation in Melbourne. During a widespread heatwave event, some of the tested scenarios with combined surface parameters could reduce the extreme values of the energy budget components- sensible heat, latent heat, and storage heat fluxes up to seasonal averages compared to the existing situation for Melbourne (control). The scenario with the highest (reasonable maximum) ground-level vegetation, green roofs, and cool roofs could reduce air temperatures up to 2.4 °C. The simulations suggest that a combined strategy with vegetative and high-albedo surfaces will deliver higher effectiveness with maximum cooling benefits and cost-effectiveness than individual strategies in cities. These results suggest the importance of collaborative strategic planning of urban surfaces to make cities healthier, sustainable, and liveable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Effect of Air Separations on Fast Pyrolysis Products for Forest Residue Feedstocks

This study investigates the intricate relationship between biomass preprocessing and pyrolysis product yields, employing the air classification technique for the treatment of loblolly pine residues with varying moisture content. A comprehensive exploration of the physicochemical properties of air-classified loblolly pine informs a sophisticated pyrolysis simulation model. Given the complex and multifaceted nature of biomass pyrolysis, operating across diverse temporal and spatial scales, a pyrolysis kinetics-based CFD–DEM simulation method is employed to predict product yields. Results showed that the elevated moisture content amplifies particle adhesiveness, necessitating augmented air velocities for effective separation, thereby influencing the efficiency of the separation process. While carbon and hydrogen contents exhibit relative stability across diverse moisture contents and blower frequencies, the oxygen content undergoes noticeable changes. For example, the oxygen contents were measured as 29.2 and 38.6 wt% in the light fraction of 30% moisture content sample at blower frequencies of 10 and 20 Hz, respectively. An intriguing finding emerges from pyrolysis simulation, indicating that a lower blower frequency in air classification moderately enhances bio-oil yield and significantly improves its quality, particularly in terms of water content. For instance, the water content in the bio-oil was about 1.5% and 10% in the heavy and light fractions, respectively from 10% moisture sample under 15 Hz blower frequency. In summary, a detailed understanding and strategic manipulation of critical material attributes in biomass through efficient fractionation techniques are imperative for advancing fast pyrolysis as a sustainable avenue for renewable energy and chemical production.

09 BIOMASS FUELS↗
DOE OSTI2025

The Effect of Updraft Entrainment on Convective Cell Deepening in Realistic Large-Eddy Simulations

Entrainment of surrounding cooler and drier air into convective updrafts is one of the key processes that influence deep convection initiation and growth. Numerous studies have investigated the effect of entrainment on isolated convective cloud growth in idealized simulations, but the importance of this effect in realistic conditions with many interacting convective clouds remains uncertain. We examine the impact of entrainment on the depth reached by convective clouds in realistic large-eddy simulations (LES) over central Argentina during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Cloudy updrafts and their associated properties are assigned to convective cells tracked with radar reflectivity signatures. Several thousand convective cells are tracked over two high convective available potential energy (CAPE) and two low CAPE cases that support cells of varying depths and intensities. Entrainment is calculated explicitly as the fluxes of air into the outer surface of each cloudy updraft. Single-predictor logistic regression models are used to determine the relative importance of updraft, near-updraft, and preconvective initiation atmospheric conditions in predicting whether convective cells become deep. We then build a multiple-predictor regression model pairing important updraft and meteorological metrics with fractional entrainment rate. The probability of cells transitioning to deep convection is most sensitive to ambient 600-hPa relative humidity (42% of total metric contribution to cloud depth predictability), followed by low-level CAPE (28%), cloud-base updraft width (19%), and fractional entrainment (11%). Thus, the initial width of the updraft along with potential buoyancy and its dilution through the midtroposphere collectively determine whether deep convection will result from shallower clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating the Feasibility of Measuring and Analyzing Strain Caused by Changes in Barometric Pressure in the Vadose Zone

Changes in barometric pressure propagate into the subsurface where they can affect water level measurements in wells and cause vertical strain. Previous strain sensors were developed to measure vertical strain, but they were limited to measurements at a single depth, making it difficult to evaluate strains from a migrating pressure wave. Technology known as Coherence-length-gated Microwave Photonics Interferometry (CMPI) uses optical fiber sensors to measure strain at multiple locations. It has the potential to detect strain caused by variations in barometric pressure at multiple depths in the subsurface. The goal of this study is to evaluate the feasibility of using this technology to record the changes in vertical strain in the vadose zone from small fluctuations in air pressure. Propagation of these air pressure fluctuations is affected by hydrogeological properties such as water content and permeability. A second goal of this study is to evaluate how these properties affect the air pressure distribution, strain, and gas-phase diffusivity with a long-term goal of using strain to monitor the vadose zone. The focus of this study is on a suite of laboratory experiments that used a sand-filled column with an open head space. The CMPI fiber was packaged and embedded along the axis of the column to measure the strain at multiple locations while air pressure transducers were installed through the wall of the column. Small periodic fluctuations in air pressure, similar to barometric pressure fluctuations, were created in the head space of the column using an audio driven speaker driven at 4 Hz using a sinusoidal signal created by a function generator. Then, the fluctuations propagated along the column where they were measured by the transducers and recorded as functions of time. Initial tests were conducted using dry sand, but then additional tests were conducted after injecting water into the column, which changed the hydrologic and mechanical properties. Another experiment was conducted after creating a thin barrier to air flow, which was designed to simulate a thin layer of saturated ground during rainfall. A final experiment was conducted by changing the grain size of the material to create a heterogeneity in the upper half of the column. The suite of six experiments was conducted to highlight effects of variations in water content and heterogeneities on air pressure and vertical strain. Numerical simulations were also created to evaluate the laboratory data and to provide a baseline analysis for the results. The numerical simulations used the governing equations to linear poroelasticity and two-phase flow with boundary conditions representing the experiments. The conceptual model for strain caused by barometric pressure recognizes that a fluctuating barometric pressure causes fluctuations in the air pore pressure that decrease in amplitude and lag in time with increase depth. The pressure distribution results in two different loads that cause strain: 1.) a mechanical load as the barometric pressure acts on the ground surface; and 2.) pressure loading in the pore space. The results show that the total vertical strain is a contribution of both of these loads, and this is the strain that was observed in the laboratory. The lab experiments show that both air pressure and strain vary as sinusoids that propagate with depth. The amplitude of the air pressure and strain decrease as a function of depth in all the conditions that were evaluate. For example, in the dry sand, the amplitude of the air pressure decreases from 25 Pa at the head space to 5 Pa at a depth of 70 cm, whereas the amplitude of the strain decreases from 0.15 to 0.06 µε over the same depth interval. In many cases the amplitudes decrease is approximately negative exponential functions of depth. When water was injected into the sand and the saturation increased, the amplitude of the air pressure increased in the partially saturated sand but decreased sharply to zero where the soil was saturated. The strain decreased with depth, but the transition from partially to fully saturated conditions had little effect on the strain profile, even though it had a major effect on the air pressure. This apparently is because strain is caused by variations in both water and air pressure. The fluctuating air pressure caused the water pressure to fluctuate, which caused strain throughout the column. Phase delay of the air pressure and the strain both increase as linear to bi-linear functions of depth. The phase delays indicate that pressure propagates at a velocity of 7 to 8 m/s in dry or partially saturated sand, but it drops by more than an order of magnitude in the vicinity of the capillary fringe. The phase delays indicate the strain propagates faster than the pressure. The data indicate that the strain propagates from 1.2 to more than 4x faster than the pressure. Switching from sand to silt in the column caused the velocity of the pressure to drop by roughly half for both pressures, but curiously, it had little effect on the velocity of the strain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗