Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computing Methodologies”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantum Circuit Cutting for Classical Shadows

Classical shadow tomography is a sample-efficient technique for characterizing quantum systems and predicting many of their properties. Circuit cutting is a technique for dividing large quantum circuits into smaller fragments that can be executed more robustly using fewer quantum resources. We introduce a divide-and-conquer circuit cutting method for estimating the expectation values of observables using classical shadows. We derive a general formula for making predictions using the classical shadows of circuit fragments from arbitrarily cut circuits and provide the sample complexity analysis for the case when observables factorize across fragments. Then, we numerically show that our divide-and-conquer method outperforms traditional uncut shadow tomography when estimating high-weight observables that act non-trivially on many qubits and discuss the mechanisms for this advantage.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2016

BIGHORN Computational Fluid Dynamics Theory, Methodology, and Code Verification & Validation Benchmark Problems

This document presents the theoretical background for a hybrid finite-element / finite-volume fluid flow solver, namely BIGHORN, based on the Multiphysics Object Oriented Simulation Environment (MOOSE) computational framework developed at the Idaho National Laboratory (INL). An overview of the numerical methods used in BIGHORN are discussed and followed by a presentation of the formulation details. The document begins with the governing equations for the compressible fluid flow, with an outline of the requisite constitutive relations. A second-order finite volume method used for solving the compressible fluid flow problems is presented next. A Pressure-Corrected Implicit Continuous-fluid Eulerian (PCICE) formulation for time integration is also presented. The multi-fluid formulation is being developed. Although multi-fluid is not fully-developed, BIGHORN has been designed to handle multi-fluid problems. Due to the flexibility in the underlying MOOSE framework, BIGHORN is quite extensible, and can accommodate both multi-species and multi-phase formulations. This document also presents a suite of verification & validation benchmark test problems for BIGHORN. The intent for this suite of problems is to provide baseline comparison data that demonstrates the performance of the BIGHORN solution methods on problems that vary in complexity from laminar to turbulent flows. Wherever possible, some form of solution verification has been attempted to identify sensitivities in the solution methods, and suggest best practices when using BIGHORN.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗
DOE OSTI2018

A mixed-integer conic programming formulation for computing the flexibility index under multivariate gaussian uncertainty

We present a methodology for computing the flexibility index when uncertainty is characterized using multivariate Gaussian random variables. Our approach computes the flexibility index by solving a mixed-integer conic program (MICP). This methodology directly characterizes ellipsoidal sets to capture correlations in contrast to previous methodologies that employ approximations. We also show that, under a Gaussian representation, the flexibility index can be used to obtain a lower bound for the so-called stochastic flexibility index (i.e., the probability of having feasible operation). Our results also show that the methodology can be generalized to capture different types of uncertainty sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Crystal c-axis mapping of hcp metals by conventional reflected polarized light microscopy: Application to untextured and textured cp-Titanium

Highlights: • A computer-aided methodology based on reflected PLM was developed for c-axis orientation mapping of hcp structures. • The methodology was applied to two alpha cp-Ti bulk samples with different texture characteristics. • Calibration and validation of the methodology were performed by using EBSD. - Abstract: Nowadays the crystallographic characterization of polycrystalline metals and metallic alloys is usually performed by means of the electron backscatter diffraction (EBSD) technique. However, when dealing with coarse structures presenting intrinsic optical anisotropy, polarized light microscopy (PLM) imaging might be an alternative or a complementary time-efficient technique. Accordingly, a new computer-aided methodology has been developed in the present study and successfully applied on the crystal c-axis characterization of a hexagonal close-packed (hcp) structure, alpha commercial pure-Titanium (cp-Ti). In this sense, the correlation between EBSD-based crystallographic parameters and PLM-based parameters under certain light microscope settings has been first established by direct comparison using a representative amount of random orientations. Finally, PLM-based c-axis mapping and subsequent c-axis texture plotting have been performed and validated with EBSD for untextured and textured cp-Ti.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Scalable Energy Efficiency with Resilience for High Performance Computing Systems: A Quantitative Methodology

Energy efficiency and resilience are two crucial challenges for HPC systems to reach exascale. While energy efficiency and resilience issues have been extensively studied individually, little has been done to understand the interplay between energy efficiency and resilience for HPC systems. Decreasing the supply voltage associated with a given operating frequency for processors and other CMOS-based components can significantly reduce power consumption. However, this often raises system failure rates and consequently increases application execution time. In this work, we present an energy saving undervolting approach that leverages the mainstream resilience techniques to tolerate the increased failures caused by undervolting.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Scalable Energy Efficiency with Resilience for High Performance Computing Systems: A Quantitative Methodology

Energy efficiency and resilience are two crucial challenges for HPC systems to reach exascale. While energy efficiency and resilience issues have been extensively studied individually, little has been done to understand the interplay between energy efficiency and resilience for HPC systems. Decreasing the supply voltage associated with a given operating frequency for processors and other CMOS-based components can significantly reduce power consumption. However, this often raises system failure rates and consequently increases application execution time. In this work, we present an energy saving undervolting approach that leverages the mainstream resilience techniques to tolerate the increased failures caused by undervolting.

energy, resilience↗
DOE OSTI2015

Second-order adjoint sensitivity analysis methodology (2ND-ASAM) for computing exactly and efficiently first- and second-order sensitivities in large-scale systems - 14356

