Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gaussian processes”

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

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗

File Specification for GEOS-5 FP-IT (Forward Processing for Instrument Teams)

The GEOS-5 FP-IT Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP-IT incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP-IT uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The analysis is performed at a horizontal resolution of 0.625-degree longitude by 0.5-degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational GEOS-5 FP-IT processing in support of the EOS instrument teams. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary.

3D Var↗

An efficient system for reliably transmitting image and video data over low bit rate noisy channels

This research project is intended to develop an efficient system for reliably transmitting image and video data over low bit rate noisy channels. The basic ideas behind the proposed approach are the following: employ statistical-based image modeling to facilitate pre- and post-processing and error detection, use spare redundancy that the source compression did not remove to add robustness, and implement coded modulation to improve bandwidth efficiency and noise rejection. Over the last six months, progress has been made on various aspects of the project. Through our studies of the integrated system, a list-based iterative Trellis decoder has been developed. The decoder accepts feedback from a post-processor which can detect channel errors in the reconstructed image. The error detection is based on the Huber Markov random field image model for the compressed image. The compression scheme used here is that of JPEG (Joint Photographic Experts Group). Experiments were performed and the results are quite encouraging. The principal ideas here are extendable to other compression techniques. In addition, research was also performed on unequal error protection channel coding, subband vector quantization as a means of source coding, and post processing for reducing coding artifacts. Our studies on unequal error protection (UEP) coding for image transmission focused on examining the properties of the UEP capabilities of convolutional codes. The investigation of subband vector quantization employed a wavelet transform with special emphasis on exploiting interband redundancy. The outcome of this investigation included the development of three algorithms for subband vector quantization. The reduction of transform coding artifacts was studied with the aid of a non-Gaussian Markov random field model. This results in improved image decompression. These studies are summarized and the technical papers included in the appendices.

Costello, Daniel J., Jr.↗

State estimation Kalman filter using optical processings Noise statistics known

Reference is made to a study by Casasent et al. (1983), which gave a description of a frequency-multiplexed acoustooptic processor and showed how it was capable of performing all the individual operations required in Kalman filtering. The data flow and organization of all required operations however, were not detailed in that study. Consideration is given here to a simpler Kalman filter state estimation problem. Equally spaced time-sampled intervals (k times T sub s, with k the iterative time index) are assumed. It is further assumed that the system noise vector w and the measurement noise vector v are uncorrelated and Gaussian distributed and that the noise statistics (Q and R) and the system model (Phi, Gamma, H) are known. The error covariance matrix P and the extrapolated error covariance matrix M can thus be precomputed and the Kalman gain matrix K sub k can be precomputed and stored for each input time sample.

Jackson, J.↗

Rossby wavetrains in the stratosphere forced by localised disturbances in the troposphere

Studies have linked elements of the stratospheric circulation with particular localized features in the tropospheric circulation. This suggests a study of the response of the stratosphere to forcing by localized disturbances in the troposphere. A multilevel, primitive equation model is used of the stratosphere and mesosphere the height of whose lower boundary at 300 mb can be prescribed. Localized height disturbances which grow to steady amplitude are applied at this lower boundary, and the response of the initially axially symmetric stratosphere is studied. The perturbation, centered at 45 deg N, has a Gaussian distribution with a half width of 15 deg., corresponding roughly to the size of persistent troughs and blocking ridges in the troposphere. The forcing is nearly at fully amplitude 10 days after being switched on, and thereafter remains steady. Two types of experiments are conducted: in one the forcing is of small amplitude (100 gpm) and in the other it is of large amplitude (600 gpm). These pairs of experiments are compared to determine how nonlinear processes affect the perturbation fields. This is done locally by defining the perturbation response to the forcing as an anomaly, i.e., as a departure from the response in a controlled experiment in which no asymmetric forcing was applied at the lower boundary of the model. Experiments have been conducted for a number of atmospheric states obtained as zonal means of observations made by a stratospheric sounding unit (SSU). In this summary, results for the zonal mean flow on January 19, 1982 are outlined.

