Search NASA⌕ Search

SEARCH · Search NASA

Results for “Experimental and numerical techniques”

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 235 records · Page 13

Micromechanical response of SiC-OPyC layers in TRISO fuel particles

Tristructural isotropic (TRISO)–coated particle fuel is a proposed fuel for multiple advanced reactor concepts. The performance of the particle depends on whether the silicon carbide (SiC) layer remains intact to prevent the release of metallic and gaseous fission products. Mechanical fracture of the SiC layer is a potential failure mode under various fuel configurations and operating environments, including the potential transmission of matrix-originating cracks through TRISO particles. Furthermore, this study uses instrumented indentation techniques on cross-sectioned surrogate particles to examine the mechanical stability of the critical interface between SiC and the outer pyrolytic carbon (OPyC) layer. The observed behavior at the interface is rationalized by examining the radially dependent fracture behavior of the SiC layer and performing a numerical analysis to quantify the residual stresses that develop during the processing and cross-sectioning of the as-fabricated particle. Characterizing the SiC-OPyC interface of surrogate TRISO particles using nanoindentation provides unique insight into the interface's room-temperature residual stress and mechanical stability. The modeling efforts were used to investigate the experimental procedure further, and the results are presented herein to validate this fuel form's potential mechanical failure modes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of flexible proton beam imaging energy spectrometers (PROBIES)

A simple repeating step-filtered structure has been designed and deployed for the detection of MeV-energy proton beams produced from high-intensity short-pulse laser-driven experiments. Due to the diagnostic's ability to recover energy-resolved images of proton beams, it is called the proton beam imager and energy spectrometer (or PROBIES). The flexible design has two useful applications. First, it can be fielded in a configuration with a scintillator and imaging setup in order to collect proton beam information in high-repetition rate experiments. Here we show a configuration that collects the equivalent of a nine-film radiochromic film (RCF) stack, and a technique for using a neural network (NN) to rapidly analyze the images on millisecond time-scales. Second, when using RCF stacks as the detector, the energy resolution capability is greatly increased compared to traditional RCF stacks, while simultaneously reducing the total number of films required to make each measurement. We demonstrate a configuration that provides up to 90 spatially- and energy-resolved datapoints using just ten films. Both designs presented here can be easily customized for anticipated experimental outputs and manufactured at low cost. Such designs allow for rapid prototyping using modern high-accuracy 3D printers or on-demand computer numerical control machining to reduce the reliance on time-consuming film media, analysis, or bulky electro-magnetically dispersing spectrometers. The analysis of data is discussed and a methodology for constructing a NN that is capable of rapidly analyzing data on-the-fly with greater than 95% accuracy is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Electron Tomography Capabilities for Shale Imaging

Despite the advantageous resolution of electron tomography (ET), reconstruction of three-dimensional (3D) images from multiple two-dimensional (2D) projections presents several challenges, including small signal-to-noise ratios, and a limited projection range. This study evaluates the capabilities of ET for thin sections of shale, a complex nanoporous medium. A numerical phantom with 1.24 nm pixel size is constructed based on the tomographic reconstruction of a Barnett shale. A dataset of 2D projection images is numerically generated from the 3D phantom and studied over a range of conditions. First, common reconstruction techniques are used to reconstruct the shale structure. The reconstruction uncertainty is quantified by comparing overall values of storage and transport metrics, as well as the misclassification of pore voxels compared to the phantom. We then select the most robust reconstruction technique and we vary the acquisition conditions to quantify the effect of artifacts. We find a strong agreement for large pores over the different acquisition workflows, while a wider variability exists for nanometer-scale features. Furthermore, the limited projection range and reconstruction are identified as the main experimental bottlenecks, thereby suggesting that sample thinning, advanced holders, and advanced reconstruction algorithms offer opportunities for improvement.

3D scanning transmission electron microscopy↗

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extending the Takagi-Taupin equations for x-ray nanobeam Bragg coherent diffraction

We present an approach for simulating x-ray nanobeam Bragg coherent diffraction patterns based on the Takagi-Taupin equations. Compared to conventional methods, the current approach can be universally applied to any weakly strained system including semi-infinite crystals that diffract dynamically. It addresses issues such as the curved wave front and redivergence of the focused incident beam. In conclusion, we show excellent agreement with experimental data for a strained La 0.7 ⁢Sr 0.3⁢ MnO 3 thin film on a SrTiO 3 substrate and a path to extracting physical information using automatic differentiation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Near surface defects: Cause of deficit between internal and external open-circuit voltage in solar cells

