Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter 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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Cummins R-SOFC System Development

The overall purpose of this project was to reduce the Reversible-Solid Oxide Fuel Cell (R-SOFC) system cost by developing two technologies, an improved cell design and the incorporation of an ejector in the fuel recycle loop instead of a blower. A Simulink model of the baseline SOFC system was developed and calibrated with experimental test data. The R-SOFC system model was built by integrating GT Suite developed models of the steam generation components into the baseline Simulink SOFC system model. The ability to run the stack in SOEC operating mode was also added to the model. The system model was used to explore the ability of the R-SOFC system to meet operational constraints on Steam/Carbon ratio and H2 concentration on the fuel side electrode. A CFD ejector model was developed and used to explore a range of ejector design parameters, leading to the final ejector design that was prototyped for testing. A prototype steam ejector was first tested in a laboratory environment using room temperature air. The steam ejector was subsequently tested using the full hot recycle loop with all relevant heat exchangers and steam generation components. The test conditions utilized temperatures, pressures, and flow rates expected in an R-SOFC application. Throughout the experimental testing work, ejector performance test data was used to improve and then validate the CFD ejector model. A CFD cell model was developed and used to optimize thermal gradients, voltage, and cost of a new cell substrate design. A CFD comparison of co-flow and cross-flow cell designs informed the decision to use a co-flow design for the new cell substrate. Multiple rounds of CFD simulation were used to improve the cell design to minimize the variation in air and fuel distribution across cell channels and to minimize the variation in air and fuel distribution across different cells in the stack. A few prototypes of the new cell substrate design were produced and validated in a laboratory environment by thermally spraying and verifying that they met established manufacturing specifications. The cell manufacturing process was adjusted in order to bring these metrics within acceptable tolerances. Cummins’ internal calculations show that the new cell design reduces cost ~50% compared to the baseline cell, while the ejector + superheater/boiler concept reduces cost of the recycle loop by ~40%. The impact of these cost reductions on the cost of producing H2 will depend on the specific system where they are applied. Therefore, a Techno-Economic Analysis was completed using system cost as a variable, and showing how the NREL Current and Future system costs translate into H2 production cost.

Henrichsen, Lars↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2025

Searching for Sterile Neutrinos In the NOvA Experiment and Measurements of Hadronic Energy Response with NOvA Test Beam Data

The NOvA (NuMI Off-Axis electron neutrino Appearance) experiment is a long-baseline neutrino oscillation experiment composed of two functionally identical detectors: a 300-ton Near Detector and a 14-kton Far Detector, separated by 809 km and placed 14 mrad off the axis of the NuMI neutrino beam created at Fermilab. This configuration enables NOvA’s rich neutrino physics program, which includes measuring neutrino mixing parameters, determining the neutrino mass hierarchy, probing CP violation in the leptonic sector, studying neutrino cross sections, searching for sterile neutrinos, and more.\\ In this thesis I will present two independent analyses developed within the NOvA experiment. The first is a search for light sterile neutrino oscillations in the 3+1 framework using $\nu_\mu$-CC and NC samples from the NOvA Near and Far Detectors. Light sterile neutrinos are hypothetical neutral leptons that do not participate in the weak interactions but can mix with the three known active neutrinos: $\nu_e$, $\nu_\mu$, and $\nu_\tau$. Anomalous results, such as $\nu_e$/$\bar{\nu}e$ appearance in a $\nu_\mu$ ($\bar{\nu}_\mu$) beam observed by the MiniBooNE and LSND experiments, can be explained by the existence of sterile neutrinos. Thus, proving the existence of this type of neutrino is essential. The latest results of these searches in the NOvA experiment, along with contributions to this effort, will be presented.\\ The second is the first measurement of the pion energy response using data from the NOvA Test Beam experiment. The NOvA Test Beam experiment, deployed at Fermilab, uses a scaled-down 30-ton NOvA detector to analyze tagged beamline particles. The beamline can select and identify electrons, muons, pions, kaons, and protons with momenta ranging from 0.4 to 1.4 GeV/c. Pions are an important component of the hadronic system in neutrino interactions, and understanding how the detector responds to these particles is crucial for the validation of the simulation and the reduction of detector calibration uncertainties, which remain one of the largest systematic uncertainties in NOvA analyses.

