Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Rotation operations”

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.

627 records · Page 35

Influence of Proton Activity Gaps between Electrodes on Open-Circuit Potential of H 2 /H 2 and H 2 /Air Cells

Polymer -electrolyte-fuel-cell open-circuit voltages (OCVs) are exactly defined by equation (1), where cathode and anode proton activities [(aH+)cathode and (aH+)anode, respectively] usually are identical, so the third term in the right-hand side of equation (2) is ignored. OCV=E0+RT/2F*ln(a1/2 O2*(a2 H+)cathode/aH2O)-RT/2F*ln((a2 H+)anode/aH2) (1) =E0+RT/2F*ln(a1/2 O2*aH2/aH2O)+RT/2F*ln((a2 H+)cathode/(a2 H+)anode) (2) Water vapor pressure is a colligative property that fundamentally correlates to electrolyte concentrations in aqueous solutions. Proton activity is a function of acid concentration, such as pH, when electrolytes are acids. In polymer-electrolyte membranes, water vapor pressure and acid concentration are understood as relative humidity (RH) and water uptake (λ), respectively, where λ represents number of water molecules per sulfonic acid molecule. Several investigations have reported the relationship between RH and λ, meaning that proton activities and associated water uptakes are intimately related to RH. In actual fuel cell operation, cathode RH is determined by ambient-atmosphere and/or humidifier RH(s), and anode RH depends on hydrogen-circulator RH. Therefore, RH is not always identical at both electrodes, and the difference between electrode RHs is considerable during dry operation of polymer electrolyte fuel cells. Therefore, the third term in the right-hand side of equation (2) may be significant for dry operation. We measured OCVs when hydrogen was supplied to both electrodes at 80°C. One electrode (A) was fixed at 30% RH, while RH at the other electrode (B) was varied (0, 5, 10, 20, and 30%). Measured OCVs varied from 0 to 75 mV. For fuel cell tests, electrode A was supplied with hydrogen at 30% RH; electrode B, oxygen at 0, 5, 10, 20, and 30% RH. OCVs deviated from that measured when RH at electrode B was 30%, increasing from 0 to 60 mV with decreasing RH at electrode B. Results are also shown in Figure 1. Proton activities of both electrodes were thermodynamically calculated. The Gibbs–Duhem relation was applied to obtain molar Gibbs free energies of water and sulfonic acid, and proton activity coefficient was calculated using the Gibbs free energy of sulfonic acid and the relationship between RH and λ1–4, assuming that protons and sulfonic anions show identical ionic-activity coefficients. OCVs were estimated using the third term in the right-hand side of equation (2). Results are shown in Figure 1. Fuel-cell current–voltage performance was poor when RHs at the anode and cathode were 30 and 20%, respectively. To determine kinetic current, we measured the oxygen-reduction reaction (ORR) using a rotating-disk electrode (RDE) in concentrated-acid aqueous solutions, which modeled catalyst-layer ionomers. Kinetic currents decreased with acid concentrations. References T. A. Zawodzinski, Jr., C. Derouin, S. Radzinski, R. J. Sherman, V. T. Smith, T. E. Springer and S. Gottesfeld , J. Electrochem. Soc., 140,1041 (1993) P. K. Das and A. Z. Weber, Proceedings of the ASME 2013 11th Fuel Cell Science, Engineering and Technology Conference, Fuel Cell 18010 (2013) V. A. Sethuraman, J. W. Weidner, A. T. Haug, S. Motupally,b and L. V. Protsailo, J. Electrochem. Soc., 155, B50 (2008) A. Kusoglu and A. Z. Weber, Chem. Rev., 117, 987 (2017) Figure 1

Yoshida, Toshihiko↗

Thermal Management for a Novel Non-Heavy Rare-Earth Interior Permanent Magnet Machine

