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 415 records · Page 23

Incompressible Navier-Stokes computations of vortical flows over double-delta wings

An implicit flux-difference splitting scheme is used to compute incompressible laminar vortical flows about a thin round-edged double-delta wing with 80 deg and 60 deg leading-edged sweep for the strake and the main wing respectively. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. It is second-order accurate spatially by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. The interaction between the two primary vortices emanating from the apex and kink leading-edges is successfully simulated. The computed trajectories of vortical cores compare well with experimental data.

Hsu, Chung-Hao↗

High-resolution studies of the structure of the solar atmosphere using a new imaging algorithm

The results of the application of a new image restoration algorithm developed by Ayers and Dainty (1988) to the multiwavelength EUV/Skylab observations of the solar atmosphere are presented. The application of the algorithm makes it possible to reach a resolution better than 5 arcsec, and thus study the structure of the quiet sun on that spatial scale. The results show evidence for discrete looplike structures in the network boundary, 5-10 arcsec in size, at temperatures of 100,000 K.

Karovska, Margarita↗

The F(N) method for the one-angle radiative transfer equation applied to plant canopies

The paper presents a semianalytical solution method, called the F(N) method, for the one-angle radiative transfer equation in slab geometry. The F(N) method is based on two integral equations specifying the intensities exiting the boundaries of the vegetation canopy; the solution is obtained through an expansion in a set of basis functions with expansion coefficients to be determined. The advantage of this method is that it avoids spatial truncation error entirely because it requires discretization only in the angular variable.

Ganapol, B. D.↗

The semi-discrete Galerkin finite element modelling of compressible viscous flow past an airfoil

A method is developed to solve the two-dimensional, steady, compressible, turbulent boundary-layer equations and is coupled to an existing Euler solver for attached transonic airfoil analysis problems. The boundary-layer formulation utilizes the semi-discrete Galerkin (SDG) method to model the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby permitting the use of a uniform finite element grid which provides high resolution near the wall and automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes, through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past NACA 0012 and RAE 2822 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack. All results show good agreement with experiment, and the coupled code proved to be a computationally-efficient and accurate airfoil analysis tool.

Meade, Andrew J., Jr.↗

Digital visual communications using a Perceptual Components Architecture

The next era of space exploration will generate extraordinary volumes of image data, and management of this image data is beyond current technical capabilities. We propose a strategy for coding visual information that exploits the known properties of early human vision. This Perceptual Components Architecture codes images and image sequences in terms of discrete samples from limited bands of color, spatial frequency, orientation, and temporal frequency. This spatiotemporal pyramid offers efficiency (low bit rate), variable resolution, device independence, error-tolerance, and extensibility.

Watson, Andrew B.↗

Silicon Micromachining in RF and Photonic Applications

Texas Instruments (TI) has developed membrane and micromirror devices since the late 1970s. An eggcrate space membrane was used as the spatial light modulator in the early years. Discrete micromirrors supported by cantilever beams created a new era for micromirror devices. Torsional micromirror and flexure-beam micromirror devices were promising for mass production because of their stable supports. TI's digital torsional micromirror device is an amplitude modulator (known as the digital micromirror device (DMD) and is in production development, discussed elsewhere. We also use a torsional device for a 4 x 4 fiber-optic crossbar switch in a 2 cm x 2 cm package. The flexure-beam micromirror device is an analog phase modulator and is considered more efficient than amplitude modulators for use in optical processing systems. TI also developed millimeter-sized membranes for integrated optical switches for telecommunication and network applications. Using a member in radio frequency (RF) switch applications is a rapidly growing area because of the micromechanical device performance in microsecond-switching characteristics. Our preliminary membrane RF switch test structure results indicate promising speed and RF switching performance. TI collaborated with MIT for modeling of metal-based micromachining.

Lin, Tsen-Hwang↗

The Thin Oil Film Equation

A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.

Brown, James L.↗

Reconfigurable Optical Spectra from Perturbations on Elliptical Whispering Gallery Resonances

