Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discretization”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

$\mathrm{REX}$: An analytical tool for reactor operating envelope expansion through fuel-clad thermo-mechanics

Here REX, an analytical reactor design tool for the identification and expansion of reactor operating envelopes, is introduced and detailed. Written in MATLAB, REX takes core geometries and axially discretized detector and depletion output from the Serpent neutron transport code and performs full-core, pin-specific thermal-hydraulic and thermo-mechanical calculations. REX identifies the operating envelopes of advanced, solid-fuel nuclear reactors on the basis of their fuel-cladding, fuel-gas, and coolant-cladding interfacial thermo-mechanics. It iteratively pushes coolant inlet temperatures and modifies fuel pin geometries to induce mechanical failure in their original cladding materials at the fullest extent of their fuel cycle lengths and permissible coolant inlet temperatures. It then attempts to expand their operating envelopes by determining the mechanical responses of alternative cladding materials under the same geometric conditions. The result of the REX calculation sequence is a set of 5-D variables describing the temperatures, pressures, geometries, and mechanics of the core as functions of assembly, fuel pin, axial zone, depletion step, and coolant inlet temperature for each candidate cladding material with the limiting fuel pin geometry for the original material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Validation of URANS and STRUCT-ε turbulence models for stratified sodium flow

Simulations of a transient stratified sodium experiment are carried out using a classic unsteady RANS model and the second-generation URANS model, STRUCT-ε. Turbulence modeling challenges and their implications to stratified flow prediction are discussed in the context of other sources of error. Input errors are discussed and addressed; discretization error is calculated to be less than 5% of the inlet velocity, for 80% of the domain; and remaining errors in temperature distributions are attributed to the turbulence model. Qualitative flow features from the simulations are presented and discussed. Compared to the experiment, the STRUCT-ε turbulence model provides a more physically accurate prediction of temperature and momentum mixing in key regions of the domain. Quantitative measures such as the L 2 norm of the temperature discrepancy demonstrate the improved performance of the STRUCT-ε approach. In conclusion, the magnitude of the temperature fluctuations is very well-predicted by the STRUCT-ε, while URANS overpredicts them by approximately 50%.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗
DOE OSTI2025

A methodology for decay heat characterization in molten salt reactors

Accurate decay heat prediction in molten salt reactors (MSRs) faces dual challenges: complex operational uncertainties and the need for interpretable models compatible with engineering workflows. This work presents a hybrid machine learning and segmented polynomial methodology that addresses both requirements through three key innovations. First, a modular data architecture encodes MSR-specific operational parameters (power density: 1-100 W cm -3 , humidity: 0-0.1 wt %, air ingress: 0-0.1 mol %) with uncertainty-aware temporal discretization spanning 15 orders of magnitude. Second, region-optimized machine learning models achieve 92.3 % root mean square error (RMSE) reduction over conventional polynomials while maintaining physical interpretability through automated piecewise equation generation. Third, dual front-end interfaces accelerate safety analyses — a Jupyter environment enables researchers to explore 10,000+ parameter combinations via interactive widgets, while a Streamlit web application reduces design iteration cycles through production-grade visualization tools. Operational deployment demonstrates prediction times of only a couple hundred milliseconds for 10 4 years decay profiles, enabling real-time optimization of spent fuel container designs.

42 - ENGINEERING↗
DOE OSTI2025

Direct local parametrization of nuclear state densities using the back-shifted Bethe formula

Level densities are often parametrized using the back-shifted Bethe formula (BBF) for nuclei that possess experimental data for s-wave neutron resonance average spacings and a complete discrete level sequence at low excitation energies. However, these parametrizations require the additional modeling of the dependence of the spin-cutoff parameter on excitation energy. Here, in this work, we avoid the need to model the spin distribution of level densities by using the experimental data to parametrize directly the state densities, for which the BBF does not depend on the spin-cutoff parameter. This approach allows for a local parameterization of state densities that is independent of the spin-cutoff parameter. We provide these parameters in a tabulated form for applications in nuclear reaction calculations and for testing microscopic approaches to state densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Techno-economic and life-cycle analysis of strategies for improving operability and biomass quality in catalytic fast pyrolysis of forest residues

