Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code comparisons”

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

Atomic forces from Dirac-Kohn-Sham equations: Implementation in flexible (APW+lo/LAPW)+LO basis set

The atomic forces formulation based on the Dirac-Kohn-Sham equation and the flexible (APW+lo/LAPW)+LO basis set is presented here. The formulation was implemented in the code FlapwMBPT and allows a user to easily switch between different basis functions of the augmentation type (APW or LAPW) and between different kind of local orbitals. Similar to the work (Phys.Rev.B 91 (2015) 035105), the implementation takes into account small discontinuities of the wave functions, density, and potential at the muffin-tin sphere boundaries. Applications to the materials with strong relativistic effects, such as $\alpha$-Uranium, PuCoGa$_{5}$, and FePt, demonstrate robustness of the method. The comparison of the calculated forces with the ones obtained by numerical differentiation of free energy shows close agreement with deviations about 0.1% or less.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

3D modeling of n = 1 RMP driven heat fluxes on the SPARC tokamak PFCs using HEAT

3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.

3D heat flux↗
DOE OSTI2022

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗
DOE OSTI2024

Developing Platinum-Group-Metal-Free Catalysts for Oxygen Reduction Reaction in Acid: Beyond the Single Metal Site

This project is to develop M (x) -N-C catalysts with dense multiple metal center (MMC) sites to meet the DOE 2025 activity target of 0.044 mA/cm 2 at 0.9 V IR-free ., as well as other goals such as durability. We have made important contributions to both catalyst development and fundamental understandings of the active sites in M-N-C catalysts in this project. We successfully made M (x) -N-C catalysts with some multiple metal center (MMC) sites by combining ionothermal carbonization with chemical vapor deposition (CVD). These catalysts, however, are not as active as the most active single-atom Fe-N-C catalysts made by the similar CVD process, owing likely to the low site density and the presence of inorganic Fe species such as iron carbides and nanoparticles. The most significant accomplishments we achieved in this project are: (1) we unraveled the formation pathway of Fe-N 4 sites during pyrolysis step-by-step and identified the trans-metalation mechanism, in collaboration with Deborah Myers and her colleagues at Argonne National Laboratory (ANL); (2) inspired by this finding, we pioneered the CVD synthesis of M-N-C catalysts (M = Mn, Fe, and Co), in which the Fe-N-C catalyst by CVD demonstrated an ORR activity of 33 mA/cm 2 at 0.9 V in H 2 -O 2 proton exchange membrane fuel cells (PEMFCs), very close to the ultimate goal of 35 mA/cm 2 of our project. This catalyst is the first Fe-N-C catalyst that contains only D1 sites without the D2 sites; whereas D1 and D2 sites have been always identified by Mossbauer in previous Fe-N-C catalysts. This finding helps to understand what the D1 and D2 sites are and their roles in catalyzing the ORR. (3) by improving the mass transport of the carbon matrix prior to the CVD process, the revised Fe-N-C catalyst made by CVD delivered a maximum power density of 0.53 W/cm 2 in H 2 -air PEMFCs. The improvement strategy was partly inspired by the computational modeling work by Adam Weber from LBNL, the Co-PI of this project, by developing, coding, and exercising a continuum level model of transport phenomena within a PGM-free catalyst layer. The model demonstrated that local resistances combined with limited site density of the PGM-free catalyst can result in limiting currents and poor polarization performance. The model also gave design guidance for impact of ECSA and overall catalyst-layer thickness. However, both Fe-N-C and Co-N-C catalysts developed by CVD showed poor durability in PEMFCS, in comparison with the traditional M-N-C catalysts synthesized via regular pyrolysis process. Consequently, we did not achieve the proposed durability targets. Despite so, we believe the FeNC-CVD catalysts with the poor durability and D1 sites only provides an excellent platform to understand the degradation mechanism of Fe-N-C in PEMFCs, the most important challenge in the development of M-N-C catalysts.

08 HYDROGEN↗
DOE OSTI2025

University of Miami G-band Vapor Radiometer Calibration (UMGVR_CAL) Field Campaign Report