Marks, C. J.↗

On the optimality of code options for a universal noiseless coder

A universal noiseless coding structure was developed that provides efficient performance over an extremely broad range of source entropy. This is accomplished by adaptively selecting the best of several easily implemented variable length coding algorithms. Custom VLSI coder and decoder modules capable of processing over 20 million samples per second are currently under development. The first of the code options used in this module development is shown to be equivalent to a class of Huffman code under the Humblet condition, other options are shown to be equivalent to the Huffman codes of a modified Laplacian symbol set, at specified symbol entropy values. Simulation results are obtained on actual aerial imagery, and they confirm the optimality of the scheme. On sources having Gaussian or Poisson distributions, coder performance is also projected through analysis and simulation.

Yeh, Pen-Shu↗

Interstellar Isomers: The Importance of Bonding Energy Differences

We present strong detections of methyl cyanide (CH3CN), vinyl cyanide (CH2CHCN), ethyl cyanide (CH3CH2CN) and cyanodiacetylene (HC4CN) molecules with the Green Bank Telescope (GBT) toward the Sgr B2(N) molecular cloud. Attempts to detect the corresponding isocyanide isomers were only successful in the case of methyl isocyanide (CH3NC) for its J(sub K) = 1(sub 0) - 0(sub 0) transition, which is the first interstellar report of this line. To determine the spatial distribution of CH3NC, we used archival Berkeley-Illinois-Maryland Association (BIMA) array data for the J(sub K) = 1(sub 0) - 0(sub 0) transitions but no emission was detected. From ab initio calculations, the bonding energy difference between the cyanide and isocyanide molecules is greater than 8500 per centimeter (greater than 12,000 K). Thus, cyanides are the more stable isomers and would likely be formed more preferentially over their isocyanide counterparts. That we detect CH3NC emission with a single antenna (Gaussian beamsize(omega(sub B))=1723 arcsec(sup 2)) but not with an interferometer (omega(sub b)=192 arcsec(sup 2)), strongly suggests that CH3NC has a widespread spatial distribution toward the Sgr B2(N) region. Other investigators have shown that CH3CN is present both in the LMH hot core of Sgr B2(N) and in the surrounding medium, while we have shown that CH3NC appears to be deficient in the LMH hot core. Thus, largescale, non-thermal processes in the surrounding medium may account for the conversion of CH3CN to CH3NC while the LMH hot core, which is dominated by thermal processes, does not produce a significant amount of CH3NC. Ice analog experiments by other investigators have shown that radiation bombardment of CH3CN can produce CH3NC, thus supporting our observations. We conclude that isomers separated by such large bonding energy differences are distributed in different interstellar environments, making the evaluation of column density ratios between such isomers irrelevant unless it can be independently shown that these species are co-spatial.

Remijan, Anthony J.↗

Using Statistical Multivariable Models to Understand the Relationship Between Interplanetary Coronal Mass Ejecta and Magnetic Flux Ropes

In-situ measurements of interplanetary coronal mass ejections (ICMEs) display a wide range of properties. A distinct subset, "magnetic clouds" (MCs), are readily identifiable by a smooth rotation in an enhanced magnetic field, together with an unusually low solar wind proton temperature. In this study, we analyze Ulysses spacecraft measurements to systematically investigate five possible explanations for why some ICMEs are observed to be MCs and others are not: i) An observational selection effect; that is, all ICMEs do in fact contain MCs, but the trajectory of the spacecraft through the ICME determines whether the MC is actually encountered; ii) interactions of an erupting flux rope (PR) with itself or between neighboring FRs, which produce complex structures in which the coherent magnetic structure has been destroyed; iii) an evolutionary process, such as relaxation to a low plasma-beta state that leads to the formation of an MC; iv) the existence of two (or more) intrinsic initiation mechanisms, some of which produce MCs and some that do not; or v) MCs are just an easily identifiable limit in an otherwise corntinuous spectrum of structures. We apply quantitative statistical models to assess these ideas. In particular, we use the Akaike information criterion (AIC) to rank the candidate models and a Gaussian mixture model (GMM) to uncover any intrinsic clustering of the data. Using a logistic regression, we find that plasma-beta, CME width, and the ratio O(sup 7) / O(sup 6) are the most significant predictor variables for the presence of an MC. Moreover, the propensity for an event to be identified as an MC decreases with heliocentric distance. These results tend to refute ideas ii) and iii). GMM clustering analysis further identifies three distinct groups of ICMEs; two of which match (at the 86% level) with events independently identified as MCs, and a third that matches with non-MCs (68 % overlap), Thus, idea v) is not supported. Choosing between ideas i) and iv) is more challenging, since they may effectively be indistinguishable from one another by a single in-situ spacecraft. We offer some suggestions on how future studies may address this.

