Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “densities”

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 55 records · Page 3

Effect of resistivity on poloidal asymmetries in electron density in the edge and scrape off layer in 3D full-F gyro-fluid simulations with FELTOR

In this contribution, we analyze the poloidal asymmetries in electron density in the edge and scrape-off layer (SOL) in a COMPASS-sized, diverted tokamak with the 3D full-F, isothermal, electromagnetic, and gyro-fluid model FELTOR. The study is performed for different simulations that span over 2 orders of magnitude in resistivity. The poloidal asymmetries of density are evaluated in the edge and the SOL relative to the outer midplane (OMP), where the highest densities are usually found. In the closed magnetic surfaces, the relative poloidal asymmetry with respect to the OMP is not larger than 20% independently of the plasma resistivity. For the open field lines in the SOL, the relative density asymmetry can range from 55% for the highest resistivity to around 40% for the lowest. The lowest densities inside the separatrix are found between the inner midplane and the top of the magnetic configuration, away from the X-point. In the SOL it is usually close to the X-point in the high field side. The observations in the closed field lines are consistent with ballooning transport but in the SOL oppose experimental evidence for high-density plasmas. In conclusion, this indicates the necessity for more complex physics to reproduce the experimental observations in the SOL for higher density plasmas, such as neutral-plasma dynamics and realistic divertor conditions.

3D full-f gyrofluid↗

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hot Hydrogen Exposure of U x Zr 1-x C y Nuclear Fuel: The Influence of Composition and Density

Refractory carbide nuclear fuel has been one of the most promising fuel candidates for space nuclear propulsion due to its high melting point, temperature stability, and compatibility in a hot hydrogen environment. In this study, U x Zr 1-x C y fuel was produced by means of a carbothermic reduction process in different UC compositions including 5,10, 20, and 30 at.% UC in the fuel compound. The powder feedstock was consolidated via direct current sintering with densities up to 97% of the theoretical density. The samples with different U x Zr 1-x C y compositions were exposed to hot hydrogen at 2600 K for a cumulative time of 300 min. The samples were characterized by SEM, XRD, density, and measured for mass losses. The high-density samples displayed improved performance in hot hydrogen by minimizing porous sites and reducing areas of direct contact with hydrogen gas, leading to reduced mass losses. Variations in sample density proved to induce large changes in mass loss rates, increasing them up to 90%. The compositions with higher UC content reported the largest mass losses in the study. The loss of uranium occurred primarily at the surfaces exposed to the hot hydrogen where changes in the lattice constant confirmed losses exceeding 50% of the initial UC content in higher compositions, specifically to U 0.3 Zr 0.7 C y . XRD analyses revealed the presence of UH 3 in U 0.3 Zr 0.7 C y suggesting that metallic uranium formed inside the sample as a product of carbon losses. High-density U x Zr x-1 C y fuel with UC concentrations at or below 20 at.% UC exhibited stability and negligible density changes in a high temperature hydrogen environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Inhomogeneous high temperature melting and decoupling of charge density waves in spin-triplet superconductor UTe 2

Charge, spin and Cooper-pair density waves have now been widely detected in exotic superconductors. Understanding how these density waves emerge — and become suppressed by external parameters — is a key research direction in condensed matter physics. Here we study the temperature and magnetic-field evolution of charge density waves in the rare spin-triplet superconductor candidate UTe 2 using scanning tunneling microscopy/spectroscopy. We reveal that charge modulations composed of three different wave vectors gradually weaken in a spatially inhomogeneous manner, while persisting to surprisingly high temperatures of 10–12 K. We also reveal an unexpected decoupling of the three-component charge density wave state. Our observations match closely to the temperature scale potentially related to short-range magnetic correlations, providing a possible connection between density waves observed by surface probes and intrinsic bulk features. Importantly, charge density wave modulations become suppressed with magnetic field both below and above superconducting T c in a comparable manner. Our work points towards an intimate connection between hidden magnetic correlations and the origin of the unusual charge density waves in UTe 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge density wave in a band insulator

