Search NASA⌕ Search

SEARCH · Search NASA

Results for “smoothing methods”

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 271 records · Page 15

Absolute electron density fluctuation reconstruction for two-dimensional hydrogen beam emission spectroscopy

Scrape-off layer (SOL) and edge plasma turbulence significantly contribute to the radial particle and heat transport, lowering the plasma confinement and increasing the heat load on the plasma facing components. SOL turbulence is predominantly intermittent, which manifests in the occurrence of isolated density filaments or blobs. Filaments propagate radially outward toward plasma facing components, limiting their lifetime by erosion and sputtering. To characterize this phenomenon in detail, few diagnostic techniques are available. Beam emission spectroscopy is a diagnostic capable of measuring plasma turbulence in both SOL and edge plasmas. Due to the finite lifetime of the excitation states during the beam–plasma interaction and the misalignment between the optics and the magnetic field, spatial smearing is introduced in the measurement. In this paper, a novel method is introduced to overcome this hindering effect by inverting the fluctuation response matrix on an optimally smoothed signal. We show that this method is fast and provides significantly more accurate absolute density fluctuation reconstruction than the direct inversion technique. Here, the presented method is usable for all types of beam emission diagnostics where the spatial resolution is higher than the combined smearing of the atomic physics and the observation.

47 OTHER INSTRUMENTATION↗

Comparison of Image Restoration Methods for Lunar Epithermal Neutron Emission Mapping

Orbital measurements of neutrons by the Lunar Exploring Neutron Detector (LEND) onboard the Lunar Reconnaissance Orbiter are being used to quantify the spatial distribution of near surface hydrogen (H). Inferred H concentration maps have low signal-to-noise (SN) and image restoration (IR) techniques are being studied to enhance results. A single-blind. two-phase study is described in which four teams of researchers independently developed image restoration techniques optimized for LEND data. Synthetic lunar epithermal neutron emission maps were derived from LEND simulations. These data were used as ground truth to determine the relative quantitative performance of the IR methods vs. a default denoising (smoothing) technique. We review and used factors influencing orbital remote sensing of neutrons emitted from the lunar surface to develop a database of synthetic "true" maps for performance evaluation. A prior independent training phase was implemented for each technique to assure methods were optimized before the blind trial. Method performance was determined using several regional root-mean-square error metrics specific to epithermal signals of interest. Results indicate unbiased IR methods realize only small signal gains in most of the tested metrics. This suggests other physically based modeling assumptions are required to produce appreciable signal gains in similar low SN IR applications.

McClanahan, T. P.↗

A variable high-order shock-capturing finite difference method with GP-WENO

Here, we present a new finite difference shock-capturing scheme for hyperbolic equations on static uniform grids. The method provides selectable high-order accuracy by employing a kernel-based Gaussian Process (GP) data prediction method which is an extension of the GP high-order method originally introduced in a finite volume framework by the same authors. The method interpolates Riemann states to high order, replacing the conventional polynomial interpolations with polynomial-free GP-based interpolations. For shocks and discontinuities, this GP interpolation scheme uses a nonlinear shock handling strategy similar to Weighted Essentially Non-Oscillatory (WENO), with a novelty consisting in the fact that nonlinear smoothness indicators are formulated in terms of the Gaussian likelihood of the local stencil data, replacing the conventional L 2 -type smoothness indicators of the original WENO method. We demonstrate that these GP-based smoothness indicators play a key role in the new algorithm, providing significant improvements in delivering high – and selectable – order accuracy in smooth flows, while successfully delivering non-oscillatory solution behavior in discontinuous flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An example of scaling MST Doppler spectra using median spectra, spectral smoothing, and velocity tracing

Although automatic, computer scaling methods appeared at the start of the MST (mesosphere stratosphere troposphere) radar technique, there is a continuing need for scaling algorithms that perform editing functions and increase the sensitivity of radar by post processing. The scaling method presented is an adaptation of the method of scaling MST Doppler spectra presented by Rastogi (1984). A brief overview of this method is as follows: a median spectrum is calculated from several sequential spectra; the median noise value is subtracted from this derived spectrum; the median spectrum is smoothed; the detection/nondetection decision is made by comparing the smoothed spectrum to the variance of the smoothed noise; and if a signal is detected, then the half-power points of the smoothed echo spectrum are used to place limits on the evaluation of the first two moments of the unsmoothed median spectrum. In all of the above steps, the algorithm is guided by tracing the expected velocity range upward from the lowest range as far as possible. The method is discussed in more detail.

Green, J. L.↗

Interpolation Methods for GTD Analysis of Shaped Reflectors

The finding of smooth analytic representations for antenna reflector surfaces which are prescribed only by discretized data obtained by various synthesis methods is examined. Frequently the data are distributed in a nonuniform grid and contain noise. The smoothness required is to C sub 1 for physical optics diffraction analysis and to C sub 2 for geometrical theory of diffraction (GTD) analysis. The GTD analysis approach requires a surface description which returns data very rapidly. Two methods of interpolation, the global and the local methods, are discussed. They each have advantages and disadvantages; characteristics are discussed and examples are presented.

Galindo-Israel, V.↗

Development of an automatic block generation algorithm

A method for automatic multiblock grid generation is described. The method combines the modified advancing front method as a predictor with an elliptic scheme as a corrector. It advances a collection of cells by one cell height in the outward direction using modified advancing front method, and then corrects newly-obtained cell positions by solving elliptic equations. This predictor-corrector type scheme is repeatedly applied until the field of interest is filled with hexahedral grid cells. Given the configuration surface grid, the scheme produces block layouts as well as grid cells with overall smoothness as its output. The method saves human-time and reduces the burden on the user in generating grids for general 3D configurations. It is used to generate multiblock grids for wings in their high-lift configuration.