Riley, P.↗

Detecting perceptual groupings in textures by continuity considerations

A generalization is presented for the second derivative of a Gaussian D(sup 2)G operator to apply to problems of perceptual organization involving textures. Extensions to other problems of perceptual organization are evident and a new research direction can be established. The technique presented is theoretically pleasing since it has the potential of unifying the entire area of image segmentation under the mathematical notion of continuity and presents a single algorithm to form perceptual groupings where many algorithms existed previously. The eventual impact on both the approach and technique of image processing segmentation operations could be significant.

Greene, Richard J.↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Detection of signals by weighted integrate-and-dump filter

A Weighted Integrate and Dump Filter (WIDF) is presented that results in reducing those losses in telemetry symbol signal to noise ratio (SNR) which occur in digital Integrate and Dump Filters (IDFs) when the samples are not phase locked to the input data symbol clock. The Minimum Mean Square Error (MMSE) criterion is used to derive a set of weights for approximating the analog integrate and dump filter, which is the matched filter for detection of signals in additive white Gaussian noise. This new digital matched filter results in considerable performance improvement compared to unweighted digital matched filters. An example is presented for a sampling rate of four times the symbol rate. As the sampling offset (or phase) varies with respect to the data symbol boundaries, the output SNR varies 1 dB for an unweighted IDF, but only 0.3 dB for the optimum WIDF, averaged over random data patterns. This improvement in performance relative to unweighted IDF means that significantly lower sampling and processing rates can be used for given telemetry symbol rates, resulting in reduced system cost.

Sadr, R.↗

The Barberplaid illusion: Even Perceived Plaid Motion is Biased by Asymmetric Apertures

Our results show that the perceived direction of motion of plaids windowed by asymmetric spatial Gaussians is biased toward the long axis of the window. The bias increases as the relative angle between the plaid motion and the window increases, peaks at a relative angle of about 40 degrees, and then decreases. The peak bias was 14 degrees for a spatial frequency of 0.6 cpd and a window aspect ratio of 4.0. The biases increase as the window is elongated and decrease as the component spatial frequency increases. We tested the predictions of several models of human motion processing (cross correlation, motion energy, intersection of constraints, and vector sum), and show that none of these can predict our data. These results suggest that spatial integration of motion signals plays a crucial role in the perception of plaid motion.

Beutter, Brent↗

Very long baseline interferometry observations of the RS Canum Venaticorum system HR 5110 at 8.4 GHz

The RS CVn binary system HR 5110 was observed with a well-calibrated four-element VLBI array at a frequency of 8.4 GHz. The total flux density, which was nearly constant at 32 + or - 2 mJy, was monitored throughout the experiment with a 20 km baseline interferometer allowing accurate visibilities to be calculated. The observed visibilities are consistent with a source unresolved on all the baselines. The maximum size of the equivalent Gaussian source is 1.4 milli-arcsec (FWHM), corresponding to a linear size of 1.1 x 10 to the 12th cm at 52 pc. This maximum size is comparable with the overall size of the binary system. The inferred lower limit of the brightness temperature is 4 x 10 to the 8th K, which is higher than the coronal temperature of 10 to the 7th K measured in a soft X-ray survey. This is consistent with a nonthermal mechanism for the radiation process of the moderate radio outburst observed.

Mutel, R. L.↗

