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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Cross-Calibration of Earth Observing System Terra Satellite Sensors MODIS and ASTER

The Advanced Spaceborne Thermal Emissive and Reflection Radiometer (ASTER) and Moderate Resolution Imaging Spectrometer (MODIS) are two of the five sensors onboard the Earth Observing System's Terra satellite. These sensors share many similar spectral channels while having much different spatial and operational parameters. ASTER is a tasked sensor and sometimes referred to a zoom camera of the MODIS that collects a full-earth image every one to two days. It is important that these sensors have a consistent characterization and calibration for continued development and use of their data products. This work uses a variety of test sites to retrieve and validate intercalibration results. The refined calibration of Collection 6 of the Terra MODIS data set is leveraged to provide the up-to-date reference for trending and validation of ASTER. Special attention is given to spatially matching radiance measurements using prelaunch spatial response characterization of MODIS. Despite differences in spectral band properties and spatial scales, ASTER-MODIS is an ideal case for intercomparison since the sensors have nearly identical views and acquisitions times and therefore can be used as a baseline of intercalibration performance of other satellite sensor pairs.

MODIS↗

NASA NTRSDate not supplied

Folded Optical Design for High Fidelity Atmospheric Emulation With A Spatial Light Modulator

Atmospheric emulators based on spatial light modulators offer the ability to test atmospheric propagation effects on a laser communication component’s performance in the laboratory setting. To create a high-fidelity atmospheric emulator, details of the optical design are key. This paper discusses the optical design choices and refinements that enabled the creation of a system that was verified to recreate multi-layer turbulence with high fidelity up to D/r0 =50. Optical design choices that affect the fidelity discussed in this paper include the characteristics of the input laser, the spatial light modulator, the holograms, the image relay optical components, and the spatial filter. Also included in this paper is a comparison of the chosen folded optical layout to an alternative angled layout.

Atmospheric turbulence↗

NASA NTRS2016

Automated Storm Tracking and the Lightning Jump Algorithm Using GOES-R Geostationary Lightning Mapper (GLM) Proxy Data

This study develops a fully automated lightning jump system encompassing objective storm tracking, Geostationary Lightning Mapper proxy data, and the lightning jump algorithm (LJA), which are important elements in the transition of the LJA concept from a research to an operational based algorithm. Storm cluster tracking is based on a product created from the combination of a radar parameter (vertically integrated liquid, VIL), and lightning information (flash rate density). Evaluations showed that the spatial scale of tracked features or storm clusters had a large impact on the lightning jump system performance, where increasing spatial scale size resulted in decreased dynamic range of the system's performance. This framework will also serve as a means to refine the LJA itself to enhance its operational applicability. Parameters within the system are isolated and the system's performance is evaluated with adjustments to parameter sensitivity. The system's performance is evaluated using the probability of detection (POD) and false alarm ratio (FAR) statistics. Of the algorithm parameters tested, sigma-level (metric of lightning jump strength) and flash rate threshold influenced the system's performance the most. Finally, verification methodologies are investigated. It is discovered that minor changes in verification methodology can dramatically impact the evaluation of the lightning jump system.

lightning jump↗

NASA NTRSDate not supplied

Folded Optical Design for High Fidelity Atmospheric Emulation With A Spatial Light Modulator

Atmospheric emulators based on spatial light modulators offer the ability to test atmospheric propagation effects on a laser communication component’s performance in the laboratory setting. To create a high-fidelity atmospheric emulator, details of the optical design are key. This paper discusses the optical design choices and refinements that enabled the creation of a system that was verified to recreate multi-layer turbulence with high fidelity up to D/r 0 =50. Optical design choices that affect the fidelity discussed in this paper include the characteristics of the input laser, the spatial light modulator, the holograms, the image relay optical components, and the spatial filter. Also included in this paper is a comparison of the chosen folded optical layout to an alternative angled layout.

Atmospheric turbulence↗

DOE OSTI2020

High Order Implicit Residual-Based Spatial Discretization Error Estimation for S N Neutron Transport

