Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Dynamical Approach Study of Spurious Numerics in Nonlinear Computations

The last two decades have been an era when computation is ahead of analysis and when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems and non-traditional fields. Ensuring a higher level of confidence in the predictability and reliability (PAR) of these numerical simulations could play a major role in furthering the design, understanding, affordability and safety of our next generation air and space transportation systems, and systems for planetary and atmospheric sciences, and in understanding the evolution and origin of life. The need to guarantee PAR becomes acute when computations offer the ONLY way of solving these types of data limited problems. Employing theory from nonlinear dynamical systems, some building blocks to ensure a higher level of confidence in PAR of numerical simulations have been revealed by the author and world expert collaborators in relevant fields. Five building blocks with supporting numerical examples were discussed. The next step is to utilize knowledge gained by including nonlinear dynamics, bifurcation and chaos theories as an integral part of the numerical process. The third step is to design integrated criteria for reliable and accurate algorithms that cater to the different multiscale nonlinear physics. This includes but is not limited to the construction of appropriate adaptive spatial and temporal discretizations that are suitable for the underlying governing equations. In addition, a multiresolution wavelets approach for adaptive numerical dissipation/filter controls for high speed turbulence, acoustics and combustion simulations will be sought. These steps are corner stones for guarding against spurious numerical solutions that are solutions of the discretized counterparts but are not solutions of the underlying governing equations.

Yee, H. C.↗
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗
NASA NTRS2014

Comparing and Reconciling Traditional Field and Photogeologic Mapping Techniques: Lessons from the San Francisco Volcanic Field, Arizona

Cartographic products and - specifically - geologic maps provide critical assistance for establishing physical and temporal frameworks of planetary surfaces. The technical methods that result in the creation of geologic maps vary depending on how observations are made as well as the overall intent of the final products [1-3]. These methods tend to follow a common linear work flow, including the identification and delineation of spatially and temporally discrete materials (units), the documentation of their primary (emplacement) and secondary (erosional) characteristics, analysis of the relative and absolute age relationships between these materials, and the collation of observations and interpretations into an objective map product. The "objectivity" of a map is critical cross comparison with overlapping maps and topical studies as well as its relevance to scientific posterity. However, the "accuracy" and "correctness" of a geologic map is very subject to debate. This can be evidenced by comparison of existing geologic maps at various scales, particularly those compiled through field- and remote-based mapped efforts. Our study focuses on comparing the fidelity of (1) "Apollo-style" geologic investigations, where typically non-geologist crew members follow static traverse routes established through pre-mission planning, and (2) "traditional" field-based investigations, where geologists are given free rein to observe without preplanned routes. This abstract summarizes the regional geology wherein our study was conducted, presents the geologic map created from traditional field mapping techniques, and offers basic insights into how geologic maps created from different tactics can be reconciled in support of exploratory missions. Additional abstracts [4-6] from this study discuss various exploration and science results of these efforts.

Skinner, J. A., Jr.↗
DOE OSTI2023

Computational Fluid Dynamics Study of a Cross-Flow Marine Hydrokinetic Turbine and the Combined Influence of Struts and Helical Blades: Preprint

A computational fluid dynamics study was performed for a cross-flow marine hydro-kinetic turbine. The analysis was done in three dimensions and used the unsteady Reynolds averaged Navier-Stokes solver in the commercial code STAR-CCM+. The base turbine configuration is the RivGen® Turbine, designed by the Ocean Renewable Power Company (ORPC). A convergence and uncertainty analysis was performed for both the spatial and temporal discretization; this was done using the base configuration which features support struts and helical foils. The proposed study aims to compare the impact of the struts on both power performance and blade loading for helical and straight blades.

CFD↗
DOE OSTI2020

Coarse-grain cluster analysis of tensors with application to climate biome identification

