Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “IMPLEMENTATION”

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 145 records · Page 8

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗

Potential reductions in fine particulate matter and premature mortality following implementation of air pollution controls on coal-fired power plants in India

Coal-fired power plants (CFPPs) account for > 70% of electricity generation in India, but < 5% of facilities have installed technologies for sulfur dioxide (SO 2 ) and nitrogen oxide (NO X ) removal. Emissions of these pollutants lead to the formation of fine particulate matter (PM 2.5 ) and an increased risk of premature mortality for exposed populations. Here, we use a nested version of the GEOS-Chem global chemical transport model (0.5° × 0.625° resolution) for India to estimate reductions in PM 2.5 concentrations that could have been achieved by implementing existing emission control technologies like flue-gas desulfurization (FGD) and/or selective catalytic reduction (SCR). We quantify the associated burden of disease using the integrated exposure response (IER) and global exposure mortality model (GEMM) functions and compare the costs of premature mortality to those for FGD installation. Model simulations for 2010 suggest installation of FGD would have reduced mean annual PM 2.5 concentrations across India by 8%, compared to 3% with SCR installation, and 11% with both FGD and SCR. A 7–28% reduction in PM 2.5 was simulated for local communities closest to CFPPs (same model grid cell), leading to up to 17% reduction in annual premature mortality. Overall, more than 0.21–0.48 million premature deaths would have been avoided over a 10-year period if FGD had been implemented on all CFPPs, compared to 0.09–0.21 million with SCR and 0.22–0.72 million with both FGD and SCR. Benefits associated with such actions are approximately $\$18.1$–$\$604$ billion USD per year, which is equivalent to ~ 0.44 to 10% of India’s GDP. These results suggest that monetary benefits from avoided premature mortality far outweigh the capital and operational costs of FGD and/or SCR installation of $\$19.5$ billion and/or $\$32.8$ billion per year, respectively. This information is essential because the high costs of installation and operation are often given as reasons for delaying installation and commissioning. Finally, we conclude that policy actions to control air pollution from CFPPs are economically justifiable.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Implementation of the riding hydrogen model in CCTBX to support the next generation of X-ray and neutron joint refinement in Phenix

A fundamental prerequisite for implementing new procedures of atomic model refinement against neutron diffraction data is the efficient handling of hydrogen atoms. The riding hydrogen model, which constrains hydrogen atom parameters to those of the non-hydrogen atoms, is a plausible parameterization for refinements. This work describes the implementation of the riding hydrogen model in the Computational Crystallography Toolbox and in Phenix. Riding hydrogen atoms can be found in several different configurations that are characterized by specific geometries. For each configuration, the hydrogen atom parameterization and the expressions for the gradients of refinement target function with respect to non-hydrogen parameters are described.

H atom parameterization↗

Lattice Physics Calculations Using the Embedded Self-Shielding Method in Polaris, Part I: Methods and Implementation

Polaris is a 2-dimensional multigroup lattice physics capability in the SCALE code system for the analysis of light water reactor fuel designs. The goal of light water reactor lattice physics codes is to generate few-group homogenized cross sections for downstream full-core nodal diffusion calculations. Additionally, lattice physics calculations contain three primary components: the cross section processing calculation, the 2D transport calculation, and the depletion calculation. This paper summarizes the calculational methods and their implementation into Polaris, with an emphasis on implementation of the embedded self-shielding method. The accuracy of the embedded self-shielding method depends on the procedure used to generate self-shielding factors on the multigroup library. Numerical benchmarks calculations reveal that the accuracy of Polaris eigenvalue predictions is enhanced by (1) using heterogeneous unit cell models to generate the self-shielding factors on the library and (2) using self-shielding factors for within-group scattering cross section.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Implementing multi-settlement decentralized electricity market design for transactive communities with imperfect communication

Recent advances in information and communication technologies and smart metering, provides strategic opportunities for ``prosumers" to reform their conventional energy practices towards more consumer-centric economies. From an operational perspective, managing power distribution networks is becoming more difficult with such active grid-edge systems providing limited to no visibility or control. Transactive Energy (TE) has been emerging as a key enabler towards effectively and efficiently integrating prosumers into competitive electricity markets. This work presents a transactive implementation of community-centric markets. A co-simulation framework is developed for evaluating the proposed market structure with high-fidelity models. Case studies on the IEEE-123 node test system demonstrate that community-centric transactive markets can enable communities of prosumers to operate collaboratively as grid-edge systems. The potential benefits of implementing community-centric TE systems are also illustrated.

Mukherjee, Monish↗

FlowPM: Distributed TensorFlow implementation of the FastPM cosmological N-body solver