The presence of interface recombination in a complex multilayered thin-film solar structure causes a disparity between the internal open-circuit voltage (V OC,in ), measured by photoluminescence, and the external open-circuit voltage (V OC,ex ) i.e. an additional V OC deficit. Higher V OC,ex value aim require a comprehensive understanding of connection between V OC deficit and interface recombination. Here, a deep near-surface defect model at the absorber/buffer interface is developed for copper indium di-selenide solar cells grown under Cu excess conditions to explain the disparity between V OC ,in and V OC,ex . The model is based on experimental analysis of admittance spectroscopy and deep-level transient spectroscopy, which show the signature of deep acceptor defect. Further, temperature-dependent current-voltage measurements confirm the presence of near surface defects as the cause of interface recombination. The numerical simulations show strong decrease in the local V OC ,in near the absorber/buffer interface leading to a V OC deficit in the device. This loss mechanism leads to interface recombination without a reduced interface bandgap or Fermi level pinning. Further, these findings demonstrate that the V OC ,in measurements alone can be inconclusive and might conceal the information on interface recombination pathways, establishing the need for complementary techniques like temperature dependent current voltage measurements to identify the cause of interface recombination in the devices.

14 SOLAR ENERGY↗

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗

Dynamic mode decomposition with core sketch

With the increase in collected data volumes, either from experimental measurements or high fidelity simulations, there is an ever-growing need to develop computationally efficient tools to process, analyze, and interpret these datasets. Modal analysis techniques have gained great interest due to their ability to identify patterns in the data and extract valuable information about the system being considered. Dynamic mode decomposition (DMD) relies on elements of the Koopman approximation theory to compute a set of modes, each associated with a fixed oscillation frequency and a decay/growth rate. Extracting these details from large datasets can be computationally expensive due to the need to implement singular value decomposition of the input data matrix. Sketching algorithms have become popular in numerical linear algebra where statistical theoretic approaches are utilized to reduce the cost of major operations. A sketch of a matrix is another matrix, which is significantly smaller, but still sufficiently approximates the original system. We put forth an efficient DMD framework, SketchyDMD, based on a core sketching algorithm that captures information about the range and corange (their mutual relationship) of input data. The proposed sketching-based framework can accelerate various portions of the DMD routines, compared to classical methods that operate directly on the raw input data. We conduct numerical experiments using the spherical shallow water equations as a prototypical model in the context of geophysical flows. In conclusion, we show that the proposed SketchyDMD is superior to existing randomized DMD methods that are based on capturing only the range of the input data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dipolar Spin Ice Regime Proximate to an All-In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2 Sn 2 O 7

The dipolar-octupolar (DO) pyrochlores, R 2 M 2 O 7 ( R = Ce , Sm , Nd ), are key players in the search for realizable novel quantum spin liquid (QSL) states as a large parameter space within the DO pyrochlore phase diagram is theorized to host QSL states of both dipolar and octupolar nature. New single crystals and powders of Ce 2 Sn 2 O 7 , synthesized by hydrothermal techniques, present an opportunity for a new characterization of the exchange parameters in Ce 2 Sn 2 O 7 using the near-neighbor X Y Z model Hamiltonian associated with DO pyrochlores. Utilizing quantum numerical linked cluster expansion fits to heat capacity and magnetic susceptibility measurements, and classical Monte Carlo calculations to the diffuse neutron diffraction of the new hydrothermally grown Ce 2 Sn 2 O 7 samples, we place Ce 2 Sn 2 O 7 ’s ground state within the ordered dipolar all-in-all-out (AIAO) Néel phase, with quantum Monte Carlo calculations showing a transition to long-range order at temperatures below those accessed experimentally. Indeed, our new neutron diffraction measurements on the hydrothermally grown Ce 2 Sn 2 O 7 powders show a broad signal at low scattering wave vectors, reminiscent of a spin ice, in striking contrast from previous powder neutron diffraction on samples grown from solid-state synthesis, which found diffuse scattering at high scattering wave vectors associated with magnetic and suggested an octupolar quantum spin ice state. We conclude that new hydrothermally grown Ce 2 Sn 2 O 7 samples host a finite-temperature proximate dipolar spin ice phase, above the expected transition to AIAO Néel order. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficient prediction of attosecond two-colour pulses from an X-ray free-electron laser with machine learning