The work presents a thermal management solution for a novel non-heavy rare-earth permanent magnet machine being developed at Oak Ridge National Laboratory. The motor has been designed to minimize losses while maximizing performance for a range of speeds and power ratings. The novel motor design reduces rare-earth magnet usage, thereby avoiding supply chain issues. The motor component heat losses are established for operating windows and desired performance. These heat losses, along with windage losses, are being used to develop cooling solutions for different components of this machine. A novel thermal management solution for stators and rotors has been developed, and progress is presented in this paper. The stator cooling is achieved with the help of water-ethylene glycol flowing over the finned aluminum stator jacket, and rotor cooling with automatic transmission fluid passing through novel channels designed in the rotor laminations. The attempt is to establish effective cooling of the stator winding, laminations, and rotor magnets. A 3D conjugate heat transfer model has been developed for overall thermal analysis to establish a down- selected thermal management solution for the machine. The model, in addition to estimated component heat losses, includes windage losses and its impact on rotor and stator cooling. Overall, the work presents a workable thermal solution for the interior permanent magnet machine with potential for further improvements. Future work will involve establishing end winding and refinement of other end parts of the machine with the aim of establishing a robust thermal management solution. The work will also focus on different shapes (e.g., round, non-round, presence of wedges) of rotor-stator gaps and investigate windage losses and their impact on thermal management for higher rotational speeds for the machine.

30 DIRECT ENERGY CONVERSION↗

Investigation of Technologies to Improve Condenser Heat Transfer and Performance in a Relevant Coal-Fired Power Plant

Improvements in thermal power-generating plant performance is correlated directly to societal benefits including lower cost of reduced fuel consumption, resulting in lower cost of electricity for the consumer and reduced carbon emissions to the atmosphere. Warm steam exhausted from low-pressure steam turbines is condensed to liquid water on the exterior of thin-walled metal condenser tubes with cooling water passing through the tube interior. The condensation of steam creates a vacuum that supports turbine rotation and the concurrent generation of electricity. This vacuum is optimized when heat transfer across the wall of condenser tube is maximized. Common hindrances to heat transfer include foulants in cooling water that may form and adhere to the interior of condenser tubes, including mineral scale, microbiological films, and particulate deposition. Flowing cooling water may also include a laminar layer at the interior metal surfaces that travels more slowly than bulk water flow, serving to impede heat transfer. On the tube exterior, condensing steam forms an insulating layer of water that flows down the tube and reduces the effectiveness of cooling. Both the interior and exterior barriers to optimal heat transfer may be alleviated to some extent by surface treatments. On the tube interior, hydrophobic coatings may be applied that can reduced the adherence of foulants and of the laminar flow layer to the tube surface. On the tube exterior, hydrophobic coatings or mechanical treatments can be applied that may result in the termination of droplet growth and the departure of droplets from the surface rather than coalescence into a continuous layer of flowing water. Fourteen surface treatments were applied to condenser tubes in this study, including eight interior coatings and six exterior treatments, five of which were coatings and one a microstructural texture. Heat transfer measurement equipment simulating conditions in the condenser of an operating power plant was used to determine heat transfer coefficients by measuring sufficient flow, temperature, water chemistry and other data. Several of the tubes with interior surface treatments showed improvement in heat transfer coefficients compared with a plain (uncoated) tube, and several of the tubes with exterior surface treatments also showed enhanced heat transfer coefficients.

01 COAL, LIGNITE, AND PEAT↗

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Consideration of Nonideal Detonation Regimes Influenced by Wave Modes in a Water-Cooled Rotating Detonation Engine Using OH* Chemiluminescence

Although inherently unstable, existing research in rotating detonation combustion supports its application in notionally steady processes resulting in greater availability compared to conventional, constant pressure combustion. Further improvements rely on a more in-depth understanding of system losses and identifying conditions which optimize device performance. Within this study, the presence and proportion of ideal and nonideal combustion regimes are compared across a variety of process conditions and wave modes. Here, large-scale data analysis seeks to summarize proportional heat release associated with commensal, parasitic, and detonative combustion averaged across individual traces of OH* chemiluminescent data acquired at the detonation plane. Means of regime partitioning based on the anatomy of the time-resolved OH* signal are proposed to ensure consistent analysis throughout the current and future studies concerning combustion regimes. Of particular interest is the possible influence of wave on the nonideal combustion in relative proportion to the desired detonation. Results showed improved percent detonation with increasing significance for the following trends: decreasing equivalence ratio, increasing wave count, decreasing wave velocity, and increasing detonation time. Increased wave number, brought on by decreased equivalence ratios and wave velocities, is thought to decrease fill region surface area, and therefore, decrease nonideal contact burning. Additional analysis is performed to consider possible trend variation due to the presence of stable galloping waves, which were found to have minimal influence on relative percent detonation values. The outcome of this study suggests operational states, which correspond to increased wave quantities for increased proportions of reactants consumed by the targeted detonative combustion regime.

42 ENGINEERING↗