Charge-density-wave implies a periodic modulation of the charge density. Typically observed in metallic systems, charge-density-wave arise from Fermi surface instabilities, resulting in the total or partial gapping of the Fermi surface. Here, we present experimental evidence for a charge-density wave state emerging in a band insulator which has no Fermi surface. The bulk and surface of our material platform, 𝛂-Bi 4 Br 4 , is gapped over the entire Brillouin zone. Through topographic and spectroscopic imaging at low temperatures, we unveil an unexpected unidirectional charge modulation in 𝛂-Bi 4 Br 4 , breaking the lattice translation symmetry. The charge-density-wave develops at temperatures below 40 K and adds an energy gap atop the existing insulating gap of 𝛂-Bi 4 Br 4 . Furthermore, our transport measurements reveal nonlinear electrical conduction, a phenomenon conventionally associated with the sliding or phason mode of incommensurate charge-density-waves. These unusual observations represent a different type of charge-density-wave

Phase transitions and critical phenomena↗

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗

Medium-density amorphous ice unveils shear rate as a new dimension in water’s phase diagram

Recent experiments revealed a new amorphous ice phase, medium-density amorphous ice (MDA), formed by ball-milling ice I h at 77 K [Rosu-Finsen et al., Science 379, 474–478 (2023)]. MDA has density between that of low-density amorphous (LDA) and high-density amorphous (HDA) ices, adding to the complexity of water’s phase diagram, known for its glass polyamorphism and two-state thermodynamics. The nature of MDA and its relation to other amorphous ices and liquid water remain unsolved. Here, we use molecular simulations under controlled pressure and shear rate at 77 K to produce and investigate MDA. Here. we find that MDA formed at constant shear rate is a steady-state nonequilibrium shear-driven amorphous ice (SDA), that can be produced by shearing ice I h , LDA, or HDA. Our results suggest that MDA could be obtained by ball-milling water glasses without crystallization interference. Increasing the shear rate at ambient pressure produces SDAs with densities ranging from LDA to HDA, revealing shear rate as a new thermodynamic variable in the nonequilibrium phase diagram of water. Indeed, shearing provides access to amorphous states inaccessible by controlling pressure and temperature alone. SDAs produced with shearing rates as high as 10 6 s −1 sample the same region of the potential energy landscape than hyperquenched glasses with identical density, pressure, and temperature. Intriguingly, SDAs obtained by shearing at ~10 8 s −1 have density, enthalpy, and structure indistinguishable from those of water “instantaneously” quenched from room temperature to 77 K over 10 ps, making them good approximants for the “true glass” of ambient liquid water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of the density field around a supersonic conical projectile using quantitative schlieren

Quantitative schlieren imaging is a flow measurement technique that is capable of measuring density fields throughout refractive flowfields. The technique was applied here to measure the density field surrounding supersonic conical projectiles in free flight. Shock waves attached to a supersonic conical projectile offer a simple geometry with the well-established Taylor–Maccoll analytical flow solution to which these experimental measurements were compared. The schlieren images recorded a projection of the index of refraction field surrounding the 10° half-angle cones which was converted to density first through an Abel inversion and then the Gladstone–Dale law. Three Abel inversion methods—two-point, three-point, and arbitrary ray axisymmetric projection (ARAP)—were applied to deconvolute the three-dimensional flow within the constrained axisymmetric flow field. The resulting reconstructed density profiles were compared to the Taylor–Maccoll solution, parameterized by cone geometry and Mach number. The experimental density fields demonstrated strong agreement with the theoretical profiles. Experimental consistency was confirmed across various projectile speeds, demonstrating quantitative schlieren’s capability to accurately reconstruct the density of the flow field, even within the resolution constraints imposed by high-speed imaging. In conclusion, an assessment of experimental uncertainties in the density reconstruction was performed.

Abel Inversion↗

High-field side scrape-off layer density profile measurements and implications for high-field side LHCD coupling in DIII-D

