Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process”

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,405 records · Page 78

IMFIT: A FAST, FLEXIBLE NEW PROGRAM FOR ASTRONOMICAL IMAGE FITTING

I describe a new, open-source astronomical image-fitting program called IMFIT, specialized for galaxies but potentially useful for other sources, which is fast, flexible, and highly extensible. A key characteristic of the program is an object-oriented design that allows new types of image components (two-dimensional surface-brightness functions) to be easily written and added to the program. Image functions provided with IMFIT include the usual suspects for galaxy decompositions (Sérsic, exponential, Gaussian), along with Core-Sérsic and broken-exponential profiles, elliptical rings, and three components that perform line-of-sight integration through three-dimensional luminosity-density models of disks and rings seen at arbitrary inclinations. Available minimization algorithms include Levenberg-Marquardt, Nelder-Mead simplex, and Differential Evolution, allowing trade-offs between speed and decreased sensitivity to local minima in the fit landscape. Minimization can be done using the standard χ{sup 2} statistic (using either data or model values to estimate per-pixel Gaussian errors, or else user-supplied error images) or Poisson-based maximum-likelihood statistics; the latter approach is particularly appropriate for cases of Poisson data in the low-count regime. I show that fitting low-signal-to-noise ratio galaxy images using χ{sup 2} minimization and individual-pixel Gaussian uncertainties can lead to significant biases in fitted parameter values, which are avoided if a Poisson-based statistic is used; this is true even when Gaussian read noise is present.

79 ASTRONOMY AND ASTROPHYSICS↗

Raytracing the long trace profiler

The optical design of the Long Trace Profiler optical system is explored with a commercial raytrace program, Zemax OpticStudio™ (ZOS)1, with the intent of finding and correcting sources of systematic error. ZOS provides both geometric raytracing tools and physical optics Gaussian beam propagation and diffraction image calculation tools, and two design modes, sequential component (SC) and non-sequential component (NSC) that are optimized for different aspects of the design process. The original LTP-II system employs a singlet lens with a 1250 mm focal length. It is optimized to provide minimum distortion over a surface slope angle range of ±5 mrad. Using the ZOS tools, we are able to simulate ghost rays that produce distortion in the beam spot image and can minimize the distortion by deliberate misalignment of the beamsplitter (BS) components. Unfortunately, the reference beam is compromised because of the component tilts. The most recent LTP500 system design simplifies the optical system and makes the reference beam usable again, even with misalignment of the polarizing beamsplitter (PBS). Two lenses are designed for the LTP500 - a cemented doublet that has been fabricated, and a singlet with one aspheric surface. Both have focal distances of 500 mm with an expanded angular measurement range of ±10 mrad. The aspheric singlet provides superior performance. ZOS allows the import of wavefront measurement data produced by commercial interferometer software. We apply the wavefront error measurement from the cemented doublet to the model to show that the 19nm RMS wavefront error needs to be improved by at least a factor of ten in order to reduce the systematic error to a level that will allow the LTP to approach its design limit of a few tens of nanoradians.

Takacs, Peter Z.↗

UQ Toolkit v 2.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterizaton and propagation of uncertainties in computational models. For the characterization of uncertainties, Bayesian inference tools are provided to infer uncertain model parameters, as well as Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, and also Karhunen-Loève expansions for representing stochastic processes. Uncertain parameters are treated as random variables and represented with Polynomial Chaos expansions (PCEs). The library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Safta, Cosmin↗

Laser amplification of incoherent radiation

The amplification of noise in a laser amplifier is treated theoretically. The model for the active medium and its description using density-matrix techniques, are taken from the theory of laser operation. The spectral behavior of the radiation in the nonlinear regime is studied and the formalism is written from the onset in the frequency domain. The statistics of the light are gradually modified by the nonlinear amplification process, and expressions are derived for the rate of change of fluctuations in intensity as a measure of statistical changes. In addition, the range of validity of Litvak's Gaussian-statistics approximation is discussed. In the homogeneous-broadening case, the evolution of initially broadband Gaussian radiation toward quasimonochromatic oscillations with laserlike statistics is explored in several numerical examples. The connections of this study with the time-domain work on self-pulsing in a ring-laser configuration, are established. Finally, spectral-narrowing and -rebroadening effects in Doppler-broadened media are discussed both analytically and with numerical examples. These examples show the distinct contribution of pulsations in the population ('Raman-type terms'), and saturation phenomena.

