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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Smoothness of Functions in Spaces of the Finite Element Method

We find necessary and sufficient conditions for the generalized smoothness of the coordinate functions obtained from the approximate relations. We show that for the coordinate functions the smoothness on their supports is equivalent to that on the boundaries of supports. We obtain conditions for the continuity of Courant type finite element approximations and conditions for the uniqueness of a linear space of such approximations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗
NASA NTRS2007

Smooth Phase Interpolated Keying

Smooth phase interpolated keying (SPIK) is an improved method of computing smooth phase-modulation waveforms for radio communication systems that convey digital information. SPIK is applicable to a variety of phase-shift-keying (PSK) modulation schemes, including quaternary PSK (QPSK), octonary PSK (8PSK), and 16PSK. In comparison with a related prior method, SPIK offers advantages of better performance and less complexity of implementation. In a PSK scheme, the underlying information waveform that one seeks to convey consists of discrete rectangular steps, but the spectral width of such a waveform is excessive for practical radio communication. Therefore, the problem is to smooth the step phase waveform in such a manner as to maintain power and bandwidth efficiency without incurring an unacceptably large error rate and without introducing undesired variations in the amplitude of the affected radio signal. Although the ideal constellation of PSK phasor points does not cause amplitude variations, filtering of the modulation waveform (in which, typically, a rectangular pulse is converted to a square-root raised cosine pulse) causes amplitude fluctuations. If a power-efficient nonlinear amplifier is used in the radio communication system, the fluctuating-amplitude signal can undergo significant spectral regrowth, thus compromising the bandwidth efficiency of the system. In the related prior method, one seeks to solve the problem in a procedure that comprises two major steps: phase-value generation and phase interpolation. SPIK follows the two-step approach of the related prior method, but the details of the steps are different. In the phase-value-generation step, the phase values of symbols in the PSK constellation are determined by a phase function that is said to be maximally smooth and that is chosen to minimize the spectral spread of the modulated signal. In this step, the constellation is divided into two groups by assigning, to information symbols, phase values that result in equal numbers of clockwise and counter-clockwise phase rotations for equally likely symbols. The purpose served by assigning phase values in this way is to prevent unnecessary generation of spectral lines and prevent net shifts of the carrier signal. In the phase-interpolation step, the smooth phase values are interpolated over a number, n, of consecutive symbols (including the present symbol) by means of an unconventional spline curve fit.

Borah, Deva K.↗
DOE OSTI2023

BeyondPlanck: VII. Bayesian estimation of gain and absolute calibration for cosmic microwave background experiments

We present a Bayesian calibration algorithm for cosmic microwave background (CMB) observations as implemented within the global end-to-end BEYONDPLANCK framework and applied to the Planck Low Frequency Instrument (LFI) data. Following the most recent Planck analysis, we decomposed the full time-dependent gain into a sum of three nearly orthogonal components: one absolute calibration term, common to all detectors, one time-independent term that can vary between detectors, and one time-dependent component that was allowed to vary between one-hour pointing periods. Each term was then sampled conditionally on all other parameters in the global signal model through Gibbs sampling. The absolute calibration is sampled using only the orbital dipole as a reference source, while the two relative gain components were sampled using the full sky signal, including the orbital and Solar CMB dipoles, CMB fluctuations, and foreground contributions. We discuss various aspects of the data that influence gain estimation, including the dipole-polarization quadrupole degeneracy and processing masks. Comparing our solution to previous pipelines, we find good agreement in general, with relative deviations of -0.67% (-0.84%) for 30 GHz, 0.12% (-0.04%) for 44 GHz and -0.03% (-0.64%) for 70 GHz, compared to Planck PR4 and Planck 2018, respectively. We note that the BEYONDPLANCK calibration was performed globally, which results in better inter-frequency consistency than previous estimates. Additionally, WMAP observations were used actively in the BEYONDPLANCK analysis, which both breaks internal degeneracies in the Planck data set and results in an overall better agreement with WMAP. Finally, we used a Wiener filtering approach to smoothing the gain estimates. We show that this method avoids artifacts in the correlated noise maps as a result of oversmoothing the gain solution, which is difficult to avoid with methods like boxcar smoothing, as Wiener filtering by construction maintains a balance between data fidelity and prior knowledge. Although our presentation and algorithm are currently oriented toward LFI processing, the general procedure is fully generalizable to other experiments, as long as the Solar dipole signal is available to be used for calibration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
NASA NTRS1979

