Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “average atom”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparison of Results from Recent NNSA and CEA Measurements of the 239 Pu(n, f) Prompt Fission Neutron Spectrum

The National Nuclear Security Administration (NNSA)/DP French Alternative Energies and Atomic Energy Commission (CEA)/DAM agreement on cooperation on fundamental science is a U.S.-French collaborative effort to combine intellectual and experimental resources and further the relevant nuclear science. Recently, both the NNSA and CEA experimental teams performed high-statistics measurements of the 239 Pu(n, f) prompt fission neutron spectrum (PFNS) at the Los Alamos Neutron Science Center, both of which were recently published in the journal Physical Review C. These separate measurements used the same experimental area and a common neutron detector array, but differ in many aspects, including background assessments, data acquisition systems and philosophies, fission detectors, and PFNS extraction techniques. Hence, some aspects of the experimental methods and associated uncertainties are highly correlated while others are independent. The results from both measurements broke new ground for PFNS measurements given their higher accuracy and more detailed study of corrections necessary for the measured quantity compared to existing literature measurements, and both will significantly impact PFNS nuclear data evaluations for the foreseeable future. Here, the focus of this work is to document a comparison of the results from these distinct measurements in terms of the acquired data, the PFNS results, and the measured average PFNS energies. While systematic differences between the PFNS results are present on the 1–3% level, the acquired data relative to each respective measurement at low incident neutron energies are in remarkable agreement, as are the conclusions regarding the magnitude and position of features in the PFNS relating to second-chance fission, third-chance fission, and pre-equilibrium neutron emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Upgrade and characterization of the fast-channels of Imaging Neutral Particle Analyzers in the DIII-D tokamak

This publication presents the characterization of the recently upgraded fast (MHz) channels of the DIII-D Imaging Neutral Particle Analyzers. The new system allows for a phase space resolution of < 10 keV and < 10 cm. The main source of noise is due to the background radiation (neutrons and gammas) present in the diagnostic lab during experiments, which can hit the photomultipliers and generate large spikes in signal. Good linear correlations (r 2 > 0.9) are found between the average rate of noise spikes in the fast channels and the average neutron rate. Additional noise is caused by the scintillation of the optical fibers used to guide the signal to the acquisition system, which is induced by neutron and gamma impacts. A correlation factor r 2 > 0.6 between the baseline noise caused by this scintillation and the signals in the neutron diagnostic was found. Possible mitigation strategies, such as the use of shielding and filters, are discussed. These allow to reduce the noise spike rate by an order of magnitude and the noise baseline by a factor 4. Despite the noise, it is shown that the upgraded fast channels provide enough signal-to-noise ratio to observe coherent fluctuations of the energetic particle population.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Inelastic scattering in isotopologues of O 2 –Ar: the effects of mass, symmetry, and density of states

The two species considered here, O 2 (oxygen molecule) and Ar (argon-atom), are both abundant components of Earth's atmosphere and hence familiar collision partners in this medium. O 2 is quite reactive and extensively involved in atmospheric chemistry, including Chapman's cycle of the formation and destruction of ozone; while Ar, like N 2 , typically plays the nevertheless crucial role of inert collider. Inert species can provide stabilization to metastable encounter-complexes through the energy transfer associated with inelastic collisions. The interplay of collision frequency and energy transfer efficiency, with state lifetimes and species concentrations, contributes to the rich and varied chemistry and dynamics found in diverse environments ranging from planetary atmospheres to the interstellar and circumstellar media. The nature and density of bound and resonance states, coupled electronic states, symmetry, and nuclear spin-statistics can all play a role. Here, we systematically investigate some of those factors by looking at the O 2 –Ar system, comparing rigorous quantum-scattering calculations for the 16 O 16 O– 40 Ar, 18 O 16 O– 40 Ar, and 18 O 18 O– 40 Ar isotope combinations. A new accurate potential energy surface was constructed for this purpose holding the O 2 bond distance at its vibrationally averaged distance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

