Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric processes”

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,171 records · Page 65

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗

Numerical Investigation of Microgravity Tank Pressure Rise Due to Boiling

The ability to control self-pressurization in cryogenic storage tanks is essential for NASAs long-term space exploration missions. Predictions of the tank pressure rise in Space are needed in order to inform the microgravity design and optimization process. Due to the fact that natural convection is very weak in microgravity, heat leaks into the tank can create superheated regions in the liquid. The superheated regions can instigate microgravity boiling, giving rise to pressure spikes during self-pressurization. In this work, a CFD model is developed to predict the magnitude and duration of the microgravity pressure spikes. The model uses the Schrage equation to calculate the mass transfer, with a different accommodation coefficient for evaporation at the interface, condensation at the interface, and boiling in the bulk liquid. The implicit VOF model was used to account for the moving interface, with bounded second order time discretization. Validation of the models predictions was carried out using microgravity data from the Tank Pressure Control Experiment, which flew aboard the Space Shuttle Mission STS-52. Although this experiment was meant to study pressurization and pressure control, it underwent boiling during several tests. The pressure rise predicted by the CFD model compared well with the experimental data. The ZBOT microgravity experiment is scheduled to fly on February 2016 aboard the ISS. The CFD model was also used to perform simulations for setting parametric limits for the Zero-Boil-Off Tank (ZBOT) Experiments Test Matrix in an attempt to avoid boiling in the majority of the test runs that are aimed to study pressure increase rates during self-pressurization. *Supported in part by NASA ISS Physical Sciences Research Program, NASA HQ, USA

pressure rise↗

Characterizing Oscillatory Bursts in Single-Trial EEG Data

Oscillatory bursts in numerous bands ranging from low (theta) to high frequencies (e.g., gamma) undoubtedly play an important role in cortical dynamics. Largely because of the inadequacy of existing analytic techniques. however, oscillatory bursts and their role in cortical processing remains poorly understood. To study oscillatory bursts effectively one must be able to isolate them and characterize them in the single trial. We describe a series of straightforward analysis techniques that produce useful indices of burst characteristics. First, stimulus-evoked responses are estimated using Differentially Variable Component Analysis (dVCA), and are subtracted from the single-trial. The single-trial characteristics of the evoked responses are stored to identify possible correlations with burst activity. Time-frequency (T-F), or wavelet, analyses are then applied to the single trial residuals. While T-F plots have been used in recent studies to identify and isolate bursts, we go further by fitting each burst in the T-F plot with a two-dimensional Gaussian. This provides a set of burst characteristics, such as, center time. burst duration, center frequency. frequency dispersion. and amplitude, all of which contribute to the accurate characterization of the individual burst. The burst phase can also be estimated. Burst characteristics can be quantified with several standard techniques (e.g.. histogramming and clustering), as well as Bayesian techniques (e.g., blocking) to allow a more parametric description analysis of the characteristics of oscillatory bursts, and the relationships of specific parameters to cortical excitability and stimulus integration.

Knuth, K. H.↗

Sensitivity and Uncertainty Quantification of Transition Scenario Simulations

