Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partitioned algorithm”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

PUMIPic: A mesh-based approach to unstructured mesh Particle-In-Cell on GPUs

Unstructured mesh particle-in-cell, PIC, simulations executing on the current and next generation of massively parallel systems require new methods for both the mesh and particles to achieve performance and scalability on GPUs. The traditional approach to implementing PIC simulations defines data structures and algorithms in terms of particles with a full copy of the unstructured mesh on every process. To effectively scale the unstructured mesh and particles, mesh-based PIC uses the unstructured mesh as the predominant data structure with the particles stored in terms of the mesh entities. Here, this paper details the PUMIPic library, a framework for developing efficient and performance-portable mesh-based PIC simulations on GPU systems. A pseudo physics simulation based on a five-dimensional gyro-kinetic code for modeling plasma physics is used to examine the performance of PUMIPic. Scaling studies of the unstructured mesh partition and number of particles are performed up to 4096 nodes of the Summit system at Oak Ridge National Laboratory. The studies show that mesh-based PIC can utilize a partitioned mesh and maintain scaling up to system limitations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Understanding the Scalability of Bayesian Network Inference using Clique Tree Growth Curves

Bayesian networks (BNs) are used to represent and efficiently compute with multi-variate probability distributions in a wide range of disciplines. One of the main approaches to perform computation in BNs is clique tree clustering and propagation. In this approach, BN computation consists of propagation in a clique tree compiled from a Bayesian network. There is a lack of understanding of how clique tree computation time, and BN computation time in more general, depends on variations in BN size and structure. On the one hand, complexity results tell us that many interesting BN queries are NP-hard or worse to answer, and it is not hard to find application BNs where the clique tree approach in practice cannot be used. On the other hand, it is well-known that tree-structured BNs can be used to answer probabilistic queries in polynomial time. In this article, we develop an approach to characterizing clique tree growth as a function of parameters that can be computed in polynomial time from BNs, specifically: (i) the ratio of the number of a BN's non-root nodes to the number of root nodes, or (ii) the expected number of moral edges in their moral graphs. Our approach is based on combining analytical and experimental results. Analytically, we partition the set of cliques in a clique tree into different sets, and introduce a growth curve for each set. For the special case of bipartite BNs, we consequently have two growth curves, a mixed clique growth curve and a root clique growth curve. In experiments, we systematically increase the degree of the root nodes in bipartite Bayesian networks, and find that root clique growth is well-approximated by Gompertz growth curves. It is believed that this research improves the understanding of the scaling behavior of clique tree clustering, provides a foundation for benchmarking and developing improved BN inference and machine learning algorithms, and presents an aid for analytical trade-off studies of clique tree clustering using growth curves.

Mengshoel, Ole Jakob↗
NASA NTRS1990

Digital halftoning methods for selectively partitioning error into achromatic and chromatic channels

A method is described for reducing the visibility of artifacts arising in the display of quantized color images on CRT displays. The method is based on the differential spatial sensitivity of the human visual system to chromatic and achromatic modulations. Because the visual system has the highest spatial and temporal acuity for the luminance component of an image, a technique which will reduce luminance artifacts at the expense of introducing high-frequency chromatic errors is sought. A method based on controlling the correlations between the quantization errors in the individual phosphor images is explored. The luminance component is greatest when the phosphor errors are positively correlated, and is minimized when the phosphor errors are negatively correlated. The greatest effect of the correlation is obtained when the intensity quantization step sizes of the individual phosphors have equal luminances. For the ordered dither algorithm, a version of the method can be implemented by simply inverting the matrix of thresholds for one of the color components.

Mulligan, Jeffrey B.↗
DOE OSTI2023

Application of the variational autoencoder to detect the critical points of the anisotropic Ising model

