Smoothing for spectral methods
Explore the source record for details and available documents.
Osher, Stanley↗
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.
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),↗
Evaluating route to impact convergence of the harmonic balance method for piecewise-smooth systems
Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.
42 ENGINEERING↗
New technique for optimal smoothing of data
Recursive method for the optimal smoothing of data has numerical superiority and is more easily understood in terms of physical reasoning than earlier methods. Using a Kalman filter, the smoothing technique, applied to a nonlinear parameter identification problem, is useful in those situations where linearization about a reference solution is valid.
Fraser, D. C.↗
Calculating the grain boundary inclination of voxelated grain structures using a smoothing algorithm
We have developed a flexible method for calculating the grain boundary (GB) inclinations of voxelated grain structure data using smoothing algorithms. We compared the performance of four algorithms: the linear interpolation, Allen–Cahn, level-set, and vertex algorithms. We assessed their accuracy using 2D and 3D cases with known inclinations. The vertex algorithm provided the best balance between accuracy and efficiency for 2D structures while the linear interpolation algorithm provided the best balance for 3D structures. We compared the GB inclinations calculated using our smoothing method on a 3D high energy X-ray diffraction microscopy (HEDM) dataset to those determined by meshing the GBs. The two approaches determined similar GB plane distributions, though they varied significantly at triple junctions. In conclusion, the smoothing method was demonstrated for two sources of 3D voxelated grain structures: HEDM data and results from Monte Carlo Potts grain growth simulations.
36 MATERIALS SCIENCE↗
Theoretical study of some methods for increasing the smoothness of flight through rough air
Explore the source record for details and available documents.
Phillips, William H↗
Score-based deterministic density sampling
We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.
97 MATHEMATICS AND COMPUTING↗
A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing
ABSTRACT Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption-line dominated, emission-line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise ratio (approaching ∼1), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.
79 ASTRONOMY AND ASTROPHYSICS↗
Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries
This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.
42 ENGINEERING↗
Surface smoothing for laser powder-bed Ti-6Al-4V by a transient liquid phase
Surface roughness is the primary driver of fatigue for additively manufactured metals. To address surface roughness, this work introduces a new method to smooth features beyond line-of-sight without material removal. The method applies a coating that triggers local surface remelting by activating a eutectic reaction during heat treatment. The associated liquid phase then wets and isothermally solidifies into a smoother surface. For Ti-6Al-4V fabricated with laser powder bed fusion, samples with and without TLP smoothing (using a Cu coating) were characterized with a suite of techniques, including mechanical testing, electron backscatter diffraction, synchrotron X-ray tomography, and fractography. TLP smoothing reduced surface roughness by 80% and amplified compressive residual stress at the surface by about 50%. With statistically equivalent virtual microstructures, crystal plasticity scrutinized the roles of phases, porosity, and surface roughness. Although the tensile strain-to-failure was reduced to 1% strain, the TLP smoothing process increased high-cycle fatigue strength by about 20% compared to control samples, pointing to future opportunities to optimize the new process through various coating compositions and heat treatment schedules. Overall, this work establishes a new paradigm for treating surfaces of materials for smoothness and compressive residual stress.
Additive manufacturing↗
Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems
The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.
97 MATHEMATICS AND COMPUTING↗
COUPLING SMOOTHED PARTICLE HYDRODYNAMICS WITH FINITE ELEMENT METHOD TO SIMULATE RESIDUAL STRESSES FROM FRICTION STIR PROCESSING
Friction stir processing (FSP) is a solid-state material processing technique that locally modifies the microstructure but also induces undesirable residual stresses. A robust numerical model for the FSP can help in mitigating these residual stresses. Heat source models within a finite element method (FEM) framework suffer from inaccuracies. In contrast, smoothed particle hydrodynamics (SPH) model that explicitly captures the material flow near the tool and the associated heat generation are accurate. However, the computational expense of SPH simulations can be prohibitive. In this work, we propose a coupled SPH-FEM framework. SPH is used to model the heat generation accurately near the tool and which is then inserted into to the FEM model as a heat source. To verify this proposed coupling approach, a test case is set up with typical FSP conditions and it is modeled in both SPH and SPH-FEM. The temperatures profiles were compared after the simulations have reached steady-state temperatures. The similarity of the temperature profiles from SPH-FEA and SPH validated the proposed coupling approach. This proposed approach achieves the accuracy of the SPH method while potentially retaining the low computational expense of FEM.
smoothed particle hydrodynamics, Finite Element Me↗
A bistatic radar method for the determination of epsilon and mu for a smooth spherical target.
Bistatic radar method for determination of permeability and permittivity of material for smooth spherical target
TARGET↗
Solution of axisymmetric and two-dimensional inviscid flow over blunt bodies by the method of lines
Comparisons with experimental data and the results of other computational methods demonstrated that very accurate solutions can be obtained by using relatively few lines with the method of lines approach. This method is semidiscrete and has relatively low core storage requirements as compared with fully discrete methods since very little data were stored across the shock layer. This feature is very attractive for three dimensional problems because it enables computer storage requirements to be reduced by approximately an order of magnitude. In the present study it was found that nine lines was a practical upper limit for two dimensional and axisymmetric problems. This condition limits application of the method to smooth body geometries where relatively few lines would be adequate to describe changes in the flow variables around the body. Extension of the method to three dimensions was conceptually straightforward; however, three dimensional applications would also be limited to smooth body geometries although not necessarily to total of nine lines.
Hamilton, H. H., II↗
Surface definition and grid generation about an Assured Crew Return Vehicle (ACRV) for Space Station Freedom
The surface definition and grid generation about an Assured Crew Return Vehicle (ACRV) for the Space Station Freedom are described. The purpose of the surface definition and grid generation is to provide the necessary geometry information for CFD calculations about the vehicle. There are two salient features in this description. The first is that the numerical model representing the ACRV configuration is obtained from the measurement of an existing wind tunnel model. The method for smoothing the measured data and obtaining the numerical model is described. The second feature is the description of the algebraic grid generation method and software to compute volume grids about the ACRV. The methods and software allow rapid computation of volume grids for a wide range of flow conditions.
Smith, R. E.↗
Sampled-Data Techniques Applied to a Digital Controller for an Altitude Autopilot
Sampled-data theory, using the Z transformation, is applied to the design of a digital controller for an aircraft-altitude autopilot. Particular attention is focused on the sensitivity of the design to parameter variations and the abruptness of the response, that is, the normal acceleration required to carry out a transient maneuver. Consideration of these two characteristics of the system has shown that the finite settling time design method produces an unacceptable system, primarily because of the high sensitivity of the response to parameter variations, although abruptness can be controlled by increasing the sampling period. Also demonstrated is the importance of having well-damped poles or zeros if cancellation is attempted in the design methods. A different method of smoothing the response and obtaining a design which is not excessively sensitive is proposed, and examples are carried through to demonstrate the validity of the procedure. This method is based on design concepts of continuous systems, and it is shown that if no pole-zero cancellations are allowed in the design, one can obtain a response which is not too abrupt, is relatively insensitive to parameter variations, and is not sensitive to practical limits on control-surface rate. This particular design also has the simplest possible pulse transfer function for the digital controller. Simulation techniques and root loci are used for the verification of the design philosophy.
Schmidt, Stanley F.↗
Solutions to Kuessner's integral equation in unsteady flow using local basis functions
The computational procedure and numerical results are presented for a new method to solve Kuessner's integral equation in the case of subsonic compressible flow about harmonically oscillating planar surfaces with controls. Kuessner's equation is a linear transformation from pressure to normalwash. The unknown pressure is expanded in terms of prescribed basis functions and the unknown basis function coefficients are determined in the usual manner by satisfying the given normalwash distribution either collocationally or in the complex least squares sense. The present method of solution differs from previous ones in that the basis functions are defined in a continuous fashion over a relatively small portion of the aerodynamic surface and are zero elsewhere. This method, termed the local basis function method, combines the smoothness and accuracy of distribution methods with the simplicity and versatility of panel methods. Predictions by the local basis function method for unsteady flow are shown to be in excellent agreement with other methods. Also, potential improvements to the present method and extensions to more general classes of solutions are discussed.
Fromme, J. A.↗
Three-dimensional phase field sintering simulations accounting for the rigid-body motion of individual grains
Sintering is a widely used powder processing technique in industrial applications. During sintering, atoms migrate to decrease the energy of the system via two main mechanisms: coarsening and densification, both of which lead to significant morphological variation of the sintered microstructure. When simulating sintering dynamics, the phase-field method has been broadly utilized because of its convenience in tracking morphology evolution. When a large number of grains is involved, it is common to use the same order parameter to describe multiple grains that are not in direct contact with one another (in order to reduce the computational memory demands). However, with this treatment it is difficult to handle the rigid-body motion of individual grains during densification. In this work, an implementation scheme is introduced to overcome the challenge of calculating individual particle motion based on existing equations. It uses a grouping algorithm and sets a cutoff radius on each grain for calculating the particle velocity during densification. This method allows for the incorporation of the densification mechanism, which has been commonly ignored in previous work, into phase-field sintering models in three-dimensional simulations with a large number of particles/grains. Moreover, through combination with the smoothed boundary method, material properties of sintered microstructures, such as the effective diffusivity and Young’s modulus, can be calculated during the sintering processes.
36 MATERIALS SCIENCE↗