SEARCH · Search NASA

Results for “Variational 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Scattering Properties and Brightness Temperatures Associated with Solid Precipitation

In the past few years, early solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. NOAA has an operational snow versus rain classifier based on AMSU-B observations. Solid precipitation retrieval algorithms reported in the literature over the past two years include those that rely on neural nets, statistics, or physical relationships. All of the algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This investigation probes the sensitivity of scattering and absorption coefficients, and hence computed brightness temperatures, resulting from variations in solid precipitation cloud profiles. The first study compares the single scattering, absorption, and asymmetry parameters associated with snow particles in clouds. Several methodologies are used to convert the physical characteristics (e.g., shape, size distributions, ice-air-water ratios) of ice particles to electromagnetic properties (e.g., absorption, scattering, and asymmetry factors). These methodologies include: conversion to solid ice particles, homogeneous dielectric mixing, or discrete dipole approximation. Changes in the conversion methodology can produce computed brightness temperature differences greater than 50 Kelvin.

Skofronick-Jackson, Gail M.↗

NASA NTRS1983

The attenuation of sound in a randomly lined duct

The kinetic theory of wave propagation in random media is applied to assess the relative merits of the peripherally segmented, axially segmented and 'checkerboard' configured acoustic duct treatments proposed by Mani (1976). It is argued that, in a first approximation, a multisegmented liner involving many discrete and distinct sections may be modelled as a liner whose impedance varies randomly with position. This results in a relatively simple analytical description of the acoustic properties of the liner in terms of a system of linear, first order ordinary differential equations, whose solution gives the distribution with axial distance of the acoustic modal intensities. The theory is expected to be valid provided that the variation in the impedance from segment to segment is small, so that the conclusions obtained here are of a preliminary nature, and must be regarded as giving a general indication of features which should be incorporated into the design of efficient acoustic liners. The principal such conclusion is that a checkerboard liner offers the best means of increasing the attenuation over and above that of a uniformly lined duct. On the other hand, only minimal gains are to be derived from the use of strictly peripherally segmented liners.

Howe, M. S.↗

NASA NTRS1989

Composite structure dynamics and space applicability

Guided wave modes in cladded or uncladded fiber-reinforced composite plates and tubes have been analyzed using a stiffness method in which the displacement variation through the thickness is approximated by polynomial interpolation functions. This allows for an arbitrary number of laminations and fiber orientations different from lamina to lamina. It is shown that dispersive behavior of guided modes depends significantly on the cladding, number of laminae, and interfaces between the adjacent laminae. A hybrid modeling technique is described in which an inner region containing cracks (or other defects) is discretized by finite elements, and the field in the exterior region is represented in terms of modes that are found using the stiffness method described above. It is found that the reflected and transmitted amplitudes of modes vary significantly with the size of a transverse or longitudinal (delamination) crack and frequency. We have also studied the impact response of a unidirectional fiber-reinforced plate. Received signals at the epicentral and other locations are shown. Strong longitudinal anisotropy of the graphite/epoxy plate causes the signal to be considerably different from that in an isotropic plate.

Datta, Subhendu K.↗

NASA NTRS2002

Efficient Reformulation of the Thermoelastic Higher-order Theory for Fgms

Functionally graded materials (FGMs) are characterized by spatially variable microstructures which are introduced to satisfy given performance requirements. The microstructural gradation gives rise to continuously or discretely changing material properties which complicate FGM analysis. Various techniques have been developed during the past several decades for analyzing traditional composites and many of these have been adapted for the analysis of FGMs. Most of the available techniques use the so-called uncoupled approach in order to analyze graded structures. These techniques ignore the effect of microstructural gradation by employing specific spatial material property variations that are either assumed or obtained by local homogenization. The higher-order theory for functionally graded materials (HOTFGM) is a coupled approach developed by Aboudi et al. (1999) which takes the effect of microstructural gradation into consideration and does not ignore the local-global interaction of the spatially variable inclusion phase(s). Despite its demonstrated utility, however, the original formulation of the higher-order theory is computationally intensive. Herein, an efficient reformulation of the original higher-order theory for two-dimensional elastic problems is developed and validated. The use of the local-global conductivity and local-global stiffness matrix approach is made in order to reduce the number of equations involved. In this approach, surface-averaged quantities are the primary variables which replace volume-averaged quantities employed in the original formulation. The reformulation decreases the size of the global conductivity and stiffness matrices by approximately sixty percent. Various thermal, mechanical, and combined thermomechanical problems are analyzed in order to validate the accuracy of the reformulated theory through comparison with analytical and finite-element solutions. The presented results illustrate the efficiency of the reformulation and its advantages in analyzing functionally graded materials.

