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 829 records · Page 46

A Monte Carlo synthetic-acceleration method for solving the thermal radiation diffusion equation

We present a novel synthetic-acceleration-based Monte Carlo method for solving the equilibrium thermal radiation diffusion equation in three spatial dimensions. The algorithm performance is compared against traditional solution techniques using a Marshak benchmark problem and a more complex multiple material problem. Our results show that our Monte Carlo method is an effective solver for sparse matrix systems. For solutions converged to the same tolerance, it performs competitively with deterministic methods including preconditioned conjugate gradient and GMRES. We also discuss various aspects of preconditioning the method and its general applicability to broader classes of problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Low frequency eddy current benchmark study for model validation

This paper presents results of an eddy current model validation study. Precise measurements were made using an impedance analyzer to investigate changes in impedance due to Electrical Discharge Machining (EDM) notches in aluminum plates. Each plate contained one EDM notch at an angle of 0, 10, 20, or 30 degrees from the normal of the plate surface. Measurements were made with the eddy current probe both scanning parallel and perpendicular to the notch length. The experimental response from the vertical and oblique notches will be reported and compared to results from different numerical simulation codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radiation transport modeling using extended quadrature method of moments

The radiative transfer equation describes the propagation of radiation through a material medium. While it provides a highly accurate description of the radiation field, the large phase space on which the equation is defined makes it numerically challenging. As a consequence, significant effort has gone into the development of accurate approximation methods. Recently, an extended quadrature method of moments (EQMOM) has been developed to solve univariate population balance equations, which also have a large phase space and thus face similar computational challenges. The distinct advantage of the EQMOM approach over other moment methods is that it generates moment equations that are consistent with a positive phase space density and has a moment inversion algorithm that is fast and efficient. The goal of the current paper is to present the EQMOM method in the context of radiation transport, to discuss advantages and disadvantages, and to demonstrate its performance on a set of standard one-dimensional benchmark problems that encompass optically thin, thick, and transition regimes. Special attention is given in the implementation to the issue of realizability—that is, consistency with a positive phase space density. Numerical results in one dimension are promising and lay the foundation for extending the same framework to multiple dimensions.

97 MATHEMATICS AND COMPUTING↗

Precise Lifetime Measurements in Light Nuclei for Benchmarking Modern Ab-initio Nuclear Structure Models

A new generation of ab-initio calculations, based on realistic two- and three-body forces, is having a profound impact on our view of how nuclei work. To improve the numerical methods, and the parameterization of 3-body forces, new precise data are needed. Electromagnetic transitions are very sensitive to the dynamics which drive mixing between configurations. We have made a series of precise (< 3%) measurements of electromagnetic transitions in the A=10 nuclei {sup 10}C and {sup 10}Be by using the Doppler Shift Attenuation method carefully. Many interesting features can be reproduced including the strong α clustering. New measurements on {sup 8}Be and {sup 12}Be highlight the interplay between the alpha clusters and their valence neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

BIGHORN Computational Fluid Dynamics Theory, Methodology, and Code Verification & Validation Benchmark Problems

This document presents the theoretical background for a hybrid finite-element / finite-volume fluid flow solver, namely BIGHORN, based on the Multiphysics Object Oriented Simulation Environment (MOOSE) computational framework developed at the Idaho National Laboratory (INL). An overview of the numerical methods used in BIGHORN are discussed and followed by a presentation of the formulation details. The document begins with the governing equations for the compressible fluid flow, with an outline of the requisite constitutive relations. A second-order finite volume method used for solving the compressible fluid flow problems is presented next. A Pressure-Corrected Implicit Continuous-fluid Eulerian (PCICE) formulation for time integration is also presented. The multi-fluid formulation is being developed. Although multi-fluid is not fully-developed, BIGHORN has been designed to handle multi-fluid problems. Due to the flexibility in the underlying MOOSE framework, BIGHORN is quite extensible, and can accommodate both multi-species and multi-phase formulations. This document also presents a suite of verification & validation benchmark test problems for BIGHORN. The intent for this suite of problems is to provide baseline comparison data that demonstrates the performance of the BIGHORN solution methods on problems that vary in complexity from laminar to turbulent flows. Wherever possible, some form of solution verification has been attempted to identify sensitivities in the solution methods, and suggest best practices when using BIGHORN.

