Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “significant quantity”

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 325 records · Page 18

Transport error estimation using residual Monte Carlo

The residual Monte Carlo (RMC) method is also known in the literature as sequential Monte Carlo and reduced-source Monte Carlo. Given a Monte Carlo method for solving a linear equation and an approximate solution to that system, the residual method enables use of essentially the same Monte Carlo algorithm to directly compute the additive error or “defect” associated with the approximate solution. As the size of the defect decreases relative to the size of the solution, the residual Monte Carlo method becomes increasingly efficient relative to the standard Monte Carlo (SMC) method. Here we present a new RMC algorithm for evaluating the space-angle error in S n radiation transport solutions, and provide computational examples demonstrating that it can be far more efficient than SMC for this purpose. Herein we also describe a particular pitfall that must be avoided if RMC is to be efficient, and explain why the performance of RMC can significantly differ between different transport problems and different quantities of interest for the same problem.

97 MATHEMATICS AND COMPUTING↗

Improved manufacturability of vacuum thermoforming molds by segmentation

Thin polymer shell components used in consumer goods, the automotive industry, and marine applications are often produced by thermoforming. Depending on the quantity of parts, molding processes can increase throughput significantly when compared to processes such as milling or additive manufacturing. One disadvantage to mold-based manufacturing methods, however, is that molds must be produced to form the required geometry. Mold production is often complex and requires a large capital investment. For large-scale molding processes, mold production is sometimes complicated by the requirement for large working volume machine tools which are capable of machining the cavities to the appropriate tolerances for the selected molding process. Here, this paper focuses on a segmented design for a large draw ratio (>3:1 molded surface area to sheet surface area) thermoforming mold and a cost-benefit analysis of segmentation versus monolithic design. Discussion of the thermoforming tests are included. The primary goals were to improve manufacturability and reduce the overall cost of the mold.

Manufacturability↗

Training quantum neural networks using the quantum information bottleneck method

Abstract We provide in this paper a concrete method for training a quantum neural network to maximize the relevant information about a property that is transmitted through the network. This is significant because it gives an operationally well founded quantity to optimize when training autoencoders for problems where the inputs and outputs are fully quantum. We provide a rigorous algorithm for computing the value of the quantum information bottleneck quantity within error ε that requires O ( log 2 ⁡ ( 1 / ϵ ) + 1 / δ 2 ) queries to a purification of the input density operator if its spectrum is supported on { 0 } ⋃ [ δ , 1 − δ ] for δ > 0 and the kernels of the relevant density matrices are disjoint. We further provide algorithms for estimating the derivatives of the QIB function, showing that quantum neural networks can be trained efficiently using the QIB quantity given that the number of gradient steps required is polynomial.

Çatlı, Ahmet Burak (ORCID:0000000152294141)↗

Role of boundary conditions in quantum computations of scattering observables

Quantum computing may offer the opportunity to simulate strongly interacting field theories, such as quantum chromodynamics, with physical time evolution. This would give access to Minkowski-signature correlators, in contrast to the Euclidean calculations routinely performed at present. However, as with present-day calculations, quantum computation strategies still require the restriction to a finite system size, including a finite, usually periodic, spatial volume. In this work, we investigate the consequences of this in the extraction of hadronic and Compton-like scattering amplitudes. Using the framework presented in Briceño et al., we estimate the volume effects for various 1+1D Minkowski-signature quantities and show that these can be a significant source of systematic uncertainty, even for volumes that are very large by the standards of present-day Euclidean calculations. We then present an improvement strategy, based in the fact that the finite volume has a reduced symmetry. This implies that kinematic points, which yield the same Lorentz invariants, may still be physically distinct in the periodic system. As we demonstrate, both numerically and analytically, averaging over such sets can significantly suppress the unwanted volume distortions and improve the extraction of the physical scattering amplitudes. As the improvement strategy is based only in kinematics, it can be applied without detailed knowledge of the system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleon tensor form factors at large N c

We investigate nucleon tensor form factors in the large- N c limit. In this picture, the nucleon emerges as a state of the N c valence quarks, which were bound by pion mean fields that were created by the presence of the valence quarks self-consistently. We find that the tensor charge ( g T u − d = 0.99 ) and the anomalous tensor magnetic moment ( κ T u + d = 7.61 ) are dominated by valence quarks, while the tensor quadrupole moment ( Q T u − d = − 7.02 ) shows significant sea quark effects. We examine how these quantities vary as the average size of the pion mean field is changed, showing interpolation between nonrelativistic quark and Skyrme limits. We also observe that g T u − d and κ T u + d depend weakly on the pion mass. In contrast, Q T u − d exhibits strong enhancement near the chiral limit. The numerical results are in good agreement with available lattice quantum chromodynamics (QCD) data and provide predictions for unmeasured quantities. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Reaction Rate Distribution Measurements in ZPPR-15