This work demonstrates our novel residual source spatial discretization error estimator (LeR/TEAD) for a DGFEM-1 discretization and assesses it along with two contemporary estimators, Ragusa and Wang's h -refinement estimator (RW) and Duo, Azmy, and Zikatanov's explicit residual-based estimator (DAZ), on a suite of Method of Manufactured Solutions (MMS) 2D problems and three realistic problem geometries. LeR/TE-AD is attractive because it directly estimates the local error in the angular flux, as opposed to a mere indicator of the error's behavior, on the same mesh and method order as the original numerical solution, thus typically being less computationally intensive than a refinement-based method. On the MMS suite, LeR/TE-AD consistently displayed a reduced performance versus its DGFEM-0 results in terms of accuracy and precision metrics, though it was not typically grossly inaccurate. This is attributed to the irregularities in the true solution across singular characteristics limiting the local accuracy of the numerical flux solution, leading to poor derivative approximations used in the residual approximations. The error transport problem then spreads the error in the residual to nearby cells, causing a greater degree of imprecision that did not afflict DAZ or RW. In testing the estimators on realistic problem geometries, however, LeR/TE-AD fared better. In practice, the true error is much larger in non-idealized geometries like in MMS, and a superlinear true solution means that RW and DAZ are not beneficially biased for DGFEM-1 error estimation. LeR/TE-AD was typically first or second in accuracy, primarily competing with RW, but the latter usually consumed 2-4 times the computational time as LeR/TE-AD, and requires a solution with four times as many unknowns. Furthermore, RW and LeR/TE-AD can be used to compute direct estimates of the error in any quantity of interest that is based on the angular ux solution, such as the fission rate density in a fuel pin, whereas DAZ requires a heuristic extension due to its norm-based nature.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

San Juan Islands Tidal Energy Characterization

The San Juan Islands are an archipelago with multiple tidal channels that produce strong tidal currents that could be harnessed for electricity generation. At present, there is limited electrical generation on the islands, with power primarily provided by a subsea cable connection with the mainland. As part of the C-MIST program, NOAA's CO-OPS program collected ADCP data at 46 stations in the San Juan Islands and adjacent waterways, from April to August of 2017. PNNL has recently developed a 3-D tidal hydrodynamic model for tidal energy resource characterization and assessment in the Salish Sea (Yang et al. 2021). Several tidal channels in the San Juan Islands were identified as top hotspots for potential tidal energy extraction. However, due to the complex geometry and inter-connected waterways, tidal currents around San Juan Islands also exhibit strong spatial and temporal asymmetry. Therefore, it is important to analyze the NOAA ADCP data to define metrics relevant to tidal energy, and further refine PNNL's tidal hydrodynamic model to accurately simulate currents in small tidal channels and account for sharp bathymetry gradients. The outcome of this technical assistance will allow developers to identify the promising opportunities for tidal energy in the San Juan Islands.

16 TIDAL AND WAVE POWER↗

DOE OSTI2024

The Scientific Case for Concurrent Neutron and X-ray Scattering and Spectroscopy

The interrogation of materials with X-rays or neutrons to determine the structure, energetics, and dynamics of materials is fundamental to advancing materials' physical and chemical science and developing innovative material technologies. A transcending challenge in developing novel materials is that progress hinges on understanding the structure and dynamics across multiple time and length scales in complex materials that feature multiple components, interfaces, and compositions. Despite the ever-growing demands on materials’ characterization, existing approaches are almost exclusively based on isolated X-ray or neutron scattering, i.e., an approach commensurate with the more narrowly defined needs of fifty years ago. A three-day workshop sponsored by the U.S. National Science Foundation (NSF) analyzed the demand for concurrent neutron and X-ray (NeX) experiments. It was held at the Spring Hill Suites, San Jose, California, from June 2 to 4, 2022. In this workshop, 70 national and international experts ascertained the crucial need to establish NeX capabilities to advance the science of complex materials and systems in the US. Here, we illustrate the need for NeX scattering and spectroscopy experiments by showcasing examples that span areas as diverse as biomaterials, energy science, soft matter, and nanomaterials. To provide NeX capability will require new instrumentation that enables concurrent experiments. Affected areas include chemistry, soft matter, quantum materials, pure and applied chemistry, bioscience, geoscience, and applied materials. NeX benefits research outcomes due to the complementarity of the two techniques, which is essential for better model refinement. While joint refinement of data from separate neutron and X-ray experiments is critical to avoid ambiguities, especially in multiphase-multicomponent materials, concurrent experiments overcome scientific and technical barriers associated with single measurements, separated by location and, thus, time. Among all the examples, these factors introduce uncertainties in the results that complicate data analysis. [1,2] [3] While models are strongly sample-dependent, the principles of joint refinement are generally applicable to these disciplines, including the development of advanced parameterization, modeling, and analysis techniques that also consider the temporal and spatial resolutions of the two methods, leading to unambiguous data interpretation. Solutions for technical barriers must be found to realize NeX experiments, including developing robust sample environments that meet the optical requirements of neutrons and X-rays.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Evaluating Variable‐Resolution CESM Over China and Western United States for Use in Water‐Energy Nexus and Impacts Modeling

