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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Development of iteration strategies for a practical implementation of a higher order transverse leakage approximation

Transverse integrated nodal diffusion methods currently represent the standard in full core neutronic simulation. The primary shortcoming in this approach is the utilization of the quadratic transverse leakage approximation. This approach, although proven to work well for typical LWR problems, is not consistent with the formulation of nodal methods and can cause accuracy and convergence problems. In previous work, an improved, consistent quadratic leakage approximation was formulated, which derived from the class of higher order nodal methods developed some years ago. In this paper a number of iteration schemes are developed around this consistent quadratic leakage approximation which yield accurate node average results in much improved calculational times. The developed consistent leakage approximation is extended in this work via a number of numerical schemes, the most promising of which results from utilizing the consistent leakage approximation as a correction method to the standard quadratic leakage approximation. Numerical results are demonstrated on a set of benchmark problems, such as the 3D IAEA LWR and MOX C5 problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Emergence of Novel Multipactor Modes Under Standing Wave Conditions in a Coaxial Line With an RF Window

The multipactor (MP) phenomenon can generate exponentially increasing electron populations in RF vacuum systems, leading to components becoming damaged. This paper uses numerical methods to analyze an RF vacuum window in a coaxial line under various standing wave conditions. Our methods are first benchmarked against the analytic and experimental results for the traveling wave and the numerical results for standing waves that were found in previous work. Second, the effect of standing waves on electron motion in the axial direction is discussed and related to electron trajectories in the presence of an RF window. It is found that standing waves with the magnetic maximum located behind the RF window lead to window collisions that can affect MP depending on the power, frequency, and conductor radii. For a low number of window collisions, MP was highly predictable by earlier work, but as the percentage of collisions occurring on the window surpassed approximately 50%, a novel mode of MP emerged to dominate the phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Parallel exponential time differencing methods for geophysical flow simulations

Two ocean models are considered for geophysical flow simulations: the multilayer shallow water equations and the multilayer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock–Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. Furthermore, these methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrates a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Convergence study of Rattlesnake solutions for the two-dimensional C5G7 MOX benchmark

This paper presents the convergence study of a specific transport scheme, self-adjoint angular flux (SAAF) formulation with the discrete ordinates (SN) method and continuous finite element method (CFEM), implemented with Rattlesnake, on solving the well known two-dimensional C5G7 MOX benchmark. Both the convergence in space and angle are studied. Numerical results show the convergence of the spatial and angular refinements.

Benchmark↗
DOE OSTI2022

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Discrete kinetic eigenmode spectra of electron plasma oscillations in weakly collisional plasma: A numerical study

It has been demonstrated that in the presence of weak collisions, described by the Lenard-Bernstein (LB) collision operator, the Landau-damped solutions become true eigenmodes of the system and constitute a complete set [C.-S. Ng et al., Phys. Rev. Lett. 83, 1974 (1999) and C. S. Ng et al., Phys. Rev. Lett. 96, 065002 (2004)]. We present numerical results from an Eulerian Vlasov code that incorporates the Lenard-Bernstein collision operator [A. Lenard and I. B. Bernstein, Phys. Rev. 112, 1456 (1958)]. The effect of collisions on the numerical recursion phenomenon seen in Vlasov codes is discussed. The code is benchmarked against exact linear eigenmode solutions in the presence of weak collisions, and a spectrum of Landau-damped solutions is determined within the limits of numerical resolution. Tests of the orthogonality and the completeness relation are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2013

Discrete Kinetic Eigenmode Spectra of Electron Plasma Oscillations in Weakly Collisional Plasma: A Numerical Study

It has been demonstrated that in the presence of weak collisions, described by the Lenard-Bernstein collision operator, the Landau-damped solutions become true eigenmodes of the system and constitute a complete set. We present numerical results from an Eulerian Vlasov code that incorporates the Lenard-Bernstein collision operator. The effect of the collisions on the numerical recursion phenomenon seen in Vlasov codes is discussed. The code is benchmarked against exact linear eigenmode solutions in the presence of weak collisions, and a spectrum of Landau-damped solutions is determined within the limits of numerical resolution. Tests of the orthogonality and the completeness relation are presented.

