Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Improved Tensor Current Limit from 8 B 𝛽 Decay Including New Recoil-Order Calculations

A precision measurement of the 𝛽 + decay of 8 B was performed using the Beta-decay Paul Trap to determine the 𝛽−𝜈 angular correlation coefficient 𝑎 𝛽⁢𝜈 . The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, 𝑎 𝛽⁢𝜈 = −0.3345 ± 0.001⁢9 stat ± 0.002⁢1 syst . This value agrees with the standard model value of −1/3 and improves uncertainties in 8 B by nearly a factor of 2. By combining results from 8 B and 8 Li , a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision 𝛽 decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3⁢𝜎. Finally, the present results are of comparable sensitivity and do not support this finding.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Precise 3D reactor core calculation using spherical harmonics and discontinuous Galerkin finite element methods

We study the use of P{sub N} method in angle and discontinuous Galerkin is space to solve 3D neutron transport problem. P{sub N} method consists in developing the angular flux on truncated spherical harmonics basic. In this paper, we couple this method with the discontinuous finite elements in space to obtain a complete discretization of the multigroup neutron transport equation. To investigate its precision, the method was applied to Takeda and C5G7 benchmark problems. These calculations point out that the proposed P{sub N}-DG method is capable of producing accurate solutions in small computational time, and that it is able to handle complex 3D geometries. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations [Slides]

TSLs for H(BeH 2 ) and Be(BeH 2 ) were evaluated in the incoherent approximation using FLASSH. Minimum critical mass estimates were performed using MC21 continuous energy Monte Carlo calculations with the BeH 2 TSLs. Researchers confirmed very low estimated critical mass for BeH 2 -moderated and Be- reflected spheres. These were the lowest critical masses NNL has studied to date. Future work includes relaxation of the incoherent approximation for Be(BeH 2 ). A small increase in critical mass is expected due to increased neutron transmission at sub-Bragg cutoff energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Theoretical and calculable dependent variables and their covariance in nuclear data libraries [Slides]

This presentation begins by defining observables, including theoretical observable, calculable observable, and measured observable. It also provides definitions and examples of experimental effects. Additionally, uncertainty in evaluated libraries and a metric to probe evaluated uncertainty is presented. In conclusion, being related to theoretical quantities such as resonance parameters, the uncertainty in ENDF libraries (e.g., ENDF/B-VIII.0) can be overestimated because it is usually evaluated under the guidance of experimental uncertainty: lim Δσexp.corr.→0 Δσ = Δσ theoretical model . The presentation states that a clear distinction of the uncertainty between nuclear theoretical models and experimental corrections should be revisited and this is important since (e.g., transport) simulations need theoretical quantities convoluted with specific operational parameters and material configurations. It also states that strong coupling between uncertainty quantification methodologies and optimization procedures exists, and it is necessary to develop methodologies to obtain uncertainty on theoretical models from physical and mathematical constraints.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Computational Materials Design for Ceramic Nuclear Waste Forms Using Machine Learning, First-Principles Calculations, and Kinetics Rate Theory

Ceramic waste forms are designed to immobilize radionuclides for permanent disposal in geological repositories. One of the principal criteria for the effective incorporation of waste elements is their compatibility with the host material. In terms of performance under environmental conditions, the resistance of the waste forms to degradation over long periods of time is a critical concern when they are exposed to natural environments. Due to their unique crystallographic features and behavior in nature environment as exemplified by their natural analogues, ceramic waste forms are capable of incorporating problematic nuclear waste elements while showing promising chemical durability in aqueous environments. Recent studies of apatite- and hollandite-structured waste forms demonstrated an approach that can predict the compositions of ceramic waste forms and their long-term dissolution rate by a combination of computational techniques including machine learning, first-principles thermodynamics calculations, and modeling using kinetic rate equations based on critical laboratory experiments. By integrating the predictions of elemental incorporation and degradation kinetics in a holistic framework, the approach could be promising for the design of advanced ceramic waste forms with optimized incorporation capacity and environmental degradation performance. Such an approach could provide a path for accelerated ceramic waste form development and performance prediction for problematic nuclear waste elements.

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

Understanding initial zeolite oligomerization steps with first principles calculations

Zeolites are porous aluminosilicate materials that find various applications in the chemical industry in separations, catalysis, ion exchange, and so forth. However, despite their widespread use, the reaction mechanisms occurring during zeolite growth are still unclear. In the study, we use density functional theory calculations to gain insights into the thermodynamics of oligomerization, which constitute the initial steps of zeolite growth. By taking into consideration solvent and temperature effects, our results demonstrate that the growth of aluminosilicate systems is significantly more exothermic than their pure silicate counterparts. Under pH neutral conditions, water prefers to dissociate on the early-growth-stage aluminosilicate complexes rather than desorb, thus generating potential Brønsted acid sites on the oligomers. Additionally, (alumino)silicate growth pathways are evaluated in the presence of Na + cation, as well as the Ca 2+ cation for the pure silicate pathway. The presence of cations increases the exothermicity of growth, with Ca 2+ exhibiting the most energetically favorable growth environment for the silicate systems. Importantly, we demonstrate through reaction extent analysis that the presence of cations modulates the speciation of the formed oligomers, with Na + favoring linear species in addition to the generally preferred cyclic ones. Overall, this work provides a fundamental understanding of the thermodynamics of complex reaction paths that occur during early stages of zeolite growth and suggests that the initial growth steps can have significant impact on the final zeolite structure.

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

