Search NASA⌕ Search

SEARCH · Search NASA

Results for “gaussian method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1985

Waveform synthesis of surface waves in a laterally heterogeneous earth by the Gaussian beam method

The present investigation is concerned with an application of the Gaussian beam method to surface waves in the laterally heterogeneous earth. The employed method has been developed for ray tracing and synthesizing seismograms of surface waves in cases involving the laterally heterogeneous earth. The procedure is based on formulations derived by Yomogida (1985). Vertical structure of the wave field is represented by the eigenfunctions of normal mode theory, while lateral variation is expressed by the parabolic equation as in two-dimensional acoustic waves or elastic body waves. It is demonstrated that a large-amplitude change can result from a slight perturbation in the phase velocity model.

Yomogida, K.↗
NASA NTRS1987

Exact image method for Gaussian beam problems involving a planar interface

Exact image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface or two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE polarized beam are correctly given as asymptotic results by the exact reflection image theory. Also, the apparent image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The present theory makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams, as well as to the surface wave.

Lindell, I. V.↗
DOE OSTI2025

Accelerating Instanton Theory with the Line Integral Nudged Elastic Band Method and Gaussian Process Regression

Quantum tunneling plays a fundamental role in many chemical reactions, particularly proton transfer processes. Ring polymer instanton theory offers a practical framework for computing tunneling rates in complex molecular systems. However, applying the ring polymer instanton method with a potential energy surface generated on-the-fly using electronic structure calculations can be computationally demanding. Here, in this work, we present a new efficient implementation of the ring polymer instanton method by combining the Line Integral Nudged Elastic Band (LI-NEB) approach with Gaussian Process Regression (GPR). We benchmarked this method on prototypical ground-state proton transfer systems, including the benchmark gas-phase hydrogen abstraction reaction H + CH 4 → H 2 + CH 3 , malonaldehyde, and Z-3-amino-propenal (aminopropenal). Our results show that this approach is an order of magnitude faster than traditional instanton algorithms while maintaining excellent agreement with their tunneling rates. This development opens the door to studying proton transfer in larger systems with improved efficiency.

chemical physics↗
NASA NTRS1986

A tutorial on the LQG/LTR method

In this paper the so-called Linear-Quadratic-Gaussian method with Loop-Transfer-Recovery is surveyed. The objective is to provide a pragmatic exposition, with special emphasis on the step-by-step characteristics for designing multivariable feedback control systems.

Athans, M.↗
NASA NTRS1972

System identification computational considerations.

System identification methods may be accurately characterized as inverse computational problems. Three approaches to the computational problem are considered. Equation error methods are based on regression analysis techniques. Attention is given to the maximum conditional likelihood estimate and the maximum unconditional likelihood (Bayesian) estimate. Gradient dependent algorithms for solving the minimization problems are discussed, taking into account Gaussian methods, quadratic model methods, and the homogeneous function model method of Jacobson and Oksman. Some advanced methods are also examined.

Bowles, R. L.↗
NASA NTRS1985

Control of an orbiting flexible square platform in the presence of solar radiation

An analysis is made of the dynamics of a thin, homogeneous, flexible square plate exposed to solar radiation disturbance. It is shown that only the first three flexible modes of the plate generate a first order net moment about the center of mass and that the solar radiation pressure does not influence the flexible modes of the plate for small-amplitude vibrations. In the absence of control, for a symmetrical homogeneous square platform, the solar radiation induces a steady angular drift about one of the body principal axes. For extremely flexible platforms, nominally oriented in the local horizontal plane, appreciable rigid modal amplitudes can be induced due to solar radiation, even in the presence of both active and passive control. Optimal control law and feedback gain values are obtained using linear quadratic Gaussian methods.

Bainum, P. M.↗
NASA NTRS1987

Constrained optimization techniques for active control of aeroelastic response

Active control of aeroelastic response is a complex problem in which the designer usually tries to satisfy many design criteria which are often conflicting in nature. To further complicate the design problem, the state space equations describing this type of control problem are usually of high order, involving a large number of states to represent the flexible structure and unsteady aerodynamics. Control laws based on the standard Linear - Quadratic - Gaussian method are of the same high order as the aeroelastic plant and may be difficult to implement in the flight computer. To overcome this disadvantage a new approach was developed for designing low-order optimized robust control laws. In this approach, a nonlinear programming algorithm is used to search for the values of control law design variables that minimize a performance index while satisfying several inequality constraints that describe the design criteria on the stability robustness and responses. The method is applied to a gust load alleviation problem and a stability robustness improvement problem of a drone aircraft.