A tensor provides a concise way to codify the interdependence of complex data. Treating a tensor as a d-way array, each entry records the interaction between the different indices. Clustering provides a way to parse the complexity of the data into more readily understandable information. Clustering methods are heavily dependent on the algorithm of choice, as well as the chosen hyperparameters of the algorithm. However, their sensitivity to data scales is largely unknown. In this work, we apply the discrete wavelet transform to analyze the effects of coarse-graining on clustering tensor data. We are particularly interested in understanding how scale affects clustering of the Earth's climate system. The discrete wavelet transform allows classification of the Earth's climate across a multitude of spatial-temporal scales. The discrete wavelet transform is used to produce an ensemble of classification estimates, as opposed to a single classification. Each element of the ensemble is a clustering at a different spatial-temporal scale. Information theoretic approaches are used to identify important scale lengths in clustering the L15 Climate Datset. We also discover a sub-collection of the ensemble that spans the majority of the variance observed, allowing for efficient consensus clustering techniques that can be used to identify climate biomes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

The stability of numerical boundary treatments for compact high-order finite-difference schemes

The stability characteristics of various compact fourth and sixth order spatial operators are assessed using the theory of Gustafsson, Kreiss and Sundstrom (G-K-S) for the semi-discrete Initial Boundary Value Problem (IBVP). These results are then generalized to the fully discrete case using a recently developed theory of Kreiss. In all cases, favorable comparisons are obtained between the G-K-S theory, eigenvalue determination, and numerical simulation. The conventional definition of stability is then sharpened to include only those spatial discretizations that are asymptotically stable. It is shown that many of the higher order schemes which are G-K-S stable are not asymptotically stable. A series of compact fourth and sixth order schemes, which are both asymptotically and G-K-S stable for the scalar case, are then developed.

Carpenter, Mark H.↗
NASA NTRS2019

Turbulence Model Implementation and Verification in the SENSEI CFD Code

This paper outlines the implementation and verification of the negative Spalart-Allmaras turbulence model into the SENSEI CFD code. The SA-neg turbulence model is implemented in a flexible, object-oriented framework where additional turbulence models can be easily added. In addition to outlining the new turbulence modeling framework in SENSEI, an overview of the other general improvements to SENSEI is provided. The results for four 2D test cases are compared to results from CFL3D and FUN3D to verify that the turbulence models are implemented properly. Several differences in the results from SENSEI, CFL3D, and FUN3D are identified and are attributed to differences in the implementation and discretization order of the boundary conditions as well as the order of discretization of the turbulence model. When a solid surface is located near or intersects an inflow or outflow boundary, higher order boundary conditions should be used to limit their effect on the forces on the surface. When the turbulence equations are discretized using second order spatial accuracy, the edge of the eddy viscosity profile seems to be sharper than when a first order discretization is used. However, the discretization order of the turbulence equation does not have a significant impact on output quantities of interest, such as pressure and viscous drag, for the cases studied.

Jackson, Charles W.↗
NASA NTRSDate not supplied

Development of Multi-Edge Slant Target for Unlocalized MTF Measurement of Airborne Imaging System Payloads

Calibration of imaging instruments is a critical step in effectively characterizing imaging data. A comprehensive characterization of an imaging system provides valuable insight into the overall uncertainties in the results derived from the raw data collected. The Scientifically Calibrated In-Flight Imagery (SCIFLI) team (based at NASA Langley Research Center) has designed, built, tested, and successfully flown the optical payload known as the SCIFLI Airborne Multispectral Imager (SAMI). SAMI is a multispectral imaging payload, complete with sensor configurations for the Ultraviolet-Visible spectrum as well as Near- Infrared, Shortwave Infrared, and Midwave Infrared bands and is fitted with various optical hardware designed to turn, split, and focus light beams to the respective in-band sensors. Radiometric calibrations are performed for SAMI such that the accuracy and uncertainty of thermal and spectral measurements are effectively characterized. However, SAMI does not regularly undergo spatial characterizations to evaluate the performance of the imager at high spatial frequencies. This is essential, as understanding the spatial sensor response of discrete sensors is just as important as the spectral sensor response. The payload currently utilizes a calibration cart for lab testing. Equipped with a modular setup and several radiation sources, it is perfect for performing high-quality radiometric calibrations. This work focuses on using the calibration cart to perform spatial calibrations. A custom designed calibration target is being used to evaluate the optical performance of SAMI. Specifically, the target is designed to measure the Modulation Transfer Function (MTF) of the imager in discrete locations around the focal plane. The results of this characterization will provide insight into the spatial resolution capabilities of the SAMI payload in practice. Future work includes using the results of the study to correct optical imperfections in imagery taken with SAMI.

