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 1,027 records · Page 57

A Combined Infrared and Microwave Technique for Studying the Diurnal Variation of Rainfall over Amazonia

Over 10 years of continuous data from the Special Sensor microwave Imager (SSM/I) aboard a series of Defense Department satellites has made it possible to construct regional rainfall climatologies at high spatial resolution. Using the Goddard Profiling Algorithm (GPROF), monthly estimates of precipitation were made over the region of northern Brazil, including the Amazon Basin, for 1987 to 1998. GPROF is a physical approach to passive microwave precipitation retrieval, which uses the Goddard Cumulus Ensemble (cloud) model to establish prior probability densities of precipitation structures. Precipitation fields from GPROF were stratified into morning and evening satellite overpasses, and accumulated at monthly intervals at 0.5 degree spatial resolution. Important diurnal effects were noted in the analysis, the most pronounced being a land/sea breeze circulation along the northern coast of Brazil and a mountain/valley circulation along the Andes. There were also indications of morning rainfall maxima along the major rivers, and evening maxima between the rivers. The addition of simultaneous geosynchronous infrared (IR) data leads to the current technique, which takes advantage of the 30 minute sampling and 4 km spatial resolution of the infrared channel and the better physics of the microwave retrieval. The resultant IR method is subsequently used to derive the diurnal variability of rainfall over the Amazon basin, and further, to investigate the relative contribution from its convective and stratiform components.

Negri, Andrew J.↗

Estimation of Bare Surface Soil Moisture and Surface Roughness Parameter Using L-Band SAR Image Data

An algorithm based on a fit of the single-scattering Integral Equation Method (IEM) was developed to provide estimation of soil moisture and surface roughness parameter (a combination of rms roughness height and surface power spectrum) from quad-polarized synthetic aperture radar (SAR) measurements. This algorithm was applied to a series of measurements acquired at L-band (1.25 GHz) from both AIRSAR (Airborne Synthetic Aperture Radar operated by the Jet Propulsion Laboratory) and SIR-C (Spaceborne Imaging Radar-C) over a well- managed watershed in southwest Oklahoma. Prior to its application for soil moisture inversion, a good agreement was found between the single-scattering IEM simulations and the L band measurements of SIR-C and AIRSAR over a wide range of soil moisture and surface roughness conditions. The sensitivity of soil moisture variation to the co-polarized signals were then examined under the consideration of the calibration accuracy of various components of SAR measurements. It was found that the two co-polarized backscattering coefficients and their combinations would provide the best input to the algorithm for estimation of soil moisture and roughness parameter. Application of the inversion algorithm to the co-polarized measurements of both AIRSAR and SIR-C resulted in estimated values of soil moisture and roughness parameter for bare and short-vegetated fields that compared favorably with those sampled on the ground. The root-mean-square (rms) errors of the comparison were found to be 3.4% and 1.9 dB for soil moisture and surface roughness parameter, respectively.

Shi, Jian-Cheng↗

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

White light modeling, algorithm development, and validation on the micro-arcsecond metrology testbed

The Space Interferometry Mission (SIM), scheduled for launch in early 2010, is an optical interferometer that will perform narrow angle and global wide angle astrometry with unprecedented accuracy, providing differential position accuracies of 1uas, and 4uas global accuracies in position, proper motion and parallax. SIM astrometric measurements are synthesized from pathlength delay measurements provided by three Michelson-type, white light interferometers. Two of the interferometers are used for making precise measurements of variations in the spacecraft attitude, while the third interferometer performs the science measurement. The ultimate performance of SIM relies on a combination of precise fringe measurements of the interfered starlight with picometer class relative distance measurements made between a set of fiducials that define the interferometer baseline vectors. The focus of the present paper is on the development and analysis of algorithms for accurate white light fringe estimation, and on the preliminary validation of these algorithms on the Micro-Arcsecond Testbed (MAM).

Shen, Janice TsaePyng↗

Preliminary results in using Deep Learning to emulate BLOB, a nuclear interaction model