This report documents the first collective attempt at developing and applying capabilities to quantify uncertainties, assess parametric sensitivities, and optimize multiple parameters and metrics in fuel cycle simulations generated by the SA&I Campaign. To do this, external codes that were designed to perform sensitivity analysis and uncertainty quantification (SA&UQ) needed to be coupled to the SA&I Campaign’s nuclear fuel cycle simulators (NFCS). In FY20, two approaches were pursued: 1) coupling Cyclus to an ORNL-internal code called MOT (Metaheuristic Optimization Tool) and 2) coupling DYMOND to the opensource SA&UQ tool kit Dakota. The primary objective of having these NFCS/SA&UQ coupled capabilities is to better inform DOE-NE and other stakeholders on the results generated from the NFCS. For a given set of fuel cycle strategies, policies, and technology assumptions that make up a fuel cycle scenario, these NFCS have traditionally been used by the SA&I Campaign to provide quantitative answers in terms of year-by-year mass flows, infrastructure requirements, costs, etc. With these newly developed coupled capabilities, the SA&I Campaign can now efficiently simulate hundreds or thousands of these scenarios, sample large ranges of parameters and assumptions, and use the unique features of the SA&UQ tools to process the data. This enables providing answers with known and propagated uncertainties, determining the sensitivity of important metrics to different parameters and assumptions, quantifying how much fuel cycle and technology parameters impact each other, and producing optimized fuel cycle strategies for single and multiple variables. To demonstrate these new capabilities, the Cyclus/MOT was used to model several scenarios ranging from simple fleet retirements to transitions to advanced reactors. Specifically, for a transition scenario from LWRs to SFRs and advanced LWRs, uncertainty quantification, sensitivity analysis, and optimization studies were applied to cases involving single and multiple parameter (input) and single and multiple metric (output) variations. In addition, a similar transition scenario was modeled to demonstrate how to optimize the reprocessing capacity parameter to minimize two performance metrics while taking into account uncertainties from two other parameters. Lastly, a depletion module based on SCALE/ORIGEN was added in Cyclus to simulate the third scenario that was designed to quantify the impact of the modeling assumption that all LWR used nuclear fuel have the same burnup. The newly developed DYMOND/Dakota capability was also applied to a transition scenario from the existing fleet to small modular reactors and fast reactors. This particular scenario involves not only explicit isotopic depletion via ORIGEN-2, but also includes multirecycling and utilizing the criticality search feature to determine the fresh fuel composition of recycled fuel, a feature unique to the DYMOND NFCS. A large database of simulations were run with 4 main parameters that were sampled: start date of reprocessing, reprocessing capacity, energy demand growth rate, and advanced reactor share of the fleet. The 4 main metrics were uranium consumption, enrichment requirements, waste generation, and levelized cost of electricity using data from the Cost Basis Report. The demonstrated SA&UQ results include those that inform on how to choose parameters to avoid “failed” scenarios, Sobol’ indices that inform on the importance of various parameters individually and synergistically, and Analysis of Variance (ANOVA) studies that decompose parameter ranges into groups and informs on whether variations are statistically significant.

Feng, B.↗

Attributes of GRB Pulses: Analysis of BATSE TTE Data

Our extensive program of modeling GRB profiles is aimed at elucidating the physical processes responsible for the burst itself, as well as possible extrinsic phenomena (e.g. time dilation) as discussed in another paper in this Symposium (Norris et al., "GRB PROFILES AS COSMIC PROBES"). We have developed special methods to extract the wealth of short time-scale information contained in the BATSE time-tag event (TTE) data. Our algorithm yields a piecewise-constant representation of the light curve -- using only the raw photon arrival times, and based on Bayesian change-point methods. This representation in effect lets the data determine the bin size and location, and avoids unwanted effects due to arbitrary choices of the bin parameters. We have determined widths, separations, and amplitudes of pulses contained in the bursts, without invoking a specific pulse model. The effect of cosmic time dilation can be easily seen in a direct plot of amplitude vs. time scale for individual pulses, without the need to lump the data into a small number of brightness classes. We are also performing noise equalization on these data (to reduce a well-known bias of pulse width as a function of signal-to-noise ratio), as well as fits of parametric pulse-shape models -- including explicit energy dependence of the pulse parameters. Such refinements are expected to improve the quality and physical significance of these results.

Norris, J. P.↗

Mid-Infrared OPO for High Resolution Measurements of Trace Gases in the Mars Atmosphere

The Martian atmosphere is composed primarily (>95%) of CO2 and N2 gas, with CO, O2, CH4, and inert gases such as argon comprising most of the remainder. It is surprisingly dynamic with various processes driving changes in the distribution of CO2, dust, haze, clouds and water vapor on global scales in the meteorology of Mars atmosphere [I]. The trace gases and isotopic ratios in the atmosphere offer important but subtle clues as to the origins of the planet's atmosphere, hydrology, geology, and potential for biology. In the search for life on Mars, an important process is the ability of bacteria to metabolize inorganic substrates (H2, CO2 and rock) to derive energy and produce methane as a by-product of anaerobic metabolism. Trace gases have been measured in the Mars atmosphere from Earth, Mars orbit, and from the Mars surface. The concentration of water vapor and various carbon-based trace gases are observed in variable concentrations. Within the past decade multiple groups have reported detection of CH4, with concentrations in the 10's of ppb, using spectroscopic observations from Earth [2]. Passive spectrometers in the mid-infrared (MIR) are restricted to the sunlit side of the planet, generally in the mid latitudes, and have limited spectral and spatial resolution. To accurately map the global distribution and to locate areas of possibly higher concentrations of these gases such as plumes or vents requires an instrument with high sensitivity and fine spatial resolution that also has global coverage and can measure during both day and night. Our development goal is a new MIR lidar capable of measuring, on global scales, with sensitivity, resolution and precision needed to characterize the trace gases and isotopic ratios of the Martian atmosphere. An optical parametric oscillator operating in the MIR is well suited for this instrument. The sufficient wavelength tuning range of the OPO can extend the measurements to other organic molecules, CO2, atmospheric water vapor, clouds, temperature, dust, and aerosols, as well as possibly polar-cap properties. Our OPO-approach may allow a new capability for active remote sensing of the outer planets and moons, where the weaker sunlight further limit passive instruments. Here we report on the OPO development effort for this lidar instrument.