Small metal foils were irradiated in ZPPR-15 to measure distributions for 239 Pu(n,f), 235 U(n,f), 238 U(n,f), 238 U(n,γ), 115 In(n,n’) and 115 In(n,γ) reaction rates in ZPPR-15. The measured data were reported as basic reaction rate distributions for individual foils, cell-averaged reaction rate distributions and spectral indices or reaction rate ratios. All of the measured data were analyzed and tabulated in such a manner that the data can be readily compared with computed values. The published uncertainties for ZPPR-15 the reaction rate measurements include the statistical uncertainties in the measurements, reproducibility uncertainty and small uncertainties related to corrections for items such as foil compositions and post-irradiation decay between reactor shutdown and foil measurement. There are additional uncertainty components related to detector calibration that are not included in the published uncertainties. A full uncertainty analysis was performed for each foil measurement and for quantities derived from the basic foil measurements. The significant uncertainties were quantified, and a total uncertainty was derived for each measured value. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models for the foil reaction rate measurements in ZPPR-15 loadings 123, 134 and 203 and for the fission chamber measurements of 237 Np and 242 Pu fission rates in ZPPR-15 loading 141. Because there were many foils distributed in the ZPPR-15 assembly and because an individual foil occupies a negligible volume in the ZPPR-15 assembly, direct Monte Carlo calculations of the measured foil reaction rates would require an impractical computational effort. The as-built Monte Carlo models and related files were prepared for conversion to deterministic models that can be calculated with realistic computational effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Feasibility demonstration of a road vehicle fueled with hydrogen-enriched gasoline

Evaluation of the concept of using hydrogen-enriched gasoline in a modified internal combustion engine in order to make possible the burning of ultralean mixtures. The use of such an engine in a road vehicle demonstrated that the addition of small quantities of gaseous hydrogen to gasoline resulted in significant reductions in exhaust emissions of carbon monoxide and nitrogen oxides as well as in thermal efficiency improvements of the engine performance.

Hoehn, F. W.↗

Very-long-baseline radio interferometry - The Mark III system for geodesy, astrometry, and aperture synthesis

Up to 112 megabit/sec from each radio telescope of an interferometric array can be recorded and processed by the Mark III VLBI system. Sample results are given for baseline lengths between three antennas in the U.S. and three in Europe, as well as for the arc lengths between the positions of six extragalactic radio sources. No significant change is detected in any of these quantities. Signals of a given polarization or of pairs of orthogonal polarizations may be recorded in up to 28 contiguous bands, each nearly 2 MHz wide, for mapping the brightness distribution of compact radio sources. The demonstrated ability to record large bandwidths, and to link many large radio telescopes, allows compact sources with flux densities below 1 millijansky to be detected and studied.

Rogers, A. E. E.↗

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

Disruption thermal load mitigation with shattered pellet injection on the Joint European Torus (JET)

Disruption mitigation remains a critical, unresolved challenge for ITER. To aid in addressing this challenge, a shattered pellet injection (SPI) system was installed on JET and experiments conducted at a range of thermal energy fractions and stored energies in excess of 7 MJ. The primary goals of these experiments were to investigate the efficacy of the SPI on JET and the ability of the plasma to assimilate multiple pellets. Single pellet injections produced a saturation in total radiated energy (Wrad) with increasing injected neon content, suggesting total radiation of stored thermal energy. Further increases in injected neon quantities resulted in reduced cooling times and current quench (CQ) durations, indicating higher impurity assimilation. No significant variation in CQ duration or W rad was observed when varying the deuterium content at fixed neon quantities. Additionally, higher assimilation, inferred by shorter CQ durations, was measured when a mechanical punch was used to launch the pellets and this was attributed to a lower pellet velocity leading to higher solid content in the pellet plume and larger fragments penetrating deeper into the plasma. Radiation asymmetries averaged over the cooling time were inferred from Emis3D and ranged from 1.6 to 1.9. Asymmetries averaged over the entire disruption sequence were found to increase at higher thermal energy fractions. The radiated energy fractions decreased with increasing thermal energy fractions but this trend was eliminated when toroidal asymmetries were accounted for with Emis3D. Pure deuterium pellets were able to produce cooling times of up to 75 ms with a gradual loss in thermal stored energy of up to 80%. Experiments with multiple pellet injection indicated W rad can be increased through pellet superposition and density can be increased with an additional D2 injection without a reduction in W rad . KPRAD modelling accurately reproduced the cooling times and the CQ duration at high thermal energies. Assimilation estimates from KPRAD indicated CQ rates scale strongly whilst W rad scales weakly and saturates with assimilated neon content. Comparable W rad can be achieved with lower assimilated neon quantities as longer cooling times are attained. Thus reduced neon content can be preferential in a thermal load mitigation scheme as it may reduce radiation asymmetries and prevent flash melting.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct numerical simulations and modeling of a spatially-evolving turbulent wake

