Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical benchmark”

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,261 records · Page 70

Power System Event Classification and Localization Using a Convolutional Neural Network

Detection and timely identification of power system disturbances are essential for situation awareness and reliable electricity grid operation. Because records of actual events in the system are limited, ensemble simulation-based events are needed to provide adequate data for building event-detection models through deep learning; e.g., a convolutional neural network (CNN). An ensemble numerical simulation-based training data set have been generated through dynamic simulations performed on the Polish system with various types of faults in different locations. Such data augmentation is proven to be able to provide adequate data for deep learning. The synchronous generators’ frequency signals are used and encoded into images for developing and evaluating CNN models for classification of fault types and locations. With a time-domain stacked image set as the benchmark, two different time-series encoding approaches, i.e., wavelet decomposition-based frequency-domain stacking and polar coordinate system-based Gramian Angular Field (GAF) stacking, are also adopted to evaluate and compare the CNN model performance and applicability. The various encoding approaches are suitable for different fault types and spatial zonation. With optimized settings of the developed CNN models, the classification and localization accuracies can go beyond 84 and 91%, respectively.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗

A Low Order and a High Order Solution for a Converging-Diverging Nozzle Problem

The Computational Aeroacoustics Workshops on Benchmark problems are conducted in order to generate solutions with a variety of methods for problems that are chosen both to be representative of significant and relevant acoustic phenomena. It is generally recognized that CAA applications require the resolution in time and in space of solution details over a significantly broader range of scales than typical computational fluid dynamics applications. These requirements are forcing the development of new methodologies for CAA applications. This paper presents numerical results from a third and a seventh order algorithm for the propagation of an acoustic signal through a converging-diverging nozzle. These two algorithms are members of a new family of high accuracy methods that have the same order accuracy in both space and time, and are an extension of previous work for linearized Euler equations to fully nonlinear time dependent problems. The simulations are all with the fully nonlinear quasi-1D Euler equations for the total solution, which includes both the fluid dynamics and the acoustics. The acoustic solution is obtained from the time dependent nonlinear solution by subtracting the steady solution.

Goodrich, John W.↗

Modeling the Benchmark Active Control Technology Wind-Tunnel Model for Application to Flutter Suppression

This paper describes the formulation of a model of the dynamic behavior of the Benchmark Active Controls Technology (BACT) wind-tunnel model for application to design and analysis of flutter suppression controllers. The model is formed by combining the equations of motion for the BACT wind-tunnel model with actuator models and a model of wind-tunnel turbulence. The primary focus of this paper is the development of the equations of motion from first principles using Lagrange's equations and the principle of virtual work. A numerical form of the model is generated using values for parameters obtained from both experiment and analysis. A unique aspect of the BACT wind-tunnel model is that it has upper- and lower-surface spoilers for active control. Comparisons with experimental frequency responses and other data show excellent agreement and suggest that simple coefficient-based aerodynamics are sufficient to accurately characterize the aeroelastic response of the BACT wind-tunnel model. The equations of motion developed herein have been used to assist the design and analysis of a number of flutter suppression controllers that have been successfully implemented.

Waszak, Martin R.↗

Benchmarking the New RESRAD-OFFSITE Source Term Model with DUST-MS and GoldSim - 13377