Coupled Experimental and High-Temperature Discrete-Element Method Modeling Studies of Aluminosilicate Particle Handling in Concentrated Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrated solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrated solar power environments. Coupled experimental and numerical studies of aluminosilicate particles in rotary kilns and dense particle curtains were performed for bulk temperatures up to 1073 K. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Temperature, spatial, and velocity profile data were extracted from experimental runs using embedded K-type thermocouple probes and particle image velocimetry techniques. Experimental and numerical studies were compared using spatial temperature profiles, velocity fields, and shape profiles of the bulk, granular flows. Numerical models were developed using commercially available discrete element method modeling software, Aspherix®. Existing Aspherix® functionality was expanded by introducing coupled radiative exchange modeling tools. The laboratory-scale rotary kiln was developed to investigate the steady-state heat and mass transfer performance of aluminosilicate particles based on particle type, bulk handling temperature, and wall roughness. The rotational speed of the rotary kiln was varied to control the relative impact of collisional and frictional effects upon the granular flow behavior. Heat and mass transfer performance was categorized based on the Froude number and the observed flow regimes of slipping, rolling, cascading, and centrifuging. Coupled discrete element method modeling studies were used to evaluate the effects of temperature-dependent, particulate mechanical properties upon bulk flow behavior and upon the relative effects of radiative, advective, and/or conductive heat transfer. A high-temperature (< 1073 K) falling particle curtain was similarly fabricated to investigate the heat and mass transfer performance of aluminosilicate particles in particle handling situations dominated by inter-particle collisions. The impact of particle type, flow preheat temperatures (< 1073K), and bulk mass flow rates were investigated upon the particle curtain shape, temperature, and velocity profiles. Coupled discrete element method modeling studies were performed to evaluate the varying impact of temperature-dependent, particulate mechanical properties on the bulk flow behavior and the temperature profile of the particle curtain.

14 SOLAR ENERGY↗

Boosted Dark Matter Search with the ICARUS Detector at the Gran Sasso Underground National Laboratory

The quest to discover dark matter is one of the most pressing topics in modern physics. Evidence of gravitational anomalies is present at different cosmological scales, spanning the stability and rotation of galaxies to the overall density profile of the entire observable Universe. Due to the discrepancy between the mass measured and the gravitational influence on the mass, if General Relativity is to correctly describe gravity at different scales, extra electromagnetically invisible matter termed dark matter must exist. In recent decades, many experiments have been performed to detect dark matter with masses that span many orders of magnitude, from Mass Compact Halo Objects (MACHOs) with masses several orders of magnitude that of our sun to very light particles on the order of $\sim\mu$eV. However, experiments have only been able to constrain the possible parameter space of dark matter models and, in general, probe dark matter models with single constituent dark matter like Weakly Interacting Massive Particles (WIMPs). In the Boosted Dark Matter (BDM) paradigm, there are four particles in the dark sector: $\chi_0$, $\chi_1$, $\chi_2$, and the dark photon $X$. WIMP-like properties are ascribed to $\chi_0$ to describe the physics of the gravitational anomalies, while relativistic components are assigned to $\chi_1$ and $\chi_2$. The inelastic boosted dark matter (iBDM) channel $\chi_1 e^-\rightarrow \chi_2 e^-\rightarrow \chi_1 Xe^-\rightarrow \chi_1 (e^+ e^- ) e^-$, where $\chi_2$ is an excited state of $\chi_1$, has a unique signature that is distinguishable from neutrino interactions, making neutrino experiments a viable search environment for BDM. Using data collected by the ICARUS collaboration with ICARUS T-600 Liquid Argon Time Projection Chamber (LArTPC) during its operation at the Gran Sasso Underground Laboratory (LNGS), the first-ever LArTPC dark matter search using the ICARUS detector is presented. With a total detector exposure of 0.13 kton$\cdot$year, a total of 4 iBDM candidate events were identified by several stages of data scanning. All were rejected due to defined selection criteria and a wire-by-wire $dE/dx$ evaluation, giving a null result. Exclusion plots at 90\%C.L. for the dark photon ($m_X$, $\epsilon$) parameter space for fixed ($m_0,m_1,m_2$) DM mass sets are presented, setting new limits on the dark photon visible decay parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Entangling-gate error from coherently displaced motional modes of trapped ions