Many of the challenges faced by the first commercial biorefineries were associated with feedstock handling, quality, and cost. Strategies are needed to enable further expansion of biorefineries and meet the growing demand for bio-based fuels and products. Here, we examine 2 key feedstock challenges and mitigation strategies in the context of a catalytic fast pyrolysis (CFP) biorefinery: (1) the operability of the feed system, which may be improved by modifying the minimum particle size fed to the reactor, and (2) the quality of the biomass, which may be improved by employing air classification to remove undesirable material and increase fuel yields. We conduct techno-economic analysis (TEA) and life-cycle analysis for these strategies, employing a discrete event simulation model for biomass preprocessing combined with a series of correlations developed from literature data and a rigorous CFP conversion model. Our results highlight the importance of balancing increased cost and material losses from preprocessing against improved operability and fuel yields. Economics and sustainability were optimized when operating at the lowest minimum particle size, emphasizing the importance of minimizing material losses while maintaining the operability of the process. Economically, additional costs and material losses from air classification could be acceptable due to improved biomass conversion, and an optimum air classification speed was identified; however, the fuel GHG emissions were minimized when air classification was not used. Valorizing material removed during preprocessing as a coproduct could improve economics and sustainability, decreasing the burden of material losses.

09 - BIOMASS FUELS↗
DOE OSTI2021

Scalable semitransparent prototype organic photovoltaic module with minimal resistance loss

Semi-transparent organic photovoltaic (ST-OPV) cells are regarded as an attractive solution to building-integrated solar energy harvesting. Both the power conversion efficiency (PCE) and average photopic transmission (APT) of ST-OPVs have shown substantial increases in recent years. However, less attention has been paid to the area scaling of ST-OPV cells. In this work, we investigate the scalability of ST-OPV cells from 4 mm 2 to 1 cm 2 , as well as 9 cm 2 active area prototype module. By integrating both top and bottom metal grids onto the transparent electrodes, the series resistance loss of 1 cm 2 ST-OPV cell is significantly reduced, and is comparable to that of the grid-free 4 mm 2 cell. Nine 1 cm 2 cells are then connected in a series-parallel circuit to realize a prototype ST-OPV module. A 100% fabrication yield with only 5% PCE deviation among discrete cells is achieved. The semitransparent module shows PCE = 7.2 ± 0.1% under simulated AM 1.5G illumination at 1 sun intensity, which exhibits no connection resistance loss compared to that of the individual cells. The ST-OPV module exhibits an APT = 38.1 ± 1.1%, which enables a light utilization efficiency, LUE = 2.74 ± 0.09%. Here, the method demonstrates a promising way for ST-OPV modules to scale without compromising performance.

14 SOLAR ENERGY↗
DOE OSTI2020

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluating MPI resource usage summary statistics

The Message Passing Interface (MPI) remains the dominant programming model for scientific applications running on today’s high-performance computing (HPC) systems. This dominance stems from MPI’s powerful semantics for inter-process communication that has enabled scientists to write applications for simulating important physical phenomena. MPI does not, however, specify how messages and synchronization should be carried out. Those details are typically dependent on low-level architecture details and the message characteristics of the application. Therefore, analyzing an application’s MPI resource usage is critical to tuning MPI’s performance on a particular platform. The result of this analysis is typically a discussion of the mean message sizes, queue search lengths and message arrival times for a workload or set of workloads. While a discussion of the arithmetic mean in MPI resource usage might be the most intuitive summary statistic, it is not always the most accurate in terms of representing the underlying data. In this paper, we analyze MPI resource usage for a number of key MPI workloads using an existing MPI trace collector and discrete-event simulator. Our analysis demonstrates that the average, while easy and efficient to calculate, is a useful metric for characterizing latency and bandwidth measurements, but may not be a good representation of application message sizes, match list search depths, or MPI inter-operation times. Additionally, we show that the median and mode are superior choices in many cases. We also observe that the arithmetic mean is not the best representation of central tendency for data that are drawn from distributions that are multi-modal or have heavy tails. Furthermore, the results and analysis of our work provide valuable guidance on how we, as a community, should discuss and analyze MPI resource usage data for scientific applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Immortal rays: Rethinking random ray neutron transport on GPU architectures

