Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity Coefficients”

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 217 records · Page 12

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗

Sensitivity study on the fin effect of thermocouple mounted on the heated surface under film boiling condition

In the reactivity-initiated accident (RIA) scenario, right after the power burst with the occurrence of departure from nucleate boiling (DNB), the cladding temperature quickly reaches its maximum and is maintained in a steady film boiling regime. Then, under limited heat transfer performance of film boiling, the local cooling effect by the thermocouple wire may not be neglected. A computational model is developed to estimate the local cooling effect of the thermocouple wire mounted on the cladding surface under film boiling condition. The model predicts the temperature profile along the cladding surface as physically expected and is verified against the benchmark problem. Sensitivity studies reveal the significant effects of the thickness of vapor film layer and the diameter of thermocouple on the temperature drop near the thermocouple, for they increase the heat transfer area to the liquid water region. The impact of film boiling heat transfer coefficient, however, is not significant, as the primary heat transfer mechanism near the thermocouple is the heat conduction to the wire. Future study will include the contribution of radiative heat transfer between the cladding surface and the vapor film layer, as well as a transient power history which represents an actual RIA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Influence of uranyl complexation on the reaction kinetics of the dodecane radical cation with used nuclear fuel extraction ligands (TBP, DEHBA, and DEHiBA)

Specialized extractant ligands – such as tri-butyl phosphate (TBP), N,N-di-(2-ethylhexyl)butyramide (DEHBA), and N,N-di-2-ethylhexylisobutryamide (DEHiBA) – have been developed for the recovery of uranium from used nuclear fuel by reprocessing solvent extraction technologies. These ligands must function in the presence of an intense multi-component radiation field, and thus it is critical that their radiolytic behaviour be thoroughly evaluated. This is especially true for their metal complexes, where there is negligible information on the influence of complexation on radiolytic reactivity, despite the prevalence of metal complexes in used nuclear fuel reprocessing solvent systems. Here we present a kinetic investigation into the effect of uranyl (UO 2 2+ ) complexation on the reaction kinetics of the dodecane radical cation (RH˙ + ) with TBP, DEHBA, and DEHiBA. Complexation had negligible effect on the reaction of RH˙ + with TBP, for which a second-order rate coefficient (k) of (1.3 ± 0.1) × 10 10 M -1 s -1 was measured. For DEHBA and DEHiBA, UO 2 2+ complexation afforded an increase in their respective rate coefficients: k(RH˙ + + [UO 2 (NO 3 ) 2 (DEHBA) 2 ]) = (2.5 ± 0.1) × 10 10 M -1 s -1 and k(RH˙ + + [UO 2 (NO 3 ) 2 (DEHiBA) 2 ]) = (1.6 ± 0.1) × 10 10 M -1 s -1 . This enhancement with complexation is indicative of an alternative RH˙ + reaction pathway, which is more readily accessible for [UO 2 (NO 3 ) 2 (DEHBA) 2 ] as it exhibited a much larger kinetic enhancement than [UO 2 (NO 3 ) 2 (DEHiBA) 2 ], 2.6× vs. 1.4×, respectively. Complementary quantum mechanical calculations suggests that the difference in reaction kinetic enhancement between TBP and DEHBA/DEHiBA is attributed to a combination of reaction pathway (electron/hole transfer vs. proton transfer) energetics and electron density distribution, wherein attendant nitrate counter anions effectively ‘shield’ TBP from RH˙ + electron transfer processes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Thickness and strain dependence of piezoelectric coefficient in BaTiO 3 thin films

We explore the thickness dependence of the converse piezoelectric coefficient (d 33 ) in epitaxial thin films of BaTiO 3 (BTO) grown on (001) SrTiO 3 substrates. Piezoresponse force microscope was performed using an atomic force microscope equipped with an interferometric displacement sensor allowing direct quantification of electromechanical coupling coefficients in BTO free from unwanted background contributions. We find that 80-nm-thick films exhibit a d 33 of ~20.5pm/V, but as the thickness is reduced, the d 33 reduces to less than 2 pm/V for a 10 nm film. To explain the atomistic origin of the effect, we performed molecular dynamics simulations with a recently developed ab initio-derived reactive force field, constructed using the ReaxFF framework. Simulations predict that under applied electric fields thin films of BaTiO 3 show an increasing thickness, with compressive strain, of the region screening the depolarization-field. This study confirms quantitatively the drop in piezoelectric performance in BTO ultrathin films and again highlights the importance of the screening mechanisms when films approach the ultrathin limits in dictating the functional behaviors.

