Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel data structures”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Interpretation of the northern boundary of Ishtar Terra from Magellan images and altimetry

Part of the controversy on the origin of western Ishtar Terra (IT) concerns the nature of Uorsar Rupes (UR), the northern boundary of IT. In the hypothesis of lithospheric convergence and underthrusting, UR is held to be the main boundary thrust fault at the toe of an accretionary wedge. A topographic rise parallel to the scarp was interpreted as a flexural bulge similar to those of terrestrial subduction zones, and quantitative models of this feature seemed broadly consistent with the expected lithospheric structure of Venus. In the alternative mantle upwelling hypothesis for western IT, the outer margins of the highland are thought to be collapsing, and UR has been interpreted as a normal fault. Herein, Magellan images and altimetry are interpreted for this region and the hypothesis that a flexural signature can be distinguished is reassessed. The Magellan images of IT show evidence of crustal shortening adjacent to UR, but extension and burial dominate northwards. Altimetric profiles display the same long wavelength trends visible in Venera data, but no clear evidence of the lithospheric flexure. A model of regional extension and burial is herein favored, but regional compression cannot be ruled out.

Mueller, S.↗
DOE OSTI2020

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↗
DOE OSTI2026

Spherical tokamak physics research in preparation for the operation of NSTX-U

The National Spherical Torus Experiment Upgrade (NSTX-U) is preparing to resume operation, representing a crucial step toward realizing compact, cost-effective fusion pilot plants. In advance of this, extensive modeling and data analysis have been conducted to advance the physics basis for low-aspect-ratio, high-performance plasma regimes, focusing on three core objectives: confinement and stability, power and particle handling, and steady-state operation. Significant progress has been made in understanding the electron temperature flattening in high-β plasmas, which is shown to be driven by a complex interplay of magnetohydrodynamic instabilities (e.g. non-resonant infernal modes), fast-ion-driven Alfvén eigenmodes, and electron and ion-scale micro-instabilities, particularly Kinetic Ballooning Modes (KBMs), whose destabilization is strongly dependent on parallel magnetic field fluctuations (δB ∥ ). Furthermore, a new gyrokinetic critical pedestal model was developed, accurately predicting pedestal structure by identifying KBMs as the primary stability limit, offering a critical constraint for future high-confinement scenarios. To address the challenge of high heat flux, novel liquid lithium plasma-facing components were modeled. The analysis confirmed that lithium vapor shielding is a self-regulating mechanism for heat mitigation, while also emphasizing that strong main ion parallel flow is essential to minimize core lithium contamination. Finally, progress toward steady-state operation was anchored by developing the required physics basis and control tools. This includes predictive modeling for reversed magnetic shear sustainment, demonstrating that magnetic island-induced bootstrap current reduction is negligible in STs, and advancing real-time control and disruption avoidance capabilities. The development of high-speed surrogate models (e.g. MMMNet) provides computationally efficient tools vital for non-inductive scenario optimization and integrated, low-disruptivity operations planned for NSTX-U.

NSTX-U↗
DOE OSTI2024

Synthetic β-sheets mimicking fibrillar and oligomeric structures for evaluation of spectral X-ray scattering technique for biomarker quantification

Archetypical cross-β spines sharpen the boundary between functional and pathological proteins including β-amyloid, tau, α-synuclein and transthyretin are linked to many debilitating human neurodegenerative and non-neurodegenerative amyloidoses. An increased focus on development of pathogenic β-sheet specific fluid and imaging structural biomarkers and conformation-specific monoclonal antibodies in targeted therapies has been recently observed. Identification and quantification of pathogenic oligomers remain challenging for existing neuroimaging modalities. We propose two artificial β-sheets which can mimic the nanoscopic structural characteristics of pathogenic oligomers and fibrils for evaluating the performance of a label free, X-ray based biomarker detection and quantification technique. Highly similar structure with elliptical cross-section and parallel cross-β motif is observed among recombinant α-synuclein fibril, Aβ-42 fibril and artificial β-sheet fibrils. We then use these β-sheet models to assess the performance of spectral small angle X-ray scattering (sSAXS) technique for detecting β-sheet structures. sSAXS showed quantitatively accurate detection of antiparallel, cross-β artificial oligomers from a tissue mimicking environment and significant distinction between different oligomer packing densities such as diffuse and dense packings. The proposed synthetic β-sheet models mimicked the nanoscopic structural characteristics of β-sheets of fibrillar and oligomeric states of Aβ and α-synuclein based on the ATR-FTIR and SAXS data. The tunability of β-sheet proportions and shapes of structural motifs, and the low-cost of these β-sheet models can become useful test materials for evaluating β-sheet or amyloid specific biomarkers in a wide range of neurological diseases. By using the proposed synthetic β-sheet models, our study indicates that the sSAXS has potential to evaluate different stages of β-sheet-enriched structures including oligomers of pathogenic proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Block-structured, equal-workload, multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Total Variation Diminishing)

