Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel coordinates”

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 127 records · Page 7

Quantization and symmetry in periodic coverage patterns with applications to earth observation

Analytical approaches based on an idealized physical model and concepts from number theory show that in periodic coverage patterns, uniquely defined by their revolution numbers R (orbital) and N (rotational), the subnodal points are earth-fixed, and they divide the equator into R equal segments of length s. The ascending subsatellite trace crosses each point once (only) each period. The descending subnodal points coincide with the ascending points if the integers N and R have like parity, and bisect the intervals between them if opposite. The interval between consecutive unidirectional crossings is Ns. Symmetries extend the equatorial results to all parallels of latitude. Complete periodic patterns of traces exhibit an overall symmetry, with trace intersections confined to discrete coordinate values which are quantized in longitude (basic s-unit) and symmetric in latitude.

King, J. C.↗

Geometric registration and rectification of spaceborne SAR imagery

This paper describes the development of automated location and geometric rectification techniques for digitally processed synthetic aperture radar (SAR) imagery. A software package has been developed that is capable of determining the absolute location of an image pixel to within 60 m using only the spacecraft ephemeris data and the characteristics of the SAR data collection and processing system. Based on this location capability algorithms have been developed that geometrically rectify the imagery, register it to a common coordinate system and mosaic multiple frames to form extended digital SAR maps. These algorithms have been optimized using parallel processing techniques to minimize the operating time. Test results are given using Seasat SAR data.

Curlander, J. C.↗

Hybrid solid element with a traction-free cylindrical surface

An eight node solid element with two parallel faces and one traction-free cylindrical surface is derived using the assumed stress hybrid method. Cylindrical coordinates are used so that the assumed stresses satisfy the equilibrium equations as well as the traction-free condition over the cylindrical boundary. In the limiting case of plane stress conditions the assumed stresses also satisfy the compatibility conditions. Example solutions have demonstrated the advantage of using this special element for analyzing solids with circular holes.

Pian, T. H. H.↗

Recognizing Patterns In Log-Polar Coordinates

Log-Hough transform is basis of improved method for recognition of patterns - particularly, straight lines - in noisy images. Takes advantage of rotational and scale invariance of mapping from Cartesian to log-polar coordinates, and offers economy of representation and computation. Unification of iconic and Hough domains simplifies computations in recognition and eliminates erroneous quantization of slopes attributable to finite spacing of Cartesian coordinate grid of classical Hough transform. Equally efficient recognizing curves. Log-Hough transform more amenable to massively parallel computing architectures than traditional Cartesian Hough transform. "In-place" nature makes it possible to apply local pixel-neighborhood processing.

Weiman, Carl F. R.↗

COHORT: Coordination of Heterogeneous Thermostatically Controlled Loads for Demand Flexibility

Demand flexibility is increasingly important for power grids. Careful coordination of thermostatically controlled loads (TCLs) can modulate energy demand, decrease operating costs, and increase grid resiliency. We propose a novel distributed control framework for the Coordination Of HeterOgeneous Residential Thermostatically controlled loads (COHORT). COHORT is a practical, scalable, and versatile solution that coordinates a population of TCLs to jointly optimize a grid-level objective, while satisfying each TCL’s end-use requirements and operational constraints. To achieve that, we decompose the grid-scale problem into subproblems and coordi- nate their solutions to find the global optimum using the alternating direction method of multipliers (ADMM). The TCLs’ local problems are distributed to and computed in parallel at each TCL, making COHORT highly scalable and privacy-preserving. While each TCL poses combinatorial and non-convex constraints, we characterize these constraints as a convex set through relaxation, thereby making COHORT computationally viable over long planning horizons. After coordination, each TCL is responsible for its own control and tracks the agreed-upon power trajectory with its preferred strategy. In this work, we translate continuous power back to discrete on/off actuation, using pulse width modulation. COHORT is generalizable to a wide range of grid objectives, which we demonstrate through three distinct use cases: generation following, minimizing ramping, and peak load curtailment. In a notable experiment, we validated our approach through a hardware-in-the-loop simulation, including a real-world air conditioner (AC) controlled via a smart thermostat, and simulated instances of ACs modeled after real-world data traces. During the 15-day experimental period, COHORT reduced daily peak loads by an average of 12.5% and maintained comfortable temperatures.

demand response↗

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗

Neural controller for adaptive movements with unforeseen payloads

