Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulator data”

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 1,333 records · Page 74

Multilevel Techniques for Compression and Reduction of Scientific Data-Quantitative Control of Accuracy in Derived Quantities

Although many compression algorithms are focused on preserving pointwise values of the data, application scientists are generally more concerned with derived quantities. Equally well, the user may even be willing to accept a high level of lossiness in the compression provided that the compressed data respect certain invariants, such as mass conservation. In the current work, we develop a mathematical framework and techniques that enable data to be adaptively compressed while maintaining a specified tolerance on a class of user-prescribed quantities. Here, the algorithm is used to augment the functionality of the data reduction package MGARD developed in previous work and the functionality is illustrated by a range of application including data from computational simulation of autocatalytic reaction simulation, turbulent combustion simulation, experimental data obtained from magnetic confinement fusion experiment, and simulation of turbulent flow along a rectangular channel. In each case, we consider one or more relevant quantities of interest and reduce the data so as to preserve these quantities.

97 MATHEMATICS AND COMPUTING↗

NuMI Beam Muon Monitor Data Analysis and Simulation for Improved Beam Monitoring

The NuMI muon monitors (MMs) are a very important diagnostic tool for monitoring the stability of the neutrino beam used by the NOvA experiment at Fermilab. The goal of our study is to maintain the quality of the MM signal and to establish the correlations between the neutrino and muon beam profile. This study could also inform the LBNF decision on the beam diagnostic tools. We report on the progress of beam scan data analysis (beam position, spot size, and magnetic horn current scan) and comparison with the simulation outcomes.

43 PARTICLE ACCELERATORS↗

Detailed high-resolution three-dimensional simulations of OMEGA separated reactants inertial confinement fusion experiments

Here, we present results from the comparison of high-resolution three-dimensional (3D) simulations with data from the implosions of inertial confinement fusion capsules with separated reactants performed on the OMEGA laser facility. Each capsule, referred to as a “CD Mixcap,” is filled with tritium and has a polystyrene (CH) shell with a deuterated polystyrene (CD) layer whose burial depth is varied. In these implosions, fusion reactions between deuterium and tritium ions can occur only in the presence of atomic mix between the gas fill and shell material. The simulations feature accurate models for all known experimental asymmetries and do not employ any adjustable parameters to improve agreement with experimental data. Simulations are performed with the RAGE radiation-hydrodynamics code using an Implicit Large Eddy Simulation (ILES) strategy for the hydrodynamics. We obtain good agreement with the experimental data, including the DT/TT neutron yield ratios used to diagnose mix, for all burial depths of the deuterated shell layer. Additionally, simulations demonstrate good agreement with converged simulations employing explicit models for plasma diffusion and viscosity, suggesting that the implicit sub-grid model used in ILES is sufficient to model these processes in these experiments. In our simulations, mixing is driven by short-wavelength asymmetries and longer-wavelength features are responsible for developing flows that transport mixed material towards the center of the hot spot. Mix material transported by this process is responsible for most of the mix (DT) yield even for the capsule with a CD layer adjacent to the tritium fuel. Consistent with our previous results, mix does not play a significant role in TT neutron yield degradation; instead, this is dominated by the displacement of fuel from the center of the implosion due to the development of turbulent instabilities seeded by long-wavelength asymmetries. Through these processes, the long-wavelength asymmetries degrade TT yield more than the DT yield and thus bring DT/TT neutron yield ratios into agreement with experiment. Finally, we present a detailed comparison of the flows in 2D and 3D simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of a Smart Signal Detection Method for Cyclic Voltammetry via Artificial Neural Intelligence

The electrorefiner (ER) is the heart of pyroprocessing technology which is a high-temperature method for separating uranium from Used Nuclear Fuel (UNF). It is important to improve this technology with respect to nuclear materials detection and accountability. Artificial Neural Intelligence (ANI) is a novel data analysis and simulation method that can be applied to electrochemical data sets. This computational code, which has been performed using the commercial software MATLAB, can be trained to generalize adequate electrical current and potential simulated data sets for the unseen data with a high accuracy of prediction. For this purpose, a massive collection of cyclic voltammetry (CV) data sets by Hoover (2014), for 0.5, 1, 2.5, and 5 wt% of zirconium chloride in LiCl-KCl molten salt with different scan rates at 773 K has been considered. The computer is trained via ANI to predict the unseen data after providing suitable hidden layers and validation numbers. In addition, this work can trace the CV plot for a blind condition by interpolating between two simulated data set. The different hidden layers with various neurons (from 5 to 30) at several validation numbers (from 5 to 30) has been studied and the average percent error between experimental and theoretical data for 0.5 wt% with 200 and 450 mV/s has been calculated. Preliminary results demonstrate that if the number of hidden layers increases from one to three, the average error falls down from 44% to around 8%. The best condition which gives a minimum average percent error has been discovered and the simulated CV graph has been compared with the experimental data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transfer learning of neural surrogates on multifidelity groundwater simulations