Abstract. We describe the development of a block-structured, equal-CPU-load (central processing unit), multi-grid-nesting interface for the Boussinesq wave model FUNWAVE-TVD (Fully Nonlinear Boussinesq Wave Model with Total Variation Diminishing Solver). The new model framework does not interfere with the core solver, and thus the core program, FUNWAVE-TVD, is still a standalone model used for a single grid. The nesting interface manages the time sequencing and two-way nesting processes between the parent grid and child grid with grid refinement in a hierarchical manner. Workload balance in the MPI-based (message passing interface) parallelization is handled by an equal-load scheme. A strategy of shared array allocation is applied for data management that allows for a large number of nested grids without creating additional memory allocations. Four model tests are conducted to verify the nesting algorithm with assessments of model accuracy and the robustness in the application in modeling transoceanic tsunamis and coastal effects.

Choi, Young-Kwang↗
NASA NTRS1986

Subcritical and supercritical interplanetary shocks - Magnetic field and energetic particle observations

A study of 34 forward interplanetary shocks observed by ISEE 3 during 1978 and 1979 has been conducted. Magnetic field and high-energy particle data have been used, and for each shock the first critical Mach number has been determined. The first surprising result is that the majority of the observed shocks appear to be supercritical, and consistent with their supercritical character, many shocks have a foot and/or an overshoot in the magnetic field structure. Large-amplitude low-frequency waves (period of about 20 s in the spacecraft frame) are commonly observed upstream of all supercritical shocks (except for a few quasi-perpendicular shocks) and also upstream of the few subcritical shocks. Intense particle events are frequently observed at many shocks: spikes at quasi-perpendicular shocks and energetic storm particle events associated with quasi-parallel shocks can be comparably intense. The correlation of the high-energy particle peak flux with various shock parameters is in agreement with the acceleration mechanisms proposed by previous studies.

Bavassano-Cattaneo, M. B.↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2022

Co-flow injection for serial crystallography at X-ray free-electron lasers

Serial femtosecond crystallography (SFX) is a powerful technique that exploits X-ray free-electron lasers to determine the structure of macromolecules at room temperature. Despite the impressive exposition of structural details with this novel crystallographic approach, the methods currently available to introduce crystals into the path of the X-ray beam sometimes exhibit serious drawbacks. Samples requiring liquid injection of crystal slurries consume large quantities of crystals (at times up to a gram of protein per data set), may not be compatible with vacuum configurations on beamlines or provide a high background due to additional sheathing liquids present during the injection. Proposed and characterized here is the use of an immiscible inert oil phase to supplement the flow of sample in a hybrid microfluidic 3D-printed co-flow device. Co-flow generation is reported with sample and oil phases flowing in parallel, resulting in stable injection conditions for two different resin materials experimentally. A numerical model is presented that adequately predicts these flow-rate conditions. The co-flow generating devices reduce crystal clogging effects, have the potential to conserve protein crystal samples up to 95% and will allow degradation-free light-induced time-resolved SFX.

3D printing↗
NASA NTRS1993

Analytic model of aurorally coupled magnetospheric and ionospheric electrostatic potentials

This paper describes modest but significant improvements on earlier studies of electrostatic potential structure in the auroral region, using the adiabatic auroral arc model. With certain simplifying assumptions, new analytic nonlinear solutions fully exhibiting the parametric dependence of potentials on magnetospheric (e.g., cross-tail potential) and ionospheric (e.g., recombination rate) parameters are found. No purely phenomenological parameters are introduced. The results are in reasonable agreement with observed average auroral potential drops, inverted-V scale sizes, and dissipation rates. The dissipation rate is quite comparable to tail energization and transport rates and should have a major effect on tail and magnetospheric dynamics. Various relations between the cross-tail potential and auroral parameters (e.g., total parallel currents and potential drops) are given which can be studied with existing data sets.