Purpose: A reliable model to simulate nuclear interactions is fundamental for Ion-therapy. We already showed how BLOB (“Boltzmann-Langevin One Body”), a model developed to simulate heavy ion interactions up to few hundreds of MeV/u, could simulate also 12 C reactions in the same energy domain. However, its computation time is too long for any medical application. For this reason we present the possibility of emulating it with a Deep Learning algorithm. Methods: The BLOB final state is a Probability Density Function (PDF) of finding a nucleon in a position of the phase space. We discretised this PDF and trained a Variational Auto-Encoder (VAE) to reproduce such a discrete PDF. As a proof of concept, we developed and trained a VAE to emulate BLOB in simulating the interactions of 12 C with 12 C at 62 MeV/u. To have more control on the generation, we forced the VAE latent space to be organised with respect to the impact parameter (b) training a classifier of b jointly with the VAE. Results: In this work, the distributions obtained from the VAE are similar to the input ones and the computation time needed to use the VAE as a generator is negligible. Conclusions: We show that it is possible to use a Deep Learning approach to emulate a model developed to simulate nuclear reactions in the energy range of interest for Ion-therapy. We foresee the implementation of the generation part in C++ and to interface it with the most used Monte Carlo toolkit: Geant4.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantitative interpretation of Great Lakes remote sensing data

The paper discusses the quantitative interpretation of Great Lakes remote sensing water quality data. Remote sensing using color information must take into account (1) the existence of many different organic and inorganic species throughout the Great Lakes, (2) the occurrence of a mixture of species in most locations, and (3) spatial variations in types and concentration of species. The radiative transfer model provides a potential method for an orderly analysis of remote sensing data and a physical basis for developing quantitative algorithms. Predictions and field measurements of volume reflectances are presented which show the advantage of using a radiative transfer model. Spectral absorptance and backscattering coefficients for two inorganic sediments are reported.

Shook, D. F.↗

A preliminary study of relaxation methods for the inviscid conservative gasdynamics equations using flux splitting

Plus minus flux vector split schemes are combined with the multigrid relaxation method to obtain fast rates of iterative (or steady state) convergence. On model nonlifting biconvex airfoil flow, steady state convergence is obtained using multigrid in 40 iteractions for subsonic cases and in 120 iterations for transonic cases. A 65 x 33 stretched grid is used. However, the multigrid scheme used here (and there are numerous variations) sometimes could not be made to work for transonic cases with a large number of supersonic points. Moreover, the multigrid method is only several times faster than the conventional algorithm which itself could perhaps be accelerated by other means. Nethertheless, with proper adjustment the multigrid method will ultimately lead to very fast rates of iterative convergence.

Steger, J. L.↗

Land use in the Paraiba Valley through remotely sensed data

A methodology for land use survey was developed and land use modification rates were determined using LANDSAT imagery of the Paraiba Valley (state of Sao Paulo). Both visual and automatic interpretation methods were employed to analyze seven land use classes: urban area, industrial area, bare soil, cultivated area, pastureland, reforestation and natural vegetation. By means of visual interpretation, little spectral differences are observed among those classes. The automatic classification of LANDSAT MSS data using maximum likelihood algorithm shows a 39% average error of omission and a 3.4% error of inclusion for the seven classes. The complexity of land uses in the study area, the large spectral variations of analyzed classes, and the low resolution of LANDSAT data influenced the classification results.

Dejesusparada, N.↗

A study of the effects of a cubic nonlinearity on a modern modal identification technique

The effect of a geometric nonlinearity on the Ibrahim Time Domain (ITD) modal data analysis technique has been studied using two analytically derived models and one laboratory model. Response data for the three models were analyzed by the ITD method. Indicators of nonlinear response were found which include harmonically related frequencies with repetitive mode shapes, clusters of frequencies in a narrow band around each harmonic, and variations in frequency with amplitude of oscillation. Also,for the cases studied, the presence of a nonlinearity has no detrimental effect on identifying linear responses. A potential for applying the algorithm to the identification of a certain class of nonlinear system was indicated.

Horta, L. G.↗

Safe microburst penetration techniques: A deterministic, nonlinear, optimal control approach

A relatively large amount of computer time was used for the calculation of a optimal trajectory, but it is subject to reduction with moderate effort. The Deterministic, Nonlinear, Optimal Control algorithm yielded excellent aircraft performance in trajectory tracking for the given microburst. It did so by varying the angle of attack to counteract the lift effects of microburst induced airspeed variations. Throttle saturation and aerodynamic stall limits were not a problem for the case considered, proving that the aircraft's performance capabilities were not violated by the given wind field. All closed loop control laws previously considered performed very poorly in comparison, and therefore do not come near to taking full advantage of aircraft performance.

