Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “k mean”

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 523 records · Page 29

A Python Tool for Aqueous Plutonium Nitrate Density Law Input Preprocessing in MCNP6

Here, this work develops a predictive density tool in Python, named Plutonium Nitrate Solutions (PuNS), to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer method and an empirical method were implemented into the PuNS tool to generate atom densities for use in MCNP6 material cards. These material cards are directly prepared into an MCNP6 input text file and are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, and plutonium isotope weight percentages. The PuNS tool is validated and verified against the International Criticality Safety Benchmark Evaluation Project Handbook experiments and is observed to predict densities within a root mean square error of 0.89% for the Pitzer method and 1.82% for the empirical method. These errors in density lead to up to 1569 pcm difference in MCNP6 calculated k eff for the Pitzer method and up to a 1751 pcm difference for the empirical method when compared to experimental benchmarks. Simultaneous work is also being performed at Los Alamos National Laboratory and the University of New Mexico to create a similar tool for plutonium chloride solutions, named Plutonium Chloride Solution, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also facilitate an improved understanding of solution systems and a potential relaxation in the conservatism of current aqueous plutonium processing criticality safety limits.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Intrinsic electric quadrupole moment of the K π = 8 - isomeric state in Hf 178

The lifetime of the 9 - state in the rotational band based on the 4.0 s, K π = 8 - , isomeric state ( 178 Hf m 1 ) from the decay of the 31-yr isomer ( 178 Hf m 2 ) was determined to be 99(2) ps by means of the fast-timing technique using two LaBr 3 (Ce) scintillators. The δ (E2/M1) mixing ratios of the ΔI = 1 γ rays depopulating levels in this band were deduced from γ-γ angular correlations by using a 178 Hf m 2 radioactive source located at the center of the Gammasphere HPGe detector array. The new results, together with previous spectroscopic information, provide a different way to extract the intrinsic quadrupole moment of Q 0 = 6.45 (14) eb for the 178 Hf m 1 band. A possible explanation for the reduction of the 178 Hf m 1 nuclear charge radius is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Including Chi-Nu 235 U PFNS Experimental Data into an ENDF/B-VIII.1 Release Candidate Evaluation

This report documents an evaluation of 235 U prompt fission neutron spectra (PFNS) that is a release candidate for the upcoming U.S. nuclear data library, ENDF/B-VIII.1. This evaluation differs from its predecessor, ENDF/B-VIII.0, mainly by the inclusion of 235 U PFNS measured by the Chi-Nu team of LANL and LLNL. This data set is the first one that covers the 235 U PFNS for continuous incident-neutron energies of 1⁻20 MeV and outgoing-neutron energies from 10 keV⁻10 MeV with high precision. Previous data sets were either measured in a limited energy range or with less precision. Hence, these new Chi-Nu data provide decisive information for the evaluation. The resulting evaluated data correspond well to the new experimental PFNS. The evaluated PFNS also produce average mean energies and 239 Pu/ 235 U PFNS in agreement with associated Chi-Nu data. If one uses the new evaluated data to predict the neutron multiplication factor, k eff , of selected ICSBEP critical assemblies, the differences of simulated values compared to those using ENDF/B-VIII.0 is modest (less than 55 pcm). This difference in k eff can be easily accommodated by changes in the 235 U average prompt fission neutron multiplicity that is currently being re-evaluated. In addition to that, the new PFNS predict on average 235 U LLNL pulsed-sphere neutron-leakage spectra better than ENDF/B-VIIII.0 PFNS.

235U↗

Influence of hydrogenation on the vibrational density of states of magnetocaloric LaFe 11.4 Si 1.6 H 1.6