SAMI↗
DOE OSTI2025

$2+1$ dimensional Floquet systems and lattice fermions: Exact bulk spectral equivalence

A connection has recently been proposed between periodically driven systems known as Floquet insulators in continuous time and static fermion theories in discrete time. This connection has been established in a (1+1) ( 1 + 1 ) -dimensional free theory, where an explicit mapping between the spectra of a Floquet insulator and a discrete-time Dirac fermion theory has been formulated. Here we investigate the potential of static discrete-time theories to capture Floquet physics in higher dimensions, where so-called anomalous Floquet topological insulators can emerge that feature chiral edge states despite having bulk bands with zero Chern number. Starting from a particular model of an anomalous Floquet system, we provide an example of a static discrete-time theory whose bulk spectrum is an exact analytic match for the Floquet spectrum. The spectra with open boundary conditions in a particular strip geometry also match up to finite-size corrections. However, the models differ in several important respects. The discrete-time theory is spatially anisotropic, so that the spectra do not agree for all lattice terminations, e.g. other strip geometries or on half spaces. This difference can be attributed to the fact that the static discrete-time model is quasi-one-dimensional in nature and therefore has a different bulk-boundary correspondence than the Floquet model.

Iadecola, Thomas (ORCID:0000000251456441)↗
NASA NTRS2019

Retrieving Aerosol Characteristics From the PACE Mission, Part 1: Ocean Color Instrument

NASA’s Plankton, Aerosol, Clouds, ocean Ecosystem (PACE) satellite mission is scheduled to launch in 2022, with the Ocean Color Instrument (OCI) on board. For the first time reflected sunlight from the Earth across a broad spectrum from the ultraviolet (UV: 350 nm) to the short wave infrared (SWIR: 2260 nm) will be measured from a single instrument at 1 km spatial resolution. While seven discrete bands will represent the SWIR, the spectrum from 350 to 890 nm will be continuously covered with a spectral resolution of 5 nm. OCI will thus combine in a single instrument (and at an enhanced spatial resolution for the UV) the heritage capabilities of the Moderate resolution Imaging Spectroradiometer (MODIS) and the Ozone Monitoring Instrument (OMI), while covering the oxygen A-band (O2A). Designed for ocean color and ocean biology retrievals, OCI also enables continuation of heritage satellite aerosol products and the development of new aerosol characterization from space. In particular the combination of MODIS and OMI characteristics allows deriving aerosol height, absorption and optical depth along with a measure of particle size distribution. This is achieved by using the traditional MODIS visible-to-SWIR wavelengths to constrain spectral aerosol optical depth and particle size. Extrapolating this information to the UV channels allows retrieval of aerosol absorption and layer height. A more direct method to derive aerosol layer height makes use of O2A absorption methods, despite the relative coarseness of the nominal 5 nm spectral resolution of OCI. Altogether the PACE mission with OCI will be an unprecedented opportunity for aerosol characterization that will continue climate data records from the past decades and propel aerosol science forward toward new opportunities.

PACE↗
NASA NTRS2017

On the Relationship Between Sub-Daily Instantaneous and Daily Total Gross Primary Production: Implications for Interpreting Satellite-Based SIF Retrievals