The high-field side (HFS) scrape-off layer (SOL) is an often under-diagnosed region of tokamak plasmas. Situated in a region with favorable curvature, the HFS SOL has minimal turbulence-induced radial transport compared to the low-field side. Using profile reflectometry, the HFS SOL density profile is measured with high temporal resolution for a wide range of DIII-D plasma discharges for the first time. The magnetic configuration, particularly the location of the secondary separatrix (SS) largely determines the HFS SOL density profile. Additionally, density perturbations induced by edge-localized modes (ELMs) are observed in the HFS SOL, and the impact of ELMs on the SOL density profile is characterized for a range of magnetic configurations. It is found that ELM-induced HFS SOL perturbations are generally localized to within the SS and ELM-induced changes to the HFS wall density are minimized in near double-null configurations. The characterization of the HFS SOL density profile on DIII-D is crucial for the effective coupling of the planned HFS lower hybrid current drive (LHCD) launcher. Experimental HFS SOL measurements are used as input to full-wave simulations of LHCD coupling. Furthermore, this work simulates the effect of the magnetic configuration on LHCD coupling. It also simulates LHCD coupling in ELMs and high q min discharges. During such discharges, the LHCD coupling is predicted to be resilient to ELM-induced SOL density perturbations.

high-field side↗

Experimental investigation of PDI bifurcation of lower hybrid waves during electron density ramp-up in EAST

Abstract The effect of parametric decay instability (PDI) on the current drive efficiency of 4.6 GHz lower hybrid (LH) waves in EAST is investigated experimentally, showing the PDI channel bifurcation of LH waves for the first time in EAST. First, experiments with three platforms of LH power were performed, achieving the LH power required for the PDI occurrence. Second, PDI bifurcation experiments were further carried out by ramping up the plasma electron density. The loop voltage increases with an increase in density, implying a decrease in the current driven by the LH wave. PDI bifurcation during electron density ramp-up was studied by analyzing the parallel refractive index ( n ∥ ) and the frequency spectrum broadening, which is measured with a radio frequency magnetic probe array recently installed close to the LH antenna. It is observed for the first time that they both first increase with density, then there is not much variation and a clear sideband in the frequency spectrum is also observed when the density is up to 4 × 10 19 m − 3 , suggesting a change in the PDI channel. Calculation of the mode growth rate driven by PDI shows that when the edge electron density is up to 1.9 × 10 18 m − 3 , the growth rate of the ion cyclotron quasi-mode (ICQM) will exceed that of the ion sound quasi-mode (ISQM), quantitively explaining that with an increase in density, the PDI channel partly transits from the ISQM to the ICQM channel. Studies provide a possible way to reduce the power deposition in the edge region and improve drive capability by means of mitigating PDI behavior.

Physics↗

Impact of ionization and transport on pedestal density structure in DIII-D and Alcator C-Mod

Abstract This paper investigates the role of ionization on the pedestal structure using both measurements and modeling for H-mode plasmas on DIII-D and Alcator C-Mod to enhance our ability to predict pedestal behavior in future pilot plants. The impact of the neutral penetration depth on the pedestal density is investigated using dimensionally matching hydrogen and deuterium DIII-D H-mode discharges at low and high electron density. The DIII-D Lyman- α diagnostic measurements show that hydrogen neutrals penetrate deeper inside the plasma on both the high field and low field side, while the pedestal electron density structure is similar for both isotopes. However, as the opaqueness increases we observe that the pedestal density gradient becomes stiff, similar to prior observations on DIII-D and C-Mod (Mordijck 2020 Nuclear Fusion 60 082006). In addition, these results also confirm prior measured and modeled poloidal asymmetries in neutral densities, indicating that to make transport predictions, 2D neutral modeling is necessary. The first direct validation of SOLPS-ITER for the measured brightness, emissivity and neutral densities for three different confinement regimes on C-Mod is introduced. The SOLPS-ITER model shows good agreement, within the constrains of the model for all regimes. In addition, a comparison of SOLPS-ITER modeling for DIII-D and C-Mod shows that as opaqueness increases, the role of divertor fueling and thus poloidal asymmetries in the neutral density profiles decreases. Based on these experimental and modeling results we estimate the size of a potential particle pinch using typical values for the diffusion coefficient for both DIII-D and C-Mod H-mode discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-linear dependence of ion heat flux on plasma density at the L–H transition of JET NBI-heated deuterium–tritium plasmas