This work presents the second-order adjoint sensitivity analysis methodology (2.-ASAM) for computing exactly and efficiently the second-order functional derivatives ('sensitivities') of system responses to the system's model parameters. For a physical system comprising N{sub α} parameters and N{sub r} responses, forward methods require a total of at least (N{sup 2}{sub α}/2+3N{sub α}/2) large-scale computations for obtaining all of the first- and second-order sensitivities, for all N{sub r} system responses. On the other hand, for one functional-type system response, the 2.-ASAM requires one large-scale computation using the first-level adjoint sensitivity system (1.-LASS) for obtaining all of the first-order sensitivities, followed by at most (2N{sub α}+1) large-scale computations using the second-level adjoint sensitivity systems (2.-LASS) for obtaining all of the second-order sensitivities. The second-order sensitivities contribute decisively to causing asymmetries in the response distribution, since they are the leading contributors to the third-order response correlations (skewness). They also cause the 'expected value of the response' to differ from the 'computed nominal value of the response'. The implementation of the 2.-ASAM requires very little additional effort beyond the construction of the adjoint sensitivity system needed for computing the first-order sensitivities. Only the sources on the right-sides of the diffusion (differential) operator needed to be modified; the left-side of the differential equations (and hence the 'solver' in large-scale practical applications) remain unchanged. The application of the 2.-ASAM is illustrated on a benchmark heat conduction/convection problem, which makes transparent the underlying mathematical derivations. For this illustrative problem, 4 'large-scale' adjoint computations suffice for computing exactly all of the 6 first- and 21 distinct second-order derivatives. The 2.-ASAM presented in this work should enable the hitherto very difficult, if not intractable, exact computation of all of the second-order response sensitivities for large-systems involving many parameters, which is expected to affect significantly other fields that need efficiently computed second-order response sensitivities, such as optimization, data assimilation/adjustment, model calibration, and predictive modeling. (author)

42 ENGINEERING↗
DOE OSTI2014

Computing model independent perturbations in dark energy and modified gravity

We present a methodology for computing model independent perturbations in dark energy and modified gravity. This is done from the Lagrangian for perturbations, by showing how field content, symmetries, and physical principles are often sufficient ingredients for closing the set of perturbed fluid equations. The fluid equations close once ''equations of state for perturbations'' are identified: these are linear combinations of fluid and metric perturbations which construct gauge invariant entropy and anisotropic stress perturbations for broad classes of theories. Our main results are the proof of the equation of state for perturbations presented in a previous paper, and the development of the required calculational tools.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Computing nucleon EDM on a lattice

I will discuss briefly recent changes in the methodology of computing the baryon EDM on a lattice. The associated correction substantially reduces presently existing lattice values for the proton and neutron theta-induced EDMs, so that even the most precise previous lattice results become consistent with zero. On one hand, this change removes previous disagreements between these lattice results and the phenomenological estimates of the nucleon EDM. On the other hand, the nucleon EDM becomes much harder to compute on a lattice. In addition, I will review the progress in computing quark chromo-EDM-induced nucleon EDM using chiral quark action.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

QCOR; A Language Extension Specification for the Heterogeneous Quantum-Classical Model of Computation

Quantum computing (QC) is an emerging computational paradigm that leverages the laws of quantum mechanics to perform elementary logic operations. Existing programming models for QC were designed with fault-tolerant hardware in mind, envisioning stand-alone applications. However, the susceptibility of near-term quantum computers to noise limits their stand-alone utility. To better leverage limited computational strengths of noisy quantum devices, hybrid algorithms have been suggested whereby quantum computers are used in tandem with their classical counterparts in a heterogeneous fashion. This modus operandi calls out for a programming model and a high-level programming language that natively and seamlessly supports heterogeneous quantum-classical hardware architectures in a single-source-code paradigm. Motivated by the lack of such a model, we introduce a language extension specification, called QCOR, which enables single-source quantum-classical programming. Programs written using the QCOR library–based language extensions can be compiled to produce functional hybrid binary executables. After defining QCOR’s programming model, memory model, and execution model, we discuss how QCOR enables variational, iterative, and feed-forward QC. Additionally, QCOR approaches quantum-classical computation in a hardware-agnostic heterogeneous fashion and strives to build on best practices of high-performance computing. The high level of abstraction in the language extension is intended to accelerate the adoption of QC by researchers familiar with classical high-performance computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Effects of Second-Order Hydrodynamic Forces on Floating Offshore Wind Turbines

Relative to first-order, second-order wave-excitation loads are known to cause significant motions and additional loads in offshore oil and gas platforms. The design of floating offshore wind turbines was partially inherited from the offshore oil and gas industry. Floating offshore wind concepts have been studied with powerful aero-hydro-servo-elastic tools; however, most of the existing work on floating offshore wind turbines has neglected the contribution of second-order wave-excitation loads. As a result, this paper presents a computationally efficient methodology to consider these loads within FAST, a wind turbine computer-aided engineering tool developed by the National Renewable Energy Laboratory. The method implemented was verified against the commercial OrcaFlex tool, with good agreement, and low computational time. A reference floating offshore wind turbine was studied under several wind and wave load conditions, including the effects of second-order slow-drift and sum-frequency loads. Preliminary results revealed that these loads excite the turbine's natural frequencies, namely the surge and pitch natural frequencies.

13 HYDRO ENERGY↗
DOE OSTI2020

Multiscale modeling of microwave-heated multiphase systems

Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