RESRAD-OFFSITE is a computer code developed by Argonne National Laboratory under the sponsorship of U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC). It is designed on the basis of RESRAD (onsite) code, a computer code designated by DOE and NRC for evaluating soil-contaminated sites for compliance with human health protection requirements pertaining to license termination or environmental remediation. RESRAD-OFFSITE has enhanced capabilities of modeling radionuclide transport to offsite locations and calculating potential radiation exposure to offsite receptors. Recently, a new source term model was incorporated into RESRAD-OFFSITE to enhance its capability further. This new source term model allows simulation of radionuclide releases from different waste forms, in addition to the soil sources originally considered in RESRAD (onsite) and RESRAD-OFFSITE codes. With this new source term model, a variety of applications can be achieved by using RESRAD-OFFSITE, including but not limited to, assessing the performance of radioactive waste disposal facilities. This paper presents the comparison of radionuclide release rates calculated by the new source term model of RESRAD-OFFSITE versus those calculated by DUST-MS and GoldSim, respectively. The focus of comparison is on the release rates of radionuclides from the bottom of the contaminated zone that was assumed to contain radioactive source materials buried in soil. The transport of released contaminants outside of the primary contaminated zone is beyond the scope of this paper. Overall, the agreement between the RESRAD-OFFSITE results and the DUST-MS and GoldSim results is fairly good, with all three codes predicting identical or similar radionuclide release profiles over time. Numerical dispersion in the DUST-MS and GoldSim results was identified as potentially contributing to the disagreement in the release rates. In general, greater discrepancy in the release rates was found for short-lived, fast-moving radionuclides than for long-lived, slow-moving radionuclides. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗

On the Approximability of Random-Hypergraph MAX-3-XORSAT Problems with Quantum Algorithms

Constraint satisfaction problems are an important area of computer science. Many of these problems are in the complexity class NP which is exponentially hard for all known methods, both for worst cases and often typical. Fundamentally, the lack of any guided local minimum escape method ensures the hardness of both exact and approximate optimization classically, but the intuitive mechanism for approximation hardness in quantum algorithms based on Hamiltonian time evolution is poorly understood. We explore this question using the prototypically hard MAX-3-XORSAT problem class. We conclude that the mechanisms for quantum exact and approximation hardness are fundamentally distinct. We qualitatively identify why traditional methods such as quantum adiabatic optimization are not good approximation algorithms. We propose a new spectral folding optimization method that does not suffer from these issues and study it analytically and numerically. We consider random rank-3 hypergraphs including extremal planted solution instances, where the ground state satisfies an anomalously high fraction of constraints compared to truly random problems. We show that, if we define the energy to be $E = N_{unsat}-N_{sat}$, then spectrally folded quantum optimization will return states with energy $E \leq A E_{GS}$ (where $E_{GS}$ is the ground state energy) in polynomial time, where conservatively, $A \simeq 0.6$. We thoroughly benchmark variations of spectrally folded quantum optimization for random classically approximation-hard (planted solution) instances in simulation, and find performance consistent with this prediction. We do not claim that this approximation guarantee holds for all possible hypergraphs, though our algorithm's mechanism can likely generalize widely. These results suggest that quantum computers are more powerful for approximate optimization than had been previously assumed.

Kapit, Eliot↗

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical Investigation of Radiative Heat Transfer in Laser Induced Air Plasmas

Radiative heat transfer is one of the most important phenomena in the laser induced plasmas. This study is intended to develop accurate and efficient methods for predicting laser radiation absorption and plasma radiative heat transfer, and investigate the plasma radiation effects in laser propelled vehicles. To model laser radiation absorption, a ray tracing method along with the Beer's law is adopted. To solve the radiative transfer equation in the air plasmas, the discrete transfer method (DTM) is selected and explained. The air plasma radiative properties are predicted by the LORAN code. To validate the present nonequilibrium radiation model, several benchmark problems are examined and the present results are found to match the available solutions. To investigate the effects of plasma radiation in laser propelled vehicles, the present radiation code is coupled into a plasma aerodynamics code and a selected problem is considered. Comparisons of results at different cases show that plasma radiation plays a role of cooling plasma and it lowers the plasma temperature by about 10%. This change in temperature also results in a reduction of the coupling coefficient by about 10-20%. The present study indicates that plasma radiation modeling is very important for accurate modeling of aerodynamics in a laser propelled vehicle.

Liu, J.↗

Performance of fully-coupled algebraic multigrid preconditioners for large-scale VMS resistive MHD

