Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

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 1,261 records · Page 70

Cluster Analysis of Spectroscopic Line Profiles and EUV Emission in RMHD Simulations and Observations of the Solar Atmosphere

Spatially-resolved observations from the IRIS, SDO/AIA, and other space mission and ground-based telescopes, coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar atmosphere. To better understand the dynamical and thermodynamic properties in the simulation data and their connection to observations, it is essential to determine similarities in the behaviors of the synthesized and observed emission. However, the complexity of observational data and physical processes makes comparison of observations and modeling results difficult. In this work, we show the initial results of application of K-Means clustering (unsupervised machine learning) algorithm to two different problems: 1) recognition of the typical spectroscopic line profiles observed by IRIS during solar flares and their typical dynamic behavior; 2) recognition of shocks and heating events in synthetic AIA emission data obtained from StellarBox quiet-Sun simulations. The average silhouette width technique for the KMeans algorithm is utilized in different ways to obtain optimal numbers of clusters. We discuss application of the emission clustering to visualizations of the computational volume, understanding its evolutionary trends and behavior patterns, and inversion (reconstruction) of physical properties of the solar atmosphere from synthesizes emission data.

Sadykov, Viacheslav↗

Bayesian inference in band excitation scanning probe microscopy for optimal dynamic model selection in imaging

The universal tendency in scanning probe microscopy (SPM) over the last two decades is to transition from simple 2D imaging to complex detection and spectroscopic imaging modes. The emergence of complex SPM engines brings forth the challenge of reliable data interpretation, i.e., conversion from detected signals to descriptors specific to tip–surface interactions and subsequently to material’s properties. In this work, we implemented a Bayesian inference approach for the analysis of the image formation mechanisms in band excitation SPM. Compared to the point estimates in classical functional fit approaches, Bayesian inference allows for the incorporation of extant knowledge of materials and probe behavior in the form of corresponding prior distribution and return the information on the material functionality in the form of readily interpretable posterior distributions. We explore the nonlinear mechanical behaviors spatially in a classical ferroelectric material, PbTiO 3 . We observe the non-trivial evolution of the Duffing stiffness term and the nonlinearity of the sample surface, determine spatial clustering of the nonlinear response, and perform a Landau analysis on predicting the nonlinear coefficient, which indicates that ferroelectric behavior can be a cause of the observed results. These observations suggest that the spectrum of anomalous behaviors at the ferroelectric domain walls may be broader than previously believed and can extend to non-conventional mechanical properties in addition to static and microwave conductance.

36 MATERIALS SCIENCE↗

New Software for Predicting Charging of Spacecraft

Charging System Analyzer Program (Nascap-2K) is a comprehensive update, revision, and extension of several NASA and Air Force codes for predicting electrical charging of spacecraft. Nascap-2K integrates the capabilities and models included in four independent programs: NASCAP/LEO for low-Earth orbits, NASCAP/GEO for geosynchronous orbits, POLAR for auroral charging in polar orbits, and DynaPAC (Dynamic Plasma Analysis Code) for time-dependent plasma interactions. While each of the earlier codes works well for the range of problems for which it was designed, by today s standards these codes are difficult to learn, cumbersome to use, and overly restrictive in their geometric modeling capabilities. Nascap-2K incorporates these models into a single software package that includes spacecraft surface modeling, spatial gridding, environmental specifications, calculating scripting, and post-processing analysis and visualization. The provided material properties database includes values from earlier programs as well as values from recent measurements. Development of Nascap-2K continues with future capabilities to include interactions with dense plasma such as those produced by electric propulsion.

Mandell, Myron↗

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗

Regional tropical rainfall shifts under global warming: an energetic perspective

Abstract Future climate simulations feature pronounced spatial shifts in the structure of tropical rainfall. We apply a novel atmospheric energy flux analysis to diagnose late 21st century tropical rainfall shifts in a large ensemble of simulations of 21st century climate. The method reconstructs 2D spatial changes in rainfall based on horizontal shifts in the lines of zero meridional and zonal divergent energy flux, called the energy flux equator (EFE) and energy flux prime meridian (EFPM), respectively. Two main sources of future atmospheric energy flux changes, and hence rainfall shifts, are identified by the analysis: the high-latitude North Atlantic due to a weakened Atlantic Meridional Overturning Circulation that shifts tropical rainfall southwards over the greater Tropical Atlantic sector and eastern Pacific; and the eastern tropical Pacific due to a permanent El-Niño-like response that produces zonal shifts over the Maritime Continent and South America. To first order, the shifts in the EFE and EFPM mirror gross distributional changes in tropical precipitation, with a southward shift in rainfall over the tropical Atlantic, West Africa, and eastern tropical Pacific and an eastward shift over the Maritime Continent and western Pacific. When used to reconstruct future rainfall shifts in the tropical Atlantic and Sahel, the method reasonably represents the simulated meridional structure of rainfall shifts but does not do so for the zonal structures.

54 ENVIRONMENTAL SCIENCES↗

X-ray computed tomography-based porosity analysis: Algorithms and application for porous woody biomass

