SEARCH · Engineering Papers

Results for “parallel methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Report to NCSP on 2008 DANCE measurements of 233 U($\eta$,$\gamma$)

Uranium-233 plays an important role in the Th-U fuel cycle, with substantial production in the cycle. This cycle has been proposed as an alternative to the U-Pu fuel cycle due to its reduced production of transuranic elements. An accurate measurement of the 233 U($\eta$,$\gamma$) cross section is required by the National Criticality Safety Program (NCSP) to complete the neutron-induced cross section data, where experimental capture cross section data are scarce and were measured decades ago. The most recent capture cross section data available in the literature were measured in 2007 at the n_TOF facility (CERN); in the 60s measurements were performed at Rensselaer Polytechnic Institute (RPI) and at LANL. Finally, as reported by ORNL, a new evaluation with a revised (renormalized) fission cross section is needed on 233 U. The challenge for this measurement lies in the difficulty of measuring the capture cross section data in the competing fission background, as the fission cross section is around one order of magnitude larger than the capture cross section for 233 U. The accuracy of a capture cross section measurement depends on discrimination between $\gamma$’s produced in capture and fission reactions, for which an experimental setup combining capture and fission detectors is needed. For the ($\eta$,$\gamma$) cross section measurement at LANSCE, this discrimination is achieved by combining the Detector for Advanced Neutron Capture Experiments (DANCE), to measure $\gamma$’s from capture reactions, with a Parallel Plate Avalanche Counter (PPAC) to tag the $\gamma$’s produced by fission. This method was successfully used to measure 235 U and 239 Pu capture cross sections. In these measurements, the neutron capture cross section was determined in a large fission background well above 100 keV. As part of the NCSP nuclear data effort, we have looked at past DANCE measurements on 233 U($\eta$,$\gamma$) and evaluated whether existing data is adequate to apply this technique.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2014

Store-operate-coherence-on-value

A system, method and computer program product for performing various store-operate instructions in a parallel computing environment that includes a plurality of processors and at least one cache memory device. A queue in the system receives, from a processor, a store-operate instruction that specifies under which condition a cache coherence operation is to be invoked. A hardware unit in the system runs the received store-operate instruction. The hardware unit evaluates whether a result of the running the received store-operate instruction satisfies the condition. The hardware unit invokes a cache coherence operation on a cache memory address associated with the received store-operate instruction if the result satisfies the condition. Otherwise, the hardware unit does not invoke the cache coherence operation on the cache memory device.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Systems and methods for photovoltaic string protection

A system and method includes a circuit for protecting a photovoltaic string. A bypass switch connects in parallel to the photovoltaic string and a hot spot protection switch connects in series with the photovoltaic string. A first control signal controls opening and closing of the bypass switch and a second control signal controls opening and closing of the hot spot protection switch. Upon detection of a hot spot condition the first control signal closes the bypass switch and after the bypass switch is closed the second control signal opens the hot spot protection switch.

14 SOLAR ENERGY↗

DOE OSTI2013

Lattice study of quark distribution amplitudes in the pion and its excitations

Lattice QCD serves as a computational framework capable of predicting the spectrum of hadronic excitations from first principles. Our desire to describe the wealth of existing experimental data on the spectrum and to predict the outcomes of future experiments poses numerous challenges. Thus, obtaining an accurate resolution of excited states using methods of LQCD is complicated due to the faster decay of excited state correlation functions in Euclidean space in comparison with those of ground states, which we overcome through the use of anisotropic lattices with a finer temporal than spatial discretization. The aim of this project is to go beyond the spectrum to discern the structure of the states through the computation of the quark distribution amplitudes for both the ground and excited pion states on improved anisotropic lattices developed by the Hadron Spectrum Collaboration. Application of variational method allows us to extract the excited-state spectrum. When combined with undergoing parallel perturbative study of renormalization coefficients for quark bilinear operators, this work will enable us to explore the internal structure of the excited states, and to investigate the approach to a quark and gluon description of hadrons when probed at high-momentum transfers.

