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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Simulation of Inflated Pahoehoe Lava Flows

A new stochastic model simulates late-stage pahoehoe lobes where random processes dominate emplacement. The model prescribes probabilistic rules for determining where and when parcels of lava move within the lobe. Unlike a classical Brownian motion random walk, the model allows individual parcels to remain dormant, but fluid, for multiple time steps. The randomness of parcel volume transfers within the lobe interior as well as at the margins qualitatively reflects inflation processes observed in the field. The fraction of inflated volume to total volume increases with the total volume, with greater than 75% of the lobe volume contributed through inflation for typical lobes. The influence on planform shape and topographic cross-sectional profiles of total volume, source area and shape, topographic confinement, and sequential breakouts at the lobe margins, are all explored with the stochastic model. Each of these factors influences the overall lobe thickness and width. The model provides a means for assessing the relative importance of these processes through comparisons with field data. For the first time, Gaussian and parabolic functions are quantitatively fit to field measurements of pahoehoe lobes. Both functional forms provide adequate description of the cross-sectional flow shapes. When comparing simulated lobes to field data, sequential breakouts at the lobe margins are found to be an important process controlling the final topographic distribution of observed pahoehoe lobes.

modeling↗
DOE OSTI2018

Geometric Analysis Based Double Closed-loop Iterative Learning Control of Output PDF Shaping of Fiber Length Distribution in Refining Process

In order to improve the pulp quality and to reduce the energy consumption, the fiber length distribution (FLD) is generally employed as one of the important technological indexes in the refining process. Considering that the traditional mean and variance of fiber length are unable to adequately characterize the non-Gaussian distribution properties, this paper proposes a novel geometric analysis based double closed-loop iterative learning control (ILC) method for probability density function (PDF) shaping of output FLD in the refining process. Primarily, a RBF neural network (NN) with Gaussian-type is utilized to approximate the square root PDF in the inner loop, where the RBF basis function parameters (center and width) are tuned between any two adjacent batches by using an ILC law, and the subspace identification method can be applied to establish the state-space model of weight vector. Then, for the sake of accelerating the convergence rate of the closed-loop sysyem, a geometric analysis based ILC method is adopted in the outer loop. Lastly, both simulation and experiments demonstrate the effectiveness and practicability of the proposed approach.

42 ENGINEERING↗
NASA NTRS2004

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗
DOE OSTI2025

FNCL Enhancements Implementation (FY25 Annual Report)

The FNCL investigations team at Lawrence Livermore National Laboratory (LLNL) has completed research and development of hardware, signal processing, and analysis tools to enhance the measurement capabilities of both the current CAEN SyS VeryFuel Fast Neutron Collar (FNCL) instrument and a next-generation FNCL prototype. The team successfully built and commissioned the LLNL Demonstrator System: a fully integrated, three-panel detector system featuring higher segmentation, plastic scintillators (EJ-276D), Silicon Photomultipliers (SiPMs), no high-voltage requirement, a reduced electronic footprint, and the LLNL-developed Gaussian Mixture Model Pulse Shape Discrimination (GMM-PSD) signal processing. An extensive experimental campaign was conducted at LLNL’s Inherently Safe Subcritical Assembly (ISSA) facility using both the baseline FNCL and the LLNL Demonstrator. The campaign results validated system performance, calibration stability, and the effectiveness of advanced signal processing and analysis algorithms in a relevant environment.

and physical protection↗
DOE OSTI2017

Asymmetry in serial femtosecond crystallography data

Serial crystallography is an increasingly important approach to protein crystallography that exploits both X-ray free-electron laser (XFEL) and synchrotron radiation. Serial crystallography recovers complete X-ray diffraction data by processing and merging diffraction images from thousands of randomly oriented non-uniform microcrystals, of which all observations are partial Bragg reflections. Random fluctuations in the XFEL pulse energy spectrum, variations in the size and shape of microcrystals, integrating over millions of weak partial observations and instabilities in the XFEL beam position lead to new types of experimental errors. The quality of Bragg intensity estimates deriving from serial crystallography is therefore contingent upon assumptions made while modeling these data. Here it is observed that serial femtosecond crystallography (SFX) Bragg reflections do not follow a unimodal Gaussian distribution and it is recommended that an idealized assumption of single Gaussian peak profiles be relaxed to incorporate apparent asymmetries when processing SFX data. The phenomenon is illustrated by re-analyzing data collected from microcrystals of the Blastochloris viridis photosynthetic reaction center and comparing these intensity observations with conventional synchrotron data. The results show that skewness in the SFX observations captures the essence of the Wilson plot and an empirical treatment is suggested that can help to separate the diffraction Bragg intensity from the background.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Spin dynamics in driven composite multiferroics

