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 829 records · Page 46

High-resolution shock-capturing schemes for inviscid and viscous hypersonic flows

A class of implicit Total Variation Diminishing (TVD) type algorithms suitable for transonic and supersonic multidimensional Euler and Navier-Stokes equations was extended to hypersonic computations. The improved conservative shock-capturing schemes are spatially second- and third-order, and are fully implicit. They can be first- or second-order accurate in time and are suitable for either steady or unsteady calculations. Enhancement of stability and convergence rate for hypersonic flows is discussed. With the proper choice of the temporal discretization and suitable implicit linearization, these schemes are fairly efficient and accurate for very complex two-dimensional hypersonic inviscid and viscous shock interactions. This study is complimented by a variety of steady and unsteady viscous and inviscid hypersonic blunt-body flow computations. Due to the inherent stiffness of viscous flow problems, numerical experiments indicated that the convergence rate is in general slower for viscous flows than for inviscid steady flows.

Yee, H. C.↗

Coherence Analysis of the Space Launch System using Unsteady Pressure Sensitive Paint

Transonic buffet forces are a major source of unsteady loading on launch vehicles, thus requiring accurate estimation for efficient vehicle design. The state of the art in modeling these unsteady loads utilizes wind-tunnel tests where the fluctuating pressures are measured by pressure transducers (PTs) at discrete locations on a rigid buffet model. These pressures are then integrated over the surface of the vehicle to yield a series of orthogonal centerline loads called buffet forcing functions (BFFs). Typically, the PT layout aims at resolving the pressure correlation along the longitudinal axis of the vehicle. As a result, the distribution of azimuthal correlation and its impact on the estimated BFFs are not well known. To fill these gaps, extremely high-spatial-resolution uPSP data were collected for three different configurations of the Space Launch System in the NASA Ames Research Center 11-Foot Transonic Unitary Plan Wind Tunnel. The spatio-temporal behavior of the pressure correlation on these vehicles is analyzed and flow features of interest are investigated. It is shown that terminal shocks interacting with turbulence are a source of increased azimuthal coherence, especially when the shock develops at a junction. Vortex shedding off the forward attachment hardware that connects the core stage to the solid rocket boosters (SRBs) is the most severe buffet environment on the vehicle. The associated fluctuating pressures are shown to be highly coherent as far as the vehicle tail and up to 40 degrees away from the boosters. For selected areas of the vehicle, factoring the azimuthal coherence into the attenuation of discrete-measurements-based BFFs results in under prediction relative to the BFFs obtained from full integration of the uPSP data.

buffet↗

Release of liquid water from the Space Shuttle

Groundbased and onboard video images of a sunlit Shuttle Orbiter water dump are interpreted as showing that the continuous 1-mm-diameter liquid stream quickly breaks up in near-vacuum to form ice/snow particles of two characteristic sizes. Discrete large droplets are most evident in the close-in photographs, and unresolved submicron 'fog' from recondensation of overexpanded evaporated water appears to dominate the ground-telescope photographs of the 2.5 km long optically detectable trail. The mean diameter of the smaller particles was estimated from the spatial distribution of visible radiance using a model of their energy balance, (small) surface roughening as they sublime, and Mie scattering of pre-dawn sunlight. The results are consistent with those from recent space-tank simulations.

Pike, C. P.↗

Vision based obstacle detection and grouping for helicopter guidance

Electro-optical sensors can be used to compute range to objects in the flight path of a helicopter. The computation is based on the optical flow/motion at different points in the image. The motion algorithms provide a sparse set of ranges to discrete features in the image sequence as a function of azimuth and elevation. For obstacle avoidance guidance and display purposes, these discrete set of ranges, varying from a few hundreds to several thousands, need to be grouped into sets which correspond to objects in the real world. This paper presents a new method for object segmentation based on clustering the sparse range information provided by motion algorithms together with the spatial relation provided by the static image. The range values are initially grouped into clusters based on depth. Subsequently, the clusters are modified by using the K-means algorithm in the inertial horizontal plane and the minimum spanning tree algorithms in the image plane. The object grouping allows interpolation within a group and enables the creation of dense range maps. Researchers in robotics have used densely scanned sequence of laser range images to build three-dimensional representation of the outside world. Thus, modeling techniques developed for dense range images can be extended to sparse range images. The paper presents object segmentation results for a sequence of flight images.

Sridhar, Banavar↗

Analysis of Vegetation Within A Semi-Arid Urban Environment Using High Spatial Resolution Airborne Thermal Infrared Remote Sensing Data