Elastic strain, electrical bias, and localized geometric deformations were applied to elliptical whispering-gallery-mode resonators fabricated with lithium niobate. The resultant perturbation of the mode spectrum is highly dependant on the modal indices, resulting in a discretely reconfigurable optical spectrum. Breaking of the spatial degeneracy of the whispering-gallery modes due to perturbation is also observed.

filters↗

A Linear-Elasticity Solver for Higher-Order Space-Time Mesh Deformation

A linear-elasticity approach is presented for the generation of meshes appropriate for a higher-order space-time discontinuous finite-element method. The equations of linear-elasticity are discretized using a higher-order, spatially-continuous, finite-element method. Given an initial finite-element mesh, and a specified boundary displacement, we solve for the mesh displacements to obtain a higher-order curvilinear mesh. Alternatively, for moving-domain problems we use the linear-elasticity approach to solve for a temporally discontinuous mesh velocity on each time-slab and recover a continuous mesh deformation by integrating the velocity. The applicability of this methodology is presented for several benchmark test cases.

higher-order↗

DNS of Low-Pressure Turbine Cascade Flows with Elevated Inflow Turbulence Using a Discontinuous-Galerkin Spectral-Element Method

Recent progress towards developing a new computational capability for accurate and efficient high-fidelity direct numerical simulation (DNS) and large-eddy simulation (LES) of turbomachinery is described. This capability is based on an entropy- stable Discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in a computationally efficient manner on a modern high performance computer architecture. An inflow turbulence generation procedure based on a linear forcing approach has been incorporated in this framework and DNS conducted to study the effect of inflow turbulence on the suction- side separation bubble in low-pressure turbine (LPT) cascades. The T106 series of airfoil cascades in both lightly (T106A) and highly loaded (T106C) configurations at exit isentropic Reynolds numbers of 60,000 and 80,000, respectively, are considered. The numerical simulations are performed using 8th-order accurate spatial and 4th-order accurate temporal discretization. The changes in separation bubble topology due to elevated inflow turbulence is captured by the present method and the physical mechanisms leading to the changes are explained. The present results are in good agreement with prior numerical simulations but some expected discrepancies with the experimental data for the T106C case are noted and discussed.

Garai, Anirban↗

Cross-Cutting Computational Modeling Project: Integrative Modeling Approach

A wide range of computational models and analyses have been applied to spaceflight risk assessment and countermeasure development. The benefits of using computational modeling to enhance Human Research Program (HRP) goals include the ability to mathematically represent physiological systems, integrate multiple, discrete experimental measures, span multiple temporal and spatial scales, determine important factors within the system and provide estimates of unmeasurable quantities. In the area of application, computational models provide a means of developing simulations to test hypotheses, determining key factors of the system to aid experimental design and bridging gaps in sparse data by mathematically simulating large populations. Specifically, computational models and their supporting analysis tools have the proven potential to integrate analyses of risk factors to enhance mission planning and preparation capabilities and to inform spacecraft design and countermeasure development. Appropriately applied, computational models may allow intelligent, unbiased physiological parameter assessment to enable hypothesis testing and model based design of experiments. HRP recently formed the Computational Modeling Project (CMP), managed out of Glenn Research Center, as a cross-cutting activity aimed at leveraging the growing power and acceptance of computational modeling in informing clinical, physiological, and biological studies. This presentation will provide an overview of the challenges and opportunities in implementing various forms of computational models in support of the HRPs path to risk reduction.

biomechanics↗

Applicability of Spectral Decomposition by Varimax-Rotated, Principal Component Analysis to the Surface Biology and Geology (SBG) VNIR Mission Concept