Mastropas, Ekaterina↗

DOE OSTI2014

Evaluation of Low Power ARM Processors for Monte Carlo Particle Transport: MCNP6 on the ARM Cortex-A8 - Paper 83

In the push towards exascale computing, the limiting factors are power requirements and heat dissipation. Using current technology, the power required to drive an exascale supercomputer would near 400-MW. Clearly, future computing hardware will have to be more energy efficient in order to achieve exascale performance. Most of the current research in advanced computing architectures is in the use of standard CPUs with accelerators. However, these accelerators are not generally supported by the traditional scientific programming languages (C, C++, and Fortran) and the scientific software must be tailored specifically to take advantage of the accelerators. Large scale scientific software (hundreds-of-thousand to millions of lines of code) can not be easily ported to the computing platform du jour. The widespread use of low power ARM based processors in the mobile computing market offer a low power alternative to accelerators that are also supported by the traditional scientific programming languages and programming paradigms such as MPI and OpenMP. Computing architectures for scientific computing are often assessed using standard benchmarks such as LINPACK benchmark or using suites of mini-applications that are greatly simplified stand-ins for scientific software. Often the real applications that will deployed on the scientific supercomputers are not ported or tested until delivery of the machines. In order to evaluate the use of ARM processors for scientific computing the Monte Carlo particle transport code, MCNP version 6 was compiled on the ARM Cortex-A8 and a set of benchmark problems were used to assess performance in comparison to an Intel Xeon E5-2670 and the low power Intel Atom N270. MCNP is a continuous energy Monte Carlo particle transport code in development at Los Alamos National Laboratory since 1977. MCNP is used to transport neutrons, photon, electrons, and charged particles using continuous energy interaction data. It is commonly used in fields of nuclear power, nuclear medicine, radiation protection, and high-energy physics through-out the world. MCNP version 6 is approximately 450 thousand lines of code, written primary in Fortran90, with some C and C++. MCNP supports distributed multiprocessing using MPI and threading with OpenMP. The ARM Cortex-A8 executes MCNP6 using about twice as much power as the Intel Xeon E5-2670 and about the same amount of power as the Intel Atom N270. However, these results are encouraging for the use of low-power ARM chips in future supercomputers. The ARM Cortex-A8 is a 32-bit CPU and it's NEON general-purpose SIMD (single-instruction multiple data) engine does not support double-precision floating-point arithmetic. More recent ARM processors (ARM Cortex-A53 and Cortex-A57) support 64-bit addressing and have NEON co-processors that support double-precision floating-point arithmetic. ARM claims that the A50 series delivers a performance increase of 15 times the Cortex-A8 performance for similar power usage. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2024

Reductive Analysis with Compiler-Guided Large Language Models for Input-Centric Code Optimizations

Input-centric program optimization aims to optimize code by considering the relations between program inputs and program behaviors. Despite its promise, a long-standing barrier for its adoption is the difficulty of automatically identifying critical features of complex inputs. This paper introduces a novel technique, reductive analysis through compiler-guided Large Language Models (LLMs), to solve the problem through a synergy between compilers and LLMs. It uses a reductive approach to overcome the scalability and other limitations of LLMs in program code analysis. The solution, for the first time, automates the identification of critical input features without heavy instrumentation or profiling, cutting the time needed for input identification by 44× (or 450× for local LLMs), reduced from 9.6 hours to 13 minutes (with remote LLMs) or 77 seconds (with local LLMs) on average, making input characterization possible to be integrated into the workflow of program compilations. Optimizations on those identified input features show similar or even better results than those identified by previous profiling-based methods, leading to optimizations that yield 92.6% accuracy in selecting the appropriate adaptive OpenMP parallelization decisions, and 20-30% performance improvement of serverless computing while reducing resource usage by 50-60%.

Input-Centric Optimization↗

DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗

DOE OSTI2018