Bansal, Yogesh↗

NASA NTRS1972

Theoretical study of corrugated plates: Shear stiffness of a trapezoidally corrugated plate with discrete attachments to a rigid flange at the ends of the corrugations

Analysis and numerical results are presented for the elastic shear stiffness of a corrugated shear web with a certain type of discrete attachments at the ends of the trough lines of the corrugations, namely point attachments to a rigid flange which interferes with the deformations of the end cross sections by preventing downward movement but permitting upward (lifting off) movement. The analysis is based on certain assumed modes of deformation of the cross sections in conjunction with the method of minimum total potential energy and the calculus of variations in order to obtain equations for the manner in which the assumed modes of deformation vary along the length of the corrugation. The numerical results are restricted to the case of equal-width crests and troughs but otherwise apply to a wide variety of geometries. They are in the form of graphs which give the overall shear stiffness as a fraction of the overall shear stiffness that could be obtained by having continuous attachment at the ends of the corrugations.

Hsiao, C.↗

NASA NTRS1972

Anomalous type III bursts explained in terms of intermittent solar wind streams observed by Pioneer spacecraft.

Occasional solar radio bursts drift downward in frequency (like type III bursts), but the drift slows or ceases or it may even reverse. These anomalous forms are called L-bursts, U-bursts, etc., and are observed as low as 1 MHz by spacecraft. The conventional explanation involves guidance of relativistic exciting particles by solar magnetic fields that do not have the spiral shape expected in undisturbed conditions. A reasonable alternative explanation is possible and, indeed, it follows directly from the inferences which have been drawn from the Stanford dual-frequency experiment flown on four Pioneer spacecraft. Analysis of temporal variations in these data leads to the conclusion that narrow streams of enhanced concentration have lifetimes which are typically about one day and the streams seem to be intermittently ejected by long-enduring sources at the sun. Each enhancement is thought to lie along a spiral and to have a beginning and an end. It may happen that relativistic exciters encounter the tail of such a discrete streamer; the concentration encountered by such exciters would initially decline but, at the streamer tail, it would level off and perhaps increase. The plasma frequency versus time (the burst frequency) would then cease to drop in the typical type III profile and would level off or increase (as in a U-burst).

Croft, T. A.↗

NASA NTRS1999

Some Characteristics of the Martian Aphelion Global Cloud Belt

The presence of discrete condensate clouds on Mars is certainly not a new discovery, having been observed through most of the documented history of telescopic monitoring. Furthermore, spacecraft data have been used to study discrete cloud features in the Martian atmosphere in greater detail, e.g., morphology, seasonal occurrence. Condensate clouds, specifically discrete water ice clouds, appeared to be regarded as fairly common but, with the possible exception of the polar regions, generally uninteresting from a climatological point of view. However, recent observations indicate that in addition to their large spatial scale, the water ice clouds may in fact play a more prominent role in the Martian climate. In this paper, we wish to examine the spatial and temporal variations of the cloud belt optical depth, as well as the microphysical characteristics of the water ice particles themselves.

Wolff, M. J.↗

NASA NTRS1993

Two-dimensional radiative transfer in cloudy atmospheres - The spherical harmonic spatial grid method

A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels of approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres.

Evans, K. F.↗

NASA NTRS2001

Fiber Grating Coupled Light Source Capable of Tunable, Single Frequency Operation

Fiber Bragg grating coupled light sources can achieve tunable single-frequency (single axial and lateral spatial mode) operation by correcting for a quadratic phase variation in the lateral dimension using an aperture stop. The output of a quasi-monochromatic light source such as a Fabry Perot laser diode is astigmatic. As a consequence of the astigmatism, coupling geometrics that accommodate the transverse numerical aperture of the laser are defocused in the lateral dimension, even for apsherical optics. The mismatch produces the quadratic phase variation in the feedback along the lateral axis at the facet of the laser that excites lateral modes of higher order than the TM 00 . Because the instability entails excitation of higher order lateral submodes, single frequency operation also is accomplished by using fiber Bragg gratings whose bandwidth is narrower than the submode spacing. This technique is particularly pertinent to the use of lensed fiber gratings in lieu of discrete coupling optics. Stable device operation requires overall phase match between the fed-back signal and the laser output. The fiber Bragg grating acts as a phase-preserving mirror when the Bragg condition is met precisely. The phase-match condition is maintained throughout the fiber tuning range by matching the Fabry-Perot axial mode wavelength to the passband center wavelength of the Bragg grating.