Kinetic Eigenmode Spectra↗
DOE OSTI2011

Benchmarking of the mono-energetic transport coefficients-results from the International Collaboration on Neoclassical Transport in Stellarators (ICNTS)

Numerical results for the three mono-energetic transport coefficients required for a complete neoclassical description of stellarator plasmas have been benchmarked within an international collaboration. These transport coefficients are flux-surface-averaged moments of solutions to the linearized drift kinetic equation which have been determined using field-line-integration techniques, Monte Carlo simulations, a variational method employing Fourier-Legendre test functions and a finite-difference scheme. The benchmarking has been successfully carried out for past, present and future devices which represent different optimization strategies within the extensive configuration space available to stellarators. A qualitative comparison of the results with theoretical expectations for simple model fields is provided. The behaviour of the results for the mono-energetic radial and parallel transport coefficients can be largely understood from such theoretical considerations but the mono-energetic bootstrap current coefficient exhibits characteristics which have not been predicted.

Beidler, C. D.↗
DOE OSTI2010

Power counting regime of chiral effective field theory and beyond.

Chiral effective field theory ({chi}EFT) complements numerical simulations of quantum chromodynamics (QCD) on a space-time lattice. It provides a model-independent formalism for connecting lattice simulation results at finite volume and a variety of quark masses to the physical world. The asymptotic nature of the chiral expansion places the focus on the first few terms of the expansion. Thus, knowledge of the power-counting regime (PCR) of {chi}EFT, where higher-order terms of the expansion may be regarded as negligible, is as important as knowledge of the expansion itself. Through the consideration of a variety of renormalization schemes and associated parameters, techniques to identify the PCR where results are independent of the renormalization scheme are established. The nucleon mass is considered as a benchmark for illustrating this general approach. Because the PCR is small, the numerical simulation results are also examined to search for the possible presence of an intrinsic scale which may be used in a nonperturbative manner to describe lattice simulation results outside of the PCR. Positive results that improve on the current optimistic application of chiral perturbation theory ({chi}PT) beyond the PCR are reported.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Large eddy simulation of a turbulent non-reacting spray jet

We performed Large Eddy Simulation (LES) of a turbulent non-reacting n-Heptane spray jet, referred to as Spray H in the Engine Combustion Network (ECN), and executed a data analysis focused on key LES metrics such as fraction of resolved turbulent kinetic energy and similarity index. In the simulation, we used the dynamic structure model for the sub-grid stress, and the Lagrangian-based spray-parcel models coupled with the blob-injection model. The finest mesh-cell size used was characterized by an Adaptive Mesh Refinement (AMR) cell size of 0.0625 mm. To obtain ensemble statistics, we performed 28 numerical realizations of the simulation. Demonstrated by the comparison with experimental data in a previous study [7], this LES has accurately predicted global quantities, such as liquid and vapor penetrations. The analysis in this work shows that 14 realizations of LES are sufficient to provide a reasonable representation of the average flow behavior that is benchmarked against the 28-realization ensemble. With the current mesh, numerical schemes, and sub-grid scale turbulence model, more than 95% of the turbulent kinetic energy is directly resolved in the flow regions of interest. The large-scale flow structures inferred from a statistical analysis reveal a region of disorganized flow around the peripheral region of the spray jet, which appears to be linked to the entrainment process.

Hu, Bing↗
DOE OSTI2025

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2010

Computational Fluid Dynamics Analysis of Very High Temperature Gas-Cooled Reactor Cavity Cooling System