Minimum reflux calculation for multicomponent distillation in multi‐feed, multi‐product columns: Mathematical model

Abstract Multi‐feed, multi‐product distillation columns are ubiquitous in multicomponent distillation systems. The minimum reflux ratio of a distillation column is directly related to its energy consumption and capital cost. Thus, it is a key parameter for distillation systems design, operation, and comparison. In this series, we present the first accurate shortcut based algorithmic method to determine the minimum reflux condition for any general multi‐feed, multi‐product (MFMP) distillation column separating any ideal multicomponent mixture. The classic McCabe‐Thiele or Underwood method is a special case of this general approach. Compared with existing techniques, this method does not involve any rigorous tray‐by‐tray calculation, nor does it require guessing of key components. In this first part of the series, we present the mathematical model for a general MFMP column, derive constraints for feasible separation and minimum reflux condition, discuss their geometric interpretations, and present an illustrative example to demonstrate the effectiveness of our approach.

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

Composite Correlated Molecular Orbital Theory Calculations of Ring Strain for Use in Predicting Polymerization Reactions

Strained ring systems play an important role in synthesis and can be characterized by the ring strain energy (RSE). The RSE of 3, 4, 5, and 6 membered saturated and unsaturated ring systems containing N, O, P, and S heteroatoms and H, F, SiMe 3 , and SO 2 Me substituents were calculated at the G3(MP2) composite correlated molecular orbital theory level using up to 5 models to predict the RSE. Generally, the RSE decreased as ring size increased with a substantial decrease from 4 to 5 membered rings. Replacement of a ring CH 2 with P or S reduced the RSE, consistent with less angle strain. The RSE for unsaturated systems were generally greater than for saturated systems due to increased angle strain. No general trends were found with respect to substituent effects. In conclusion, the RSE values suggest that 3-pyrroline and 2-pyrroline and their derivatives may be able to support ring opening metathesis polymerization and warrant further study.

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

Updating Critical Temperature Calculations for Several Secondary High Explosives

In designing experiments where high explosives (HEs) are heated, it is important to have an understanding of where thermal runaway may occur. This determination is often done by using the Frank-Kamenetskii (FK) equation. For several highly-studied HEs, the parameters necessary for FK calculations are typically referenced from decades-old literature, and more recent experimental data have shown that these values require adjustment. For example, some of the legacy parameters dangerously overpredict the critical temperature by tens of degrees relative to more recent observed values. Herein, we first summarize historical results and highlight the importance of insulation effects when estimating critical temperatures for small samples, which is especially relevant to small-scale aging, compatibility, and characterization experiments. Here, we then present updated FK parameters for several common secondary CHNO HEs like RDX, HMX, PETN, and TATB, as well as the first reported values for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105). Our updated parameters produce critical-temperature curves that are consistent with both mm-scale differentical scanning calorimetry (DSC) measurements and cm-scale data from the Lawrence Livermore National Laboratory one-dimensional time-to-explosion (ODTX) experiment. In analyzing the critical temperature versus sample size curves derived from the transcendental FK equation, we discovered that these curves are remarkably well described by a simple power-law function with a universal power of 0.142 and a HE-specific scaling factor.

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

Calculating the angle between jet axes

We study the angle between i) the standard jet axis, ii) the axis of a jet which has been groomed using soft drop, with reduced sensitivity to soft radiation, iii) the jet axis obtained with the winner-take-all recombination scheme, which is insensitive to soft radiation at leading power. We calculate the distributions for these angles at next- to-leading logarithmic accuracy, including non-global logarithms. The angle between the standard and groomed jet axis directly probes soft wide-angle radiation, leading to a novel factorization formula. This angle is also very sensitive to nonperturbative physics, which is directly connected to nonperturbative contribution to the rapidity anomalous dimension for transverse momentum distributions. Comparing our predictions to Pythia we find good agreement, and we foresee applications to jet substructure in proton-proton and heavy ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Lattice QCD calculation of the pion distribution amplitude with domain wall fermions at physical pion mass

We present a direct lattice QCD calculation of the x-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing a ≈ 0.084 fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction x ϵ [0.25, 0.75] at the largest pion momentum P z ≈ 1.85 GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and a = 0.076 fm. By comparison we find with 2σ confidence level that the DA obtained from chiral fermions is flatter and lower near x = 0.5.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An analysis of Bayesian estimates for missing higher orders in perturbative calculations

