Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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 667 records · Page 37

Synthesis of the land carbon fluxes of the Amazon region between 2010 and 2020

The Amazon is the largest continuous tropical forest in the world and plays a key role in the global carbon cycle. Human-induced disturbances and climate change have impacted the Amazon carbon balance. Here we conduct a comprehensive synthesis of existing state-of-the-art estimates of the contemporary land carbon fluxes in the Amazon using a set of bottom-up methods (i.e., dynamic vegetation models and bookkeeping models) and a top-down inversion (atmospheric inversion model) over the Brazilian Amazon and the whole Biogeographical Amazon domain. Over the whole biogeographical Amazon region bottom-up methodologies suggest a small average carbon sink over 2010-2020, in contrast to a small carbon source simulated by top-down inversion (2010-2018). However, these estimates are not significantly different from one another when accounting for their large individual uncertainties, highlighting remaining knowledge gaps, and the urgent need to reduce such uncertainties. Nevertheless, both methodologies agreed that the Brazilian Amazon has been a net carbon source during recent climate extremes and that the south-eastern Amazon was a net land carbon source over the whole study period (2010-2020). Overall, our results point to increasing human-induced disturbances (deforestation and forest degradation by wildfires) and reduction in the old-growth forest sink during drought.

54 ENVIRONMENTAL SCIENCES↗

Genome-Resolved Metagenomics and Detailed Geochemical Speciation Analyses Yield New Insights into Microbial Mercury Cycling in Geothermal Springs

Geothermal systems emit substantial amounts of aqueous, gaseous, and methylated mercury, but little is known about microbial influences on mercury speciation. Here, we report results from genome-resolved metagenomics and mercury speciation analysis of acidic warm springs in the Ngawha Geothermal Field (< 55°C, pH < 4.5), Northland Region, Aotearoa New Zealand. Our aim was to identify the microorganisms genetically equipped for mercury methylation, demethylation, or Hg(II) reduction to volatile Hg(0) in these springs. Dissolved total and methylated mercury concentrations in two adjacent springs with different mercury speciation ranked among the highest reported from natural sources (250 to 16,000 ng liter -1 and 0.5 to 13.9 ng liter -1 , respectively). Total solid mercury concentrations in spring sediments ranged from 1,274 to 7,000 μg g -1 . In the context of such ultrahigh mercury levels, the geothermal microbiome was unexpectedly diverse and dominated by acidophilic and mesophilic sulfur- and iron-cycling bacteria, mercury- and arsenic-resistant bacteria, and thermophilic and acidophilic archaea. By integrating microbiome structure and metagenomic potential with geochemical constraints, we constructed a conceptual model for biogeochemical mercury cycling in geothermal springs. The model includes abiotic and biotic controls on mercury speciation and illustrates how geothermal mercury cycling may couple to microbial community dynamics and sulfur and iron biogeochemistry.

59 BASIC BIOLOGICAL SCIENCES↗

TCO Analysis Approach and Regional Analysis of dWPT for Class 8 Tractors

Dynamic Wireless Power Transfer (dWPT) is a method by which battery electric vehicles (BEVs) can charge their battery while traveling on the road without the need for a physical conductive connection to the power source. dWPT has been proposed as a strategy to enable a reduction in vehicle battery capacity and associated mass and cost. In this slide deck presented at the EVs@Scale Consortium - Wireless Power Transfer Pillar Deep-Dive Meeting on November 11th, 2023, NREL provides results from an evaluation of dWPT using data from Class 8 tractors driving in the Atlanta Metro Area. NREL selected data for archetypal days representing local, regional, and long-haul trips, defined according to trip length, that included travel on primary roadways. EVI-InMotion (Electric Vehicle Infrastructure - InMotion), a systems planning and optimization tool developed at NREL, was used to evaluate dWPT performance assuming dWPT charging on 120 road segments for a total roadway lane distance of 2,365 miles. The EVI-InMotion results and representative day drive cycles were analyzed with NREL's T3CO (Transportation Technology Total Cost of Ownership) tool to estimate the total cost of ownership (TCO) for scenarios comprising two model years - 2030 and 2040 - and two technology progress cases. TCO was calculated for diesel, fuel cell electric, BEVs with batteries sized assuming no dWPT capabilities, and 200kWh BEVs with dWPT installed. This analysis finds that en-route stationary charging frequency and downtime when not on electrified roadways are the main contributors to TCO for the dWPT vehicles and that these vehicles can achieve cost parity with FCEVs at low electricity costs. Based on the scenario assumptions used here, low electricity costs would further help the cost parity with diesel vehicles in regional and long-haul cases due to stationary fueling downtime. This presentation also concludes that key factors affecting the parity potential of dWPT-capable vehicles include more extensive dWPT road coverage, higher en-route charging power, less expensive power batteries, and higher hydrogen or diesel fuel costs.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

