Search NASA⌕ Search

SEARCH · Search NASA

Results for “asymptotic preserving”

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 19 records

Consensus-Based Formation Control of a Class of Multi-Agent Systems

This paper presents a consensus-based formation control scheme for autonomous multi-agent systems represented by double integrator dynamics. Assuming that the information graph topology consists of an undirected connected graph, a leader-based consensus-type control law is presented and shown to provide asymptotic formation stability when subjected to piecewise constant formation velocity commands. It is also shown that global asymptotic stability is preserved in the presence of (0, infinity)- sector monotonic non-decreasing actuator nonlinearities.

Joshi, Suresh↗

On the Stability of Collocated Controllers in the Presence or Uncertain Nonlinearities and Other Perils

Robustness properties are investigated for two types of controllers for large flexible space structures, which use collocated sensors and actuators. The first type is an attitude controller which uses negative definite feedback of measured attitude and rate, while the second type is a damping enhancement controller which uses only velocity (rate) feedback. It is proved that collocated attitude controllers preserve closed loop global asymptotic stability when linear actuator/sensor dynamics satisfying certain phase conditions are present, or monotonic increasing nonlinearities are present. For velocity feedback controllers, the global asymptotic stability is proved under much weaker conditions. In particular, they have 90 phase margin and can tolerate nonlinearities belonging to the (0,infinity) sector in the actuator/sensor characteristics. The results significantly enhance the viability of both types of collocated controllers, especially when the available information about the large space structure (LSS) parameters is inadequate or inaccurate.

Joshi, S. M.↗

Robustness properties of collocated controllers for flexible spacecraft

Robustness properties are investigated for two types of controllers for large flexible space structures, which use collocated sensors and actuators. The first type is an attitude controller which uses negative definite feedback of measured attitude and rate, while the second type is a damping enhancement controller which uses only velocity (rate) feedback. It is proved that collocated attitude controllers preserve closed loop global asymptotic stability when linear actuator/sensor dynamics satisfying certain phase conditions are present, or monotonic increasing nonlinearities are present. For velocity feedback controllers, the global asymptotic stability is proved under much weaker conditions. In particular, they have 90 phase margin and can tolerate nonlinearities belonging to the (0, infinity) sector in the actuator/sensor characteristics. The results significantly enhance the viability of both types of collocated controllers, especially when the available information about the large space structure (LSS) parameters is inadequate or inaccurate.

Joshi, S. M.↗

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup -alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state toward which a complete self-preserving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

Tracing Presolar Grain Populations Through Matrix Alteration Histories

Primitive meteorites preserve remnants of evolved stars as presolar grains, from asymptotic and red giant branch stars, novae and supernovae. These grains retain records of the formation mechanisms of the stellar environment from which they condensed, as well as physical and chemical processes in the interstellar medium (ISM), solar nebula and the parent bodies in which they are found. Consequently, their study can refine our understanding of the solar system and external stellar and planetary systems. However, presolar grains and their records can be altered, re-equilibrated with the surrounding matrix, and eventually lost through aqueous and thermal processing on the meteorite parent body. NanoSIMS is the most rapid and effective method for in situ identification of presolar grains via their highly anomalous isotopic ratios. Previous studies have measured variations in presolar grain populations between different matrix regions within the same meteorite thin section. However, most of these studies do not identify if there are any ties between these variations in presolar grain populations and the chemical and mineralogical differences in measured matrix regions. Therefore, the goal of this study is to determine if there is measurable variations in the degree of alteration of different matrix regions within the same thin section, and if mineralogical and chemical changes in matrix regions can be used as indicators for presolar grain preservation. Such a correlation could improve the efficiency of presolar grain identification during searches across large surface areas. To do this, we carried out coordinated scanning electron microscope (SEM) and NanoSIMS analyses, using a new SEM method which measures the ppm concentration of 100 nm to 5 μm sized Mg-, S-, Al-, C- and Ca-rich matrix grains.

presolar grains↗

A wave envelope finite element scheme for acoustical radiation

The aeroacoustic problem associated with the radiation of fan noise from the inlet of a turbofan aircraft engine, the dimensions of which are generally many times larger than the acoustical wavelengths of the major energy-carrying frequencies, is considered. In the present approach, a conventional finite element solution in the inner region is compatibly matched to a 'wave envelope' finite element solution in a large but finite outer region. The inclusion of a wavelike variation with the correct asymptotic decay in the shape functions for the outer region preserves the correct behavior of the solution at large distances. The method is initially presented for a simple one-dimensional model based on the solution of Webster's horn equation. Results are presented for the specific case of a uniform cylindrical section joined to a conical expansion, and also for a simple axisymmetric test case of the calculation of acoustical pressure generated by a vibrating circular piston located at the center of an infinite rigid baffle.

