Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code comparisons”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

MCNP® Code Version 6.3.1: Verification & Validation Testing (Rev.1)

This report describes the verification and validation testing performed on MCNP® code version 6.3.1. The purpose of this report is to act as a compendium of test suite descriptions and results. Accordingly the document is divided into two parts. Part I (this part) describes each test suite in the following sections. Part II provides the results of testing each suite and comparisons to experimental and/or alternative computational results, as appropriate.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

Abstract The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s −1 and diffusion coefficient of 0.05 m 2 s −1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: actuator design for particle and divertor condition control *

Abstract This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1 ⋅ 10 − 3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Perturbative model for the saturation of energetic-particle-driven modes limited by self-generated zonal modes

We present a simplified energy-conserving approach to incorporate wave–wave nonlinear effects within the framework commonly used to describe wave–particle nonlinearities. In particular, the effects of zonal mode (ZM) generation on the determination of the saturation amplitude of energetic particle (EP)-driven Alfvénic instabilities is studied. The model assumes that the zonal perturbations grow at a rate twice that of the original (pump) wave, consistent with a beat-driven (or force-driven) generation mechanism. The evolution and saturation of the mode amplitude are investigated both analytically and numerically within our reduced model assumptions, in both the collisionless and scattering-dominated regimes. These studies underscore the crucial role of sources and sinks in capturing the impact and the role of beat-driven zonal perturbations on mode evolution. In the realistic case of saturation set by sources and sinks, we discuss the role of a finite amplitude ZM in reducing microturbulent particle scattering, thus limiting the energy source for the resonant mode. We then discuss comparisons between the model’s predictions and simulation results. The model reproduces key features observed in gyrokinetic simulations as the reduction in saturated mode amplitude and the onset of wave–wave nonlinear effects as functions of mode growth rate and amplitude. Thanks to its simplicity, it can be readily implemented into codes based on reduced models, thereby improving their predictive capability for strongly driven instabilities.

Alfvén eigenmodes↗
DOE OSTI2022

Spiner-EOSPAC Comparison: performance and accuracy on Power9 CPU and GPU

The goal of this report is to concisely compare EOSPAC to Spiner EOS for a selection of commonly-used materials over as much of their domain as is possible. Importantly, this comparison does not address the way that the Euler equations may amplify interpolation errors and we take as given that the rational function interpolator gives the best possible representation of the EOS represented by the Sesame data. Neither does the report address all, or even majority of, the capabilities of each code nor does it explore performance or accuracy over multiple compilers and platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Evaluation and Comparison of a New Robust Waveform Against Direct Sequence Spread Spectrum for HF

This paper provides a comprehensive performance comparison between a filter bank multicarrier spread-spectrum (FBMC-SS) waveform and a current robust military waveform; namely, MIL-STD-188-110D, Waveform 0, proposed for communications through ionospheric/skywave HF channels. Waveform 0 is effectively a direct sequence spread spectrum waveform that uses the Walsh multi-codes to enhance the information transmission rate. It may thus be referred to as Walsh-DSSS. FBMC-SS, on the other hand, makes use of filter banks to provide excellent performance when the received signal is subject to partial band interference. Successful application of FBMC-SS for communications across skywave HF channels has been previously demonstrated, both theoretically and through experimental work. However, very little has been done to contrast FBMC-SS against Walsh-DSSS. The goal of this paper is to first add new features to FBMC-SS to bring it on par with Walsh- DSSS. These features include: (i), introduction of multi-codes that achieve a comparable (or better) data rate to the Walsh-DSSS; and (ii), addition of a scrambling step applied to the multi-codes to make the receiver detection robust against widely spread multipaths. With this established, in the second part of the paper, we examine the performance of the developed FBMC-SS against Walsh-DSSS when both are applied for communications across skywave HF channels. The two waveforms are compared both through a theoretical study and through experimental works across several skywave channels ranging from hundreds to thousands of kilometers.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

A scalable compressible volume of fluid solver using a stratified flow model

Here, in this paper, a Volume of Fluid method which uses a Stratified Flow model for flux calculation is proposed to solve a compressible multiphase flow in a way that has high parallel efficiency. The solver is part of an open-source computing tool to study fuel entrainment and combustion in the application of paraffin-based hybrid rocket motors. This paper focuses on the compressible Volume of Fluid solver for the liquid and gaseous phases in this environment. This solver uses perfect gas and stiffened gas models for equation of state, with the ability to easily add other models as desired. Additionally a comparison of a Riemann solver versus an AUSM+up scheme is conducted, showing improvement with the later in the solution of multiphase flows. This solver successfully reproduces shock tube cases with the expected accuracy, as well as simulating two-dimensional shear and gravity-driven flows. Results also show that the code is able to utilize a Stratified Flow model to evaluate shear flow without the need for interface reconstruction or gradient calculations. Finally, it is shown that the solver has near-ideal scaling under strong scaling tests as well as good performance in static scaling, giving improved performance over current options and implying future advancement of high-performance computing efficiency for multiphase flow solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2021