Here, we present FlowPM, a Particle-Mesh (PM) cosmological N-body code implemented in Mesh-TensorFlow for GPU-accelerated, distributed, and differentiable simulations. We implement and validate the accuracy of a novel multi-grid scheme based on multiresolution pyramids to compute large-scale forces efficiently on distributed platforms. We explore the scaling of the simulation on large-scale supercomputers and compare it with corresponding Python based PM code, finding on an average 10x speed-up in terms of wallclock time. We also demonstrate how this novel tool can be used for efficiently solving large scale cosmological inference problems, in particular reconstruction of cosmological fields in a forward model Bayesian framework with hybrid PM and neural network forward model. We provide skeleton code for these examples and the entire code is publicly available at https://github.com/modichirag/flowpm[Formula presented].

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation of an extensible property modeling framework in ESPEI with applications to molar volume and elastic stiffness models

Property models are becoming more widely adopted by commercial Calphad databases, but they are not nearly as common in non-commercial or traditional academic Calphad databases. A primary driver is that user-friendly Calphad modeling tools that support property models are not widely available. Here we present new property modeling capabilities that have been implemented in ESPEI (the Extensible, Self-optimizing Phase Equilibrium Infrastructure). These capabilities include both generating property model parameters from data and improvements to the algorithmic selection of the most appropriate model from a series of candidates. Additionally, two illustrative examples are given that use ESPEI to fit different property models. First, we generate molar volume model parameters for Group IV, V, and VI refractory BCC alloys based on the model by Lu et al. (2005). Second, we demonstrate the extensibility of ESPEI’s property modeling capabilities by implementing a custom PyCalphad model for BCC elastic stiffness parameters to generate and compare parameters to the ones assessed by Marker et al. (2018) using the same data. Property models generated by ESPEI can be used in PyCalphad or further optimized with uncertainty quantification using ESPEI.

36 MATERIALS SCIENCE↗

Engineering intricacies of implementing single-atom alloy catalysts for low-temperature electrocatalytic CO 2 reduction

Catalysts research for electrocatalytic CO 2 reduction reactions (CO 2 R) has undergone rapid growth in the last decade. Single-atom alloy catalysts (SAAs) featuring atomically dispersed metal dopants on host metal surfaces have shown promises in boosting CO 2 R yield by optimizing the structure and synergy of the catalytic metals at the atomic scale. Despite the exciting development of SAAs for CO 2 R in fundamental science, dedicated studies for its engineering implementation have been absent. We use this perspective to discuss our non-exhaustive engineering considerations for implementing SAAs for CO 2 R. Here, the perspective starts with a brief overview of the current research status for SAAs in CO 2 R, followed by focal points on structure uncertainties associated with catalyst manufacturing, catalyst layer degradation during reaction, and possibilities for SAAs to mitigate the salt precipitation issue at the device level. We hope our opinions will engage increasing attention toward the engineering catalysis research for applying SAAs to CO 2 R at scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward implementation of carbon-conservation networks in nonmodel organisms

Decarboxylation — the release of carbon dioxide (CO 2 ) from a substrate — reduces the carbon yield of bioproduced chemicals. When overlaid onto central carbon metabolism, carbon-conservation networks (CCNs) that reroute flux around CO 2 release can theoretically achieve higher carbon yields for products derived from intermediates that traditionally require CO 2 release, such as acetyl-CoA. Recently, CCNs have started to be implemented in model organisms to produce compounds at higher carbon yields. However, implementation of CCNs in nonmodel hosts may have the greatest impact given their ability to assimilate a larger array of feedstocks, greater environmental tolerance, and unique biosynthetic pathways, ultimately enabling access to a wider range of products. Here, we review recent advances in CCNs with a focus on their application to nonmodel organisms. The differences in central carbon metabolism among different nonmodel hosts reveal opportunities to engineer and apply new CCNs.

59 BASIC BIOLOGICAL SCIENCES↗

Scalable implementation of polynomial filtering for density functional theory calculation in PARSEC

In this work, we present an efficient implementation of polynomial filtering methods in PARSEC, a real-space pseudopotential based Kohn–Sham density functional theory solver. The implementation described here improves upon a Chebyshev-filtered subspace iteration algorithm used in the previous version of PARSEC. We present a hybrid polynomial filtering scheme that combines Chebyshev-filtered subspace iteration and a spectrum slicing method that partitions the spectrum into several spectral slices and uses bandpass-filtered subspace iteration to compute approximate eigenpairs within each interior slice simultaneously. We describe a procedure to partition a spectrum and construct polynomial filters. We also discuss a number of practical issues such as the use of appropriate data layouts for carrying out the computation on a two-dimensional process grid and how to achieve good load balance by allocating an appropriate number of process groups to each spectral slice. Numerical examples are presented to demonstrate the effectiveness of the hybrid polynomial filtering method as well as the superior parallel scalability of spectrum slicing in comparison to that of Chebyshev-filtered subspace iteration.

