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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗
DOE OSTI2011

Petascale direct numerical simulation of turbulent combustion—fundamental insights towards predictive models

The advent of petascale computing applied to direct numerical simulation (DNS) of turbulent combustion has transformed our ability to interrogate fine-grained ‘turbulence-chemistry’ interactions in canonical and laboratory configurations. In particular, three-dimensional DNS, at moderate Reynolds numbers and with complex chemistry, is providing unprecedented levels of detail to isolate and reveal fundamental causal relationships between turbulence, mixing and reaction. This information is leading to new physical insight, providing benchmark data for assessing model assumptions, suggesting new closure hypotheses, and providing interpretation of statistics obtained from lower-dimensional measurements. In this paper the various roles of petascale DNS are illustrated through selected examples related to lifted flame stabilization, premixed and stratified flame propagation in intense turbulence, and extinction and reignition in turbulent non-premixed jet flames. Extending the DNS envelope to higher Reynolds numbers, higher pressures, and greater chemical complexity will require exascale computing in the next decade. The future outlook of DNS in terms of challenges and opportunities in this regard are addressed.

Energy & Fuels↗
DOE OSTI2021

A Big Problem for Small Earthquakes: Benchmarking Routine Magnitudes and Conversion Relationships with Coda Envelope-Derived Mw in Southern Kansas and Northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude (M w ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w <4) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w . Therefore, to assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L , m bLg , and M w for M ~2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ~3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b-values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b-values, should be treated with caution.

58 GEOSCIENCES↗
DOE OSTI2022

The 2022 solar fuels roadmap

Abstract Renewable fuel generation is essential for a low carbon footprint economy. Thus, over the last five decades, a significant effort has been dedicated towards increasing the performance of solar fuels generating devices. Specifically, the solar to hydrogen efficiency of photoelectrochemical cells has progressed steadily towards its fundamental limit, and the faradaic efficiency towards valuable products in CO 2 reduction systems has increased dramatically. However, there are still numerous scientific and engineering challenges that must be overcame in order to turn solar fuels into a viable technology. At the electrode and device level, the conversion efficiency, stability and products selectivity must be increased significantly. Meanwhile, these performance metrics must be maintained when scaling up devices and systems while maintaining an acceptable cost and carbon footprint. This roadmap surveys different aspects of this endeavor: system benchmarking, device scaling, various approaches for photoelectrodes design, materials discovery, and catalysis. Each of the sections in the roadmap focuses on a single topic, discussing the state of the art, the key challenges and advancements required to meet them. The roadmap can be used as a guide for researchers and funding agencies highlighting the most pressing needs of the field.

08 HYDROGEN↗
DOE OSTI2021

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗
NASA NTRS1983

Evaluation of the synoptic and mesoscale predictive capabilities of a mesoscale atmospheric simulation system

The overall performance characteristics of a limited area, hydrostatic, fine (52 km) mesh, primitive equation, numerical weather prediction model are determined in anticipation of satellite data assimilations with the model. The synoptic and mesoscale predictive capabilities of version 2.0 of this model, the Mesoscale Atmospheric Simulation System (MASS 2.0), were evaluated. The two part study is based on a sample of approximately thirty 12h and 24h forecasts of atmospheric flow patterns during spring and early summer. The synoptic scale evaluation results benchmark the performance of MASS 2.0 against that of an operational, synoptic scale weather prediction model, the Limited area Fine Mesh (LFM). The large sample allows for the calculation of statistically significant measures of forecast accuracy and the determination of systematic model errors. The synoptic scale benchmark is required before unsmoothed mesoscale forecast fields can be seriously considered.

Koch, S. E.↗
NASA NTRS2000

Numerical Study of Wave Propagation in a Non-Uniform Flow