36 MATERIALS SCIENCE↗

Effect of Solvent on the Local Structure, Dynamics, and Vibrational Density of States in Sn-BEA Zeolite

Lewis acid zeolites are attractive catalysts for epoxidation and biomass valorization, as they are highly active and selective in the liquid phase and can operate at or near ambient conditions. While a rich experimental literature exists on liquid-phase Lewis acid zeolite catalysis, our understanding of the molecular organization and solvent dynamics in the vicinity of Lewis acid sites with differing metal site speciation remains limited. In this work, we investigate the molecular coordination and diffusion of two common solvents (methanol and water) around the closed and open Sn-BEA zeolite active sites using molecular dynamics simulations with a machine-learned interatomic potential trained on ab initio molecular dynamics trajectories. Molecular dynamics simulations reveal that introducing active sites significantly enhances local order in the first and second solvation shells compared to the pure silica case. For methanol, both closed and open active sites are singly coordinated, while more than two water molecules coordinate the open site. In contrast to methanol, we observed that water molecules dissociate, leading to the formation of additional Sn-OH and silanol groups away from the active site. The diffusion coefficients of water and methanol are functions of the solvent population in the pore. Here, our work provides insights into how active site speciation in Lewis acid zeolites affects solvent coordination, diffusion, and vibrational signature. This information is foundational for catalyst design and optimization of liquid-phase catalytic processes in zeolites. It also demonstrates the suitability of machine-learned interatomic potentials for modeling reactive systems, enabling sufficiently long trajectories for appropriate statistical averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of carbamate formation on the surface tension of epoxy‐amine curing systems

Abstract Epoxy‐amine reactive systems, including their surface tension, have been studied due to their use and importance in composites and adhesives in laboratory and industrial applications. Understanding how the system's surface tension changes during curing could help target applications for tailored composites. Using the Wilhelmy plate method, surface tensions for epoxy‐curing agent systems were investigated under different environments. While the epoxy exhibited an expected negative temperature coefficient of −0.096 mN/m/°C, the curing agent exhibited a positive temperature coefficient of 0.073 mN/m/°C, which is due to the curing agent reacting with atmospheric CO 2 to form carbamate through a zwitterion pathway and functional groups accumulating at the surface of the reacting systems. Carbamate formation results in a surface tension below that of the individual constituents in the epoxy‐amine curing agent system. Fourier transform infrared (FTIR) spectroscopy was used to verify in situ formation of carbamate, and by reducing the system's temperature, it was demonstrated that the system surface tension profiles displayed micellar surfactant‐like characteristics. FTIR was used to explore the difference in composition between the surface and the bulk, in particular carbamate formation at the surface, resulting in a lowering of the surface tension of reacting systems.

Shenk, Timothy M.↗

Accelerated Zymonic Acid Formation from Pyruvic Acid at the Interface of Aqueous Nanodroplets

To explore the role of the liquid interface in mediating reactivity in small compartments, the formation kinetics of zymonic acid (ZA) is measured in submicron aerosols (average radius = 240 nm) using mass spectrometry. The formation of ZA, from a condensation reaction of two pyruvic acid (PA) molecules, proceeds over days in bulk solutions, while in submicron aerosols, it occurs in minutes. The experimental results are replicated in a kinetic model using an apparent interfacial reaction rate coefficient of krxn = (0.9 ± 0.2) × 10-3 M-1 s-1. The simulation reveals that surface activity of PA coupled with an enhanced interfacial reaction rate drives accelerated ZA formation in aerosols. Experimental and simulated results provide compelling evidence that the condensation reaction of PA occurs exclusively at the aerosol interface with a reaction rate coefficient that is enhanced by 4 orders of magnitude (∼104) relative to what is estimated for macroscale solutions.

Kim, Pyeongeun↗