Entangling gates in trapped-ion quantum computers are most often applied to stationary ions with initial motional distributions that are thermal and close to the ground state, while those demonstrations that involve transport generally use sympathetic cooling to reinitialize the motional state prior to applying a gate. Future systems with more ions, however, will face greater nonthermal excitation due to increased amounts of ion transport and exacerbated by longer operational times and variations over the trap array. In addition, pregate sympathetic cooling may be limited due to time costs and laser access constraints. In this paper, we analyze the impact of such coherent motional excitation on entangling-gate error by performing simulations of Mølmer-Sørenson (MS) gates on a pair of trapped-ion qubits with both thermal and coherent excitation present in a shared motional mode at the start of the gate. Here, we quantify how a small amount of coherent displacement erodes gate performance in the presence of experimental noise, and we demonstrate that adjusting the relative phase between the initial coherent displacement and the displacement induced by the gate or using Walsh modulation can suppress this error. We then use experimental data from transported ions to analyze the impact of coherent displacement on MS-gate error under realistic conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modelling pulsed field magnetization of iron-based bulk superconductors

Abstract Bulk superconductors can be used as super-strength quasi-permanent magnets capable of providing magnetic flux densities considerably superior to conventional permanent magnets. This makes them attractive for several engineering applications that rely on strong magnetic fields like rotating machines, NMR/MRI and magnetic drug delivery systems. Recently, the authors reported a record trapped magnetic field in an iron-based bulk superconductor: 2.83 T was trapped in potassium-doped barium iron arsenide (Ba, K)Fe 2 As 2 (or Ba122) at 5 K. Of particular significance is that the strength and temporal stability of this magnetic field exceeds the requirements of MRI machines, indicating iron-based bulks can now perform at levels demanded by engineering applications. One crucial challenge for their practical use, however, is the need to apply and remove an external magnetic field to magnetize them. Pulsed field magnetization (PFM) shows great promise as a practical method of magnetizing bulks, but the process generates heat in the bulk that is detrimental to its superconducting performance and ability to act as a super-strength magnet. In this paper, coupled electromagnetic–thermal numerical models are used to simulate the PFM of iron-based bulk superconductors. Here we focus on the recent-record-breaking, fine-grain polycrystalline K-doped Ba122 bulks. The impact that the specific J c ( B ) characteristics and thermal properties of the Ba122 material—all of which have been experimentally measured from state-of-the-art samples—have on the magnetic flux dynamics and thermal behaviour during PFM, including the final trapped field, is investigated. We show that because the thermal properties are similar to those of REBa 2 Cu 3 O 7 −δ bulks, a similar response to pulsed fields is obtained. A maximum trapped field of ∼0.81 T (∼43.3% of the maximum trapped field capability under ideal, field-cooling conditions) was simulated at 5 K, with a magnetization efficiency of ∼54%. The modelling framework provides a fast and flexible tool for optimising the practical PFM process at different operating temperatures to maximise the trapped field in state-of-the-art Ba122 bulks and to guide the design of future experiments.

bulk superconductors↗

In-Situ Blade Strain Measurements of a Crossflow Turbine Operating in a Tidal Flow

This data was collected between October 25 and December 12 of 2022 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). The goal was to collect blade strain data from a crossflow turbine operating in a tidal flow. A table in ('Deployment Schedule.PNG') outlines the various dates when each instrument was operational, and more details can be found via literature listed in 'Related Publications'.txt. This dataset includes zipped folders for each instrument containing data in .csv files for the relevant duration specific to each instrument, along with separate README file for each measurement. Some instrument files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. All data contained in this submission is unfiltered/unprocessed data unless otherwise noted in the README file. Blade strain was measured using 8 foil-based strain gauges along the span of a single turbine blade. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) both upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal flow with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Chancey 2019. The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. A single blade of this turbine was further modified to accommodate 8 full-bridge strain gauges (Bharath et al 2023, Bichanich et al 2024). For power performance and other relevant details on the turbine and its characteristics, see O'Byrne 2022.

16 TIDAL AND WAVE POWER↗

Charge per Micro-Pulse Calculation for Ideal Lujan and WNR beams