Eberhardt, Scott↗

Neptune radio emission in dipole and multipole magnetic fields

We study Neptune's smooth radio emission in two ways: we simulate the observations and we then consider the radio effects of Neptune's magnetic multipoles. A procedure to deduce the characteristics of radio sources observed by the Planetary Radio Astronomy experiment minimizes limiting assumptions and maximizes use of the data, including quantitative measurement of circular polarization. Study of specific sources simulates time variation of intensity and apparent polarization of their integrated emission over an extended time period. The method is applied to Neptune smooth recurrent emission (SRE). Time series are modeled with both broad and beamed emission patterns, and at two frequencies which exhibit different time variation of polarization. These dipole-based results are overturned by consideration of more complex models of Neptune's magnetic field. Any smooth emission from the anticipated auroral radio source is weak and briefly observed. Dominant SRE originates complex fields at midlatitude. Possible SRE source locations overlap that of 'high-latitude' emission (HLE) between +(out) and -(in) quadrupoles. This is the first identification of multipolar magnetic structure with a major source of planetary radio emission.

Sawyer, C. B.↗

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge–Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6–0.7 bohr in OPAW-TDDFT rather than the 0.4–0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe–Salpeter equation.

Chemistry↗

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse to fine. The coarser grids are used to diminish the magnitude of the smooth part of the residuals. The method was applied to the solution of the transonic small disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is studied with meshes of both constant and variable step size.

South, J. C., Jr.↗

Q-Method Extended Kalman Filter

A new algorithm is proposed that smoothly integrates non-linear estimation of the attitude quaternion using Davenport s q-method and estimation of non-attitude states through an extended Kalman filter. The new method is compared to a similar existing algorithm showing its similarities and differences. The validity of the proposed approach is confirmed through numerical simulations.

Zanetti, Renato↗

Characterizing Wheel-Soil Interaction Loads Using Meshfree Finite Element Methods: A Sensitivity Analysis for Design Trade Studies

A wheel experiencing sinkage and slippage events poses a high risk to planetary rover missions as evidenced by the mobility challenges endured by the Mars Exploration Rover (MER) project. Current wheel design practice utilizes loads derived from a series of events in the life cycle of the rover which do not include (1) failure metrics related to wheel sinkage and slippage and (2) performance trade-offs based on grouser placement/orientation. Wheel designs are rigorously tested experimentally through a variety of drive scenarios and simulated soil environments; however, a robust simulation capability is still in development due to myriad of complex interaction phenomena that contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree nite element approaches enable simulations that capture su cient detail of wheel-soil interaction while remaining computationally feasible. This study implements the JPL wheel-soil benchmark problem in the commercial code environment utilizing the large deformation modeling capability of Smooth Particle Hydrodynamics (SPH) meshfree methods. The nominal, benchmark wheel-soil interaction model that produces numerically stable and physically realistic results is presented and simulations are shown for both wheel traverse and wheel sinkage cases. A sensitivity analysis developing the capability and framework for future ight applications is conducted to illustrate the importance of perturbations to critical material properties and parameters. Implementation of the proposed soil-wheel interaction simulation capability and associated sensitivity framework has the potential to reduce experimentation cost and improve the early stage wheel design proce

terramechanics↗

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reductions in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

Asymmetric hypersonic flow

A general method for the analysis of the inviscid asymmetric hypersonic flow fields enveloping smooth bodies of general shape is given. The method is based on the assumption of a thin shock layer which yields an expression for pressure in generalized Mises coordinates. Numerical results for elliptic cones at angle of attack are shown to compare well with experiments and other theories. The computing logic for a blunt body is described, and a limiting solution at the stagnation point is presented.

Maslen, S. H.↗

Feasibility of SPF/DB titanium sandwich for LFC wings

The feasibility of fabricating SPF/DB titanium structures of sufficient smoothness to be used for laminar flow wing surfaces was demonstrated. Two methods of fabricating panels which meet the surface smoothness criteria were demonstrated. The first consists of superplasticity forming/diffusion bonding a panel using steel dies, then machining the surface to the required flatness and finish after forming. This approach, however, was estimated to be more costly than the second approach, in which the panel is formed against ceramic platens which produce the desired surface smoothness without subsequent finishing. Four panels were fabricated, and their surface quality was evaluated. The specified criteria for waviness compared with actual measured values on the panels is shown. The acceptable surface quality as well as feasibility of the laminar flow control surface design, in which separate strips incorporating the boundary layer bleed provisions are bonded into slots on the surface, was also demonstrated.

Mcquilkin, F. T.↗

The method of complex characteristics for design of transonic blade sections

A variety of computational methods were developed to obtain shockless or near shockless flow past two-dimensional airfoils. The approach used was the method of complex characteristics, which determines smooth solutions to the transonic flow equations based on an input speed distribution. General results from fluid mechanics are presented. An account of the method of complex characteristics is given including a description of the particular spaces and coordinates, conformal transformations, and numerical procedures that are used. The operation of the computer program COMPRES is presented along with examples of blade sections designed with the code. A user manual is included with a glossary to provide additional information which may be helpful. The computer program in Fortran, including numerous comment cards is listed.

Bledsoe, M. R.↗

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