XFEL SASE pulses can enhance time-dependent observables

Abstract X-ray free electron lasers (XFELs) have emerged as powerful sources of short and intense x-ray pulses. We propose a simple and robust procedure which takes advantage of the inherent stochasticity of self-amplified stimulated emission (SASE) pulses to enhance the time-resolution and signal strength of the recorded data. Notably, the proposed method is able to enhance the average signal without knowledge of the signal strength of individual shots. Simple metrics for the probe pulses are introduced, such as an effective pulse duration applicable to SASE pulses characterised in the time domain using e.g. an X-band transverse cavity. The approach is evaluated using simulated and real pulse data in the context of ultrafast electron dynamics in a molecule. Utilising H 2 as a model system, we demonstrate the efficacy of the method theoretically, successfully enhancing the predicted nonresonant ultrafast x-ray scattering signal associated with electron dynamics. The method presented is broadly applicable and offers a general strategy for enhancing time-dependent observables at XFELs.

data analysis↗
DOE OSTI2024

Resilience of the Aurivillius structure upon La and Cr doping in a Bi 5 Ti 3 FeO 15 multiferroic

Here, combining the experimental techniques of high-resolution X-ray diffraction, magnetometry, specific heat measurement, and X-ray photoelectron, Raman and dielectric spectroscopy techniques, we have studied the influence of La and Cr doping on the crystal structure and magnetism of the room temperature Aurivillius multiferroic Bi 5 Ti 3 FeO 15 by investigating the physical properties of (Bi 4 La)Ti 3 FeO 15 and Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 . The parent (Bi 5 Ti 3 FeO 15 ) and the doped ((Bi 4 La)Ti 3 FeO 15 and Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 ) compounds crystallize in the A2 1 am space group, which is confirmed through our analysis of high-resolution synchrotron X-ray diffraction data obtained on phase-pure polycrystalline powders. We determined the oxidation states of the metal atoms in the studied compounds as Fe 3+ , Ti 4+ , Cr 3+ , and La 3+ through the analysis of X-ray photoelectron spectroscopy data. The magnetic susceptibilities of the three compounds are marked by the absence of a long-range ordered ground state, but dominated by superparamagnetic clusters with dominant antiferromagnetic interactions. This signature of short-range magnetism is also seen in specific heat as a low temperature enhancement which is suppressed upon the application of external magnetic fields up to 8 T. Our dielectric spectroscopy experiments showed that the three studied compounds display similar features in the dielectric constant measured as a function of frequency. However, upon doping La at the Bi site, the width of the ferroelectric hysteresis loop increases for (Bi 4 La)Ti 3 FeO 15 compared to that of the parent compound Bi 5 Ti 3 FeO 15 , and with Cr doping, Bi 5 Ti 3 (Fe 0.5 Cr 0.5 )O 15 becomes a leaky dielectric. The resilience of the Aurivillius crystal structure towards doping of La at the Bi site and Cr at the Fe site is clearly seen in the bulk properties of magnetic susceptibility, specific heat and the average crystal structure. The relevance of changes in the local structure is evident from our Raman spectroscopy and X-ray pair distribution function studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

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

Genetic-algorithm retrieval of the molecular alignment distribution with high-order harmonics generated from transiently aligned CO 2 molecules

Using genetic algorithm, we propose a method to retrieve the alignment distribution of transiently aligned CO 2 molecules from the high-order harmonic generation (HHG) spectra. The retrieval method is based on the quantitative rescattering (QRS) model where averaged photorecombination transition dipole can be factored out from the parallel (or perpendicular) harmonic spectra after the propagation of the harmonic fields in the gas medium. We examine how the retrieved alignment distributions are affected by uncertainty in alignment dependent ionization probability and on multiple orbital contribution to the HHG. We further confirm that alignment distribution is more accurately retrieved by using the minima in the HHG spectra driven by a long-wavelength laser. In addition, we show that earlier experimental data on the ratios between the perpendicular and the parallel HHG components of aligned CO 2 molecules are in better agreement with the QRS model if the macroscopic propagation and multiple orbital interference are included in the theoretical calculation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Origins of heat transport anisotropy in MoTe2 and other bulk van der Waals materials