We generalize the previous study on the application of variational autoencoders to the two-dimensional Ising model to a system with anisotropy. Due to the self-duality property of the system, the critical points can be located exactly for the entire range of anisotropic coupling. This presents an excellent test bed for the validity of using a variational autoencoder to characterize an anisotropic classical model. Furthermore, we reproduce the phase diagram for a wide range of anisotropic couplings and temperatures via a variational autoencoder without the explicit construction of an order parameter. Considering that the partition function of ($d$ + 1)-dimensional anisotropic models can be mapped to that of the $d$-dimensional quantum spin models, the present study provides numerical evidence that a variational autoencoder can be applied to analyze quantum systems via the quantum Monte Carlo method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Distributionally Robust Decentralized Volt-Var Control With Network Reconfiguration

Here, this paper presents a decentralized volt-var optimization (VVO) and network reconfiguration strategy to address the challenges arising from the growing integration of distributed energy resources, particularly photovoltaic (PV) generation units, in active distribution networks. To reconcile control measures with different time resolutions and empower local control centers to handle intermittency locally, the proposed approach leverages a two-stage distributionally robust optimization; decisions on slow-responding control measures and set points that link neighboring subnetworks are made in advance while considering all plausible distributions of uncertain PV outputs. We present a decomposition algorithm with an acceleration scheme for solving the proposed model. Numerical experiments on the IEEE 123 bus distribution system are given to demonstrate its outstanding out-of-sample performance and computational efficiency, which suggests that the proposed method can effectively localize uncertainty via risk-informed proactive timely decisions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2006

Observing Ice in Clouds from Space

There are many satellite observations of cloud top properties and the liquid and rain content of clouds, however, we do not yet quantitatively understand the processes that control the water budget of the upper troposphere where ice is the predominant phase, and how these processes are linked to precipitation processes and the radiative energy budget. The ice in clouds in the upper troposphere either melts into rain or is detrained, and persists, as cirrus clouds affecting the hydrological and energy cycle, respectively. Fully modeling the Earth's climate and improving weather and climate forecasts requires accurate satellite measurements of various cloud properties at the temporal and spatial scales of cloud processes. These properties include cloud horizontal and vertical structure, cloud water content and some measure of particle sizes and shapes. The uncertainty in knowledge of these ice characteristics is reflected in the large discrepancies in model simulations of the upper tropospheric water budget. Model simulations are sensitive to the partition of ice between precipitation and outflow processes, i.e., to the parameterization of ice clouds and ice processes. One barrier to achieving accurate global ice cloud properties is the lack of adequate observations at millimeter and submillimeter wavelengths (183-874 GHz). Recent advances in instrumentation have allowed for the development and implementation of an airborne submillimeter-wave radiometer. The brightness temperatures at these frequencies are especially sensitive to cirrus ice particle sizes (because they are comparable to the wavelength). This allows for more accurate ice water path estimates when multiple channels are used to probe into the cloud layers. Further, submillimeter wavelengths offer simplicity in the retrieval algorithms because they do not probe into the liquid and near surface portions of clouds, thus requiring only one term of the radiative transfer equation (ice scattering) to relate brightness temperatures to ice. The next step is a satellite mission designed to acquire global Earth radiance measurements in the submillimeter-wave region, thus bridging the measurement gap between microwave sounders and shorter-wavelength infrared and visible sensors. This presentation provides scientific justification and an approach to measuring ice water path and particle size from a satellite platform that spans a range encompassing both the hydrologically active and radiatively active components of cloud systems.

Ackerman, S.↗
NASA NTRS2026

PERSIANN-Unet: A Global Deep Learning Framework for Near-Real-Time Precipitation Estimation Using Infrared Data

