Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational challenges”

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 739 records · Page 41

Spectral Ratio Imaging with Hyperion Satellite Data for Geological Mapping

Since the advent of LANDSAT I in 1972, many different multispectral satellites have been orbited by the U.S. and other countries. These satellites have varied from 4 spectral bands in LANDSAT I to 14 spectral bands in the ASTER sensor aboard the TERRA space platform. Hyperion is a relatively new hyperspectral sensor with over 220 spectral bands. The huge increase in the number of spectral bands offers a substantial challenge to computers and analysts alike when it comes to the task of mapping features on the basis of chemical composition, especially if little or no ground truth is available beforehand from the area being mapped. One approach is the theoretical approach of the modeler, where all extraneous information (atmospheric attenuation, sensor electronic gain and offset, etc.) is subtracted off and divided out, and laboratory (or field) spectra of materials are used as training sets to map features in the scene of similar composition. This approach is very difficult to keep accurate because of variations in the atmosphere, solar illumination, and sensor electronic gain and offset that are not always perfectly recorded or accounted for. For instance, to apply laboratory or field spectra of materials as data sets from the theoretical approach, the header information of the files must reflect the correct, up-to-date sensor electronic gain and offset and the analyst must pick the exact atmospheric model that is appropriate for the day of data collection in order for classification procedures to accurately match pixels in the scene with the laboratory or field spectrum of a desired target on the basis of the hyperspectral data. The modeling process is so complex that it is difficult to tell when it is operating well or determine how to fix it when it is incorrect. Recently RSI has announced that the latest version of their ENVI software package is not performing atmospheric corrections correctly with the FLAASH atmospheric model. It took a long time to determine that it was wrong, and may take an equally long time (or longer) to fix.

Vincent, Robert K.↗

Three-Dimensional Application of DAMAS Methodology for Aeroacoustic Noise Source Definition

At the 2004 AIAA/CEAS Aeroacoustic Conference, a breakthrough in acoustic microphone array technology was reported by the authors. A Deconvolution Approach for the Mapping of Acoustic Sources (DAMAS) was developed which decouples the array design and processing influence from the noise being measured, using a simple and robust algorithm. For several prior airframe noise studies, it was shown to permit an unambiguous and accurate determination of acoustic source position and strength. As a follow-on effort, this paper examines the technique for three-dimensional (3D) applications. First, the beamforming ability for arrays, of different size and design, to focus longitudinally and laterally is examined for a range of source positions and frequency. Advantage is found for larger array designs with higher density microphone distributions towards the center. After defining a 3D grid generalized with respect to the array s beamforming characteristics, DAMAS is employed in simulated and experimental noise test cases. It is found that spatial resolution is much less sharp in the longitudinal direction in front of the array compared to side-to-side lateral resolution. 3D DAMAS becomes useful for sufficiently large arrays at sufficiently high frequency. But, such can be a challenge to computational capabilities, with regard to the required expanse and number of grid points. Also, larger arrays can strain basic physical modeling assumptions that DAMAS and all traditional array methodologies use. An important experimental result is that turbulent shear layers can negatively impact attainable beamforming resolution. Still, the usefulness of 3D DAMAS is demonstrated by the measurement of landing gear noise source distributions in a difficult hard-wall wind tunnel environment.

Brooks, Thomas F.↗

Turbulent Output-Based Anisotropic Adaptation

Controlling discretization error is a remaining challenge for computational fluid dynamics simulation. Grid adaptation is applied to reduce estimated discretization error in drag or pressure integral output functions. To enable application to high O(10(exp 7)) Reynolds number turbulent flows, a hybrid approach is utilized that freezes the near-wall boundary layer grids and adapts the grid away from the no slip boundaries. The hybrid approach is not applicable to problems with under resolved initial boundary layer grids, but is a powerful technique for problems with important off-body anisotropic features. Supersonic nozzle plume, turbulent flat plate, and shock-boundary layer interaction examples are presented with comparisons to experimental measurements of pressure and velocity. Adapted grids are produced that resolve off-body features in locations that are not known a priori.

Park, Michael A.↗

Technique for Calculating Solution Derivatives With Respect to Geometry Parameters in a CFD Code

