SEARCH · Search NASA

Results for “Spatial refinement”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows: Minimization of Delta * B Numerical Error - II

The generalization of a class of low-dissipative high order filter finite difference schemes for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been developed. The new scheme consists of a divergence free preserving high order spatial base scheme with a filter approach which can be divergence-free preserving depending on the type of filter operator being used, the method of applying the filter step, and the type of flow problem to be considered. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta * B) numerical error in the sense that no standard divergence cleaning is required. Performance evaluation of these variants, and the key role that the proper treatment of their corresponding numerical boundary conditions can play will be illustrated. Many levels of grid refinement and detailed comparison with several commonly used compressible MHD shock-capturing schemes will be sought. For certain MHD 2-D test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Sjoegreen, Bjoern↗

NASA NTRS1991

Dynamic edge warping - An experimental system for recovering disparity maps in weakly constrained systems

Dynamic edge warping (DEW), a technique for recovering reasonably accurate disparity maps from uncalibrated stereo image pairs, is presented. No precise knowledge of the epipolar camera geometry is assumed. The technique is embedded in a system including structural stereopsis on the front end and robust estimation in digital photogrammetry on the other for the purpose of self-calibrating stereo image pairs. Once the relative camera orientation is known, the epipolar geometry is computed and the system can use this information to refine its representation of the object space. Such a system will find application in the autonomous extraction of terrain maps from stereo aerial photographs, for which camera position and orientation are unknown a priori, and for online autonomous calibration maintenance for robotic vision applications, in which the cameras are subject to vibration and other physical disturbances after calibration. This work thus forms a component of an intelligent system that begins with a pair of images and, having only vague knowledge of the conditions under which they were acquired, produces an accurate, dense, relative depth map. The resulting disparity map can also be used directly in some high-level applications involving qualitative scene analysis, spatial reasoning, and perceptual organization of the object space. The system as a whole substitutes high-level information and constraints for precise geometric knowledge in driving and constraining the early correspondence process.

Boyer, K. L.↗

NASA NTRS1994

Effects of spatial order of accuracy on the computation of vortical flowfields

The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Ekaterinaris, J. A.↗

NASA NTRS2021

Detection and Tracking of Aircraft in the Far-Field from Small Unmanned Aerial Systems

Onboard far-field aircraft detection is needed for safe non-cooperative traffic mitigation in autonomous small Unmanned Aerial System (sUAS) operations. Machine vision systems, based on standard optics and visible light detectors, possess the ideal size, weight, and power (SWaP) requirements for sUAS. This work presents the design and analysis of a novel aircraft detection and tracking pipeline based on optical sensing alone. Key contributions of the work include a refined range inequality model based on sensing and detection with FAA well-clear separation assurance distances between aircraft in mind, a detector fusion method to maximize the benefit of two image detectors, and a comparative analysis of Linear Kalman-filtering and Extended Kalman-filtering to seek optimal tracking performance. The pipeline is evaluated offline against multiple intruder platforms, using two types of flight encounters: multirotor sUAS vs. fixed-wing sUAS and multirotor sUAS vs. general aviation(GA)plane. Analysis is restricted to the rate-limiting head-on and departing collision volume cases vertically separated for safety. Results indicate that it is feasible to use the proposed optical spatial-temporal tracking algorithm to provide adequate alerting time to prevent penetration of well-clear separation volumes for both sUAS and GA aircraft.

Unmanned Aerial System↗

NASA NTRSDate not supplied

Isolator Shock Train Response to Dual-Mode Scramjet Throttling-Like Forcing