We report on the impact of magnetoelastic coupling on the magnetocaloric properties of LaFe 11.4 Si 1.6 H 1.6 in terms of the vibrational (phonon) density of states (VDOS), which we determined with 57 Fe nuclear resonant inelastic x-ray scattering (NRIXS) measurements and with density functional theory (DFT) based first-principles calculations in the ferromagnetic (FM) low-temperature and paramagnetic (PM) high-temperature phase. In experiments and calculations, we observe pronounced differences in the shape of the Fe-partial VDOS between nonhydrogenated and hydrogenated samples. Furthermore, this shows that hydrogen not only shifts the temperature of the first-order phase transition, but also affects the elastic response of the Fe subsystem significantly. In turn, the anomalous redshift of the Fe VDOS, observed by going to the low-volume PM phase, survives hydrogenation. As a consequence, the change in the Fe-specific vibrational entropy ΔS lat across the phase transition has the same sign as the magnetic and electronic contribution. DFT calculations show that the same mechanism, which is a consequence of the itinerant electron metamagnetism associated with the Fe subsystem, is effective in both the hydrogenated and the hydrogen-free compounds. Although reduced by 50% as compared to the hydrogen-free system, the measured change ΔS lat of 3.2 ± 1.9)J/ kg K across the FM-to-PM transition contributes with ~35 % significantly and cooperatively to the total isothermal entropy change ΔS iso . Hydrogenation is observed to induce an overall blueshift of the Fe VDOS with respect to the H-free compound; this effect, together with the enhanced Debye temperature observed, is a fingerprint of the hardening of the Fe sublattice by hydrogen incorporation. In addition, the mean Debye velocity of sound of LaFe 11.4 Si 1.6 H 1.6 was determined from the NRIXS and the DFT data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High temperature oxidation of corrosion resistant alloys from machine learning

Parabolic rate constants, k p , were collected from published reports and calculated from corrosion product data (sample mass gain or corrosion product thickness) and tabulated for 75 alloys exposed to temperatures between ~800 and 2000 K (~500–1700 °C; 900–3000°F). Data were collected for environments including lab air, ambient and supercritical carbon dioxide, supercritical water, and steam. Materials studied include low- and high-Cr ferritic and austenitic steels, nickel superalloys, and aluminide materials. A combination of Arrhenius analysis, simple linear regression, supervised and unsupervised machine learning methods were used to investigate the relations between composition and oxidation kinetics. The supervised machine learning techniques produced the lowest mean standard errors. The most significant elements controlling oxidation kinetics were Ni, Cr, Al, and Fe, with Mo and Co composition also found to be significant features. The activation energies produced from the machine learning analysis were in the correct distributions for the diffusion constants for the oxide scales expected to dominate in each class.

Materials Science↗

Thermodynamics of diamond formation from hydrocarbon mixtures in planets

Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the melting line of diamond, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where it is thermodynamically possible for diamond to form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, here we show a depletion zone at pressures above 200 GPa and temperatures below 3000 K-3500 K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electromagnetic pinned solitons for space debris detection

Electromagnetic “pinned” solitons in the form of stationary nonlinear waves are studied within the framework of an inertial magneto-hydrodynamic model. These structures, that can arise when a charged source moves in a magnetized plasma, have a velocity that is equal to the source velocity and, hence, appear as “pinned” structures that envelope the source. We investigate the excitation of such solitons in the Low Earth Orbit region due to the passage of charged orbital debris objects. The spatial size of these electromagnetic solitons, typically of the order of a few ion skin depths, can be very large in this region. Such solitons can be detected using a variety of ground- or orbit-based radio sounding techniques and may provide a convenient additional means of tracking small sized orbital debris objects that are difficult to spot optically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Cauchy problem for Boltzmann equation modeling a polyatomic gas

In the present article, we consider the Boltzmann equation that models a polyatomic gas by introducing one additional continuous variable, referred to as microscopic internal energy. We establish existence and uniqueness theory in the space homogeneous setting for the full non-linear case, under an extended Grad-type assumption on transition probability rates, which comprises hard potentials for both the relative speed and internal energy with the rate in the interval 0,2, multiplied by an integrable angular part and integrable partition functions. The Cauchy problem is resolved by means of an abstract ordinary differential equation (ODE) theory in Banach spaces for the initial data with finite and strictly positive gas mass and energy, finite momentum, and additionally finite K* polynomial moment, with K* depending on the rate of the transition probability and the structure of a polyatomic molecule or its internal degrees of freedom. Moreover, we prove that polynomially and exponentially weighted Banach space norms associated with the solution are both generated and propagated uniformly in time.

Physics↗

Suppression of low temperature magnetic ordering in samarium nanoparticles

The role of finite size effects on magnetic order has been investigated in samarium nanoparticles prepared by physical vapor deposition. In this study, a dense layer composed of distinct nanoparticles with a mean particle diameter of 26 nm was deposited on a diamagnetic substrate. M(T) measurements identify the expected pair of antiferromagnetic ordering temperatures in the bulk Sm precursor, at 113 K and 14 K, where the magnetic unit cell for the lower ordering temperature is 10.36 nm along the c-axis. The high temperature ordering of the hexagonal sites in the Sm nanocrystals is slightly decreased with respect to that of bulk Sm, while the low temperature transition associated with the cubic sites is significantly suppressed. The observed changes are attributed to finite size effects, with ordering suppressed as the particle radius approaches the length of the magnetic unit cell, and surface moments become more prominent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fathoming the anisotropic magnetoelasticity and magnetocaloric effect in GdNi