A complete redesign of the LANSCE front-end is currently under development for the Los Alamos Modernization Project (LAMP). This includes the replacement of the traditional Cockroft-Walton injection system to the newer radiofrequency quadrupole (RFQ) standard. LANSCE accelerates both H- and H+ beams and therefore requires an unconventional application for RFQ systems. A new facility could consider using an independent RFQ for each beam species. However, with the tight space of the LANSCE injector facility, it would prove difficult to instrument. Our initial studies have shown that a single RFQ could be used. The negatively- and positively-charged direct current (DC) beams are shaped, bunched, and accelerated out of phase to each other, with very little interaction. This would work for all possible H+ beam gates and most H- Beam gates. However, a challenge arises with the Weapons Neutron Research (WNR) beam structure that requires the acceleration of bunched beams at higher charge density per pulse. The Long Bunch Enable Gate (LBEG) used for the Lujan center is the best example of a standard injector beam gate. The beam gate, or Macro-Pulse (MP), is 625-us long, separated by at least 8.3 ms. A chopper is used to dice the MP into 1750 mini-Pulses (mP), which are 290-ns wide. After chopping, the MP contains 100,000 micro-Pulses (uP), which are generated by the linac acceleration structure modulated at 201.25 MHz; that is, the micro-Pulses are separated 5 ns apart from each other. The H+ beam structure will be very similar once a chopper is added to its transport. However, the Multi-Pulse Enable Gate (MPEG) used for WNR is significantly different in structure. In this case, the mP is chopped to be only 35-ns wide. This width is then rotated into a single 201.25 MHz RF bucket using the Low Frequency Buncher (LFB). This is proving difficult for the RFQ design in the LAMP upgrade. We began investigating the amount of current sent to WNR for two reasons. Firstly, we wanted to find a solution to avoid using the Low Frequency Buncher to simplify the LAMP RFQ design. It is simple to inject and accelerate a 35 mA DC beam in the LAMP RFQ (this is done in numerous accelerator facilities). However, injecting a short bunch with a larger current (once compressed to 5 ns to meet the RFQ frequency) may be difficult. Second, we have a technique called PSR2WNR in which we do not employ the Low Frequency Buncher, but rather accumulate 5 LBEG uP (5 ns) in the PSR and deliver the accumulated beam to WNR every 1.8 us. As a result, we could provide to WNR with five accumulated LBEG uPs every shot. The goal here is to create a rapid kicker with a frequency response of 555 kHz. DARHT-II and other accelerator facilities have constructed kickers of this type for other uses, and our engineers consider this technology as feasible.

43 PARTICLE ACCELERATORS↗

Second Target Station High-Fidelity Target Activation Comparison

The development of the Second Target Station (STS) target system at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory (ORNL) is well underway. The target system at STS consists of a rotating target disk that contains 21 segments of tungsten clad in tantalum clad in steel. A key aspect of the design of the target system is to account for the delayed heating and material damage caused by the delayed dose from decaying radionuclides. These radionuclides are a product of either spallation reactions or transmutation of the nuclei in the target system. These radionuclides build up in the target system components over the lifetime of the facility, and the radiation that is emitted can deposit energy in the components causing significant component heating and material damage. Monte Carlo N-Particle (MCNP) Version 6.2 transports the various particle species and calculates the spallation products and neutron fluxes throughout the target system. These spallation products and neutron fluxes along with the material definition of each component are relayed to the CINDER2008 transmutation code to calculate the radionuclide inventories and the corresponding decay gamma emission spectra. MCNP6.2 coupled with CINDER2008 is the computational method-of-choice for the analysis discussed in the following sections of this report. The analysis focuses on validating major assumptions in calculating the radionuclide inventory in the STS target system: all of the target segments are fresh, unirradiated material when the protons are incident on the segment, the average of the 21 segments of the target is sufficient to represent a single segment, and that averaging the proton pulse structure over time does not significantly affect the radionuclide inventory. The position-averaged, high-fidelity, and single-tally computational methods are used to validate the assumptions and provide a point of comparison to evaluate how the assumptions impact the radionuclide inventories. A more detailed explanation of the three computational methods is provided in Section 2. The position-averaged and single-tally methods are less computationally expensive when compared with the high-fidelity method where 54,000 individual calculations are needed to calculate 1 hr of STS operation. Section 3 details the comparison of the three methods to show that the assumptions made in the position-averaged method do not significantly impact the radionuclide inventory after 1 hr of operation. The discussions and results in this report are for 1 hr of operation. Due to the computational cost associated with calculating the transmutation and activation using the high-fidelity method, only 1 hr of operation has been calculated. The discrepancies observed after 1 hr of operation are not extrapolated out to longer operational times, and this report does not address how the discrepancies between the computational methods may manifest for longer operational periods.

43 PARTICLE ACCELERATORS↗

Gripper Test Assembly Status Report (FY2020)