Cyanobacterial and Harmful Algal Blooms (CyanoHABs) are a growing concern in coastal and inland waters. But, spectral interference from multiple constituents in optically complex waters can hamper application of remote sensing using traditional image processing methods. The Kent State University (KSU) spectral decomposition method can be applied to multispectral and hyperspectral remote sensing images (e.g. HICO and the NASA Glenn HSI2) to partition and identify signals related to cyanobacteria, algae, pigment degradation products and suspended sediment in each pixel. Fundamental to the use of remote sensing data is the ability to extract independent signals from correlated hyperspectral VNIR data cubes. The Kent State University varimax-rotated, principal component analysis method (VPCA) is important to integrate into the SBG VNIR mission concept because it provides greater specificity, a software-based SNR boost relative to hardware performance, and can assist with Cal/Val, Modeling and Applications. We present examples of the hyperspectral application of the KSU VPCA method with relevance to SBG. The information extracted by VPCA can be validated spectrally or spatially with laboratory and/or in situ sensors, which capture spatial or time series of information at discrete points within remote sensing images. Comparisons show hyperspectral sensors extract more components than multispectral ones, but more independent information can be extracted from multispectral sensors by VPCA than traditional band ratio approaches. The spectral decomposition method is capable of enhancing the signal to noise ratio (SNR) of the NASA Glenn, second-generation hyperspectral imager by a factor of 7x to 20x, with a spectral reproducibility of ±3%. The spectral decomposition method, when compared against existing remote sensing monitoring methods exhibits both greater specificity and a lower detection limit. The method has been validated with multispectral images in Lake Erie to quantify the Microcystis CyanoHAB and from the Indian River Lagoon, Florida to quantify the Brown Tide resulting from A. lagunesnsis. Field operations in the Western Basin of Lake Erie were conducted using a bbe Fluoroprobe to collect vertical profiles and horizontal tows along a transect from the Toledo to the Detroit Lighthouse during coincident satellite overpasses. Extraction of pixel values from the MODIS Aqua sensor yields agreement between in situ field and lab-based measures of cyanobacterial, cryptophyte, diatoms and green algae, suspended sediment and pigment degradation products with R2>0.8.

hyperspectral↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗

Visual sensitivity to spatially sampled modulation in human observers

Thresholds were measured for detecting spatial luminance modulation in regular lattices of visually discrete dots. Thresholds for modulation of a lattice are generally higher than the corresponding threshold for modulation of a continuous field, and the size of the threshold elevation, which depends on the spacing of the lattice elements, can be as large as a one log unit. The largest threshold elevations are seen when the sample spacing is 12 min arc or greater. Theories based on response compression cannot explain the further observation that the threshold elevations due to spatial sampling are also dependent on modulation frequency: the greatest elevations occur with higher modulation frequencies. The idea that this is due to masking of the modulation frequency by the spatial frequencies in the sampling lattice is considered.

Mulligan, Jeffrey B.↗

A numerical extension of the spatially-filtered Euler equations for contact discontinuities

Solving the spatially-filtered Euler equations with a kinetic energy-preserving centered discretization is an efficient framework for simulating shock-dominated flows. Here, the present work describes a numerical treatment of the spatially-filtered Euler equations to minimize oscillations when simulating contact discontinuities. To counteract dispersive errors inherent to centered schemes, a WENO-like correction term is applied to the enthalpy transport in the energy equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Intercalibration of the reflective solar bands of MODIS and MISR instruments on the Terra platform

The multispectral imaging sensors on the Terra platform have been operating for over two decades facilitating a variety of scientific applications. The MODIS sensor provides the largest spectral coverage of 0.41 to 14.2 μm, acquiring data at three different spatial resolutions, 250 m, 500 m, and 1 km. The MISR instrument views the Earth using nine discrete cameras pointed at fixed angles including viewing the nadir direction at a spatial resolution of 275 m and covering a wavelength range from 0.44 to 0.86 μm. Being on the same platform, the two sensors complement each other in terms of spatial coverage (and target viewing geometry) and facilitate synergistic applications using multispectral data. A consistent radiometric calibration between these sensors is a prerequisite for creating high quality science products from their observations. Both instruments underwent intensive prelaunch characterization, with calibration monitored on-orbit using their onboard calibrators. In this paper, we perform a calibration inter-comparison of the spectrally matching bands of the two instruments using vicarious techniques. Vicarious techniques include multiyear simultaneous views of the North African desert, North Atlantic Ocean and Dome Concordia, thereby covering the different parts of the dynamic range. Also included in this work are the near-simultaneous top-of-atmosphere (TOA) reflectance measurements from Railroad Valley, USA, as provided by the RadCalNet (converted to TOA), that are used as a calibration reference to compare the on-orbit observations between MODIS and MISR.

Intercalibration↗