Determination of Seven Organic Impurities in FD&C Yellow No. 6 by Ultra-High-Performance Liquid Chromatography

The U.S. Food and Drug Administration (FDA) batch certifies FD&C Yellow No. 6 (Y6) to ensure that the color additive meets requirements published in the Code of Federal Regulations (CFR), including specifications for seven organic manufacturing impurities consisting of two intermediates, a reaction by-product, an impurity originating from an intermediate, and three subsidiary colors. An ultra-high-performance liquid chromatography (UHPLC) method was developed and validated for determining seven organic impurities in Y6. For comparison, the currently used high-performance liquid chromatography (HPLC) method was also validated. The new UHPLC method uses a 1.7 µm particle size biphenyl column with aqueous ammonium formate and methanol as eluants. Analytes are identified by comparing their retention times and UV-visible spectra to those of reference standards. Calibration is performed in the presence of the dye matrix and analyte levels are determined from their peak areas. UHPLC and HPLC validation studies obtained linear calibration curves and excellent values for limits of detection, limits of quantitation, recovery, repeatability, and intermediate precision for all analytes. Survey analyses of 30 samples using the UHPLC and HPLC methods yielded results that were consistent within experimental error; however, the UHPLC method provided more accurate results for two analytes that coelute using the HPLC method. The UHPLC method satisfies the accuracy and precision requisites for routine certification of Y6. The UHPLC method is faster and generates less waste than the HPLC method while delivering better sensitivity, separation, and accuracy.

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

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗
DOE OSTI2024

A full-scope, high-fidelity simulator-based hardware-in-the-loop testbed for comprehensive nuclear power plant cybersecurity research

Nuclear power plant (NPP) cybersecurity research often relies on hardware-in-the-loop (HIL) testbeds that integrate real hardware components into simulated environments. These testbeds allow researchers to identify vulnerabilities, evaluate attack impacts, and test security measures in a controlled setting. Furthermore, previous HIL testbeds lacked fidelity to accurately represent real nuclear systems, limiting the scope of cybersecurity analysis. This study presents the creation of a HIL testbed, devised upon a full-scope, high-fidelity NPP simulator, to facilitate realistic and comprehensive cybersecurity research. To demonstrate its capabilities, the control logic for the steam generator water level was migrated from the simulator to an external programmable logic controller. As a practical application of the developed testbed, supply chain attack scenarios were simulated by injecting malicious code into the controller logic, and the effects of manipulating sensor inputs and control commands were observed. While this HIL testbed provides more detailed simulations, enhanced realism, and wider applicability compared to other options utilizing a less complex simulator, it is also more intricate and costly. For this reason, we include a detailed comparison with some alternative architectures to aid fellow researchers and practitioners in the selection of a suitable HIL architecture based on specific research objectives.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

tonalli : an asexual genetic code to characterize APOGEE-2 stellar spectra. I. Validation with synthetic and solar spectra

ABSTRACT We present tonalli, a spectroscopic analysis python code that efficiently predicts effective temperature, stellar surface gravity, metallicity, $\alpha$-element abundance, and rotational and radial velocities for stars with effective temperatures between 3200 and 6250 K, observed with the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2). tonalli implements an asexual genetic algorithm to optimize the finding of the best comparison between a target spectrum and the continuum-normalized synthetic spectra library from the Model Atmospheres with a Radiative and Convective Scheme (MARCS), which is interpolated in each generation. Using simulated observed spectra and the APOGEE-2 solar spectrum of Vesta, we study the performance, limitations, accuracy, and precision of our tool. Finally, a Monte Carlo realization was implemented to estimate the uncertainties of each derived stellar parameter.

Adame, Lucía (ORCID:0000000263286099)↗
DOE OSTI2022

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗
DOE OSTI2024

Using Likwid and Byfl to Benchmark Hardware Performance