Cornwall, John M.↗
NASA NTRS2000

Numerical Propulsion System Simulation: An Overview

The cost of implementing new technology in aerospace propulsion systems is becoming prohibitively expensive and time consuming. One of the main contributors to the high cost and lengthy time is the need to perform many large-scale hardware tests and the inability to integrate all appropriate subsystems early in the design process. The NASA Glenn Research Center is developing the technologies required to enable simulations of full aerospace propulsion systems in sufficient detail to resolve critical design issues early in the design process before hardware is built. This concept, called the Numerical Propulsion System Simulation (NPSS), is focused on the integration of multiple disciplines such as aerodynamics, structures and heat transfer with computing and communication technologies to capture complex physical processes in a timely and cost-effective manner. The vision for NPSS, as illustrated, is to be a "numerical test cell" that enables full engine simulation overnight on cost-effective computing platforms. There are several key elements within NPSS that are required to achieve this capability: 1) clear data interfaces through the development and/or use of data exchange standards, 2) modular and flexible program construction through the use of object-oriented programming, 3) integrated multiple fidelity analysis (zooming) techniques that capture the appropriate physics at the appropriate fidelity for the engine systems, 4) multidisciplinary coupling techniques and finally 5) high performance parallel and distributed computing. The current state of development in these five area focuses on air breathing gas turbine engines and is reported in this paper. However, many of the technologies are generic and can be readily applied to rocket based systems and combined cycles currently being considered for low-cost access-to-space applications. Recent accomplishments include: (1) the development of an industry-standard engine cycle analysis program and plug 'n play architecture, called NPSS Version 1, (2) A full engine simulation that combines a 3D low-pressure subsystem with a 0D high pressure core simulation. This demonstrates the ability to integrate analyses at different levels of detail and to aerodynamically couple components, the fan/booster and low-pressure turbine, through a 3D computational fluid dynamics simulation. (3) Simulation of all of the turbomachinery in a modern turbofan engine on parallel computing platform for rapid and cost-effective execution. This capability can also be used to generate full compressor map, requiring both design and off-design simulation. (4) Three levels of coupling characterize the multidisciplinary analysis under NPSS: loosely coupled, process coupled and tightly coupled. The loosely coupled and process coupled approaches require a common geometry definition to link CAD to analysis tools. The tightly coupled approach is currently validating the use of arbitrary Lagrangian/Eulerian formulation for rotating turbomachinery. The validation includes both centrifugal and axial compression systems. The results of the validation will be reported in the paper. (5) The demonstration of significant computing cost/performance reduction for turbine engine applications using PC clusters. The NPSS Project is supported under the NASA High Performance Computing and Communications Program.

Lytle, John K.↗
DOE OSTI2017

Insight into 3D micro-CT data: exploring segmentation algorithms through performance metrics

Three-dimensional (3D) micro-tomography (µ-CT) has proven to be an important imaging modality in industry and scientific domains. Understanding the properties of material structure and behavior has produced many scientific advances. An important component of the 3D µ-CT pipeline is image partitioning (or image segmentation), a step that is used to separate various phases or components in an image. Image partitioning schemes require specific rules for different scientific fields, but a common strategy consists of devising metrics to quantify performance and accuracy. The present article proposes a set of protocols to systematically analyze and compare the results of unsupervised classification methods used for segmentation of synchrotron-based data. Additionally, the proposed dataflow for Materials Segmentation and Metrics (MSM) provides 3D micro-tomography image segmentation algorithms, such as statistical region merging (SRM), k-means algorithm and parallel Markov random field (PMRF), while offering different metrics to evaluate segmentation quality, confidence and conformity with standards. Both experimental and synthetic data are assessed, illustrating quantitative results through the MSM dashboard, which can return sample information such as media porosity and permeability. The main contributions of this work are: (i) to deliver tools to improve material design and quality control; (ii) to provide datasets for benchmarking and reproducibility; (iii) to yield good practices in the absence of standards or ground-truth for ceramic composite analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A Communication-Avoiding 3D LU Factorization Algorithm for Sparse Matrices