Reversible Reactions, Mesh Size, and Segmental Dynamics Control Penetrant Diffusion in Ethylene Vitrimers

The diffusion of two aromatic dyes with nearly identical sizes was measured in ethylene vitrimers with precise linker lengths and borate ester cross-links using fluorescence recovery after photobleaching (FRAP). One dye possessed a reactive hydroxyl group, while the second was inert. The reaction of the hydroxyl group with the network is slow relative to the hopping times of the dye, resulting in a large slowdown by a factor of 50 for a reactive probe molecule. A kinetic model was fit to the fluorescence intensity data to determine rate constants for the reversible reaction of the dye from the network, which confirms the role of slow reaction kinetics. A second network cross-linker was also investigated with a substituted boronic ester showing ∼10,000 times faster exchange kinetics. In this system, the two dyes show the same diffusion coefficient, as the reaction is no longer the rate-limiting step. The role of dense meshes on small and large dyes is also discussed in the context of the existing theories. Finally, these results highlight the potential of dynamic networks to control penetrant transport through synergistic effects of the mesh size, dynamic bond kinetics, and penetrant–network interactions.

confinement↗

Uncertainty quantification and sensitivity analysis for SPERT III E-core reactivity measurement benchmarking

The Special Power Excursion Reactor Test (SPERT) III E-core experiment is important because it provides critical data on reactor behavior under significant reactivity insertions, which is essential for validating computer simulations and ensuring the safety of modern light water reactors. Its design similarities to contemporary reactors make it a valuable resource for understanding and mitigating extreme hazards in nuclear operations. The current study details the application of formal parametric uncertainty quantification and sensitivity analysis to a model of the SPERT-III E-core model for zero power reactivity benchmarking. Additionally, the reactivity impact from various modeling assumptions is quantified. Overall, an conservative estimate for an uncertainty in k$_{\text{eff}}$ of $\pm$1257 was observed. A less conservative, more realistic, uncertainty estimate of $\pm$1096 pcm can be justified by the potential for various parametric uncertainties to become negligible when sampled independently across the ~1400 pins in the core. The experimental results fall within both of these uncertainty bounds. Standardized regression coefficient as well as Sobol indices are used to identify the guide tube thicknesses as the primary contributors to the uncertainty in k$_{\text{eff}}$. Overall, this study provides information on how the uncertainties in input parameters and modeling methods impact simulated k$_{\text{eff}}$ values and can be used to aid model building efforts for future code validation with the SPERT-III E-core experiment.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Kinetics of Radiation-induced Cr( II ) and Cr( III ) Redox Chemistry in Molten LiCl-KCl Eutectic

Chromium (Cr) is a frequent constituent of the metal alloys proposed for molten salt nuclear reactor (MSR) applications, and is typically the least noble metal ion present. Consequently, chromium is preferentially corroded into molten salt solutions. Here, the redox poise and redox cycling of chromium ions in the salt can greatly influence its corrosivity towards structural alloys, ultimately impacting the longevity of MSR systems. Radiation-induced chemistry is expected to play a significant role in determining the chromium oxidation state distribution during MSR operations. In the present research, electron pulse radiolysis techniques were employed to characterize the reactivity of Cr(II) and Cr(III) ions with primary radiolysis products in molten lithium chloride–potassium chloride (LiCl–KCl) eutectic over a temperature range of 400–600 °C. Both chromium oxidation states were found to rapidly react with the primary products of molten chloride salt radiolysis, i.e., the solvated electron (e S - ) and the dichlorine radical anion (Cl 2 ˙ - ). For reactions with the e S - , second-order rate coefficients (k) of k = (4.1 ± 0.2) and (6.1 ± 0.3) × 10 10 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively, were determined. Temperature-dependent measurements allowed for the derivation of activation parameters for electron capture by Cr(II) and Cr(III). Both chromium ions also react with Cl 2 ˙ - , k = (7.2 ± 0.3) and (1.4 ± 0.1) × 10 9 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Model predictive control of mixing controlled compression ignition operation for low reactivity fuels