Menegozzi, L. N.↗

A fiber-optic probe for particle sizing in concentrated suspensions

A fiber-optic probe employing two monomode optical fibers, one for transmitting a Gaussian laser beam to the scattering volume and the second, positioned at some backscatter angle, for receiving the scattered light is described. Performance and suitability of the system for a process control environment is assessed by studying a suspension of polystyrene latex particles over a wide range of sizes and concentrations. The results show that the probe is ideal for a process control environment in industrial and laboratory applications. Particle size is recovered, without any additional corrections for multiple light scattering, in concentrations containing up to 10 percent solids of 39-nm polystyrene latex spheres.

Dhadwal, Harbans S.↗

Preliminary structural control results from the Middeck Active Control Experiment (MACE)

Results are presented of on-going closed-loop ground experiments on the MACE test article, the objective of which is to investigate the extent to which closed-loop behavior of flexible spacecraft in zero gravity can be predicted, as well as to examine orbit system identification and control reconfiguration. The MACE hardware consists of three torque wheels, a two-axis gimballing payload, inertial sensors, and a flexible support structure. With the acquisition of a second payload, this is to represent a multiple payload platform with significant structural flexibility. When linear quadratic Gaussian control is used, payload pointing accuracy is improved by an order of magnitude when disturbed by a broadband torque disturbance. The successes and failures of the design and implementation process are discussed.

Miller, David W.↗

A Comparison of Three Random Number Generators for Aircraft Dynamic Modeling Applications

Three random number generators, which produce Gaussian white noise sequences, were compared to assess their suitability in aircraft dynamic modeling applications. The first generator considered was the MATLAB (registered) implementation of the Mersenne-Twister algorithm. The second generator was a website called Random.org, which processes atmospheric noise measured using radios to create the random numbers. The third generator was based on synthesis of the Fourier series, where the random number sequences are constructed from prescribed amplitude and phase spectra. A total of 200 sequences, each having 601 random numbers, for each generator were collected and analyzed in terms of the mean, variance, normality, autocorrelation, and power spectral density. These sequences were then applied to two problems in aircraft dynamic modeling, namely estimating stability and control derivatives from simulated onboard sensor data, and simulating flight in atmospheric turbulence. In general, each random number generator had good performance and is well-suited for aircraft dynamic modeling applications. Specific strengths and weaknesses of each generator are discussed. For Monte Carlo simulation, the Fourier synthesis method is recommended because it most accurately and consistently approximated Gaussian white noise and can be implemented with reasonable computational effort.

Grauer, Jared A.↗

Gaussian membership functions are most adequate in representing uncertainty in measurements

In rare situations, like fundamental physics, we perform experiments without knowing what their results will be. In the majority of real-life measurement situations, we more or less know beforehand what kind of results we will get. Of course, this is not the precise knowledge of the type 'the result will be between alpha - beta and alpha + beta,' because in this case, we would not need any measurements at all. This is usually a knowledge that is best represented in uncertain terms, like 'perhaps (or 'most likely', etc.) the measured value x is between alpha - beta and alpha + beta.' Traditional statistical methods neglect this additional knowledge and process only the measurement results. So it is desirable to be able to process this uncertain knowledge as well. A natural way to process it is by using fuzzy logic. But, there is a problem; we can use different membership functions to represent the same uncertain statements, and different functions lead to different results. What membership function do we choose? In the present paper, we show that under some reasonable assumptions, Gaussian functions mu(x) = exp(-beta(x(exp 2))) are the most adequate choice of the membership functions for representing uncertainty in measurements. This representation was efficiently used in testing jet engines to airplanes and spaceships.