97 MATHEMATICS AND COMPUTING↗

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗

DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks

Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.

97 MATHEMATICS AND COMPUTING↗

Electromagnetic properties of indium isotopes illuminate the doubly magic character of 100 Sn

Understanding the nuclear properties in the vicinity of 100 Sn, which has been suggested to be the heaviest doubly magic nucleus with proton number Z equal to neutron number N, has been a long-standing challenge for experimental and theoretical nuclear physics. In particular, contradictory experimental evidence exists regarding the role of nuclear collectivity in this region of the nuclear chart. Here, we provide further evidence for the doubly magic character of 100 Sn by measuring the ground-state electromagnetic moments and nuclear charge radii of indium (Z = 49) isotopes as N approaches 50 from above using precision laser spectroscopy. Our results span almost the complete range between the two major closed neutron shells at N = 50 and N = 82 and reveal parabolic trends as a function of the neutron number, with a clear reduction towards these two closed neutron shells. Furthermore, a detailed comparison between our experimental results and numerical results from two complementary nuclear many-body frameworks (density functional theory and ab initio methods) exposes deficiencies in nuclear models and establishes a benchmark for future theoretical developments.

100Sn↗

Radiation from particles with arbitrary energy falling into higher-dimensional black holes

We consider point particles with arbitrary energy per unit mass E that fall radially into a higher-dimensional, nonrotating, asymptotically flat black hole. We compute the energy and linear momentum radiated in this process as functions of E and of the spacetime dimensionality D=n+2 for n=2,...,9 (in some cases we go up to 11). We find that the total energy radiated increases with n for particles falling from rest (E=1). For fixed particle energies 1<E{<=}2 we show explicitly that the radiation has a local minimum at some critical value of n, and then it increases with n. We conjecture that such a minimum exists also for higher particle energies. The present point-particle calculation breaks down when n=11, because then the radiated energy becomes larger than the particle mass. Quite interestingly, for n=11 the radiated energy predicted by our calculation would also violate Hawking's area bound. This hints at a qualitative change in gravitational radiation emission for n > or approx. 11. Our results are in very good agreement with numerical simulations of low-energy, unequal-mass black hole collisions in D=5 (that will be reported elsewhere) and they are a useful benchmark for future nonlinear evolutions of the higher-dimensional Einstein equations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Benchmarking a new closed-form thermal analysis technique against a traditional lumped parameter, finite-difference method

A benchmarking effort was conducted to determine the accuracy of a new analytic generic geology thermal repository model developed at LLNL relative to a more traditional, numerical, lumped parameter technique. The fast-running analytical thermal transport model assumes uniform thermal properties throughout a homogenous storage medium. Arrays of time-dependent heat sources are included geometrically as arrays of line segments and points. The solver uses a source-based linear superposition of closed form analytical functions from each contributing point or line to arrive at an estimate of the thermal evolution of a generic geologic repository. Temperature rise throughout the storage medium is computed as a linear superposition of temperature rises. It is modeled using the MathCAD mathematical engine and is parameterized to allow myriad gridded repository geometries and geologic characteristics [4]. It was anticipated that the accuracy and utility of the temperature field calculated with the LLNL analytical model would provide an accurate 'birds-eye' view in regions that are many tunnel radii away from actual storage units; i.e., at distances where tunnels and individual storage units could realistically be approximated as physical lines or points. However, geometrically explicit storage units, waste packages, tunnel walls and close-in rock are not included in the MathCAD model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigation of Microbunching Instabilities in Modern Recirculating Accelerators