Acceleration of transonic potential flow calculations on arbitrary meshes by the multiple grid method

A multiple grid method for transonic flow calculations is developed. The proposed scheme incorporates a generalized alternating direction method as the smoothing algorithm. Numerical experiments indicate that this multigrid alternating direction method converges rapidly and reliably for a range of cases typical of the cruising regime up to the onset of drag rise. It also appears that the method can be readily generalized to treat three-dimensional flows.

Jameson, A.↗
DOE OSTI2020

Topological Regularization via Persistence-Sensitive Optimization

Optimization, a key tool in machine learning and statistics, relies on regularization to reduce overfitting. Traditional regularization methods control a norm of the solution to ensure its smoothness. Recently, topological methods have emerged as a way to provide a more precise and expressive control over the solution, relying on persistent homology to quantify and reduce its roughness. All such existing techniques back-propagate gradients through the persistence diagram, which is a summary of the topological features of a function. Their downside is that they provide information only at the critical points of the function. We propose a method that instead builds on persistence-sensitive simplification and translates the required changes to the persistence diagram into changes on large subsets of the domain, including both critical and regular points. This approach enables a faster and more precise topological regularization, the benefits of which we illustrate with experimental evidence.

Nigmetov, Arnur↗
NASA NTRS1993

A variational method for finite element stress recovery and error estimation

A variational method for obtaining smoothed stresses from a finite element derived nonsmooth stress field is presented. The method is based on minimizing a functional involving discrete least-squares error plus a penalty constraint that ensures smoothness of the stress field. An equivalent accuracy criterion is developed for the smoothing analysis which results in a C sup 1-continuous smoothed stress field possessing the same order of accuracy as that found at the superconvergent optimal stress points of the original finite element analysis. Application of the smoothing analysis to residual error estimation is also demonstrated.

Tessler, A.↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗
DOE OSTI2020

Impact of isolation and fiducial cuts on qT and N-jettiness subtractions

Kinematic selection cuts and isolation requirements are a necessity in experimental measurements for identifying prompt leptons and photons that originate from the hard-interaction process of interest. We analyze how such cuts affect the application of the q T and N-jettiness subtraction methods for fixed-order calculations. We consider both fixed-cone and smooth-cone isolation methods. We find that kinematic selection and isolation cuts both induce parametrically enhanced power corrections with considerably slower convergence compared to the standard power corrections that are already present in inclusive cross sections without additional cuts. Using analytic arguments at next-to-leading order we derive their general scaling behavior as a function of the subtraction cutoff. We also study their numerical impact for the case of gluon-fusion Higgs production in the H → γγdecay mode and for pp → γγ direct diphoton production. We find that the relative enhancement of the additional cut-induced power corrections tends to be more severe forq T , where it can reach an order of magnitude or more, depending on the choice of parameters and subtraction cutoffs. We discuss how all such cuts can be incorporated without causing additional power corrections by implementing the subtractions differentially rather than through a global slicing method. We also highlight the close relation of this formulation of the subtractions to the projection-to-Born method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Analysis of a finite composite plate with smooth rigid pin

An elastic continuum method of analysis developed for a smooth rigid pin in a finite composite plate by using complex potentials that satisfy the governing orthotropic differential equation is presented. The plate is subjected to uniaxial loading and the nonlinear problems of partial contact in the pin-plate are studied. Three types of pin fits are examined: interference, push, and clearance. The problems of stresses due to pure interference, load contact behavior, and stresses variations are addressed. The effects of the rectangular geometry and the orthotropic property of the plate are demonstrated through numerical results. It is found that the use of large hole diameter/plate width ratios have deleterious effects and the use of interference has beneficial effects. The solutions obtained will help in the development of special finite elements for fastener joints. Finally, the boundary conditions on the finite outer boundaries are found to be satisfied by a successive integration technique.