The propagation of acoustic waves originating from cylindrical and spherical pulses, in a non-uniform mean flow, and in the presence of a reflecting wall is investigated by Hardin and Pope approach using compact approximation of spatial derivatives. The 2-D and 3-D stagnation flows and a flow around a cylinder are taken as prototypes of real world flows with strong gradients of mean pressure and velocity. The intensity and directivity of acoustic wave patterns appear to be quite different from the benchmark solutions obtained in a static environment for the same geometry. The physical reasons for amplification and weakening of sound are discussed in terms of dynamics of wave profile and redistribution of acoustic energy and its potential and kinetic components. For an acoustic wave in the flow around a cylinder, the observed mean acoustic pressure is approximately doubled (upstream pulse position) and halved (downstream pulse position) in comparison with the sound propagation in static ambient conditions.

Povitsky, Alex↗
DOE OSTI2024

Exact enforcement of temporal continuity in sequential physics-informed neural networks

The use of deep learning methods in scientific computing represents a potential paradigm shift in engineering problem solving. One of the most prominent developments is Physics-Informed Neural Networks (PINNs), in which neural networks are trained to satisfy partial differential equations (PDEs). While this method shows promise, the standard version has been shown to struggle in accurately predicting the dynamic behavior of time-dependent problems. To address this challenge, methods have been proposed that decompose the time domain into multiple segments, employing a distinct neural network in each segment and directly incorporating continuity between them in the loss function of the minimization problem. In this work we introduce a method to exactly enforce continuity between successive time segments via a solution ansatz. This hard constrained sequential PINN (HCS-PINN) method is simple to implement and eliminates the need for any loss terms associated with temporal continuity. The method is tested for a number of benchmark problems involving both linear and non-linear PDEs. Examples include various first order time dependent problems in which traditional PINNs struggle, namely advection, Allen–Cahn, and Korteweg–de Vries equations. Furthermore, second and third order time-dependent problems are demonstrated via wave and Jerky dynamics examples, respectively. Notably, the Jerky dynamics problem is chaotic, making the problem especially sensitive to temporal accuracy. Finally, the numerical experiments conducted with the proposed method demonstrated superior convergence and accuracy over both traditional PINNs and the soft-constrained counterparts.

42 ENGINEERING↗
DOE OSTI2017

MC 2 -3: Multigroup Cross Section Generation Code for Fast Reactor Analysis

This paper presents the methods and performance of the MC2 -3 code, which is a multigroup cross-section generation code for fast reactor analysis, developed to improve the resonance self-shielding and spectrum calculation methods of MC2 -2 and to simplify the current multistep schemes generating region-dependent broad-group cross sections. Using the basic neutron data from ENDF/B data files, MC2 -3 solves the consistent P1 multigroup transport equation to determine the fundamental mode spectra for use in generating multigroup neutron cross sections. A homogeneous medium or a heterogeneous slab or cylindrical unit cell problem is solved in ultrafine (2082) or hyperfine (~400 000) group levels. In the resolved resonance range, pointwise cross sections are reconstructed with Doppler broadening at specified temperatures. The pointwise cross sections are directly used in the hyperfine group calculation, whereas for the ultrafine group calculation, self-shielded cross sections are prepared by numerical integration of the pointwise cross sections based upon the narrow resonance approximation. For both the hyperfine and ultrafine group calculations, unresolved resonances are self-shielded using the analytic resonance integral method. The ultrafine group calculation can also be performed for a two-dimensional whole-core problem to generate region-dependent broad-group cross sections. Verification tests have been performed using the benchmark problems for various fast critical experiments including Los Alamos National Laboratory critical assemblies; Zero-Power Reactor, Zero-Power Physics Reactor, and Bundesamt für Strahlenschutz experiments; Monju start-up core; and Advanced Burner Test Reactor. Verification and validation results with ENDF/B-VII.0 data indicated that eigenvalues from MC2 -3/DIF3D agreed well with Monte Carlo N-Particle5 MCNP5 or VIM Monte Carlo solutions within 200 pcm and regionwise one-group fluxes were in good agreement with Monte Carlo solutions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1998