Here, we explore the current performance and scaling of a fully-implicit stabilized unstructured finite element (FE) variational multiscale (VMS) capability for large-scale simulations of 3D incompressible resistive magnetohydrodynamics (MHD). The large-scale linear systems that are generated by a Newton nonlinear solver approach are iteratively solved by preconditioned Krylov subspace methods. The efficiency of this approach is critically dependent on the scalability and performance of the algebraic multigrid preconditioner. Our study considers the performance of the numerical methods as recently implemented in the second-generation Trilinos implementation that is 64-bit compliant and is not limited by the 32-bit global identifiers of the original Epetra-based Trilinos. The study presents representative results for a Poisson problem on 1.6 million cores of an IBM Blue Gene/Q platform to demonstrate very large-scale parallel execution. Additionally, results for a more challenging steady-state MHD generator and a transient solution of a benchmark MHD turbulence calculation for the full resistive MHD system are also presented. These results are obtained on up to 131,000 cores of a Cray XC40 and one million cores of a BG/Q system.

97 MATHEMATICS AND COMPUTING↗

Benchmark experiment for electron-impact ionization of argon: Absolute triple-differential cross sections via three-dimensional electron emission images

Single ionization of argon by 195-eV electron impact is studied in an experiment, where the absolute triple-differential cross sections are presented as three-dimensional electron emission images for a series of kinematic conditions. Thereby a comprehensive set of experimental data for electron-impact ionization of a many-electron system is produced to provide a benchmark for comparison with theoretical predictions. Theoretical models using a hybrid first-order and second-order distorted-wave Born plus R-matrix approach are employed to compare their predictions with the experimental data. While the relative shape of the calculated cross section is generally in reasonable agreement with experiment, the magnitude appears to be the most significant problem with the theoretical treatment for the conditions studied in the present work. This suggests that the most significant challenge in the further development of theory for this process may lie in the reproduction of the absolute scale rather than the angular dependence of the cross section.

74 ATOMIC AND MOLECULAR PHYSICS↗

Benchmarking Quantum Simulators Using Ergodic Quantum Dynamics

Here we propose and analyze a sample-efficient protocol to estimate the fidelity between an experimentally prepared state and an ideal target state, applicable to a wide class of analog quantum simulators without advanced spatiotemporal control. Our protocol relies on universal fluctuations emerging from generic Hamiltonian dynamics, which we discover in the present work. It does not require fine-tuned control over state preparation, quantum evolution, or readout capability, while achieving near optimal sample complexity: a percent-level precision is obtained with ~103 measurements, independent of system size. Furthermore, the accuracy of our fidelity estimation improves exponentially with increasing system size. We numerically demonstrate our protocol in a variety of quantum simulator platforms, including quantum gas microscopes, trapped ions, and Rydberg atom arrays. We discuss applications of our method for tasks such as multiparameter estimation of quantum states and processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evolutionary Optimization of Quadrifilar Helical and Yagi-Uda Antennas

We present optimization results obtained for two type of antennas using evolutionary algorithms. A quadrifilar helical UHF antenna is currently flying aboard NASA's Mars Odyssey spacecraft and is due to reach final Martian orbit insertion in January, 2002. Using this antenna as a benchmark, we ran experiments employing a coevolutionary genetic algorithm to evolve a quadrifilar helical design in-situ - i.e., in the presence of a surrounding structure. Results show a 93% improvement at 400 MHz and a 48% improvement at 438 MHz in the average gain. The evolved antenna is also one-fourth the size. Yagi-Uda antennas are known to be difficult to design and optimize due to their sensitivity at high gain and the inclusion of numerous parasitic elements. Our fitness calculation allows the implicit relationship between power gain and sidelobe/backlobe loss to emerge naturally, a technique that is less complex than previous approaches. Our results include Yagi-Uda antennas that have excellent bandwidth and gain properties with very good impedance characteristics. Results exceeded previous Yagi-Uda antennas produced via evolutionary algorithms by at least 7.8% in mainlobe gain.

Lohn, Jason D.↗

The Effects of Different Lattice Codes on Results of Sensitivity and Uncertainty Analysis