Dueñas Tonguino, David Francisco [Cincinnati U.] (↗
DOE OSTI2025

Deciphering Reaction Mechanisms of Molecular Proton Reduction Catalysts with Cyclic Voltammetry: Kinetic vs Thermodynamic Control

The kinetics and thermodynamics of elementary reaction steps involved in the catalytic reduction of protons to hydrogen define the reaction landscape for catalysis. The mechanisms can differ in the order of the elementary proton transfer, electron transfer, and bond-forming steps and can be further differentiated by the sites at which protons and electrons localize. Access to fully elucidated mechanistic, kinetic, and thermochemical details of molecular catalysts is crucial to facilitate the development of new catalysts that operate with optimal efficiency, selectivity, and durability. The mechanism by which a catalyst operates, as well as the kinetics and thermodynamics associated with the individual steps, can often be accessed through electroanalytical studies. Here, this Account details the application of cyclic voltammetry to interrogate reaction mechanisms and quantify the kinetics and thermodynamics of elementary reaction steps for a series of molecular catalysts that mediate electrochemical proton reduction. I distinguish the limiting scenarios wherein a catalyst operates under kinetic control vs thermodynamic control, with a focus on detecting how cyclic voltammetry features shift with proton source strength and concentration, as well as scan rate. For systems that operate under kinetic control, catalytic currents are observed at, or slightly positive toward, the formal potential for the redox process that triggers catalysis. Under thermodynamic control, catalytic responses shift as a function of the proton source pKa and effective pH of the solution. After drawing this distinction, we introduce the appropriate voltammetry experiments and accompanying analytical expressions for extracting key metrics from the data. To illustrate analytical strategies to quantify elementary reaction steps of catalysts operating under kinetic control, I describe our studies of proton reduction catalysts Co(dmgBF 2 ) 2 (CH 3 CN) 2 (dmgBF 2 = difluoroboryl-dimethylglyoxime) and [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ (P 2 Ph N 2 Ph = 1,5-phenyl-3,7-phenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, peak shift analysis, foot-of-the-wave analysis, and plateau current analysis are applied to data sets wherein voltammetric response are recorded as a function of catalyst concentration, proton source concentration, proton source strength, and scan rate to quantify rate constants for elementary proton transfer and bond-forming steps in a catalytic cycle. Further, the case study of [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ illustrates how complementary spectroscopic methods can bolster the mechanistic assignment. Collectively, these two studies showcase how detailed mechanistic studies inform on rate-limiting elementary steps in catalysis and other key processes underpinning catalysis. Second, I present analytical strategies to interrogate catalysts operating under thermodynamic control, centered on the case study of [Ni II (P 2 Ph N 2 Bn ) 2 ] 2+ (P 2 Ph N 2 Bn = 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, the application of nonaqueous Pourbaix theory to extract thermodynamic information is introduced, and the construction of a coupled Pourbaix diagram is detailed. This study identifies ligand-based protonation as the key process that places catalysis under thermodynamic control and influences the reaction mechanism. Together, the work detailed in this Account showcases the utility of electroanalytical methods to disentangle complex reaction mechanisms and extract key thermochemical and kinetic parameters for elementary steps of catalysis. Through detailed presentation of the key analytical expressions that underpin these analyses, this Account seeks to facilitate the adoption of cyclic voltammetry by the community to fully extract kinetic, thermochemical, and mechanistic information on electrochemical small-molecule activation.

catalysts↗
DOE OSTI2025

Design and full core fuel performance assessment of high burnup cores for 4-loop PWRs

Increasing the fuel discharge burnup of current light water reactors (LWRs) promises reductions in fuel cycle and/or operations costs. By assuming a constant core power density, the economic gain is enabled by better fuel utilization and/or an increased capacity factor. In this effort to investigate greater than 62 MWd/kgU maximum rod average burnup for 110+ kW/l core power density, two core designs have been developed for a standard 17x17, 193 fuel assemblies pressurized water reactor (PWR). The levelized unit cost methodology is employed to evaluate fuel cycle, operation and maintenance, and capital cost impacts and to examine the economic viability of both core design pathways. Core design and optimization are performed using the commercial STUDSVIK code package. Fuel performance analysis is realized in full core configuration via auditing FRAPCON4.1, FAST1.2, and the high-fidelity code BISON. To provide a realistic assessment, the core design process takes into consideration best practices in current PWR core design. It features acceptable performance in terms of various core design constraints on maximum allowable peaking and boron concentration. Gadolinia (Gd2O3) is used as a burnable poison with a maximum of 9 wt% concentration while feeding 89 or 77 fuel assemblies in a 3-batch refueling scheme. Full core fuel performance simulation, which allows for characterization of relevant fuel temperatures, plenum pressures, stresses, and strains, is performed with respect to two bounding burnup levels. Such performance is potentially licensable for the 18-month high burnup core (<68 MWd/kgU peak pin), while it is more challenging for the 24-month high burnup core design pathway (<75 MWd/kgU peak pin). Maximum rod plenum pressure is identified as the most limiting fuel performance parameter. Here, while the scope of the present study focuses on the steady-state plus overpower conditions, the acceptability of the new discharge burnup has to be further assessed by considering uncertainties and impacts under accident scenarios in the future.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

An infeasible-start framework for convex quadratic optimization, with application to constraint-reduced interior-point and other methods

A framework is proposed for solving general convex quadratic programs (CQPs) from an infeasible starting point by invoking an existing feasible-start algorithm tailored for inequality-constrained CQPs. The central tool is an exact penalty function scheme equipped with a penalty-parameter updating rule. The feasible-start algorithm merely has to satisfy certain general requirements, and so is the updating rule. Under mild assumptions, the framework is proved to converge on CQPs with both inequality and equality constraints and, at a negligible additional cost per iteration, produces an infeasibility certificate, together with a feasible point for an (approximately) ℓ 1 -least relaxed feasible problem, when the given problem does not have a feasible solution. The framework is applied to a feasible-start constraint-reduced interior-point algorithm previously proved to be highly performant on problems with many more inequality constraints than variables (“imbalanced”). Numerical comparison with popular codes (OSQP, qpOASES, MOSEK) is reported on both randomly generated problems and support-vector machine classifier training problems. The results show that the former typically outperforms the latter on imbalanced problems. Finally, application of the proposed infeasible-start framework to other feasible-start algorithms is briefly considered, and is tested on a simplex iteration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Particle size and shape effect of Crumbler® rotary shear-milled granular woody biomass on the performance of Acrison® screw feeder: A computational and experimental investigation

Physical experiments and discrete element model (DEM) simulations are conducted to evaluate particle characteristics and operation parameter effects on screw feeding performance for rotary shear-milled Douglas fir. Three performance metrics are used: mass flow rate, shaft driving torque, and specific energy consumption. The impact of particle size, particle size distribution (PSD), shaft rotational speed (rpm), and hopper dimensions on the performance are investigated. All employed performance metrics reveal the superior flowability of the 2-mm particles in contrast to the larger 6-mm counterpart. Remarkably, wider PSD results in poorer flowability than the two mono-sized particles, proving the flowability enhancement achieved by narrower PSD of the woody feedstock. More importantly, DEM simulations unveil PSD-induced degradation in flowability is attributed to mechanical interlocking and particle segregation effects. Furthermore, higher shaft rpm causes higher mass flow rate at the cost of higher specific energy consumption due to viscous dissipation and changes in flow pattern.

09 BIOMASS FUELS↗
DOE OSTI2025

Modifications in the charge trap landscape in Hf 0.5 Zr 0.5 O 2 as a function of oxygen vacancy concentration observed with photoemission electron microscopy

Oxygen vacancies in Hf x Zr (1-x) O 2 (HZO) both contribute to stabilization of the ferroelectric orthorhombic phase and promote leakage pathways that limit the endurance of devices based on the material. Here, for this reason, the defect states of oxygen vacancies were investigated using photoemission electron microscopy (PEEM) and photoluminescence spectroscopy (PL), as their concentration was varied via ex situ laser exposure. Following a controlled oxygen vacancy reduction via visible (2.54 eV) laser dosing of HZO, deep-ultraviolet (DUV, 5.82 eV) PEEM was used to spatially probe the resulting mid-gap defect states and work function. Work function was found to increase monotonically with the laser-induced reduction in oxygen vacancy concentration culminating in a total increase near 70 meV. The change implies a Fermi level shift toward the valence band as the total available electron-filled charge states are reduced with the removal of oxygen vacancies. A reduction in charge states is corroborated by the observed lessening of both photoemission and photoluminescence intensities after laser dosing. The deduced position of the Fermi level is within a band of near-conduction band defect states produced by oxygen vacancies that are linked to endurance limiting leakage currents. Together, these results directly identify the primary role of oxygen vacancies on the defect states in HZO while demonstrating that laser exposure can be used for their modification.

42 ENGINEERING↗
DOE OSTI2024

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Impact of metal diffusion, lattice distortions, native defects, and ambient on dielectric breakdown in Ni–Ga 2 O 3 Schottky diodes

Ga 2 O 3 unipolar devices are of high interest due to their ~8 MV/cm predicted breakdown fields, which have not yet been achieved due to premature device failure. Pre- and post-failure defect analysis of Ni–Ga 2 O 3 Schottky diodes in ultrahigh vacuum (UHV) and air were performed using depth-resolved cathodoluminescence, high angle annular dark field scanning transmission electron microscopy, and energy dispersive x-ray analysis to understand the physical mechanisms that precede premature breakdown. The breakdown voltage in UHV was dramatically reduced by nearly 40% compared with the breakdown in air. This reduction in the breakdown voltage correlated with post-breakdown differences in Ni distribution, indicating that the coordination and bonding of Ni contribute strongly to electrical behavior in Ni–Ga 2 O 3 Schottky diodes. Breakdown studies in UHV revealed that Ni diffuses away more from the metal–semiconductor interface than with air breakdown, where Ni localizes more near the interface, indicative of the preferential formation of a Ni oxide under O-poor conditions. These measurements also identified the formation of divacancy-interstitial complexes and their characteristic luminescence signature ~150 nm from the interface, the densities of which correlated with breakdown fields. These findings show that electric-field-induced degradation occurs via the rearrangement of native point defects, which act as an additional precursor to device failure. Macroscopically, they show the impact of both vacuum conditions and metal reactivity on Ga 2 O 3 device fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of fast and thermal neutron irradiation on Ga-polar and N-polar GaN diodes

Studies of the radiation tolerance and electrical behavior of gallium nitride (GaN) based devices are important for the next generation of high-power and high-voltage electronics that may be subjected to harsh environments such as nuclear reactor and fusion facilities, particle accelerators, and post-denotation environments. Here, in this work, we study the behavior of Ga-polar and N-polar GaN Schottky diodes before and after exposure to fast and thermal + fast neutrons. Temperature-dependent current–voltage (I–V) and circular transmission line method (CTLM) measurements were used to study the electrical characteristics. A strong reduction in reverse leakage current and an increase in differential resistance in forward bias were observed after neutron irradiation. Thermionic emission (TE), Frenkel–Poole (FP) emission, and Fowler–Nordheim (FN) tunneling models were used to explain the forward and reverse I–V characteristics pre- and post-irradiation. The study confirms that Ga-polar and N-polar GaN Schottky diodes exhibit different electrical responses to fast and thermal neutron irradiations. The reverse bias characteristics of N-polar diodes are less affected after the fast neutron irradiation compared to Ga-polar diodes, while in the forward bias region, the electrical behavior after fast and thermal neutron irradiations is similar in Ga-polar and N-polar diodes. The results indicate that the role of orientation should be considered in the design of GaN-based radiation-tolerant electronics.

42 ENGINEERING↗
DOE OSTI2020

Implementation of Ion Exchange Processes for Carbon Dioxide Mineralization Using Industrial Waste Streams

Sequestration of CO 2 within stable mineral carbonates (e.g., CaCO 3 ) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO 2 in an environmentally benign form. However, the pH of aqueous streams equilibrated with gaseous streams containing CO 2 (pH < 4) are typically lower than that which is required for carbonate precipitation (pH > 8). Traditionally, alkalinity is provided by a stoichiometric reagent (e.g., NaOH) which renders these processes environmentally hazardous and economically unfeasible. This work investigates the use of regenerable ion-exchange materials to induce alkalinity in CO 2 -saturated aqueous solutions such that the pH shift required for mineralization occurs without the need for stoichiometric reagents. Na + -H + exchange isotherms (at [H + ] = 10 −8 –10 −1 M) and rates were measured for 13X and 4A zeolites and TP-207 and TP-260 organic exchange resins in batch equilibrium and fixed-bed exchange experiments, respectively. At solutions equilibrated with CO 2 at 1.0 atm (pH = 3.9), H + exchange capacities for the materials were similar (1.7–2.4 mmol H + /g material) and resulted in pH increases from 3.9 to greater than 8.0. Multi-component mixtures using Ca 2+ and Mg 2+ cations (at 10 −3 –10 −1 M) in CO 2 -saturated water were used to probe competitive ion exchange. The presence of divalent cations in solution inhibited H + exchange, reducing capacities to as low as 0.2 mmol H + /g for both resins and zeolites. Dynamic H + exchange capacities in fixed-bed ion exchange columns were similar to equilibrium values for resins (∼1.5 mmol/g) and zeolites (∼0.8 mmol/g) using inlet solutions that were equilibrated with gaseous streams of CO 2 at 1.0 atm. However, exchange kinetics were limited by intraparticle diffusion as indicated by the increased rate parameters with increasing inlet flow rates (20–160 cm 3 min −1 ). Experimental calcite precipitation from mixing the alkaline CO 3 2− -rich water solution obtained from the ion-exchange column with a simulated liquid waste stream solution achieved thermodynamic maximum yields. The results from these studies indicate that ion exchange processes can be used as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO 3 , thereby opening a pathway toward sustainable and economic mineralization processes.

42 ENGINEERING↗
DOE OSTI2022

Effects of a CFD-improved dimple stepped-lip piston on thermal efficiency and emissions in a medium-duty diesel engine

Diesel piston-bowl shape is a key design parameter that affects spray-wall interactions and turbulent flow development, and in turn affects the engine’s thermal efficiency and emissions. It is hypothesized that thermal efficiency can be improved by enhancing squish-region vortices as they are hypothesized to promote fuel-air mixing, leading to faster heat-release rates. However, the strength and longevity of these vortices decrease with advanced injection timings for typical stepped-lip (SL) piston geometries. Dimple stepped-lip (DSL) pistons enhance vortex formation at early injection timings. Previous engine experiments with such a bowl show 1.4% thermal efficiency gains over an SL piston. However, soot was increased dramatically [SAE 2022-01-0400]. In a previous study, a new DSL bowl was designed using non-combusting computational fluid dynamic simulations. This improved DSL bowl is predicted to promote stronger, more rotationally energetic vortices than the baseline DSL piston: it employs shallower, narrower, and steeper-curved dimples that are placed further out into the squish region. In the current experimental study, this improved bowl is tested in a medium-duty diesel engine and compared against the SL piston over an injection timing sweep at low-load and part-load operating conditions. No substantial thermal efficiency gains are achieved at the early injection timing with the improved DSL design, but soot emissions are lowered by 45% relative to the production SL piston, likely due to improved air utilization and soot oxidation. However, these benefits are lost at late injection timings, where the DSL piston renders a lower thermal efficiency than that of the SL piston. Energy balance analyses show higher wall heat transfer with the DSL piston than with the SL piston despite a 1.3% reduction in the piston surface area. Vortex enhancement may not necessarily lead to improved efficiency as more energetic squish-region vortices can lead to higher convective heat transfer losses.

42 ENGINEERING↗
DOE OSTI2023

Quantum simulations in effective model spaces: Hamiltonian-learning variational quantum eigensolver using digital quantum computers and application to the Lipkin-Meshkov-Glick model

Quantum simulations offer the potential to predict the structure and dynamics of nuclear many-body systems that are beyond the capabilities of classical computing. Generally, preparing the ground state of strongly-interacting many-body systems relevant to nuclear physics is however inefficient, even using ideal quantum computers. In addition, currently available NISQ-era quantum devices possess modest numbers of qubits, limiting the size of quantum many-body systems that can be simulated. In this context, a reformulation of the quantum many-body problems using truncated model spaces and Hamiltonians is desirable to make them more amenable to near-term quantum computers. The importance of symmetries in low-energy theories, including effective field theories (EFTs), lattice quantum chromodynamics (QCD), and effective model spaces for nuclear systems, in particular their interplay with the reduction of active Hilbert spaces, is well known. Lesser known is the fact that the non-commutivity of some symmetries and truncations of the model space can be profitably combined with variational calculations to rearrange the entanglement into localized structures and enable more efficient simulations. Here, the goal of the present study is to explore and utilize the non-commutivity of symmetries and model-space truncations of quantum many-body systems important to nuclear physics, particularly in combination with variational algorithms for quantum simulations and effective Hamiltonian learning. We introduce an iterative hybrid classical-quantum algorithm, Hamiltonian learning variational quantum eigensolver (HL-VQE), that simultaneously optimizes an effective Hamiltonian, thereby rearranging entanglement into the effective model space, and the associated ground-state wavefunction. Quantum simulations, using classical computers and IBM's superconducting-qubit quantum computers, are performed to demonstrate the HL-VQE algorithm, in the context of the Lipkin-Meshkov-Glick (LMG) model of interacting fermions, where the Hamiltonian transformation corresponds to an orbital rotation. We use a mapping where the number of qubits scales with the $\log$ of the size of the effective model space, rather than the particle number. HL-VQE is found to provide an exponential improvement in LMG-model calculations of the ground-state energy and wavefunction, compared to naive truncations without Hamiltonian learning, throughout a significant fraction of the Hilbert space. In the context of EFT, this corresponds to counterterms scaling exponentially with the cut-off as opposed to power law. Implementations on IBM's QExperience quantum computers and simulators for 1- and 2-qubit effective model spaces are shown to provide accurate and precise results, reproducing classical predictions. For a range of parameters defining the LMG model, the HL-VQE algorithm is found to have better scaling of quantum resources requirements than previously explored algorithms. In particular, the HL-VQE scales efficiently over a large fraction of the model space, in contrast to VQE alone. This work constitutes a step in the development of entanglement-driven quantum algorithms for descriptions of nuclear many-body systems. This, in part, leverages the potential of noisy intermediate-scale quantum (NISQ) devices. The exponential scaling of counterterms observed in this study suggests the possibility of more general applicability to other non-perturbative EFTs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Updates to a Preliminary Probabilistic Performance Assessment Model for Radiological and Chemical Contamination at the West Valley Site, New York - 20420

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States at the WNYNSC. NFS built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (640 metric tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations halted in 1972 and never restarted, leaving behind radioactive and chemical wastes. The U.S. Department of Energy (DOE) was required to complete certain waste management activities under the West Valley Demonstration Project (WVDP) Act of 1980 including decommissioning of WVDP facilities. As collaborating agencies, NYSERDA and the DOE are tasked with making decisions about decommissioning and risk reduction for the West Valley Site. Neptune and Company, Inc. (Neptune) was contracted to develop a probabilistic performance assessment (PPA) model to assist the agencies in their decision making process for decommissioning the WVDP and WNYNSC. One important tool that is needed in order to inform the decision-making process is a science-based model of the West Valley Site that evaluates potential future consequences for human health and the environment. This forms the core of the spatial domain of the West Valley PPA Model. The PPA Model, developed using the GoldSim system modeling software, is a tool intended to provide support for decision making that evaluates uncertainty, in a manner that is transparent, defensible, and robust. The PPA Model includes contaminant transport and health effects components, and is organized around geographically-grouped contaminated facilities. These include the waste disposal areas licensed by the U.S. Nuclear Regulatory Commission and the State of New York, a waste tank farm for storage of high level radioactive waste resulting from reprocessing operations, and several areas contaminated with radioactive and chemical constituents. The PPA Model evaluates contaminant transport from these sources to points of exposure across the site and into receiving surface waters and sediments. Hypothetical people and wildlife could be exposed to contamination at these locations, and the effects of these exposures are evaluated. Contaminant transport processes to be evaluated in the PPA Model include groundwater and surface water transport, contaminant translocation by plants and animals, diffusion, and erosion. The evaluation of exposures to people in this preliminary model is limited to a resident farmer scenario, and ecological assessment is performed at the level of a screening analysis. The results of these preliminary evaluations inform future model developments. PPA Model results are subjected to sensitivity analysis in order to determine those pathways and parameters that are most significant in influencing the results. This information allows analysts and decision makers to focus on those aspects of Site behavior and processes. With this information, the decision makers can drive informed, defensible decisions regarding decommissioning of the Site. This paper includes an update of the information presented at WM2019 [1]. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Structure-Activity Relationships of Hydrothermally Aged Titania-Supported Vanadium-Tungsten Oxide Catalysts for SCR of NOx Emissions with NH3

Supported V2O5-WO3/TiO2 materials employed as selective catalytic reduction (SCR) catalysts for NOx emission control from power plants. Fresh SCR catalysts usually receive exposure to harsh treatments in industry to accelerate catalyst activation (calcination in air at 650oC) and catalyst aging (hydrothermal aging at 650oC) in a way that represents various points in the catalyst/product lifetime. The present study investigates the catalyst structural and chemical changes occurring during such harsh treatments. Three series of supported V2O5-WO3/TiO2 catalysts were prepared by incipient-wetness impregnation of aqueous ammonium metavanadate and metatungstate precursors. The catalysts were subsequently dried and calcined at 550oC in O2, 650oC in O2, and hydrothermal conditions (10% O2, 8% H2O, 7% CO2 and 75% N2) at 650oC. The resulting catalysts were physically characterized by numerous techniques (in situ Raman, in situ IR, in situ High Field-High Spinning solid-state 51V MAS NMR, in situ EPR, XRD, BET surface area and ICP) and chemically probed with adsorbed ammonia, SCR-TPSR, and the SCR reaction. The surface WOx sites on the TiO2 support behave as a textural promoter that stabilizes the TiO2 (anatase) phase from sintering and transforming to the undesirable crystalline TiO2 (rutile) phase that can lead to formation of a Ti1-xVxO2 (rutile) solid solution with reduced V4+ cations (~7-15%). The surface VOx sites are mostly oligomerized as surface V+5Ox sites (~50-85% oligomers) and the extent of oligomerization tends to increase with surface WOx coverage and calcination temperature. A major difference between the calcined and hydrothermally treated catalysts was the low concentration of surface NH3* species on Lewis acid sites for the hydrothermally treated catalysts, yet the SCR activity was almost comparable for both catalysts. This finding suggests that surface NH4+*, primarily associated with the surface VOx sites, are able to efficiently perform the SCR reaction. Given that multiple catalyst parameters were simultaneously varying during these treatments, it was difficult to correlate the SCR activity with any single catalyst parameter. A correlation, however, was found between the SCR TOF/activity and the sum of the surface NH3* and NH4+* species, which is dominated by the surface NH4+* species.

SCR, NO, NH3, catalyst, V2O5, WO3, TiO2, promotion↗
DOE OSTI2022

Optimizing feed modulation for coupled methane and NO x conversion over Pd-Pt/Mn 0.5 Fe 2.5 O 4 /Al 2 O 3 monolith catalyst

Here the impacts of feed modulation (frequency, amplitude) and catalyst design (composition and architecture) parameters are reported for the conversion of methane and NO x over a dual-layer Pt+Pd/Al 2 O 3 + Mn 0.5 Fe 2.5 O 4 /Al 2 O 3 monolith. CH 4 and NO x conversion data show that the dual-layer catalyst outperforms single-layer samples having the same catalyst loadings, with and without spinel. Close proximity of the PGM and MFO functions in the mixed-layer catalyst lowers the CH 4 conversion at high temperature while separating the PGM and spinel layers with an intermediate Al 2 O 3 layer does not. Methane conversion enhancement is linked to its nonmonotonic dependence on O 2 . The performance gains are tied to a transient activity spike that occurs during the lean-to-rich feed transition when water is present in the feed. The transient spike is attributed to the removal of CO and H 2 products via reactions with stored O 2 in the spinel, eliminating inhibition of methane steam reforming.

03 NATURAL GAS↗