Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Single-gap isotropic s- wave superconductivity in single crystals AuSn 4

In this article, London, λ L (T), and Campbell, λ C (T), penetration depths were measured in single crystals of a topological superconductor candidate AuSn 4 . At low temperatures, λ L (T) is exponentially attenuated and, if fitted with the power law, λ(T) ~ T n , gives exponents n > 4, indistinguishable from the isotropic single s- wave gap Bardeen-Cooper-Schrieffer (BCS) asymptotic. The superfluid density fits perfectly in the entire temperature range to the BCS theory. The superconducting transition temperature, T c = 2.40 ± 0.05 K, does not change after 2.5 MeV electron irradiation, indicating the validity of the Anderson theorem for isotropic s- wave superconductors. Campbell penetration depth before and after electron irradiation shows no hysteresis between the zero-field cooling (ZFC) and field cooling (FC) protocols, consistent with the parabolic pinning potential. Interestingly, the critical current density estimated from the original Campbell theory decreases after irradiation, implying that a more sophisticated theory involving collective effects is needed to describe vortex pinning in this system. In general, our thermodynamic measurements strongly suggest that the bulk response of the AuSn 4 crystals is fully consistent with the isotropic s- wave weak-coupling BCS superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Octet scalars shaping LHC distributions in 4-jet final states

We study properties of a hypothetical scalar particle, Θ, which is a color octet and an electroweak singlet. At hadron colliders, Θ is pair produced through its QCD coupling to gluons, so that its mass determines the cross section. It decays at tree level into $q\bar{q}$ through dimension-5 operators, and at one loop into gluons. Thus, the main LHC signature of Θ is a pair of dijets of equal invariant mass. The CMS search in this channel shows a 3.6σ excess over the QCD background for a dijet mass M jj ≈ 0.95 TeV, which can be due to Θ : its production cross section (65 fb for a real scalar) and the acceptance of the CMS event selection applied to pp → ΘΘ → ($q\bar{q}$) ($q\bar{q}$) yield a rate consistent with the excess. Furthermore, the shape of the dσ/dM jj signal is in agreement with the CMS result. Given the data-driven background fit performed by CMS, we find that a complex scalar (whose production rate is twice as large) fits better the data than a real scalar. Besides the pair of dijets, testable LHC signals include a trijet-dijet topology, a $t\bar{t}$ pair plus a dijet resonance, as well as final states involving a Higgs, W or Z boson plus jets.