Abstract X-ray free-electron lasers are sources of coherent, high-intensity X-rays with numerous applications in ultra-fast measurements and dynamic structural imaging. Due to the stochastic nature of the self-amplified spontaneous emission process and the difficulty in controlling injection of electrons, output pulses exhibit significant noise and limited temporal coherence. Standard measurement techniques used for characterizing two-coloured X-ray pulses are challenging, as they are either invasive or diagnostically expensive. In this work, we employ machine learning methods such as neural networks and decision trees to predict the central photon energies of pairs of attosecond fundamental and second harmonic pulses using parameters that are easily recorded at the high-repetition rate of a single shot. Using real experimental data, we apply a detailed feature analysis on the input parameters while optimizing the training time of the machine learning methods. Our predictive models are able to make predictions of central photon energy for one of the pulses without measuring the other pulse, thereby leveraging the use of the spectrometer without having to extend its detection window. We anticipate applications in X-ray spectroscopy using XFELs, such as in time-resolved X-ray absorption and photoemission spectroscopy, where improved measurement of input spectra will lead to better experimental outcomes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GEN3D Experimental and Numerical Development of GEN3 Durability Models

Understanding of the high temperature durability of particles and the materials that contain them is critical to next generation of concentrating solar power (CSP) technology. Here we studied the durability of particles and their containment materials under extreme UV cycling, thermal cycling, and in low-speed high temperature mechanical wear situations. The optical stability of seven candidate particles has been determined following exposure to the high temperature conditions present in a generation 3 particle-based CSP technology. Particle solar weighted absorptance and emittance measured periodically during 10,000 high solar irradiance exposure cycles and up to 400 hours of isothermal aging has been documented. The particle aging due to repeated exposure to concentrated solar flux represented the 30-year lifetime of a power plant. Models were fit to the absorptivity and emissivity data following the isothermal aging provides the projected optical degradation of the particles as a function of temperature. Mechanical wear was studied through the use of custom developed wear testing facilities for measuring high temperature impact wear, abrasion wear, and particle attrition. Additionally, a novel technique for measuring the high temperature mechanical properties of single particles was developed. Through these tests it was observed that high nickel alloys generally showed lower wear than comparable iron based steels, particularly at elevated temperatures of 800°C. Mechanical wear at these temperatures is a highly complex phenomenon combining both mechanical wear and oxidation. Additionally, the containment materials wear rates are influenced by the particles (both hardness and roundness), making the wear mechanisms complex. The initial wear test conducted in the abrasion test rig revealed a substantial amount of oxide materials in the particle bed after testing (in relative to later tests), and substantially more wear, likely indicating a need for concern in startup operation of particle facilities to not incur high wear from the presence of oxides. Particle attrition experiments have only been conducted for a single material but increase size distribution, reduction in circularity, and particle diameter is observed. Efforts to develop predictive models was limited due to a testing campaign that prioritized testing materials for particle pathway developers over building a comprehensive design of experiments. The results discussed in this report inform future CSP developers and researchers further de risking the technology and assisting in its future development. ParticleBased CSP development provides a path to dispatchable solar power generation with storage at a price competitive in the current energy market. Lowering the cost of CSP technology provides a carbon free power generation solution that can assist in the transition from fossil fuels to renewable sources of electricity.

14 SOLAR ENERGY↗

Optimization-Inspired Pin-Fin Array for Supercritical Carbon Dioxide Recuperator

Additively manufactured heat exchangers are one possible route to cost-effective sCO 2 power cycles. In this paper, experimental results are obtained for two helical pin fin tubes that were designed following parametric optimization of the fin array. Neither the numerical optimization nor the experimental testing have been previously reported in the literature. The two tube designs were, (1) the optimization inspired design, and (2) the optimization-inspired design with a fin diameter increased by a factor of 2. To characterize the print, the designs were scanned using X-ray computed tomography to measure feature sizes and heat transfer area. The optimization was conducted in a commercial, computational fluid dynamics code. The code solved the Reynolds Averaged Navier-Stokes (RANS) and energy equations with turbulence closure provided by the shear stress transport (SST) k-ω model. In the experiments, the Nusselt number augmentation was measured using the Wilson plot technique and the friction factor was determined with mass flow and pressure drop measurements. The experimental testing indicated that the optimization-inspired design had a friction factor that was four times less than the baseline tube design at equal Nusselt number. Additionally, the optimization-inspired design had a 14% improvement in Nusselt number, at equal friction factor, relative to the best performing computational fluid dynamics (CFD) trial points. Tube design (2), with the larger diameter pin fins, had similar performance, within experimental error, as the tube with the smaller diameter pin fins (1). Both tubes achieved overall fin array efficiencies near 1. A performance factor, V/V 0 , equal to the volume of the enhanced heat exchanger divided by the volume of the baseline (no-fins) heat exchanger, is recommended to quantify internal cooling performance. In conclusion, the experimental shell and tube heat exchanger, using the additively manufactured tube, is competitive with printed circuit heat exchangers in its pressure drop class and could be further improved by optimizing a shell-and-tube heat exchanger utilizing this heat transfer enhancement feature.