High spatial resolution (5 m) remote sensing data obtained using the airborne Thermal Infrared Multispectral Scanner (TIMS) sensor for daytime and nighttime have been used to measure thermal energy responses for 2 broad classes and 10 subclasses of vegetation typical of the Salt Lake City, Utah urban landscape. Polygons representing discrete areas corresponding to the 10 subclasses of vegetation types have been delineated from the remote sensing data and are used for analysis of upwelling thermal energy for day, night, and the change in response between day and night or flux, as measured by the TIMS. These data have been used to produce three-dimensional graphs of energy responses in W/ sq m for day, night, and flux, for each urban vegetation land cover as measured by each of the six channels of the TIMS sensor. Analysis of these graphs provides a unique perspective for both viewing and understanding thermal responses, as recorded by the TIMS, for selected vegetation types common to Salt Lake City. A descriptive interpretation is given for each of the day, night, and flux graphs along with an analysis of what the patterns mean in reference to the thermal properties of the vegetation types surveyed in this study. From analyses of these graphs, it is apparent that thermal responses for vegetation can be highly varied as a function of the biophysical properties of the vegetation itself, as well as other factors. Moreover, it is also seen where vegetation, particularly trees, has a significant influence on damping or mitigating the amount of thermal radiation upwelling into the atmosphere across the Salt Lake City urban landscape. Published by Elsevier Science Ltd.

Quattrochi, Dale A.↗

Gradient-based surface nuclear magnetic resonance for groundwater investigation

In medical magnetic resonance imaging, spatial localization (imaging) is based upon the application of controlled magnetic field gradients on top of the main magnetic field to spatially modulate the frequency and/or phase of the nuclear magnetic resonance (NMR) signal across the volume of investigation. In this work, we have applied similar physical principles to produce controlled magnetic field gradients during surface NMR-based groundwater investigations. In this approach, a gradient pulse of variable amplitude or duration is applied immediately after the excitation pulse to cause predictable phase encoding of the NMR signal as a function of depth. This approach is also applicable to emerging surface NMR detection methods that use a prepolarization field with fast nonadiabatic turn-off to generate detectable NMR signals from the shallow subsurface. In this case, the gradient pulse is applied after terminating the prepolarization field and provides a heretofore unavailable means of localizing the NMR response as a function of depth. The application of gradients can also be combined with tip-angle-based modulation to yield higher imaging resolution than can be achieved through either gradient- or tip-angle-based imaging alone. We implemented this new gradient-based capability into a surface NMR gradient generation accessory that is compatible with the GMR-Flex instrument and developed surface NMR-specific forward modeling and linear inverse models. We validated the accuracy of this novel gradient-based sNMR technology using computer simulations, experiments using a small pool filled with a discrete layer of bulk water, and field experiments at well-characterized groundwater test sites along Ebey Island, WA, and Larned, KS. The gradient-based sNMR imaging observations were compared with high-resolution direct push NMR results observed at these sites. The results of computer simulations and field experiments indicate improvements in both the detection (signal-to-noise ratio) and spatial resolution of shallow subsurface water content using gradient-based surface NMR, compared with traditional surface NMR imaging methods.

Geochemistry & Geophysics↗

A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

97 MATHEMATICS AND COMPUTING↗

Exploring the MeV sky with a combined coded mask and Compton telescope: the Galactic Explorer with a Coded aperture mask Compton telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

79 ASTRONOMY AND ASTROPHYSICS↗

Physical properties of the OMC-2 and OMC-3 cores from CS and C(18)O observations

We have investigated the properties of the OMC-2 and OMC-3 cores in the Orion giant molecular cloud using high spatial spectral resolution observations of several transitions of the (13)CO, C(18)O, C(S-32) and C(S-34) molecules taken with the SEST telescope. The OMC-2 core consists of one clump (22 solar mass) with a radius of 0.11 pc surrounded by a cluster of 11 discrete infrared sources. The H2 column density and volume density in the center of this clump are 2 x 10(exp 22)/sq cm and 9 x 10(exp 5)/cu cm respectively. From a comparison between physical parameters derived from C(18)O and C(S-32) observations we conclude that the molecular envelope around the core has been completely removed by these sources and that only the very dense gas is left. OMC-3 shows a more complex elongated structure in C(18)O and CS than OMC-2. The C(S-32) and C(S-34) maps show that the denser region can be separated into at least sub-cores of roughly equal sizes (radius approximately equals 0.13 pc), with n(H2) = 6 x 10(exp 5)/cu cm, and a mass of 10 solar mass (from C(S-32)). The very different masses obtained for the central core from C(18)O and C(S-32) (55 and 12 solar mass respectively) indicate that a massive envelope is still present around the very dense sub-cores. We report the first detection of several molecular outflows in OMC-3. The presence of an IRAS source and the first detection of these outflows confirm that star formation is going on in OMC-3. Based on the different physical properties of these regions compared with OMC-1, OMC-2 appears to be in an intermediate evolutionary stage between OMC-1 and OMC-3.