According to the UAM ('Uncertainty Analysis in Modeling') benchmark, uncertainties in the LWR system calculations at all stages should be determined and quantified. As the neutronics calculations are prerequisite for the predictions of the reactor system, and the uncertainties introduced by the neutronics calculations would be propagated to the subsequent calculations, hence uncertainty analysis has been firstly focused on the neutronics phase to quantify the uncertainties of the reactor-physics predictions in UAM. For the lattice physics, the objective of uncertainty analysis is focused on uncertainty quantifications of the few-group constants; for the core physics, uncertainties of the core steady-state stand-alone neutronics calculations are focused on. Based on the previous researches, the nuclear-data uncertainty has been treated as one of the most significant source of uncertainty and received the most focus. In the previous works, the domestic code UNICORN has the capability of performing sensitivity analysis and uncertainty analysis using the direct numerical perturbation (DNP) method and statistical sampling (SS) method. In this paper, the capability of the UNICORN code has been improved. The domestic lattice code NECP-CACTI has been added in UNICORN to perform the lattice calculations, with DRAGON5.0 and NECP-CACTI available for lattice calculations. For DRAGON5.0 and NECP-CACTI, all the methods and technologies are shared in UNICORN. Therefore, the effects of different lattice codes on the results of sensitivity and uncertainty analysis have been researched and determined in this paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Langevin model of low-energy fission

Since the earliest days of fission, stochastic models have been used to describe and model the process. For a quarter century, numerical solutions of Langevin equations have been used to model fission of highly excited nuclei, where microscopic potential-energy effects have been neglected. In this paper I present a Langevin model for the fission of nuclei with low to medium excitation energies, for which microscopic effects in the potential energy cannot be ignored. I solve Langevin equations in a five-dimensional space of nuclear deformations. The macroscopic-microscopic potential energy from a global nuclear structure model well benchmarked to nuclear masses is tabulated on a mesh of approximately 10 7 points in this deformation space. The potential is defined continuously inside the mesh boundaries by use of a moving five-dimensional cubic spline approximation. Because of reflection symmetry, the effective mesh is nearly twice this size. For the inertia, I use a (possibly scaled) approximation to the inertia tensor defined by irrotational flow. A phenomenological dissipation tensor related to one-body dissipation is used. A normal-mode analysis of the dynamical system at the saddle point and the assumption of quasiequilibrium provide distributions of initial conditions appropriate to low excitation energies, and are extended to model spontaneous fission. A dynamical model of postscission fragment motion including dynamical deformations and separation allows the calculation of final mass and kinetic-energy distributions, along with other interesting quantities. The model makes quantitative predictions for fragment mass and kinetic-energy yields, some of which are very close to measured ones. Varying the energy of the incident neutron for induced fission allows the prediction of energy dependencies of fragment yields and average kinetic energies. With a simple approximation for spontaneous fission starting conditions, quantitative predictions are made for some observables which are close to measurements. In conclusion, this model is able to reproduce several mass and energy yield observables with a small number of physical parameters, some of which do not need to be varied after benchmarking to 235 U (n, f) to predict results for other fissioning isotopes.

233U↗

Ion Trapping Studies and Mitigation Strategies for the EIC ERL-Based Strong Hadron Cooler

An Energy Recovery Linac based strong hadron cooler was previously considered for the Electron-Ion Collider. The required electron beam parameters for variable-energy strong hadron cooling place significant constraints on ion trapping and collective effects. This paper presents initial studies of these constraints through a combination of analytical modelling and numerical simulations of ion production, trapping behaviour, and mitigation strategies. A multi-bunch tracking framework based on ELEGANT with the ionEffects module is used to simulate machine operation over millisecond time scales, corresponding to more than 3 × 10^5 electron bunches. The simulations include modelling of ionisation processes together with transverse electron–ion dynamics, allowing the evolution and accumulation of ions to be investigated. Analytical expressions based on Gaussian beam distributions are used to estimate ion trapping conditions and benchmark the simulation results. A bi-periodic bunch spacing scheme is also investigated as a possible mitigation method by detuning the ion oscillation frequency. These studies provide an initial assessment of ion trapping in the strong hadron cooler and demonstrate possible approaches for reducing beam–ion effects.