Yu, Anthony W.↗

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran↗

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗

Unified architecture for quantum lookup tables

Quantum access to arbitrary classical data encoded in unitary black-box oracles underlies interesting data-intensive quantum algorithms, such as machine learning or electronic structure simulation. The feasibility of these applications depends crucially on gate-efficient implementations of these oracles, which are commonly some reversible versions of the Boolean circuit for a classical lookup table. Here, we present a general parametrized architecture for quantum circuits implementing a lookup table that encompasses all prior work in realizing a continuum of optimal trade-offs between qubits, non-Clifford gates, and error resilience, up to logarithmic factors. Our architecture assumes only local 2D connectivity, yet recovers results, with the appropriate parameters, polylogarithmic error scaling. We also identify regimes, such as simultaneous sublinear scaling, in all parameters. These results enable tailoring implementations of the commonly used lookup table primitive to any given quantum device with constrained resources.

quantum circuits↗

Effective description of ajar systems with a 𝑈⁡(1) symmetry

We introduce the concept of “ajar systems” as an intermediate case between closed and open systems, where the time scale for charge exchange with the environment is parametrically larger than all other characteristic time scales. The Schwinger-Keldysh effective action for such finite-temperature systems exhibits a symmetry group G 1 × G 2 weakly broken to its diagonal subgroup G diag , which we systematically describe using spurion techniques. For G = U(1), we calculate leading-order corrections to correlation functions in both diffusive and spontaneously broken phases. Unlike systems with approximate symmetries where both real and imaginary parts of dispersion relations receive corrections, ajar systems modify only the imaginary part—preserving gapless modes while adding finite damping.

Effective field theory↗

Effects of noise on the overparametrization of quantum neural networks

Overparametrization is one of the most surprising and notorious phenomena in machine learning. Recently, there have been several efforts to study if, and how, quantum neural networks (QNNs) acting in the absence of hardware noise can be overparametrized. In particular, it has been proposed that a QNN can be defined as overparametrized if it has enough parameters to explore all available directions in state space. That is, if the rank of the quantum Fisher information matrix (QFIM) for the QNN's output state is saturated. Here, we explore how the presence of noise affects the overparametrization phenomenon. Our results show that noise can “turn on” previously zero eigenvalues of the QFIM. This enables the parametrized state to explore directions that were otherwise inaccessible, thus potentially turning an overparametrized QNN into an underparametrized one. For small noise levels, the QNN is quasioverparametrized, as large eigenvalues coexists with small ones. Then, we prove that as the magnitude of noise increases all the eigenvalues of the QFIM become exponentially suppressed, indicating that the state becomes insensitive to any change in the parameters. As such, there is a pull-and-tug effect where noise can enable new directions but also suppress the sensitivity to parameter updates. Finally, our results imply that current QNN capacity measures are ill-defined when hardware noise is present. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

A robust plasma-based laser amplifier via stimulated Brillouin scattering

Brillouin amplification in plasma is more resilient to fluctuations in the laser and plasma parameters than Raman amplification, making it an attractive alternative to Raman amplification. In this work, we focus on high plasma densities, n 0 > n cr /4, where stimulated Raman scattering is not possible and laser beam filamentation is the dominant competing process. Through analytic theory and multi-dimensional particle-in-cell simulations, we identify a parameter regime for which Brillouin amplification can be efficient while maintaining filamentation of the probe at a controlled level. We demonstrate pump-to-probe compression ratios of up to 72 and peak amplified probe fluences over 1 kJ cm -2 with ≃50% efficiency. High pulse quality is maintained through control of parasitic filamentation, enabling operation at large beam diameters. Provided the pump and probe pulse diameters can be increased to 1 mm, our results suggest that Brillouin amplification can be used to produce sub-picosecond pulses of petawatt power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Creation of a Rapid High-Fidelity Aerodynamics Module for a Multidisciplinary Design Environment