CFD Modelling of Bore Erosion in Two-Stage Light Gas Guns

A well-validated quasi-one-dimensional computational fluid dynamics (CFD) code for the analysis of the internal ballistics of two-stage light gas guns is modified to explicitly calculate the ablation of steel from the gun bore and the incorporation of the ablated wall material into the hydrogen working cas. The modified code is used to model 45 shots made with the NASA Ames 0.5 inch light gas gun over an extremely wide variety of gun operating conditions. Good agreement is found between the experimental and theoretical piston velocities (maximum errors of +/-2% to +/-6%) and maximum powder pressures (maximum errors of +/-10% with good igniters). Overall, the agreement between the experimental and numerically calculated gun erosion values (within a factor of 2) was judged to be reasonably good, considering the complexity of the processes modelled. Experimental muzzle velocities agree very well (maximum errors of 0.5-0.7 km/sec) with theoretical muzzle velocities calculated with loading of the hydrogen gas with the ablated barrel wall material. Comparison of results for pump tube volumes of 100%, 60% and 40% of an initial benchmark value show that, at the higher muzzle velocities, operation at 40% pump tube volume produces much lower hydrogen loading and gun erosion and substantially lower maximum pressures in the gun. Large muzzle velocity gains (2.4-5.4 km/sec) are predicted upon driving the gun harder (that is, upon using, higher powder loads and/or lower hydrogen fill pressures) when hydrogen loading is neglected; much smaller muzzle velocity gains (1.1-2.2 km/sec) are predicted when hydrogen loading is taken into account. These smaller predicted velocity gains agree well with those achieved in practice. CFD snapshots of the hydrogen mass fraction, density and pressure of the in-bore medium are presented for a very erosive shot.

Bogdanoff, D. W.↗
DOE OSTI2013

Solution of stochastic media transport problems using a numerical quadrature-based method

We present a new conceptual framework for analyzing transport problems in random media. We decompose such problems into stratified subproblems according to the number of material pseudo-interfaces within realizations. For a given subproblem we assign pseudo-interface locations in each realization according to product quadrature rules, which allows us to deterministically generate a fixed number of realizations. Quadrature integration of the solutions of these realizations thus approximately solves each subproblem; the weighted superposition of solutions of the subproblems approximately solves the general stochastic media transport problem. We revisit some benchmark problems to determine the accuracy and efficiency of this approach in comparison to randomly generated realizations. We find that this method is very accurate and fast when the number of pseudo-interfaces in a problem is generally low, but that these advantages quickly degrade as the number of pseudo-interfaces increases. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Developing a microwave-driven reactor for ammonia synthesis: insights into the unique challenges of microwave catalysis

Rapid development of new ammonia (NH 3 ) synthesis techniques that enable modular, intermittent production is essential to actualizing NH 3 's potential as a clean energy carrier, since contemporary methods are configured to centralized, continuous production methods with high emissions and are incompatible with renewable sources. Here, in this mission, microwave-driven catalysis is promising for its ability to enhance reaction kinetics and apply targeted heating for efficient energy use. However, owing to an incomplete understanding of the interaction between microwave fields and catalyst beds, the development of such microwave-catalysis systems remains underexplored and challenging. This paper investigates the 10× scale-up of a microwave-based NH 3 synthesis reactor using numerical and experimental approaches, achieving the largest reported microwave-driven NH 3 reactor to date. Results elucidate phenomena unique to microwave processes that are only predictable through numerical modeling, including how a catalyst's dielectric properties influence microwave field distribution by affecting penetration depth and how energy utilization can be poor even with sufficient attenuation. These dynamics change with scale, constrain reactor geometry, and potentially hamper performance. Nonetheless, we demonstrate a production rate of 56.6 g NH$_3$ per day, the highest reported NH 3 synthesis rate for laboratory-scale alternative techniques; correspondingly, the benchmark energy efficiency achieved in this paper (45.6 g NH$_3$ kW h -1 ) is the highest reported for such reactors of sufficient scale. Even with this exemplary energy efficiency, energy losses were found in excess of 50%, an issue resolvable through scale and reactor design. The efficiencies imparted by microwaves were key in these achievements, warranting further investigation toward development of microwave-driven NH 3 systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Zemach Moments and Radii of 2,3 H and 3,4 He