Cold-air outbreak (CAO) clouds in the Arctic are commonly mixed-phase (MP); however, the partitioning of the amount of ice and water in CAO clouds and precipitation is not always well observed. Understanding how cloud phases partition as a function of cloud life cycle is important for predicting snowfall rates, convective life cycle, and intensity at weather timescales. The partitioning into liquid versus ice also has radiative impacts that are consequential for climate. These concerns motivated the incorporation of an airborne G-band vapor radiometer (GVR) into an National Science Foundation-supported aircraft campaign named the Cold Air outbreak Experiment in the Sub-Arctic Region (CAESAR). The GVR is an upward-pointing passive microwave radiometer using four frequencies centered around the 183.31 GHz water vapor absorption line, displaced by +- 1, 3, 7, and 14 GHz. For context, The U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility operates a surface-based GVR at its North Slope of Alaska (NSA) site. The same GVR has been used previously for a field campaign in the southeast Pacific, where an offset was noticed between brightness temperatures (Tbs) measured under clear skies compared to those calculated from a radiative transfer model. To account for any calibration offsets to the Tbs, a request was made to DOE to allow the GVR to operate at ARM’s Southern Great Plains (SGP) observatory and enable comparisons between its measurements and those available at SGP. This request was granted, titled ‘UMGVR_CAL’, short for ‘UMGVR_Calibration’. From October 30 to November 10, 2023, the GVR was deployed to the ARM SGP site to take advantage of their regular, nearby radiosonde launches under clear-sky conditions. The latter were determined using the SGP total sky imagery data. The GVR brightness temperatures in these clear-sky conditions were compared to those calculated by a radiative transfer code (PAMTRA) based on the SGP radiosondes. During the campaign, four suitable clear-sky episodes could be used for the GVR calibration. While few in number, these proved to be enough to satisfy our goal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Numerical simulation of the hot-tail runaway electron production mechanism using CQL3D and comparison with Smith–Verwichte analytical model

Abstract The hot-tail mechanism of runaway electron (RE) production (Harvey et al 2000 Phys. Plasmas 7 4590) is the primary source of RE in the case of rapidly cooling tokamak plasma. Quantifying this mechanism is very important as it can provide most of the post-thermal-quench (TQ) current, or a seed current for the secondary source of RE through the avalanche mechanism. An analytic model which omits pitch-angle scattering is often used in literature for estimating the hot-tail RE density (Smith and Verwichte 2008 Phys. Plasmas 15 072502). In the present study, we use the CQL3D bounce-averaged Fokker–Planck code (Harvey and McCoy 1992 Proc. IAEA Technical Committee Meeting on Advances in Simulation and Modeling of Thermonuclear Plasmas p 527) to test the limits of validity of the model. In particular, we examine the cases of Z = 1 and Z = 18 ions, for sets of different initial temperature, density, electric field and the characteristic time of temperature decay. We show that for Z = 1 plasma, the ratio of RE density computed by CQL3D to that estimated from the model is within 0.6–6.0 in studied cases. For the Z = 18 case, this factor is systematically a much smaller number, typically 0.02–0.6. We suggest a simple correction to the model that narrows down the range of this ratio to 0.3–3.8 in all of the cases, including Z = 1 and Z = 18 plasmas.

Physics↗
DOE OSTI2022

SDSS-IV MaStar: [α/Fe] for the MaNGA Stellar Library from synthetic model spectra

ABSTRACT We calculate the α-enhancement ratio [α/Fe] for the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar) while also fitting for the fundamental atmospheric parameters effective temperature, surface gravity, and metallicity – Teff, log g, [Fe/H]. This approach builds upon a previous catalogue of stellar parameters, whereby only the fundamental atmospheric parameters are fit with solar-scaled models. Here, we use the same Markov Chain Monte Carlo method with the additional free parameter [α/Fe]. Using the full spectral fitting code pPXF, we are able to fit multiple lines sensitive to [α/Fe] for a more robust measurement. Quality flags based on the convergence of the sampler, errors in [α/Fe] and a cut in the χ2 of the model fit are used to clean the final catalogue, returning 17 214 spectra and values in the range of −0.25 < [α/Fe] < 0.48. Comparing our calculated [α/Fe] with literature values reveals a degeneracy in cool stars with log g ≥ ∼4; this comparison is then used to create an alternative and calibrated parameter set. We also plot the final catalogue in an [Fe/H] versus [α/Fe] diagram and recover the expected result of increasing [α/Fe] with decreasing [Fe/H] for Milky Way disc-halo stars. We apply our method to a subsample of spectra of uniform resolution and higher signal to noise that finds that our results are independent of this higher signal to noise. In the context of stellar population models, we are able to cover a parameter space for the creation of intermediate to old age models at solar-scaled [α/Fe], high [Fe/H] and enhanced [α/Fe], low [Fe/H].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

GEAR-MC and Differential-Operator Methods Applied to Electron-Photon Transport in the Integrated TIGER Series