Astley, R. J.↗

Near-field turbulence properties of single- and two-stream plane mixing layers

Detailed mean-flow and turbulence measurements have been made in the near-field of two plane mixing layers in air with a maximum velocity of 21 m/s. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Cases with both tripped and untripped initial boundary layers were studied. In all cases, it was found that the two-stream layer developed to the self-preserving state in a distance much shorter than the single-stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two-stream layer, and the development of the single-stream layer showed agreement with existing data. The results suggest that the two-stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single-stream mixing layer.

Mehta, R. D.↗

An experimental study of plane mixing layer development

Mean flow and turbulence was measured in the near field of two plane mixing layers in air with a maximum velocity of 21 m/sec. The experimental rig enabled mixing layers of velocity ratios 0 and 0.46 to be generated simultaneously. Gases with both tripped and untripped initial boundary layers were studied. It was found that the two stream layer developed to the self preserving state in a distance much shorter than did the single stream layer which followed accepted criteria for the development distance. The asymptotic levels of the turbulence quantities in the two stream layer and the development of the single stream layer showed agreement with existing data. It is suggested that the two stream mixing layer should provide a better test case for the development of turbulence models and calculation methods than the single stream mixing layer.

Mehta, R. D.↗

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

A small-scale turbulence model

A model for the small-scale structure of turbulence is reformulated in such a way that it may be conveniently computed. The model is an ensemble of randomly oriented structured two dimensional vortices stretched by an axially symmetric strain flow. The energy spectrum of the resulting flow may be expressed as a time integral involving only the enstrophy spectrum of the time evolving two-dimensional cross section flow, which may be obtained numerically. Examples are given in which a k(exp -5/3) spectrum is obtained by this method without using large wave number asymptotic analysis. The k(exp -5/3) inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the two-dimensional enstrophy spectrum. The results are insensitive to time dependence of the strain-rate, including even intermittent on-or-off strains.

Lundgren, T. S.↗

Satellite rainfall retrieval by logistic regression

The potential use of logistic regression in rainfall estimation from satellite measurements is investigated. Satellite measurements provide covariate information in terms of radiances from different remote sensors.The logistic regression technique can effectively accommodate many covariates and test their significance in the estimation. The outcome from the logistical model is the probability that the rainrate of a satellite pixel is above a certain threshold. By varying the thresholds, a rainrate histogram can be obtained, from which the mean and the variant can be estimated. A logistical model is developed and applied to rainfall data collected during GATE, using as covariates the fractional rain area and a radiance measurement which is deduced from a microwave temperature-rainrate relation. It is demonstrated that the fractional rain area is an important covariate in the model, consistent with the use of the so-called Area Time Integral in estimating total rain volume in other studies. To calibrate the logistical model, simulated rain fields generated by rainfield models with prescribed parameters are needed. A stringent test of the logistical model is its ability to recover the prescribed parameters of simulated rain fields. A rain field simulation model which preserves the fractional rain area and lognormality of rainrates as found in GATE is developed. A stochastic regression model of branching and immigration whose solutions are lognormally distributed in some asymptotic limits has also been developed.

Chiu, Long S.↗

The energy decay in self-preserving isotropic turbulence revisited