A theory and computer simulation of a neural controller that learns to move and position a link carrying an unforeseen payload accurately are presented. The neural controller learns adaptive dynamic control from its own experience. It does not use information about link mass, link length, or direction of gravity, and it uses only indirect uncalibrated information about payload and actuator limits. Its average positioning accuracy across a large range of payloads after learning is 3 percent of the positioning range. This neural controller can be used as a basis for coordinating any number of sensory inputs with limbs of any number of joints. The feedforward nature of control allows parallel implementation in real time across multiple joints.

Kuperstein, Michael↗

Dispersion relation for bianisotropic materials and its symmetry properties

The dispersion relation for an arbitrary general bianisotropic medium is derived in Cartesian coordinates, in a form well suited to imposing the boundary conditions when dealing with layered media with planar and parallel interfaces. Special cases of practical interest are also considered. Eleven fundamental coefficient families are identified by considering in detail all the symmetries present in the dispersion relation. An ad hoc expression of the determinant of the sum of two 3 x 3 matrices permits the use of a simple procedure to obtain the coefficients of the dispersion equation. The discussed symmetry properties have general validity, and this technique to evaluate the coefficients may be useful in other fields of application where dispersion relations are of importance.

Graglia, Roberto D.↗

High Temperature Anode Recycle Blower for Solid Oxide Fuel Cell, Phase II (Final Report)

Broad commercialization of solid oxide fuel cells (SOFCs) requires anode offgas recycle blowers (ARCB) that are specifically designed for handling the challenging operating conditions presented by the SOFC process gases. Otherwise, they can be susceptible to frequent maintenance, low reliability, and short life. This report presents details of a DOE-funded Phase II effort conducted by Mohawk Innovative Technology, Inc. (MiTi®) for the development and testing of an oil-free, low cost, high-temperature centrifugal ARCB based on compliant foil bearing (CFB) technology for support of a 100 kW solid oxide fuel cell power plant. This Phase II builds on the results of a successful Phase I development project that resulted in the demonstration of a low TRL-6 prototype, shown in Figure 1. This Phase II final report presents design improvements over the Phase I prototype, presents the assembled test ARCBs, and discusses the advantages of this novel technology, particularly its long life and maintenance-free operation, which are derived from the use of CFBs. Also included are preliminary techno-economic analysis considerations for cost-effective deployment of the technology. The specific objectives of this project, as stated in the statement of program objectives (SOPO) were 1) to follow the methodology of design for manufacturing (or manufacturability) and assembly for implementing improvements identified as part of the Phase I effort with the purpose of reducing cost, enabling mass production, and facilitating the commercialization of a revised ARCB design, and 2) to fabricate four complete ARCB units based on the revised design and demonstrating their performance in both a laboratory setting and in an actual SOFC power plant. Execution of this SOPO would be supported by a number of technical tasks resulting in the fabrication and testing of the prototypes. To this end, MiTi continued the teaming relationship started during Phase I of the program with subcontracting partner FuelCell Energy, Inc. (FCE), which integrated (under a parallel effort funded by DOE Award DE-FE0026199) a modular 200 kWe SOFC power plant. MiTi and FCE coordinated to ultimately incorporate one of MiTi’s ARCB prototypes into one of the 100 kWe SOFC Modular Power Blocks (MPB) that constitute the core of FCE’s SOFC power plant for in-situ long-duration testing. Additionally, MiTi would explore the scalability and extendibility of the technology to other applications, as well as conduct a basic techno-economic analysis.

03 NATURAL GAS↗

Parallel transport sweeps on two-dimensional cartesian and hexagonal grids

This paper aims to provide a proof of concept for parallel transport sweeps on two-dimensional hexagonal grids for the discrete ordinates transport equation. While the method is an extension of the popular and well-established Koch-Baker-Alcoulffe (KBA) algorithm, there are significant differences between the cartesian and hexagonal grid and thereafter sweep. The most important is the three-way connectivity of hexagons within the grid which creates greater dependencies between the elements. The KBA method in structured orthogonal grids was first implemented in the DRAGON5 code and the method is first described here. The differences in implementation for the hexagonal grid are also described. Benchmark results are also presented, showing roughly 10 times speedup in computational times with roughly 100 processors, in both cases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

General field evaluation in high-order meshes on GPUs