Multifidelity data used in the paper published in Advances in Water Resources 206 (2025) 105140, https://doi.org/10.1016/j.advwatres.2025.105140 The code used to process the data is openly available on GitHub at https://github.com/Model-Reduction-and-UQ-Group/Transfer_Learning_K_reconstruction Computationally inexpensive surrogates of process-based models, such as deep neural networks, enable ensemble-based computations used in risk assessment, data assimilation, etc. However, generation of large datasets required to train a neural network can be as expensive as the ensemble simulations themselves. We ameliorate this challenge by using data from multifidelity (MF) groundwater simulations and transfer learning (TL) to reduce data generation costs while maintaining model accuracy. As a computational example, we train a deep convolutional neural network (CNN) to reconstruct permeability fields from saturation maps derived from a multiphase flow model. Starting with very low- and low-fidelity data generated on increasingly coarse meshes, we pretrain the CNN, followed by output-layer training and fine-tuning using only a limited number of high-fidelity samples. We demonstrate the surrogate’s robustness when interpreting low-quality inputs—such as interpolated maps or data affected by noise—which has strong implications for the applicability in practical hydrogeological scenarios. This multilevel MF-TL strategy achieves a favorable trade-off between computational efficiency and predictive accuracy, significantly outperforming high-fidelity-only approaches under the same computational budget.