The Random Ray Method (TRRM) is a recently developed adaptation of the Method of Characteristics for neutral particle transport simulations. TRRM has demonstrated excellent performance on 3D nuclear reactor benchmark problems using CPU-based compute systems. When porting to GPU-based systems, however, new performance challenges arise that are unique to processors targeting massive fine-grained parallelism. For smaller problems, or for large problems that are domain decomposed across many computational nodes, the problem size per node has insufficient parallelism to saturate GPU node resources, thus greatly limiting speedup. In this study, we report on a newly developed “immortal ray” variant of TRRM. Here, the immortal ray technique exposes significantly more fine-grained parallelism by fundamentally reformulating the numerical details of ray discretization, resulting in performance tradeoffs with significant overall benefit on GPUs. For very small 2D simulation problems we found the new immortal ray variant allowed for up to a 4.4x speedup when run on a single GPU. For larger 3D simulation problems we found the new variant improved strong scaling by 3x when run on the Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

NekRS, a GPU-accelerated spectral element Navier–Stokes solver

The development of NekRS, a GPU-oriented thermal-fluids simulation code based on the spectral element method (SEM) is described. For performance portability, the code is based on the open concurrent compute abstraction and leverages scalable developments in the SEM code Nek5000 and in libParanumal, which is a library of high-performance kernels for high-order discretizations and PDE-based miniapps. Critical performance sections of the Navier–Stokes time advancement are addressed. Performance results on several platforms are presented here, including scaling to 27,648 V100s on OLCF Summit, for calculations of up to 60B gridpoints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A review of rare earth elements and yttrium in coal ash: Content, modes of occurrences, combustion behavior, and extraction methods

Rare earth elements and yttrium (REY) have attracted considerable attention over the last decade because of their vital roles in clean energy, consumer product, national defense and security applications, among other uses. Due to the retention of REY during coal burning, coal combustion ash is considered as potential alternative sources for REY. Understanding the content, speciation, retention and/or transformation behavior of REY during coal combustion not only expands our knowledge of the combustion behavior of the trace elements in coal, but also provides basis for modeling REY partitioning during coal combustion and for developing economically viable REY recovery technologies. This review makes a critical summary of recent progress in the study of REY in coal ash. The contents and the extraction potentials of REY in coal ash derived from 15 major coal-producing countries worldwide were summarized and evaluated. Various analytical methods for determining REY bulk contents and speciation, together with the solid sample pretreatment, analytical accuracy and precision, advantages and disadvantages were summarized and compared. Modern analytical approaches combined indirect methods (e.g., sequential extraction) shed light on the physical distribution, mineralogy, and the chemical state of REY in coal ash. Three types of REY occurrences in coal ash, including Si-Al glassy association, discrete minerals or compounds, and organic association (bound with unburned carbon) were defined in the review. The glassy association can be further divided into REY minerals closely bound to glass phases and dispersed throughout the glassy structure. REY partitioning in various emission streams, the size distribution, and their enrichment behavior in coal ash were discussed. Additionally, thermal behavior and transformation of various REY forms in coal during combustion process, including organic-associated REY, REY phosphates, REY carbonates, clay-bound REY and among others were summarized. Two possible retention mechanisms of REY by aluminosilicate glass at boiler temperature were proposed: the incorporation of the individual REY phases into the glass as inclusions and the diffusion of REY phases throughout Si-Al glass structures in the melting process. Feed coal mineral types, mineral-mineral associations, boiler conditions, and other factors control the retention process. After coal combustion, the speciation of REY in fly ash may be modified by the reactions of REY phases with flue gas components. Further, an overview of REY transformation mechanisms during coal combustion was deeply discussed. Finally, current extraction techniques for REY recovery from coal combustion ash were introduced. Future outlooks and research problems were also identified.

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

Coarse grained simulation of convectively driven turbulent mixing, transition, and turbulence decay