Abstract Climate models are critical tools to study earth system processes and are often further downscaled to provide refined‐resolution data sets for regional water‐energy impact analyses. In this study, the Variable‐Resolution Community Earth System Model (VR‐CESM), a new dynamical downscaling method, is employed to generate 1/8° simulations for China and the western United States over 1970–2006. Precipitation, temperature, snowpack, radiation, and wind speed are compared with two 1° global climate models, regional climate model ensemble, and reference data sets to evaluate their fidelity. Simulated precipitation skill is dependent upon the reference data sets used, a slight (<+10%) and spatially consistent wet bias in China during summer and the western United States during winter. Precipitation bias in VR‐CESM is the smallest among models and demonstrates the added value provided by resolution refinement. However, temperature biases are generally greater than precipitation due to more complicated interactions among processes. VR‐CESM simulates +2°C warm bias in the low‐elevation regions in China but −2°C cold bias in mountains. However, DJF cold bias (−2.23°C to−3.37°C) and JJA warm bias (+1.13°C to +1.87°C) exist in the western United States. The modeling skills of precipitation and temperature at various topographies are similar in the two evaluated regions. The cold bias was likely affected by downward shortwave radiation deficit (−15%) and slow bias of surface wind speeds (−16% to −29%). Despite the wet and cold bias, snow water equivalent and snow cover are underestimated in mountains. We conclude that VR‐CESM outperforms other models evaluated and provides the best estimation of downscaled climate model data for water‐energy studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗

DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2001

Effect of Under-Resolved Grids on High Order Methods

There has been much discussion on verification and validation processes for establishing the credibility of CFD simulations. Since the early 1990s, many of the aeronautical and mechanical engineering related reference journals mandated that any accepted articles in numerical simulations (without known solutions to compared with) need to perform a minimum of one level of grid refinement and time step reduction. Due to the difficulty in analysis, the effect of under-resolved grids and the nonlinear behavior of available spatial discretizations, are scarcely discussed in the literature. Here, an under-resolved numerical simulation is one where the grid spacing being used is too coarse to resolve the smallest physically relevant scales of the chosen continuum governing equations that are of interest to the numerical modeler. With the advent of new developments in fourth-order or higher spatial schemes, it has become common to regard high order schemes as more accurate, reliable and require less grid points. The danger comes when one tries to perform computations with the coarsest grid possible while still hoping to maintain numerical results sufficiently accurate for complex flows, and especially, data-limited problems. On one hand, high order methods when applies to highly coupled multidimensional complex nonlinear problems might have different stability, convergence and reliability behavior than their well studied low order counterparts, especially for nonlinear schemes such as TVD, MUSCL with limiters, ENO, WENO and discrete Galerkin. On the other hand, high order methods involve more operation counts and systematic grid convergence study can be time consuming and prohibitively expansive. At the same time it is difficult to fully understand or categorize the different nonlinear behavior of finite discretizations, especially at the limits of under-resolution when different types of bifurcation phenomena might occur, depending on the combination of grid spacings, time steps, initial conditions and numerical treatments of boundary conditions.

Yee, H. C.↗

NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗

NASA NTRS2005

Receptivity of Supersonic Boundary Layers to Acoustic Disturbances

Boundary layer receptivity to two-dimensional slow and fast acoustic waves is investigated by solving Navier-Stokes equations for Mach 4.5 flow over a flat plate with a finite-thickness leading edge. Higher order spatial and temporal schemes are employed to obtain the solution whereby the flat-plate leading edge region is resolved by providing a sufficiently refined grid. The results show that the instability waves are generated in the leading edge region and that the boundary-layer is much more receptive to slow acoustic waves (by almost a factor of 20) as compared to the fast waves. Hence, this leading-edge receptivity mechanism is expected to be more relevant in the transition process for high Mach number flows. The effect of acoustic wave incidence angle is also studied and it is found that the receptivity of the boundary layer on the windward side (with respect to the acoustic forcing) decreases by more than a factor of 4 when the incidence angle is increased from 0 to 45 deg. However, the receptivity coefficient for the leeward side is found to vary relatively weakly with the incidence angle. The effect of leading-edge thickness is also studied and bluntness is found to stabilize the boundary layer. The relative significance of fast acoustic waves is enhanced in the presence of bluntness.

Malik, Mujeeb R.↗

NASA NTRS2007

Acoustic Receptivity of Mach 4.5 Boundary Layer with Leading- Edge Bluntness