Anomaly Detection in Power Quality at Data Centers

The goal during my internship at the National Center for Critical Information Processing and Storage (NCCIPS) is to implement an anomaly detection method through the StruxureWare SCADA Power Monitoring system. The benefit of the anomaly detection mechanism is to provide the capability to detect and anticipate equipment degradation by monitoring power quality prior to equipment failure. First, a study is conducted that examines the existing techniques of power quality management. Based on these findings, and the capabilities of the existing SCADA resources, recommendations are presented for implementing effective anomaly detection. Since voltage, current, and total harmonic distortion demonstrate Gaussian distributions, effective set-points are computed using this model, while maintaining a low false positive count.

Grichine, Art↗

Results of an integrated structure-control law design sensitivity analysis

Next generation air and space vehicle designs are driven by increased performance requirements, demanding a high level of design integration between traditionally separate design disciplines. Interdisciplinary analysis capabilities have been developed, for aeroservoelastic aircraft and large flexible spacecraft control for instance, but the requisite integrated design methods are only beginning to be developed. One integrated design method which has received attention is based on hierarchal problem decompositions, optimization, and design sensitivity analyses. This paper highlights a design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changess in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient that finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗

New Insights on Tropical Ozone from SHADOZ (Southern Hemisphere Additional Ozonesondes) Profiles

The SHADOZ (Southern Hemisphere Additional Ozonesondes) ozone sounding network was initiated in 1998 to improve the coverage of tropical in-situ ozone measurements for satellite validation, algorithm development and related process studies. Over 2000 soundings have been archived at the website, , for 12 stations: Ascension Island; Nairobi and Malindi, Kenya; Irene, South Africa; Reunion Island, Watukosek, Java; Fiji; Tahiti; American Samoa; San Cristobal, Galapagos; Natal, Brazil; Paramaribo, Surinam. Key results from SHADOZ will be described from among the following: 1) By using ECC sondes with similar procedures, 5-10% accuracy and precision (1-sigma) of the sonde total ozone measurement was achieved; 2) Week-to-week variability in tropospheric ozone is so great that statistics are frequently not Gaussian; most stations vary up to a factor of 3 in tropospheric column over the course of a year; 3) Longitudinal variability in tropospheric ozone profiles is a consistent feature, with a 10-15 DU column-integrated difference between Atlantic and Pacific sites; this causes a "zonal wave-one" feature in total ozone. 4) The ozone record from Paramaribo, Surinam (6N, 55W) is a marked contrast to southern tropical ozone because Surinam is often north of the Intertropical Convergence Zone; 5) Indian Ocean region pollution may contribute up to half of the excess ozone observed in the south tropical Atlantic paradox in the December-January-February period of the year.

Thompson, Anne M.↗

Effects of correlated noise on the full-spectrum combining and complex-symbol combining arraying techniques

The process of combining telemetry signals received at multiple antennas, commonly referred to as arraying, can be used to improve communication link performance in the Deep Space Network (DSN). By coherently adding telemetry from multiple receiving sites, arraying produces an enhancement in signal-to-noise ratio (SNR) over that achievable with any single antenna in the array. A number of different techniques for arraying have been proposed and their performances analyzed in past literature. These analyses have compared different arraying schemes under the assumption that the signals contain additive white Gaussian noise (AWGN) and that the noise observed at distinct antennas is independent. In situations where an unwanted background body is visible to multiple antennas in the array, however, the assumption of independent noises is no longer applicable. A planet with significant radiation emissions in the frequency band of interest can be one such source of correlated noise. For example, during much of Galileo's tour of Jupiter, the planet will contribute significantly to the total system noise at various ground stations. This article analyzes the effects of correlated noise on two arraying schemes currently being considered for DSN applications: full-spectrum combining (FSC) and complex-symbol combining (CSC). A framework is presented for characterizing the correlated noise based on physical parameters, and the impact of the noise correlation on the array performance is assessed for each scheme.

Vazirani, P.↗

Results of an integrated structure/control law design sensitivity analysis

A design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations is discussed. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changes in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient than finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