Intermetallic GdNi adopts a CrB type of crystal structure (space group Cmcm), and it orders ferromagnetically via a second-order phase transition at 70 K, exhibiting unusually strong spontaneous striction along the three independent crystallographic axes in the ferromagnetically ordered state. Here we introduce a microscopic model to describe anisotropic changes of lattice parameters and elastic contribution to magnetocaloric effect of GdNi. In the model, results of density functional theory (DFT) calculations are used as inputs into a Hamiltonian that includes elastic energy of an anisotropic crystal lattice, exchange interactions, and Zeeman effect. The magnetic and elastic Hamiltonians are coupled through an anisotropic Bean-Rodbell model of magnetoelastic interactions. This coupling gives rise to anisotropic changes in the lattice parameters observed experimentally, and the model reveals good to reasonable agreements between the current theoretical results and earlier experimental data, thus validating the model within the limits of assumptions made. We also show that DFT calculations with 4f electrons of Gd treated as core electrons lead to a more adequate estimate of elastic constants of GdNi in comparison with the LDA + U method where 4f electrons are treated as valence electrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Isotopically resolved neutron total cross sections at intermediate energies

The neutron total cross sections σ tot of 16,18 O, 58 64 Ni, 103 Rh, and 112,124 Sn have been measured at the Los Alamos Neutron Science Center from low to intermediate energies (3 ≤ E lab ≤ 450 MeV) by leveraging wave-form-digitizer technology. The σ tot relative differences between isotopes are presented, revealing additional information about the isovector components needed for an accurate optical-model description away from stability. Digitizer-enabled σ tot -measurement techniques are discussed and a series of uncertainty-quantified dispersive optical model (DOM) analyses using these new data is presented, validating the use of the DOM for modeling light systems ( 16,18 O) and systems with open neutron shells ( 58,64 Ni and 112,124 Sn). Here, the valence-nucleon spectroscopic factors extracted for each isotope reaffirm the usefulness of high-energy proton reaction cross sections for characterizing depletion from the mean-field expectation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure of Low-spin States in \(^{45}\)Sc Studied via Coulomb Excitation

The electromagnetic structure of 45 Sc at low excitation energy was investigated via low-energy Coulomb excitation at the Heavy Ion Laboratory (HIL) of the University of Warsaw and at the Inter-University Accelerator Centre (IUAC) in New Delhi. A set of reduced E2, E3, and M1 matrix elements was extracted from the collected data using the GOSIA code. The reduced transition probability B(E2; 11/2 – → 7/2 – ) has been determined, allowing us to deduce the lifetime of the 11/2 – state at 1237 keV. In addition, the upper limit on the reduced transition probability B(E3; 7/2 – → 5/2 + ) has been determined for the first time. Finally, new large-scale shellmodel and beyond-mean-field calculations were performed to interpret the structure of this nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application of a Density Law via Python for Aqueous Plutonium Nitrate

A predictive density tool has been developed in Python to reduce bias and uncertainty in nuclear criticality safety calculations for plutonium nitrate systems. The Pitzer Method and an empirical method were implemented into the tool, allowing for plutonium nitrate density calculations. Additionally, the Python tool can generate atom densities for a MCNP6.2 material card using the density from the selected method and directly the densities into a prepared MCNP6 input text file. The material card and density are calculated based on customizable user inputs of plutonium content, nitric acid content, temperature, plutonium isotope weight percentages and impurity concentrations. The Python tool has been validated and verified against the International Handbook of Evaluated Criticality Safety Benchmark Experiments to predict densities within a root mean square error of 1.0% for the Pitzer method and 1.8% for the empirical method. These errors in density were shown to lead to a ±0.5% error in MCNP6.2 calculated k effective for the Pitzer method and a ±1.7% error for the Empirical method. Simultaneous work is also being done at the University of New Mexico and Los Alamos National Laboratory to create a similar tool for plutonium chloride solutions, which aims to provide the accreditation of the chlorine absorption. These capabilities will not only provide more accurate models but also lead the way towards a better understanding of solution systems and potential relaxation in the conservatism of the current aqueous plutonium processing limits.

97 MATHEMATICS AND COMPUTING↗

Probing FeSi, a d -electron topological Kondo insulator candidate, with magnetic field, pressure, and microwaves