Kreinovich, V.↗

Reconfigurable unitary transformations of optical beam arrays

Spatial transformations of light are ubiquitous in optics, with examples ranging from simple imaging with a lens to quantum and classical information processing in waveguide meshes. Multi-plane light converter (MPLC) systems have emerged as a platform that promises completely general spatial transformations, i.e., a universal unitary. However, until now, MPLC systems have demonstrated transformations that are far from general, e.g., converting from a Gaussian to Laguerre-Gauss mode. Here, we demonstrate the promise of an MLPC, the ability to impose an arbitrary unitary transformation that can be reconfigured dynamically. Specifically, we consider transformations on superpositions of parallel free-space beams arranged in an array, which is a common information encoding in photonics. We experimentally test the full gamut of unitary transformations for a system of two parallel beams and make a map of their fidelity. We obtain an average transformation fidelity of 0.85 ± 0.03. This high-fidelity suggests that MPLCs are a useful tool for implementing the unitary transformations that comprise quantum and classical information processing.

47 OTHER INSTRUMENTATION↗

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗

Geologic evolution of the terrestrial planets

The paper presents a geologic comparison of the terrestrial planets Mercury, Venus, Earth, the Moon and Mars, in the light of the recent photogeologic and other evidence gathered by satellites, and discusses the relationships between their regional terrain types, ages, and planetary evolution. The importance of the two fundamental processes, impact cratering and volcanism, which had formed these planets are stressed and the factors making the earth unique, such as high planetary evolution index (PEI), dynamic geological agents and the plate tectonics, are pointed out. The igneous processes which dominate earth and once existed on the others are outlined together with the planetary elevations of the earth which has a bimodal distribution, the moon which has a unimodal Gaussian distribution and Mars with a distribution intermediate between the earth and moon. Questions are raised concerning the existence of a minimum planetary mass below which mantle convection will not cause lithospheric rifting, and as to whether each planet follows a separate path of evolution depending on its physical properties and position within the solar system.

Head, J. W.↗

Robust contour decomposition using a constant curvature criterion

