Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 397 records · Page 22

A Three-Dimensional Total Odd Nitrogen (NO(y)) Simulation During SONEX using a Stretched-Grid Chemical Transport Model

The relative importance of various odd nitrogen (NOy) sources including lightning, aircraft, and surface emissions on upper tropospheric total odd nitrogen is illustrated as a first application of the three-dimensional Stretched-Grid University of Maryland/Goddard Chemical-Transport Model (SG-GCTM). The SG-GCTM has been developed to look at the effect of localized sources and/or small scale mixing processes on the large-scale or global chemical balance. For this simulation, the stretched-arid was chosen so that its maximum resolution is located over eastern North America and the North Atlantic; a region that includes most of the SONEX (the SASS (Subsonic Assessment) Ozone and Nitrogen Oxides Experiment) flight paths. The SONEX period (October-November 1997) is simulated by driving the SG-GCTM with assimilated data from the GEOS-STRAT DAS (Goddard Earth Observing System-STRAT Data Assimilation System). A new algorithm is used to parameterize the lightning, flash rates that are needed to calculate emissions of NOy by lightning. Model-calculated upper tropospheric NOy and NOy measurements from the NASA DC-8 aircraft are compared. Spatial variations in NOy were well captured especially with the stretched-grid run; however, model-calculated concentrations were often too high in the upper troposphere, particularly during the first several flights. The lightning algorithm does a reasonably good job; however, the use of emissions from observed lightning, flashes significantly improves the simulation on a few occasions, especially November 3, 1997, indicating that significant uncertainty remains in parameterizing lightning in CTMS. Aircraft emissions play a relatively minor role (about 12%) in the upper tropospheric NOY budget averaged along SONEX flight paths; however, the contribution of such emmissions is as large as about 30% during portions of some flights.

Allen, Dale↗

Comparison of soil dielectric mixing models for Soil Moisture Retrieval using SMAP Brightness Temperature over croplands in India

The accurate estimation of soil moisture (SM) using microwave remote sensing depends mostly on careful selection of retrieval parameters among which the soil dielectric mixing model is the important one. These models are often categorized into empirical, semi-empirical or volumetric based on their methodologies and input data requirements. To study in detail, the comparative performance of four dielectric mixing models -- Wang & Schmugge model, Hallikainen model, Dobson model and Mironov model were used with Soil Moisture Active Passive (SMAP) L-band brightness temperature and Single Channel Algorithm for SM retrieval over agricultural landscapes in India. The highest performance statistics combination in terms of Root Mean Square Error (RMSE), correlation coefficient (R^(2)) and percentage bias (PBIAS) against the concurrent in-situ SM measurements were calculated at the selected validation sites. The overall results indicate that the best performance was given by the Mironov model (RMSE = 0.07 cu. m/cu. m), followed by Wang & Schmugge model (RMSE = 0.08 cu. m3/cu. m), Hallikainen model (RMSE = 0.09 cu. m/cu. m), Dobson model (RMSE = 0.10 cu. m/cu. m) and original SMAP radiometer SM (RMSE = 0.12 cu. m/cu. m). Findings of this study provides important insights into application and performance of dielectric mixing models in mapping surface SM variations. This study also underlines the pivotal role of local conditions for SM retrieval which should be carefully included in the algorithms.

Swati Suman↗

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗

Nimbus-7 coastal zone color scanner - System description and initial imagery

Initial imagery from the Nimbus-7 Coastal Zone Color Scanner (CZCS) shows subtle variations in water color. Organisms, especially phytoplankton, play a major influence on the variations in water color. Processing of the visual data is described, and in particular, the use of an algorithm to remove aerosol from the image is discussed. Data on the six spectral bands (433-12,500 nm) are presented and comparisons are made between the CZCS and the Landsat-1 multispectral scanner. The implications for management of fisheries is noted.

Hovis, W. A.↗

Evaluation of Nimbus 7 SMMR sensor with airborne radiometers and surface observations in the Norwegian Sea

Measurements made by the Nimbus 7 SMMR are compared with near simultaneous observations using the airborne SMMR simulator and with surface observations. The area of the test is in the Norwegian Sea between Bear Island and Northern Norway. It is noted that during the observation period two low-pressure systems were located in the test area, giving a spatial wind variation from 3-20 m/s. It is shown that the use of the currently available brightness temperatures and algorithms for SMMR does not give universally satisfactory results for SST and wind speed under extreme weather conditions. In addition, the SMMR simulator results are seen as indicating the need for more work on calibration.

Gloersen, P.↗

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗

Development of upwind schemes for the Euler equations

