Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 919 records · Page 51

An Experimental Feasibility Study on Applying Neutron Tomography to Encapsulated Spent Nuclear Fuel - 20050

Visual inspection makes easier to ensure the integrity and safety of spent nuclear fuel (SNF) than any classical techniques. Various classical techniques have been applied but there are no reliable methods to qualitatively and quantitatively verify spent fuel in dry storage. Thus, the present authors have developed the prototype safeguards apparatus for dry storage employing the array of He-4 gas scintillation detectors (S670E, Arktis Radiation Detectors Ltd., Switzerland), newly designed to simultaneously measure thermal and fast neutrons without any moderators. The S670E detector has a cylindrical shape with a diameter of 52 mm and active length of 600 mm (total length: 875 mm). The detector is filled by He-4 gas with an approximate pressure of 180 bar for fast neutron detection, and its inner wall is coated by Li-6 for thermal neutron detection. The scintillation lights generated via Li-6 nuclear reaction and elastic scattering are collected by 24 SiPMs linearly paired at the center of the detector. The detector delivers a TTL (Transistor-Transistor Logic) output for pulse readout and UART (Universal Asynchronous Receiver Transmitter) for device control. In order to assess feasibility of the apparatus, an experimental system has been designed, built, and optimized via computational studies. Cf-252 neutron sources and linearly arrayed detectors, working as a single detector, were occupied for this study due to the difficulties in working with SNF. The laboratory scale cask (diameter: 0.67 m, height: 1.5 m), minimized by a factor of 10 compared to the actual thickness of a commercial TN-32 cask, was also manufactured. The detector array was designed to rotate the lab-scale cask and obtain 36 image profiles at every 10 degrees. All profiles were aligned in single frame image called a sinogram, and the cross-sectional image was then fabricated by the inverse radon transform algorithm. These experiments have been repeated with different configurations and numbers of sources. Some gamma-ray sources were also measured with neutron sources in order to distinguish between neutron and gamma-ray pulses. Basically, a He-4 detector is designed to run on Linux OS so it is difficult to directly apply to Windows-based equipment widely used in S. Korea. Therefore, a new data acquisition board working on Windows OS was designed and built. The board mainly consists of FPGA (Field Programmable Gate Array) and SoC (System on Chip) for TTL pulse readout, sorting measured data, and transferring data to a user interface. In conclusion, the tomographic system at lab-scale has shown considerable potential to detect a partial or gross defect of encapsulated assemblies in dry storage. Next steps of this study will be to 1) repeatedly carry out experiments to demonstrate scientific reliability and validity, and 2) numerically integrate signals with weight factors to enhance the image quality since the suggested system based on passive interrogation method requires longer measurement time. Finally, the system will apply to a commercial dry storage phased out soon in S. Korea. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Transition and Turbulence Modeling for Blunt-Body Wake Flows

Aerobraking has been proposed as an efficient means of decelerating spacecraft for planetary missions. Most current aerobrake designs feature a blunt forebody shielding the payload from the intense heat generated during atmospheric entry. Although this forebody will absorb the largest portion of the heat pulse, accurate prediction of heating in the near wake is of great importance, since large local heating values can occur at points of shear-layer impingement. In order to address the various issues associated with these blunt-body wake flowfields, the Advisory Group for Aerospace Research and Development (AGARD) formed Working Group 18 in 1992. One of the objectives of this activity was to examine real-gas effects in high-speed flow fields around a 70 deg. blunted cone. To date, many researchers have conducted experiments using this geometry in various facilities, such as the Large Energy National Shock (LENS) tunnel at Cubric/Calspan and the HEG shock tunnel at DLR-Goettingen. Several computational studies have also been conducted in concert with these tests. Many of the experimental results have indicated the possible presence of a transitional shear layer through a large increase in heat transfer downstream of the reattachment point. The presence of transition could in fact lead to much higher peak heating than if the separated flow is entirely laminar or turbulent. In the shock-tunnel tests, however, it is difficult to separate such viscous-flow phenomena from real-gas effects. In order to help make this distinction, Horvath et al. recently conducted a set of experiments in the NASA Langley 20-Inch Mach 6 Tunnel, and compared the results to laminar Navier-Stokes calculations. They found heat-transfer distributions similar to those obtained in the high-enthalpy facilities, with the measured peak heating along the sting support markedly greater than that predicted by the laminar computations. These trends point to the need to find transitional and turbulent computational solutions for these flowfields.