Access to high-quality, high-resolution, near-real-time precipitation data is essential for hydrological and meteorological research and disaster mitigation. Traditional tools such as rain gauges and radar networks, though effective, have limitations, including sparse coverage in remote areas and high operational costs. Satellite data, with its global coverage and high spatial and temporal resolutions, mitigates limitations in coverage. Satellite precipitation products like Hydro Estimator (HE), Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG), and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN) utilize both geosynchronous thermal infrared (IR) and passive microwave (PMW) data in their operation. PMW sensors offer detailed atmospheric profiles but suffer from higher latency, whereas IR sensors provide lower latency but only capture cloud-top information. Despite this constraint, IR data remains attractive for low-latency precipitation estimation. Recent advances in deep learning, particularly convolutional neural networks (CNNs), have further improved satellite precipitation retrievals. This study introduces PERSIANN-Unet (PUnet or PERSIANN V3), a quasi-global algorithm covering 60°N–60°S that combines IR data, monthly climatology, and the UNet architecture to produce half-hourly precipitation estimates at 0.04° resolution. The product is evaluated against HE, IMERG, and PDIR-Now for 2022–2023. Results show that PUnet closely matches its training target, IMERG V07 Final, at the global scale, and performance is further evaluated against Stage IV as a reference over CONUS. Training PUnet on IMERG (2016–2021) leverages a high-quality, integrated PMW IR-gauge precipitation product while developing an IR-based framework not reliant on PMW availability. By operating on a single global image, PUnet avoids tile partitioning and blending steps, reducing edge discontinuities, and produces more spatially consistent precipitation fields across hemispheres.

Phu Nguyen↗
DOE OSTI2025

Fragme∩t: An Open‐Source Framework for Multiscale Quantum Chemistry Based on Fragmentation

Fragment-based quantum chemistry offers a means to circumvent the nonlinear computational scaling of conventional electronic structure calculations, by partitioning a large calculation into smaller subsystems then considering the many-body interactions between them. Variants of this approach have been used to parameterize classical force fields and machine learning potentials, applications that benefit from interoperability between quantum chemistry codes. However, there is a dearth of software that provides interoperability yet is purpose-built to handle the combinatorial complexity of fragment-based calculations. To fill this void we introduce “Fragme∩t”, an open-source software application that provides a tool for community validation of fragment-based methods, a platform for developing new approximations, and a framework for analyzing many-body interactions. Fragme∩t includes algorithms for automatic fragment generation and structure modification, and for distance- and energy-based screening of the requisite subsystems. Checkpointing, database management, and parallelization are handled internally and results are archived in a portable database. Interfaces to various quantum chemistry engines are easy to write and exist already for Q-Chem, PySCF, xTB, Orca, CP2K, MRCC, Psi4, NWChem, GAMESS, and MOPAC. Applications reported here demonstrate parallel efficiencies around 96% on more than 1000 processors but also showcase that the code can handle large-scale protein fragmentation using only workstation hardware, all with a codebase that is designed to be usable by non-experts. Fragme∩t conforms to modern software engineering best practices and is built upon well established technologies including Python, SQLite, and Ray. The source code is available under the Apache 2.0 license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Rule groupings in expert systems using nearest neighbour decision rules, and convex hulls

Expert System shells are lacking in many areas of software engineering. Large rule based systems are not semantically comprehensible, difficult to debug, and impossible to modify or validate. Partitioning a set of rules found in CLIPS (C Language Integrated Production System) into groups of rules which reflect the underlying semantic subdomains of the problem, will address adequately the concerns stated above. Techniques are introduced to structure a CLIPS rule base into groups of rules that inherently have common semantic information. The concepts involved are imported from the field of A.I., Pattern Recognition, and Statistical Inference. Techniques focus on the areas of feature selection, classification, and a criteria of how 'good' the classification technique is, based on Bayesian Decision Theory. A variety of distance metrics are discussed for measuring the 'closeness' of CLIPS rules and various Nearest Neighbor classification algorithms are described based on the above metric.

Anastasiadis, Stergios↗
NASA NTRS1989

JPL control/structure interaction test bed real-time control computer architecture

The Control/Structure Interaction Program is a technology development program for spacecraft that exhibit interactions between the control system and structural dynamics. The program objectives include development and verification of new design concepts - such as active structure - and new tools - such as combined structure and control optimization algorithm - and their verification in ground and possibly flight test. A focus mission spacecraft was designed based upon a space interferometer and is the basis for design of the ground test article. The ground test bed objectives include verification of the spacecraft design concepts, the active structure elements and certain design tools such as the new combined structures and controls optimization tool. In anticipation of CSI technology flight experiments, the test bed control electronics must emulate the computation capacity and control architectures of space qualifiable systems as well as the command and control networks that will be used to connect investigators with the flight experiment hardware. The Test Bed facility electronics were functionally partitioned into three units: a laboratory data acquisition system for structural parameter identification and performance verification; an experiment supervisory computer to oversee the experiment, monitor the environmental parameters and perform data logging; and a multilevel real-time control computing system. The design of the Test Bed electronics is presented along with hardware and software component descriptions. The system should break new ground in experimental control electronics and is of interest to anyone working in the verification of control concepts for large structures.

