On the Harmonic Balance Method Augmented with Non-Smooth Basis Functions for Contact/Impact Problems
Explore the source record for details and available documents.
Saunders, Brian Evan↗
SEARCH · Engineering Papers
Results for “Smooth Basis Functions”
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.
On the Harmonic Balance Method Augmented with Non-Smooth Basis Functions for Contact/Impact Problems.
Abstract not provided.
Saunders, Brian E↗
Determination of discontinuities in the derivative of a piecewise-smooth function on the basis of the results of measurements
A piecewise-smooth function with discontinuity in the first derivative on a given interval is considered. The values of the function are measured at a sequence of points in the interval and a random error is included in the measurements. A method is proposed to estimate the position of the discontinuity in the derivative. Regression lines are associated with each measurement point and account for k - 1 points preceding or following the point. The estimate for the position of the discontinuity is the measurement point with the largest angle between the regression lines. The error in the estimate is analyzed and the results are verified.
Gurin, L. S.↗
A hybrid meshfree discretization to improve the numerical performance of peridynamic models
Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.
42 ENGINEERING↗
IGA-MPM: The Isogeometric Material Point Method
In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.
42 ENGINEERING↗
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.↗
Exponentially accurate approximations to piece-wise smooth periodic functions
A family of simple, periodic basis functions with 'built-in' discontinuities are introduced, and their properties are analyzed and discussed. Some of their potential usefulness is illustrated in conjunction with the Fourier series representations of functions with discontinuities. In particular, it is demonstrated how they can be used to construct a sequence of approximations which converges exponentially in the maximum norm to a piece-wise smooth function. The theory is illustrated with several examples and the results are discussed in the context of other sequences of functions which can be used to approximate discontinuous functions.
Greer, James↗
Spanning the Gap From Bulk to Bin: A Novel Spectral Microphysics Method
Microphysics methods for climate models and numerical weather prediction typically track one, two, or three moments of a droplet size distribution for various categories of liquid, ice, and aerosol. Such methods rely on conversion parameters between these categories, which introduces uncertainty into predictions. While higher-fidelity options such as bin and Lagrangian schemes exist, they require too many degrees of freedom for climate modeling applications and introduce numerical challenges. Here we introduce a flexible spectral microphysics method based on collocation of basis functions. This method generalizes to a linear bulk scheme when using few basis functions and to a smoothed bin scheme with more degrees of freedom. Tested in an idealized box setting, the method improves spectral accuracy for droplet collision-coalescence and may eliminate the need for precipitation autoconversion rates required by bulk methods; furthermore, it generalizes well to multimodal distributions with less complexity than a bin method. The potential to extend this collocation representation to multiple hydrometeor classes suggests a path forward to unify liquid, ice, and aerosol microphysics in a single, flexible, computational framework for climate modeling.
54 ENVIRONMENTAL SCIENCES↗
A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations
Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.
97 MATHEMATICS AND COMPUTING↗
Basis Selection for Wavelet Regression
A wavelet basis selection procedure is presented for wavelet regression. Both the basis and the threshold are selected using cross-validation. The method includes the capability of incorporating prior knowledge on the smoothness (or shape of the basis functions) into the basis selection procedure. The results of the method are demonstrated on sampled functions widely used in the wavelet regression literature. The results of the method are contrasted with other published methods.
Wheeler, Kevin R.↗
Forward variable selection enables fast and accurate dynamic system identification with Karhunen-Loève decomposed Gaussian processes
A promising approach for scalable Gaussian processes (GPs) is the Karhunen-Loève (KL) decomposition, in which the GP kernel is represented by a set of basis functions which are the eigenfunctions of the kernel operator. Such decomposed kernels have the potential to be very fast, and do not depend on the selection of a reduced set of inducing points. However KL decompositions lead to high dimensionality, and variable selection thus becomes paramount. This paper reports a new method of forward variable selection, enabled by the ordered nature of the basis functions in the KL expansion of the Bayesian Smoothing Spline ANOVA kernel (BSS-ANOVA), coupled with fast Gibbs sampling in a fully Bayesian approach. It quickly and effectively limits the number of terms, yielding a method with competitive accuracies, training and inference times for tabular datasets of low feature set dimensionality. Theoretical computational complexities are O ( N P 2 ) in training and O ( P ) per point in inference, where N is the number of instances and P the number of expansion terms. The inference speed and accuracy makes the method especially useful for dynamic systems identification, by modeling the dynamics in the tangent space as a static problem, then integrating the learned dynamics using a high-order scheme. The methods are demonstrated on two dynamic datasets: a ‘Susceptible, Infected, Recovered’ (SIR) toy problem, along with the experimental ‘Cascaded Tanks’ benchmark dataset. Comparisons on the static prediction of time derivatives are made with a random forest (RF), a residual neural network (ResNet), and the Orthogonal Additive Kernel (OAK) inducing points scalable GP, while for the timeseries prediction comparisons are made with LSTM and GRU recurrent neural networks (RNNs) along with the SINDy package.
Hayes, Kyle↗
Multiresolution Representation Using Biorthogonal Multiwavelets
We generalize Harten's multiresolution representation to biorthogonal multiwavelets. Several variants are considered. For example, a given array of discrete point values is transformed to point values and derivatives or point 'values and cell averages'. Compact Hermite interpolation is used in the decomposition and reconstruction algorithm. The resulting basis functions that are symmetric or skewsymmetric, compact, and smooth with optimal order accuracy. Harten's approach has several advantages: the multiresolution scheme is inherently discrete, non-periodic boundary conditions are easy to implement, and the representation can be extended to unstructured grids in bounded domains. We demonstrate the compression features of the new mutliwavelets by application to variable scale piecewise smooth functions with jump discontinuities typical of numerical solutions of nonlinear hyperbolic conservation laws.
Warming, Robert F.↗
A prescription for the asteroseismic surface correction
In asteroseismology, the surface effect refers to a disparity between the observed and the modelled frequencies in stars with solar-like oscillations. It originates from improper modelling of the surface layers. Correcting the surface effect usually requires using functions with free parameters, which are conventionally fitted to the observed frequencies. On the basis that the correction should vary smoothly across the H–R diagram, we parameterize it as a simple function of surface gravity, effective temperature, and metallicity. We determine this function by fitting a wide range of stars. The absolute amount of the surface correction decreases with luminosity, but the ratio between it and νmax increases, suggesting the surface effect is more important for red giants than dwarfs. Applying the prescription can eliminate unrealistic surface correction, which improves parameter estimations with stellar modelling. Using two open clusters, we found a reduction of scatter in the model-derived ages for each star in the same cluster. As an important application, we provide a new revision for the Δν scaling relation that, for the first time, accounts for the surface correction. The values of the correction factor, f Δν , are up to 2 percent smaller than those determined without the surface effect considered, suggesting decreases of up to 4 percent in radii and up to 8 percent in masses when using the asteroseismic scaling relations. This revision brings the asteroseismic properties into an agreement with those determined from eclipsing binaries. The new correction factor and the stellar models with the corrected frequencies are available at https://www.github.com/parallelpro/surface.
79 ASTRONOMY AND ASTROPHYSICS↗
ARPIST: Provably accurate and stable numerical integration over spherical triangles
Numerical integration on spherical triangles, including the computation of their areas, is a core computation in geomathematics. The commonly used techniques sometimes suffer from instabilities and significant loss of accuracy. We describe a new algorithm, called ARPIST, for accurate and stable integration of functions on spherical triangles. ARPIST is based on an easy-to-implement transformation to the spherical triangle from its corresponding linear triangle via radial projection to achieve high accuracy and efficiency. More importantly, ARPIST overcomes potential instabilities in computing the Jacobian of the transformation, even for poorly shaped triangles that may occur at poles in regular longitude-latitude meshes, by avoiding potential catastrophic rounding errors. We compare our proposed technique with L’Huilier’s Theorem for computing the area of spherical triangles, and also compare it with the recently developed LSQST method (Beckmann et al., 2014) and a radial-basis-function-based technique (Reeger and Fornberg, 2016) for integration of smooth functions on spherical triangulations. In conclusion, our results show that ARPIST enables better or comparable accuracy over previous methods while being easier to implement, significantly faster, and more tolerant of poor element shapes.
97 MATHEMATICS AND COMPUTING↗
Complex Volume Grid Generation Through the Use of Grid Reusability
This paper presents a set of surface and volume grid generation techniques which reuse existing surface and volume grids. These methods use combinations of data manipulations to reduce grid generation time, improve grid characteristics, and increase the capabilities of existing domain discretization software. The manipulation techniques utilize physical and computational domains to produce basis function on which to operate and modify grid character and smooth grids using Trans-Finite Interpolation, a vector interpolation method and parametric re-mapping technique. With these new techniques, inviscid grids can be converted to viscous grids, multiple zone grid adaption can be performed to improve CFD solver efficiency, and topological changes to improve modeling of flow fields can be done simply and quickly. Examples of these capabilities are illustrated as applied to various configurations.
Alter, Stephen J.↗
A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries
Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.
97 MATHEMATICS AND COMPUTING↗
Spacesuit and Mobility Performance Changes
The complex interactions between the human body and spacesuit lead to changes inmovement patterns and mobility performancesof the wearer. In general, factors including the geometric properties, such as shape and size, mechanical properties of the suit, and pressurization of the suit are known to be associated with altered movement patterns as compared to an unsuited human.However, their relative contributions have not been explicitly quantified from the mobility performance perspectives. The goal of this study wasthus to assessthe effects fromthe different types of mobility constraintconditions, namely by wearing either a 3D printed hard upper torso (HUT) assembly orfully pressurized spacesuit. The outcome was also compared against the unsuited motions. For this study, an xEMU (exploration Extravehicular Mobility Unit) suitwas considered, which is the next generation spacesuit developed by NASA. In each test condition, the subject was asked to move the arm and hand as prescribedfor different task types, and the corresponding body segment locations were recordedusing a 3D motion capture system. The following three tasks were performedand analyzed: 1) Outward one-handed reaches:the subject in a standing pose made sweeping motions with the extended right arm from the extreme end-to-end positions, including side-to-side at different elevations and top-to-bottom at different azimuths. 2) Outward two-handed reaches:similarto the previous task, howeverthe subject kept the hands together during the motions in order to assessthe areas that can be reached by both hands. 3) Inward one-handed reaches:the subject made right-hand reach motions to the surface of the HUT.The hand traces collected from each task were modeled by a template shape parametrically deformed with a radial basis function. This process enabled foran abstraction of the hand traces into a smooth surface envelope representing the maximally reachable area of the test subject, of which the shapes and sizes were compared across the different test conditions. The preliminary analysis has shown that theoverall size of reachenvelopes decreases in a pressurized suit compared to 3D printed mockup HUT and unsuited conditions. The specific shape of the envelopes, which were determined by the reachable and unreachable zones, alsovary with the testconditions, and the differences werepronounced with the inward reaches to the HUT surface. The latter observation ispotentially relatedto the increased demandfor shoulder and elbow flexions.Overcomingthe resistance from the pressurizedsoft goods and mechanical constraints of the shoulder assemblywas seen to be associated with the difference in motion patterns between the suited and unsuited conditions. Overall, the information quantified from this study is expected to provide structured metrics for spacesuit mobility, which can improve design optimization and human-system integration.
K Han Kim↗
Molecular and functional analyses of the contractile apparatus in lymphatic muscle
Lymphatics are necessary for the generation and regulation of lymph flow. Lymphatics use phasic contractions and extrinsic compressions to generate flow; tonic contractions alter resistance. Lymphatic muscle exhibits important differences from typical vascular smooth muscle. In this study, the thoracic duct exhibited significant functional differences from mesenteric lymphatics. To understand the molecular basis for these differences, we examined the profiles of contractile proteins and their messages in mesenteric lymphatics, thoracic duct, and arterioles. Results demonstrated that mesenteric lymphatics express only SMB smooth muscle myosin heavy chain (SM-MHC), whereas thoracic duct and arterioles expressed both SMA and SMB isoforms. Both SM1 and SM2 isoforms of SM-MHC were detected in arterioles and mesenteric and thoracic lymphatics. In addition, the fetal cardiac/skeletal slow-twitch muscle-specific beta-MHC message was detected only in mesenteric lymphatics. All four actin messages, cardiac alpha-actin, vascular alpha-actin, enteric gamma-actin, and skeletal alpha-actin, were present in both mesenteric lymphatics and arterioles. However, in thoracic duct, predominantly cardiac alpha-actin and vascular alpha-actin were found. Western blot and immunohistochemical analyses corroborated the mRNA studies. However, in arterioles only vascular alpha-actin protein was detected. These data indicate that lymphatics display genotypic and phenotypic characteristics of vascular, cardiac, and visceral myocytes, which are needed to fulfill the unique roles of the lymphatic system.
Non-NASA Center↗