A spin dynamics approach has been used to study the behavior of the magnetic spins and the electric pseudo-spins in a 1-D composite multiferroic chain with a linear magneto-electric coupling at the interface. The response is investigated with either external magnetic or electric fields driving the system. The spin dynamics is based on the Landau-Lifshitz-Gilbert equation. A Gaussian white noise is later added into the dynamic process to include the thermal effects. The interface requires a closer inspection of the magneto-electric effects. Thus, we construct a 2-D ladder model to describe the behavior of the magnetic spins and the electric pseudo-spins with different magneto-electric couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Predicted performance of a tangential viewing hard x-ray camera for the DIII-D high field side lower hybrid current drive experiment

High field side launch of lower hybrid current drive (LHCD) has improved accessibility and penetration over low field side launch on DIII-D. Simulations predict single pass absorption under a wide range of plasma conditions. Hard x-ray (HXR) measurement of LHCD generated fast electron bremsstrahlung (50–250 keV) will validate wave propagation and absorption. Emissivity profiles are recovered from one-dimensional inversion of HXR brightness to determine LH damping location, fast electron slowing down time, and some indication of the fast electron energy. The camera will be implemented by populating 32 tangential sightlines of the existing Gamma Ray Imager with Kromek SPEARTM Cadmium Zinc Telluride (CZT) detectors sensitive to 10–1000 keV photons with 10 keV energy resolution. Expected count rates allow for <0.5 ms time resolution. Pulses are processed using 50 ns shaping time Cremat CR-200 Gaussian shaping modules and are digitized by 25 MHz D-TACQ ACQ216 digitizers. The performance of the HXR camera is evaluated by comparing predicted fast electron density profiles and inverted synthetic brightnesses obtained from the ray-tracing/Fokker–Planck codes GENRAY/CQL3D. Inversions closely matched predicted fast electron profiles for a range of experimental parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗
DOE OSTI2018

A Comparison of Different Methods of Separation of Continuous Overlapping Spectral Lines

Several methods of separation of continuous overlapping spectral lines are compared. In the majority of methods, the profile of each line is modeled by a Gaussian or Lorentzian, and total measured spectrum z is processed. Number of lines N and their parameters are usually estimated by the method of derivatives; however, the differentiation of noisy spectrum z causes large errors. To improve the differentiation accuracy, it is proposed to use smoothing splines. In the Fourier-self-deconvolution method, apodization (artificial truncation of the interferogram) is used to resolve overlapping lines, which makes it possible to resolve lines, but at the expense of a significant decrease in their widths. In this work, reducing the widths of lines is not used for their resolution, but rather, true line profiles are reconstructed by minimizing the residual functional with the modified coordinate-descent method using the decremental-constraint technique, and also, for comparison, with the Nelder–Mead simplex method. In the Manoilov method, the parameters of lines (peaks) are determined from convolutions of spectrum derivatives with individual peaks. In that method, the notion of the degree of overlap has been also introduced. In this work, we introduce a generalized degree of overlap for the case in which the amplitudes, widths, and spacings between neighboring lines are different. Numerical illustrations are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Damping of surface pressure fluctuations in hypersonic turbulent flow past expansion corners

Surface pressure fluctuations of Mach 8 turbulent flow past a 2.5- and a 4.25-deg expansion corner maintained a Gaussian distribution but were severely attenuated by the expansion process. The pressure fluctuations did not recover to those of an equilibrium turbulent flow even though the mean pressures reached downstream inviscid values in four to six boundary-layer thicknesses. The fluctuations were convected with a velocity comparable to that on a flat plate, and they maintained their identities longer for the stronger expansion. The damping of pressure fluctuations at hypersonic Mach numbers, even by small corner angles, may be exploited in fatigue design.

Chung, Kung-Ming↗
DOE OSTI2015

Gaussian variational ansatz in the problem of anomalous sea waves: Comparison with direct numerical simulation