Psiaki, Mark L.↗

Numerical study of the vortex burst phenomenon for delta wings

A flux-difference splitting scheme is employed to compute low-speed flows over a delta wing for angles of attack from 0 to 40 deg as steady-state solutions to the three-dimensional, Reynolds-averaged Navier-Stokes equations in their thin-layer approximation. The finite-difference scheme is made spatially second-order accurate by applying a total variation diminishing-like discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Using first-order accurate Euler backward-time differencing, an efficient implicit algorithm is constructed, which combines approximate factorization in cross planes with a symmetric planar Gauss-Seidel relaxation in the remaining third spatial direction. The geometry of the thin (maximum thickness is 0.021), slender (aspect ratio is unity), sharp-edged delta wing is taken from Hummel's (1967, 1978) wind tunnel model. Over the entire angle-of-attack range, the computed values of lift and pitching moment are in good agreement with the experimental data. Also details of the flow-fieldlike spanwise surface pressure distributions compare well with the experiment. Computed flow-field results with a bubble-type vortex burst are analyzed in detail.

Hartwich, PETER-M.↗

Data indexing techniques for the EUVE all-sky survey

This poster describes techniques developed for manipulating large full-sky data sets for the Extreme Ultraviolet Explorer project. The authors have adapted the quatrilateralized cubic sphere indexing algorithm to allow us to efficiently store and process several types of large data sets, such as full-sky maps of photon counts, exposure time, and count rates. A variation of this scheme is used to index sparser data such as individual photon events and viewing times for selected areas of the sky, which are eventually used to create EUVE source catalogs.

Lewis, J.↗

NOAA AVHRR Land Surface Albedo Algorithm Development

The primary objective of this research is to develop a surface albedo model for the National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVHRR). The primary test site is the Konza prairie, Kansas (U.S.A.), used by the International Satellite Land Surface Climatology Project (ISLSCP) in the First ISLSCP Field Experiment (FIFE). In this research, high spectral resolution field spectrometer data was analyzed to simulate AVHRR wavebands and to derive surface albedos. Development of a surface albedo algorithm was completed by analysing a combination of satellite, field spectrometer, and ancillary data. Estimated albedos from the field spectrometer data were compared to reference albedos derived using pyranometer data. Variations from surface anisotropy of reflected solar radiation were found to be the most significant albedo-related error. Additional error or sensitivity came from estimation of a shortwave mid-IR reflectance (1.3-4.0 micro-m) using the AVHRR red and near-IR bands. Errors caused by the use of AVHRR spectral reflectance to estimate both a total visible (0.4-0.7 micro-m) and near-IR (0.7-1.3 micro-m) reflectance were small. The solar spectral integration, using the derived ultraviolet, visible, near-IR and SW mid-IR reflectivities, was not sensitive to many clear-sky changes in atmospheric properties and illumination conditions.

Toll, D. L.↗

Application of an Ensemble Smoother to Precipitation Assimilation

Assimilation of precipitation in a global modeling system poses a special challenge in that the observation operators for precipitation processes are highly nonlinear. In the variational approach, substantial development work and model simplifications are required to include precipitation-related physical processes in the tangent linear model and its adjoint. An ensemble based data assimilation algorithm "Maximum Likelihood Ensemble Smoother (MLES)" has been developed to explore the ensemble representation of the precipitation observation operator with nonlinear convection and large-scale moist physics. An ensemble assimilation system based on the NASA GEOS-5 GCM has been constructed to assimilate satellite precipitation data within the MLES framework. The configuration of the smoother takes the time dimension into account for the relationship between state variables and observable rainfall. The full nonlinear forward model ensembles are used to represent components involving the observation operator and its transpose. Several assimilation experiments using satellite precipitation observations have been carried out to investigate the effectiveness of the ensemble representation of the nonlinear observation operator and the data impact of assimilating rain retrievals from the TMI and SSM/I sensors. Preliminary results show that this ensemble assimilation approach is capable of extracting information from nonlinear observations to improve the analysis and forecast if ensemble size is adequate, and a suitable localization scheme is applied. In addition to a dynamically consistent precipitation analysis, the assimilation system produces a statistical estimate of the analysis uncertainty.

Zhang, Sara↗

Verification of a New NOAA/NSIDC Passive Microwave Sea-Ice Concentration Climate Record