A solution has been developed to the challenges of computation of derivatives with respect to geometry, which is not straightforward because these are not typically direct inputs to the computational fluid dynamics (CFD) solver. To overcome these issues, a procedure has been devised that can be used without having access to the mesh generator, while still being applicable to all types of meshes. The basic approach is inspired by the mesh motion algorithms used to deform the interior mesh nodes in a smooth manner when the surface nodes, for example, are in a fluid structure interaction problem. The general idea is to model the mesh edges and nodes as constituting a spring-mass system. Changes to boundary node locations are propagated to interior nodes by allowing them to assume their new equilibrium positions, for instance, one where the forces on each node are in balance. The main advantage of the technique is that it is independent of the volumetric mesh generator, and can be applied to structured, unstructured, single- and multi-block meshes. It essentially reduces the problem down to defining the surface mesh node derivatives with respect to the geometry parameters of interest. For analytical geometries, this is quite straightforward. In the more general case, one would need to be able to interrogate the underlying parametric CAD (computer aided design) model and to evaluate the derivatives either analytically, or by a finite difference technique. Because the technique is based on a partial differential equation (PDE), it is applicable not only to forward mode problems (where derivatives of all the output quantities are computed with respect to a single input), but it could also be extended to the adjoint problem, either by using an analytical adjoint of the PDE or a discrete analog.

Mathur, Sanjay↗

Spin-Up and Tuning of the Global Carbon Cycle Model Inside the GISS ModelE2 GCM

Planetary carbon cycle involves multiple phenomena, acting at variety of temporal and spacial scales. The typical times range from minutes for leaf stomata physiology to centuries for passive soil carbon pools and deep ocean layers. So, finding a satisfactory equilibrium state becomes a challenging and computationally expensive task. Here we present the spin-up processes for different configurations of the GISS Carbon Cycle model from the model forced with MODIS observed Leaf Area Index (LAI) and prescribed ocean to the prognostic LAI and to the model fully coupled to the dynamic ocean and ocean biology. We investigate the time it takes the model to reach the equilibrium and discuss the ways to speed up this process. NASA Goddard Institute for Space Studies General Circulation Model (GISS ModelE2) is currently equipped with all major algorithms necessary for the simulation of the Global Carbon Cycle. The terrestrial part is presented by Ent Terrestrial Biosphere Model (Ent TBM), which includes leaf biophysics, prognostic phenology and soil biogeochemistry module (based on Carnegie-Ames-Stanford model). The ocean part is based on the NASA Ocean Biogeochemistry Model (NOBM). The transport of atmospheric CO2 is performed by the atmospheric part of ModelE2, which employs quadratic upstream algorithm for this purpose.

Climate model↗

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗

Model Predictions and Observed Performance of JWST's Cryogenic Position Metrology System

The James Webb Space Telescope cryogenic testing requires measurement systems that both obtain a very high degree of accuracy and can function in that environment. Close-range photogrammetry was identified as meeting those criteria. Testing the capability of a close-range photogrammetric system prior to its existence is a challenging problem. Computer simulation was chosen over building a scaled mock-up to allow for increased flexibility in testing various configurations. Extensive validation work was done to ensure that the actual as-built system meet accuracy and repeatability requirements. The simulated image data predicted the uncertainty in measurement to be within specification and this prediction was borne out experimentally. Uncertainty at all levels was verified experimentally to be less than 0.1 millimeters.

Lunt, Sharon R.↗

Model Predictions and Observed Performance of JWST's Cryogenic Position Metrology System

The James Webb Space Telescope cryogenic testing requires measurement systems that both obtain a very high degree of accuracy and can function in that environment. Close-range photogrammetry was identified as meeting those criteria. Testing the capability of a close-range photogrammetric system prior to its existence is a challenging problem. Computer simulation was chosen over building a scaled mock-up to allow for increased flexibility in testing various configurations. Extensive validation work was done to ensure that the actual as-built system meet accuracy and repeatability requirements. The simulated image data predicted the uncertainty in measurement to be within specification and this prediction was borne out experimentally. Uncertainty at all levels was verified experimentally to be less than 0.1 millimeters.

Lunt, Sharon R.↗

Photoionization Calculations of the Radiation Force Due To Spectral Lines in AGNs