Naidu, A. C. B.↗
DOE OSTI2020

Embedded smoothed particle hydrodynamics

In this paper, a hybrid method for embedding a meshfree smoothed particle hydrodynamics into mesh-based ALE hydro schemes is presented. The method relies on the embedded mesh framework presented in Puso et al. (2015, 2012) and utilizes SPH (Owen, 2010) in the solid foreground; when an appropriate damage or failure model is used with SPH, it becomes a natural way to model material fracture and cracking in high-explosive driven experiments. Overall, we provide various examples which validate the coupling between SPH and the embedded mesh method.

42 ENGINEERING↗
NASA NTRS1989

The 3DGRAPE book: Theory, users' manual, examples

A users' manual for a new three-dimensional grid generator called 3DGRAPE is presented. The program, written in FORTRAN, is capable of making zonal (blocked) computational grids in or about almost any shape. Grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. The smoothness for which elliptic methods are known is seen here, including smoothness across zonal boundaries. An introduction giving the history, motivation, capabilities, and philosophy of 3DGRAPE is presented first. Then follows a chapter on the program itself. The input is then described in detail. A chapter on reading the output and debugging follows. Three examples are then described, including sample input data and plots of output. Last is a chapter on the theoretical development of the method.

Sorenson, Reese L.↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
NASA NTRS1969

Towards an optimal inversion method for remote atmospheric sensing

The inference of atmospheric structure from satellite radiometric observations requires an inversion algorithm. A variety of techniques was spawned to meet these demands. One class, the nonlinear inversion methods, copes with the problem of data noise. Unlike linear techniques which require a priori data smoothing, the nonlinear method can be applied directly to raw data. The algorithm discriminates the noise input by resolving the inferences into two types of solution, associating the real roots with atmospheric structure while ascribing the imaginary roots to noise.

King, J. I. F.↗
DOE OSTI2024

Autonomous Electrochemistry Platform with Real-Time Normality Testing of Voltammetry Measurements Using ML

Electrochemistry workflows utilize various instruments and computing systems to execute workflows consisting of electrocatalyst synthesis, testing and evaluation tasks. The heterogeneity of the software and hardware of these ecosystems makes it challenging to orchestrate a complete workflow from production to characterization by automating its tasks. We propose an autonomous electrochemistry computing platform for a multi-site ecosystem that provides the services for remote experiment steering, real-time measurement transfer, and AI/ML-driven analytics. We describe the integration of a mobile robot and synthesis workstation into the ecosystem by developing custom hub-networks and software modules to support remote operations over the ecosystem’s wireless and wired networks. We describe a workflow task for generating I-V voltammetry measurements using a potentiostat, and a machine learning framework to ensure their normality by detecting abnormal conditions such as disconnected electrodes. We study a number of machine learning methods for the underlying detection problem, including smooth, non-smooth, structural and statistical methods, and their fusers. We present experimental results to illustrate the effectiveness of this platform, and also validate the proposed ML method by deriving its rigorous generalization equations.

Alnajjar, Anees↗
NASA NTRS1994

Intercomparison of single-frequency methods for retrieving a vertical rain profile from airborne or spaceborne radar data

This paper briefly reviews several single-frequency rain profiling methods for an airborne or spaceborne radar. The authors describe the different methods from a unified point of view starting from the basic differential equation. This facilitates the comparisons between the methods and also provides a better understanding of the physical and mathematical basis of the methods. The application of several methods to airborne radar data taken during the Convective and Precipitation/Electrification Experiment is shown. Finally, the authors consider a hybrid method that provides a smooth transition between the Hitschfeld-Bordan method, which performs well at low attenuations, and the surface reference method, for which the relative error decreases with increasing path attenuation.

Iguchi, Toshio↗