Nance, Robert P.↗

Effect of Apex Flap Deflection on Vertical Tail Buffeting

A computational study of the effect of vortex breakdown location on vertical tail buffeting is conducted. The position of the breakdown is modified by employing an apex flap deflected by an experimentally determined optimal angle. The delayed breakdown flow and buffeting response is then compared to the nominal undeflected case. This multidisciplinary problem is solved sequentially for the fluid flow, the elastic tail deformations and the grid displacements. The fluid flow is simulated by time accurately solving the unsteady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, upwind, flux-difference splitting finite volume scheme. The elastic vibrations of the tails are modeled by uncoupled bending and torsion beam equations. These equations are solved accurately in time using the Galerkin method and a five-stage Runge-Kutta-Verner scheme. The grid for the fluid dynamics calculations is continuously deformed using interpolation functions to disperse the displacements smoothly throughout the computational domain. An angle-of-attack of 35 deg.is chosen such that the wing primary-vortex cores experience vortex breakdown and the resulting turbulent wake flow impinges on tile vertical tails. The dimensions and material properties of the vertical tails are chosen such that the deflections are large enough to insure interaction with the flow, and the natural frequencies are high enough to facilitate a practical computational solution. Results are presented for a baseline uncontrolled buffeting case and a delayed breakdown case in which the apex flap has been deflected 15 deg. The flap was found to be very effective in delaying the breakdown, increasing the location from 50%c to 94%c, which resulted in a 6% increase in lift coefficient and pitching moment. However, the integrated buffet loads and tip responses were roughly equivalent for the two cases.

Massey, Steven J.↗

Grand canonically optimized grain boundary phases in hexagonal close-packed titanium

Grain boundaries (GBs) profoundly influence the properties and performance of materials, emphasizing the importance of understanding the GB structure and phase behavior. As recent computational studies have demonstrated the existence of multiple GB phases associated with varying the atomic density at the interface, we introduce a validated, open-source GRand canonical Interface Predictor (GRIP) tool that automates high-throughput, grand canonical optimization of GB structures. While previous studies of GB phases have almost exclusively focused on cubic systems, we demonstrate the utility of GRIP in an application to hexagonal close-packed titanium. We perform a systematic high-throughput exploration of tilt GBs in titanium and discover previously unreported structures and phase transitions. In low-angle boundaries, we demonstrate a coupling between point defect absorption and the change in the GB dislocation network topology due to GB phase transformations, which has important implications for the accommodation of radiation-induced defects.

36 MATERIALS SCIENCE↗

Mechanistic Insights into Molecular Copper Hydride Catalysis: the Kinetic Stability of CuH Monomers toward Aggregation is a Critical Parameter for Catalyst Performance

The activity of molecular copper hydride (CuH) complexes towards the selective insertion of unsaturated hydrocarbons under mild conditions has contributed significantly to versatile methodologies for upgrading these feedstocks. However, these catalysts are particularly susceptible to deleterious aggregation, leading to the depletion of active CuH species. Little is known about the mechanisms of CuH aggregation, how it influences overall catalyst performance, and how it can be controlled. We address these challenges with mechanistic studies on a model reaction of unactivated alkene hydroboration catalyzed by (IPr*CPh 3 )CuH (LCuH). Here, we report a comprehensive mechanistic investigation of this system, identifying an aggregation pathway that continuously depletes catalytically active LCuH to form inactive CuH clusters during turnover. Deactivation of LCuH is controlled primarily by the competition between the kinetics of the initial LCuH dimerization step and that of alkene insertion. We therefore propose that a more comprehensive understanding of CuH catalyst performance must account for the kinetics of the initial LCuH dimerization step, revising a previously explored thermodynamic understanding of CuH aggregation, where the concentration of active species is controlled by equilibria established between CuH dimers and monomers. With a series of (NHC)CuH congeners (NHC = N-heterocyclic carbene), we demonstrate that ostensibly minor structural modifications to the ligand peripheries can drastically affect the LCuH dimerization kinetics, while maintaining reactivity towards on–cycle alkene insertion. We employed a computational approach based on molecular dynamics simulations to provide an in-depth understanding of how specific structural ligand modifications can substantially increase the kinetic stability of monomeric CuH catalysts. Our combined experimental and computational studies suggest strategies for rational ligand design that can be broadly applied to molecular catalyst systems that are susceptible to deactivation via aggregation pathways.