Briggs, Hugh C.↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2021

Seawater Debye Model Function at L-Band and Its Impact on Salinity Retrieval From Aquarius Satellite Data

A model function of seawater, which specifies the dielectric constant of seawater as a function of salinity, temperature, and frequency, is important for the retrieval of sea surface salinity using satellite data. In 2017, a model function has been developed based on measurement data at 1.4134 GHz using a third-order polynomial expression in salinity ( S ) and temperature ( T ). Although the model showed improvements in salinity retrieval, it had an inconsistent behavior between partitioned salinities. To improve the stability of the model, new dielectric measurements of seawater have been made recently over a broad range of salinities and temperatures to expand the data set used for developing the model function. The structure of the model function has been changed from a polynomial expansion in S and T to a physics-based model consisting of a Debye molecular resonance term plus a conductivity term. Each unknown parameter is expressed in S and T based on the expanded measurement data set. Physical arguments have been used to limit the number of unknown coefficients in these expressions to improve the stability of the model function. The new model function has been employed in the retrieval algorithm of the Aquarius satellite mission to obtain a global salinity map. The retrieved salinity using a different model function is compared with in situ data collected by Argo floats to evaluate the impact and the performance of model functions. The results indicate that the new model function has significant improvements in salinity retrieval compared with other existing models.

Seawater↗
NASA NTRS2024

Accurate Infrared Line Lists for 20 Isotopologues of CS 2 at Room Temperature

To facilitate atmospheric and spectroscopic studies of carbon disulfide, or CS 2 , in both planetary and exoplanetary atmospheres, we adopt the “Best Theory + Reliable High-resolution Experiment” algorithm to generate semiempirical IR line lists for the 20 most abundant CS 2 isotopologues, denoted as Ames-296K. The IR lists are computed using the Ames-1 potential energy surface, refined using the experimental transition set and an ab initio dipole moment surface fitted from CCSD(T)/aug-cc-pV(T/Q/5+d)Z dipoles extrapolated to a one-particle basis set limit. The IR lists cover the range of 0–10,000 cm -1 , with an S 296K cutoff at 10 -31 cm -1 /molecule·cm -2 (abundance included). A “natural” IR line list at 296 K includes about 10 million lines of the 20 isotopologues, with their intensities scaled by the corresponding abundances. The zero-point energy, partition functions, and abundances are reported for each isotopologue. The energy levels in the global effective Hamiltonian model for 12 C 32 S 2 are adopted to improve the line position accuracy. This new IR list for the main isotopologue is denoted as A+I.296K. Reliable HITRAN2020 line positions are also utilized to improve the accuracy of the 32 S 12 C 34 S, 32 S 12 C 33 S, and 32 S 13 C 32 S isotopologue line lists. The final composite line list is validated against Pacific Northwest National Laboratory experimental cross sections, showing excellent agreement. The agreement supports the quality of the composite line list and the power of synergy between experiment and theory. The new data are proposed for use in updating and expanding the CS 2 data in HITRAN and other high-resolution IR databases. Supplementary files are available in Zenodo and AHED.

Exoplanet atmospheric composition↗
NASA NTRS2017

Applying Graph Theory to Problems in Air Traffic Management

Graph theory is used to investigate three different problems arising in air traffic management. First, using a polynomial reduction from a graph partitioning problem, it isshown that both the airspace sectorization problem and its incremental counterpart, the sector combination problem are NP-hard, in general, under several simple workload models. Second, using a polynomial time reduction from maximum independent set in graphs, it is shown that for any fixed e, the problem of finding a solution to the minimum delay scheduling problem in traffic flow management that is guaranteed to be within n1-e of the optimal, where n is the number of aircraft in the problem instance, is NP-hard. Finally, a problem arising in precision arrival scheduling is formulated and solved using graph reachability. These results demonstrate that graph theory provides a powerful framework for modeling, reasoning about, and devising algorithmic solutions to diverse problems arising in air traffic management.