The sensitivity analysis algorithms that have been developed by the radiation transport community in multiple neutron transport codes, such as MCNP and SCALE, are extensively used by fields such as the nuclear criticality community. However, these techniques have seldom been considered for electron transport applications. In the past, the differential-operator method with the single scatter capability has been implemented in Sandia National Laboratories’ Integrated TIGER Series (ITS) coupled electron-photon transport code. This work is meant to extend the available sensitivity estimation techniques in ITS by implementing an adjoint-based sensitivity method, GEAR-MC, to strengthen its sensitivity analysis capabilities. To ensure the accuracy of this method being extended to coupled electron-photon transport, it is compared against the central-difference and differential-operator methodologies to estimate sensitivity coefficients for an experiment performed by McLaughlin and Hussman. Energy deposition sensitivities were calculated using all three methods, and the comparison between them has provided confidence in the accuracy of the newly implemented method. Unlike the current implementation of the differential-operator method in ITS, the GEAR-MC method was implemented with the option to calculate the energy-dependent energy deposition sensitivities, which are the sensitivity coefficients for energy deposition tallies to energy-dependent cross sections. The energy-dependent cross sections could be the cross sections for the material, elements in the material, or reactions of interest for the element. Further, these sensitivities were compared to the energy-integrated sensitivity coefficients and exhibited a maximum percentage difference of 2.15%.

42 ENGINEERING↗
DOE OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗
DOE OSTI2023

Convergence of small scale Ly α structure at high- z under different reionization scenarios

ABSTRACT The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc and 5 − 80 h−1 kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

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

Current aqueous plutonium processing models for criticality safety often contain significant bias due to material modeling assumptions. These solutions include plutonium chloride solutions, which are modeled as fictitious plutonium metal-water mixtures because little is known about the actual density of the solution. Furthermore, there is no current predictive capability for modeling plutonium metal-water mixtures that is approved for use at Los Alamos National Laboratory (LANL). Recent density measurements for aqueous plutonium chloride systems (PuCl 3 -HCl-H 2 O) now allow for the development and application of a more realistic density law. This work develops a Python-based density law for this ternary solution using an empirical method. This code can be used in conjunction with an MCNP6 input to determine the density and composition of a solution based upon user inputs of plutonium concentration, hydrochloric acid concentration, and temperature. The tool allows users to input plutonium and acid content of a solution in terms of molality, molarity, or concentration, and predicts density within the current data range within 1.4% of experimental data. Current preliminary MCNP6 calculations utilizing this tool have demonstrated a minimum decrease in system reactivity of 5% in comparison to the current modeling conventions. Thus, this tool enables more accurate criticality safety operational limits by better crediting chorine content while still maintaining necessary conservatism.

Density Law↗
DOE OSTI2020

Hanging String Cuts in SPR Caverns: Modeling Investigation and Comparison with Sonar Data

Investigation of leaching for oil sales includes looking closely at cavern geometries. Anomalous cavern "features" have been observed near the foot of some caverns subsequent to partial drawdowns. One potential mitigation approach to reducing further growth of preexisting features is based on the hypothesis that reducing the brine string length via a "string cue' would serve to move the zone associated with additional leaching to a location higher up in the cavern and thus away from the preexisting feature. Cutting of the hanging string is expected to provide a control of leaching depth that could be used to "smooth" existing features and thus reduce geomechanical instability in that region of the cavern. The SANSMIC code has been used to predict cavern geometry changes (i.e., the extent of cavern growth with depth) based on variable input parameters for four caverns: West Hackberry 11 (WH11), West Hackberry 113 (WH113), Big Hill 104 (BH104), Big Hill 114 (BH114). By comparing the initial sonar geometry with resultant geometries calculated by the SANSMIC code, conclusions may be drawn about the potential impact of these variables on future cavern growth. Ultimately, these conclusions can be used to assess possible mitigation strategies such as the potential advantage of cutting versus not cutting a brine string. This work has resulted in a recommendation that a hanging string cut of 80 ft in WH11 would be beneficial to future cavern geometry, while there would be little to no benefit to string cuts in the other three caverns investigated here. The WH11 recommendation was followed in 2019, resulting in an operational string cut. A sonar performed after the string cut showed no adverse leaching in the area of the preexisting flare, as expected from the results of the preliminary SANSMIC runs described in this report. Additional SANSMIC modeling of the actual amount of injected raw water resulted in good agreement with the post-cut sonar.

02 PETROLEUM↗
DOE OSTI2020

Neutral leakage, power dissipation and pedestal fueling in open vs closed divertors