Particle accelerators are machines to accelerate and store charged particles, such as electrons or protons, to the energy levels for various scientific applications. A collection of charged particles usually forms a particle beam. There are three basic types of particle accelerators: linear accelerators (linac), storage-ring (or circular) accelerators, and recirculating accelerators. In a linac, particles are accelerated and pass through once along a linear or straight beamline. Storage-ring accelerators propel particles around a circular track and repetitively append the energy to the stored beam. The third type, also the most recent one in chronology, the recirculating accelerator, is designed to accelerate the particle beam in a short section of linac, circulate the beam, and then either continue to accelerate for energy boost or decelerate it for energy recovery. The beam properties of a linac machine are set at best by the initial particle sources. For storage rings, the beam equilibria are instead determined by the overall machine design. The modern recirculating machines share with linacs the advantages to both accelerate and preserve the beam with high beam quality, as well as efficiently reuse the accelerating components. The beamline design in such a machine configuration can however be much more complicated than that of linacs. As modern accelerators push toward the high-brightness or high-intensity frontier by demanding particles in a highly charged bunch (about nano-Coulomb per bunch) to concentrate in an ever-decreasing beam phase space (transverse normalized emittance about 1 μm and relative energy spread of the order of 10^-5 in GeV beam energy), the interaction amongst particles via their self-generated electromagnetic fields can potentially lead to coherent instabilities of the beam and thus pose significant challenges to the machine design and operation. In the past decade and a half, microbunching instability (MBI) has been one of the most challenging issues for such high-brightness or high-intensity beam transport, as it would degrade lasing performance in the fourth-generation light sources, reduce cooling efficiency in electron cooling facilities, and compromise the luminosity of colliding beams in lepton or lepton-hadron colliders. The dissertation work will focus on the MBI in modern recirculating electron accelerators. It has been known that the collective interactions, the coherent synchrotron radiation (CSR) and the longitudinal space charge (LSC) forces, can drive MBI. The CSR effect is a collective phenomenon in which the electrons in a curved motion, e.g. a bending dipole, emit radiation at a scale comparable to the micro-bunched structure of the bunch distribution. The LSC effect stems from non-uniformity of the charge distribution, acts as plasma oscillation, and can eventually accumulate an amount of energy modulation when the beam traverses a long section of a beamline. MBI can be seeded by non-uniformity or shot noise of the beam, which originates from granularity of the elementary charge. Through the aforementioned collective effects, the modulation of the bunch sub-structure can be amplified and, once the beam-wave interaction formed a positive feedback, can result in MBI. The problem of MBI has been intensively studied for linac-based facilities and for storage-ring accelerators. However, systematic studies for recirculation machines are still very limited and form a knowledge gap. Because of the much more complicated machine configuration of the recirculating accelerators than that of linacs, the existing MBI analysis needs to be extended to accommodate the high-brightness particle beam transport in modern recirculating accelerators. This dissertation is focused on theoretical investigation of MBI in such machine configuration in the following seven themes: (1) Development and generalization of MBI theory The theoretical formulation has been extended so as to be applicable to a general linear beamline lattice including horizontal and vertical transport bending elements, and beam acceleration or deceleration. These featured generalizations are required for MBI analysis in recirculation accelerators. (2) Construction of CSR impedance models In addition to the steady-state CSR interaction, it has been found that the exit transient effect (or CSR drift) can even result in more serious MBI in high-brightness recirculation arcs. The onedimensional free-space CSR impedances, especially the exit transients, are derived. The steady-state CSR impedance is also extended to non-ultrarelativistic beam energy for MBI analysis of low-energy merger sections in recirculating accelerators. (3) Numerical implementation of the derived semi-analytical formulation This includes the development of a semi-analytical Vlasov solver for MBI analysis, and also benchmarking of the solver against massive particle tracking simulations. (4) Exploration of multistage amplification behavior of CSR microbunching development The CSR-induced MBI acts as an amplifier, which amplifies the sub-bunch modulation of a beam. The amplification is commonly quantified by the amplification gain. A beam transport system can be considered as a cascaded amplifier. Unlike the two-stage amplification of four-dipole bunch compressor chicanes employed in linacs, the recirculation arcs, which are usually constituted by several tens of bending magnets, show a distinguishing feature of up to six-stage microbunching amplification for our example arc lattices. That is, the maximal CSR amplification gain can be proportional to the peak bunch current up to sixth power. A method to compare lattice performance has been developed in terms of gain coefficients, which nearly depend on the lattice properties only. This method has also proven to be an effective way to quantify the current dependence of the maximal (5) Control of CSR MBI in multibend transport or recirculation arcs The existing mitigation schemes of MBI mostly aim to linac-based accelerators and may not be practical to the recirculating accelerator facilities. Thus a set of conditions for suppression of CSR MBI was proposed and examined for example lattices from low (~100 MeV) to high (~1 GeV) energies. (6) Study of more aspects of microbunched structures in beam phase spaces For a cascaded amplifier in circuit electronics, the total amplification gain can be estimated as the product of individual gains. In a beam transport line of an accelerator, the (scalar) gain multiplication was examined and found to under-estimate the overall microbunching amplification. The concept of gain matrix was developed, which includes the density, energy and transverse-longitudinal modulations in a beam phase space, and used to analyze MBI for a proposed recirculating machine. Throughout the gain matrix approach, it reasonably gives the upper limit of spectral MBI gain curves. This extended analysis can be employed to study multi-pass recirculation. (7) Study of MBI for magnetized beams Driven by a recent energy-recovery-linac based cooler design for electron cooling at Jefferson Lab Electron-Ion Collider Project, the generalized theoretical formulation for MBI to a transversely coupled beam has been developed and applied to this study. A magnetized beam in general features non-zero canonical angular momentum, thus considered to be a transversely coupled beam. A novel idea of utilizing magnetized beam transport was proposed for improvement of cooling efficiency and possible mitigation of collective effects. A concern of MBI regarding this design was studied and excluded. The large transverse beam size associated with the beam magnetization is found to help suppress MBI via the transverse-longitudinal correlation.