We propose a new algorithm to improve the strong scalability of right-looking sparse LU factorization on distributed memory systems. Our 3D sparse LU algorithm uses a three-dimensional MPI process grid, aggressively exploits elimination tree parallelism and trades off increased memory for reduced per-process communication. We also analyze the asymptotic improvements for planar graphs (e.g., from 2D grid or mesh domains) and certain non-planar graphs (specifically for 3D grids and meshes). For planar graphs with n vertices, our algorithm reduces communication volume asymptotically in n by a factor of O{log n} and latency by a factor of O{log n}. For non-planar cases, our algorithm can reduce the per-process communication volume by 3× and latency by O{n^1/3} times. In all cases, the memory needed to achieve these gains is a constant factor. We implemented our algorithm by extending the 2D data structure used in superLU. Our new 3D code achieves speedups up to 27× for planar graphs and up to 3.3× for non-planar graphs over the baseline 2D superLU when run on 24,000 cores of a Cray XC30.

Sao, Piyush↗
NASA NTRS1981

Investigation from Japanese MAGSAT Team. Part A. Crustal structure near Japan and in Antarctic station. Part B. Electric currents and hydromagnetic waves in the ionosphere and the magnetosphere

Preliminary results of MAGSAT data analysis are described. Regional anomaly maps (deviations from the MGST model field) for X,Y,Z, and F in the area of 115 to 155 deg E and 20 to 60 deg N were obtained. A similar map for the geomagnetic total force anomaly in the vicinity of Japan showed that the observed anomaly can be explained by the difference in crustal magnetization between the Japan Sea and the Japan Island, which reflects a difference of 25 km in the thickness of the magnetized layer. The MAGSAT record of a sudden commencement of a magnetic storm above the South Atlantic Ocean showed a reverse impulse particularly in the D-component. Results relating to toroidal currents in the ionosphere, transverse and parallel perturbations over the polar regions, the relationship between field aligned currents and precipitating electrons, and the calculation of the subsatellite electric field are also discussed.

Fukushima, N.↗
NASA NTRS1973

Structural lineaments of Gaspe from ERTS imagery

A test study was conducted to assess the value of ERTS images for mapping geologic features of the Gaspe Peninsula, Quebec. The specific objectives of the study were: 1) to ascertain the best procedure to follow in order to obtain valuable geologic data as a result of interpretation; and 2) to indicate in which way these data could relate to mineral exploration. Of the four spectral bands of the Multispectral scanner, the band from 700 to 800 nanometers, which seems to possess the best informational content for geologic study, was selected for analysis. The original ERTS image at a scale of 1:3,700,000 was enlarged about 15 times and reproduced on film. Geologically meaningful lines, called structural lineaments, were outlined and classified according to five categories: morpho-lithologic boundaries, morpho-lithologic lineaments, fault traces, fracture zones and undefined lineaments. Comparison with the geologic map of Gaspe shows that morpho-lithologic boundaries correspond to contacts between regional stratigraphic units. Morpholithologic lineaments follow bedding trends, whereas fracture traces appear as sets of parallel lineaments, intersecting at high angles the previous category of lineaments. Fault traces mark more precisely the location of faults already mapped and spot the presence of presumable faults, not indicated on the geologic map.

Steffensen, R.↗
NASA NTRS1988

Interiors of the giant planets

This theoretical/observational project constrains interior structure of Jovian planets through observational data. Researchers continue to concentrate on Neptune in support of the 1989 Voyager encounter. Occultations of stars by Neptune are observed from the Tucson area and from Chile to obtain information about Neptune's atmosphere and to continue to search for Neptune arcs. Occultations by other solar system objects are also observed as part of collaborative efforts from time to time. New results on the structure of scintillations in the central flash occultation by Neptune on 20 August 1985 were derived. Analysis shows that scintillations are present throughout the lightcurve, both near the half-intensity points (at a pressure of 1 microbar) and near the central flash (at 0.4 mbar). Near the planetary limb, the scintillations are extended parallel to the limb; near the shadow center, they are extended parallel to the limb; near the shadow center, they are extended in a radial direction. Researchers collaborated with Ramesh Narayan to derive a theory relating the scintillations to density fluctuations in Neptune's atmosphere. The theory will ultimately enable researchers to test whether the scintillations are caused by internal gravity waves in Neptune's upper atmosphere.

Hubbard, W. B.↗
DOE OSTI2025

LATTE: open-source, high-performance traveltime computation, tomography and source location in acoustic and elastic media