Spatially and temporally continuous estimation of plant photosynthetic carbon fixation (or gross primary production, GPP) is crucial to our understanding of the global carbon cycle and the impact of climate change. Besides spatial, seasonal and interannual variations, GPP also exhibits strong diurnal variations. Satellite retrieved solar-induced chlorophyll fluorescence (SIF) provides a spatially continuous, but temporally discrete measurement of plant photosynthesis, and has the potential to be used to estimate GPP at global scale. However, it remains unclear whether the seasonal time series of SIF snapshots taken at a fixed time of the day can be used to infer daily total GPP variation at spatial and seasonal scales. In this study, we first used GPP estimates from 135 eddy covariance flux sites, covering a wide range of geographic locations and biome types, to investigate the relationship between the instantaneous GPP (GPP(sub inst)) and daily GPP (GPP(sub daily)) on the seasonal course for different times of the day. Latitudinal and diurnal patterns were found to correspond to variations in photosynthetically active radiation (PAR) and light use efficiency (LUE), respectively. We then used the Soil-Canopy Observation Photosynthesis and Energy Balance (SCOPE) model and the FluxCom GPP product to investigate the instantaneous and daily SIF-GPP relationships at five flux tower sites along a latitudinal gradient and at a global scale for different biome types. The results showed that daily SIF had a stronger linear correlation with daily GPP than instantaneous SIF at the seasonal scale, with an instantaneous to daily SIF conversion factor following the latitudinal and seasonal pattern driven by PAR. Our study highlights the necessity to take the latitudinal and diurnal factors into consideration for SIF-GPP relationship analyses or for physiological phenology analyses based on SIF.

photosynthesis model↗
NASA NTRS1991

Rotational wind indicator enhances control of rotated displays

Rotation by 108 deg of the spatial mapping between a visual display and a manual input device produces large spatial errors in a discrete aiming task. These errors are not easily corrected by voluntary mental effort, but the central nervous system does adapt gradually to the new mapping. Bernotat (1970) showed that adding true hand position to a 90 deg rotated display improved performance of a compensatory tracking task, but tracking error rose again upon removal of the explicit cue. This suggests that the explicit error signal did not induce changes in the neural mapping, but rather allowed the operator to reduce tracking error using a higher mental strategy. In this report, we describe an explicit visual display enhancement applied to a 108 deg rotated discrete aiming task. A 'wind indicator' corresponding to the effect of the mapping rotation is displayed on the operator-controlled cursor. The human operator is instructed to oppose the virtual force represented by the indicator, as one would do if flying an airplane in a crosswind. This enhancement reduces spatial aiming error in the first 10 minutes of practice by an average of 70 percent when compared to a no enhancement control condition. Moreover, it produces adaptation aftereffect, which is evidence of learning by neural adaptation rather than by mental strategy. Finally, aiming error does not rise upon removal of the explicit cue.

Cunningham, H. A.↗
DOE OSTI2025

Shearing approach to gauge-invariant Trotterization

Universal quantum simulations of gauge field theories are exposed to the risk of gauge symmetry violations when it is not known how to compile the desired operations exactly using the available gate set. In this article, we show how time evolution can be compiled in an Abelian gauge theory—if only approximately—without compromising gauge invariance, by graphically motivating a block-diagonalization procedure. When gauge-invariant interactions are associated with a “spatial network” in the space of discrete quantum numbers, it is seen that cyclically shearing the spatial network converts simultaneous updates to many quantum numbers into conditional updates of a single quantum number; ultimately, this eliminates any need to pass through (and acquire overlap onto) unphysical intermediate configurations. Shearing is explicitly applied to gauge-matter and magnetic interactions of lattice quantum electrodynamics. The features that make shearing successful at preserving Abelian gauge symmetry may also be found in non-Abelian theories, bringing one closer to gauge-invariant simulations of quantum chromodynamics.

Gauge theories↗
DOE OSTI2022

Spatiotemporal Plasma-Particle Characterization of Dry Aerosols Using Nanosecond, Femtosecond, and Filament Laser-Produced Plasmas