Ryan, David E. [Pacific Northwest National Laborat↗

Preparation of Neptunyl and Plutonyl Acetates To Access Nonaqueous Transuranium Coordination Chemistry

Uranyl diacetate dihydrate is a useful reagent for the preparation of uranyl (UO 2 2+ ) coordination complexes, as it is a well-defined stoichiometric compound featuring moderately basic acetates that can facilitate protonolysis reactivity, unlike other anions commonly used in synthetic actinide chemistry such as halides or nitrate. Despite these attractive features, analogous neptunium (Np) and plutonium (Pu) compounds are unknown to date. Here, in this study, a modular synthetic route is reported for accessing stoichiometric neptunyl(VI) and plutonyl(VI) diacetate compounds that can serve as starting materials for transuranic coordination chemistry. The new NpO 2 2+ and PuO 2 2+ complexes, as well as a corresponding molecular UO 2 2+ complex, are isomorphous in the solid state, and in solution show similar solubility properties that facilitate their use in synthesis. In both solid and solution state, the +VI oxidation state (O.S.) is maintained, as demonstrated by vibrational and optical spectroscopy, confirming that acetate anions stabilize the oxidizing, high-valent +VI states of Np and Pu as they do for the more stable U(VI). All three acetate salts readily react with a model diprotic ligand, affording incorporation of U(VI), Np(VI), and Pu(VI) cores into molecular coordination compounds that occurs concomitantly with elimination of acetic acid; the new complexes are high-valent, yet overall charge neutral, facilitating entry into nonaqueous chemistry by rational synthesis. Computational studies reveal that the dianionic ligand framework assists in stabilizing the +VI O.S. via donation to the 5f shells of the actinides, highlighting the potential usefulness of protonolysis reactivity toward preparation of stabilized high-valent transuranic species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Separation and Transition on the ROTEX-T Cone-Flare

As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.

CFD Validation↗

Separation and Transition on the ROTEX-T Cone-Flare

As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.

CFD Validation↗

A Pd III Sulfate Dimer Initiates Rapid Methane Monofunctionalization by H Atom Abstraction

An electrogenerated Pd III 2 species in fuming sulfuric acid is competent for rapid and concurrent methane monohydroxylation to methyl bisulfate (CH 3 OSO 3 H) and methane sulfonation to methanesulfonic acid (CH 3 SO 3 H). In situ NMR at 50 °C is used to track methane transformation exclusively to CH 3 OSO 3 H and CH 3 SO 3 H at high conversions. Integrating a set of kinetic and computational studies, the mechanism of methane monofunctionalization by Pd III 2 is examined. Here, experimental rate laws and common kinetic isotope effects for CH 3 OSO 3 H and CH 3 SO 3 H formation suggest that both transformations proceed via a common rate-limiting C-H activation step. Introduction of O 2 or Pd II,III 2 suppresses CH 3 SO 3 H generation, indicating a radical chain sequence. Although the metal-metal bonded Pd III 2 complex is a net two-electron oxidant, our aggregate kinetic data point to a mechanistic model that features rate-limiting H atom abstraction by the Pd III 2 complex to generate a methyl radical intermediate. The CH 3 • intermediate then recombines with Pd II,III 2 to furnish a CH 3 Pd III 2 intermediate that reductively eliminates CH 3 OSO 3 H. Alternatively, the CH 3 intermediate can enter a chain reaction with SO 3 to generate CH 3 SO 3 H. DFT computations support the radical-based C-H activation by Pd III 2 and delineate H atom abstraction pathways with computed reaction barriers and kinetic isotope effects (KIEs) that are consistent with experimental data. These mechanistic investigations challenge the paradigm of electrophilic C-H activation and highlight H atom abstraction as a potent pathway for selective methane C-H oxidative functionalization at high reaction rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Building Load Control Using Distributionally Robust Chance-Constrained Programs with Right-Hand Side Uncertainty and the Risk-Adjustable Variants

Aggregation of heating, ventilation, and air conditioning (HVAC) loads can provide reserves to absorb volatile renewable energy, especially solar photo-voltaic (PV) generation. In this paper, we decide HVAC control schedules under uncertain PV generation, using a distributionally robust chance-constrained (DRCC) building load control model under two typical ambiguity sets: the moment-based and Wasserstein ambiguity sets. We derive mixed integer linear programming (MILP) reformulations for DRCC problems under both sets. Especially, for the Wasserstein ambiguity set, we use the right-hand side (RHS) uncertainty to derive a more compact MILP reformulation than the commonly known MILP reformulations with big-M constants. All the results also apply to general individual chance constraints with RHS uncertainty. Furthermore, we propose an adjustable chance-constrained variant to achieve tradeoff between the operational risk and costs. We derive MILP reformulations under the Wasserstein ambiguity set and second-order conic programming (SOCP) reformulations under the moment-based set. Using real-world data, we conduct computational studies to demonstrate the efficiency of the solution approaches and the effectiveness of the solutions. Summary of Contribution: The problem studied in this paper is motivated by a building load control problem that uses the aggregation of heating, ventilation, and air conditioning (HVAC) loads as flexible reserves to absorb uncertain solar photovoltaic (PV) generation. The problem is formulated as distributionally robust chance-constrained (DRCC) programs with right-hand side (RHS) uncertainty. In addition, we propose a risk-adjustable variant of the DRCC programs, where the risk level, instead of being predetermined, is treated as a decision variable. The paper aims to provide tractable reformulations and solution algorithms for both the (general) DRCC and the (general) adjustable DRCC models with RHS uncertainty.

97 MATHEMATICS AND COMPUTING↗

Diastereoselectivity Controlled by the Hydrogenation Mechanisms during the Electrochemical Reduction of a Carbonyl Group

Stereocontrol is of critical importance in organic synthesis. In this study, we demonstrate how heterogeneous electrochemical hydrogenation enables diastereocontrol simply by tuning electrochemical hydrogenation mechanisms without altering the adsorption conformation of a reactant on the electrode. We use 4-hydroxy-1-tetralone (4-OH-tetralone) as a model reactant, where diastereomers can be produced during the hydrogenation of the carbonyl group. In traditional thermocatalytic hydrogenation, H 2 first dissociates on the catalyst surface to form surface-adsorbed hydrogen (H*), and therefore, H* is always added to the organic reactant from the catalyst side via hydrogen atom transfer (HAT). Thus, in order to flip the diastereoselectivity, the adsorbed reactant itself must be physically flipped. In contrast, electrochemical hydrogenation can occur either via HAT, where H is added from the electrode surface, or via proton-coupled electron transfer (PCET), where H is added from the solution side of the adsorbed reactant. Thus, without changing the adsorption conformation of the reactant, opposite diastereomers can be obtained by switching the hydrogenation mechanism (HAT vs PCET). In this work, using a combination of experimental and computational studies, we demonstrate two examples of flipping diastereoselectivity by different electrochemical hydrogenation mechanisms. In the first case, we achieve opposite diastereoselectivities using metals that adopt different hydrogenation mechanisms (HAT vs PCET). In the second case, we flip the diastereoselectivity by varying the applied potential, which switches one hydrogenation mechanism to the other on the same metal electrode. In each case, our results offer an atomic-level understanding of the preferred hydrogenation mechanism that enables the corresponding diastereoselectivity.

adsorption↗

Effects of detailed geometry and real fluid thermodynamics on Spray G atomization

We present recent results toward the quantification of spray characteristics at engine conditions for an eight-hole counter-bored (stepped) GDI injector-Spray G in the ECN denomination. This computational study is characterized by two novel features: the detailed description of a real injector's internal surfaces via tomographic reconstruction; and a general equation of state that represents the thermodynamic properties of homogeneous liquid-vapor mixtures. The combined level-set moment-of-fluid approach, coupled to an embedded boundary formulation for moving solid walls, makes it possible to seamlessly connect the injector's internal flow to the spray. The Large Eddy Simulation (LES) discussed here presents evidence of partial hydraulic flipping and, during the closing transient, string cavitation. Results are validated by measurements of spray density profiles and droplet size distribution.

42 ENGINEERING↗

The Effect of Wake Passing on Turbine Blade Film Cooling

The effect of upstream blade row wake passing on the showerhead film cooling performance of a downstream turbine blade has been investigated through a combination of experimental and computational studies. The experiments were performed in a steady-flow annular turbine cascade facility equipped with an upstream rotating row of cylindrical rods to produce a periodic wake field similar to that found in an actual turbine. Spanwise, chordwise, and temporal resolution of the blade surface temperature were achieved through the use of an array of nickel thin-film surface gauges covering one unit cell of showerhead film hole pattern. Film effectiveness and Nusselt number values were determined for a test matrix of various injectants, injectant blowing ratios, and wake Strouhal numbers. Results indicated a demonstratable reduction in film effectiveness with increasing Strouhal number, as well as the expected increase in film effectiveness with blowing ratio. An equation was developed to correlate the span-average film effectiveness data. The primary effect of wake unsteadiness was found to be correlated well by a chordwise-constant decrement of 0.094-St. Measurable spanwise film effectiveness variations were found near the showerhead region, but meaningful unsteady variations and downstream spanwise variations were not found. Nusselt numbers were less sensitive to wake and injection changes. Computations were performed using a three-dimensional turbulent Navier-Stokes code which was modified to model wake passing and film cooling. Unsteady computations were found to agree well with steady computations provided the proper time-average blowing ratio and pressure/suction surface flow split are matched. The remaining differences were isolated to be due to the enhanced mixing in the unsteady solution caused by the wake sweeping normally on the pressure surface. Steady computations were found to be in excellent agreement with experimental Nusselt numbers, but to overpredict experimental film effectiveness values. This is likely due to the inability to match actual hole exit velocity profiles and the absence of a credible turbulence model for film cooling.

Heidmann, James David↗

Nuclear Fuel and Pu Redox Studies from The Glenn T. Seaborg Institute at Idaho National Laboratory

The Glenn T. Seaborg Institute at Idaho National Laboratory (INL-GTSI) focuses on advancing fundamental research in the actinide sciences by providing unique opportunities to early career scientists and engineers to gain experience studying the actinide elements and their associated systems. The INL-GTSI is built from three focus areas that are based on the expertise and supporting infrastructure at INL and include solid state chemistry and physics, solution phase chemistry and physics, and forensic and isotope science. INL is the lead Laboratory for nuclear energy research and development in the U. S. and the research on nuclear fuels performed under the INL-GTSI gives good examples of solid state studies. Uranium-Molybdenum (U-Mo) alloys are leading fuel candidates for conversion of high performance research and test reactors to low-enriched fuels. During irradiation, generated fission gas accumulates into bubbles and self-organizes into a gas bubble superlattice (GBS) that effectively stores fission gases and inhibits fuel swelling. A study on the early self- organizing behavior of the GBS shows that not only grain boundaries but the interfaces between the U-Mo matrix and uranium carbide (UC) impurities are important to GBS formation.[1] In solution, understanding the complex redox behavior of plutonium in aqueous environments is critical for establishing optimized nuclear waste reprocessing solvent systems and storage tank environments. INL-GTSI researchers have produced an experimentally validated multi-scale model of the gamma radiation induced behavior of plutonium ions in concentrated aqueous HNO3 solutions.[2] Here, gamma radiation effected only minimal steady state changes in the redox distribution of the plutonium oxidation states. The redox cycling between Pu(IV) and Pu(III) is demonstrated to be mediated by the •OH/NO3• radical oxidation of Pu(III) and the H2O2/HNO3 driven reduction of Pu(IV). The INL-GTSI offers young researchers the unique chance to work directly with actinide bearing materials in a U. S. National Laboratory environment. Further topical areas of interest to the INL-GTSI include, but are not limited to, fundamental actinide properties, structure/property (electronic, magnetic, thermal) relations, actinide quantum criticality, f- electron interactions, electron correlations, computational studies, new phases, defect effects, interface interactions, isotope production and separation, forensic analytical chemistry, structure and dynamic properties of actinides in non-aqueous media, separations chemistry and kinetics for advanced nuclear fuel cycles, radiation effects, and innovative and advanced ligand design for complexation of the actinides.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗

Aerothermodynamic Measurements in Hypersonic Non-Equilibrium Flows

High enthalpy arc-jets are unique facilities particularly suited for simulating complex flows in the aerospace field, such as the aerothermodynamics of a re-entry vehicle. Arc-jets are often used to evaluate important design factors such as heat shield materials and vehicle design. Characterization of these facilities is important, as studies often match specific in-flight environments during experiments. Due to the complex environment produced by an arc-jet, with effects such as thermodynamic and chemical non-equilibrium occurring in the flow, characterization experiments are significantly more difficult than in traditional blow-down wind tunnels. The current work aims to characterize an arc-jet facility in terms of flow temperature, velocity and unsteadiness, in a spatially-resolved analysis. To achieve this goal, a non-intrusive imaging technique called “planar laser-induced fluorescence” was performed in the NASA Langley Hypersonic Materials Environ-mental Test System arc-jet facility. The experimental data was analyzed to produce the quantitative measurements in multiple regions of the flow around a blunt body sample. A quasi-1D low fidelity solver was created to numerically simulate the flow to investigate non-equilibrium effects occur-ring outside the imaging region in the arc-jet nozzle. Unsteadiness in the test section was minimized by analyzing a subset of data assessing the gas injection location. Radial velocity, rotational temperature and translational temperature measurements are provided that can be used to validate future high-fidelity computational studies. The temperature measurements revealed rotational non-equilibrium potentially occurring in the flow. This work provides a methodology for characterizing arc-jet facilities, identifying and correcting sources of systematic error, and estimating the uncertainty associated with the quantitative measurements.

laser induced fluorescence↗

Reticular Materials and AI-Driven Computer Simulations for Seawater Mining of Valuable Metals (Final Technical Report)

This Final Technical Report describes our exploratory efforts that combine reticular materials synthesis (hydrolytically robust metal–organic frameworks, MOFs) with AI‑enabled molecular simulations to develop mechanistic, quantitative design rules for recovering lithium and other alkali-metal ions from highly dilute, competitive aqueous resources (e.g., seawater). The central outcome is a joint experimental–computational study of ion uptake in MOF‑808 (Chemical Science, 2025) that quantifies both thermodynamics and kinetics of Li + , Na + , and K + uptake and identifies how pore size, pore hydration state, dehydration penalties, and pore-window transport barriers govern selectivity. Guided by these insights, we synthesized and tested functionalized MOF‑808 and multivariate MOFs incorporating ion-recognition motifs (including carboxylates and crown-ether linkers) and evaluated uptake in synthetic seawater, highlighting framework topology and pore chemistry as levers for improved Li + /Na + discrimination. We also developed transferable simulation models, enhanced-sampling protocols, and automated workflows that enable systematic screening of porous sorbents.

42 ENGINEERING↗

Enhanced resistance to decay of imprinted nanopatterns in thin films by bare nanoparticles compared to polymer-grafted nanoparticles

We extend a previous study on the influence of nanoparticles on the decay of nanoimprinted polymer film patterns to compare the effects of “bare” silica (SiO 2 ) nanoparticles and SiO 2 nanoparticles with grafted polymer layers having the same chemical composition as the polymer matrix. This method involves nanoimprinting substrate-supported polymer films using a pattern replicated from a digital versatile disc (DVD), and then annealing the patterned polymer nanocomposite films at elevated temperatures to follow the decay of the topographic surface pattern with time by atomic force microscopy imaging after quenching. We quantified the relaxation of the pattern height (“slumping”) and determined the relaxation time τ for this pattern decay process as a function of nanoparticle filler type and concentration to determine how nanoparticle additives influence relative film stability. Attractive interactions between the bare nanoparticles and the polymer matrix significantly enhance the thermal resilience of the nanopatterns to decay, compared to those of the particle brushes, wherein the particle core interactions are screened from the matrix via the brush layer. A novel aspect of this method is that it readily lends itself to in situ film relaxation measurements in a manufacturing context. We observe that the relaxation time of the pattern relaxation exhibits entropy–enthalpy compensation in the free energy parameters governing the pattern relaxation process as a function of temperature, irrespective of the NP system used, consistent with our previous experimental and computational studies.

36 MATERIALS SCIENCE↗