Knowledge of how heat flows anisotropically in van der Waals (vdW) materials is crucial for thermal management of emerging 2D materials devices and design of novel anisotropic thermoelectric materials. Despite the importance, anisotropic heat transport in vdW materials is yet to be systematically studied and is often presumably attributed to anisotropic speeds of sound in vdW materials due to soft interlayer bonding relative to covalent in-plane networks of atoms. In this work, we investigate the origins of the anisotropic heat transport in vdW materials, through time-domain thermoreflectance (TDTR) measurements and first-principles calculations of anisotropic thermal conductivity of three different phases of MoTe 2 . MoTe 2 is ideal for the study due to its weak anisotropy in the speeds of sound. We find that even when the speeds of sound are roughly isotropic, the measured thermal conductivity of MoTe 2 along the c-axis is 5–8 times lower than that along the in-plane axes. We derive meaningful characteristic heat capacity, phonon group velocity, and relaxation times from our first principles calculations for selected vdW materials (MoTe 2 , BP, h-BN, and MoS 2 ), to assess the contributions of these factors to the anisotropic heat transport. Interestingly, we find that the main contributor to the heat transport anisotropy in vdW materials is anisotropy in heat capacity of the dominant heat-carrying phonon modes in different directions, which originates from anisotropic optical phonon dispersion and disparity in the frequency of heat-carrying phonons in different directions. Further, the discrepancy in frequency of the heat-carrying phonons also leads to ~2 times larger average relaxation times in the cross-plane direction, and partially explains the apparent dependence of the anisotropic heat transport on the anisotropic speeds of sound. This work provides insight into understanding of the anisotropic heat transport in vdW materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Combined Machine Learning and Molecular Dynamics Reveal Two States of Hydration of a Single Functional Group of Cationic Polymeric Brushes

The state of hydration of a macromolecular system regulates a plethora of different properties of such a system. In this article, we develop a novel machine learning (ML) approach, based on the unsupervised clustering algorithm, for probing the hydration behavior of the {N(CH 3 ) 3 } + functional group of the PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] polyelectrolyte (PE) brush system. The PE brushes and the brush-supported water molecules and counterions (chloride ions) are first described using all-atom molecular dynamics (MD) simulations. The simulation data is subsequently used in our ML framework to identify that (1) the {N(CH 3 ) 3 } + functional groups of the PMETAC brushes have two distinct hydration states with one state (state 1) being characterized by less structured water molecules and the other state (state 2) being characterized by more structured water molecules and (2) an enhancement in the brush grafting density leads to the progressive dissapparenace of state 2. An increase in the grafting density increases the number of chloride counterions in a given volume around the {N(CH 3 ) 3 } + functional group and increases the number of shared water molecules between the {N(CH 3 ) 3 } + and Cl - . The chloride counterions are associated with a hydration layer with much less structured water molecules. Therefore, with an increase in the grafting density, an increase in the percentage of shared water molecules leads to the prevalence of the hydration state [of the {N(CH 3 ) 3 } + moiety] with less structured water molecules. Finally, we explain how the present findings are commensurate with two key previous related results, namely a significantly large chloride ion mobility inside the PMETAC brush layer and the {N(CH 3 ) 3 } + -Cl - average distance remaining independent of the PMETAC brush grafting density. Furthermore, we anticipate that the combined ML-MD-simulation approach proposed in this study can be adapted to probe other soft matter systems to reveal new insights of the underlying mechanisms of emergent phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Assessment of BPM options for the EIC Beam Accumulator Ring