Dobrescu, Bogdan A. [Fermi National Accelerator La↗
DOE OSTI2020

A Systems Biology Approach to Energy Flow in H 2 Producing Microbial Communities: Microbial Mats

This report describes the results of work conducted at NASA Ames Research Center, in collaboration with investigators at Lawrence Livermore National Laboratory, during the time period 12/15/2014 to 12/15/2017. The work was an extension of a collaboration with LLNL that began in 2008. Work this reporting period focused primarily on tasks related to microbial mats from both the Elkhorn Slough and Baja California field sites, including the analysis and publication of data produced prior to 12/15/2014. During this reporting period, in addition to the microbial mat work, we also conducted a series of experiments with a particular algal-bacterial association (C. sorokiniana and A. brasilense). A number of aspects of this well-characterized algal-bacterial association enabled us not only to further our understanding of interactions between these particular organisms, but to apply that knowledge and techniques to our collection of microorganisms isolated from microbial mats. Microbial mat experiments focused both on intact mat communities and on cultures of organisms isolated from the mats (the Microbial Mat Living Library). In a series of co-culture experiments, we documented both positive and negative interactions amongst heterotrophic and autotrophic isolates of microbial mat microorganisms. One particularly important and novel mat cyanobacterium (ESFC-1) was studied in detail. We developed techniques to pursue an investigation of functional redundancy (the idea that closely related organisms perform identical roles in a community) in cyanobacterial autotrophs in the Elkhorn Slough microbial mat. In those experiments, we documented numerous instances of increases in growth rate between co-cultures of two organisms (relative to their growth rate in isolation) performing what would seem to be identical functional roles. Simplified microbial mat communities consisting of organisms from the Microbial Mat Living Library were constructed and studied in laboratory experiments. The “constructed mats” were grown in standard laboratory glassware and plasticware, but cultures were also combined with artificial materials in an effort to re-create the three-dimensional structure of naturally occurring microbial mats. These constructed mats re-created a number of functions attributable to natural microbial mats (N-fixation, H2-production) and simplified their study using molecular ecological techniques due to the smaller number of microbial “species” present. Cyanobacteria which migrated in response to vertical variations in light intensity (irradiance) obtained characteristic positions in constructed mats that were similar to their positions in natural mats. A detailed investigation of nitrogen cycling in both the Baja California and Elkhorn Slough microbial mats was initiated during this reporting period. The nitrogen cycle in the microbial mats has not been studied in depth, despite the fact that nitrogen usually limits growth in marine environments. The organisms conducting nitrogen cycling in both microbial mats were studied using the genomes available from earlier work on these mats, as well as with a targeted approach using molecular ecological methods developed here. Quantitative measurements of rates of nitrogen transformations were made using stable isotopically labeled nitrogen. The processes quantified and examined using molecular ecological methods included: nitrogen fixation, ammonification, ammonium and nitrate assimilation, denitrification, anaerobic oxidation of ammonium, and dissimilatory nitrate reduction to ammonium. In a series of experiments designed to examine a potential application for microbial mats, we demonstrated that both Elkhorn Slough and Baja California microbial mats were very effective in removing large quantities of both nitrate and ammonium from waste streams containing these compounds at levels where they would be considered to be pollutants. Considerable effort was expended during this reporting period to analyze, and prepare for publication, metagenomes and metatranscriptomes collected from two experimental manipulations of microbial mats. The first, a diel study, was conducted in November 2011, while the second, an investigation of the effects of sulfate on community composition, was conducted during this reporting period. For the earlier study, reference-based and reference-free methods were used to assess microbial ecology and genetic partitioning in these complex microbial systems, and that work is available as a pre-print. Four metagenomes from the same experiment were published. In the second, longer (nearly 200 day) manipulation, rRNA iTag libraries for cDNA and DNA were analyzed to determine the effects of sulfate concentration on microbial mat community composition over the course of the experiment. Work on the algal-bacterial association (Chlorella sorokiniana and Azospirillum brasilense) focused initially on demonstrating the transfer of carbon and nitrogen between the partners, using stable isotopic labeling and the LLNL NanoSIMS. Subsequently, the potential role of diazotrophy (nitrogen fixation) by the bacterial partner in supporting the growth of the association was assessed in a series of experiments. It appears to be the case that, although the bacterial partner demonstrably increases the fitness of the algae under nitrogen limitation, and cells of C. sorokiniana receive nitrogen originally fixed from atmospheric nitrogen by A. brasilense, the amount of nitrogen obtained by diazotrophy is too small to explain the increase in fitness.

08 HYDROGEN↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Machine learning applied to classifying neutron resonances

The performance of nuclear reactors and other nuclear systems depends on a precise understanding of the neutron interaction cross sections for materials used in these systems. These cross sections exhibit resonance structure whose shape is determined in part by the angular momentum quantum numbers of the resonances. The correct assignment of the quantum numbers of neutron resonances is therefore of paramount importance. In this project, we apply a machine learning technique, namely decision trees, to automate the quantum number assignments. The tree is trained from simulated data generated to mimic the errors found in real data. We explore the use of several physics-motivated features for training our trees, including the nearest neighbor spacing distribution, cumulative level distribution, and channel width distributions. Initial results using random matrix theory motivated fits which demonstrated that we can determine resonance spin groups somewhat reliably. If we use these fits as features in our trees, we can train them to spot outliers corresponding to misassigned resonances. We found that with the large number of features used in this project that the decision tree tended to over t training data resulting in poor performance with respect to the test data. By reducing the number of features, we can achieve nearly perfect assignment of quantum numbers with our training data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Broadband Characterization and Circuit Model Development of Transmission-Scale Transformers

This report describes broadband measurements of transmission-scale transformers typical in the electric power grid. This work was performed as part of the EMP Resilient Grid LDRD project at Sandia National Laboratories to generate circuit models that can be used for high-altitude electromagnetic pulse (HEMP) coupling simulations and response predictions. The objective of the work was to obtain characterization data of substation yard equipment across a frequency range relevant to HEMP. Vector network analyzer measurements up to 100 MHz were performed on two power transformers at ABB-Hitachi and a single ITEC potential transformer. Custom cable breakouts were designed to interface with the transformer terminals and provide ground connections to the chassis at the base of the transformer bushings. The three-phase terminals of the power transformers were measured as a common mode impedance using a parallel resistive splitter, and the single-phase terminals of the potential transformer were measured directly. A vector fitting algorithm was used to empirically fit circuit models to the resulting two-port networks and input impedances of the measured objects. Simplified circuit representations of the input impedances were also generated to assess the degree of precision needed for high-altitude electromagnetic pulse response predictions, which were performed in Sandia's XYCE circuit simulator platform. HEMP coupling simulations using the transformer models showed significant reduction in the voltage peak and broadening in the pulse width seen at the power transformer compared to the traveling wave voltage. This indicated the importance of the load condition when defining the coupled insult in an electric power substation. Simplified circuit models showed a similar voltage at the transformer with a smoothed waveform. The presence of potential transformers in the simulation did not significantly change the simulated voltage at the power transformer. Single-port input impedance models were also developed to define load conditions when transfer characteristics were not necessary.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Unitary Bethe-Salpeter Methods in Two- and Three-Body Systems

The study of resonances, unstable particles formed in particle scattering, is motivated by the questions about the strong interaction and the composition of hadrons. In order to advance this study in the intermediate energy region, a variety of techniques are used. In this thesis, we present work done on two different analyses that use the principle of unitarity and the Bethe-Salpeter equation to study the resonances Lambda(1405), Sigma(1385), and a1(1260). Firstly, we perform a simultaneous analysis of Sand P-waves of the Strangeness = -1 two-body meson-baryon scattering amplitude using all low-energy data. For the first time, differential cross section data are included for chiral unitary coupled-channel models. From this model, S- and P-wave amplitudes are extracted and we observe both well-known I(JP)=0(1/2^-) S-wave states as well as a new I(JP)=1(1/2^+) state absent in quark models and lattice QCD results. Multiple statistical and phenomenological tests suggest that, while the data clearly require an I =1 P-wave resonance, the new state just accounts for the absence of the decuplet Sigma(1385)3/2+ in the model which is subsequently included in the parameterization. Secondly, we formulate the final state interaction of the a1(1260) resonance decay in a manifestly three-body-unitary parameterization and fit it to the a1(1260) lineshape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model is a good test case because it can also be applied to other meson decays including exotic states and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity. The understanding of the resonances being analyzed in both the two-body work and the three-body work gives a better understanding of intermediate energy QCD.

Sadasivan, Daniel↗
DOE OSTI2021

Spectroscopic and Dynamical Studies of Highly Energized Small Polyatomic Molecules (End-of-Grant Report)

The aim of this research program has been to develop models for qualitative and quantitative understanding of complex processes in small molecules, such as vibronic coupling, predissociation, photodissociation, and isomerization. Simple models capture physical insights and lead to intuitive ideas that can enlighten a broad array of problems. Numerical fits can provide predictive power, but rarely impart insights into the underlying mechanisms responsible for the observed phenomena. This DOE supported program has prioritized a quantum mechanical understanding over a perfect numerical fit of spectra. Finding few-parameter models, when those parameters have physical meaning, allows the insights gained from one problem to be easily applied to related problems in other systems. The systems that have been under investigation include atmospherically relevant sulfur-containing compounds, various unimolecular dynamical processes of electronically excited acetylene, and the rovibrational dynamics of formaldehyde and SO 2 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

National Ignition Campaign Data Analysis by an ICF Outsider

A simple DT burn model, with no alpha-particle warming feedback, using a combination of ICF inputs and the measured ion temperature, T, volume compression, C, ablator velocity, v, and burn width time, Δt, can reproduce the national-ignition campaign (NIC) fusion yields with a rms percent difference of ~40%. Analysis suggestion the fusion-yield reproducibility of NIC events with nominally identically measured T, C, v, and Δt is <30%. The highest fusion yields, at a fixed laser energy are obtained with depleted-uranium (DU) hohlraums and high-density carbon (HDC) ablators. Small diameter fit tubes (≤5 μm) enhance fusion yields by an average of 26±15% relative to 10-μm diameter fit tubes. A dependence of fusion yield on the diameter of the capsule relative to the diameter of the hohlraum is inferred for fixed T, C, v, and Δt. No evidence of conditions close to ignition, or alpha-particle warming is found in the existing NIC data with LEH< 2MJ. To obtain a predictive model for extrapolation proposes we map the T, C, and v into power laws as a function of LEH. The above mentioned extrapolatable model has a rms percent difference from the yields of the second half of the NIC of (56 ± 4)%. Higher (and significantly lower) DT temperatures are obtained with HDC ablators (and low-foot drives). The highest compressions are obtained via low-foot drives, but this does not compensate for the corresponding lower temperatures. An extrapolation of the presented model suggests significant alpha-particle heating at an LEH of ~5 MJ. It is hoped that others more familiar it the conditions needed for ignition will use to presented model to estimate the LEH required for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

In Situ Inference for Earth System Predictability

An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

In Situ Inference for Earth System Predictability

Focal Area: Focal Area 3: Insight gleaned from complex simulated data using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

PBX9012 Tracker Gauge Data

Five shock-to-detonation transition (SDT) gas gun experiments for PBX 9012 were performed by Burns and Chiquete [2020]. The data from the experiments is available on the small-scale database. Here we focus on the tracker gauge data for the lead shock trajectory, x(t). In order to determine the shock speed U s (t) and the detonation transition point, a least squared fit to the shock trajectory is used based on a pair of ODEs with a fitting form for the shock acceleration developed by Hill and Gustavsen [2002].

42 ENGINEERING↗
DOE OSTI2021

Theoretical justification for heat flux limiter 0.15

In this section, we describe a semi-nonlocal kinetic model for the electron heat flux. The model has no free parameters. It takes the more physically motivated approach of placing a limiter on the perturbation to the electron distribution function, rather than the usual ad-hoc limiter on the heat flux. If the kinetic model is fitted by a heat flux limiter model, we find that the higher value of f e = 0.15 is a much closer fit than other models with f e ≃ 0.02.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SURFplus model calibration for PBX 9012

PBX 9012 is a plastic bonded explosive composed of 90.2 wt % HMX (cyclo-tetramethylene-tetranitramine) and 9.8 % Viton A (vinylidene-hexafluoropropylene copolymer). It has nearly the same composition as LX-07 (90% HMX and 10% Viton A) but with a slightly lower density. Hence it has a higher porosity and is slightly more shock sensitive. Here we describe the SURFplus model calibration for PBX 9012 using data from five shock-to-detonation transition (SDT) gas gun experiments by Burns and Chiquete, detonation front shape data from rate stick experiments and cylinder test experiment reported in Anderson et al., Jackson et al. The SDT experiments provide Pop plot data points (run distance-to-detonation) which characterize shock initiation, Lagrangian velocity time histories from embedded magnetic velocity gauges which give information on the reaction behind the lead shock, and Hugoniot data points for the reactants EOS. The rate stick experiments provide data on the curvature effect which characterize propagating detonation waves. The cylinder test experiments provide data on the detonation release isentrope used to calibrate the products EOS. The HE model uses a Davis reactants EOS with parameters that Burns and Chiquete [2020] fit to Hugoniot data, and Davis products EOS with parameters that Anderson et al. [2020] fit to cylinder test data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Supercritical Treatment Technology for Water Purification (Annual Final Report)

This project continued the development of the Supercritical Water Extraction – Enhanced Targeted Recovery (SWEETR™) technology, a novel desalination process for treating hypersaline brines. The overall objective is to demonstrate the technical and economic feasibility of an innovative, energy efficient, and robust supercritical desalination technology to treat hypersaline solutions and separating saltwater into a pure water stream and valuable recovered solids, resulting in zero liquid discharge (ZLD). The proposed technology, focuses on innovative methods of applying supercritical water to treat highly concentrated brine solutions without incurring a high energy penalty. Integrating the technology with solar energy reduces the energy cost for the system. The results obtained during the proof-of-concept testing demonstrate that the SWEETR™ concept is technically feasible and that supercritical water treatment is a viable option with strong merit for treating high salinity waters. The strategic design of the system provides the opportunity to localize the supercritical zone, which helps reduce the overall energy cost associated with the desalination process. The process can be tuned to remove only the amount of salts required to make the treated water “fit for purpose”, providing further opportunities to minimize energy costs. The ability to destroy organic compounds associated with the water while simultaneously being able to desalinate makes SWEETR™ uniquely fitted for treating aqueous streams such as produced water from oil and gas extraction which contain small amounts of organics. Such mixed contaminant streams pose an extreme challenge for other treatment options. Future work will further develop the SWEETR™ technology to take it from the bench-scale to the pilot-scale, and ready the technology for commercial application.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Fast Emulation of Expensive Simulations using Approximate Gaussian Processes [Slides]

Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