42 ENGINEERING↗

Excitation functions and isomeric cross-section ratios of (d,xn) reactions on 86 Sr

Excitation functions of the 86 Sr(d,n) 87m,87g Y, 86 Sr(d,2n) 86m,86g Y and 86 Sr(d,3n) 85m,85g Y reactions on enriched 86 Sr target were measured by the activation technique up to deuteron energies of 49 MeV. The isomeric cross-section ratios as a function of projectile energy were deduced from the measured data for 87m Y, 87g Y(cum), 86m Y, 86g Y(cum), 85m Y and 85g Y pairs for the same energy range. All measurements are reported for the first time. The experimental data were compared with the data from the TENDL library which is based on TALYS calculation with default parameters. No satisfactory agreement was observed. Nuclear model calculations were then performed using the codes TALYS and EMPIRE with some parameter adjustments, and compared with the experimental data. The quality of the agreement between experimental data and model calculations was numerically quantified. In general, the data as well as the isomeric cross-section ratios are partially reproduced by the model calculations, provided the input model parameters are properly chosen and the level structure of the product nucleus is thoughtfully considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Quantifying thermal barrier coating delamination through luminescence modeling

Thermal barrier coatings (TBCs) are widely used in the hot sections of gas turbine systems as they are remarkably efficient in insulating the underlying superalloys, leading to higher operating temperatures and therefore enhancing the performance of the engines. However, this benefit is only possible if the integrity of the TBC under aggressive thermo-mechanical environments is ensured. Delamination is a common but hard to detect failure mode. We present, in this work, supported by experimentation, the implementation of a modeling approach applying the Kubelka-Munk theory to provide numerical quantification of luminescence contrast and intensity due to top coat delamination in TBCs. The method relies on the drastic change in reflectivity when a delamination forms, exploiting it for high-contrast luminescence mapping. Two distinct TBC configurations containing an erbium-doped yttria-stabilized zirconia (YSZ:Er 3+ ) layer for delamination sensing were used to validate this model. A delamination zone induced by Rockwell indentation was successfully tracked by measuring an increase of the intensity of the erbium emission line at 562 nm. Luminescence-based methods for delamination detection can provide a revolutionary non-invasive technique, with potential for both off-line and on-line engine monitoring.

36 MATERIALS SCIENCE↗

Bi-fidelity modeling of uncertain and partially unknown systems using DeepONets

Recent advances in modeling large-scale, complex physical systems have shifted research focuses towards data-driven techniques. However, generating datasets by simulating complex systems can require significant computational resources. Similarly, acquiring experimental datasets can prove difficult. For these systems, often computationally inexpensive, but in general inaccurate models, known as the low-fidelity models, are available. Here in this paper, we propose a bi-fidelity modeling approach for complex physical systems, where we model the discrepancy between the true system's response and a low-fidelity response in the presence of a small training dataset from the true system's response using a deep operator network, a neural network architecture suitable for approximating nonlinear operators. We apply the approach to systems that have parametric uncertainty and are partially unknown. Three numerical examples are used to show the efficacy of the proposed approach to model uncertain and partially unknown physical systems.

MATHEMATICS AND COMPUTING↗

Simulation-Based Inference for Neutrino Interaction Model Parameter Tuning

High-energy physics experiments studying neutrinos rely heavily on simulations of their interactions with atomic nuclei. Limitations in the theoretical understanding of these interactions typically necessitate ad hoc tuning of simulation model parameters to data. Traditional tuning methods for neutrino experiments have largely relied on simple algorithms for numerical optimization. While adequate for the modest goals of initial efforts, the complexity of future neutrino tuning campaigns is expected to increase substantially, and new approaches will be needed to make progress. In this paper, we examine the application of simulation-based inference (SBI) to the neutrino interaction model tuning for the first time. Using a previous tuning study performed by the MicroBooNE experiment as a test case, we find that our SBI algorithm can correctly infer the tuned parameter values when confronted with a mock data set generated according to the MicroBooNE procedure. This initial proof-of-principle illustrates a promising new technique for next-generation simulation tuning campaigns for the neutrino experimental community.

Tame-Narvaez, Karla Maria [Fermilab]↗