Using gasoline or other low reactivity fuels with a pilot injection or port fuel injection in a compression ignition engine has shown great potential in reducing NOx emissions while keeping high thermal efficiency compared to diesel. However, excessive combustion noise is caused by a high maximum pressure rise rate in the cylinder due to the higher fractions of premixed charge of the low-reactivity fuel. This noise can result in structural damage to engine components and as such, combustion noise limits the range of the operating parameters and makes the control of such engines challenging. In this study, a simulation environment was built up in MATLAB/Simulink leveraging a physics-based zero-dimension combustion model to capture the in-cylinder pressure time traces as well as metrics relevant to thermal efficiency and combustion noise. Here, in order to also facilitate the control of emissions, machine learning models were investigated to capture NOx emissions. A kernel-based extreme learning machine (K-ELM) performed best and had a coefficient of correlation (R-squared) of 0.998. The combustion and NOx emission models are valid for not only conventional gasoline fuel but also oxygenated alternative fuel blends at three different pilot injection strategies. In order to track key combustion metrics while keeping noise and emissions within constraints, a model predictive control (MPC) was applied for a compression ignition engine operating with a range of potential fuels and fuel injection strategies. The MPC is validated under different scenarios, including a load step change, fuel type change, and injection strategy change, with proportional–integral (PI) control as the baseline. The simulation results show that MPC reduces about 26% of ringing intensity in the transient process and 17% at the steady state for E30. Generally, MPC can optimize the overall performance through modifying the main injection timing, pilot fuel mass, and exhaust gas recirculation (EGR) fraction.

42 ENGINEERING↗

Horizontal Split Table Conceptual Design for Advanced Reactor Validation

Oak Ridge National Laboratory and Lawrence Livermore National Laboratory are collaborating to develop a conceptual design for a horizontal split table for use in performing critical experiments. The goal of this design effort is to provide nuclear data testing and validation capabilities for advanced reactors such as pebble-bed high-temperature gas-cooled reactors, molten salt reactors, and heat pipe microreactor, but it could also be used for the current generation of reactors. The first concept being explored for the horizontal split table, a pebble-bed design based on the HTR-10 reactor, is described in this paper. A critical configuration corresponding to a footprint of about 4.5 m 2 was determined with SCALE/KENO-VI to fit the planned dimensions of the horizontal split table. The similarity of the pebble-bed design and the HTR-10 reactor application was assessed using SCALE/TSUNAMI, and a similarity coefficient $c_k$ of 0.9982 was obtained, proving that the concept will be useful for nuclear data validation and assimilation of pebble-bed type advanced reactors. In the proposed design, the materials with the highest $k_{eff}$ sensitivity are graphite and uranium, demonstrating that particular care must be given to carbon-related cross-section data. The effect of mechanical uncertainties between the fixed and moving tables was also assessed by calculating the reactivity change caused by vertical and horizontal gaps, as well as angular and torsion offsets between the two sides of the horizontal split table concept. The highest relative changes on the concept’s reactivity were caused by angular perturbations. The same analysis process is currently being used to create a molten salt advanced reactor type horizontal split table concept based on the Molten Salt Reactor Experiment (MSRE).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ferroelectricity of wurtzite Al 1−x Hf x N heterovalent alloys

Thin films of aluminum hafnium nitride (Al1−xHfxN) were synthesized via reactive magnetron sputtering for Hf contents up to x = 0.13. X-ray diffraction showed a single c-axis oriented wurtzite phase for all films. Hard x-ray photoelectron spectroscopy demonstrated homogeneous Al:Hf distribution through the thin films and confirmed their insulating character. A collection of complementary tests showed unambiguous polarization inversion, and thus ferroelectricity in multiple samples. Current density vs electric field hysteresis measurements showed distinct ferroelectric switching current peaks, the piezoelectric coefficient d 33,f,meas measured using a double beam laser interferometer (DBLI) showed a reversal in sign with similar magnitude, and anisotropic wet etching confirmed field-induced polarization inversion. This demonstrates the possibility of using tetravalent–and not just trivalent–alloying elements to enable ferroelectricity in AlN-based thin films, highlighting the compositional flexibility of ferroelectricity in wurtzites and greatly expanding the chemistries that can be considered for future devices.

36 MATERIALS SCIENCE↗