43 PARTICLE ACCELERATORS↗

Efficient and accurate numerical integration of exchange-correlation density functionals

We present a method for efficient calculation of exchange correlation in a density-functional theory calculation with Gaussian basis functions. It seamlessly combines the treatments of smooth (low-energy cutoff) and nonsmooth parts of the electron density such that computational performance can be maximized without approximation. Benchmark calculations show that the scheme can speed up the calculation of the exchange-correlation by as much as ten times or more without error.

74 ATOMIC AND MOLECULAR PHYSICS↗

Numerical Simulation of Textile Composite Stamping On Double Dome

Stamping is one of the most effective ways to form textile composites in industry for providing high-strength, low-weight and cost-effective products. This paper presents a fully continuum mechanics-based approach for stamping simulation of textile fiber reinforced composites by using finite element (FE) method. A previously developed non-orthogonal constitutive model is used to represent the anisotropic mechanical behavior of textile composites under large deformation during stamping. Simulation are performed on a balanced plain weave composite with 0 deg./90 deg. and {+-}45 deg. as initial yarn orientation over a benchmark double dome device. Simulation results show good agreement with experimental output in terms of a number of parameters selected for comparison.

36 MATERIALS SCIENCE↗

Performance Characterization for Fusion Co-design Applications

Magnetic fusion is a long-term solution for producing electrical power for the world, and the large thermonuclear international device (ITER) being constructed will produce net energy and a path to fusion energy provided the computer modeling is accurate. To effectively address the requirements of the high-end fusion simulation community, application developers, algorithm designers, and hardware architects must have reliable simulation data gathered at scale for scientifically valid configurations. This paper presents detailed benchmarking results for a set of magnetic fusion applications with a wide variety of underlying mathematical models including both particle-in-cell and Eulerian codes using both implicit and explicit numerical solvers. Our evaluation on a petascale Cray XE6 platform focuses on profiling these simulations at scale identifying critical performance characteristics, including scalability, memory/network bandwidth limitations, and communication overhead. Overall results are a key in improving fusion code design, and are a critical first step towards exascale hardware-software co-design — a process that tightly couples applications, algorithms, imple- mentation, and computer architecture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flux in a 1D-Reactor Cell by Response Matrix