Optimizing Renewable Ammonia Production for a Sustainable Fertilizer Supply Chain Transition

Local renewable ammonia production using electrolytic hydrogen is an emerging approach to alleviate emissions attributed to synthetic nitrogen fertilizer production while also insulating against fluctuations in fertilizer prices and mitigating transportation costs and emissions. However, replacing ammonia currently produced using fossil fuels will not be immediate. To this end, we develop a supply chain transition model, which first optimizes the design and hourly operation of new renewable ammonia facilities to minimize production costs and then optimizes the annual installation timing, production scale, and location of these new renewable facilities along with ammonia transportation to meet county resolution demands. The objective is to augment and eventually replace conventional ammonia market imports in an economically competitive manner. We performed a case study for Minnesota's ammonia supply chain and found that a full transition to in-state renewable production by 2032 is optimal. This is incentivized by the U.S. federal government's clean hydrogen production credits. This transition results in 99 % reduction in carbon intensity along with stable supply costs below $475 per metric tonne. New renewable production facilities are an order of magnitude smaller than existing conventional plants. They use both wind and solar resources and operate dynamically to minimize expensive battery and hydrogen storage capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polymer principles behind solubilizing lignin with organic cosolvents for bioenergy

Lignin solubilization is key to a viable biorefinery because its removal leads to facile deconstruction of biomass and because the isolated lignin can serve to derive precursors of novel high-value materials. The mixing of organic solvents with water has been shown to improve biomass fractionation and lignin conversion reactions. However, generally-applicable solubilization strategies are lacking because of the remarkable variability of lignin across plant feedstocks. Here, to obtain a predictive understanding of lignin solvation, we perform molecular dynamics simulations of model lignin polymers in two mixtures of water with polar aprotic solvents: tetrahydrofuran (THF) : water and γ-valerolactone (GVL) : water. The model lignins include H-, G- and S-only homopolymers and lignins with an S : G 1 : 1 ratio. We find that a well-established theory of self-avoiding polymers in a “good” solvent describes accurately the physical conformations of all types of lignin in both solvents. As the degree of methoxy substitution increases in the homopolymers, the distributions of the lignin radius of gyration and the Flory exponent ν , which describes the lignin-solvent interactions, do not change in THF : water, while ν shifts to slightly higher values in GVL : water. We attribute this increase to the interaction between the methyl group of GVL with the lignin methoxy groups. We also find that the reduction in the lignin radius of gyration due to branching is accurately described by the Zimm–Stockmayer theory for both THF : water and GVL : water. The above findings validate the applicability of polymer physics concepts to lignin and suggest that GVL : water may have the most favorable interaction with S-lignin, whereas the interactions of THF : water with lignin are independent of lignin monomeric content.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A rapid embodied carbon assessment tool for priority materials

Embodied carbon limits within building materials are a driving factor in global trade, generating new research and analysis tools in industry. These product assessments, which require utilizing life-cycle assessment (LCA) across broad supply chains, can be expensive, time and data-intensive, and subject to significant variations. Existing methods and tools, such as specific environmental product declarations, typically do not capture these variations and dynamics in supply and manufacturing. Moreover, models and tools must enable stakeholders to assess customized supply chains and future scenarios. In this study, we present the Rapid Embodied Carbon Assessment and Target-setting for Emissions-intensive Materials (REDuCE) tool for building materials. We developed a tool that allows users to select production technologies, transportation mode and distances, concrete carbonation, fuel sources, and regional electricity mixes to supply customization for cement and concrete produced and consumed in California. We generated and integrated a material demand model using residential building stock projections. We provide the user with a wide range of mitigation alternatives through low-carbon production pathways, material use efficiency, transportation modes, and projected electricity grid mixes. In a case study application through four mitigation scenarios, we find emission savings up to 80% by maximizing user mitigation alternatives, primarily driven by reductions in material use intensities. This work represents a foundation for expanding LCA and embodied carbon tools to better enable stakeholders to rapidly and accurately assess customized supply chains while meeting trade requirements.

building materials↗

"Designing, simulating, and performing the 100-AV field test for the CIRCLES consortium: Methodology and Implementation of the Largest mobile traffic control experiment to date"