Understanding of turbulent free shear flows (wakes, jets, and mixing layers) is important, not only for scientific interest, but also because of their appearance in numerous practical applications. Turbulent wakes, in particular, have recently received increased attention by researchers at NASA Langley. The turbulent wake generated by a two-dimensional airfoil has been selected as the test-case for detailed high-resolution particle image velocimetry (PIV) experiments. This same wake has also been chosen to enhance NASA's turbulence modeling efforts. Over the past year, the author has completed several wake computations, while visiting NASA through the 1993 and 1994 ASEE summer programs, and also while on sabbatical leave during the 1993-94 academic year. These calculations have included two-equation (K-omega and K-epsilon) models, algebraic stress models (ASM), full Reynolds stress closure models, and direct numerical simulations (DNS). Recently, there has been mutually beneficial collaboration of the experimental and computational efforts. In fact, these projects have been chosen for joint presentation at the NASA Turbulence Peer Review, scheduled for September 1994. DNS calculations are presently underway for a turbulent wake at Re(sub theta) = 1000 and at a Mach number of 0.20. (Theta is the momentum thickness, which remains constant in the wake of a two dimensional body.) These calculations utilize a compressible DNS code written by M. M. Rai of NASA Ames, and modified for the wake by J. Cimbala. The code employs fifth-order accurate upwind-biased finite differencing for the convective terms, fourth-order accurate central differencing for the viscous terms, and an iterative-implicit time-integration scheme. The computational domain for these calculations starts at x/theta = 10, and extends to x/theta = 610. Fully developed turbulent wake profiles, obtained from experimental data from several wake generators, are supplied at the computational inlet, along with appropriate noise. After some adjustment period, the flow downstream of the inlet develops into a fully three-dimensional turbulent wake. Of particular interest in the present study is the far wake spreading rate and the self-similar mean and turbulence profiles. At the time of this writing, grid resolution studies are underway, and a code is being written to calculate turbulence statistics from these wake calculations; the statistics will be compared to those from the ongoing PIV wake measurements, those of previous experiments, and those predicted by the various turbulence models. These calculations will lead to significant long-term benefits for the turbulence modeling effort. In particular, quantities such as the pressure-strain correlation and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. Present turbulence models do a very good job at prediction of the shape of the mean velocity and Reynolds stress profiles in a turbulent wake, but significantly underpredict the magnitude of the stresses and the spreading rate of the wake. Thus, the turbulent wake is an ideal flow for turbulence modeling research. By careful comparison and analysis of each term in the modeled Reynolds stress equations, the DNS data can show where deficiencies in the models exist; improvements to the models can then be attempted.

Cimbala, John M.↗

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗

Characterization of Spent Purolite A530E Resin with Implications for Long-Term Radioactive Contaminant Removal

Direct removal of pertechnetate anion (TcO4-) from the subsurface contaminated plumes at legacy nuclear weapons production sites, is effective remediation strategy to prevent their spread to adjacent aquifers but is also challenging because TcO4- is a trace component of the contaminated groundwater typically containing large access of other anionic constituents, such nitrate, sulfate, chloride and others. It can be achieved through ion exchange treatment using Purolite A530E resin which to date however has only been evaluated under laboratory conditions and many questions regarding its long-term operational performance remain. To address this need, here we report comprehensive characterization of the spent Purolite A530E resin which processed over 5.38x109 L of contaminated groundwater and successfully removed about 3.78 Ci of Tc-99 during four years of operation at the 200 West Pump & Treat facility at the US DOE Hanford site. This Tc-99 loading constitutes however only about 0.034% of the theoretical capacity of the resin which retained significant amounts of sulfate. Among other radioactive contaminants, small quantities of U-238, Co-60, and I-129 were retained by the resin. The total loading of iodine (combined 1-127 and I-129) on the spent resin exceeded that of Tc-99. To elucidate the mechanism of iodine retention, ion exchange behavior of iodide and iodate was investigated. Purolite A530E resin exhibited highly efficient uptake of iodide and only moderate affinity for iodate in accord with their Gibbs energy of hydration. Sorption isotherms for both anions obeyed Freundlich model.

Levitskaia, Tatiana G.↗

Comparison of Results from Recent NNSA and CEA Measurements of the 239 Pu(n, f) Prompt Fission Neutron Spectrum