Michael A Krainak↗

NASA NTRS2010

Cryogenic Carbonate Formation on Mars: Clues from Stable Isotope Variations Seen in Experimental Studies

Discoveries of large deposits of sedimentary materials on the planet Mars by landers and orbiters have confirmed the widely held hypothesis that water has played a crucial role in the development of the martian surface. Recent studies have indicated that both water ice and liquid water may have been present and in the case of water ice perhaps is still present on or near the surface of Mars. However, there remains much controversy about the prevailing atmospheric conditions and climate of Mars during its history and whether liquid water existed on the martian surface simply during discrete geological events or whether this water was present over relatively much longer geologic time periods. The recent identification of Ca-rich carbonate by the Phoenix lander as well as its measurement of the isotopic composition of atmospheric CO2 has shown the importance of understanding the carbonates on Mars as an important sink of atmospheric carbon. This work compliments that of our past experiments where we produced cryogenic calcite in open containers, as analogs for terrestrial aufeis formation, and as a means for evaluating the fractionation of C-13 in CO2 during bicarbonate freezing [13]. Unlike our previous experiments in which carbonates were grown in ambient laboratory condition in open containers (atmospheric pressure and composition), this work attempts to quantify the amount of delta C-13 enrichment possible in both fluids and secondary carbonates formed from freezing of bicarbonate fluids under martian-like atmospheric conditions. Morphologic textures of produced carbonates in these experiments are also examined under SEM in order to identify the effect that the cryogenic freezing process has on the mineral's mineralogy. Understanding the role of kinetic isotope fractionation during formation of carbonates under martian-like conditions will aid in our ability to quantify the isotopic composition of the carbonate sink furthering our ability to model the climate history of Mars.

Socki, Richard A.↗

NASA NTRS1972

Analysis of photometric data obtained on the 1969 NASA airborne auroral expedition

The analysis of the 1969 airborne expedition data is centered around two main topics: dayside auroral characteristics and hydroxial airglow variations. Some of the conclusions drawn from the analysis include the following: (1) The midday sector of the auroral oval is bound on the equatorward side by a region of greater than or = to 3 keV electrons. (2) The dayside oval is associated with the cusp region in the magnetosphere which has particle energies usually less than or = to 1 keV and lambda 6300 (OI)/lambda 4278 N2(+) ratios 10. (3) Discrete earth-sun aligned auroras appear in the polar cap just poleward of the oval superimposed on a 1 keV particle energy background. (4) The discrete earth-sun aligned auroras seen in the open field line region on the dayside have a similar appearance to some rayed auroras seen on the nightside both in terms of charactertistic particle energy spectra as well as visual appearance. (5) Airglow OH measurements show no appreciable change with latitude or longitude but do show diurnal changes. Two types of enhancements are observed: one associated with visible aurora, and the other uncorrelated with visible aurora.

Romick, G. J.↗

NASA NTRS2011

Superpixel-Augmented Endmember Detection for Hyperspectral Images

Superpixels are homogeneous image regions comprised of several contiguous pixels. They are produced by shattering the image into contiguous, homogeneous regions that each cover between 20 and 100 image pixels. The segmentation aims for a many-to-one mapping from superpixels to image features; each image feature could contain several superpixels, but each superpixel occupies no more than one image feature. This conservative segmentation is relatively easy to automate in a robust fashion. Superpixel processing is related to the more general idea of improving hyperspectral analysis through spatial constraints, which can recognize subtle features at or below the level of noise by exploiting the fact that their spectral signatures are found in neighboring pixels. Recent work has explored spatial constraints for endmember extraction, showing significant advantages over techniques that ignore pixels relative positions. Methods such as AMEE (automated morphological endmember extraction) express spatial influence using fixed isometric relationships a local square window or Euclidean distance in pixel coordinates. In other words, two pixels covariances are based on their spatial proximity, but are independent of their absolute location in the scene. These isometric spatial constraints are most appropriate when spectral variation is smooth and constant over the image. Superpixels are simple to implement, efficient to compute, and are empirically effective. They can be used as a preprocessing step with any desired endmember extraction technique. Superpixels also have a solid theoretical basis in the hyperspectral linear mixing model, making them a principled approach for improving endmember extraction. Unlike existing approaches, superpixels can accommodate non-isometric covariance between image pixels (characteristic of discrete image features separated by step discontinuities). These kinds of image features are common in natural scenes. Analysts can substitute superpixels for image pixels during endmember analysis that leverages the spatial contiguity of scene features to enhance subtle spectral features. Superpixels define populations of image pixels that are independent samples from each image feature, permitting robust estimation of spectral properties, and reducing measurement noise in proportion to the area of the superpixel. This permits improved endmember extraction, and enables automated search for novel and constituent minerals in very noisy, hyperspatial images. This innovation begins with a graph-based segmentation based on the work of Felzenszwalb et al., but then expands their approach to the hyperspectral image domain with a Euclidean distance metric. Then, the mean spectrum of each segment is computed, and the resulting data cloud is used as input into sequential maximum angle convex cone (SMACC) endmember extraction.