One of the main mechanisms that could drive mass outflows in active galactic nuclei (AGNs) is radiation pressure due to spectral lines. Although straightforward to understand, the actual magnitude of the radiation force is challenging to compute because the force depends on the physical conditions in the gas, as well as the strength, spectral energy distribution (SED), and geometry of the radiation field. We present results from our photoionization and radiation transfer calculations of the force multiplier, M(ξ, t), using the same radiation field to compute the gas photoionization and thermal balance. We assume low gas density (n = 104 per cu.cm) and column density (N ≤ 1017 per sq.cm), a Boltzmann distribution for the level populations, and the Sobolev approximation. Here, we describe results for two SEDs corresponding to an unobscured and obscured AGN in NGC 5548. Our main results are the following: (1) although M(ξ, t) starts to decrease with ξ for ξ gsim 1 as shown by others, this decrease in our calculations is relatively gradual and could be nonmonotonic as M(ξ, t) can increase by a factor of a few for ξ ≈ 10–1000; (2) at these same ξ for which the multiplier is higher than in previous calculations, the gas is thermally unstable by the isobaric criterion; (3) non-local thermodynamic equilibrium effects reduce M(t, ξ) by over two orders of magnitude for ξ gsim 100. The dynamical consequence of result (1) is that line driving can be important for ξ as high as 1000 when the LTE approximation holds, while result (2) provides a natural cloud formation mechanism that may account for the existence of narrow line regions. Result (3) suggests that line driving may not be important for ξ gsim 100 in tenuous plasma.

Dannen, Randall C.↗

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗

Overview of Mars Sample Return – Earth Entry System Woven Roughness Heating Augmentation Test in NASA Langley’s Mach 6 Wind Tunnel

The Mars Sample Return Mission (MSR) is a planned NASA flagship mission in which a sample retrieval lander (SRL) with a rover will be flown to Mars to obtain sample tubes on the surface that were dropped by the Mars 2020 rover [1]. After obtaining the sam-ples, the rover will return and ascend back to Martian orbit onboard the Mars Ascent Vehicle (MAV). Upon return to Earth orbit, the samples will perform Entry, Descent, and Landing (EDL) with the Earth Entry Sys-tem (EES) architecture, and land in Utah. The EES vehicle will utilize a HEEET-variant as its TPS, which will be the first time a woven TPS will be used on a flagship NASA mission [2]. This TPS offers a unique challenge for Computational Fluid Dynamics (CFD) modeling of the aerothermal envi-ronment of the vehicle, as woven roughness heating augmentation has not been extensively investigated experimentally. As a result, in order to validate com-putational models for woven roughness heating aug-mentation, a wind tunnel test campaign at NASA Langley Research Center’s Mach 6 wind tunnel was performed in April of 2023. This test campaign consisted of over a hundred runs with Reynolds numbers spanning from 1-7 mil-lion 1/ft and with six separate wind tunnel models used. A second campaign with a suite of new models will be conducted in Summer 2023 as well as a cam-paign with a flat plate model, both of which are of great interest to the MSR-EES project. The data obtained from this test are extremely vital for the MSR mission, as they will validate CFD roughness heating models which will be directly used to design the TPS of the EES portion of MSR and characterize the heating environment that the entry ve-hicle will experience. Further extensions of the MSR-EES test campaign will continue to provide validation data for developing more effective computational tools.

Jonathan Cheatwood↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

3D Cloud Masking Across A Broad Swath Using Multi-Angle Polarimetry and Deep Learning