Abstract Recent JET D–T campaigns opened the possibility of unique isotope studies to investigate the L–H transition physics in view of reactor plasmas and to study the origin of the observed power threshold minimum. In the present paper, we characterise L–H transitions in the low and high-density branches of JET NBI-heated D–T plasmas. As discussed in the paper, L–H transition has been hypothesised to be determined by the transport power losses of plasma ions, i.e. the so-called ion heat flux ( Q i ). We present the first power balance analysis of JET NBI-heated D–T plasmas to evaluate the ion heat flux at the transition. Due to the experimental setting being similar to previous JET D experiments, we also directly compare the results, discussing the isotope effect and similarities between datasets. First, we find an isotope effect between D and D–T Q i , with a lower Q i in D–T plasmas. We confirm that the ion heat flux deviates from density linearity compared to the linear trend observed in wave-heated D plasmas of other tokamaks. The deviation we observe in NBI-heated L–H transitions happens at an isotope-dependent density. Plasma edge rotation correlates with Q i deviation from density linearity in the low-density branch. However, further investigations would be required to assess the role of rotation on Q i and the power threshold minimum at JET. At low plasma density, NBI power dominates Q i , while increasing the density makes the equipartition power dominant. We finally compare our results with hypotheses proposed from evidence in other tokamaks to present a complete overview of ion heat flux analyses in D and D–T NBI-heated plasmas at JET.

Vincenzi, P. (ORCID:0000000251564354)↗

Exploring the transition from continuous turbulence fluctuations to bursting ELMs in high SOL density regimes

BOUT++ turbulence simulations of the DIII-D reveal that the density profile between the separatrix and pedestal plays a crucial role in the dynamics of edge localized modes (ELMs) and edge plasma turbulent transport. Nonlinear simulations demonstrate that small ELMs in the DIII-D hybrid scenario under high SOL density conditions are predominantly driven by local ballooning modes near the separatrix, stabilizing global instabilities while enhancing localized pressure fluctuations. A key control parameters for ELM dynamics is the separatrix-to-pedestal density ratio, n e,sep /n e,ped . A high ratio indicates a shallow gradient, favoring small ELMs, while a lower ratio signals a steep gradient, which increases the likelihood of large ELMs. Comprehensive parameter scans, including n e,sep /n e,ped , density gradient profiles near the separatrix, and resistivity, reveal the critical role of these parameters in shaping transitions between turbulence-driven transport and ELM bursting. The scans demonstrate that in high SOL density regimes, small ELMs can result from either global resistive MHD instabilities or local ballooning modes near the separatrix, depending on the steepness of the separatrix density gradient. These findings also highlight the transition from continuous turbulence to small ELMs. The post-crash peak in pressure fluctuations, δP rms serves as a critical metric for identifying transition from continuous turbulence fluctuations to ELM bursting. Larger δP rms values correlate with ELM bursts driven by local or global instabilities, whereas smaller values indicate turbulence-dominated transport. Drift-Alfvén and resistive ballooning turbulence enhance the entrainment of fluctuations from the pedestal to the SOL, contributing to the complex interplay of dynamics in this regime. These findings emphasize the importance of separatrix density shaping and pedestal gradient control for optimizing ELM behavior in ITER and future fusion devices.

