Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Comparison of the Tangent Linear Properties of Tracer Transport Schemes Applied to Geophysical Problems.

A number of geophysical applications require the use of the linearized version of the full model. One such example is in numerical weather prediction, where the tangent linear and adjoint versions of the atmospheric model are required for the 4DVAR inverse problem. The part of the model that represents the resolved scale processes of the atmosphere is known as the dynamical core. Advection, or transport, is performed by the dynamical core. It is a central process in many geophysical applications and is a process that often has a quasi-linear underlying behavior. However, over the decades since the advent of numerical modelling, significant effort has gone into developing many flavors of high-order, shape preserving, nonoscillatory, positive definite advection schemes. These schemes are excellent in terms of transporting the quantities of interest in the dynamical core, but they introduce nonlinearity through the use of nonlinear limiters. The linearity of the transport schemes used in Goddard Earth Observing System version 5 (GEOS-5), as well as a number of other schemes, is analyzed using a simple 1D setup. The linearized version of GEOS-5 is then tested using a linear third order scheme in the tangent linear version.

4DVAR↗
NASA NTRS1972

Decay of isotropic turbulence generated by a mechanically agitated grid.

Experimental study of weak isotropic turbulence, created by a mechanically agitated grid, has indicated that in the absence of large linear-momentum wakes the energy of turbulence relaxes very quickly into a stable self-preserving structure, which, depending on the initial Reynolds number of turbulence, decays at different constant inverse powers of time. Both the longitudinal correlation coefficients and the corresponding spectral distributions, except for the difference in the parametric constants, are of the same functional type as those found previously for a passive grid.

Ling, S. C.↗
NASA NTRS1976

Feature combinations and the divergence criterion

Classifying large quantities of multidimensional remotely sensed agricultural data requires efficient and effective classification techniques and the construction of certain transformations of a dimension reducing, information preserving nature. The construction of transformations that minimally degrade information (i.e., class separability) is described. Linear dimension reducing transformations for multivariate normal populations are presented. Information content is measured by divergence.

Decell, H. P., Jr.↗
NASA NTRS1978

Linear dimension reduction and Bayes classification

An explicit expression for a compression matrix T of smallest possible left dimension K consistent with preserving the n variate normal Bayes assignment of X to a given one of a finite number of populations and the K variate Bayes assignment of TX to that population was developed. The Bayes population assignment of X and TX were shown to be equivalent for a compression matrix T explicitly calculated as a function of the means and covariances of the given populations.

Decell, H. P., Jr.↗
NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRS2011

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗
NASA NTRS2006

Sulfur during the Transition from Anoxic to Oxic Atmospheres

The invention of oxygenic photosynthesis was likely accompanied by the introduction of large amounts of O2 and complementary reduced gases (chiefly CH4) into the atmosphere. To first approximation the venting of O2 and CH4 are stochiometrically linked. We therefore present a suite of numerical photochemical models that address the anoxic-oxic transition in an atmosphere driven by large linked inputs of biogenic 02 and CH4. We find in general that, in steady state, there are two solutions, one oxic and the other anoxic. The anoxic solution appears to be linearly stable. If volcanic SO2 fluxes are large, S disproportionates into oxidized (H2S04) and reduced (S8) exit channels. As elemental sulfur is insoluble it provides a means of preserving photochemical mass-independent fractionation (MIF). On the other hand, if the source of volcanic SO2 is smaller than today, all S can leave the atmosphere as S8. Under these conditions there would be no MIF signal. The oxic solution appears to be linearly unstable. In the oxic solutions S is invariably oxidized to sulfate, and the MIF signal would be absent. The transitional atmosphere is relatively unstable and is also the most photochemically active. Consequently it is the transitional atmosphere, not the oxic or anoxic atmospheres, that has the lowest CH4 levels and weakest greenhouse warming. As a practical matter we expect the transitional atmospheres to vary strongly in response to diurnal and seasonal biological forcing.

Zahnle, Kevin↗
NASA NTRS1999

The Reconstruction Problem Revisited

The role of reconstruction in avoiding oscillations in upwind schemes is reexamined, with the aim of providing simple, concise proofs. In one dimension, it is shown that if the reconstruction is any arbitrary function bounded by neighboring cell averages and increasing within a cell for increasing data, the resulting scheme is monotonicity preserving, even though the reconstructed function may have overshoots and undershoots at the cell edges and is in general not a monotone function. In the special case of linear reconstruction, it is shown that merely bounding the reconstruction between neighboring cell averages is sufficient to obtain a monotonicity preservinc,y scheme. In two dimensions, it is shown that some ID TVD limiters applied in each direction result in schemes that are not positivity preserving, i.e. do not give positive updates when the data are positive. A simple proof is given to show that if the reconstruction inside the cell is bounded by the neighboring cell averages (including corner neighbors), then the scheme is positivity preserving. A new limiter that enforces this condition but is not as dissipative as the Minmod limiter is also presented.

Suresh, Ambaby↗
NASA NTRS1991

Model order reduction applied to a hot-bench simulation of an aeroelastic wind-tunnel model