Described are many algorithmic and computational aspects of upwind schemes and their second-order accurate formulations based on Total-Variation-Diminishing (TVD) approaches. An operational unification of the underlying first-order scheme is first presented encompassing Godunov's, Roe's, Osher's, and Split-Flux methods. For higher order versions, the preprocessing and postprocessing approaches to constructing TVD discretizations are considered. TVD formulations can be used to construct relaxation methods for unfactored implicit upwind schemes, which in turn can be exploited to construct space-marching procedures for even the unsteady Euler equations. A major part of the report describes time- and space-marching procedures for solving the Euler equations in 2-D, 3-D, Cartesian, and curvilinear coordinates. Along with many illustrative examples, several results of efficient computations on 3-D supersonic flows with subsonic pockets are presented.

Chakravarthy, Sukumar R.↗

Artificial dissipation models for hypersonic external flow

Four artificial dissipation models which augment central difference schemes were examined for hypersonic external flows. The models were a first and third order dissipation model, a directionally scaled first and third order dissipation model, a flux limited dissipation model, and a flux difference split dissipation model. Each model was implemented in the lower-upper symmetric-Gauss-Seidel (LU-SGS) algorithm to solve the full Navier-Stokes equations. The latter two models can be regarded as total variation diminishing (TVD) schemes. Test results for model problems showed that the flux limited dissipation model was robust enough to predict a high speed blunt body flow with strong shock and expansion waves. The flux difference split dissipation model was capable of shock capturing with higher resolution, but was less robust. First and third order dissipation models turned out to be neither accurate nor robust enough for high Mach number flow computations.

Yoon, Seokkwan↗

Liquid crystal point diffraction interferometer

A new instrument, the liquid crystal point diffraction-interferometer (LCPDI), has been developed for the measurement of phase objects. This instrument maintains the compact, robust design of Linnik's point diffraction interferometer (PDI) and adds to it phase stepping capability for quantitative interferogram analysis. The result is a compact, simple to align, environmentally insensitive interferometer capable of accurately measuring optical wavefronts with very high data density and with automated data reduction. This dissertation describes the theory of both the PDI and liquid crystal phase control. The design considerations for the LCPDI are presented, including manufacturing considerations. The operation and performance of the LCPDI are discussed, including sections regarding alignment, calibration, and amplitude modulation effects. The LCPDI is then demonstrated using two phase objects: defocus difference wavefront, and a temperature distribution across a heated chamber filled with silicone oil. The measured results are compared to theoretical or independently measured results and show excellent agreement. A computer simulation of the LCPDI was performed to verify the source of observed periodic phase measurement error. The error stems from intensity variations caused by dye molecules rotating within the liquid crystal layer. Methods are discussed for reducing this error. Algorithms are presented which reduce this error; they are also useful for any phase-stepping interferometer that has unwanted intensity fluctuations, such as those caused by unregulated lasers.

Mercer, Carolyn R.↗

HEAO 3 observations of the Galactic center 511 keV line

The existence of an intense narrow electron/position annihilation line from the direction of the Galactic center is unquestioned; however, the source of the positrons remains a mystery. It is believed the line emisssion has two components: one compact at or near the Galactic center, and the other extended in the equatorial plane of the Galaxy. A compact component is suggested by several reports of intensity variations, including one from our group based on HEAO 3 data. Since the original HEAO 3 analysis, our algorithms have been greatly improved and we have gained a much better understanding of systematic effects which could influence the results. Given the importance of the result and the significantly improved analysis techniques now avaiable, the data have been reanalyzed. For a single point source at the Galactic center, the revised line intensities are (1.25 +/- 0.18) x 10(exp -3) photons cm(exp -2) s(exp -1) for the fall of 1979 and (0.99 +/- 0.18) x 10(exp 3) photons cm(exp -2) x(exp -1) for the spring of 1980.

Mahoney, W. A.↗

A Fully Implicit Time Accurate Method for Hypersonic Combustion: Application to Shock-induced Combustion Instability

A new fully implicit, time accurate algorithm suitable for chemically reacting, viscous flows in the transonic-to-hypersonic regime is described. The method is based on a class of Total Variation Diminishing (TVD) schemes and uses successive Gauss-Siedel relaxation sweeps. The inversion of large matrices is avoided by partitioning the system into reacting and nonreacting parts, but still maintaining a fully coupled interaction. As a result, the matrices that have to be inverted are of the same size as those obtained with the commonly used point implicit methods. In this paper we illustrate the applicability of the new algorithm to hypervelocity unsteady combustion applications. We present a series of numerical simulations of the periodic combustion instabilities observed in ballistic-range experiments of blunt projectiles flying at subdetonative speeds through hydrogen-air mixtures. The computed frequencies of oscillation are in excellent agreement with experimental data.

Yungster, Shaye↗

The feasibility of using time-dependent photochemical calculations to infer radical species concentrations from solar occultation absorption measurements

