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 307 records · Page 17

High Spatial Resolution Airborne Multispectral Thermal Infrared Remote Sensing Data for Analysis of Urban Landscape Characteristics

We have used airborne multispectral thermal infrared (TIR) remote sensing data collected at a high spatial resolution (i.e., 10m) over several cities in the United States to study thermal energy characteristics of the urban landscape. These TIR data provide a unique opportunity to quantify thermal responses from discrete surfaces typical of the urban landscape and to identify both the spatial arrangement and patterns of thermal processes across the city. The information obtained from these data is critical to understanding how urban surfaces drive or force development of the Urban Heat Island (UHI) effect, which exists as a dome of elevated air temperatures that presides over cities in contrast to surrounding non-urbanized areas. The UHI is most pronounced in the summertime where urban surfaces, such as rooftops and pavement, store solar radiation throughout the day, and release this stored energy slowly after sunset creating air temperatures over the city that are in excess of 2-4'C warmer in contrast with non-urban or rural air temperatures. The UHI can also exist as a daytime phenomenon with surface temperatures in downtown areas of cities exceeding 38'C. The implications of the UHI are significant, particularly as an additive source of thermal energy input that exacerbates the overall production of ground level ozone over cities. We have used the Airborne Thermal and Land Applications Sensor (ATLAS), flown onboard a Lear 23 jet aircraft from the NASA Stennis Space Center, to acquire high spatial resolution multispectral TIR data (i.e., 6 bandwidths between 8.2-12.2 (um) over Huntsville, Alabama, Atlanta, Georgia, Baton Rouge, Louisiana, Salt Lake City, Utah, and Sacramento, California. These TIR data have been used to produce maps and other products, showing the spatial distribution of heating and cooling patterns over these cities to better understand how the morphology of the urban landscape affects development of the UHI. In turn, these data have been used by government officials, urban planners, and other decision-makers, to make more informed decisions on how to mitigate the UHI and its subsequent impacts.

Quattrochi, Dale A.↗

Aliasing and blurring in 2-D sampled imagery

The quality of image reconstructions from discrete data suffers not only from the blurring of spatial detail caused by limitations in the spatial frequency response of electrooptical systems, but also from the aliasing generated if spatial detail has been undersampled. Previous studies have observed that reasonable spot intensity profiles and photosensor aperture shapes of equivalent size result in about equal blurring but that some profiles and shapes suppress aliasing better than others. This paper presents quantitative results of the magnitude of aliasing and blurring as a function of random radiance fields typical for natural scenes and of spatial responses and sampling intervals typical for TV cameras and optical-mechanical scanners. These results indicate that aliasing may often be a larger source of degradation than either blurring or electronic noise.

Huck, F. O.↗

Do clusters of galaxies affect the spectrum of the microwave background?

The distortion, averaged over the sky, due to the Compton scattering of background photons with electrons in the hot gas in clusters of galaxies is calculated. Using an existing sample of X-ray clusters, various values of the density parameter Omega, and plausible models for cluster evolution, Monte Carlo realizations of the microwave sky are generated. The spatial structure of these simulations shows a network of discrete sources whose properties can be a strong function of both evolution and Omega. The amount of spectral distortion in the models is greatest for models characterized by self-similar cluster evolution in an open universe and is within an order of magnitude of the current upper limits. Thus, improved observational sensitivity must inevitably detect some deviation from a blackbody spectrum.

Markevitch, M.↗

Remote Observations of Reentering Spacecraft Including the Space Shuttle Orbiter

Flight measurement is a critical phase in development, validation and certification processes of technologies destined for future civilian and military operational capabilities. This paper focuses on several recent NASA-sponsored remote observations that have provided unique engineering and scientific insights of reentry vehicle flight phenomenology and performance that could not necessarily be obtained with more traditional instrumentation methods such as onboard discrete surface sensors. The missions highlighted include multiple spatially-resolved infrared observations of the NASA Space Shuttle Orbiter during hypersonic reentry from 2009 to 2011, and emission spectroscopy of comparatively small-sized sample return capsules returning from exploration missions. Emphasis has been placed upon identifying the challenges associated with these remote sensing missions with focus on end-to-end aspects that include the initial science objective, selection of the appropriate imaging platform and instrumentation suite, target flight path analysis and acquisition strategy, pre-mission simulations to optimize sensor configuration, logistics and communications during the actual observation. Explored are collaborative opportunities and technology investments required to develop a next-generation quantitative imaging system (i.e., an intelligent sensor and platform) with greater capability, which could more affordably support cross cutting civilian and military flight test needs.

Horvath, Thomas J.↗

Indigenous Organic Matter Associated with Unusual Low Temperature Aqueous Secondary Phases in the Mars Meteorite Nakhla

In freshly fractured chips of the Mars meteorite Nakhla, we found carbonaceous phases intimately associated with secondary aqueous alteration phases, both mineral and amorphous, interpreted to have formed through the low-temperature, mineral-water dissolution of the host Martian basalt.We have characterized the nature of carbonaceous matter and spatially associated secondary phases within Nakhla using a multidisciplinary suite of optical and electron beam instrumentation(Thomas-Keprta et al., 2022). Secondary Phases: The preserved secondary mineral phases are best described as ‘juvenile’ hydrated amorphous silica formed when water interacts with mafic rocks. Although amorphous silica on Earth rarely persists through geologic time due to alteration into other phases, e.g.,microcrystalline opal-CT and eventually quartz, the lack of maturation of primarily amorphous silica with embedded nano phase ferrihydrite in Nakhla infers these phases formed in a limited water environment by low temperature, circumneutral pH fluids stabilized by the presence of coexisting amorphous silica.Additionally, we report the presence of the ferrous hydroxycarbonate mineral chukanovite (Fe2(OH)2CO3)within some of the secondary mineral aggregates studied. Organic Matter:Indigenous Nakhla carbonaceous matter is present both in condensed phases and in a dispersed state spatially associated with secondary alteration phases. In the former, discrete refractory micron to submicron assemblages were identified that appear macromolecular in nature and, in several cases, associated with F and, in one case, significant N. While organically bound Fhas yet to be detected remotely on Mars, the identification of organic N within one of the carbonaceous phases is considered a prerequisite for Martian astrobiological potential. Summary: Neither the intimate association of carbonaceous matter with secondary phases nor the identification of chukanovite have been previously reported in any of the Martian meteorites.A variety of habitability-related metastable solid-phase, oxidation-reduction, and organic-molecular sample attributes in Nakhla formed hundreds of millions of years ago near Mars’ surface; these alteration products persisted on Mars until very recentlyand resemble the amorphous materials observed at Gale crater by Mars Science Laboratory CuriosityRover (Rampe et al., 2020).Similar phases may be expected to occur elsewhere on Mars (e.g., Jezerocrater). References: Rampe E. B. et al.(2020) Mineralogy and geochemistry of sedimentary rocks and eolian sediments in Gale crater, Mars: A review after six Earth years of exploration with Curiosity.Geochemistry80, 125605; Thomas-Keprta K.L. et al. (2022) Indigenous carbon-rich matter associated with unusual aqueous alteration features in Nakhla: Complex formation and preservation history. Geochim. et Cosmochim. Acta320, 41-78 (online).

Mars meteorite↗

The Need for High Spatial Resolution Multispectral Thermal Remote Sensing Data In Urban Heat Island Research

Although the study of the Urban Heat Island (UHI) effect dates back to the early 1800's when Luke Howard discovered London s heat island, it has only been with the advent of thermal remote sensing systems that the extent, characteristics, and impacts of the UHI have become to be understood. Analysis of the UHI effect is important because above all, this phenomenon can directly influence the health and welfare of urban residents. For example, in 1995, over 700 people died in Chicago due to heat-related causes. UHI s are characterized by increased temperature in comparison to rural areas and mortality rates during a heat wave increase exponentially with the maximum temperature, an effect that is exacerbated by the UHI. Aside from the direct impacts of the UHI on temperature, UHI s can produce secondary effects on local meteorology, including altering local wind patterns, increased development of clouds and fog, and increasing rates of precipitation either over, or downwind, of cities. Because of the extreme heterogeneity of the urban surface, in combination with the sprawl associated with urban growth, thermal infrared (TIR) remote sensing data have become of significant importance in understanding how land cover and land use characteristics affect the development and intensification of the UHI. TIR satellite data have been used extensively to analyze the surface temperature regimes of cities to help observe and measure the impacts of surface temperatures across the urban landscape. However, the spatial scales at which satellite TIR data are collected are for the most part, coarse, with the finest readily available TIR data collected by the Landsat ETM+ sensor at 60m spatial resolution. For many years, we have collected high spatial resolution (10m) data using an airborne multispectral TIR sensor over a number of cities across the United States. These high resolution data have been used to develop an understanding of how discrete surfaces across the urban environment (e.g., rooftops, pavements) interact from a surface-lower atmosphere energy flux perspective, to force the development of the UHI. Moreover, the airborne TIR sensor we used in our UHI studies was a multispectral sensor that had six channels in the 8-12pm range. The advantages of collecting multispectral TIR data became readily evident as a valuable tool for better calculation of unique surface thermal energy responses for urban materials over the 8-12 micrometer region, and also for getting a better handle on surface emissivity characteristics for these discrete surfaces. In this presentation, we will provide evidence on the virtues of how high spatial resolution multispectral TIR data can provide for better analysis of the UHI that cannot now be attained via TIR data obtained from satellites. Furthermore, we wish to provide compelling evidence on why future TIR satellite sensors should collect data at fine spatial resolutions (e.g. less than or equal to 30m) to better allow for measurement of surface thermal energy fluxes from discrete urban surfaces, and to better understand how surface fluxes from different urban materials in cities around the world in different climatic regimes, affect development of the UHI characteristics.

Quattrochi, Dale A.↗

Modeling of supersonic reacting flow fields

A detailed understanding of the scramjet combustor flow field is critical to the achievement of a successful design. Even though the combustor flow field is quite complex, it can be realistically viewed as a collection of spatially developing and reacting supersonic mixing layers that are initially discrete, but that ultimately merge into larger more complex zones. These mixing layers begin downstream of a set of fuel injectors that introduce gaseous hydrogen in both a parallel and transverse direction into a supersonic air stream entering from the engine inlet. The behavior of the initial portion of the combustor flow, in the mixing layers near the fuel injectors, appears to be most critical, since this is where the mechanism for efficient high speed mixing must be established to achieve the required degree of combustion downstream. Because of the structure of the flow field in this initial portion of the combustor, a single supersonic, spatially developing and reacting mixing layer serves as an excellent physical model for the overall flow field. Even though this reacting mixing layer flow is geometrically simple, it can still be made to retain all of the fluid mechanical and chemical complexities present in the actual combustor flow field.

Drummond, J. P.↗

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗

A Stable and Conservative Interface Treatment of Arbitrary Spatial Accuracy

Stable and accurate interface conditions are derived for the linear advection-diffusion equation. The conditions are functionally independent of the spatial order of accuracy and rely only on the form of the discrete operator. We focus on high-order finite-difference operators that satisfy the summation-by-parts (SBP) property. We prove that stability is a natural consequence of the SBP operators used in conjunction with the new boundary conditions. In addition, we show that the interface treatments are conservative. New finite-difference operators of spatial accuracy up to sixth order are constructed: these operators satisfy the SBP property. Finite-difference operators are shown to admit design accuracy (p(sup th)-order global accuracy) when (p - 1)(sup th)-order stencil closures are used near the boundaries if the physical boundary conditions are implemented to at least p(sup th)-order accuracy. Stability and accuracy are demonstrated on the nonlinear Burgers' equation for an twelve-subdomain problem with randomly distributed interfaces.

Carpenter, Mark H.↗

Aliased noise in radiometric measurements

The magnitude of aliased noise that degrades the accuracy of continuous reconstructions of discrete radiometric measurements was evaluated as a function of the spatial response and sampling intervals of the radiometer, and of the resolution of the reconstructed measurements. A Wiener spectrum, representative of a wide range of scenes, was used to characterize the radiance fluctuations.

Huck, F. O.↗

Structure of the midcontinent basement. Topography, gravity, seismic, and remote sensing

Some 600,000 discrete Bouguer gravity estimates of the continental United States were spatially filtered to produce a continuous tone image. The filtered data were also digitally painted in color coded form onto a shaded relief map. The resultant image is a colored shaded relief map where the hue and saturation of a given image element is controlled by the value of the Bouguer anomaly. Major structural features (e.g., midcontinent gravity high) are readily discernible in these data, as are a number of subtle and previously unrecognized features. A linear gravity low that is approximately 120 to 150 km wide extends from southeastern Nebraska, at a break in the midcontinent gravity high, through the Ozark Plateau, and across the Mississippi embayment. The low is also aligned with the Lewis and Clark lineament (Montana to Washington), forming a linear feature of approximately 2800 km in length. In southeastern Missouri the gravity low has an amplitude of 30 milligals, a value that is too high to be explained by simple valley fill by sedimentary rocks.

Guinness, E. A.↗

Three-dimensional incompressible Navier-Stokes simulations of slender-wing vortices

The application of an implicit flux-difference splitting scheme to leading-edge vortex flows is described. The numerical scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. Second-order spatial accuracy is achieved by applying TVD-like upwind discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Good agreement between computed results and experimental data has been obtained for three different low-aspect-ratio wings.

Hsu, Chung-Hao↗

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

Matrix-based, finite-difference algorithms for computational acoustics

A compact numerical algorithm is introduced for simulating multidimensional acoustic waves. The algorithm is expressed in terms of a set of matrix coefficients on a three-point spatial grid that approximates the acoustic wave equation with a discretization error of O(h exp 5). The method is based on tracking a local phase variable and its implementation suggests a convenient coordinate splitting along with natural intermediate boundary conditions. Results are presented for oblique plane waves and compared with other procedures. Preliminary computations of acoustic diffraction are also considered.

Davis, Sanford↗

User's manual for three dimensional FDTD version D code for scattering from frequency-dependent dielectric and magnetic materials

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code version D is a 3-D numerical electromagnetic scattering code based upon the finite difference time domain technique (FDTD). The manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into 14 sections: introduction; description of the FDTD method; operation; resource requirements; version D code capabilities; a brief description of the default scattering geometry; a brief description of each subroutine; a description of the include file; a section briefly discussing Radar Cross Section computations; a section discussing some scattering results; a sample problem setup section; a new problem checklist; references and figure titles. The FDTD technique models transient electromagnetic scattering and interactions with objects of arbitrary shape and/or material composition. In the FDTD method, Maxwell's curl equations are discretized in time-space and all derivatives (temporal and spatial) are approximated by central differences.

Beggs, John H.↗

Contributions of the NASA Langley Research Center to the DARPA/AFRL/NASA/ Northrop Grumman Smart Wing Program

An overview of the contributions of the NASA Langley Research Center (LaRC) to the DARPA/AFRL/NASA/ Northrop Grumman Corporation (NGC) Smart Wing program is presented. The overall objective of the Smart Wing program was to develop smart** technologies and demonstrate near-flight-scale actuation systems to improve the aerodynamic performance of military aircraft. NASA LaRC s roles were to provide technical guidance, wind-tunnel testing time and support, and Computational Fluid Dynamics (CFD) analyses. The program was divided into two phases, with each phase having two wind-tunnel entries in the Langley Transonic Dynamics Tunnel (TDT). This paper focuses on the fourth and final wind-tunnel test: Phase 2, Test 2. During this test, a model based on the NGC Unmanned Combat Air Vehicle (UCAV) concept was tested at Mach numbers up to 0.8 and dynamic pressures up to 150 psf to determine the aerodynamic performance benefits that could be achieved using hingeless, smoothly-contoured control surfaces actuated with smart materials technologies. The UCAV-based model was a 30% geometric scale, full-span, sting-mounted model with the smart control surfaces on the starboard wing and conventional, hinged control surfaces on the port wing. Two LaRC-developed instrumentation systems were used during the test to externally measure the shapes of the smart control surface and quantify the effects of aerodynamic loading on the deflections: Videogrammetric Model Deformation (VMD) and Projection Moire Interferometry (PMI). VMD is an optical technique that uses single-camera photogrammetric tracking of discrete targets to determine deflections at specific points. PMI provides spatially continuous measurements of model deformation by computationally analyzing images of a grid projected onto the model surface. Both the VMD and PMI measurements served well to validate the use of on-board (internal) rotary potentiometers to measure the smart control surface deflection angles. Prior to the final entry, NASA LaRC also performed three-dimensional unstructured Navier Stokes CFD analyses in an attempt to predict the potential aerodynamic impact of the smart control surface on overall model forces and moments. Eight different control surface shapes were selected for study at Mach = 0.6, Reynolds number = 3.25 x 10(exp 6), and + 2 deg., 3 deg., 8 deg., and 10 deg.model angles-of-attack. For the baseline, undeflected control surface geometry, the CFD predictions and wind-tunnel results matched well. The agreement was not as good for the more complex aero-loaded control surface shapes, though, because of the inability to accurately predict those shapes. Despite these results, the NASA CFD study served as an important step in studying advanced control effectors.

Florance, Jennifer P.↗

Stabilization of a Swept-Wing Boundary Layer by Discrete Roughness Elements at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to study potential stabilizing ef- fect of spanwise periodic discrete roughness elements (DREs) on cross ow instabilities in a spatially developing three-dimensional boundary layer over an in nite-swept natural- laminar- ow wing at a freestream Mach number of 0:75 and a chord Reynolds number of approximately 25 million. In the DNS, both the spanwise periodic DREs and distributed roughness in the leading-edge region are implemented to simulate a typical experimen- tal scenario in which multiple steady cross ow modes including the most unstable mode (i.e., the \target" mode) emerge because of the presence of naturally distributed surface roughness in the leading edge region and spanwise periodic control cylinders of subcritical wavelength are used to force small-wavelength disturbances (i.e., the control mode) for damping the target mode. The DNS results show that the e ectiveness of DRE control is sensitive to roughness diameter, height, and chordwise placement. For the DRE parame- ters considered in this study, the stabilizing e ect on the target mode is small within the computational domain that ended at about 35% of the chord.

Nicholson, Gary↗

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