Linking big models to big data: efficient ecosystem model calibration through Bayesian model emulation

Data-model integration plays a critical role in assessing and improving our capacity to predict ecosystem dynamics. Similarly, the ability to attach quantitative statements of uncertainty around model forecasts is crucial for model assessment and interpretation and for setting field research priorities. Bayesian methods provide a rigorous data assimilation framework for these applications, especially for problems with multiple data constraints. However, the Markov chain Monte Carlo (MCMC) techniques underlying most Bayesian calibration can be prohibitive for computationally demanding models and large datasets. We employ an alternative method, Bayesian model emulation of sufficient statistics, that can approximate the full joint posterior density, is more amenable to parallelization, and provides an estimate of parameter sensitivity. Analysis involved informative priors constructed from a meta-analysis of the primary literature and specification of both model and data uncertainties, and it introduced novel approaches to autocorrelation corrections on multiple data streams and emulating the sufficient statistics surface. We report the integration of this method within an ecological workflow management software, Predictive Ecosystem Analyzer (PEcAn), and its application and validation with two process-based terrestrial ecosystem models: SIPNET and ED2. In a test against a synthetic dataset, the emulator was able to retrieve the true parameter values. A comparison of the emulator approach to standard brute-force MCMC involving multiple data constraints showed that the emulator method was able to constrain the faster and simpler SIPNET model's parameters with comparable performance to the brute-force approach but reduced computation time by more than 2 orders of magnitude. The emulator was then applied to calibration of the ED2 model, whose complexity precludes standard (brute-force) Bayesian data assimilation techniques. Both models are constrained after assimilation of the observational data with the emulator method, reducing the uncertainty around their predictions. Performance metrics showed increased agreement between model predictions and data. Our study furthers efforts toward reducing model uncertainties, showing that the emulator method makes it possible to efficiently calibrate complex models.

58 GEOSCIENCES↗

DOE OSTI2021

Efficient exascale discretizations: High-order finite element methods

Efficient exploitation of exascale architectures requires rethinking of the numerical algorithms used in many large-scale applications. These architectures favor algorithms that expose ultra fine-grain parallelism and maximize the ratio of floating point operations to energy intensive data movement. One of the few viable approaches to achieve high efficiency in the area of PDE discretizations on unstructured grids is to use matrix-free/partially assembled high-order finite element methods, since these methods can increase the accuracy and/or lower the computational time due to reduced data motion. In this paper we provide an overview of the research and development activities in the Center for Efficient Exascale Discretizations (CEED), a co-design center in the Exascale Computing Project that is focused on the development of next-generation discretization software and algorithms to enable a wide range of finite element applications to run efficiently on future hardware. CEED is a research partnership involving more than 30 computational scientists from two US national labs and five universities, including members of the Nek5000, MFEM, MAGMA and PETSc projects. We discuss the CEED co-design activities based on targeted benchmarks, miniapps and discretization libraries and our work on performance optimizations for large-scale GPU architectures. We also provide a broad overview of research and development activities in areas such as unstructured adaptive mesh refinement algorithms, matrix-free linear solvers, high-order data visualization, and list examples of collaborations with several ECP and external applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Hierarchical calibration and validation for modeling bench-scale solvent-based carbon capture. Part 1: Non-reactive physical mass transfer across the wetted wall column