Previous controlled experiments on single-lane ring roads have shown that a single partially autonomous vehicle (AV) can effectively mitigate traffic waves. This naturally prompts the question of how these findings can be generalized to field operational, high-density traffic conditions. To address this question, the Congestion Impacts Reduction via CAV-in-the-loop Lagrangian Energy Smoothing (CIRCLES) Consortium conducted MegaVanderTest (MVT), a live traffic control experiment involving 100 vehicles near Nashville, TN, USA. This article is a tutorial for developing analytical and simulation-based tools essential for designing and executing a live traffic control experiment like the MVT. It presents an overview of the proposed roadmap and various procedures used in designing, monitoring, and conducting the MVT, which is the largest mobile traffic control experiment at the time. The design process is aimed at evaluating the impact of the CIRCLES AVs on surrounding traffic. The article discusses the agent-based traffic simulation framework created for this evaluation. A novel methodological framework is introduced to calibrate this microsimulation, aiming to accurately capture traffic dynamics and assess the impact of adding 100 vehicles to existing traffic. The calibration model's effectiveness is verified using data from a six-mile section of Nashville's I-24 highway. The results indicate that the proposed model establishes an effective feedback loop between the optimizer and the simulator, thereby calibrating flow and speed with different spatiotemporal characteristics to minimize the error between simulated and real-world data. Finally, We simulate AVs in multiple scenarios to assess their effect on traffic congestion. This evaluation validates the AV routes, thereby contributing to the execution of a safe and successful live traffic control experiment via AVs.

Ameli, Mostafa↗

Noncentrosymmetric Triangular Magnet CaMnTeO 6 : Strong Quantum Fluctuations and Role of s 0 versus s 2 Electronic States in Competing Exchange Interactions

Abstract Noncentrosymmetric triangular magnets offer a unique platform for realizing strong quantum fluctuations. However, designing these quantum materials remains an open challenge attributable to a knowledge gap in the tunability of competing exchange interactions at the atomic level. Here, a new noncentrosymmetric triangular S = 3/2 magnet CaMnTeO 6 is created based on careful chemical and physical considerations. The model material displays competing magnetic interactions and features nonlinear optical responses with the capability of generating coherent photons. The incommensurate magnetic ground state of CaMnTeO 6 with an unusually large spin rotation angle of 127°(1) indicates that the anisotropic interlayer exchange is strong and competing with the isotropic interlayer Heisenberg interaction. The moment of 1.39(1) µB , extracted from low‐temperature heat capacity and neutron diffraction measurements, is only 46% of the expected value of the static moment 3 µB . This reduction indicates the presence of strong quantum fluctuations in the half‐integer spin S = 3/2 CaMnTeO 6 magnet, which is rare. By comparing the spin‐polarized band structure, chemical bonding, and physical properties of AMnTeO 6 ( A = Ca, Sr, Pb), how quantum‐chemical interpretation can illuminate insights into the fundamentals of magnetic exchange interactions, providing a powerful tool for modulating spin dynamics with atomically precise control is demonstrated.

36 MATERIALS SCIENCE↗

Modelling and analysis of salt-convection effect on oxide reduction process for uranium oxides using smoothed particle hydrodynamics

Oxygen ion transport via molten salt continuum is a key kinetic feature of electrochemical reduction process for uranium oxides in molten salt. During the electrolysis, the evolution of oxygen ion concentration field in molten salt continuum appears from porous metallic layer to open bulk salt. Despite the important role affecting process throughput, its effect coupled with controlled salt movement was not adequately investigated. To capture the reactive oxygen ion transport phenomena, a multidimensional reduction model based on smoothed particle hydrodynamics (SPH) was combined with the advection-diffusion model. A pelletized oxide fuel bed capturing interconnected salt continuum is considered as a reactant material geometry, which allows forced convective delivery of molten salt through the reactant materials with a relatively low pressure head. Here, in this study, reduction time of a pelletized-fuel bed was quantified with respect to key design parameters of the reactant material configuration. The simulation results demonstrated the advantage of salt convection which effectively removes oxygen ion products from the reactant materials and consequently accelerates the reduction process. Moreover, the reduction-time correlation of a single pelletized fuel was developed as a function of relevant dimensionless parameters and was extended to the application for a pelletized-fuel bed. Finally, several insights and implications in the use of a pelletized oxide fuel were presented. The developed SPH-based framework enabled multi-physics analysis of governing dynamics in appreciating and screening candidate reactant material configurations in a cost-effective manner.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Perturbative model for the saturation of energetic-particle-driven modes limited by self-generated zonal modes

We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.

Alfvén eigenmodes↗

2023 Roadmap on molecular modelling of electrochemical energy materials

New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.