The National Nuclear Security Administration (NNSA)/DP French Alternative Energies and Atomic Energy Commission (CEA)/DAM agreement on cooperation on fundamental science is a U.S.-French collaborative effort to combine intellectual and experimental resources and further the relevant nuclear science. Recently, both the NNSA and CEA experimental teams performed high-statistics measurements of the 239 Pu(n, f) prompt fission neutron spectrum (PFNS) at the Los Alamos Neutron Science Center, both of which were recently published in the journal Physical Review C. These separate measurements used the same experimental area and a common neutron detector array, but differ in many aspects, including background assessments, data acquisition systems and philosophies, fission detectors, and PFNS extraction techniques. Hence, some aspects of the experimental methods and associated uncertainties are highly correlated while others are independent. The results from both measurements broke new ground for PFNS measurements given their higher accuracy and more detailed study of corrections necessary for the measured quantity compared to existing literature measurements, and both will significantly impact PFNS nuclear data evaluations for the foreseeable future. Here, the focus of this work is to document a comparison of the results from these distinct measurements in terms of the acquired data, the PFNS results, and the measured average PFNS energies. While systematic differences between the PFNS results are present on the 1–3% level, the acquired data relative to each respective measurement at low incident neutron energies are in remarkable agreement, as are the conclusions regarding the magnitude and position of features in the PFNS relating to second-chance fission, third-chance fission, and pre-equilibrium neutron emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Role of Nuclear Data Sensitivities in Prompt α-Eigenvalue Predictions of Delayed Critical Benchmarks

Alpha (α) eigenvalues, which describe the logarithmic time derivative of the neutron population in a multiplying system, are integral to time-dependent behavior and diagnostic applications. However, uncertainties in the evaluated nuclear data can significantly impact the accuracy of transport simulations for such quantities. This work explores the use of machine learning models to predict two key outputs, α-eigenvalues and keff bias, using input features derived from α-eigenvalue sensitivities to nuclear data. The criticality safety benchmark models used in this study come from the International Handbook of Evaluated Criticality Safety Benchmark Experiments. Three models, random forest, XGBoost, and NGBoost, are trained on both energy-resolved and energy-summed α sensitivities. For the α-eigenvalue bias prediction, NGBoost achieved the highest R 2 (0.9476) using energy-resolved features, while XGBoost performed best using summed sensitivities. In contrast, when predicting the keff bias, all the models showed moderate predictive capability (best R 2 ≈ 0.72), as the mapping from the static α-sensitivities to the static keff bias was less direct. SHAP (SHapley Additive exPlanations) analysis was used to interpret the model predictions. Across both prediction tasks, the features associated with neutron capture [H-1 (n, γ)], uranium scattering reactions (such as 235 U elastic/inelastic), and actinide capture/fission reactions (such as 239 Pu and 234 U) were consistently identified as the most impactful. This highlights the key role of specific nuclear reactions and energy ranges in shaping both time-dependent and steady-state criticality behavior. These results demonstrated that α-sensitivities, despite being computed for time-dependent metrics, can provide valuable insights for predicting both α-eigenvalues and the keff bias. Moreover, machine learning models offer a promising pathway for uncovering important nuclear data dependencies and guiding future data evaluation efforts.

Nuclear data↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Direct Measurement of the Be 7 L / K Capture Ratio in Ta-Based Superconducting Tunnel Junctions

We report a high-statistics measurement of the L / K orbital electron capture ratio in 7 Be embedded in cryogenic Ta. The thin Ta film formed part of a high-resolution superconducting tunnel junction radiation detector that was used to identify the signals from different decay channels. The measured L/K capture ratio of 0.070(7) is significantly larger than the only previous measurement of this quantity and the theoretical predictions that include in-medium effects. This value is a uniquely sensitive probe of the 1s and 2s orbital overlaps with the nucleus and is of relevance to nuclear and atomic physics, as well as Li production in novae and other astrophysical scenarios. Furthermore, this is the first experiment that uses superconducting tunnel junctions for nuclear-recoil detection, opening a new experimental avenue for low-energy precision measurements with rare isotopes.

6 ≤ A ≤ 19↗

Measuring ground movement in geothermal areas of Imperial Valley, California

Significant ground movement may accompany the extraction of large quantities of fluids from the subsurface. In Imperial Valley, California, one of the potential hazards of geothermal development is the threat of both subsidence and horizontal movement of the land surface. Regional and local survey nets are being monitored to detect and measure possible ground movement caused by future geothermal developments. Precise measurement of surface and subsurface changes will be required to differentiate man-induced changes from natural processes in this tectonically active region.

Lofgren, B. E.↗