Chiofalo, Alessia [University of Bologna] (ORCID:0↗

Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV

Improved jet energy scale corrections, based on a data sample corresponding to an integrated luminosity of 19.7 fb(-1) collected by the CMS experiment in proton-proton collisions at a center-of-mass energy of 8 TeV, are presented. The corrections as a function of pseudorapidity η and transverse momentum p(T) are extracted from data and simulated events combining several channels and methods. They account successively for the effects of pileup, uniformity of the detector response, and residual data-simulation jet energy scale differences. Further corrections, depending on the jet flavor and distance parameter (jet size) R, are also presented. The jet energy resolution is measured in data and simulated events and is studied as a function of pileup, jet size, and jet flavor. Typical jet energy resolutions at the central rapidities are 15–20% at 30 GeV, about 10% at 100 GeV, and 5% at 1 TeV. The studies exploit events with dijet topology, as well as photon+jet, Z+jet and multijet events. Several new techniques are used to account for the various sources of jet energy scale corrections, and a full set of uncertainties, and their correlations, are provided.The final uncertainties on the jet energy scale are below 3% across the phase space considered by most analyses (p(T),>30 GeV and 0| η| <5.). In the barrel region (| η| <1.3) an uncertainty below 1% for p(T)>30 GeV is reached, when excluding the jet flavor uncertainties, which are provided separately for different jet flavors. A new benchmark for jet energy scale determination at hadron colliders is achieved with 0.32% uncertainty for jets with \pt of the order of 165–330\GeV, and | η| <0.8.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Light hydrocarbon and noble gas migration as an analogue for potential CO 2 leakage: numerical simulations and field data from three hydrocarbon systems

Abstract Sites where anthropogenic CO 2 captured from industrial sources is stored in deep geological formations for climate change mitigation are required to show secure retention of the injected CO 2 . Monitoring, reporting, and verification (MRV) plans are needed to indicate that no CO 2 release has occurred. We explored the degree to which direct comparison between a surface anomaly and reservoir geochemistry using various geochemical parameters can be used for attribution. We used data collected on light hydrocarbons and noble gases throughout the sedimentary column at three CO 2 enhanced oil recovery (EOR) sites to understand the processes that may cause fluid evolution. Light hydrocarbon and noble gases were sampled from reservoirs, gas‐bearing intervals above reservoirs, and groundwater. Vertical profiles indicated that lighter components are relatively enriched during migration (i.e. chromatographic separation). Static and numerical models were designed to simulate episodic gas migration and geochemical alteration of these geochemical parameters from solubility and sorption. The effects of hydrocarbon solubility were minimal (Bernard ratio changes within 5.2%) although field data were within the range of expected alterations from sorption. Forward models of CO 2 migration and noble gas interactions showed that CO 2 stripping causes an enrichment of crustal noble gases. In areas where natural fluxes of CO 2 from depth are non‐existent, the occurrence of crustal noble gas signature may distinguish fugitive CO 2 from the reservoir from natural near‐surface sources, and could be considered to explain apparent fluid anomalies. © 2019 Society of Chemical Industry and John Wiley & Sons, Ltd.

Anderson, Jacob↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Strong Lensing Parameter Estimation on Ground-Based Imaging Data Using Simulation-Based Inference

Current ground-based cosmological surveys, such as the Dark Energy Survey (DES), are predicted to discover thousands of galaxy-scale strong lenses, while future surveys, such as the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will increase that number by 1-2 orders of magnitude. The large number of strong lenses discoverable in future surveys will make strong lensing a highly competitive and complementary cosmic probe. To leverage the increased statistical power of the lenses that will be discovered through upcoming surveys, automated lens analysis techniques are necessary. We present two Simulation-Based Inference (SBI) approaches for lens parameter estimation of galaxy-galaxy lenses. We demonstrate the successful application of Neural Posterior Estimation (NPE) to automate the inference of a 12-parameter lens mass model for DES-like ground-based imaging data. We compare our NPE constraints to a Bayesian Neural Network (BNN) and find that it outperforms the BNN, producing posterior distributions that are for the most part both more accurate and more precise; in particular, several source-light model parameters are systematically biased in the BNN implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Assessment of geophysical monitoring methods for detection of brine and CO 2 leakage in drinking water aquifers

Here, we evaluated and compared the effectiveness of four surface-based geophysical monitoring methods to downhole pressure and chemical methods to detect brine and CO 2 leakage in underground sources of drinking water overlying a CO 2 storage reservoir. This assessment uses synthetic monitoring data generated from 400 simulated aquifer impact data sets. Here, the six monitoring techniques can detect impacted groundwater once 20,000 tonnes of CO 2 have leaked into the drinking water resource. Geophysical methods are most effective at detecting shallow plumes. Although downhole monitoring methods outperform geophysical methods in detecting deep plumes, geophysical methods may help reduce the false negative during the post injection site care because they detect impacted groundwater where downhole sensors may be absent.

58 GEOSCIENCES↗

Comparison of residual dose rates from MARS and DORIAN simulations with experimental data from the SINBAD database

This benchmarking study aims to compare residual dose rate predictions from the MARS and DORIAN codes with experimental measurements taken from samples irradiated at CERN's CERF facility. The research evaluates the residual dose rates of these samples over varied cooling intervals. Additionally, it delves into the influence of parameter choices, such as specific physics processes and energy thresholds, on the production and transport of particles and their effects on simulation outcomes. While both MARS and DORIAN largely align with experimental findings, discrepancies were observed in the copper sample. This prompted an in-depth examination of the elevated dose rates from MARS.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analysis of the Impact of Radiation on Cameras and Robotic Perception - 19242

Robotic systems are increasingly used for inspection, maintenance, and disaster assistance for nuclear power plants. These robots utilize increasingly high levels of autonomy to operate in environments that are unsafe for humans to operate the robot in or require more precise control over the robot than a human can offer. In these cases, the robot must include perception and computation components to autonomously sense and navigate through the environment to fulfill its mission. This difficult task is complicated by radiation interfering with sensors' operation. The impact of this interference is largely unexplored, as existing works on radiation interference on sensors focus on estimating the incident radiation, rather than how it impacts a robotic system using the sensor. In this work, we seek to quantify and describe the impact of radiation on cameras, and how the radiation affects robotic performance. In this paper we study cameras that are suitable for small, lightweight unmanned aerial vehicles (UAVs) that are increasingly deployed to fly in confined areas in emergencies, where the surrounding environment is unsafe for workers to operate. These small UAVs lack the payload capacity to carry significant radiation shielding, so any onboard sensors will be affected by radiation. In the case of cameras, this radiation appears as noise in the camera imagery. This noise will interfere with the UAV's operation, by obscuring obstacles and interfering with localization and optical flow algorithms that detect detailed features in the environment. To date, there is a lack of information on how the radiation noise in the camera image impacts the overall robot system performance Existing work focuses on estimating the amount of incident radiation by measuring the corruption in the images. However, these works suffer several shortcomings for robotic systems. First, a robot must be able to successfully navigate its environment during a successful mission, and merely quantifying the amount of radiation noise in an image does not adequately explain how it may interfere with higher level navigation and perception tasks. Second, radiation estimation procedures rely on focusing the camera on a known and static background. This is unsuitable for robotic exploration tasks, where the robot frequently moves around an unknown environment that may be dynamically changing. Our work contains two contributions that address these shortcomings. First, we conduct a detailed analysis of the noise radiation causes in cameras by exposing cameras that are commonly used on UAVs [e.g., small form factor complementary metal-oxide-semiconductor (CMOS) cameras] to varying levels of radiation that may be present in emergency scenarios. From this exposure, we quantify not only how much noise is present, but quantify its appearance. By quantifying the appearance and amount of noise in a camera image, researchers can design mitigation strategies for the corruption, so that their robotic systems can operate in these challenging environments. Second, we demonstrate how to leverage this experimental data in robotic simulations. Using data from the first set of experiments allows us to test and verify the robots in simulation. Using simulators allows us to verify the robots will perform as expected during critical times without needing expensive test equipment or damaging valuable robot prototypes. Other researchers using robots for exploration of irradiated environments will benefit from our research by incorporating similar noise models in their research. 1 Southwest Research Institute{sup R} is a registered trademark of Southwest Research Institute in the United States and/or other countries. This will increase the reliability of robotic systems prior to deployment in emergency situations, which will in turn increase their effectiveness and speed the development of robots that would otherwise need expensive and time-consuming test procedures to validate. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗

Experimental validation of the TOPAS Monte Carlo system for passive scattering proton therapy

Purpose: TOPAS (TOol for PArticle Simulation) is a particle simulation code recently developed with the specific aim of making Monte Carlo simulations user-friendly for research and clinical physicists in the particle therapy community. The authors present a thorough and extensive experimental validation of Monte Carlo simulations performed with TOPAS in a variety of setups relevant for proton therapy applications. The set of validation measurements performed in this work represents an overall end-to-end testing strategy recommended for all clinical centers planning to rely on TOPAS for quality assurance or patient dose calculation and, more generally, for all the institutions using passive-scattering proton therapy systems. Methods: The authors systematically compared TOPAS simulations with measurements that are performed routinely within the quality assurance (QA) program in our institution as well as experiments specifically designed for this validation study. First, the authors compared TOPAS simulations with measurements of depth-dose curves for spread-out Bragg peak (SOBP) fields. Second, absolute dosimetry simulations were benchmarked against measured machine output factors (OFs). Third, the authors simulated and measured 2D dose profiles and analyzed the differences in terms of field flatness and symmetry and usable field size. Fourth, the authors designed a simple experiment using a half-beam shifter to assess the effects of multiple Coulomb scattering, beam divergence, and inverse square attenuation on lateral and longitudinal dose profiles measured and simulated in a water phantom. Fifth, TOPAS’ capabilities to simulate time dependent beam delivery was benchmarked against dose rate functions (i.e., dose per unit time vs time) measured at different depths inside an SOBP field. Sixth, simulations of the charge deposited by protons fully stopping in two different types of multilayer Faraday cups (MLFCs) were compared with measurements to benchmark the nuclear interaction models used in the simulations. Results: SOBPs’ range and modulation width were reproduced, on average, with an accuracy of +1, −2 and ±3 mm, respectively. OF simulations reproduced measured data within ±3%. Simulated 2D dose-profiles show field flatness and average field radius within ±3% of measured profiles. The field symmetry resulted, on average in ±3% agreement with commissioned profiles. TOPAS accuracy in reproducing measured dose profiles downstream the half beam shifter is better than 2%. Dose rate function simulation reproduced the measurements within ∼2% showing that the four-dimensional modeling of the passively modulation system was implement correctly and millimeter accuracy can be achieved in reproducing measured data. For MLFCs simulations, 2% agreement was found between TOPAS and both sets of experimental measurements. The overall results show that TOPAS simulations are within the clinical accepted tolerances for all QA measurements performed at our institution. Conclusions: Our Monte Carlo simulations reproduced accurately the experimental data acquired through all the measurements performed in this study. Thus, TOPAS can reliably be applied to quality assurance for proton therapy and also as an input for commissioning of commercial treatment planning systems. This work also provides the basis for routine clinical dose calculations in patients for all passive scattering proton therapy centers using TOPAS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A program to generate simulated radioxenon beta–gamma data for concentration verification and validation and training exercises

Abstract Several hundred simulated radioxenon beta-gamma data files were developed to assist in evaluating the performance and results from radioxenon concentration calculation analysis at the International Data Center (IDC) and other National Data Centers (NDC). PNNL developed a Beta-Gamma Simulator (BGSim) that incorporated GEANT-modeled data sets from radioxenon decay chains, as well as functionality to use nuclear detector-acquired data sets to create new beta-gamma spectra with varying amounts of background, 133Xe, 131mXe, 133mXe, 135Xe, and 222Rn and its decay products. The program has been implemented on a web-based applications platform and allows the user to create very specific data sets that incorporate most of the operational parameters for the current beta-gamma systems deployed in the International Monitoring System (IMS) and the On-site Inspection (OSI) equipment. After an initial beta-gamma simulations program was developed, additional uses began to be identified for the program output: training sets of two-dimensional spectra for data analysts at the IDC and other NDC, spectra for exercises such as the Integrated Field Exercise 2014 (IFE14) held in Jordan at the Dead Sea, and testing new analysis methods and algorithms

custom datasets↗

Understanding and Modeling Lossy Compression Schemes on HPC Scientific Data

Scientific simulations generate large amounts of floating-point data, which are often not very compressible using the traditional reduction schemes, such as deduplication or lossless compression. The emergence of lossy floating-point compression holds promise to satisfy the data reduction demand from HPC applications; however, lossy compression has not been widely adopted in science production. We believe a fundamental reason is that there is a lack of understanding of the benefits, pitfalls, and performance of lossy compression on scientific data. In this paper, we conduct a comprehensive study on state-of-the-art lossy compression, including ZFP, SZ, and ISABELA, using real and representative HPC datasets. Our evaluation reveals the complex interplay between compressor design, data features and compression performance. The impact of reduced accuracy on data analytics is also examined through a case study of fusion blob detection, offering domain scientists with the insights of what to expect from fidelity loss. Furthermore, the trial and error approach to understanding compression performance involves substantial compute and storage overhead. To this end, we propose a sampling based estimation method that extrapolates the reduction ratio from data samples, to guide domain scientists to make more informed data reduction decisions.

Lu, Tao↗

Effects of Confinement and Pressure on Vibrational Behavior of Nano-confined Propane

Fluids confined in nanopores exhibit significant deviations in their structure and dynamics from the bulk behavior. Although phase, structural and diffusive behaviors of confined fluids have been investigated and reported extensively, confinement effects on the vibrational properties are less understood. Here, we study the vibrational behavior of propane confined in 1.5 nm nanopores of MCM-41-S using inelastic neutron scattering (INS) and molecular dynamics (MD) simulations. Vibrational spectra have been obtained from INS data as functions of temperature and pressure. At ambient pressure, a strong quasielastic signal observed in the INS spectrum at 80 K suggests that confined propane remains liquid below the bulk phase melting point of 85 K. The quasielastic signal is heavily suppressed when either the pressure is increased to 1 kbar, or the temperature is lowered to 30 K, indicating solidification of pore-confined propane. Confinement in MCM-41-S pores results in a glass-like state of propane that exhibits a relatively featureless low-energy vibrational spectrum compared to the bulk crystalline propane. Increasing the pressure to 3 kbar results in hardening of the intermolecular and methyl torsional modes. The INS data are used for estimating the isochoric specific heat of confined propane, which is compared with the specific heat of bulk propane reported in literature. Data from MD simulations is used to calculate the vibrational power spectra that agree qualitatively with the experimental data. Simulation data also suggests a reduction in the structural ordering (positional, orientational and intramolecular) of propane under confinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

PFLOTRAN-E4D: A parallel open source PFLOTRAN module for simulating time-lapse electrical resistivity data

Time-lapse electrical resistivity tomography (ERT) is finding increased application for remotely monitoring processes occurring in the near subsurface in three-dimensions (i.e. 4D monitoring). However, there are few codes capable of simulating the evolution of subsurface resistivity and corresponding tomographic measurements arising from a particular process, particularly in parallel and with an open source license. Herein we describe and demonstrate an electrical resistivity tomography module for the PFLOTRAN subsurface simulation code, named PFLOTRAN-E4D. The PFLOTRAN-E4D module operates in parallel using a dedicated set of compute cores in a master-slave configuration. At each time step, the master processes receives subsurface states from PFLOTRAN, converts those states to bulk electrical conductivity, and instructs the slave processes to simulate a tomographic data set. The resulting multi-physics simulation capability enables accurate feasibility studies for ERT imaging, the identification of the ERT signatures that are unique to a given process, and facilitates the joint inversion of ERT data with hydrogeological data for subsurface characterization. PFLOTRAN-E4D is demonstrated herein using a field study of stage-driven groundwater/river water interaction ERT monitoring along the Columbia River, Washington, USA. Results demonstrate the complex nature of changes subsurface electrical conductivity, in both the saturated and unsaturated zones, arising from water table changes and from river water intrusion into the aquifer. The results also demonstrate the sensitivity of surface based ERT measurements to those changes over time. PFLOTRAN-E4D is available with the PFLOTRAN development version with an open-source license at https://bitbucket.org/pflotran/pflotran-dev .

Environmental Molecular Sciences Laboratory↗