Recently, evidence for a conducting surface state (CSS) below 19 K was reported for the correlated d -electron small gap semiconductor FeSi. In the work reported herein, the CSS and the bulk phase of FeSi were probed via electrical resistivity ρ measurements as a function of temperature T , magnetic field B to 60 T, and pressure P to 7.6 GPa, and by means of a magnetic field-modulated microwave spectroscopy (MFMMS) technique. The properties of FeSi were also compared with those of the Kondo insulator SmB 6 to address the question of whether FeSi is a d -electron analogue of an f -electron Kondo insulator and, in addition, a “topological Kondo insulator” (TKI). The overall behavior of the magnetoresistance of FeSi at temperatures above and below the onset temperature T S = 19 K of the CSS is similar to that of SmB 6 . The two energy gaps, inferred from the ρ( T ) data in the semiconducting regime, increase with pressure up to about 7 GPa, followed by a drop which coincides with a sharp suppression of T S . Several studies of ρ( T ) under pressure on SmB 6 reveal behavior similar to that of FeSi in which the two energy gaps vanish at a critical pressure near the pressure at which T S vanishes, although the energy gaps in SmB 6 initially decrease with pressure, whereas in FeSi they increase with pressure. The MFMMS measurements showed a sharp feature at T S ≈ 19 K for FeSi, which could be due to ferromagnetic ordering of the CSS. However, no such feature was observed at T S ≈ 4.5 K for SmB 6 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Comparing Pr 3+ and Nd 3+ for deactivating the Er 3+ : 4 I 13/2 level in lanthanum titanate glass

Erbium lanthanum titanate glasses were prepared by levitation melting for the spectroscopic study of ways to promote the mid-infrared fluorescence. Two series of heavily erbium doped glasses (15 wt%) were prepared with the addition of either Pr 3+ or Nd 3+ in amounts relative to Er 3+ of 0.05, 0.1, and 0.2. Both ions quench the lower Er 3+ laser level with the Pr 3+ doing so more rapidly. Although high co-dopant concentrations result in higher energy transfer, as clearly evidenced in upconversion and downconversion fluorescence measurements, the mid-infrared lifetime also suffers a reduction and, therefore, a balance must be struck in the co-dopant concentration. Lifetime and spectral measurements indicate that, at a fixed relative co-dopant amount, Pr 3+ is more effective than Nd 3+ at removing the bottleneck of the Er 3+ 4 I 13/2 level. Moreover, consideration of the lifetimes alongside the absorption data of the individual ions indicates that despite the large absorption cross-section of Nd 3+ at 808 nm, the concentration needed to yield more absorbed power than utilizing direct 976 nm excitation of Er 3+ results in unfavorable lifetimes of the mid-infrared transition. In the end, Pr 3+ prevails as the superior co-dopant in terms of the effects on fluorescence lifetimes as well as potential laser system design considerations. In a unique self-doping approach, a reducing melt atmosphere of Ar instead of O 2 creates a small fraction of Ti 3+ . In 5Er 2 O 3 -12La 2 O 3 -83TiO 2 glass, the presence of Ti 3+ quenches the 4 I 13/2 emission about 2.6 times more than the 4 I 11/2 when lifetimes are compared to an O 2 melt environment. As an additional means of increasing the mid-infrared emission, the effect of temperature on the mid- and near- infrared lifetimes of a lightly doped lanthanum titanate composition is investigated between 77-300 K. The mid-infrared lifetime increases by ∼30% while the near-infrared lifetime increases by ∼10%, which suggests in addition to co-doping, active cooling of the gain media will further enhance performance.

Materials Science↗

Data-driven estimation of energy consumption for electric bus under real-world driving conditions

Reliable and accurate estimation of an electric bus’s instantaneous energy consumption is critical in evaluating energy impacts of planning and control of electric bus operations. In this study, we developed machine learning-based long short-term memory (LSTM) and artificial neural network (ANN) models to estimate 1 Hz energy consumption of electric buses based on continuous monitoring data of electric buses in Chattanooga, Tennessee, in 2019 and 2020. We propose a data-partitioning algorithm to separate energy charging and discharging modes before applying data-driven estimation models. Here, a K-fold cross-validation-based model selection process was conducted to identify the optimal model structure and input variables in terms of prediction accuracy. The estimation results show the predicted mean absolute percentage error rates of LSTM and ANN models were 3% and 5%, respectively. We compared the proposed models with existing models in the literature based on the same testing data to demonstrate the predictability of our models.

Artificial neural network↗