In the traditional aerospace vehicle design process, each successive design phase is accompanied by an increment in the modeling fidelity of the disciplinary analyses being performed. This trend follows a corresponding shrinking of the design space as more and more design decisions are locked in. The correlated increase in knowledge about the design and decrease in design freedom occurs partly because increases in modeling fidelity are usually accompanied by significant increases in the computational expense of performing the analyses. When running high fidelity analyses, it is not usually feasible to explore a large number of variations, and so design space exploration is reserved for conceptual design, and higher fidelity analyses are run only once a specific point design has been selected to carry forward. The designs produced by this traditional process have been recognized as being limited by the uncertainty that is present early on due to the use of lower fidelity analyses. For example, uncertainty in aerodynamics predictions produces uncertainty in trajectory optimization, which can impact overall vehicle sizing. This effect can become more significant when trajectories are being shaped by active constraints. For example, if an optimal trajectory is running up against a normal load factor constraint, inaccuracies in the aerodynamic coefficient predictions can cause a feasible trajectory to be considered infeasible, or vice versa. For this reason, a trade must always be performed between the desired fidelity and the resources available. Apart from this trade between fidelity and computational expense, it is very desirable to use higher fidelity analyses earlier in the design process. A large body of work has been performed to this end, led by efforts in the area of surrogate modeling. In surrogate modeling, an up-front investment is made by running a high fidelity code over a Design of Experiments (DOE); once completed, the DOE data is used to create a surrogate model, which captures the relationships between input variables and responses into regression equations. Depending on the dimensionality of the problem and the fidelity of the code for which a surrogate model is being created, the initial DOE can itself be computationally prohibitive to run. Cokriging, a modeling approach from the field of geostatistics, provides a desirable compromise between computational expense and fidelity. To do this, cokriging leverages a large body of data generated by a low fidelity analysis, combines it with a smaller set of data from a higher fidelity analysis, and creates a kriging surrogate model with prediction fidelity approaching that of the higher fidelity analysis. When integrated into a multidisciplinary environment, a disciplinary analysis module employing cokriging can raise the analysis fidelity without drastically impacting the expense of design iterations. This is demonstrated through the creation of an aerodynamics analysis module in NASA’s OpenMDAO framework. Aerodynamic analyses including Missile DATCOM, APAS, and USM3D are leveraged to create high fidelity aerodynamics decks for parametric vehicle geometries, which are created in NASA’s Vehicle Sketch Pad (VSP). Several trade studies are performed to examine the achieved level of model fidelity, and the overall impact to vehicle design is quantified.

Aerodynamic Properties↗

Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river

Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗

Chamber Flux and Porewater Concentration of CH4, CO2 and N2O, 2018, Columbia River bank at the Hanford site, WA, USA

This is the corresponding observations data for the paper: Jorge A Villa; Garrett J Smith; Yang Ju; Lupita Renteria; Jordan C Angle; Evan Arntzen; Samuel F Harding; Huiying Ren; Xingyuan Chen; Audrey H Sawyer; Emily B Graham; James C Stegen; Kelly C Wrighton; Gil Bohrer (2020) Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river. Science of the Total Environment.Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow, intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during the low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Compact Augmented Spark Igniters for Liquid Rocket Engines

The objective of this dissertation is to characterize augmented spark impinging (torch) igniters. Torch igniters inject propellants into a small prechamber where a spark igniter imparts the activation energy to initiate combustion. The flame stabilizes within the prechamber and exits into the main combustion chamber where it ignites the main flow directly or by mixing with a coaxial flow to create a secondary flame. The underlying phenomena that promote reliable ignition are multifaceted and required a broad scope of investigations. Initial efforts focused on spark discharges across an annular geometry. Gaseous propellent projects the electrical arc toward the prechamber, where the arc and the thermal energy is imparted to the surrounding fluid then meets with a combustible mixture. Discharges are parametrically examined against pressure, spark-gap size, and exciter types to emulate transient engine startup. The energy imparted to the flow, approximate velocity, and distance traveled by the exhaust plume are determined. Spark quenching (no energy deposition) and its effects on perceived versus actual measurements are shown. Results outline the prechamber target zone where a combustible mixture is necessary for ignition. Next, three-dimensional, time-accurate, and nonreacting computational fluid dynamics simulations assess effects of geometric and mass flow differences on the size, location, and composition of a combustible mixture. Fluid dynamic phenomena are elucidated. Correlations between nondimensional variables are found. Results are coupled with a conservation-law analysis to produce novel objective functions. Lastly, a full-scale, modular torch igniter is tested. Injector sizes, mixture ratios, and injector momentum ratios are systematically varied to map ignition probability and ignition delay. Results show the importance of local mixture ratios, provide injector momentum ratios guidelines, and demonstrate reliable ignition of core mixture ratios typically outside flammability limits. The deliverable of this project is not a torch igniter or derivative thereof. Instead, this dissertation produces novel objective-functions for igniter design. These equations are functions of variables that may be constrained by engine-level requirements, physics, and experimental results from this work. They provide a foundation for future igniter design and will minimize iterative processes from the design cycle for next-generation ignition hardware.

Torch Igniter↗