Robust and scalable function evaluation at any arbitrary point in the finite/spectral element mesh is required for querying the partial differential equation solution at points of interest, comparison of solution between different meshes, and Lagrangian particle tracking. This is a challenging problem, particularly for high-order unstructured meshes partitioned in parallel with MPI, as it requires identifying the element that overlaps a given point and computing the corresponding reference space coordinates. Here, we present a robust and efficient technique for general field evaluation in large-scale high-order meshes with quadrilaterals and hexahedra. In the proposed method, a combination of globally partitioned and processor-local maps are used to first determine a list of candidate MPI ranks, and then locally candidate elements that could contain a given point. Next, element-wise bounding boxes further reduce the list of candidate elements. Finally, Newton’s method with trust region is used to determine the overlapping element and corresponding reference space coordinates. Since GPU-based architectures have become popular for accelerating computational analyses using meshes with tensor-product elements, specialized kernels have been developed to utilize the proposed methodology on GPUs. The method is also extended to enable general field evaluation on surface meshes. The paper concludes by demonstrating the use of the proposed method in various applications ranging from mesh-to-mesh transfer during r-adaptivity to Lagrangian particle tracking.

97 MATHEMATICS AND COMPUTING↗

Software Control Program For Transportable Microgrid State-of-charge Balancing And Frequency Stability Controls

A deterministic state-of-charge (SOC) balancing approach software control code is introduced as an integral secondary management to primary control layer of an islanded small microgrid or nanogrid system made up of multiple grid-forming inverter/battery/solar combination systems, where each set of batteries with each inverter are on independent DC buses (i.e. non-paralleled on the DC sides). A DERMS-level control approach, algorithm and automation controller program was developed to improve coordination and enable microgrid asset compliance and SOC balancing, enabling provision of a system-level power stability support architecture, load support, and asset scalability. The architecture is configured to treat each unit or micro/nano-grid as a node in a microgrid network, allowing for autonomous DERMS control regarding load and SOC balancing and power stability. As the network grows with the addition of units, greater coordination efforts may be required. The ideal small network microgrid ranges from 2-10 inverter/battery units before additional control parameters must be considered in the existing architecture. The control approach focuses on a deterministic state-of-charge analysis as the primary level control process followed by a secondary control loop using a forced frequency-watt droop strategy to conform off-the-shelf components into behaving under a leader-follower configuration. Adopting this control scheme has been shown to allow for a balanced, unit-coordinated microgrid network, enabling stable power flow. The deterministic state-of-charge approach is introduced as an integral primary control layer of an islanded small network microgrid. A standard strategy for SOC balancing is implementing a battery management system (BMS) to control SOC on the DC side. An alternative approach is to determine how to coordinate sending and receiving power on the AC side with multiple units. The latter approach assesses all the integrated units in the microgrid network. Once the individual units are identified, further system data is required to calculate each unit's total kWh, provided information about its capability to supply or consume kWh and availability. The secondary control layer in the multi-layered small network microgrid methodology uses the primary layer’s decision to initiate frequency setpoint changes, initializing the SOC balancing. The secondary control layer considers numerous system-dependent variables to enable a charging and discharging profile based on adjustable frequency setpoints. The combined architecture will result in stable, coordinated power flow enhancing an AC microgrid's functionalities.