With current high precision collider data, the reliable estimation of theoretical uncertainties due to missing higher orders (MHOs) in perturbation theory has become a pressing issue for collider phenomenology. Traditionally, the size of the MHOs is estimated through scale variation, a simple but ad hoc method without probabilistic interpretation. Bayesian approaches provide a compelling alternative to estimate the size of the MHOs, but it is not clear how to interpret the perturbative scales, like the factorisation and renormalisation scales, in a Bayesian framework. Recently, it was proposed that the scales can be incorporated as hidden parameters into a Bayesian model. In this paper, we thoroughly scrutinise Bayesian approaches to MHO estimation and systematically study the performance of different models on an extensive set of high-order calculations. We extend the framework in two significant ways. First, we define a new model that allows for asymmetric probability distributions. Second, we introduce a prescription to incorporate information on perturbative scales without interpreting them as hidden model parameters. We clarify how the two scale prescriptions bias the result towards specific scale choice, and we discuss and compare different Bayesian MHO estimates among themselves and to the traditional scale variation approach. Finally, we provide a practical prescription of how existing perturbative results at the standard scale variation points can be converted to 68%/95% credibility intervals in the Bayesian approach using the new public code MiHO.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

CGC for ultra-peripheral Pb+Pb collisions at the Large Hadron Collider: a more realistic calculation

We provide the first calculation of two-gluon production at mid-rapidity in ultra-peripheral collisions in the Color Glass Condensate framework. To estimate systematic uncertainty associated with poor understanding of the wave function of the nearly real photon, we consider two diametrically different models: the dilute quark-antiquark dipole approximation and a vector meson, in which color charge density is approximated by McLerran-Venugopalan model. In the experimentally relevant range, the target nucleus can be faithfully approximated by a highly saturated state. This simplification enables us to perform efficient numerical simulations and extract the two-gluon correlation functions and the associated azimuthal harmonics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Method for calculating radon activity and radon rejection using a beta-gamma detector

Radon interferes with concentration measurements used by atmospheric radioxenon systems. We demonstrate a method to quantify the amount of radon that is present in the detectors, the impact of radon activity on the minimum-detectable-concentrations, and how to determine the needed radon rejection levels. An example calculation shows a radon rejection level of 10 5 is sufficient to limit impact on the detector sensitivity. We anticipate this method will give analysist a better understanding of radon present in their measurements and allow system designers to tailor their systems’ radon rejection better for its location.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Interlaboratory Reproducibility of Contour Method Data Analysis and Residual Stress Calculation

While the contour method for residual stress assessment has developed rapidly, no published study documents its interlaboratory reproducibility. Objective: Here we report an initial reproducibility experiment focused on contour method data analysis and residual stress calculation. The experiment uses surface topography data from a physical process simulation of elastic-plastic beam bending. The simulation provides surface topography, for input to the contour method data analysis, as well as a known residual stress field with 130 MPa peak magnitude. Additionally, to increase realism, noise and specific artifacts are added to the topography data. A group of participants received the topography data (without the known residual stress), independently analyzed the data, and submitted results as a two-dimensional residual stress field. Analysis of submissions provides a group average residual stress field and the spatial distribution of reproducibility standard deviation. The group average residual stress agrees with the known stress in magnitude and spatial trend. The reproducibility standard deviation ranges from 2 to 54 MPa over the measurement plane, with an average of 5.4 MPa. Reproducibility standard deviation is smaller in the cross-section interior (≤ 5 MPa), modest near local extrema in the stress field (5 to 10 MPa), and larger near the cross-section boundaries (10 to 30 MPa). Overall, the largest values of reproducibility standard deviation (up to 54 MPa) occur in limited areas where artifacts had been added to the topography data; while some participants identified and removed these artifacts, some did not, leading to systematic differences that elevated the standard deviation.

42 ENGINEERING↗
DOE OSTI2020

Anharmonic multi-phonon nonradiative transition: An ab initio calculation approach

Nonradiative carrier recombinations at deep centers in semiconductors are of great importance for both fundamental physics and device engineering. Here, we provide a revised analysis of Huang’s original nonradiative multi-phonon (NMP) theory with ab initio calculations. First, we confirmed at the first-principles level that Huang’s concise formula gives the same results as the matrix-based formula, and that Huang’s high-temperature formula provides an analytical expression for the coupling constant in Marcus theory. Secondly, we correct for anharmonic effects by taking into account local phonon-mode variations for different charge states of a defect. The corrected capture rates for defects in GaN and SiC agree well with experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Temperature-Dependent Configurational Entropy Calculations for Refractory High-Entropy Alloys

The cluster expansion formalism for alloys is used to construct surrogate models for three refractory high-entropy alloys (NbTiVZr, HfNbTaTiZr, and AlHfNbTaTiZr). These cluster expansion models are then used along with Monte Carlo methods and thermodynamic integration to calculate the configurational entropy of these refractory high-entropy alloys as a function of temperature. Many solid solution alloy design guidelines are based on the ideal entropy of mixing, which increases monotonically with N , the number of elements in the alloy. However, our results show that at low temperatures, the configurational entropy of these materials is largely independent of N , and the assumption described above only holds in the high-temperature limit. This suggests that alloy design guidelines based on the ideal entropy of mixing require further examination.

36 MATERIALS SCIENCE↗