Mukhopadhyay, Vivekananda↗
NASA NTRS1989

Holographic signal processing for control of a flexible beam

A robust control system is presented for damping vibrations in a cantilevered beam. The distributed deflections of the beam are sensed using real-time holography. A hybrid processor which uses optical and digital processors is proposed to implement a distributed control algorithm. The controller is designed using the linear-quadratic-Gaussian method to add damping to a finite number of vibrational modes.

Barsky, Michael F.↗
NASA NTRS1993

Improving the reconstruction of the velocity potential and primordial density fluctuations by choice of smoothing windows

Recently considerable attention has been focused on improving algorithms for restoring primordial density fluctuations in the universe or investigating large-scale velocity fields by going to higher order approximations relating nonlinear states to their initial conditions. An alternative approach is to investigate the use of alternative smoothing windows within the context of first-order approximations. I present evidence that the universally used Gaussian smoothing window is far from optimal. A sharp truncation P(k) = 0 for k greater than k(c) leads to a much more direct connection to initial conditions for either 'Gaussianization' methods or use of the Zel'dovich approximation, as discussed in the text.

Melott, Adrian L.↗
NASA NTRS2001

Atmospheric Multiple Scattering Effects on GLAS Altimetry: Analysis of Expected Errors in Antarctic Altitude Measurements - Part 2

The altimetry bias in GLAS (Geoscience Laser Altimeter System) or other laser altimeters resulting from atmospheric multiple scattering is studied in relationship to current knowledge of cloud properties over the Antarctic Plateau. Estimates of seasonal and interannual changes in the bias are presented. Results show the bias in altitude from multiple scattering in clouds would be a significant error source without correction. The selective use of low optical depth clouds or cloudfree observations, as well as improved analysis of the return pulse such as by the Gaussian method used here, are necessary to minimize the surface altitude errors. The magnitude of the bias is affected by variations in cloud height, cloud effective particle size and optical depth. Interannual variations in these properties as well as in cloud cover fraction could lead to significant year-to-year variations in the altitude bias. Although cloud-free observations reduce biases in surface elevation measurements from space, over Antarctica these may often include near-surface blowing snow, also a source of scattering-induced delay. With careful selection and analysis of data, laser altimetry specifications can be met.

Mahesh, Ashwin↗
NASA NTRS1964

The Effect of Information Feedback Upon Psychophysical Judgments

An analysis was made of the role of presentation schedules and information feedback on performance in a forced-choice signal detection task. The experimental results indicate that information feedback facilitates performance, but only for certain presentation schedules. The present study was designed to assess performance in a signal detection task under two conditions of information feedback. In the I-condition, S was told on each trial whether his detection response was correct or incorrect; in the !-condition S was given no feedback regarding the correctness of his response. The task involved a 2-response, forced-choice auditory detection problem. On each trial 2 temporal intervals were defined and S was required to report which interval he believed contained the signal; i. e., in one interval a tone burst in a background of white noise was presented, while the other interval contained only white noise. A trial will be denoted as s1 or s2, depending on whether the signal was embedded in the 1st or 2nd interval; the S's response will be denoted A1 or A2 to indicate which interval he reported contained the signal. The probability of an s1 trial will be denoted as y. In this study two values of y were used (.50 and.75) and, as indicated above, two conditions of information feedback. Thus there were 4 experimental conditions (501, · 50I, 751, 75I); each S was run under all 4 conditions. Method Gaussian noise was presented binaurally in S's headphones throughout a test session and the signal was a 1000-cps sinusoid tone; the tone was presented for 100 msec. including equal fall and rise times of 20 msec. The ratio of signal energy to noise power in a unit bandwidth was 2.9, and was constant throughout the study. The. S was seated before a stimulus display board. On each trial a red warning light was flashed for 100 msec. Two amber lights then came on successively each for 1 sec.; these lights defined the 2 observation intervals. The onset of the signal occurred 500 msec. after the onset of one of the observation intervals. After the second amber light went off, S indicated his response by pressing 1 of 2 wand switches under cards reading "1st interval" and "2nd interval." For the !-condition a green light flashed on above the correct response key after S's response; the green light was omitted in the !-condition. Each trial lasted 6 sec. The S's were 12 male college students with normal hearing. They were run for two practice sessions followed by 20 test sessions. Test sessions were run on consecutive days, 350 trials/day. Each day S ran on 1 of the 4 experimental conditions; in successive 4-day blocks S ran one day on each of the 4 experimental conditions in a random order. Thus, over 20 days each of the experimental conditions was repeated 5 times.