Thompson, David R.↗

NASA NTRS1992

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗

NASA NTRS1983

Nearfield observations of tones generated from supersonic jet flows

The physics of discrete tone generation (screech) from improperly expanded supersonic jets have been investigated experimentally. Both nearfield acoustic survey and phase-averaged schlieren flow visualization were conducted. The dominant mode of flow oscillation that occurred during intense screech and its coupling with the internal shock structure and external sound field were established. Deductions made from the observed coupling mechanism revealed the importance of Helmholtz number and shock cell spacing in the maintenance of screech. Measured amplitude variation of screech are compared with computed amplitude variation of the dominant instability. The results suggest that jet instability plays a direct role in the generation of screech.

Yu, J. C.↗

NASA NTRS1989

Supersonic rectangular jet impingement noise experiments

The discrete frequency sound produced by jets issuing from a convergent, rectangular nozzle of aspect ratio 4.24 was investigated. Experiments were performed both with the free jet and with the jet impinging on a hard ground surface. The impingement tones that dominate the impinging jet spectra show a definite staging behavior which appears to be biased toward the free jet screech frequency once the separation distance exceeds the region of substantial shock cell development. The frequency variation of the impingement tone stages fit the details of a feedback cycle if the disturbance convection velocity is chosen to be 20 percent higher than that necessary to satisfy the screech feedback loop. Phase locked optical records show a flapping mode of jet oscillation with tones at or near the screech frequency, with superimposed symmetric oscillations when a second dominant tone of unrelated frequency appears in the spectrum.

Norum, Thomas D.↗

NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

NASA NTRS2012

Reynolds-Averaged Turbulence Model Assessment for a Highly Back-Pressured Isolator Flowfield

The use of computational fluid dynamics in scramjet engine component development is widespread in the existing literature. Unfortunately, the quantification of model-form uncertainties is rarely addressed with anything other than sensitivity studies, requiring that the computational results be intimately tied to and calibrated against existing test data. This practice must be replaced with a formal uncertainty quantification process for computational fluid dynamics to play an expanded role in the system design, development, and flight certification process. Due to ground test facility limitations, this expanded role is believed to be a requirement by some in the test and evaluation community if scramjet engines are to be given serious consideration as a viable propulsion device. An effort has been initiated at the NASA Langley Research Center to validate several turbulence closure models used for Reynolds-averaged simulations of scramjet isolator flows. The turbulence models considered were the Menter BSL, Menter SST, Wilcox 1998, Wilcox 2006, and the Gatski-Speziale explicit algebraic Reynolds stress models. The simulations were carried out using the VULCAN computational fluid dynamics package developed at the NASA Langley Research Center. A procedure to quantify the numerical errors was developed to account for discretization errors in the validation process. This procedure utilized the grid convergence index defined by Roache as a bounding estimate for the numerical error. The validation data was collected from a mechanically back-pressured constant area (1 2 inch) isolator model with an isolator entrance Mach number of 2.5. As expected, the model-form uncertainty was substantial for the shock-dominated, massively separated flowfield within the isolator as evidenced by a 6 duct height variation in shock train length depending on the turbulence model employed. Generally speaking, the turbulence models that did not include an explicit stress limiter more closely matched the measured surface pressures. This observation is somewhat surprising, given that stress-limiting models have generally been developed to better predict shock-separated flows. All of the models considered also failed to properly predict the shape and extent of the separated flow region caused by the shock boundary layer interactions. However, the best performing models were able to predict the isolator shock train length (an important metric for isolator operability margin) to within 1 isolator duct height.

Baurle, Robert A.↗

NASA NTRS1988

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.↗