Typical dual-mode scramjet designs employ an isolator component situated between the inlet and combustor. The isolator fluidically buffers the inlet from the combustor via the formation of a shock train, an interconnected series of shock-boundary layer interactions. The isolator is required during ramjet-engine operation and at the lower edges of the scramjet operational envelope. In an operational engine, transient behavior imposes fluidic dynamics due to attitude adjustments, combustor throttling, aberrant resonances within the combustor, etc. To accommodate these transients, the shock train adjusts on both temporal and spatial scales. Sufficiently large downstream pressure increases in the combustor play a critical role during unstart, whereby strong compression waves propagate upstream in the isolator to ultimately force the shock train out the inlet, if it cannot stabilize to the pressure increase. This critically reduces the engine mass flow rate, potentially causing a catastrophic reduction in thrust and loss of the vehicle. A nearly instantaneous unstart can, in principle, be predicted by a ratio of the normal shock pressure jump in an inviscid isolator. However, lower-magnitude perturbations, below this bound, also induce unstarts. Understanding the transient behavior inherent in unsteady isolator operation is critical to reducing safety margin in isolator length and refining engine control algorithms to improve engine performance. The Isolator Dynamics Research Lab (IDRL) is a cold-flow, direct-connect isolator facility ideal for assessing the fundamental fluid dynamics of shock trains. The internal flow of a scramjet isolator is modeled using a 2D, symmetric Mach 2.5 nozzle followed by a 695 mm long, 25.4x50.8 mm rectangular isolator test section. A blockage cone on a high-speed ball screw actuator downstream of the isolator is pulsed to simulate the throttling behavior in an operational combustor. Wall flush high-speed static pressure transducers detect the streamwise shock train within the isolator. Time delays between the onset of the throttling and the response motion of the shock train are analyzed to assess the controlling scales and parameters. From this analysis, the response time of the isolator shock train system is found to depend primarily on the isolator back pressure time-rate-of-change but quickly saturates to a maximum propagation speed along the isolator. These results inform dual-mode scramjet engine designers of the critical parameters and limits of an unstart resistant engine.

scramjet↗

DOE OSTI2025

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Localizing tetrahedral aluminum in nitrate-bearing gibbsite to constrain defect-impurity coupling

The enhanced radiolytic stability of gibbsite (α-Al(OH) 3 ) containing trace nitrate (NO 3 − ) is a phenomenon in nuclear waste management, but its structural origins remain unresolved. Motivated by the detection of minority tetrahedral aluminum (T d ) defects in synthetic gibbsite, we hypothesized that these sites may participate in NO 3 − retention or mediate H 2 suppression. To evaluate this, we combined orthogonal techniques comprised of spatially selective solid-state 27 Al MAS NMR, comparative spectroscopy, and density functional theory (DFT) modeling. Paramagnetic editing and dynamic nuclear polarization (DNP) MAS NMR confirm that T d defects are confined to the particle interior. DFT calculations reveal no energetic stabilization of NO 3 − near T d sites. Comparative NMR analysis shows that T d is also present in chloride-bearing gibbsite, which exhibits high radiolytic hydrogen yields. These three independent disqualifications rule out T d as a structural contributor to nitrate-mediated suppression and narrow the scope of defect-driven explanations. The findings redirect mechanistic attention away from coordination defects and toward redox-active impurity pathways, providing a refined foundation for understanding radiation tolerance in Al(OH) 3 .