Castets, A.↗

Implementation of Hybrid V-Cycle Multilevel Methods for Mixed Finite Element Systems with Penalty

The goal of this paper is the implementation of hybrid V-cycle hierarchical multilevel methods for the indefinite discrete systems which arise when a mixed finite element approximation is used to solve elliptic boundary value problems. By introducing a penalty parameter, the perturbed indefinite system can be reduced to a symmetric positive definite system containing the small penalty parameter for the velocity unknown alone. We stabilize the hierarchical spatial decomposition approach proposed by Cai, Goldstein, and Pasciak for the reduced system. We demonstrate that the relative condition number of the preconditioner is bounded uniformly with respect to the penalty parameter, the number of levels and possible jumps of the coefficients as long as they occur only across the edges of the coarsest elements.

Lai, Chen-Yao G.↗

Higher-Order Theory for Functionally Graded Materials

Functionally graded materials (FGM's) are a new generation of engineered materials wherein the microstructural details are spatially varied through nonuniform distribution of the reinforcement phase(s). Engineers accomplish this by using reinforcements with different properties, sizes, and shapes, as well as by interchanging the roles of the reinforcement and matrix phases in a continuous manner (ref. 1). The result is a microstructure that produces continuously or discretely changing thermal and mechanical properties at the macroscopic or continuum scale. This new concept of engineering the material's microstructure marks the beginning of a revolution both in the materials science and mechanics of materials areas since it allows one, for the first time, to fully integrate the material and structural considerations into the final design of structural components. Functionally graded materials are ideal candidates for applications involving severe thermal gradients, ranging from thermal structures in advanced aircraft and aerospace engines to computer circuit boards. Owing to the many variables that control the design of functionally graded microstructures, full exploitation of the FGM's potential requires the development of appropriate modeling strategies for their response to combined thermomechanical loads. Previously, most computational strategies for the response of FGM's did not explicitly couple the material's heterogeneous microstructure with the structural global analysis. Rather, local effective or macroscopic properties at a given point within the FGM were first obtained through homogenization based on a chosen micromechanics scheme and then subsequently used in a global thermomechanical analysis.

Aboudi, J.↗

Exploring the MeV Sky With A Combined Coded Mask and Compton Telescope: the Galactic Explorer With A Coded Aperture Mask Compton Telescope (GECCO)

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept that builds upon the heritage of past and current missions improving the sensitivity and, very importantly, the angular resolution. This consists in combining a Compton telescope and a coded-mask telescope. We delineate the motivation for such a concept and we define the scientific goals for such a mission. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that bridges the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO can disentangle the discrete sources from the truly diffuse emission. Individual Galactic and extragalactic sources are detected. This also allows to understand the gamma-ray Galactic center excess and the Fermi Bubbles, and to trace the low-energy cosmic rays, and their propagation in the Galaxy. Nuclear and annihilation lines are spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. Such an instrument also detects explosive transient gamma-ray sources, which, in turn, enables identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. By looking at a poorly explored energy band it also allows discoveries of new astrophysical phenomena.

Elena Orlando↗

Multi-Scale Fractal Analysis of Image Texture and Pattern

Fractals embody important ideas of self-similarity, in which the spatial behavior or appearance of a system is largely independent of scale. Self-similarity is defined as a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. An ideal fractal (or monofractal) curve or surface has a constant dimension over all scales, although it may not be an integer value. This is in contrast to Euclidean or topological dimensions, where discrete one, two, and three dimensions describe curves, planes, and volumes. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution. However, most geographical phenomena are not strictly self-similar at all scales, but they can often be modeled by a stochastic fractal in which the scaling and self-similarity properties of the fractal have inexact patterns that can be described by statistics. Stochastic fractal sets relax the monofractal self-similarity assumption and measure many scales and resolutions in order to represent the varying form of a phenomenon as a function of local variables across space. In image interpretation, pattern is defined as the overall spatial form of related features, and the repetition of certain forms is a characteristic pattern found in many cultural objects and some natural features. Texture is the visual impression of coarseness or smoothness caused by the variability or uniformity of image tone or color. A potential use of fractals concerns the analysis of image texture. In these situations it is commonly observed that the degree of roughness or inexactness in an image or surface is a function of scale and not of experimental technique. The fractal dimension of remote sensing data could yield quantitative insight on the spatial complexity and information content contained within these data. A software package known as the Image Characterization and Modeling System (ICAMS) was used to explore how fractal dimension is related to surface texture and pattern. The ICAMS software was verified using simulated images of ideal fractal surfaces with specified dimensions. The fractal dimension for areas of homogeneous land cover in the vicinity of Huntsville, Alabama was measured to investigate the relationship between texture and resolution for different land covers.

Emerson, Charles W.↗

Three-state nematicity in the triangular lattice antiferromagnet Fe 1/3 NbS 2

Nematic order is the breaking of rotational symmetry in the presence of translational invariance. While originally defined in the context of liquid crystals, the concept of nematic order has arisen in crystalline matter with discrete rotational symmetry, most prominently in the tetragonal Fe-based superconductors where the parent state is four-fold symmetric. In this case the nematic director takes on only two directions, and the order parameter in such 'Ising-nematic' systems is a simple scalar. In this study, using a spatially resolved optical polarimetry technique, we show that a qualitatively distinct nematic state arises in the triangular lattice antiferromagnet Fe 1/3 NbS 2 . The crucial difference is that the nematic order on the triangular lattice is a Z 3 or three-state Potts-nematic order parameter. As a consequence, the anisotropy axes of response functions such as the resistivity tensor can be continuously reoriented by external perturbations. This discovery lays the groundwork for devices that exploit analogies with nematic liquid crystals.

36 MATERIALS SCIENCE↗

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗

The cosmic X-ray background

As the sky in the microwave band is dominated by a cosmic background, so too is the X-ray sky. In this presentation the observational situation regarding the extragalactic X-ray background is reviewed, emphasizing data obtained from HEAO-1 and the Einstein Observatory (HEAO-2). Spectral characteristics and spatial variations are described and discussed within the context of what is currently known about individual extragalactic sources. It is concluded that the bulk of the cosmic X-ray background is yet to be understood. Possibilities range from genuinely diffuse emission to an 'as yet' unknown large population of unresolved discrete sources. The role of cosmological effects is examined and could be important.

Boldt, Elihu↗

Observation of pressure ridges in SAR images of sea ice: Scattering theory and comparison with observations

Ridges and keels (hummocks and bummocks) in sea ice flows are important in sea ice research for both scientific and practical reasons. Sea ice movement and deformation is driven by internal and external stresses on the ice. Ridges and keels play important roles in both cases because they determine the external wind and current stresses via drag coefficients. For example, the drag coefficient over sea ice can vary by a factor of several depending on the fluid mechanical roughness length of the surface. This roughness length is thought to be strongly dependent on the ridge structures present. Thus, variations in ridge and keel structure can cause gradients in external stresses which must be balanced by internal stresses and possibly fracture of the ice. Ridging in sea ice is also a sign of fracture. In a practical sense, large ridges form the biggest impediment to surface travel over the ice or penetration through sea ice by ice-strengthened ships. Ridges also play an important role in the damage caused by sea ice to off-shore structures. Hence, observation and measurement of sea ice ridges is an important component of sea ice remote sensing. The research reported here builds on previous work, estimating the characteristics of ridges and leads in sea ice from SAR images. Our objective is to develop methods for quantitative measurement of sea ice ridges from SAR images. To make further progress, in particular, to estimate ridge height, a scattering model for ridges is needed. Our research approach for a ridge scattering model begins with a survey of the geometrical properties of ridges and a comparison with the characteristics of the surrounding ice. For this purpose we have used airborne optical laser (AOL) data collected during the 1987 Greenland Sea Experiment. These data were used to generate a spatial wavenumber spectrum for height variance for a typical ridge - the typical ridge is the average over 10 large ridges. Our first-order model radar scattering includes both the quasi-specular and Bragg resonant scatter mechanisms. This model is extended to include contributions from volume scatter and scatter from discrete objects. Geometrical characteristics from the AOL survey and model calculations imply that for radar wavelengths and observation geometries that are dominated by the quasi-specular scattering mechanism radar backscatter from a ridge is a measure of peak ridge height. We present scattering model results and compare them with ridges observed during the LEADEX experiment of March-April 1992 when both X, C, and L-band aircraft SAR and the ERS-1 satellite SAR observed a region in the Beaufort Sea near 86 deg N, 10 deg W. Data were also collected documenting ridge characteristics on the surface. The surface data are used to generate a SAR signature via the scattering model described above. The predicted SAR signatures compare well with the SAR observations.

Vesecky, J. F.↗