Bi, R. [Lancaster University, Cockcroft Institute]↗

Sulfur Pellets Responses to a Bare and Steel Reflected Pulse of the Oak Ridge National Laboratory Health Physics Research Reactor

The experiments analyzed in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the $\textit{Fast Burst Reactor}$. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. The HPRR was an unmoderated, unshielded fast reactor that used highly enriched uranium and molybdenum alloy as fuel. The reactor was initially sent to the Nevada Test Site in 1962, where it was used to evaluate radiation doses received as a result of the Hiroshima and Nagasaki bombings during World War II. A few years later, the reactor was sent back to ORNL to be part of the Dosimetry Application Research (DOSAR) facility shown in Figure 1, which included a reactor building shown on the left (west) of the picture and a control and laboratory building in the upper right corner (northeast). The critical assembly was used for numerous technical studies, including systems calibration, dosimetry, radiobiology of plants and animals, testing of radiation alarms, as well as teaching and training in radiation dosimetry and nuclear engineering. Between 1963 and 1987, the HPRR was operated for thousands of hours, achieving criticality close to 10,000 times and motivating many publications. The HPRR was decommissioned in 1987. The goal of this effort was to use historical data from operation of the HPRR to create a criticality accident alarm system (CAAS) benchmark to be included in the $\textit{International Handbook of Evaluated Criticality Safety Benchmark Experiments}$ (ICSBEP Handbook). A thorough inspection was performed of all available documentation and information available. The most promising experiments that were selected for evaluation were those described in the 1987 ORNL report entitled $\textit{Health Physics Research Reactor Reference Dosimetry}$, ORNL-6240. The report includes reference dosimetry results of the shielded and unshielded configurations of the HPRR after burst operations. Because of changes to the reactor positioning and storage systems that were made in 1985, the previous dosimetry reports became obsolete, and the newly designed experiments were needed to create the HPRR’s adjusted dosimetry data. The various results reported in ORNL-6240 include reference doses and dose equivalents from different conventions at different distances and elevations as determined using the detected neutron fluence and conversion factors. The HPRR neutron fluence was obtained through different methods, including sulfur pellet analysis and threshold detector unit data. Information about the HPRR spectrum was also obtained through Bonner sphere measurements. This benchmark is focused on a part of the measured sulfur fluences reported in Appendix H of ORNL-6240. Standard commercial sulfur pellets were placed at different distances from the HPRR centerline during burst operation and were activated due to the 32 S(n,p) 32 P reaction. The resulting 32 P activity was then measured and the information about the corresponding sulfur fluence and/or neutron dose could be extracted. Many of those measurements have 7 been performed with the HPRR in its bare configuration or with different shields (combinations of Lucite, concrete, steel). All the necessary, precise information about material and/or dimensions of the different shields was not found, so it was decided to focus only on the unshielded and steel-shielded configurations to minimize the benchmark uncertainty. A total of 31 cases (24 unshielded and 7 shielded cases at different positions) of sulfur fluence were selected before evaluation to develop the benchmark.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mean Estimate Distances for Galaxies with Multiple Estimates in NED-D

Numerous research topics rely on an improved cosmic distance scale (e.g., cosmology, gravitational waves) and the NASA/IPAC Extragalactic Database of Distances (NED-D) supports those efforts by tabulating multiple redshift-independent distances for 12,000 galaxies (e.g., Large Magellanic Cloud (LMC) zero-point). Six methods for securing a mean estimate distance (MED) from the data are presented (e.g., indicator and Decision Tree). All six MEDs yield surprisingly consistent distances for the cases examined, including for the key benchmark LMC and M106 galaxies. The results underscore the utility of the NED-D MEDs in bolstering the cosmic distance scale and facilitating the identification of systematic trends.

79 ASTRONOMY AND ASTROPHYSICS↗