A new satellite-based passive microwave sea-ice concentration product developed for the National Oceanic and Atmospheric Administration (NOAA)Climate Data Record (CDR) programme is evaluated via comparison with other passive microwave-derived estimates. The new product leverages two well-established concentration algorithms, known as the NASA Team and Bootstrap, both developed at and produced by the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC). The sea ice estimates compare well with similar GSFC products while also fulfilling all NOAA CDR initial operation capability (IOC) requirements, including (1) self describing file format, (2) ISO 19115-2 compliant collection-level metadata,(3) Climate and Forecast (CF) compliant file-level metadata, (4) grid-cell level metadata (data quality fields), (5) fully automated and reproducible processing and (6) open online access to full documentation with version control, including source code and an algorithm theoretical basic document. The primary limitations of the GSFC products are lack of metadata and use of untracked manual corrections to the output fields. Smaller differences occur from minor variations in processing methods by the National Snow and Ice Data Center (for the CDR fields) and NASA (for the GSFC fields). The CDR concentrations do have some differences from the constituent GSFC concentrations, but trends and variability are not substantially different.

Passive Microwave↗

GOES-16 ABI Navigation Assessment

The US Geostationary Operational Environmental Satellite – R Series (GOES-R) was launched on November 19, 2016and was designated GOES-16 upon reaching geostationary orbit ten days later. After checkout and calibration, GOES-16 was relocated to its operational location of 75.2 degrees west and officially became GOES East on December 18, 2017. The Advanced Baseline Imager (ABI) is the primary instrument on the GOES-R series for imaging Earth's surface and atmosphere to significantly improve the detection and observation of severe environmental phenomena. A team supporting the GOES-R Flight Project at NASA's Goddard Space Flight Center developed algorithms and software for independent verification of ABI Image Navigation and Registration (INR), which became known as the INR Performance Assessment Tool Set (IPATS). In this paper, we will briefly describe IPATS on top concept level, and then introduce the Landsat chips, chip registration algorithms, and how IPATS measurements are filtered. We present GOES-16 navigation (NAV) errors from flight data from January 2017 to May 2018. The results show a) IPATS characterized INR variations throughout the post-launch test phase; and b) ABI INR has improved over time as post-launch tests were performed and corrections applied. Finally, we will describe how estimated NAV errors have been used to assess and understand satellite attitude anomalies and scale errors etc. This paper shows that IPATS is an effective tool for assessing and improving GOES-16 ABI INR and is also useful for INR long-term monitoring.

GOES-R↗

ANC of Radiated Sound from a Small UAS Rotor using Virtual Error Sensors

Small unmanned aerial systems (UAS) are becoming increasingly common for private, military, and commercial use, increasing community noise exposure. Reducing the noise produced by UAS could help improve community acceptance. Active noise control (ANC) might be used to attenuate noise produced by UAS, however, traditional ANC systems require a physical sensor in the far-field, which is not feasible. This paper assesses a virtual error sensing (VES) method that eliminates the need for a far-field sensor. This paper describes the proposed VES strategy, and presents numerical simulations and experimental results that highlight the benefits and limitations of the approach. Results for the VES system alone and with an ANC approach are presented and discussed. Experimental testing focused on attenuating the tonal noise produced by one 2-bladed rotor with a tip radius of 11.9 cm. Pressure variations caused by blade rotation were measured in the near and far-field using electret microphones and externally polarized condenser microphones, respectively. The filtered-x least mean squares algorithm was used in conjunction with the VES system to attenuate the far-field response. Experimental results show reductions between 6-13 dB at varying far-field locations and rotation rates.

virtual error sensor↗

Short timescale stability monitoring of GOES-16 visible channels based on daily inter-calibration events with VIIRS and MODIS

The CERES project uses geostationary (GEO) imager derived broadband fluxes to estimate the diurnal flux variations between the Terra (10:30 LT) and Aqua (1:30 LT) flux obs. The GEO derived cloud properties are utilized for narrowband to broadband and angular directional model (ADM) algorithms. The GEO imager channels are first inter-calibrated to the Aqua-MODIS C6.1 calibration reference to ensure that the retrieved clouds and fluxes will be consistent. The CERES project has a 2-month latency, all GEO, MODIS and VIIRS calibration anomalies must be identified and mitigated. Convert the monthly GEO ray-matching algorithm into a daily monitoring algorithm to catch any daily calibration anomalies.

geostationary sensor calibration↗