The assumption of self-preservation allows for an analytical determination of the energy decay in isotropic turbulence. Here, the self-preserving isotropic decay problem is analyzed, yielding a more complete picture of self-serving isotropic turbulence. It is proven rigorously that complete self-serving isotropic turbulence admits two general types of asymptotic solutions: one where the turbulent kinetic energy K approximately t (exp -1) and one where K approximately t (sup alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one-point equations, it is demonstrated that the K approximately t (exp -1) and where K approximately t (sup -alpha) with an exponent alpha greater than 1 that is determined explicitly by the initial conditions. By a fixed point analysis and numerical integration of the exact one point equations, it is demonstrated that the K approximately t (exp -1) power law decay is the asymptotically consistent high Reynolds number solution; the K approximately 1 (sup - alpha) decay law is only achieved in the limit as t yields infinity and the turbulence Reynolds number vanishes. Arguments are provided which indicate that a K approximately t (exp -1) power law decay is the asymptotic state towards which a complete self-preseving isotropic turbulence is driven at high Reynolds numbers in order to resolve the imbalance between vortex stretching and viscous diffusion.

Speziale, Charles G.↗

Self-preservation of turbulent wakes

The present experiment has ascertained the development of the wake flow behind NASA GA(W)-1 airfoils, showing that, in the far wake, the mean velocity profiles exhibit self-similar behavior irrespective of the upstream boundary layer's character. It is noted, however, that the processes by means of which different wakes reach the asymptotic stage can be very different for different types of wake generators.

Mehta, Jayesh M.↗

Design of adaptive control systems by means of self-adjusting transversal filters

The design of closed-loop adaptive control systems based on nonparametric identification was addressed. Implementation is by self-adjusting Least Mean Square (LMS) transversal filters. The design concept is Model Reference Adaptive Control (MRAC). Major issues are to preserve the linearity of the error equations of each LMS filter, and to prevent estimation bias that is due to process or measurement noise, thus providing necessary conditions for the convergence and stability of the control system. The controlled element is assumed to be asymptotically stable and minimum phase. Because of the nonparametric Finite Impulse Response (FIR) estimates provided by the LMS filters, a-priori information on the plant model is needed only in broad terms. Following a survey of control system configurations and filter design considerations, system implementation is shown here in Single Input Single Output (SISO) format which is readily extendable to multivariable forms. In extensive computer simulation studies the controlled element is represented by a second-order system with widely varying damping, natural frequency, and relative degree.

Merhav, S. J.↗

Wavelength dependence in radio-wave scattering and specular-point theory

Radio-wave scattering from natural surfaces contains a strong quasispecular component that at fixed wavelengths is consistent with specular-point theory, but often has a strong wavelength dependence that is not predicted by physical optics calculations under the usual limitations of specular-point models. Wavelength dependence can be introduced by a physical approximation that preserves the specular-point assumptions with respect to the radii of curvature of a fictitious, effective scattering surface obtained by smoothing the actual surface. A uniform low-pass filter model of the scattering process yields explicit results for the effective surface roughness versus wavelength. Interpretation of experimental results from planetary surfaces indicates that the asymptotic surface height spectral densities fall at least as fast as an inverse cube of spatial frequency. Asymptotic spectral densities for Mars and portions of the lunar surface evidently decrease more rapidly.

Tyler, G. L.↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

Blocking reduction of Landsat Thematic Mapper JPEG browse images using optimal PSNR estimated spectra adaptive postfiltering

Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.

Linares, Irving↗

Solutions of the benchmark problems by the dispersion-relation-preserving scheme

The 7-point stencil Dispersion-Relation-Preserving scheme of Tam and Webb is used to solve all the six categories of the CAA benchmark problems. The purpose is to show that the scheme is capable of solving linear, as well as nonlinear aeroacoustics problems accurately. Nonlinearities, inevitably, lead to the generation of spurious short wave length numerical waves. Often, these spurious waves would overwhelm the entire numerical solution. In this work, the spurious waves are removed by the addition of artificial selective damping terms to the discretized equations. Category 3 problems are for testing radiation and outflow boundary conditions. In solving these problems, the radiation and outflow boundary conditions of Tam and Webb are used. These conditions are derived from the asymptotic solutions of the linearized Euler equations. Category 4 problems involved solid walls. Here, the wall boundary conditions for high-order schemes of Tam and Dong are employed. These conditions require the use of one ghost value per boundary point per physical boundary condition. In the second problem of this category, the governing equations, when written in cylindrical coordinates, are singular along the axis of the radial coordinate. The proper boundary conditions at the axis are derived by applying the limiting process of r approaches 0 to the governing equations. The Category 5 problem deals with the numerical noise issue. In the present approach, the time-independent mean flow solution is computed first. Once the residual drops to the machine noise level, the incident sound wave is turned on gradually. The solution is marched in time until a time-periodic state is reached. No exact solution is known for the Category 6 problem. Because of this, the problem is formulated in two totally different ways, first as a scattering problem then as a direct simulation problem. There is good agreement between the two numerical solutions. This offers confidence in the computed results. Both formulations are solved as initial value problems. As such, no Kutta condition is required at the trailing edge of the airfoil.

Tam, Christopher K. W.↗