computational complexity↗
NASA NTRS2006

A Mission to Observe Ice in Clouds from Space

To date there have been multiple satellite missions to observe and retrieve cloud top properties and the liquid in, and precipitation from, clouds. There are currently a few missions that attempt to measure cloud ice properties as a byproduct of other observations. However, we do not yet quantitatively understand the processes that control the water budget of the upper troposphere where ice is the predominant phase, and how these processes are linked to precipitation processes and the radiative energy budget. The ice in clouds either melts into rain or is detrained, and persists, as cirrus clouds affecting the hydrological and energy cycle, respectively. Fully modeling the Earth's climate and improving weather and climate forecasts requires accurate satellite measurements of various cloud properties at the temporal and spatial scales of cloud processes. The uncertainty in knowledge of these ice characteristics is reflected in the large discrepancies in model simulations of the upper tropospheric water budget. Model simulations are sensitive to the partition of ice between precipitation and outflow processes, i.e., to the parameterization of ice clouds and ice processes. This presentation will describe the Submillimeter-wave InfraRed Ice Cloud Experiment (SIRICE) concept, a satellite mission designed to acquire global Earth radiance measurements in the infrared and submillimeter-wave region (183-874 GHz). If successful, this mission will bridge the measurement gap between microwave sounders and shorter-wavelength infrared and visible sensors. The brightness temperatures at submillimeter-wave frequencies are especially sensitive to cirrus ice particle sizes (because they are comparable to the wavelength). This allows for more accurate ice water path estimates when multiple channels are used to probe into the cloud layers. Further, submillimeter wavelengths offer simplicity in the retrieval algorithms because they do not probe into the liquid and near surface portions of clouds, thus requiring only one term of the radiative transfer equation (ice scattering) to relate brightness temperatures to ice. Scientific justification and the SIRICE approach to measuring ice water path and particle size that span a range encompassing both the hydrologically active and radiatively active components of cloud systems will be presented.

Ackerman, S.↗
NASA NTRS2017

Applying Graph Theory to Problems in Air Traffic Management

Graph theory is used to investigate three different problems arising in air traffic management. First, using a polynomial reduction from a graph partitioning problem, it is shown that both the airspace sectorization problem and its incremental counterpart, the sector combination problem are NP-hard, in general, under several simple workload models. Second, using a polynomial time reduction from maximum independent set in graphs, it is shown that for any fixed e, the problem of finding a solution to the minimum delay scheduling problem in traffic flow management that is guaranteed to be within n1-e of the optimal, where n is the number of aircraft in the problem instance, is NP-hard. Finally, a problem arising in precision arrival scheduling is formulated and solved using graph reachability. These results demonstrate that graph theory provides a powerful framework for modeling, reasoning about, and devising algorithmic solutions to diverse problems arising in air traffic management.

precision arrival scheduling↗
NASA NTRS2005

Studying the Diurnal Cycle of Convection Using a TRMM-Calibrated Infrared Rain Algorithm

The development of a satellite infrared (IR) technique for estimating convective and stratiform rainfall and its application in studying the diurnal variability of rainfall on a global scale is presented. The Convective-Stratiform Technique (CST), calibrated by coincident, physically retrieved rain rates from the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR), is applied over the global tropics. The technique makes use of the IR data from the TRMM Visible/Infrared Scanner (VIRS) before application to global geosynchronous satellite data. The calibrated CST technique has the advantages of high spatial resolution (4 km), filtering of nonraining cirrus clouds, and the stratification of the rainfall into its convective and stratiform components, the last being important for the calculation of vertical profiles of latent heating. The diurnal cycle of rainfall, as well as the division between convective and Stratiform rainfall will be presented. The technique is validated using available data sets and compared to other global rainfall products such as Global Precipitation Climatology Project (GPCP) IR product, calibrated with TRMM Microwave Imager (TMI) data. Results from five years of PR data will show the global-tropical partitioning of convective and stratiform rainfall.

Negri, Andrew J.↗