A hierarchical model calibration and validation are proposed here for quantifying the confidence level of mass transfer prediction using a computational fluid dynamics (CFD) model, where the solvent-based carbon dioxide (CO 2 ) capture is simulated and the predicted results are compared to the corresponding bench-scale experimental data. Two unit problems with an increasing level of complexity are proposed to break down the complex physical/chemical processes of solvent-based CO 2 capture into relatively simpler problems to separate the effects of physical transport and chemical reaction. This paper focuses on the calibration and validation of the first unit problem, i.e., the CO 2 mass transfer across a falling ethanolamine (MEA) film in the absence of chemical reaction. This problem is investigated both experimentally and numerically using non-reactive nitrous oxide (N 2 O) as a surrogate for CO 2 . To capture motion of the gas-liquid interface, a volume of fluid method is employed with a one-fluid formulation to compute the mass transfer between the two phases. Parallel bench-scale experiments are designed and conducted to validate and calibrate the CFD models using a general Bayesian calibration approach. Two important transport parameters, Henry's constant and gas diffusivity, are calibrated to produce the posterior distributions, which will be used as input for the second unit problem to address the chemical absorption of CO 2 across the MEA falling film, where both mass transfer and chemical reaction are involved. Finally, mass transfer coefficients predicted by CFD are also compared with those predicted by the traditional/empirical correlations under both steady and wavy falling film conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Parallel accelerated vector similarity calculations for genomics applications

The surge in availability of genomic data holds promise for enabling determination of genetic causes of observed individual traits, with applications to problems such as discovery of the genetic roots of phenotypes, be they molecular phenotypes such as gene expression or metabolite concentrations, or complex phenotypes such as diseases. However, the growing sizes of these datasets and the quadratic, cubic or higher scaling characteristics of the relevant algorithms pose a serious computational challenge necessitating use of leadership scale computing. In this study we describe a new approach to performing vector similarity metrics calculations, suitable for parallel systems equipped with graphics processing units (GPUs) or Intel Xeon Phi processors. Our primary focus is the Proportional Similarity metric applied to Genome Wide Association Studies (GWAS) and Phenome Wide Association Studies (PheWAS). We describe the implementation of the algorithms on accelerated processors, methods used for eliminating redundant calculations due to symmetries, and techniques for efficient mapping of the calculations to many-node parallel systems. Finally, results are presented demonstrating high per-node performance and parallel scalability with rates of more than five quadrillion (5 × 10 15 ) elementwise comparisons achieved per second on the ORNL Titan system. In a companion paper we describe corresponding techniques applied to calculations of the Custom Correlation Coefficient for comparative genomics applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

An O(N logN) multi-layer boundary element method for direct computation of sound propagation in shallow water environments

Direct three-dimensional (3D) numerical simulations of acoustic fields in range-dependent shallow water environments remains a challenge due to environmental complexities and large computational cost.We develop an efficient 3D boundary element method (BEM) for shallow water acoustic propagation which utilizes a Pre-corrected Fast Fourier Transform (PFFT) approach to reduce the computational effort from O ( N 2~3 ) to O ( N log N ) where N is the total number of boundary unknowns. To account for inhomogeneous media, the method allows for arbitrary number of coupled multi-layer BEM sub-domains. With O ( N log N ) efficiency and the use of massively parallel high-performance computing platforms, we are able to conduct multi-layer 3D direct simulations of low-mid frequency acoustics over kilometer ranges. Furthermore we perform extensive validations of the method and provide two shallow water waveguide examples benchmarked against theoretical solutions. To illustrate the efficacy and usefulness of the PFFT-BEM method, we perform 3D large-scale direct simulations to assess the performance of two established canonical models: axisymmetric coupled mode model for 3D seamount; and Kirchhoff approximation and perturbation theory for 3D rough surface scattering.

42 ENGINEERING↗

DOE OSTI2017

Anisotropic Thermal Response of Packed Copper Wire

The apparent thermal conductivity of packed copper wire test specimens was measured parallel and perpendicular to the axis of the wire using laser flash, transient plane source, and transmittance test methods. Approximately 50% wire packing efficiency was produced in the specimens using either 670- or 925-μm-diameter copper wires that both had an insulation coating thickness of 37 μm. The interstices were filled with a conventional varnish material and also contained some remnant porosity. The apparent thermal conductivity perpendicular to the wire axis was about 0.5–1 W/mK, whereas it was over 200 W/mK in the parallel direction. The Kanzaki model and an finite element analysis (FEA) model were found to reasonably predict the apparent thermal conductivity perpendicular to the wires but thermal conductivity percolation from nonideal wire-packing may result in their underestimation of it.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2018