Graham, Trent R. [Pacific Northwest National Labor↗

NASA NTRS1997

On Approximate Factorization Schemes for Solving the Full Potential Equation

An approximate factorization scheme based on the AF2 algorithm is presented for solving the three-dimensional full potential equation for the transonic flow about isolated wings. Two spatial discretization variations are presented, one using a hybrid first-order/second-order-accurate scheme and the second using a fully second-order-accurate scheme. The present algorithm utilizes a C-H grid topology to map the flow field about the wing. One version of the AF2 iteration scheme is used on the upper wing surface and another slightly modified version is used on the lower surface. These two algorithm variations are then connected at the wing leading edge using a local iteration technique. The resulting scheme has improved linear stability characteristics and improved time-like damping characteristics relative to previous implementations of the AF2 algorithm. The presentation is highlighted with a grid refinement study and a number of numerical results.

Holst, Terry L.↗

NASA NTRS2011

Land Surface Verification Toolkit (LVT) - A Generalized Framework for Land Surface Model Evaluation

Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it also supports hydrological data products from other, non-LIS environments. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

Kumar, Sujay V.↗

NASA NTRS2020

Entropy Stable h/p-Nonconforming Discretization with the Summation-by-Parts Property for the Compressible Euler and Navier–Stokes Equations

In this paper, we extend the entropy conservative/stable algorithms presented by Del Rey Fernandez and coauthors for the compressible Euler and Navier-Stokes equations on nonconforming p-refined/coarsened curvilinear grids to h/p refinement/coarsening. The main difficulty in developing nonconforming algorithms is the construction of appropriate coupling procedures across nonconforming interfaces. Here, we utilize a computationally simple and efficient approach based upon using decoupled interpolation operators. The resulting scheme is entropy conservative/stable and elementwise conservative. Numerical simulations of the isentropic vortex and viscous shock propagation con firm the entropy conservation/stability and accuracy properties of the method (achieving ~ p + 1 convergence), which are comparable to those of the original conforming scheme. Simulations of the Taylor{Green vortex at R(e) = 1,600 and turbulent flow past a sphere at R(e(infinity)) = 2,000 show the robustness and stability properties of the overall spatial discretization for unstructured grids. Finally, to demonstrate the entropy conservation property of a fully-discrete explicit entropy stable algorithm with h=p refinement/coarsening, we present the time evolution of the entropy function obtained by simulating the propagation of the isentropic vortex using a relaxation Runge-Kutta scheme.

Nonconforming interfaces↗

NASA NTRSDate not supplied

Development of a Consistent MODIS and VIIRS Cloud Detection Approach for CERES

A consistent cloud fraction record across various satellite platforms is essential for maintaining a long-term and stable climate data record of Earth's energy budget. With the Aqua satellite nearing the end of its operational lifetime, the continuation of this record relies on utilizing VIIRS observations from NOAA20 for cloud detection in NASA’s Clouds and Earth’s Radiant Energy System (CERES) project. However, integrating data from VIIRS and MODIS instruments poses challenges due to their distinct characteristics, such as varying spatial resolutions and different spectral channels. As a result, deriving consistent cloud properties from these two sensors without introducing artificial discontinuities in the time series remains a complex and challenging task. This paper will present progress toward developing a unified MODIS and VIIRS cloud mask using common channels to produce consistent cloud properties for CERES next edition (Ed5) Earth radiation budget data products. The fundamental approach taken in the CERES cloud mask is to compare the observed radiances to the expected background clear sky radiances. Therefore, one vital step is to compute clear sky radiances with a radiative transfer model that accurately accounts for satellite-specific, spectrally dependent surface reflectance, surface emission, and atmospheric absorption. Refined radiative transfer models and updated ancillary data inputs including surface emissivity maps, IGBP, snow and ice maps are incorporated into the processing framework to improve cloud detection consistency and accuracy. Pixel level cloud mask results and monthly global cloud fraction comparisons between MODIS and VIIRS will be presented to evaluate their consistency. Remaining challenges will be discussed. It is expected that this work will contribute consistent cloud properties for CERES that adequately bridges the MODIS and VIIRS imager data records.

CERES↗

DOE OSTI2026

Performance and Stability Characterization of the ICARUS Light Detection System

The ICARUS detector, a key component of the Short Baseline Neutrino (SBN) Program at Fermilab, consists of two identical T300 modules filled with liquid argon. It is equipped with a Light Detection System (LDS) based on 360 8-inch Hamamatsu R5912-MOD photomultiplier tubes (PMTs) arranged behind the wire planes to collect Vacuum Ultraviolet (VUV, $\sim$ 128 nm) scintillation light. Operating under cryogenic conditions ( $\sim$ 87 K), the LDS is essential for determining the event start time (t0) with nanosecond precision for beam spill synchronization, improving longitudinal spatial resolution, and contributing to the event trigger and cosmic-ray mitigation. The performance of the LDS was investigated addressing both hardware and data analysis aspects. Following a progressive degradation in PMT gain observed during operations at FNAL, systematic gain measurements were first carried out from room temperature down to low temperatures. The results show stable performance at room temperature but a significant, irreversible reduction in gain at low temperatures. Based on these findings, a series of mitigation strategies were implemented in the ICARUS detector to preserve PMT performance and ensure reliable cryogenic operation. Currently, ongoing waveform analysis of the PMT signals is being performed to characterize signal shape, charge integration, and timing properties, aiming to refine and improve the agreement between experimental data and Monte Carlo simulations.

Saia, Clara [U. Catania (main); INAF, Catania; Cat↗

DOE OSTI2025

Uncertainties in Simulating Flooding During Hurricane Harvey Using 2D Shallow Water Equations

Abstract Flooding is one of the most impactful weather‐related natural hazards. Numerical models that solve the two dimensional (2D) shallow water equations (SWE) represent the first‐principles approach to simulate all types of spatial flooding, such as pluvial, fluvial, and coastal flooding, and their compound dynamics. High spatial resolution (e.g., () m) is needed in 2D SWE simulations to capture flood dynamics accurately, resulting in formidable computational challenges. Thus, relatively coarser spatial resolutions are used for large‐scale simulations of flooding, which introduce uncertainties in the results. It is unclear how the uncertainty associated with the model resolution compares to the uncertainties in precipitation data sets and assumptions regarding boundary conditions when channelized flows interact with other water bodies. In this study, we compare these three sources of uncertainties in 2D SWE simulations for the 2017 Houston flooding event. Our results show that precipitation uncertainty and mesh resolution have more significant impacts on the simulated streamflow and inundation dynamics than the choice of the downstream boundary condition at the watershed outlet. We point out the viability to confine the uncertainty of coarsening mesh resolution by using the variable resolution mesh (VRM) which refines critical topographic features with far fewer grid cells. Specifically, in simulations with VRM, the simulated inundation depths over the refined region are comparable to that use the finest uniform mesh. This study contributes to understanding the challenges and pathways for applying 2D SWE models to improve the realism of flood simulations over large scales.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

DOE OSTI2021

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗

DOE OSTI2025

Steady-state mixing state of black carbon aerosols from a particle-resolved model

Abstract. Black carbon (BC) exerts a notable warming effect due to its strong light absorption, largely influenced by its “mixing state”. However, due to computational constraints, the mixing state is challenging to accurately represent in large-scale models. In this study, we employ a particle-resolved model to simulate the evolution of BC mixing state based on field observation. Our result shows that aerosol compositions, coating thickness (CT) distribution, and optical properties of BC aerosols all exhibit a tendency toward a steady state with a characteristic timescale of less than 1 d, considerably shorter than the BC atmospheric lifetime. The rapid attainment of a steady state suggests that it is reasonable to disregard this pre-steady-state period and instead concentrate on the average properties of BC across extensive spatial and temporal scales. The distribution of CT follows an exponential linear distribution and can be characterized by a single slope parameter k. This distribution is independent of the BC core's distribution. In the model simulation, the mean CT, equivalent to the 1/k, is 62 nm, which is consistent with the statistical results indicating a mean CT of 63 nm. Utilizing the slope parameter k, which effectively characterizes the CT distribution under the steady-state simplifying assumption, the BC absorption enhancement closely corresponds to the results obtained via the particle-resolved method. This study simplifies the BC mixing state description and yields a precise evaluation of the BC optical properties, which has the potential utility for modeling efforts in the refinement of the assessment of BC's radiative effects.

Zhang, Zhouyang↗

DOE OSTI2023

Deformation mechanisms in an additively manufactured dual-phase eutectic high-entropy alloy

Nanostructured metals and alloys often exhibit high strengths but at the expense of reduced ductility. Through harnessing the far-from-equilibrium processing conditions of laser powder-bed fusion (L-PBF) additive manufacturing, we develop a dual-phase nanolamellar structure comprised of FCC/L1 2 and BCC/B2 phases in a Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 eutectic high-entropy alloy (EHEA), which exhibits a combination of ultrahigh yield strength (>1.4 GPa) and large tensile ductility (~17%). The deformation mechanisms of the additively manufactured EHEA are studied via in-situ synchrotron X-ray diffraction and high-resolution transmission electron microscopy. The high yield strength mainly results from effective blockage of dislocation motion by the high density of lamellar interfaces. The refined nanolamellar structures and low stacking fault energy (SFE) promote stacking fault (SF)-mediated deformation in FCC/L1 2 nanolamellae. The accumulation of abundant dislocations and SFs at lamellar interfaces can effectively elevate local stresses to promote dislocation multiplication and martensitic transformation in BCC/B2 nanolamellae. The cooperative deformation of the dual phases, assisted by the semi-coherent lamellar interfaces, gives rise to the large ductility of the as-printed EHEA. In addition, here we also demonstrate that post-printing heat treatment allows us to tune the non-equilibrium microstructures and deformation mechanisms. After annealing, the significantly reduced SFE and thicknesses of the FCC nanolamellae facilitate the formation of massive SFs. The dissolution of nano-precipitates in the BCC/B2 nanolamellae reduces spatial confinement and further promotes martensitic transformation to enhance work hardening. Our work provides fundamental insights into the rich variety of deformation mechanisms underlying the exceptional mechanical properties of the additively manufactured dual-phase nanolamellar EHEAs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Local Resource Assessment: Beetle-Killed Spruce in Alaska

This report synthesizes forest inventory and spatial data from federal and state sources to present a snapshot of beetle-killed spruce across Southcentral and Interior Alaska. This natural resource is under investigation for use in cellulose-based products, particularly building insulation. Research funded by the Department of Energy's ARPA-E program, has successfully developed insulation - Celium - using beetle-killed spruce combined with mycelium, the root network of fungi. The resulting insulation product can be manufactured locally and used across all Alaskan climate zones for residential, commercial, and shipping applications. Removing standing dead wood also contributes to wildfire risk reduction and mitigates associated public health impacts. This Local Resource Assessment, funded by the Denali Commission, contextualizes and reconciles unaligned forestry data sources, establishing a pathway to refined, replicable analysis and, ultimately, presenting estimates of beetle-killed spruce biomass to inform the potential scale and development of a Celium industry.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