This paper outlines a benchmarking study conducted during my internship at LANL, focusing on CPU (Computer Processing Unit) and program performance assessment. The primary goal was to gather memory access data using three methods across five polybench kernels The data gathered would then be used to compare and contrast to one another and calculate operational intensity for performance comparisons. Benchmarking tools like Byfl and Likwid were employed, with Byfl offering hardware-independent data through LLVM compiler communication and Likwid directly interacting with computer hardware. The study considered various benchmarking factors, including optimization levels, Big O notation ((n)), CPU diversity and specific kernel equations. Big O notation was utilized to simplify code complexity, with detailed breakdwons of operations and memory components for each polybench application. Specific O(n) equations enabled nuanced kernel compariosns, facilitating the identification of performance variations. CPU efficiency assessments were conducted using Likwid tests on two CPUs. The central focus on code optimization aimed at achieving higher speeds and reduced memory usage through streamlined code. Future work propsoes creating a roofline model, synthesizing benchmarking data into a comprehensive data graph to assist in optimizing code and improving hardware performance. The potential impact on the laboratory or national mission was underscored, emphasizing the importance of optimizing applications and hardware to conserve resources and accelerate program execution. The specific relevance to LANL’s operations in math-intensive fields such as Nuclear Fission, Space Exploration, and Nanotechnology highlights the necessity of efficient benchmarking for resource conservation and proram speed. Overall, this study contributes to the understanding of CPU and program performance, providing insights for future optimization efforts in a laboratory setting

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Higher Order Multipole Analysis of Beam Crabbing Mode at 197 MHz for Electron-Ion Collider

In the local crabbing scheme for Electron Ion Collider interaction region, beam rotation is initiated by one set of crab cavities at a high beta region with near 90 deg phase advance to the interaction point (IP), and ceased by an identical set of cavities on the symmetrically opposite side. The fundamental mode of all crab cavities is designed to apply a momentum kick to the particle bunches to complete this half cycle of rotation in the horizontal plane, meanwhile the electromagnetic (EM) field can be expanded into a series of harmonics. Despite the lowest dipole mode, higher harmonics cause instabilities during beam operation and loss of particles. We used SimTrack, which is the 6-d symplectic element-by-element tracking code developed at Brookhaven National Laboratory in 2015, to study the threshold of each higher-order harmonic multipolar mode with the criteria of maintaining a dynamic aperture of 6 σ or greater for the Hadron Storage Ring (HSR). In this paper, we showed the threshold for each higher harmonic multipole and their comparison with the multipole strength for the current crab cavity design.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

CNRS Molten Salt Reactor Benchmark Analysis Using Griffin-Pronghorn Coupled Multi-physics Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time, since the fuel salt is circulating through the whole primary loop. This movement of fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and the corresponding redistribution of the DNPs in the active region of the core. To capture this effect, neutronics and thermal-hydraulics (T/H) codes need to be able to handle the movement of DNPs and their decay [1, 2]. Idaho National Laboratory (INL) is actively working on developing the neutronics and T/H codes to model the MSRs with flowing fuel. As one of the efforts, the neutronics code Griffin [3] and the T/H code Pronghorn [4] that are built upon the MOOSE framework [5] have been extended to handle the flowing fuel involving the drift of DNPs, and this capability is being tested and verified with MSR benchmark problems including the CNRS benchmark problem [6]. Additionally, a multi-physics analysis of the molten salt reactor experiment (MSRE) is being conducted to verify and validate the codes against available experimental data [7]. This work presents the analysis of the CNRS benchmark problem using the Griffin-Pronghorn coupled multi-physics code system of INL. Both steady-state and transient problems are analyzed, and the results are compared with those of other participants presented in [6]. The following section provides a description of the benchmark, along with input parameters and observables for each step. Then, the results for each step are presented in comparison with other participants’ results, followed by summary and conclusion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Evaluation of the Mechanical Properties of Cast and Wrought CF8C-Plus Relevant to ASME Code Case Qualification

The mechanical behavior of a cast form of an advanced austenitic stainless steel, CF8C-Plus, is compared with that of its wrought equivalent in terms of both tensile and creep-rupture properties and estimated allowable stress values for pressurized service at temperatures up to about 850°C. A traditional Larson-Miller parametric model is used to analyze the creep-rupture data and to predict long-term lifetimes for comparison of the two alloy types. The cast CF8C-Plus exhibited lower yield and tensile strengths, but higher creep strength compared to its wrought counterpart. Two welding methods, shielded-metal-arc welding (SMAW) and gas-metal-arc welding, met the weld qualification acceptance criteria in ASME BPVC Section IX for the cast CF8C-Plus. However, for the wrought CF8C-Plus, while SMAW and gas-tungsten-arc welding passed the tensile acceptance criteria, they failed the side bend tests due to lack of fusion or weld metal discontinuities.

Chen, Xiang↗