The properties of lignocellulosic biomass, such as bulk physical, thermal, and mechanical properties, as well as the mobility of enzymes or catalysts, are largely affected by porosity, pore structures, and pore size distribution. While X-ray computed tomography (CT) has been introduced to effectively produce 3D volumetric images of material microstructure, a quantitative porosity analysis from the 3D images has never been done. This work introduces a first-of-its-kind X-ray CT-based quantitative porosity analysis method and applies the developed toolkit to characterize the internal porosity distribution of loblolly pine. A sample loblolly pine chip was scanned by a nano-resolution X-ray CT system and digitally reconstructed after a sequence of image processing operations. Here, a comprehensive porosity analysis was then performed to quantify the envelope porosity, the spatial distribution of local porosity, and the directional porosity. Solutions to the various challenges in image processing, porosity calculation, and large data handling are provided.

42 ENGINEERING↗

Manifold projection image segmentation for nano-XANES imaging

As spectral imaging techniques are becoming more prominent in science, advanced image segmentation algorithms are required to identify appropriate domains in these images. We present a version of image segmentation called manifold projection image segmentation (MPIS) that is generally applicable to a broad range of systems without the need for training because MPIS uses unsupervised machine learning with a few physically motivated hyperparameters. We apply MPIS to nanoscale x-ray absorption near edge structure (XANES) imaging, where XANES spectra are collected with nanometer spatial resolution. We show the superiority of manifold projection over linear transformations, such as the commonly used principal component analysis (PCA). Moreover, MPIS maintains accuracy while reducing computation time and sensitivity to noise compared to the standard nano-XANES imaging analysis procedure. Finally, we demonstrate how multimodal information, such as x-ray fluorescence data and spatial location of pixels, can be incorporated into the MPIS framework. We propose that MPIS is adaptable for any spectral imaging technique, including scanning transmission x-ray microscopy, where the length scale of domains is larger than the resolution of the experiment.

36 MATERIALS SCIENCE↗

Multi-scale variation in spatial heterogeneity for microbial community structure in an eastern Virginia agricultural field

To better understand the distribution of soil microbial communities at multiple spatial scales, a survey was conducted to examine the spatial organization of community structure in a wheat field in eastern Virginia (USA). Nearly 200 soil samples were collected at a variety of separation distances ranging from 2.5 cm to 11 m. Whole-community DNA was extracted from each sample, and community structure was compared using amplified fragment length polymorphism (AFLP) DNA fingerprinting. Relative similarity was calculated between each pair of samples and compared using geostatistical variogram analysis to study autocorrelation as a function of separation distance. Spatial autocorrelation was found at scales ranging from 30 cm to more than 6 m, depending on the sampling extent considered. In some locations, up to four different correlation length scales were detected. The presence of nested scales of variability suggests that the environmental factors regulating the development of the communities in this soil may operate at different scales. Kriging was used to generate maps of the spatial organization of communities across the plot, and the results demonstrated that bacterial distributions can be highly structured, even within a habitat that appears relatively homogeneous at the plot and field scale. Different subsets of the microbial community were distributed differently across the plot, and this is thought to be due to the variable response of individual populations to spatial heterogeneity associated with soil properties. c2003 Federation of European Microbiological Societies. Published by Elsevier Science B.V. All rights reserved.

NASA Discipline Life Support Systems↗

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Doppler lidar data

Dual Doppler lidar analyses of data taken by pulsed lidars demonstrated feasibility of deriving wind fields from coordinated lidar scans. Limited case histories of thunderstorm outflows were obtained. Co-located comparison between Marshall Space Flight Center lidar and NCAR 5.5 cm radar demonstrated desirability of lidar in cases of marginal radar reflectivity in clear air and low-elevation scans. Analysis continued on backscattered intensity and velocity measurements made from April 1983 to February 1984. A slant path method was used to calculate vertical profiles of volumetric backscatter and adsorption in the lower troposphere. High-quality VAD scans were identified as candidates for investigating feasibility of calculating horizontal motion fields using single Doppler lidar. Activities during FY-85 also included participation in Fall 1984 airborne Doppler lidar flight experiments. Preliminary data review was begun using McIdas system. Analysis of backscatter and absorpiton profiles continues. Focus is on understanding spatial and temporal variations, as well as frequency distribution, of backscatter at several tropospheric levels. Results from this study provide input to evaluation of clean/dirty airmass hypothesis of aerosol distribution. Assistance is being given to preparation of a comprehensive, global backscatter measurement plan. Analysis of data from Fall 1984 flight experiments is just beginning. Work has begun on preprocessing data to minimize errors due to electro-optic modulator malfunction during flights.

Rothermel, J.↗

Potential Multi-Component Structure of the Debris Disk Around HIP 17439 Revealed by Herschel DUNES