25 ENERGY STORAGE↗

Planning Optimization for Inductively Charged On-demand Automated Electric Shuttles Project at Greenville, South Carolina

Wireless charging technology presents an ideal fit for autonomous electric vehicles for realizing a fully automated system (vehicle and charger). This paper presents a planning optimization analysis for a pilot project of in-route wireless charging infrastructure serving fixed-route on-demand shared automated electric shuttles (SAESs) at Greenville, South Carolina, USA. A single-objective non-linear mixed integer system planning optimization problem is formulated. A comprehensive cost model representing the overall inductively charged SAESs system is developed, considering road construction, power electronics and materials, traction battery, and installation costs. The optimization problem is solved to determine the best combination of the system key design parameters (number and allocations of wireless chargers, charging power level, track length and on-board battery capacity) that show the most cost-effective solution and allow the SAESs to achieve charge-sustaining operation. The planning platform incorporates representative simulated traffic data (driving speed and routes) for four SAESs in the Greenville using the Simulation of Urban Mobility tool. These data are fed to a vehicle powertrain model and a wireless charger power model to predict the battery power, energy and state-of-charge profiles, which are provided to the search algorithm to assess the design objectives under specific constraints. Finally, the results indicate that implementing high-power (100 kW) wireless charger at a few designated stops for fixed-route SAESs with the proper track length allows the vehicles to realize charge-sustaining operation, infinite range, and zero recharge downtime, with a significant reduction in the on-board battery (36%) and road coverage (69%), at minimum cost.

33 ADVANCED PROPULSION SYSTEMS↗

Accounting for the effect of dislocation climb-mediated flow on the anisotropy and texture evolution of Mg alloy, AZ31B

The current study explores the role that dislocation climb has in mediating plasticity a Mg alloy at moderately elevated temperatures. Interrupted tensile tests were performed on samples of Mg alloy, AZ31B, sheet in the O temper condition over a range of strain rates (10 -5 to 10 -1 s -1 ) and temperatures (20–350°C) along the rolling and transverse directions. Experimental measurements of the resulting strain anisotropy and texture evolution were used as constraints during a parametric study employing a new crystal plasticity model (VPSC-CLIMB), which explicitly accounts for the kinematics of dislocation climb. The results reveal that the climb of basal (a) dislocations is not only important for dislocation recovery, but also demonstrate that climb accommodates a significant fraction of the strain in conditions where a power-law creep-type constitutive response prevails. This work does not discredit the notion that non-basal slip of (a) and (c+a) dislocations is important over a wide range of temperatures and strain rates. However, it demonstrates that the activation of dislocation climb as the mechanistic change within the power law regime provides an explanation for a wide range of observations, including the simultaneous reduction in strain anisotropy, slowed texture evolution, and rapid increase in strain rate sensitivity. Finally, it is hypothesized that these conclusions may even apply to cases in which grain boundary sliding and/or dynamic recrystallization are observed.

36 MATERIALS SCIENCE↗

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

Floating solar panels on reservoirs impact phytoplankton populations: A modelling experiment

Floating solar photovoltaic (FPV) deployments are increasing globally as the switch to renewable energy intensifies, representing a considerable water surface transformation. FPV installations can potentially impact aquatic ecosystem function, either positively or negatively. However, these impacts are poorly resolved given the challenges of collecting empirical data for field or modelling experiments. In particular, there is limited evidence on the response of phytoplankton to changes in water body thermal dynamics and light climate with FPV. Given the importance of understanding phytoplankton biomass and species composition for managing ecosystem services, we use an uncertainty estimation approach to simulate the effect of FPV coverage and array siting location on a UK reservoir. FPV coverage was modified in 10% increments from a baseline with 0% coverage to 100% coverage for three different FPV array siting locations based on reservoir circulation patterns. Results showed that FPV coverage significantly impacted thermal properties, resulting in highly variable impacts on phytoplankton biomass and species composition. The impacts on phytoplankton were often dependent on array siting location as well as surface coverage. Changes to phytoplankton species composition were offset by the decrease in phytoplankton biomass associated with increasing FPV coverage. We identified that similar phytoplankton biomass reductions could be achieved with less FPV coverage by deploying the FPV array on the water body's faster-flowing area than the central or slower flowing areas. The difference in response dependent on siting location could be used to tailor phytoplankton management in water bodies. Simulation of water body-FPV interactions efficiently using an uncertainty approach is an essential tool to rapidly develop understanding and ultimately inform FPV developers and water body managers looking to minimise negative impacts and maximise co-benefits.

14 SOLAR ENERGY↗

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