In connection with fast chemical reactions, short-lived stratospheric species experience rapid concentration variations at sunset and sunrise. For solar occultation absorption measurements, these rapid concentration variations may introduce significant errors with respect to the inference of atmospheric abundances for some species due to asymmetrical concentration distributions. Most retrieval algorithms assume that concentration distributions are spherically symmetric. The effect of this assumption on the accuracy of inferred concentrations has been studied by Kerr et al. (1977). The present investigation considers the feasibility of using a time-dependent one-dimensional photochemical model to provide detailed information about the asymmetrical distribution for use in the retrieval procedure. As shown by Boughner et al. (1980), diurnal effects can be represented by an inhomogeneity factor. It is found that the NO retrieval improves considerably with the inclusion of a correction factor containing the asymmetrical variations.

Larsen, J. C.↗

Airbreathing Propulsion System Analysis Using Multithreaded Parallel Processing

In this paper, parallel processing is used to analyze the mixing, and combustion behavior of hypersonic flow. Preliminary work for a sonic transverse hydrogen jet injected from a slot into a Mach 4 airstream in a two-dimensional duct combustor has been completed [Moon and Chung, 1996]. Our aim is to extend this work to three-dimensional domain using multithreaded domain decomposition parallel processing based on the flowfield-dependent variation theory. Numerical simulations of chemically reacting flows are difficult because of the strong interactions between the turbulent hydrodynamic and chemical processes. The algorithm must provide an accurate representation of the flowfield, since unphysical flowfield calculations will lead to the faulty loss or creation of species mass fraction, or even premature ignition, which in turn alters the flowfield information. Another difficulty arises from the disparity in time scales between the flowfield and chemical reactions, which may require the use of finite rate chemistry. The situations are more complex when there is a disparity in length scales involved in turbulence. In order to cope with these complicated physical phenomena, it is our plan to utilize the flowfield-dependent variation theory mentioned above, facilitated by large eddy simulation. Undoubtedly, the proposed computation requires the most sophisticated computational strategies. The multithreaded domain decomposition parallel processing will be necessary in order to reduce both computational time and storage. Without special treatments involved in computer engineering, our attempt to analyze the airbreathing combustion appears to be difficult, if not impossible.

Schunk, Richard Gregory↗

Latent Heating Retrievals Using the TRMM Precipitation Radar: A Multi-Seasonal Study

The Goddard Convective-Stratiform Heating (CSH) algorithm is used to retrieve profiles of latent heating over the global tropics for a period of several months using TRMM precipitation radar data. The seasonal variation of heating over the tropics is then examined. The period of interest also coincides with several TRMM field campaigns that recently occurred over the South China Sea in 1998 (SCSMEX), Brazil in 1999 (TRMM-LBA), and in the central Pacific in 1999 (KWAJEX). Sounding diagnosed Q1 budgets from these experiments could provide a means of validating the retrieved profiles of latent heating from the CSH algorithm.

Tao, W.-K.↗

Optimum propellant usage for reaction jet systems of space vehicles

The on-off type control for reaction jet systems is proven to be the optimal fuel scheme. However, due to the nonlinear characteristics of this type control, no direct method of solution is known for this optimal process being applied to the attitude control of space vehicles. This paper will discuss the optimization of fuel usage for an attitude control of space vehicles. A computation technique is developed for the calculation of optimal control law. The method of calculus of variations is applied to the estimation of the changes of performance index as well as terminal constraints. Thus an algorithm is obtained by the steepest descent method. A numerical example is given in the paper.

Liu, T. C.↗

Thoughts on standardization of parameters for image evaluation

Images received for image processing and analysis are obtained from a wide variety of sources and with a wide variety of sensors. Because it is desirable to have image processing algorithms be as universally applicable as possible, they should be designed, where possible, to be insensitive to the parametric variations of the source material. Where this is not possible, these variations must be taken into account. Therefore, it is necessary to consider what parameters may be defined in common across a suite of image types. Objective parameters or measurements of images which, in the proper combinations, may serve as surrogates for real images may be pixel-specific, location dependent, or combinations thereof. Parameters which have proven useful in defining the characteristics of images include the gray scale linearity, granularity of the quantization, spectral content, geometrical fidelity, resolution of the system expressed as either the point spread function or the modulation transfer function, and the spatial frequency content and characteristics of the data itself.

Billingsley, F. C.↗

Determination of sea ice parameters with the Nimbus 7 SMMR

A method of determining sea ice parameters using Nimbus 7 polarized multispectral radiance data obtained with the Nimbus 7 Scanning Multichannel Microwave Radiometer (SMMR) is presented. Observed radiances from selected areas in the Arctic region for the period February 3-7, 1979 were used in computing algorithm coefficients. Polar maps of sea ice concentration, multiyear fraction, and ice temperature are illustrated for this period. The variation of the mean and standard deviation of ice concentration and multiyear ice fraction for a region of perennial ice cover over the first 11 months of SMMR operation is also presented. Comparisons are made between the calculated sea ice parameters and information obtained from previous studies using aircraft, submarine and surface observations. The absolute accuracy of the SMMR parameters remains uncertain.

Cavalieri, D. J.↗