Context. The dust observed in debris disks is produced through collisions of larger bodies left over from the planet/planetesimal formation process. Spatially resolving these disks permits to constrain their architecture and thus that of the underlying planetary/planetesimal system. Aims. Our Herschel open time key program DUNES aims at detecting and characterizing debris disks around nearby, sun-like stars. In addition to the statistical analysis of the data, the detailed study of single objects through spatially resolving the disk and detailed modeling of the data is a main goal of the project. Methods. We obtained the first observations spatially resolving the debris disk around the sun-like star HIP 17439 (HD 23484) using the instruments PACS and SPIRE on board the Herschel Space Observatory. Simultaneous multi-wavelength modeling of these data together with ancillary data from the literature is presented. Results. A standard single component disk model fails to reproduce the major axis radial profiles at 70 μm, 100 μm, and 160 μm simultaneously. Moreover, the best-fit parameters derived from such a model suggest a very broad disk extending from few au up to few hundreds of au from the star with a nearly constant surface density which seems physically unlikely. However, the constraints from both the data and our limited theoretical investigation are not strong enough to completely rule out this model. An alternative, more plausible, and better fitting model of the system consists of two rings of dust at approx. 30 au and 90 au, respectively, while the constraints on the parameters of this model are weak due to its complexity and intrinsic degeneracies. Conclusions. The disk is probably composed of at least two components with different spatial locations (but not necessarily detached), while a single, broad disk is possible, but less likely. The two spatially well-separated rings of dust in our best-fit model suggest the presence of at least one high mass planet or several low-mass planets clearing the region between the two rings from planetesimals and dust.

sun-like stars↗

Self-consistent, high-order spatial profiles in a model for two-fluid turbulent mixing

A Reynolds-averaged Navier-Stokes model is presented with the property that it admits self-consistent, high-order spatial profiles in simulations of two-fluid turbulent mixing layers. Whereas previous models have been limited by the assumption of a linear mixing profile, the present paper relaxes this assumption and, as a result, is shown to achieve much better agreement with experimental profiles. Similarity analysis is presented to derive constraints on model coefficients to enforce desired self-similar growth rates that are fully consistent with the high-order spatial profiles. Through this similarity analysis, it is shown that care must be taken in model construction, as it is possible to construct certain terms in such a way as to leave growth rates unconstrained. This model, termed the k-Φ-L-a-V model, is then applied in simulations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz mixing layers. Overall, these simulations confirm that the expected growth parameters are recovered and high-order spatial profiles are maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Chandra Interactive Analysis of Observations (CIAO)

The Chandra (formerly AXAF) telescope, launched on July 23, 1999, provides X-rays data with unprecedented spatial and spectral resolution. As part of the Chandra scientific support, the Chandra X-ray Observatory Center provides a new data analysis system, CIAO ("Chandra Interactive Analysis of Observations"). We will present the main components of the system: "First Look" analysis; SHERPA: a multi-dimensional, multi-mission modeling and fitting application; Chandra Imaging and Plotting System; Detect package-source detection algorithms; and DM package generic data manipulation tools, We will set up a demonstration of the portable version of the system and show examples of Chandra Data Analysis.

Dobrzycki, Adam↗

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g., ELM filaments) in the scrape-off layer and at the plasma edge with 2.5 μs temporal and 10 mm optical resolution. A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics and then extended to a database of 159 ELM events. Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average, an ELM crash consists of a filament with a circular cross section, which is propelled outward with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100 μs until the radial propulsion settles back to the pre-ELM level. Furthermore, the experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current–filament interaction model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crack Detection with Lamb Wave Wavenumber Analysis

In this work, we present our study of Lamb wave crack detection using wavenumber analysis. The aim is to demonstrate the application of wavenumber analysis to 3D Lamb wave data to enable damage detection. The 3D wavefields (including vx, vy and vz components) in time-space domain contain a wealth of information regarding the propagating waves in a damaged plate. For crack detection, three wavenumber analysis techniques are used: (i) two dimensional Fourier transform (2D-FT) which can transform the time-space wavefield into frequency-wavenumber representation while losing the spatial information; (ii) short space 2D-FT which can obtain the frequency-wavenumber spectra at various spatial locations, resulting in a space-frequency-wavenumber representation; (iii) local wavenumber analysis which can provide the distribution of the effective wavenumbers at different locations. All of these concepts are demonstrated through a numerical simulation example of an aluminum plate with a crack. The 3D elastodynamic finite integration technique (EFIT) was used to obtain the 3D wavefields, of which the vz (out-of-plane) wave component is compared with the experimental measurement obtained from a scanning laser Doppler vibrometer (SLDV) for verification purposes. The experimental and simulated results are found to be in close agreement. The application of wavenumber analysis on 3D EFIT simulation data shows the effectiveness of the analysis for crack detection. Keywords: : Lamb wave, crack detection, wavenumber analysis, EFIT modeling

Tian, Zhenhua↗

Spatial distribution of earth flux density from unrestricted field of view radiation measurements

Results of an analysis to determine if unrestricted field of view measurements can be used to estimate flux densities on scales smaller than the detector's Earth field of view are reported. Simple geometric considerations were employed to reduce simulated irradiances to Earth flux densities. A spatial distribution of emitted flux density is also presented and discussed which was inferred from unrestricted field of view measurements of Nimbus 6.

Weaver, W. L.↗