Legacy computational reactor physics software tools and protocols currently used in support of the Advanced Test Reactor (ATR) core fuel management and safety assurance, and to some extent, experiment management, are inconsistent with the state of modern nuclear engineering practice. Thus, they are difficult, if not impossible, to verify and validate (V and V) according to modern standards. Furthermore, the legacy staff knowledge required for effective application of their application and protocols from the 1960's and 1970's is rapidly being lost due to staff turnover and retirements. In late 2009, the Idaho National Laboratory (INL) initiated a focused effort, the ATR Core Modeling Update Project, to address this situation through the introduction of modern high-fidelity transport theory based computational software and models to replace the legacy methods, with appropriate (V and V) protocols according to applicable national standards. The new computational tools include: - The 2-dimensional HELIOS deterministic transport code for integrated ATR core fuel and experiment management and Core Safety Analysis Package (CSAP) physics data generation. The MCNP5 3-dimensional stochastic simulation (Monte Carlo) code for verification and validation applications and for specialized ATR core analysis along with its longstanding traditional role in experiment support - The 2-dimensional deterministic NEWT transport code and the 3-dimensional KENO6 stochastic simulation code, with their SCALE-based support infrastructure for use in V and V applications and as supporting modules for sensitivity-uncertainty analysis via the TSUNAMI system. The ATTILA 3-dimensional discrete-ordinates transport code for specialized experiment support and possibly additional applications. Regardless of application, every neutron transport code, new or legacy, requires verification. The aim of this presentation is to provide a neutron transport benchmark for a simple 1D-reactor cell in plane geometry. The benchmark finds use as a standard of comparison for the codes listed above intended for the ATR upgrade. The method developed here is in contrast to the finite Fourier transform solution found in Ref. 1, which is limited to isotropic scattering and particularly difficult to numerically implement unless the total cross section is constant over the cell.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multidisciplinary benchmarks of a conservative spectral solver for the nonlinear Boltzmann equation

The Boltzmann equation describes the evolution of the phase-space probability distribution of classical particles under binary collisions. Approximations to it underlie the basis for several scholarly fields, including aerodynamics and plasma physics. While these approximations are appropriate in their respective domains, they can be violated in niche but diverse applications which require direct numerical solution of the original nonlinear Boltzmann equation. An expanded implementation of the Galerkin–Petrov conservative spectral algorithm is employed to study a wide variety of physical problems. Enabled by distributed precomputation, solutions of the spatially homogeneous Boltzmann equation can be achieved in seconds on modern personal hardware, while spatially-inhomogeneous problems are solvable in minutes. Here, several benchmarks are presented focusing on accuracy compared to both analytic theoretical predictions and other Boltzmann solvers. These benchmarks span several physical domains including weakly ionized plasma, gaseous fluids, and atomic-plasma interaction.

97 MATHEMATICS AND COMPUTING↗

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗

Thermal Transport by First-Principles Anharmonic Lattice Dynamics

Predictive methodologies for calculating lattice thermal transport and related properties have been recently developed and benchmarked successfully against measurements in a wide range of materials, from bulk to the nanoscale. This chapter presents an overview of the theoretical underpinnings of typical numerical methods for solving the Peierls-Boltzmann equation (PBE) for phonon transport when coupled with density functional theory (DFT). We discuss how these numerical methods may be employed to advance fundamental understanding of thermal transport, particularly as applied to measurements of anharmonic properties (e.g., linewidths, lattice expansion) in materials. In particular, calculations of interatomic force constants, integration methods, and construction of transport properties will be explained in detail with reference to relevant research in the literature. Furthermore, this chapter will provide an overview of applications and methods for calculations of thermal conductivities beyond perfect infinite crystals (including phonon-defect scattering) and of other measured observables, including mean free path accumulation, linewidths, and densities of states from scattering experiments.

Lindsay, Lucas↗