Li, Nami [Lawrence Livermore National Laboratory (↗

CryoTEN: efficiently enhancing cryo-EM density maps using transformers

Abstract Motivation Cryogenic electron microscopy (cryo-EM) is a core experimental technique used to determine the structure of macromolecules such as proteins. However, the effectiveness of cryo-EM is often hindered by the noise and missing density values in cryo-EM density maps caused by experimental conditions such as low contrast and conformational heterogeneity. Although various global and local map-sharpening techniques are widely employed to improve cryo-EM density maps, it is still challenging to efficiently improve their quality for building better protein structures from them. Results In this study, we introduce CryoTEN—a 3D UNETR++ style transformer to improve cryo-EM maps effectively. CryoTEN is trained using a diverse set of 1295 cryo-EM maps as inputs and their corresponding simulated maps generated from known protein structures as targets. An independent test set containing 150 maps is used to evaluate CryoTEN, and the results demonstrate that it can robustly enhance the quality of cryo-EM density maps. In addition, automatic de novo protein structure modeling shows that protein structures built from the density maps processed by CryoTEN have substantially better quality than those built from the original maps. Compared to the existing state-of-the-art deep learning methods for enhancing cryo-EM density maps, CryoTEN ranks second in improving the quality of density maps, while running >10 times faster and requiring much less GPU memory than them. Availability and implementation The source code and data are freely available at https://github.com/jianlin-cheng/cryoten.

Biochemistry & Molecular Biology↗

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗

A measurement of $H_0$ from DESI DR1 using energy densities

We present a new measurement of the Hubble constant, independent of standard rulers and robust to pre-recombination modifications such as Early Dark Energy (EDE), obtained by calibrating the total energy density of the Universe. We start using the present-day photon density as an anchor, and use the baryon-to-photon ratio from Big Bang Nucleosynthesis based measurements and the baryon-to-matter ratio from the baryons' imprint on galaxy clustering to translate to a physical matter density at present day. We then compare this to measurements of the ratio of the matter density to the critical density ($Ω_{\mathrm{m}}$), calculated using the relative positions of the baryon acoustic oscillations, to measure the critical density of the universe and hence $H_0$. The important measurements of the evolution of the energy density all happen at low redshift, so we consider this a low-redshift measurement. We validate our method both on a suite of $N$-body mocks and on noiseless theory vectors generated across a wide range of Hubble parameters in both $Λ$CDM and EDE cosmologies. Using DESI DR1 data combined with the angular CMB acoustic scale and the latest BBN constraints, we find $H_0 = 69.0 \pm 2.5$ km s$^{-1}$ Mpc$^{-1}$, consistent with existing early and late-time determinations of the Hubble constant. We consider the impact of non-standard dark energy evolution on our measurement. Future data, including that from further iterations of DESI and from Euclid, will add to these results providing a powerful test of the Hubble tension.

Krolewski, Alex [U. Waterloo (main); Waterloo U., ↗

Intermolecular Interactions in Direct Air Capture Materials: Insights from Charge Density Analysis

Direct air capture (DAC) materials enable the removal of CO 2 from the atmosphere, but improving their efficiency requires a detailed understanding of the intermolecular interactions that govern CO 2 sorption and release. Here, we present an experimental electron density study of methylglyoxal-bis(iminoguanidine) (MGBIG), a promising DAC material, using high-resolution X-ray and neutron diffraction data combined with quantum crystallographic analysis. This approach bridges theoretical and experimental data by quantifying electron density distributions and revealing how hydrogen bonds stabilize CO 2 -derived carbonate phases and may influence the desorption behavior. We identify distinct hydrogen-bonding environments in two crystalline carbonate phases: P1, a transient kinetic product, and P3, a thermodynamically stable phase. Multipolar refinement and electrostatic potential and multipole moment calculations precisely map electron density distributions, revealing key hydrogen bonds involved in CO 2 capture. Topological analysis of electron density highlights a cooperative hydrogen-bonding network in the thermodynamically favored P3 phase, where enhanced electron density delocalization and water-mediated interactions contribute to a more stable lattice. Energetic analyses confirm that stronger hydrogen bonding networks enhance the stability of P3 with a binding energy of −607.0 kJ/mol and greater lattice stability (−847.3 kJ/mol) compared to P1 (−302.5 and −571.0 kJ/mol, respectively). Electrostatic potential maps further illustrate polarization patterns that may influence the stability of the binding of CO 2 and release conditions. These findings establish a direct experimental framework for linking electron density distributions to intermolecular interactions in DAC materials, providing a rational design strategy for optimizing sorbents with improved CO 2 capture efficiency and reduced energy demands.

Electron density↗

Application of an Empirical Density Law via Python for Aqueous Plutonium Chloride Systems for MCNP6

Criticality safety models for aqueous plutonium chloride systems often contain a significant bias due to assumptions in material compositions. These systems are currently modeled as a fictitious metal-water mixture because little is known about the true solution density. Furthermore, no predictive density tools or capabilities for modeling aqueous plutonium chloride systems are approved for use at Los Alamos National Laboratory. Recent density measurements of this ternary system (PuCl 3 -HCl-H 2 O) have allowed for the development of a more realistic density law, which is applied in this work via an empirical method based in Python. This tool, entitled PuCS (Plutonium Chloride Solution tool) may be used to determine solution density and composition based on the plutonium content, acid content, and temperature for MCNP6 inputs. PuCS has been found to predict density within 2% of experimental data. In conclusion, MCNP6 calculations have found that crediting minimal amounts of free acid (0.5 M) may correspond to a ~12% decrease in peak reactivity in comparison to current modeling methods.

42 ENGINEERING↗