The electron injection system for the Electron-Ion Collider (EIC) at BNL is designed to provide a beam of polarized electrons, which is crucial for studying the structure of protons and atomic nuclei. The Beam Accumulator Ring (BAR) is a part of the injection chain between the 750 MeV linear accelerator and the Rapid Cycling Synchrotron (RCS), which accelerates the beam up to full energy (5–18 GeV). The functional role of the BAR is to accumulate the charge injected from the linear accelerator in order to achieve the high intensity of the polarized electron beam required by the specifications for injection into the RCS. This objective will be realized through the sequential injection of bunches with a charge of 1.1 nC at a repetition rate of 30 Hz. Once the charge of a single bunch reaches 28 nC, the beam will be extracted from the BAR and injected into the RCS. The beam instrumentation needs to provide reliable measurements with the required accuracy over the dynamic range from 0.1 nC (one tenth of a typical injected bunch charge) to 32 nC – the maximum expected accumulated current. There are 10 Beam Position Monitors (BPMs) in the ring, 1 extra button-electrode assembly for the RF system, and 8 BPMs in the beam transport lines. The BPM locations are shown in Fig. 1, marked by blue rectangles. The ring BPMs will be capable of both average orbit and turn-by-turn measurement modes. Accurate measurement of the beam position with a large horizontal offset requires polynomial correction of the BPM nonlinearity.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Temperature controlled transformations of giant unilamellar vesicles of amphiphilic triblock copolymers synthesized via microfluidic mixing

We report on a simple approach for synthesis of temperature-responsive giant unilamellar vesicles (GUVs) from poly(N-vinylcaprolactam) 15 -block-poly(dimethylsiloxane) 65 -block-poly(N-vinylcaprolactam) 15 (PVCL 15 -PDMS65-PVCL 15 ) triblock copolymer and non-temperature responsive small and giant vesicles from novel poly(N-vinylpyrrolidone)-block-poly(dimethylsiloxane)-block-poly(N-vinylpyrrolidone) (PVPON 15 -PDMS 65 -PVPON 15 and PVPON 6- PDMS 30- PVPON 6 ) triblock copolymers using microfluidic mixing at 25 °C. We show that temperature-responsive PVCL 15 -PDMS 65- PVCL 15 GUVs with the average diameter of 1.4 ± 0.2 µm while being stable at room temperature for at least 14 days, transformed irreversibly into small vesicles of 168 ± 40 nm after incubation of their aqueous solution at 42 °C for 24 h. We hypothesized that this transformation is induced by local compressive stresses of the vesicle membrane due to the collapse of PVCL blocks above the copolymer lower critical solution temperature (LCST) leading to the decrease of the vesicle membrane thickness. Consequently, we found that the temperature-induced size transformation of the PVCL-based GUVs at 42°C can be suppressed by substituting PVCL with its hydrophilic homologue PVPON, or by suppressing the PVCL's LCST behavior through hydrogen-bonding with tannic acid molecules. In the former case, novel PVPON n -PDMS m -PVPON n triblock copolymers (n = 15, m = 65 and n = 6, m = 30) assemble into vesicles stable from 25 °C to 55 °C as confirmed by optical and electron microscopy, dynamic light scattering (DLS) and small-angle neutron scattering (SANS). In the latter, hydrogen bonding interactions of PVCL with the polyphenol tannic acid (TA) at room temperature resulted in stable PVCL-based GUVs at 42°C as confirmed by optical, electron, and atomic force microscopies. We also found that physical crosslinking of the PVCL corona through hydrogen bonding with TA in PVCL 15- PDMS 65- PVCL 15 GUVs will delay their low pH-induced degradation at 37°C by 48 h compared to non-modified GUVs. Our findings open opportunities for the development of temperature-regulated stable micro-vehicles that would change their structural characteristics in the physiologically relevant temperature range from 25 to 42°C and can be utilized for cell mimicking studies. The developed GUVs also have potential in theranostic drug delivery as substitutes for polymer microcapsules and lipid microbubbles as well as for stimuli-triggered sensing, protection, and rapid response in an aqueous environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