The nonlinear dynamics of an obliquely oriented wave packet on a sea surface is analyzed analytically and numerically for various initial parameters of the packet in relation to the problem of the so-called rogue waves. Within the Gaussian variational ansatz applied to the corresponding (1+2)-dimensional hyperbolic nonlinear Schrödinger equation (NLSE), a simplified Lagrangian system of differential equations is derived that describes the evolution of the coefficients of the real and imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative description of the process of nonlinear spatiotemporal focusing, which is one of the most probable mechanisms of rogue wave formation in random wave fields. The system of equations is integrated in quadratures, which allows one to better understand the qualitative differences between linear and nonlinear focusing regimes of a wave packet. Predictions of the Gaussian model are compared with the results of direct numerical simulation of fully nonlinear long-crested waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

Parallel Gaussian elimination of a block tridiagonal matrix using multiple microcomputers

The solution of a block tridiagonal matrix using parallel processing is demonstrated. The multiprocessor system on which results were obtained and the software environment used to program that system are described. Theoretical partitioning and resource allocation for the Gaussian elimination method used to solve the matrix are discussed. The results obtained from running 1, 2 and 3 processor versions of the block tridiagonal solver are presented. The PASCAL source code for these solvers is given in the appendix, and may be transportable to other shared memory parallel processors provided that the synchronization outlines are reproduced on the target system.

Blech, Richard A.↗
DOE OSTI2022

Dynamics of bottlebrush polymers

Bottlebrushes are an interesting class of polymers which shows intriguing material properties often associated with dynamics. While dynamical phenomena in linear polymers are well understood and existing theories can describe them in a good way, bottlebrush dynamics have only rarely been investigated. Therefore, we performed dielectric spectroscopy and quasi-elastic neutron scattering to study the dynamics of polydimethylsiloxane-based bottlebrush polymers, PDMS-g-PDMS focusing mostly on the segmental dynamics of the side chains. Comparing the relaxation times of the α – relaxation, tracked with dielectric spectroscopy, of bottlebrush polymers with those of their respective linear side chains show a slowing down once the side chains are attached to the backbone. This effect diminishes and finally vanishes with increasing side chain length. The time and length scale, offered by quasi-elastic neutron scattering, fits for the segmental dynamics together with faster processes. The Q -dependence of the segmental relaxation times allows to classify bottlebrush polymers as heterogenous including a non-Gaussian character. For such a dynamical system, the mean square displacement needs to be separated into single processes before an overall mean square displacement can be generated by applying the time temperature superposition principle.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Hardware Implementation of a Bilateral Subtraction Filter

A bilateral subtraction filter has been implemented as a hardware module in the form of a field-programmable gate array (FPGA). In general, a bilateral subtraction filter is a key subsystem of a high-quality stereoscopic machine vision system that utilizes images that are large and/or dense. Bilateral subtraction filters have been implemented in software on general-purpose computers, but the processing speeds attainable in this way even on computers containing the fastest processors are insufficient for real-time applications. The present FPGA bilateral subtraction filter is intended to accelerate processing to real-time speed and to be a prototype of a link in a stereoscopic-machine- vision processing chain, now under development, that would process large and/or dense images in real time and would be implemented in an FPGA. In terms that are necessarily oversimplified for the sake of brevity, a bilateral subtraction filter is a smoothing, edge-preserving filter for suppressing low-frequency noise. The filter operation amounts to replacing the value for each pixel with a weighted average of the values of that pixel and the neighboring pixels in a predefined neighborhood or window (e.g., a 9 9 window). The filter weights depend partly on pixel values and partly on the window size. The present FPGA implementation of a bilateral subtraction filter utilizes a 9 9 window. This implementation was designed to take advantage of the ability to do many of the component computations in parallel pipelines to enable processing of image data at the rate at which they are generated. The filter can be considered to be divided into the following parts (see figure): a) An image pixel pipeline with a 9 9- pixel window generator, b) An array of processing elements; c) An adder tree; d) A smoothing-and-delaying unit; and e) A subtraction unit. After each 9 9 window is created, the affected pixel data are fed to the processing elements. Each processing element is fed the pixel value for its position in the window as well as the pixel value for the central pixel of the window. The absolute difference between these two pixel values is calculated and used as an address in a lookup table. Each processing element has a lookup table, unique for its position in the window, containing the weight coefficients for the Gaussian function for that position. The pixel value is multiplied by the weight, and the outputs of the processing element are the weight and pixel-value weight product. The products and weights are fed to the adder tree. The sum of the products and the sum of the weights are fed to the divider, which computes the sum of products the sum of weights. The output of the divider is denoted the bilateral smoothed image. The smoothing function is a simple weighted average computed over a 3 3 subwindow centered in the 9 9 window. After smoothing, the image is delayed by an additional amount of time needed to match the processing time for computing the bilateral smoothed image. The bilateral smoothed image is then subtracted from the 3 3 smoothed image to produce the final output. The prototype filter as implemented in a commercially available FPGA processes one pixel per clock cycle. Operation at a clock speed of 66 MHz has been demonstrated, and results of a static timing analysis have been interpreted as suggesting that the clock speed could be increased to as much as 100 MHz.