Myers, KurtS [Idaho National Laboratory (INL), Ida↗

Ignition of confined gaseous mixtures by hot surfaces and hot wires

Ignition times and spatial and temporal variations of temperature and concentration in gaseous mixtures confined between two infinite parallel walls or two infinite cylinders have been obtained by numerical integration of the appropriate conservation equations written in Lagrangian coordinates. Ignition times and ignition energies are presented for the case of an isothermal wall in terms of the initial mixture pressure and equivalence ratio for both one and two-step chemical reaction mechanisms. The numerical results indicate that there is a critical mixture pressure for which the ignition time is minimum. The values of this critical pressure are larger (smaller) than 1 atm for the one- (two-) step reaction mechanism. The critical pressure for the ignition time is not equal to the critical pressure for the ignition energy. The ignition time and energy decrease with the equivalence ratio within a certain range and then remain constant.

Ramos, J. I.↗

Supercomputing 2002: NAS Demo Abstracts

The hyperwall is a new concept in visual supercomputing, conceived and developed by the NAS Exploratory Computing Group. The hyperwall will allow simultaneous and coordinated visualization and interaction of an array of processes, such as a the computations of a parameter study or the parallel evolutions of a genetic algorithm population. Making over 65 million pixels available to the user, the hyperwall will enable and elicit qualitatively new ways of leveraging computers to accomplish science. It is currently still unclear whether we will be able to transport the hyperwall to SC02. The crucial display frame still has not been completed by the metal fabrication shop, although they promised an August delivery. Also, we are still working the fragile node issue, which may require transplantation of the compute nodes from the present 2U cases into 3U cases. This modification will increase the present 3-rack configuration to 5 racks.

Parks, John↗

Software for Verifying Image-Correlation Tie Points

A computer program enables assessment of the quality of tie points in the image-correlation processes of the software described in the immediately preceding article. Tie points are computed in mappings between corresponding pixels in the left and right images of a stereoscopic pair. The mappings are sometimes not perfect because image data can be noisy and parallax can cause some points to appear in one image but not the other. The present computer program relies on the availability of a left- right correlation map in addition to the usual right left correlation map. The additional map must be generated, which doubles the processing time. Such increased time can now be afforded in the data-processing pipeline, since the time for map generation is now reduced from about 60 to 3 minutes by the parallelization discussed in the previous article. Parallel cluster processing time, therefore, enabled this better science result. The first mapping is typically from a point (denoted by coordinates x,y) in the left image to a point (x',y') in the right image. The second mapping is from (x',y') in the right image to some point (x",y") in the left image. If (x,y) and(x",y") are identical, then the mapping is considered perfect. The perfect-match criterion can be relaxed by introducing an error window that admits of round-off error and a small amount of noise. The mapping procedure can be repeated until all points in each image not connected to points in the other image are eliminated, so that what remains are verified correlation data.

Klimeck, Gerhard↗

A GPU accelerated mixed-precision Smoothed Particle Hydrodynamics framework with cell-based relative coordinates

Smoothed Particle Hydrodynamics (SPH) is essential for modeling complex large-deformation problems across various applications, requiring significant computational power. A major portion of SPH computation time is dedicated to the Nearest Neighboring Particle Search (NNPS) process. While advanced NNPS algorithms have been developed to enhance SPH efficiency, the potential efficiency gains from modern computation hardware remain underexplored. Here, this study investigates the impact of GPU parallel architecture, low-precision computing on GPUs, and GPU memory management on NNPS efficiency. Our approach employs a GPU-accelerated mixed-precision SPH framework, utilizing low precision float-point 16 (FP16) for NNPS while maintaining high precision for other components. To ensure FP16 accuracy in NNPS, we introduce a Relative Coordinated-based Link List (RCLL) algorithm, storing FP16 relative coordinates of particles within background cells. Our testing results show three significant speedup rounds for CPU-based NNPS algorithms. The first comes from parallel GPU computations, with up to a 1000x efficiency gain. The second is achieved through low-precision GPU computing, where the proposed FP16-based RCLL algorithm offers a 1.5x efficiency improvement over the FP64-based approach on GPUs. By optimizing GPU memory bandwidth utilization, the efficiency of the FP16 RCLL algorithm can be further boosted by 2.7x, as demonstrated in an example with 1 million particles. Our code is released at https://github.com/pnnl/lpNNPS4SPH.

97 MATHEMATICS AND COMPUTING↗

Analysis of impingement heat transfer for two parallel liquid-metal slot jets

An analytical method is developed for determining heat transfer by impinging liquid-metal slot jets. The method involves mapping the jet flow region, which is bounded by free streamlines, into a potential plane where it becomes a uniform flow in a channel of constant width. The energy equation is transformed into potential plane coordinates and is solved in the channel flow region. Conformal mapping is then used to transform the solution back into the physical plane and obtain the desired heat-transfer characteristics. The analysis given here determines the heat-transfer characteristics for two parallel liquid-metal slot jets impinging normally against a uniformly heated flat plate. The liquid-metal assumptions are made that the jets are inviscid and that molecular conduction is dominating heat diffusion. Wall temperature distributions along the heated plate are obtained as a function of spacing between the jets and the jet Peclet number.

Siegel, R.↗

Implementation of a Parallel Kalman Filter for Stratospheric Chemical Tracer Assimilation

A Kalman filter for the assimilation of long-lived atmospheric chemical constituents has been developed for two-dimensional transport models on isentropic surfaces over the globe. An important attribute of the Kalman filter is that it calculates error covariances of the constituent fields using the tracer dynamics. Consequently, the current Kalman-filter assimilation is a five-dimensional problem (coordinates of two points and time), and it can only be handled on computers with large memory and high floating point speed. In this paper, an implementation of the Kalman filter for distributed-memory, message-passing parallel computers is discussed. Two approaches were studied: an operator decomposition and a covariance decomposition. The latter was found to be more scalable than the former, and it possesses the property that the dynamical model does not need to be parallelized, which is of considerable practical advantage. This code is currently used to assimilate constituent data retrieved by limb sounders on the Upper Atmosphere Research Satellite. Tests of the code examined the variance transport and observability properties. Aspects of the parallel implementation, some timing results, and a brief discussion of the physical results will be presented.

Chang, Lang-Ping↗