The effect of divertor closure on neutral leakage, power dissipation and pedestal fueling is explored using the SOLPS code for different densities including attached and detached divertor conditions for DIII-D open and closed divertors. Here, a metric for divertor closure in terms of the percentage of neutrals escaping the divertor is presented. At the detachment onset the population of neutrals escaping the divertor is 12% for the closed divertor and 35% for the open divertor. This results in a lower contribution of the convective term to the total heat flux and increased power losses in the closed divertor. The comparison between the two geometries reveals that for the same upstream density, the upstream temperature is 23% higher in the closed divertor. The physical mechanisms for detachment at lower upstream density in the closed divertor are also identified with the reduction of the convective term directly resulting in 40% decrease and the increased power losses contributing with a factor of 2. The baffling also lowers the neutral flux refueling the pedestal resulting in further reduction of the core density for achieving detachment. The results introduced here show that the divertor recycling is the primary fueling mechanism in DIII-D. It is shown that fueling occurring in the region between the X-point and the outer midplane is at least four times lower in the closed divertor due the baffling structure which enhances target neutral trapping and consequently reduced upstream leakage. Due to the openness of the inner leg to the common flux and HFS regions both open and closed divertor experience almost identical fueling at the HFS highlighting the need to also consider the inboard divertor for shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

iMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs – I. Construction and analysis of the mock data cubes

ABSTRACT Galaxy formation and evolution simulations are essential tools to probe poorly known astrophysics processes, but particular care is needed to compare simulations with galaxy observations, as observed data need to be modelled as well. We present a method to generate mock galaxies from the hydro-dynamical IllustrisTNG simulations which are suited to compare with integral field spectroscopic observation of galaxies from the SDSS-IV/Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. First, we include the same instrumental effects and procedures as adopted in the acquisition and analysis of real data. Furthermore, we generate the galaxy spectra from the simulations using new stellar population models based on the MaNGA stellar library (MaStar). In this way, our mock data cubes have the same spatial sampling, cover the same wavelength range (3600–10 300 Å), and share the same spectral resolution (R ≈ 1800) and flux calibration of real MaNGA galaxy spectra. In this first paper, we demonstrate the method over an early- and a late-type simulated galaxy from TNG50. We analyse the correspondent mock MaNGA-like data cubes with the same full spectral fitting code, firefly, which was used for the observed spectra. We find that the intrinsic and recovered age and metallicity gradients are consistent within 1σ, with residuals over all tassels consistent with 0 at the 68 ${{\ \rm per\ cent}}$ confidence level. We also perform the challenging test at comparing intrinsic and recovered star formation histories, finding a close resemblance between input and output. In follow-up papers, we will present a full simulated MaNGA-like catalogue (≈10 000 galaxies) with a comprehensive comparison of TNG50 simulations to MaNGA observational results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling Enhancements, Cross-Section Generation Updates, and Benchmarking with Shift

This technical report documents the modeling enhancements, cross-section generation updates, and bench marking with the Shift Monte Carlo code performed under the US Department of Energy Nuclear Energy Advanced Modeling and Simulation Program in FY 2024. The work performed included several modeling enhancements, such as integration of cross-section generation in Titan and the ability to produce microscopic multigroup cross sections with Shift. Benchmarking of the cross sections produced by Shift and the two-step workflow with Griffin was performed for three problems: the Advanced Breeder Test Reactor, a generic pebble bed reactor, and a TRISO heat pipe microreactor. Comparisons of results from these benchmark problems were done with Serpent, OpenMC, and Griffin. These enhancements provide a robust foundation for applying Shift for both reference and two-step neutronics analysis for advanced reactor simulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

2025 Review and Revision of Federal Guidance Report 15