Simulations of an aeroelastically scaled wind-tunnel model were developed for hot-bench testing of a digital controller. The digital controller provided active flutter-suppression, rolling-maneuver-load alleviation, and plant estimation. To achieve an acceptable time scale for the hot-bench application, the mathematical model of the wind-tunnel model was reduced from 220 states to approximately 130 states while assuring that the required accuracy was preserved for all combinations of 10 inputs and 56 outputs. The reduction was achieved by focussing on a linear, aeroelastic submodel of the full mathematical model and by applying a method based on the internally balanced realization of a dynamic system. The error-bound properties of the internally balanced realization significantly contribute to its utility in the model reduction process. The reduction method and the results achieved are described.

Buttrill, Carey S.↗
NASA NTRS1992

Coherent-backscatter effect - A vector formulation accounting for polarization and absorption effects and small or large scatterers

Previous theoretical work on the coherent-backscatter effect in the context of speckle time autocorrelation has gone beyond the diffusion approximation and the assumption of isotropic (point) scatterers. This paper extends the theory to include the effects of polarization and absorption, and to give the angular line shape. The results are expressions for angular variations valid for small and large scatterers and linear and circular polarizations, in lossless or lossy media. Calculations show that multiple anisotropic scattering results in the preservation of incident polarization. Application to a problem in radar astronomy is considered. It is shown that the unusual radar measurements (high reflectivity and polarization ratios) of Jupiter's icy Galilean satellites can be explained by coherent backscatter from anisotropic (forward) scatterers.

Peters, Kenneth J.↗
NASA NTRS1997

Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology

The compact form of the discontinuous Galerkin method allows for a detailed local analysis of the method in the neighborhood of the shock for a non-linear model problem. Insight gained from the analysis leads to new flux formulas that are stable and that preserve the compactness of the method. Although developed for a model equation, the flux formulas are applicable to systems such as the Euler equations. This article presents the analysis for methods with a degree up to 5. The analysis is accompanied by supporting numerical experiments using Burgers' equation and the Euler equations.

Atkins, H. L.↗
NASA NTRS1974

Physical conditions for conducting radiobiological experiments in beams of accelerated particles with high linear energy transfer

The design and construction of an accelerator to conduct radiobiological experiments is reported that uses aluminum filters to control the accelerated ion beam while preserving its stability, and a vacuum chamber to conduct the ion beam with the help of a collector through a lavsan exit port to the target. Depth distribution of the absorbed dose from a monodirectional ion beam is practically completely represented by the change in the energy spectrum of the biological object.

Kudryashov, Y. I.↗
NASA NTRS1991

Applications of singular value analysis and partial-step algorithm for nonlinear orbit determination

An adaptive method in which cruise and nonlinear orbit determination problems can be solved using a single program is presented. It involves singular value decomposition augmented with an extended partial step algorithm. The extended partial step algorithm constrains the size of the correction to the spacecraft state and other solve-for parameters. The correction is controlled by an a priori covariance and a user-supplied bounds parameter. The extended partial step method is an extension of the update portion of the singular value decomposition algorithm. It thus preserves the numerical stability of the singular value decomposition method, while extending the region over which it converges. In linear cases, this method reduces to the singular value decomposition algorithm with the full rank solution. Two examples are presented to illustrate the method's utility.

Ryne, Mark S.↗
NASA NTRS1997

An Investigation of High-Order Shock-Capturing Methods for Computational Aeroacoustics

This project is motivated by the desire to develop numerical methods that will be useful in the study of compressible flows that exhibit aeroacoustic phenomena. Solutions to linear problems have been investigated through the development of a computer code based on the recent dispersion-relation-preserving (DRP) methodology. In regard to nonlinear problems, the class of essentially nonoscillatory (ENO) schemes have been considered as the primary candidates for solving aeroacoustic problems in which discontinuities are involved. Discontinuities in the solution itself (e.g. shocks) as well as in the geometry on which the problem is defined have been studied. Two-dimensional nonlinear problems were considered in order to determine if the one-dimensional results obtained in the first phase of this project were extendable to more realistic problems. Conclusions have been drawn in regard to the ability to numerically predict solutions of nonlinear problems with shocks to high-order accuracy.

Casper, Jay↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
NASA NTRS1991

New bounds on parameter uncertainties for robust stability

The problem of robust stability in linear systems with parametric uncertainties is considered. A method is developed to determine bounds on uncertain parameters in the state-space model of the system, so that stability is preserved. Both continuous- and discrete-time systems are considered. Unlike previous results, the stability bounds are derived in terms of actual parameters instead of their absolute values, and consequently the stable region in the parameter space is not necessarily symmetric with respect to the origin. Furthermore, the uncertain parameters do not have to be linearly independent to each other.

Gao, Zhiqiang↗
NASA NTRS1989

A nearly-linear computational-cost scheme for the forward dynamics of an N-body pendulum

The dynamic equations of motion of an n-body pendulum with spherical joints are derived to be a mixed system of differential and algebraic equations (DAE's). The DAE's are kept in implicit form to save arithmetic and preserve the sparsity of the system and are solved by the robust implicit integration method. At each solution point, the predicted solution is corrected to its exact solution within given tolerance using Newton's iterative method. For each iteration, a linear system of the form J delta X = E has to be solved. The computational cost for solving this linear system directly by LU factorization is O(n exp 3), and it can be reduced significantly by exploring the structure of J. It is shown that by recognizing the recursive patterns and exploiting the sparsity of the system the multiplicative and additive computational costs for solving J delta X = E are O(n) and O(n exp 2), respectively. The formulation and solution method for an n-body pendulum is presented. The computational cost is shown to be nearly linearly proportional to the number of bodies.

Chou, Jack C. K.↗