Traveltime-based tomography and source location are fundamental approaches for imaging subsurface structures and understanding the spatiotemporal distribution of seismicity from local to global scales. We present an open-source, high-performance framework integrating eikonal equation solvers and adjoint-state theory for traveltime computation, velocity tomography, source location and joint tomography-location in 2-D/3-D acoustic and elastic media. We introduce novel regularization schemes based on total generalized p-variation, structural similarity and multitask machine learning to enhance the fidelity and interpretability of inverted models and source locations. Key features of our implementation also include the ability to leverage both absolute-difference and double-difference traveltime misfits for high-fidelity velocity tomography and source parameter estimation; support for traveltime computation and inversion in diverse 2-D/3-D scenarios with arbitrary source and receiver distributions; and a perturbation-based optimal step-size estimation method to reduce computational costs. In addition, our implementation employs shared-memory and distributed-memory parallelization to provide an efficient solution for traveltime computation, tomography, and source location. In conclusion, we validate the efficacy and accuracy of our approach through multiple synthetic data examples.

58 GEOSCIENCES↗
DOE OSTI2017

Experimental study of shock-accelerated inclined heavy gas cylinder

An experimental study examines shock acceleration with an initially diffuse cylindrical column of sulfur hexafluoride surrounded by air and inclined with respect to the shock front. Three-dimensional vorticity deposition produces flow patterns whose evolution is captured with planar laser-induced fluorescence in two planes. Both planes are thus parallel to the direction of the shock propagation. The first plane is vertical and passes through the axis of the column. The second visualization plane is normal to the first plane and passes through the centerline of the shock tube. Vortex formation in the vertical and centerline planes is initially characterized by different rates and morphologies due to differences in initial vorticity deposition. In the vertical plane, the vortex structure manifests a periodicity that varies with Mach number. The dominant wavelength in the vertical plane can be related to the geometry and compressibility of the initial conditions. At later times, the vortex interaction produces a complex and irregular three-dimensional pattern suggesting transition to turbulence. We present highly repeatable experimental data for Mach numbers 1.13, 1.4, 1.7, and 2.0 at column incline angles of 0, 20, and 30 degrees for about 50 nominal cylinder diameters (30 cm) of downstream travel.

42 ENGINEERING↗
NASA NTRS2005

Recent Advances in Radar Polarimetry and Polarimetric SAR Interferometry

The development of Radar Polarimetry and Radar Interferometry is advancing rapidly, and these novel radar technologies are revamping Synthetic Aperture Radar Imaging decisively. In this exposition the successive advancements are sketched; beginning with the fundamental formulations and high-lighting the salient points of these diverse remote sensing techniques. Whereas with radar polarimetry the textural fine-structure, target-orientation and shape, symmetries and material constituents can be recovered with considerable improvements above that of standard amplitude-only Polarization Radar ; with radar interferometry the spatial (in depth) structure can be explored. In Polarimetric-Interferometric Synthetic Aperture Radar (POL-IN-SAR) Imaging it is possible to recover such co-registered textural plus spatial properties simultaneously. This includes the extraction of Digital Elevation Maps (DEM) from either fully Polarimetric (scattering matrix) or Interferometric (dual antenna) SAR image data takes with the additional benefit of obtaining co-registered three-dimensional POL-IN-DEM information. Extra-Wide-Band POL-IN-SAR Imaging - when applied to Repeat-Pass Image Overlay Interferometry - provides differential background validation and measurement, stress assessment, and environmental stress-change monitoring capabilities with hitherto unattained accuracy, which are essential tools for improved global biomass estimation. More recently, by applying multiple parallel repeat-pass EWB-POL-D(RP)-IN-SAR imaging along stacked (altitudinal) or displaced (horizontal) flight-lines will result in Tomographic (Multi- Interferometric) Polarimetric SAR Stereo-Imaging , including foliage and ground penetrating capabilities. It is shown that the accelerated advancement of these modern EWB-POL-D(RP)-IN-SAR imaging techniques is of direct relevance and of paramount priority to wide-area dynamic homeland security surveillance and local-to-global environmental ground-truth measurement and validation, stress assessment, and stress-change monitoring of the terrestrial and planetary covers. In addition, various closely related topics of (i) acquiring additional and protecting existing spectral windows of the Natural Electromagnetic Spectrum (NES) pertinent to Remote Sensing; (ii) mitigating against common "Radio Frequency Interference (RFI)" and intentional Directive Jamming of Airborne & Space borne POL-IN-SAR Imaging Platforms are appraised.

Boerner, Wolfgang-Martin↗