Huertas, Andres↗
DOE OSTI2019

An Adaptive Bayesian Parameter Estimation of a Synchronous Generator Under Gross Errors

Polynomial-chaos-expansion-based surrogate models have recently been advocated in the literature for power system dynamic parameter estimation. Regarding the estimation of the uncertain generator parameters, a Bayesian inference framework has been proposed based on a polynomial-based reduced-order representation of the synchronous machines using assumed parameter values. Then, the non-Gaussian posterior probability distribution functions (pdfs) of these parameters are recovered through the stochastic sampling approach efficiently. However, facing very large parameter errors, the reliability of the surrogate model decreases, yielding biased estimation results. To overcome this problem, here we develop hierarchical Bayesian inference framework that processes the measurements provided by phasor measurement units, while making use of multifidelity surrogates together with the importance sampling method. The latter allows us to estimate in an efficient manner the posterior pdfs of the uncertain parameters through the normalized weights of the prior samples. To improve the accuracy of the posterior pdfs, an adaptive procedure is further adopted in the importance sampling for the gradual evolution of its proposal functions. The new proposals assist in fine-tuning the sample space and thereby help to construct surrogates with higher fidelity. Through an iterative process, this approach is able to estimate accurately and efficiently non-Gaussian posterior pdfs of the uncertain generator parameters subject to gross errors.

42 ENGINEERING↗
NASA NTRS1977

Study of photon correlation techniques for processing of laser velocimeter signals

The objective was to provide the theory and a system design for a new type of photon counting processor for low level dual scatter laser velocimeter (LV) signals which would be capable of both the first order measurements of mean flow and turbulence intensity and also the second order time statistics: cross correlation auto correlation, and related spectra. A general Poisson process model for low level LV signals and noise which is valid from the photon-resolved regime all the way to the limiting case of nonstationary Gaussian noise was used. Computer simulation algorithms and higher order statistical moment analysis of Poisson processes were derived and applied to the analysis of photon correlation techniques. A system design using a unique dual correlate and subtract frequency discriminator technique is postulated and analyzed. Expectation analysis indicates that the objective measurements are feasible.

Mayo, W. T., Jr.↗
DOE OSTI2015

PROBABILISTIC MODEL OF BEAM–PLASMA INTERACTION IN RANDOMLY INHOMOGENEOUS PLASMA

We propose a new model that describes beam–plasma interaction in the presence of random density fluctuations with a known probability distribution. We use the property that, for the given frequency, the probability distribution of the density fluctuations uniquely determines the probability distribution of the phase velocity of waves. We present the system as discrete and consisting of small, equal spatial intervals with a linear density profile. This approach allows one to estimate variations in wave energy density and particle velocity, depending on the density gradient on any small spatial interval. Because the characteristic time for the evolution of the electron distribution function and the wave energy is much longer than the time required for a single wave–particle resonant interaction over a small interval, we determine the description for the relaxation process in terms of averaged quantities. We derive a system of equations, similar to the quasi-linear approximation, with the conventional velocity diffusion coefficient D and the wave growth rate γ replaced by the average in phase space, by making use of the probability distribution for phase velocities and by assuming that the interaction in each interval is independent of previous interactions. Functions D and γ are completely determined by the distribution function for the amplitudes of the fluctuations. For the Gaussian distribution of the density fluctuations, we show that the relaxation process is determined by the ratio of beam velocity to plasma thermal velocity, the dispersion of the fluctuations, and the width of the beam in the velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