Federal Guidance Report No. 15 (FGR 15), External Exposure to Radionuclides in Air, Water and Soil, published in 2019 and referred to below as FGR 2019, provides age-specific effective dose rate coefficients for reference persons externally exposed to each of 1252 radionuclides homogenously distributed in environmental media. Soon after completion of FGR 2019, the International Commission on Radiological Protection (ICRP) published a similar report (ICRP Publication 144, 2020) addressing the same radionuclides and many of the external exposure scenarios addressed in FGR 2019. In 2023 Argonne National Laboratory (ANL) published a review of dose coefficients in FGR 2019, concluding that “The external effective dose from beta radiation is not appropriately accounted for in FGR 15 for both low- and high-energy beta emitters in air, in soil, and on soil surfaces.” That conclusion was based largely on comparisons of dose coefficients in FGR 2019 with values in ICRP Publication 144 but also on consideration of some unexpected patterns of equivalent dose rate coefficients across tissues and of effective dose rate coefficients across different soil depths, for a selected set of low-energy beta emitters. In response to the ANL report, the Center for Radiation Protection Knowledge (CRPK) at Oak Ridge National Laboratory (ORNL) performed an extensive reexamination of the methods and published values of FGR 2019. CRPK found that coding errors had resulted in inaccurate estimates, primarily overestimates, in many of the dose coefficients tabulated in FGR 2019 but that many of the differences between results in FGR 2019 and ICRP Publication 144 could be traced to differences in methodology. CRPK has corrected all errors in FGR 2019 associated with the coding errors and has taken the opportunity to improve a major portion of the remaining dose coefficients in FGR 2019, primarily through additional Monte Carlo calculations resulting in improved statistics for tissue equivalent dose rate coefficients for exposure to monoenergetic sources. This has eliminated the need for extrapolation of results from high and medium energies to low energies, as was done in FGR 2019. In addition, inconsistencies in the methodology identified in the review, such as computational representation of the adult female, were eliminated. The revised effective dose rate coefficients are consistent with values in ICRP Publication 144 except for differences in values clearly arising from differences in methodology; and differences in values for a relatively small set of very low-energy radionuclides with highly uncertain dose coefficients regardless of methodology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cluster-dynamics-based parameterization for sulfuric acid–dimethylamine nucleation: comparison and selection through box and three-dimensional modeling

Clustering of gaseous sulfuric acid (SA) enhanced by dimethylamine (DMA) is a major mechanism for new particle formation (NPF) in polluted atmospheres. However, uncertainty remains regarding the SA–DMA nucleation parameterization that reasonably represents cluster dynamics and is applicable across various atmospheric conditions. This uncertainty hinders accurate three-dimensional (3-D) modeling of NPF and the subsequent assessment of its environmental and climatic impacts. Here we extensively compare different cluster-dynamics-based parameterizations for SA–DMA nucleation and identify the most reliable one through a combination of box model simulations, 3-D modeling, and in situ observations. Results show that the parameterization derived from Atmospheric Cluster Dynamic Code (ACDC) simulations, incorporating the latest theoretical insights (DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97X-D/6-311++G(3df,3pd) level of theory) and adequate representation of cluster dynamics, exhibits dependable performance in 3-D NPF simulation for both winter and summer conditions in Beijing and shows promise for application in diverse atmospheric conditions. Another ACDC-derived parameterization, replacing the level of theory with RI-CC2/aug-cc-pV(T+d)Z//M06-2X/6–311++G(3df,3pd), also performs well in NPF modeling at relatively low temperatures around 280 K but exhibits limitations at higher temperatures due to inappropriate representation of SA–DMA cluster thermodynamics. Additionally, a previously reported parameterization incorporating simplifications is applicable for simulating NPF in polluted atmospheres but tends to overestimate particle formation rates under conditions of elevated temperature (>∼300 K) and low-condensation sink (<∼3×10 -3 s −1 ). Our findings highlight the applicability of the new ACDC-derived parameterization, which couples the latest SA–DMA nucleation theory and holistic cluster dynamics, in 3-D NPF modeling. The ACDC-derived parameterization framework provides a valuable reference for developing parameterizations for other nucleation systems.

Shen, Jiewen↗
DOE OSTI2021

Numerical Simulation of Container Breach and Airborne Release of Solids Due to Mechanical Insults

Throughout U.S. Department of Energy (DOE) complexes, safety engineers employ the five-factor formula to calculate the source term (ST) that includes parameters of airborne release fraction (ARF), respirable fraction (RF) and damage ratio (DR). Limited experimental data on fragmentation of solids, such as ceramic pellets (i.e., PuO2), and container breach due to mechanical insults (i.e., drop and forklift impact), can be supplemented by modeling and simulation using high fidelity computational tools to estimate these parameters. This paper presents the use of Sandia National Laboratories' SIERRA solid mechanics (SM) finite element code to investigate the behavior of the widely utilized waste container (such as 7A Drum) subject to a range of free fall impact and puncture scenarios. The resulting behavior of the container is assessed, and the estimates are presented for bounding DRs from calculated breach areas for the various accident conditions considered. This paper also describes a novel multiscale constitutive model recently implemented in SIERRA/SM that simulates the fracture of brittle materials such as PuO 2 and determines ARF during the fracture process. Furthermore, comparisons are made between model predictions and simple bench-top experiments.

61 RADIATION PROTECTION AND DOSIMETRY↗