Accurate predictions with quantifiable uncertainties are essential to many practical turbulent flow applications exhibiting extreme geometrical complexity and broad ranges of length and time scales. Under-resolved computer simulations are typically unavoidable in such applications, and implicit large-eddy simulation (ILES) often becomes the effective strategy. In this work, we focus on assessing ILES initialized with well-characterized 25 6 3 homogeneous isotropic turbulence datasets generated with direct numerical simulation (DNS). ILES solutions based on the LANL xRAGE code are studied as function of resolution for 6 4 3 , 12 8 3 , 25 6 3 , and 51 2 3 grids. ILES performance of new directionally-unsplit high-order numerical hydrodynamics algorithms in xRAGE is examined in this context. Compared to the initial 25 6 3 DNS, we find longer inertial subranges and higher turbulence Reynolds number (Re) for 25 6 3 and 51 2 3 xRAGE — attributed to having linked DNS (Navier–Stokes based) solutions to nominally inviscid (higher Re) Euler based ILES solutions. For fixed 25 6 3 resolution, we find that significantly larger simulated turbulence Re can be achieved with the higher-order unsplit (vs. split) discretizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Direct ab initio calculation of the 4 He nuclear electric dipole polarizability

The calculation of nuclear electromagnetic sum rules by directly diagonalizing the nuclear Hamiltonian in a large basis is numerically challenging and has not been performed for A>2 nuclei. With the significant progress of high performance computing, we show that calculating sum rules using numerous discretized continuum states obtained by directly diagonalizing the ab initio no-core shell model Hamiltonian is achievable numerically. Specifically, we calculate the 4 He electric dipole (E1) polarizability, that is an inverse energy weighted sum rule, employing the Daejeon16 NN interaction. We demonstrate that the calculations are numerically tractable as the dimension of the basis increases and are convergent. Our results for the 4 He electric dipole polarizability are consistent with the most recent experimental data and are compared with those of other theoretical studies employing different techniques and various interactions.

Astronomy & Astrophysics↗
DOE OSTI2025

Point cloud-based diffusion models for the Electron-Ion Collider

At high-energy collider experiments, generative models can be used for a wide range of tasks, including fast detector simulations, unfolding, searches of physics beyond the Standard Model, and inference tasks. In particular, it has been demonstrated that score-based diffusion models can generate high-fidelity and accurate samples of jets or collider events. This work expands on previous generative models in three distinct ways. First, our model is trained to generate entire collider events, including all particle species with complete kinematic information. We quantify how well the model learns event-wide constraints such as the conservation of momentum and discrete quantum numbers. We focus on the events at the future Electron-Ion Collider, but we expect that our results can be extended to proton-proton and heavy-ion collisions. Second, previous generative models often relied on image-based techniques. The sparsity of the data can negatively affect the fidelity and sampling time of the model. We address these issues using point clouds and a novel architecture combining edge creation with transformer modules called Point Edge Transformers. Third, we adapt the foundation model OmniLearn, to generate full collider events. This approach may indicate a transition toward adapting and fine-tuning foundation models for downstream tasks instead of training new models from scratch.

Araz, Jack Y. [Stony Brook Univ., NY (United State↗
DOE OSTI2026

Variational neural network approach to QFT in the field basis

We present a variational neural network approach for solving quantum field theories in the field basis, focusing on the free Klein-Gordon model formulated in momentum space. While recent studies have explored neural-network-based variational methods for scalar field theory in position space, a systematic benchmark of the analytically solvable Klein-Gordon ground state—particularly in the momentum-space field basis—has been lacking. In this work, we represent the ground-state wavefunctional as a neural network defined on a discretized set of field configurations and train it by minimizing the Hamiltonian expectation value. This framework enables direct comparison to exact analytic results for a range of key observables, including the ground-state energy, two-point correlators, expectation value of the field, and the structure of the learned wavefunctional itself. Our results provide quantitative diagnostics of accuracy and establish a validated foundation for extending neural-network wavefunctional methods to interacting field theories and position-space formulations.

Klein-Gordon model↗