97 MATHEMATICS AND COMPUTING↗

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. As a result, the utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed.

36 MATERIALS SCIENCE↗

Dataset describing two reference models for full-spectral lighting and daylight simulations together with implementations for two software systems

A dataset of two spectral lighting simulation reference models - one office and one factory hall - is presented. It aims to demonstrate and support full-spectral daylight and electric lighting simulations and facilitate evaluation of non-visual effects of light. The dataset includes Rhino CAD geometry, comprehensive spectral material and light source data and window system BSDF data. Example implementations in the two software tools, Radiance and OWL, enable reproducible workflows and support adoption in other software. The dataset is openly available on Zenodo. The office model reproduces Room 518 at the University of Innsbruck, including a west-facing façade and interior furnishings. The factory hall model follows the proposed geometry in the European standard 15193 for building energy performance. Interior reflectances in the office were measured in-situ using a handheld spectrometer. Exterior spectra and factory hall materials matching specified reflectances were obtained from an online spectral materials database. Glazing transmittance was derived from IGDB data using LBNL Optics/WINDOW. BSDFs for venetian blinds at various tilt angles, and for a diffusing pane adapted from the Complex Glazing Database, were generated in WINDOW. Luminaires in both models are specified with photometric files (Eulumdat/IES) and lamp spectra (Fluorescent 840, 4000 K LED). The provided example implementations (Radiance, OWL) include prepared input data and scripts to run first spectral simulations; example results are also included. The dataset is prepared to support reuse by researchers, designers and software developers for method validation, software engineering and comparison, and development of spectral metrics and controls.

Geisler-Moroder, David↗

Implementation and test of an automated control hunting fault correction algorithm in a fault detection and diagnostics tool

Control hunting due to improper proportional–integral–derivative (PID) parameters in the building automation system (BAS) is one of the most common faults identified in commercial buildings. It can cause suboptimal performance and early failure of heating, ventilation, and air conditioning (HVAC) equipment. Commercial fault detection and diagnostics (FDD) software represents one of the fastest growing market segments in smart building technologies in the United States. Implementation of PID retuning procedures as an auto-correction algorithm and integration into FDD software has the potential to mitigate control hunting across a heterogeneous portfolio of buildings with different BAS in a scalable way. This paper presents the development, implementation, and field testing of an automated control hunting fault correction algorithm based on lambda tuning open-loop rules. The algorithm was developed in a commercial FDD software and successfully tested among nine variable air volume boxes in an office building in the United States. The paper shows the feasibility of using FDD tools to automatically correct control hunting faults, discusses scalability considerations, and proposes a path forward for the HVAC industry and academia to further improve this technology.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.

54 ENVIRONMENTAL SCIENCES↗

QoS-aware edge AI placement and scheduling with multiple implementations in FaaS-based edge computing

Resource constraints on the computing continuum require that we make smart decisions for serving AI-based services at the network edge. AI-based services typically have multiple implementations (e.g., image classification implementations include SqueezeNet, DenseNet, and others) with varying trade-offs (e.g., latency and accuracy). The question then is how should AI-based services be placed across Function-as-a-Service (FaaS) based edge computing systems in order to maximize total Quality-of-Service (QoS). To address this question, we propose a problem that jointly aims to solve (i) edge AI service placement and (ii) request scheduling. These are done across two time-scales (one for placement and one for scheduling). Here we first cast the problem as an integer linear program. We then decompose the problem into separate placement and scheduling subproblems and prove that both are NP-hard. We then propose a novel placement algorithm that places services while considering device-to-device communication across edge clouds to offload requests to one another. Our results show that the proposed placement algorithm is able to outperform a state-of-the-art placement algorithm for AI-based services, and other baseline heuristics, with regard to maximizing total QoS. Additionally, we present a federated learning-based framework, FLIES, to predict the future incoming service requests and their QoS requirements. Our results also show that our FLIES algorithm is able to outperform a standard decentralized learning baseline for predicting incoming requests and show comparable predictive performance when compared to centralized training.

97 MATHEMATICS AND COMPUTING↗

Code-verification techniques for the method-of-moments implementation of the magnetic-field integral equation

For computational physics simulations, code verification plays a major role in establishing the credibility of the results by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, surface integral equations, such as the method-of-moments implementation of the magnetic-field integral equation, are frequently used to solve Maxwell's equations on the surfaces of electromagnetic scatterers. These electromagnetic surface integral equations yield many code-verification challenges due to the various sources of numerical error and their possible interactions. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Code-verification techniques for the method-of-moments implementation of the combined-field integral equation

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. Here, we demonstrate the effectiveness of these approaches for cases with and without coding errors.

97 MATHEMATICS AND COMPUTING↗

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