The problem of decomposing an extended boundary or contour into simple primitives is addressed with particular emphasis on Laplacian-of-Gaussian (LoG) zero-crossing contours. A technique is introduced for partitioning such contours into constant curvature segments. A nonlinear `blip' filter matched to the impairment signature of the curvature computation process, an overlapped voting scheme, and a sequential contiguous segment extraction mechanism are used. This technique is insensitive to reasonable changes in algorithm parameters and robust to noise and minor viewpoint-induced distortions in the contour shape, such as those encountered between stereo image pairs. The results vary smoothly with the data, and local perturbations induce only local changes in the result. Robustness and insensitivity are experimentally verified.

Wuescher, Daniel M.↗

Is the Non-Dipole Magnetic Field Random?

Statistical modelling of the Earth's magnetic field B has a long history. In particular, the spherical harmonic coefficients of scalar fields derived from B can be treated as Gaussian random variables. In this paper, we give examples of highly organized fields whose spherical harmonic coefficients pass tests for independent Gaussian random variables. The fact that coefficients at some depth may be usefully summarized as independent samples from a normal distribution need not imply that there really is some physical, random process at that depth. In fact, the field can be extremely structured and still be regarded for some purposes as random. In this paper, we examined the radial magnetic field B(sub r) produced by the core, but the results apply to any scalar field on the core-mantle boundary (CMB) which determines B outside the CMB.

Walker, Andrew D.↗

Energy Measurement Studies for CO2 Measurement with a Coherent Doppler Lidar System

The accurate measurement of energy in the application of lidar system for CO2 measurement is critical. Different techniques of energy estimation in the online and offline pulses are investigated for post processing of lidar returns. The cornerstone of the techniques is the accurate estimation of the spectrum of lidar signal and background noise. Since the background noise is not the ideal white Gaussian noise, simple average level estimation of noise level is not well fit in the energy estimation of lidar signal and noise. A brief review of the methods is presented in this paper.

Beyon, Jeffrey Y.↗

Extended Line Spread Function of TES Microcalorimeters with Au/Bi Absorbers

Microcalorimeters have the potential to provide line shapes well-described by a single Gaussian broadening term of few-eV width. This attribute makes the detectors especially well suited for x-ray astrophysics observations; however, low-level non-Gaussian broadening terms are expected and must be characterized. These terms depend on the composition of the x-ray absorber, the detailed x-ray absorption physics, the device thermalization processes, and the incident x-ray energy. Here we present the first measurements targeted at understanding the extended line-spread function (LSF) of x-ray microcalorimeter pixels under development for the X-ray Integral Field Unit (XIFU) on the Athena X-ray Observatory. These pixels are comprised of Mo/Au transition-edge sensors (TESs) with overhanging electroplated Au/Bi absorbers. We measured the line shapes using monochromatic x-ray sources with < 1 eV width at several x-ray energies (0.85, 0.93, 1.25, 1.5, 5.4, and 8.0 keV) across the instrument bandpass (0.3-12 keV) and modeled the line profiles. These results are compared to the extended LSF of the Hitomi microcalorimeter pixels that used HgTe absorbers.

Superconducting devices↗

Investigation of post-breakup Coulomb acceleration using a trajectory model

Intermediate mass fragments ejected during the deexcitation of excited projectilelike fragments may promptly decay following ejection; the daughter particles that are subsequently produced are subject to interactions with the residual nucleus that affect final-state observables, a process herein referred to as post-breakup Coulomb acceleration. A simple classical Coulomb interaction model was used to study modification of 8 Be (2 + ), 5 Li (3/2 – ), 7 Li (7/2 – ), 7 Be (7/2 – ), and states in 12 B emitted following heavy-ion collisions of 28 Si + 12 C at 35 MeV/nucleon. Here, in contrast to previous work studying 8 Be (2 + ), excellent agreement between simulation and experiment was obtained using only Coulomb forces when either a Lorentzian or R-matrix line shape was used to describe the initial relative energy rather than a Gaussian. In consideration of the obtained results, improvements to the model and evaluation of experimental data are discussed as future directions, but it was concluded that the effects observed in the present data can be accurately described using only elements of classical mechanics and that the process is largely understood for a wide range of state lifetimes and mass and charge (a)symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of the EMCCD point-source response using x-rays

Electron Multiplying Charge Coupled Devices, EMCCD are used as x-ray detectors. The NSLS-II Soft Inelastic x-ray Scattering (SIX) beam line has two EMCCDs for x-ray detection in the spectrometer arm. The spectrometer with high resolving power disperses x-rays vertically. The x-ray vertical position on the sensor plane is related to its energy. This allows for very accurate x-ray energy measurements through x-ray coordinates. X-rays interact with silicon and create a number of electron-hole pairs proportional to the x-ray energy. Electrons drift and diffuse toward pixel gates and are collected there. The diffused electrons form a charge cloud distributed over several neighboring pixels. This charge sharing enables coordinate measurements with accuracy better than the pixel pitch. The charge distribution shape has to be taken into account to achieve ultimate accuracy in coordinate measurements. In this paper, we present a method of the charge distribution shape analysis and demonstrate its applications. The drift and diffusion of electrons ftom the point of generation to pixel gates results in the bell-shaped electron cloud usually approximated by Gaussian shape. The number of electrons collected under a pixel is proportional to the shape function integral. These electron packets get transferred to the sense node of the output amplifier. The transfer process could introduce distortions to the original charge distribution. For example, during transfers, electrons in the packet could be exposed to traps if they are present in the sensor. The trapping and later the release processes distort the apparent shape of the charge distribution. Therefore, deviations of the charge distribution shape ftom the originally symmetrical form can indicate the presence of trap centers in the sensor and can be used for sensor diagnostics.

47 OTHER INSTRUMENTATION↗