Growth and characterization of HgBa 2 CaCu 2 O 6+δ and HgBa 2 Ca 2 Cu 3 O 8+δ crystals

We report the successful growth of sizable crystals of the cuprate superconductors HgBa 2 CaCu 2 O 6+δ and HgBa 2 Ca 2 Cu 3 O 8+δ . These compounds are well known for their high optimal superconducting critical temperatures of T c = 128 K and 134 K at ambient pressure, respectively, and for their challenging synthesis. Using a conventional quartz-tube encapsulation method and a two-layer encapsulation method that utilizes custom-built high-pressure furnaces, we are able to grow single crystals with linear dimensions up to several millimeters parallel to the CuO 2 planes. Extended postgrowth annealing is shown to lead to sharp superconducting transitions, indicative of high macroscopic homogeneity. X-ray diffraction and polarized Raman spectroscopy are identified as viable nondestructive methods to help separate the two compounds from synthesis products, as the latter often contain intergrowths of the former. Furthermore, our work helps to remove obstacles toward the study of these model cuprate systems with experimental probes that require sizable high-quality crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Spin-Wave Frequency Division Multiplexing in an Yttrium Iron Garnet Microstripe Magnetized by Inhomogeneous Field.

Spin waves are promising candidates for information processing and transmission in a broad frequency range. In the realization of magnonic devices, the frequency related division of the spin waves is a critical function for parallel information processing. In this work, we demonstrate a proof-of-concept spin-wave frequency division multiplexing method by magnetizing a homogeneous magnetic microstripe with an inhomogeneous field. The symmetry breaking additional field is introduced by a Permalloy stripe simply placed in lateral proximity to an yttrium iron garnet waveguide. Spin waves with different frequencies can propagate independently, simultaneously, and separately in space along the shared waveguide. This work demonstrates one potential way for parallel information transmission and processing in magnonics.

Zhang, Zhizhi↗

DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Efficient smoothed particle radiation hydrodynamics I: Thermal radiative transfer

This work presents efficient solution techniques for radiative transfer in the smoothed particle hydrodynamics discretization. Two choices that impact efficiency are how the material and radiation energy are coupled, which determines the number of iterations needed to converge the emission source, and how the radiation diffusion equation is solved, which must be done in each iteration. The coupled material and radiation energy equations are solved using an inexact Newton iteration scheme based on nonlinear elimination, which reduces the number of Newton iterations needed to converge within each time step. During each Newton iteration, the radiation diffusion equation is solved using Krylov iterative methods with a multigrid preconditioner, which abstracts and optimizes much of the communication when running in parallel. The code is verified for an infinite medium problem, a one-dimensional Marshak wave, and a two and three-dimensional manufactured problem, and exhibits first-order convergence in time and second-order convergence in space. For these problems, the number of iterations needed to converge the inexact Newton scheme and the diffusion equation is independent of the number of spatial points and the number of processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Overlapping localized exponential time differencing methods for diffusion problems

The localized exponential time differencing (ETD) based on overlapping domain decomposition has been recently introduced for extreme-scale phase field simulations of coarsening dynamics, which displays excellent parallel scalability in supercomputers. This paper serves as the first step toward building a solid mathematical foundation for this approach. We study the overlapping localized ETD schemes for a model time-dependent diffusion equation discretized in space by the standard central difference. Two methods are proposed and analyzed for solving the fully discrete localized ETD systems: the first one is based on Schwarz iteration applied at each time step and involves solving stationary problems in the subdomains at each iteration, while the second one is based on the Schwarz waveform relaxation algorithm in which time-dependent subdomain problems are solved at each iteration. The convergences of the associated iterative solutions to the corresponding fully discrete localized ETD solution and to the exact semidiscrete solution are rigorously proved. Finally, numerical experiments are also carried out to confirm theoretical results and to compare the performance of the two methods.

58 GEOSCIENCES↗