Refueling systems for fast reactors are designed to handle fresh and used core assemblies (fuel, reflector, and shield core assemblies) within the reactor vessel in an opaque coolant environment without visual reference. These refueling machines are designed to work in a sodium (or other fast reactor coolant) and argon vapor space environment and are engineered with the rotatable plug system to allow for the movement of fresh and spent fuel into and out of the reactor core. The refueling machines are a critical component in any reactor and thus need to undergo extensive testing in a prototypic environment to ensure that they will meet all of the system functions and requirements. Argonne has developed an innovative compact refueling system design for the Advanced Fast Reactor-100 that is based upon some mechanisms used in previous reactor designs, such as the U.K.’s Prototype Fast Reactor (PFR) and some mechanisms that have not been used in sodium. This compact refueling machine supports the reduction in size of the AFR-100’s reactor vessel, and if fully developed, would support and inform the development of the in-vessel refueling machines for such commercial reactors as the GEH PRISM reactor plant, the ARC Clean Energy’s ARC-100 reactor, and the Natrium reactor, among others. This refueling system is a vital component of a fast reactor that supports reducing the cost of the reactor and increasing its reliability. During the development of the compact fuel handling machine conceptual design for the AFR-100, a lack of testing data for many mechanical components in sodium under typical in-reactor loads and conditions was discovered. The reduction in lifetime of the various mechanical components in the liquid sodium environment needs to be quantified versus the calculated component lifetimes under normal conditions in the testing while they are subjected to typical loading profiles experienced in the past. The Gripper Test Assembly discussed here includes a full-size gripper device with appropriate mechanical features that will be tested in sodium to provide this testing data. The Gripper Test Assembly is used to test various mechanical fuel handling components submerged in high temperature liquid sodium. These mechanisms are gears, bearings, gripper jaws and head, universal joints and shafts, ball screws, among others. These components will be tested under the typical sodium environmental conditions experienced during refueling operations with appropriate loading conditions that simulate the removal and insertion of core assemblies into a fast reactor grid plate structure. This Gripper Test Assembly is the second in a series of refueling system mechanisms developed for testing in sodium. The first test assembly is the Gear Test Assembly (GTA) which was used to test the performance of gears and bearings operating in sodium. Because of the successful testing conducted with the GTA, it was decided to continue with the development of the Gripper Test Assembly which uses the same gears and bearings tested in GTA. Using the data collected during operation of this gripper test assembly, lifetime reduction factors of the various mechanical components can be calculated for the material combinations selected. These lifetime reduction factors can be used in the design of future mechanical systems which operate in these environments to accurately predict component end of life. In addition, an understanding of the ability of these components and mechanisms to operate under-sodium with the chosen materials will be accomplished. Again, this Gripper Test Assembly is the follow-on test article to the Gear Test Assembly which was testing the ability of gears and bearings operating in a sodium environment. Once fabricated and qualified, it will be tested in the Mechanism Engineering Test Loop (METL) facility located in Building 308 at Argonne National Laboratory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coherent anti-Stokes Raman scattering with squeezed light: CARS for quantum-enhanced spectroscopy and imaging

We theoretically investigate quantum-enhanced coherent anti-Stokes Raman scattering (CARS) using squeezed light to amplify vibrational transition rates at low photon flux. Quantum sensing approaches are needed for nondestructive nanometrology such as in bioimaging where reduced photodamage is desired while retaining resolution and sensitivity. We analyze both single-mode squeezing applied to the pump field and two-mode squeezing between the pump and Stokes fields. We also show that the ordering of displacement and squeezing operations—whether displacement precedes squeezing or squeezing precedes displacement—has an impact on the resulting CARS transition amplitudes due to a difference in the photon number and the quantum-enhancement coefficients, with the latter offering a stronger enhancement in the case of two modes squeezing of the pump and Stokes under experimentally accessible conditions. Furthermore, our calculations capture these quantum enhancements through the intrinsic photon-number correlations of squeezed light, eliminating the need for interferometric detection or higher pump powers that are otherwise required to reach comparable sensitivities in classical CARS. Finally, we outline a quantum plasmonic extension of our model in which local field enhancements caused by surface plasmon excitation in metallic nanoparticles can be incorporated via mode-selective field amplification factors, offering a pathway toward combining squeezed-light quantum optics with surface-enhanced nanoscale spectroscopy and imaging.

Atomic & molecular structure↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