High-resolution infrared spectroscopy of jet cooled CH 2 Br radicals: The symmetric CH stretch manifold and absence of nuclear spin cooling

Direct laser absorption of a slit supersonic discharge expansion provides the first high-resolution spectroscopic results on the symmetric CH stretch excitation (ν 1 ) of the bromomethyl (CH 2 Br) radical in the ground electronic state. Narrowband (<1 MHz) mid-infrared radiation is produced by difference-frequency generation of two visible laser beams, with the open shell halohydrocarbon radical generated by electron dissociative attachment of CH 2 Br 2 in a discharge and rapidly cooled to T rot = 18 ± 1 K in the subsequent slit-jet supersonic expansion. A rovibrational structure in the radical spectrum is fully resolved, as well as additional splittings due to spin-rotation effects and 79 Br/ 81 Br isotopologues in natural abundance. Spectroscopic constants and band origins are determined by fitting the transition frequencies to a non-rigid Watson Hamiltonian, yielding results consistent with a vibrationally averaged planar radical and an unpaired electron in the out-of-plane p π orbital. Additionally, extensive satellite band structure from a vibrational hot band is observed and analyzed. The hot band data is compared to CFOUR/VPT2 (CCSD(T)cc-pVQZ) ab initio anharmonic predictions of the vibration rotation alpha matrix, which permits unambiguous assignment to CH 2 symmetric-stretch excitation built on the singly excited CH 2 out-of-plane bending mode (ν 1 + ν 4 ← ν 4 ). Longitudinal cooling of the Doppler width in the slit-jet expansion geometry also reveals partially resolved hyperfine structure on transitions out of the lowest angular momentum states in excellent agreement with predictions based on microwave studies. High level ab initio MOLPRO calculations [CCSD(T)-f12b/VnZ-f12 (n = 3, 4, CBS)] are also performed with explicitly correlated f12 electron methods for the out-of-plane CH 2 bending mode over the halogen series CH 2 X (X = F, Cl, Br, I), which clearly reveals a non-planar geometry for X = F (with a ΔE ≈ 0.3 kcal/mol barrier) and yet planar equilibrium geometries for X = Cl, Br, and I. Finally, a detailed Boltzmann analysis of the transition intensities provides support for negligible collisional equilibration of the entangled H atom nuclear spin states on the few hundred microsecond time scale and high collision densities of a slit supersonic expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Au 76 (SC 6 H 4 - p -CH 3 ) 42 Square Quantum Platelet: One-Dimensional Growth of Quantum Rods Turns 90 Degrees

The growth pattern of atomically precise nanoclusters (NCs) is of fundamental interest, and the structural effect on their photoluminescence (PL) is important owing to their PL in the second near-infrared (NIR-II) range (900–1700 nm optoelectronically or 1000–1700 nm biologically) that holds great promise for optoelectronic and biomedical applications. Because of the small energy gaps required for NIR-II emission, the PL performance of NIR-II luminophores is largely limited by nonradiative processes. In this work, we discovered a Au 76 (p-MBT) 42 (Au 76 ) (p-MBTH = p-methylbenzenethiol) nanocluster featuring a face-centered cubic (fcc) core in a square shape (edge length: 1 nm). This square quantum platelet can be viewed as a side-facet (010) growth of the Au 52 (p-MBT) 32 (Au 52 ) rod, as opposed to the (001) facet growth. We found that Au 76 exhibits bright emission centered at 970 nm with a PL quantum yield (PLQY) of 30% in solution under ambient conditions, which can be further enhanced to 40% when the solution is deaerated. X-ray crystallography analysis coupled with time-resolved spectroscopy revealed that the nearly doubled PLQY compared to Au 52 (18.3%) was resulted from shorter Au–Au bond lengths in Au 76 (average 2.835 Å) than that in Au 52 (3.04 Å). This work provides important insights into the design of highly luminescent NCs, which are promising for photovoltaics, photocatalysis, and optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗
DOE OSTI2022

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