Understanding the 3-dimensional structure of clouds is of crucial importance to modeling our changing climate. Active sensors, such as radar and lidar, provide accurate vertical cloud profiles, but are mostly restricted to along-track sampling. Passive sensors can capture a wide swath, but struggle to see beneath cloud tops. In essence, both types of products are restricted to two dimensions: as a cross-section in the active case, and an image in the passive case. However, multi-angle sensor configurations contain implicit information about 3D structure, due to parallax and atmospheric path differences. Extracting that implicit information can be challenging, requiring computationally expensive radiative transfer techniques. Machine learning, as an alternative, may be able to capture some of the complexity of a full 3D radiative transfer solution with significantly less computational expense. In this work, we make three contributions towards understanding 3D cloud structure from multi-angle polarimetry. First, we introduce a large-scale, publicly available dataset that fuses existing cloud products into a format more amenable to machine learning. This dataset treats multi-angle polarimetry as an input, and radar-based vertical cloud profiles as an output. Second, we describe and evaluate strong baseline machine learning models based that predict these profiles from the passive imagery. Notably, these models are trained only on center-swath labels, but can predict cloud profiles over the entire passive imagery swath. Third, we leverage the information-theoretic nature of machine learning to draw conclusions about the relative utility of various sensor configurations, including spectral channels, viewing angles, and polarimetry. These findings have implications for Earth-observing missions such as NASA’s Plankton, Aerosol, Cloud-ocean Ecosystem (PACE) and Atmosphere Observing System (AOS) missions, as well as in informing future applications of computer vision to atmospheric remote sensing.

clouds↗

Transversity from First Principles in QCD

Transversity observables, such as the T-odd Sivers single-spin asymmetry measured in deep inelastic lepton scattering on polarized protons and the distributions which are measured in deeply virtual Compton scattering, provide important constraints on the fundamental quark and gluon structure of the proton. In this talk I discuss the challenge of computing these observables from first principles; i.e.; quantum chromodynamics, itself. A key step is the determination of the frame-independent light-front wavefunctions (LFWFs) of hadrons - the QCD eigensolutions which are analogs of the Schroedinger wavefunctions of atomic physics. The lensing effects of initial-state and final-state interactions, acting on LFWFs with different orbital angular momentum, lead to T-odd transversity observables such as the Sivers, Collins, and Boer-Mulders distributions. The lensing effect also leads to leading-twist phenomena which break leading-twist factorization such as the breakdown of the Lam-Tung relation in Drell-Yan reactions. A similar rescattering mechanism also leads to diffractive deep inelastic scattering, as well as nuclear shadowing and non-universal antishadowing. It is thus important to distinguish 'static' structure functions, the probability distributions computed the target hadron's light-front wavefunctions, versus 'dynamical' structure functions which include the effects of initial- and final-state rescattering. I also discuss related effects such as the J = 0 fixed pole contribution which appears in the real part of the virtual Compton amplitude. AdS/QCD, together with 'Light-Front Holography', provides a simple Lorentz-invariant color-confining approximation to QCD which is successful in accounting for light-quark meson and baryon spectroscopy as well as hadronic LFWFs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

PCM-Based Durable Write Cache for Fast Disk I/O

Flash based solid-state devices (FSSDs) have been adopted within the memory hierarchy to improve the performance of hard disk drive (HDD) based storage system. However, with the fast development of storage-class memories, new storage technologies with better performance and higher write endurance than FSSDs are emerging, e.g., phase-change memory (PCM). Understanding how to leverage these state-of-the-art storage technologies for modern computing systems is important to solve challenging data intensive computing problems. In this paper, we propose to leverage PCM for a hybrid PCM-HDD storage architecture. We identify the limitations of traditional LRU caching algorithms for PCM-based caches, and develop a novel hash-based write caching scheme called HALO to improve random write performance of hard disks. To address the limited durability of PCM devices and solve the degraded spatial locality in traditional wear-leveling techniques, we further propose novel PCM management algorithms that provide effective wear-leveling while maximizing access parallelism. We have evaluated this PCM-based hybrid storage architecture using applications with a diverse set of I/O access patterns. Our experimental results demonstrate that the HALO caching scheme leads to an average reduction of 36.8% in execution time compared to the LRU caching scheme, and that the SFC wear leveling extends the lifetime of PCM by a factor of 21.6.

99 GENERAL AND MISCELLANEOUS↗

Big Data Processing for Power Grid Event Detection

In this paper we present the application of big data processing for the development of machine learning (ML) models to detect relevant events in power grid operations. This is based on almost 20TB of phasor measurement unit data corresponding to up to two years of operation of three grid interconnections which provide power to most of the United States. A significant aspect of the work consists in having all data processing performed on a single standard GPU server, from pre-processing to ML model training and testing. We describe the data and computational infrastructure, challenges faced and methods used in data processing, main findings and results. The ML approach employed for best utilization of the big data is also discussed, including sample results.

Paes Leao, Bruno↗