Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
NASA NTRS1976

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.↗
NASA NTRS2012

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↗
NASA NTRS2013

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↗
NASA NTRS2004

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.↗
NASA NTRS2006

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.↗
NASA NTRS1972

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.↗
NASA NTRS1982

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.↗
NASA NTRS1986

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.↗
NASA NTRS1987

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.↗
NASA NTRS1987

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.↗
NASA NTRS2008

Determining Concentration of Nanoparticles from Ellipsometry

A method of using ellipsometry or polarization analysis of light in total internal reflection of a surface to determine the number density of gold nanoparticles on a smooth substrate has been developed. The method can be modified to enable determination of densities of sparse distributions of nanoparticles in general, and is expected to be especially useful for measuring gold-nanoparticle-labeled biomolecules on microarrays. The method is based on theoretical calculations of the ellipsometric responses of gold nanoparticles. Elements of the calculations include the following: For simplicity, the gold nanoparticles are assumed to be spherical and to have the same radius. The distribution of gold nanoparticles is assumed to be a sub-monolayer (that is, sparser than a monolayer). The optical response of the sub-monolayer is modeled by use of a thin-island-film theory, according to which the polarizabilities parallel and perpendicular to the substrate are functions of the wavelength of light, the dielectric functions (permittivities expressed as complex functions of frequency or wavelength) of the gold and the suspending medium (in this case, the suspending medium is air), the fraction of the substrate area covered by the nanoparticles, and the radius of the nanoparticles. For the purpose of the thin-island-film theory, the dielectric function of the gold nanoparticles is modeled as the known dielectric function of bulk gold plus a correction term that is necessitated by the fact that the mean free path length for electrons in gold decreases with decreasing radius, in such a manner as to cause the imaginary part of the dielectric function to increase with decreasing radius (see figure). The correction term is a function of the nanoparticle radius, the wavelength of light, the mean free path and the Fermi speed of electrons in bulk gold, the plasma frequency of gold, and the speed of light in a vacuum. These models are used to calculate ellipsometric responses for various concentrations of gold nanoparticles having an assumed radius. The modeled data indicates distinct spectral features for both the real and the imaginary part of the dielectric function. An ellipsometric measurement would determine this distinct feature and thus can be used to measure nanoparticle concentration. By "ellipsometric responses" is meant the intensities of light measured in various polarization states as functions of the angle of incidence and the polarization states of the incident light. These calculated ellipsometric responses are used as calibration curves: Data from subsequent ellipsometric measurements on real specimens are compared with the calibration curves. The concentration of the nanoparticles on a specimen is assumed to be that of the calibration curve that most closely matches the data pertaining to that specimen.

Venkatasubbarao, Srivatsa↗
NASA NTRS1992

Recursive Filtering And Smoothing In Robot Dynamics

Techniques developed originally for electronic systems also useful for multibody mechanical systems. Report summarizes methods developed to solve nonlinear forward-dynamics problem for robot of multiple-link arms connected by joints. Primary objective to show equivalence between recursive methods of dynamical analysis and some filtering and smoothing techniques from state-estimation theory.

Rodriguez, Guillermo↗
NASA NTRS1984

SASS wind ambiguity removal by direct minimization. II - Use of smoothness and dynamical constraints

A variational analysis method (VAM) is used to remove the ambiguity of the Seasat-A Satellite Scatterometer (SASS) winds. The VAM yields the best fit to the data by minimizing an objective function S which is a measure of the lack of fit. The SASS data are described and the function S and the analysis procedure are defined. Analyses of a single ship report which are analogous to Green's functions are presented. The analysis procedure is tuned and its sensitivity is described using the QE II storm. The procedure is then applied to a case study of September 6, 1978, south of Japan.

Hoffman, R. N.↗
NASA NTRS1974

Aerodynamic influence coefficient method using singularity splines

A numerical lifting surface formulation, including computed results for planar wing cases is presented. This formulation, referred to as the vortex spline scheme, combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of loading function methods. The formulation employes a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfied all the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise. Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral. The current formulation uses the elementary horseshoe vortex as the basic singularity and is therefore restricted to linearized potential flow. As part of the study, a non planar development was considered, but the numerical evaluation of the lifting surface concept was restricted to planar configurations. Also, a second order sideslip analysis based on an asymptotic expansion was investigated using the singularity spline formulation.

Mercer, J. E.↗
DOE OSTI2021

Fast & accurate emulation of two-body scattering observables without wave functions

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