Atkinson, Richard C.↗
DOE OSTI2025

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗
NASA NTRS1987

Exact-image method for Gaussian-beam problems involving a planar interface

The exact-image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface between two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian-beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE-polarized beam are correctly given as asymptotic results by the exact-reflection-image theory. Also, the apparent-image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The theory described here makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams as well as to the surface wave.

Lindell, Ismo V.↗
NASA NTRS1987

Probabilistic load model development and validation for composite load spectra for select space propulsion engines

A major task of the program to develop an expert system to predict the loads on selected components of a generic space propulsion engine is the design development and application of a probabilitic loads model. This model is being developed in order to account for the random nature of the loads and assess the variable load ranges' effect on the engine performance. A probabilistic model has been developed. The model is based primarily on simulation methods, but also has a Gaussian algebra method (if all variables are near normal), a fast probability integrator routine (for the calculation of low probability events), and a separate, stand alone program for performing barrier crossing calculations. Each of these probabilistic methods has been verified with theoretical calculations using assumed distributional forms.

Kurth, R.↗
DOE OSTI2025

Development of a forced advection sampling technique (FAST) for quantification of methane emissions from orphaned wells

Abstract. Orphaned wells, meaning unplugged and non-producing wells lacking responsible owners, pose a significant and undersampled environmental challenge due to their vast number and unknown associated emissions. We propose, develop and test an alternative method for estimating emissions from orphaned wells using a forced advection sampling technique (FAST) that can overcome many of the limitations in current methods (cost, accuracy, safety). In contrast to existing ambient Gaussian plume methods, our approach uses a fan-generated flow to force advection between the emission source and a point methane (CH4) sensor. The fan flow field is characterized using a colocated sonic anemometer to measure the 3D wind profile generated by the fan. Using time-series measurements of CH4 concentration and wind, a simple estimate of the CH4 emission rate of the source can be inferred. The method was calibrated using outdoor controlled-release experiments and then tested on four orphaned wells in Lufkin, TX, and Osage County, OK. Our results suggest that the FAST method can provide a low-cost, portable, fast and safe alternative to existing methods with reasonable estimates of orphaned well emissions over a range of leak rates below 40 g h−1 and within certain geometric and atmospheric constraints.

Dubey, Mohit L↗
NASA NTRS2004

Improved Discrete Approximation of Laplacian of Gaussian

An improved method of computing a discrete approximation of the Laplacian of a Gaussian convolution of an image has been devised. The primary advantage of the method is that without substantially degrading the accuracy of the end result, it reduces the amount of information that must be processed and thus reduces the amount of circuitry needed to perform the Laplacian-of- Gaussian (LOG) operation. Some background information is necessary to place the method in context. The method is intended for application to the LOG part of a process of real-time digital filtering of digitized video data that represent brightnesses in pixels in a square array. The particular filtering process of interest is one that converts pixel brightnesses to binary form, thereby reducing the amount of information that must be performed in subsequent correlation processing (e.g., correlations between images in a stereoscopic pair for determining distances or correlations between successive frames of the same image for detecting motions). The Laplacian is often included in the filtering process because it emphasizes edges and textures, while the Gaussian is often included because it smooths out noise that might not be consistent between left and right images or between successive frames of the same image.

Shuler, Robert L., Jr.↗
NASA NTRS2000

Higher Order Modeling In the BEM/FEM Hybrid Formulation

Hybrid formulations using low order curl-conforming bases to represent the total electric field within a finite element region and low order divergence-conforming bases to represent equivalent electric and magnetic currents on the boundary are well known. However, higher-order divergence and curl-conforming bases have been shown to provide significant benefits in convergence rates and accuracy when employed in strictly integral equation and strictly finite element formulations. In this paper, a hybrid electric field formulation employing higher order bases is presented. The paper addresses benefits and issues associated with using higher order divergence-and curl-conforming bases in the hybrid finite element/boundary element electric field formulation. The method of singularity subtraction may be used to compute the self terms of the boundary integral when the bases are of low order. But this method becomes laborious and requires great care when the divergence conforming bases are of higher order. In order to handle these singularities simply and accurately, a generalized Gaussian quadrature method is employed in which the expansion functions account for the singularity. In preliminary tests of the higher order hybrid formulation, the equivalent electric current induced by scattering of a plane wave from a square dielectric cylinder is examined. Accurate results are obtained using only a two-triangle mesh when the current basis is of order 4 or 5. Additional results are presented comparing the error obtained using higher order bases to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained with higher order bases are compared to those obtained with lower order bases for selected sample problems.

Fink, Patrick W.↗