The design of passive heat removal systems is one of the main concerns for the modular very high temperature gas-cooled reactors (VHTR) vessel cavity. The reactor cavity cooling system (RCCS) is a key heat removal system during normal and off-normal conditions. The design and validation of the RCCS is necessary to demonstrate that VHTRs can survive to the postulated accidents. The computational fluid dynamics (CFD) STAR-CCM+/V3.06.006 code was used for three-dimensional system modeling and analysis of the RCCS. A CFD model was developed to analyze heat exchange in the RCCS. The model incorporates a 180-deg section resembling the VHTR RCCS experimentally reproduced in a laboratory-scale test facility at Texas A&M University. All the key features of the experimental facility were taken into account during the numerical simulations. The objective of the present work was to benchmark CFD tools against experimental data addressing the behavior of the RCCS following accident conditions. Two cooling fluids (i.e., water and air) were considered to test the capability of maintaining the RCCS concrete walls' temperature below design limits. Different temperature profiles at the reactor pressure vessel (RPV) wall obtained from the experimental facility were used as boundary conditions in the numerical analyses to simulate VHTR transient evolution during accident scenarios. Mesh convergence was achieved with an intensive parametric study of the two different cooling configurations and selected boundary conditions. To test the effect of turbulence modeling on the RCCS heat exchange, predictions using several different turbulence models and near-wall treatments were evaluated and compared. The comparison among the different turbulence models analyzed showed satisfactory agreement for the temperature distribution inside the RCCS cavity medium and at the standpipes walls. For such a complicated geometry and flow conditions, the tested turbulence models demonstrated that the realizable k-epsilon model with two-layer all y+ wall treatment performs better than the other k-epsilon and k-omega turbulence models when compared to the experimental results and the Reynolds stress transport turbulence model results. A scaling analysis was developed to address the distortions introduced by the CFD model in simulating the physical phenomena inside the RCCS system with respect to the full plant configuration. The scaling analysis demonstrated that both the experimental facility and the CFD model achieve a satisfactory resemblance of the main flow characteristics inside the RCCS cavity region, and convection and radiation heat exchange phenomena are properly scaled from the actual plant.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Enhancement of Distribution System State Estimation Using Pruned Physics-Aware Neural Networks: Preprint

Realizing complete observability in the three-phase distribution system remains a challenge that hinders the implementation of classical state estimation algorithms. In this paper, a new method so-called pruned physics-aware neural network (P2N2) is developed to improve the voltage estimation accuracy in the distribution system. The method relies on the physical grid topology, which is used to design the connections between different hidden layers of a neural network model. To verify the proposed method, a numerical simulation based on one-year smart meter data of load consumptions for threephase power flow is developed to generate the measurement and voltage state data. The IEEE 123 node system is selected as the test network to benchmark the proposed algorithm against the classical weighted least squares (WLS). Numerical results show that P2N2 outperforms WLS, in terms of data redundancy and estimation accuracy.

distribution systems state estimation↗
DOE OSTI2021

Enhancement of Distribution System State Estimation Using Pruned Physics-Aware Neural Networks

Realizing complete observability in the three-phase distribution system remains a challenge that hinders the implementation of classic state estimation algorithms. In this paper, a new method, called the pruned physics-aware neural network (P2N2), is developed to improve the voltage estimation accuracy in the distribution system. The method relies on the physical grid topology, which is used to design the connections between different hidden layers of a neural network model. To verify the proposed method, a numerical simulation based on one- year smart meter data of load consumptions for three-phase power flow is developed to generate the measurement and voltage state data. The IEEE 123-node system is selected as the test network to benchmark the proposed algorithm against the classic weighted least squares (WLS). Numerical results show that P2N2 outperforms WLS in terms of data redundancy and estimation accuracy.

distribution system state estimation↗
DOE OSTI2019

A Leray regularized ensemble-proper orthogonal decomposition method for parameterized convection-dominated flows