The ability to rapidly characterize dry aerosols in air using laser-induced breakdown spectroscopy (LIBS) with femtosecond laser pulses promises advancement towards real-time atmospheric sampling and standoff capabilities. Of particular interest is the ability to apply LIBS in the context of low-particle loaded environments where discrete particle interactions must be observed within the sampling volume of the laser-produced plasma (LPP). In this study, dry nanoparticles in suspension are generated from a standard solution and sampled in air using Q-switched nanosecond (ns-) pulses, short-focus (SF) femtosecond (fs-) pulses, and filaments. Short time-gated plasma images are captured to observe spatially and temporally varying discrete plasma-particle interactions, which is shown to influence early air breakdown behavior and subsequent plasma evolution. Along with images, photo-multiplier tube (PMT) measurements are conducted where strong spatiotemporal dependencies are exhibited by the collected emission signal on particle proximity and plasma expansion behavior. Finally, conditional analysis is performed on LIBS measurements to determine associated sampling probabilities and filter out spectra with poor or absent emission peaks with an adaptive threshold algorithm.

42 ENGINEERING↗
NASA NTRS2020

Spectral Spherical Harmonics Discrete Ordinate Method

A new method for modeling the radiative transfer in inhomogeneous three-dimensional media illuminated by a Gaussian beam is described. This approach, called the Spectral Spherical Harmonics Discrete Ordinate Method (SSHDOM), uses the Fourier expansion method to transform the three-dimensional radiative transfer into an one-dimensional equation in the spectral domain, and the Spherical Harmonics Discrete Ordinate Method (SHDOM) for its solution. Specifically, (i) the source function is represented in the spectral domain through a spherical harmonic expansion, (ii) the spectral one-dimensional radiative transfer equation is integrated along discrete ordinates through a spatial grid, and (iii) the solution method is based on the Picard iteration. Both SSHDOM and SHDOM algorithms are implemented in a common computer code.

Gaussian beam↗
NASA NTRS2008

The Foggy EUV Corona and Coronal Heating by MHD Waves from Explosive Reconnection Events

In 0.5 arcsec/pixel TRACE coronal EUV images, the corona rooted in active regions that are at the limb and are not flaring is seen to consist of (1) a complex array of discrete loops and plumes embedded in (2) a diffuse ambient component that shows no fine structure and gradually fades with height. For each of two not-flaring active regions, found that the diffuse component is (1) approximately isothermal and hydrostatic and (2) emits well over half of the total EUV luminosity of the active-region corona. Here, from a TRACE Fe XII coronal image of another not-flaring active region, the large sunspot active region AR 10652 when it was at the west limb on 30 July 2004, we separate the diffuse component from the discrete loop component by spatial filtering, and find that the diffuse component has about 60% of the total luminosity. If under much higher spatial resolution than that of TRACE (e. g., the 0.1 arcsec/pixel resolution of the Hi-C sounding-rocket experiment proposed by J. W. Cirtain et al), most of the diffuse component remains diffuse rather being resolved into very narrow loops and plumes, this will raise the possibility that the EUV corona in active regions consists of two basically different but comparably luminous components: one being the set of discrete bright loops and plumes and the other being a truly diffuse component filling the space between the discrete loops and plumes. This dichotomy would imply that there are two different but comparably powerful coronal heating mechanisms operating in active regions, one for the distinct loops and plumes and another for the diffuse component. We present a scenario in which (1) each discrete bright loop or plume is a flux tube that was recently reconnected in a burst of reconnection, and (2) the diffuse component is heated by MHD waves that are generated by these reconnection events and by other fine-scale explosive reconnection events, most of which occur in and below the base of the corona where they are seen as UV explosive events, EUV blinkers, and type II spicules. These MHD waves propagate across field lines and dissipate, heating the plasma in the field between the bright loops and plumes.

Moore, Ron L.↗
NASA NTRS1992

Development of an on-line parameter estimation system using the discrete modal filter

Ongoing development of an automated system for monitoring the frequency and damping of time varying structural systems is discussed. A spatial filtering technique called the discrete modal filter is used to uncouple multiple response measurements into modal coordinate responses. Since the modal coordinate responses are the responses of single-degree-of-freedom systems, the frequency and damping may be accurately estimated with short data records. An on-line monitoring system is discussed which may be useful in a variety of applications. The focus in this paper is on the application to aircraft flight flutter testing.

Shelley, S. J.↗