Boundary layer receptivity to two-dimensional slow and fast acoustic waves is investigated by solving Navier-Stokes equations for Mach 4.5 flow over a flat plate with a finite-thickness leading edge. Higher order spatial and temporal schemes are employed to obtain the solution whereby the flat-plate leading edge region is resolved by providing a sufficiently refined grid. The results show that the instability waves are generated in the leading edge region and that the boundary-layer is much more receptive to slow acoustic waves (by almost a factor of 20) as compared to the fast waves. Hence, this leading-edge receptivity mechanism is expected to be more relevant in the transition process for high Mach number flows where second mode instability is dominant. Computations are performed to investigate the effect of leading-edge thickness and it is found that bluntness tends to stabilize the boundary layer. Furthermore, the relative significance of fast acoustic waves is enhanced in the presence of bluntness. The effect of acoustic wave incidence angle is also studied and it is found that the receptivity of the boundary layer on the windward side (with respect to the acoustic forcing) decreases by more than a factor of 4 when the incidence angle is increased from 0 to 45 deg. However, the receptivity coefficient for the leeward side is found to vary relatively weakly with the incidence angle.

Malik, Mujeeb R.↗

NASA NTRS2019

A New Model of the Crustal Magnetic Field of Mars Using MGS and MAVEN

While devoid of an active magnetic dynamo field today, Mars possesses a remanent magnetic field that may reach several thousand nanoteslas locally. The exact origin and the events that have shaped the crustal magnetization remain largely enigmatic. Three magnetic field data sets from two spacecraft collected over 13 cumulative years have sampled the Martian magnetic field over a range of altitudes from 90 up to 6,000 km: (a) Mars Global Surveyor (MGS) magnetometer (1997–2006), (b) MGS Electron Reflectometer (1999–2006), and (c) Mars Atmosphere and Volatile EvolutioN (MAVEN) magnetometer (2014 to today). In this paper we combine these complementary data sets for the first time to build a new model of the Martian internal magnetic field. This new model improves upon previous ones in several aspects: comprehensive data coverage, refined data selection scheme, modified modeling scheme, discrete-to-continuous transformation of the model, and increased model resolution. The new model has a spatial resolution of ∼160 km at the surface, corresponding to spherical harmonic degree 134. It shows small scales and well-defined features, which can now be associated with geological signatures.

Magnetic↗

NASA NTRS2018

Numerical Simulation of Dynamic Stall Using Near-Body Adaptive Mesh Refinement

Time-dependent Reynolds-averaged Navier-Stokes simulations have been carried out for a NACA 0012 airfoil and a flexible UH-60A rotor undergoing dynamic stall. The OVERFLOW computational fluid dynamics (CFD) code is used to explore the use of near-body adaptive mesh refinement (NB-AMR) for the first time on a flexible helicopter rotor in forward flight. Emphasis is placed on understanding the temporal and spatial convergence of the solutions, their numerical efficiency and stability, and establishing a grid-converged solution. Some of the remaining differences between CFD and flight-test measurements are discussed. Time-dependent flow visualization is used to provide an improved understanding of the physical mechanisms involved with two-dimensional and three-dimensional dynamic stall, and the NB-AMR process.

Near-Body↗

NASA NTRS2004

The X-ray halo of an extremely luminous LSB disk galaxy

We are continuing to refine our upper limit on emission from halo gas in Malin 2. The upper limit is, of course, below the detected flux, but is made more difficult to quantify by the disk and possible AGN sources. We are also exploring spectral and spatial-size constraints to help separate the sources of emission. On the theory side, more recent work on the X-ray halo luminosity from halo gas leftover from galaxy formation has lowered the prediction for disk galaxies (e.g. Toft et al. 2002, MNRAS, 335, 799). While our upper limit is well below the original prediction, refinements in model have moved the theoretical goalposts, so that the observation may be consistent with newer models. A recent theoretical development, which our observations of Malin 2 appear to support, is that a substantial amount of mass can be accreted onto galaxies without being heated at a virial shock. The previous standard theory was that gas accreting into a halo hits a virial shock and is heated to high temperatures, which could produce X-ray halos in massive galaxies. Recent models show that "smooth accretion" of matter bypasses the virial shocking (Murali e t al. 2002, ApJ, 571, 1; Birnboim & Dekel 2003, MNRAS, 345, 349). Additionally, new hydrodynamical simulations of galaxy mergers by UCSC graduate student T. J. Cox show that hot gas halos can be created by gas blown out from the merger, taking up orbital energy of the merging galaxies (Cox et al. 2004, ApJ, 607, L87). If mergers rather than virial shocking are the origin of hot gas halos, the existence of an X-ray halo should depend more on past merger activity than halo mass. Then it makes sense that elliptical galaxies and poor groups with ellipticals, which are probably formed in mergers, have X-ray gas halos; while a giant, quiescent LSB disk galaxy like Malin 2, which has never suffered a major merger, does not have an X-ray halo. While both the observational expectations and theoretical models have changed since we began this project, which has forced us to re-evaluate the goals, we are pressing forward to firm up the observations and put them in context of the current models for X-ray halos of massive galaxies.

Weiner, Benjamin J.↗