Abstract Partial differential equations (PDEs) are often dependent on input quantities that are uncertain. To quantify this uncertainty PDEs must be solved over a large ensemble of parameters. Even for a single realization this can be a computationally intensive process. In the case of flows governed by the Navier–Stokes equations, an efficient method has been devised for computing an ensemble of solutions. To further reduce the computational cost of this method, an ensemble-proper orthogonal decomposition (POD) method was recently proposed. The main contribution of this work is the introduction of POD spatial filtering for ensemble-POD methods. The POD spatial filter makes possible the construction of the Leray ensemble-POD model, which is a regularized-reduced order model for the numerical simulation of convection-dominated flows of moderate Reynolds number. The Leray ensemble-POD model employs the POD spatial filter to smooth (regularize) the convection term in the Navier–Stokes equations, and diminishes the numerical inaccuracies produced by the ensemble-POD method in the numerical simulation of convection-dominated flows. Specifically, for the numerical simulation of a convection-dominated two-dimensional flow between two offset cylinders, we show that the Leray ensemble-POD method better reflects the dynamics of the benchmark results than the ensemble-POD scheme. The second contribution of this work is a new numerical discretization of the variable viscosity ensemble algorithm in which the average viscosity is replaced with the maximum viscosity. It is shown that this new numerical discretization is significantly more stable than those in current use. Furthermore, error estimates for the novel Leray ensemble-POD algorithm with this new numerical discretization are also proven.

Mathematics↗
NASA NTRS2005

Benchmark Problems Used to Assess Computational Aeroacoustics Codes

The field of computational aeroacoustics (CAA) encompasses numerical techniques for calculating all aspects of sound generation and propagation in air directly from fundamental governing equations. Aeroacoustic problems typically involve flow-generated noise, with and without the presence of a solid surface, and the propagation of the sound to a receiver far away from the noise source. It is a challenge to obtain accurate numerical solutions to these problems. The NASA Glenn Research Center has been at the forefront in developing and promoting the development of CAA techniques and methodologies for computing the noise generated by aircraft propulsion systems. To assess the technological advancement of CAA, Glenn, in cooperation with the Ohio Aerospace Institute and the AeroAcoustics Research Consortium, organized and hosted the Fourth CAA Workshop on Benchmark Problems. Participants from industry and academia from both the United States and abroad joined to present and discuss solutions to benchmark problems. These demonstrated technical progress ranging from the basic challenges to accurate CAA calculations to the solution of CAA problems of increasing complexity and difficulty. The results are documented in the proceedings of the workshop. Problems were solved in five categories. In three of the five categories, exact solutions were available for comparison with CAA results. A fourth category of problems representing sound generation from either a single airfoil or a blade row interacting with a gust (i.e., problems relevant to fan noise) had approximate analytical or completely numerical solutions. The fifth category of problems involved sound generation in a viscous flow. In this case, the CAA results were compared with experimental data.

Dahl, Milo D.↗
DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Application of Jacobian-free Newton–Krylov method in implicitly solving two-fluid six-equation two-phase flow problems: Implementation, validation and benchmark

This work represents a first-of-its-kind successful application to employ advanced numerical methods in solving realistic two-phase flow problems with two-fluid six-equation two-phase flow model. These advanced numerical methods include high-resolution spatial discretization scheme with staggered grids (high-order) fully implicit time integration schemes, and Jacobian-free Newton–Krylov (JFNK) method as the nonlinear solver. The computer code developed in this work has been extensively validated with existing experimental flow boiling data in vertical pipes and rod bundles, which cover wide ranges of experimental conditions, such as pressure, inlet mass flux, wall heat flux and exit void fraction. Additional code-to-code benchmark with the RELAP5-3D code further verifies the correct code implementation. The combined methods employed in this work exhibit strong robustness in solving two-phase flow problems even when phase appearance (boiling) and realistic discrete flow regimes are considered. Transitional flow regimes used in existing system analysis codes, normally introduced to overcome numerical difficulty, were completely removed in this work. As a result, this in turn provides the possibility to utilize more sophisticated flow regime maps in the future to further improve simulation accuracy.

42 ENGINEERING↗