We present benchmark calculations of Zemach moments and radii of 2,3 H and 3,4 He using various few-body methods. Zemach moments are required to interpret muonic atom data measured by the CREMA collaboration at the Paul Scherrer Institute. Conversely, radii extracted from spectroscopic measurements can be compared with ab initio computations, posing stringent constraints on the nuclear model. For a given few-body method, different numerical procedures can be applied to compute these quantities. A detailed analysis of the numerical uncertainties entering the total theoretical error is presented. Uncertainties from the few-body method and the calculational procedure are found to be smaller than the dependencies on the dynamical modeling and the single nucleon inputs, which are found to be ≲2%. Furthermore when relativistic corrections and two-body currents are accounted for, the calculated moments and radii are in very good agreement with the available experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2010

Fundamental High-Pressure Calibration from All-Electron Quantum Monte Carlo Calculations

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high-pressure calibration based on the equation of state of cubic boron nitride. We compute the static contribution to the free energy with the QMC method and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. We compute the anharmonic Raman frequency shift with QMC simulations as a function of pressure and temperature, allowing optical pressure calibration. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron method is applicable to first-row solids, providing a new reference for ab initio calculations of solids and benchmarks for pseudopotential accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Porting AMG2013 to Heterogeneous CPU+GPU Nodes

LLNL's future advanced technology system SIERRA will feature heterogeneous compute nodes that consist of IBM PowerV9 CPUs and NVIDIA Volta GPUs. Conceptually, the motivation for such an architecture is quite straightforward: While GPUs are optimized for throughput on massively parallel workloads, CPUs strive to minimize latency for rather sequential operations. Yet, making optimal use of heterogeneous architectures raises new challenges for the development of scalable parallel software, e.g., with respect to work distribution. Porting LLNL's parallel numerical libraries to upcoming heterogeneous CPU+GPU architectures is therefore a critical factor for ensuring LLNL's future success in ful lling its national mission. One of these libraries, called HYPRE, provides parallel solvers and precondi- tioners for large, sparse linear systems of equations. In the context of this intern- ship project, I consider AMG2013 which is a proxy application for major parts of HYPRE that implements a benchmark for setting up and solving di erent systems of linear equations. In the following, I describe in detail how I ported multiple parts of AMG2013 to the GPU (Section 2) and present results for di erent experiments that demonstrate a successful parallel implementation on the heterogeneous ma- chines surface and ray (Section 3). In Section 4, I give guidelines on how my code should be used. Finally, I conclude and give an outlook for future work (Section 5).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗
DOE OSTI2020

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

Here, a conservative discontinuous Galerkin scheme for a non-linear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Using kernel-based statistical distance to study the dynamics of charged particle beams in particle-based simulation codes

Measures of discrepancy between probability distributions (statistical distance) are widely used in the fields of artificial intelligence and machine learning. We describe how certain measures of statistical distance can be implemented as numerical diagnostics for simulations involving charged-particle beams. Related measures of statistical dependence are also described. The resulting diagnostics provide sensitive measures of dynamical processes important for beams in nonlinear or high-intensity systems, which are otherwise difficult to characterize. Here, the focus is on kernel-based methods such as maximum mean discrepancy, which have a well-developed mathematical foundation and reasonable computational complexity. Several benchmark problems and examples involving intense beams are discussed. While the focus is on charged-particle beams, these methods may also be applied to other many-body systems such as plasmas or gravitational systems.

47 OTHER INSTRUMENTATION↗