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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Prediction of sonic boom at a focus

The behavior of sonic boom at a focus has been reviewed for the purpose of extending present sonic boom computational methods to include focal zones. The geometry of a focal zone - whether a smooth caustic, a cusped caustic, or a perfect focus to a point - determines the character of focused signatures. The seeming contradiction of various experimental data can be resolved by noting these differences. A ray acoustic analysis has been developed for quantitative determination of caustic geometry. The only reliable theory presently available for signatures at a focus is for a smooth caustic. There has been some controversy between theoretical and experimental values of a constant in the scaling law for this case. It has been found that this discrepancy can be resolved by accounting for the finite thickness of real sonic boom shock waves. These findings have been incorporated into an existing sonic boom computer program.

Plotkin, K. J.↗
NASA NTRS1980

Design of slotted transonic wind tunnels for supersonic flow development

A simple analytical method to predict the supersonic flow development in a two-dimensional slotted transonic wind tunnel has been developed and validated. While doing that, an insight has been gained into the flow mechanism which causes overexpansion with some slot shapes. As a consequence of the success of the analysis method, a design method has been developed on similar lines, to obtain slot shapes for prescribed smooth supersonic flow development at a design Mach number. The performance of these slot shapes at off-design conditions have been studied to enable choice of slot shapes to be made, based on acceptable performance over the whole supersonic range of interest.

Ramaswamy, M. A.↗
NASA NTRS1984

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗
NASA NTRS1987

Color graphics techniques for shaded surface displays of aerodynamic flowfield parameters

A variety of algorithms are presented that are useful in the color graphics displays of large sets of three-dimensional aerodynamic data. These methods augment established techniques for rapidly and interactively reviewing an assortment of smooth-shaded, color-coded illustrations for experimental and computational distributions of both scalar and vector quantities.

Weston, Robert P.↗
NASA NTRS1988

Spectral element multigrid. Part 2: Theoretical justification

A multigrid algorithm is analyzed which is used for solving iteratively the algebraic system resulting from tha approximation of a second order problem by spectral or spectral element methods. The analysis, performed here in the one dimensional case, justifies the good smoothing properties of the Jacobi preconditioner that was presented in Part 1 of this paper.

Maday, Yvon↗
NASA NTRS1989

SAGE: A 2-D self-adaptive grid evolution code and its application in computational fluid dynamics

SAGE is a user-friendly, highly efficient, two-dimensional self-adaptive grid code based on Nakahashi and Deiwert's variational principles method. Grid points are redistributed into regions of high flowfield gradients while maintaining smoothness and orthogonality of the grid. Efficiency is obtained by splitting the adaption into 2 directions and applying one-sided torsion control, thus producing a 1-D elliptic system that can be solved as a set of tridiagonal equations.

Davies, Carol B.↗
NASA NTRS1993

Integration of Design, Thermal, Structural, and Optical Analysis, Including Thermal Animation

In many industries there has recently been a concerted movement toward 'quality management' and the issue of how to accomplish work more efficiently. Part of this effort is focused on concurrent engineering; the idea of integrating the design and analysis processes so that they are not separate, sequential processes (often involving design rework due to analytical findings) but instead form an integrated system with smooth transfers of information. Presented herein are several specific examples of concurrent engineering methods being carried out at Langley Research Center (LaRC): integration of thermal, structural and optical analyses to predict changes in optical performance based on thermal and structural effects; integration of the CAD design process with thermal and structural analyses; and integration of analysis and presentation by animating the thermal response of a system as an active color map -- a highly effective visual indication of heat flow.

Amundsen, Ruth M.↗
NASA NTRSDate not supplied

Single-Scattering Properties of Melting Precipitation for GPM Passive Microwave and Radar Remote Sensing Applications

Over the past two decades, detailed computational simulations of the intricate three dimensional structures of ice-phase crystals and aggregates of those crystals have been developed. The microwave single-scattering properties of these simulated ice particles have been computed and used to improve quantitative estimates of snow rates and to better deRne the vertical structure of snow water contents in deep convective systems, as derived from satellite-borne passive microwave and/or radar remote sensing measurements. The same icephase particles have more recently been used as the starting point for simulations of melting precipitation using computational melting methods. In the current study, a heuristic melting method, as well as a physically-based melting procedure based on smoothed-particle hydrodynamics, are applied to ice particle models to describe the full evolution of the particles from dry snow to liquid drops. The discrete dipole approximation is utilized to calculate the single-scattering properties of the mixed-phase particles throughout the melting process. Then, the properties of the particles are “mapped” into simpliRed microphysical simulations of particle size spectra in the melting layers of stratiform, precipitating clouds. The bulk single-scattering properties of the melting layers and the sensitivity of their properties to modeling assumptions are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗
DOE OSTI2021

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗
NASA NTRS1994

A spectral boundary integral equation method for the 2-D Helmholtz equation

In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the nonsmoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the nonsmoothness of the integral kernels in the spectral implementation. The present method is robust for a general boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. A numerical example of wave scattering is given in which the exponential accuracy of the present numerical method is demonstrated.

Hu, Fang Q.↗
DOE OSTI2018

An Evaluation of Bang-Time Measurements from a Multichannel, Triaxial, nTOF Diagnostic for MagLIF Experiments at the Z facility

Neutron bang times for a series of MagLIF (Magnetic Liner Inertial Fusion) experiments with D 2 -filled targets have been measured at the Z facility. The emitted neutrons were detected as current-mode pulses in a multichannel, neutron time-of-flight (nTOF) diagnostic with conventional, scintillator-photomultiplier-tube (PMT) detectors. In these experiments, the detectors were fielded at known, fixed distances L (690-2510 cm) from the target, and on three, non-coplanar (but convergent) lines-of-sight (LOS). The primary goal of this diagnostic was to estimate a fiducial time (bang time) relative to an externally generated time-base for synchronizing all the diagnostics in an experiment. Recorded arrival times (A7) of the pulses were characterized experimentally by three numerical methods: a first-moment estimate (centroid) and two nodal measures — Savitzky-Golay (SG) smoothing and a single point peak estimate of the raw data. These times were corrected for internal detector time delays (transit and impulse-response function) — an adjustment that linked the recorded ATs to the corresponding arrival of uncollided neutrons at each detector. The bang time was then estimated by linearly regressing the arrival times against the associated distances to the source; t bang (on the system timescale) was taken as the temporal intercept of the regression equation at distance L = 0. This article reports the analysis for a representative shot #2584 for which (a) the recorded ATs — even without detector corrections — agreed by method in each channel to within 1-2 ns; (b) internal corrections were each ~3 — 5 ns; and (c) a 95% uncertainty (confidence) interval for t bang in this shot was estimated at ±3 ns with 4 degrees of freedom. A secondary goal for this diagnostic was to check that the bang time measurements corresponded to neutrons emitted by the D(d,n) 3 He reaction in a thermalized DD plasma. According to the theoretical studies by Brysk, such neutrons should be emitted with an isotropic Gaussian distribution of mean kinetic energy $ \overline{E}$ of 2.449 MeV; this energy translates to a mean neutron speed $ \overline{u}$ of 2.160 cm/ns [D. H. Munro, Nuclear Fusion, 56(3) 036001 (2016)]. In the MagLIF series of shots there was no evidence of spatial asymmetry in the time-distance regressions, and it was possible to extract the mean neutron speed from the slope of these fits. In shot 2584 $ \overline{u}$ was estimated at 2.152 cm/ns ± 0.010 cm/ns [95 % confidence, 4 dof] and the mean kinetic energy $ \overline{E}$ (with relativistic corrections) was 2.431 MeV ± 0.022 MeV [95 % confidence, 4 dof] — results supporting the assumption that D-D neutrons were, in fact, measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Variational projection methods for gradient crystal plasticity using Lie algebras

Summary A computational method is developed for evaluating the plastic strain gradient hardening term within a crystal plasticity formulation. While such gradient terms reproduce the size effects exhibited in experiments, incorporating derivatives of the plastic strain yields a nonlocal constitutive model. Rather than applying mixed methods, we propose an alternative method whereby the plastic deformation gradient is variationally projected from the elemental integration points onto a smoothed nodal field. Crucially, the projection utilizes the mapping between Lie groups and algebras in order to preserve essential physical properties, such as orthogonality of the plastic rotation tensor. Following the projection, the plastic strain field is directly differentiated to yield the Nye tensor. Additionally, an augmentation scheme is introduced within the global Newton iteration loop such that the computed Nye tensor field is fed back into the stress update procedure. Effectively, this method results in a fully implicit evolution of the constitutive model within a traditional displacement‐based formulation. An elemental projection method with explicit time integration of the plastic rotation tensor is compared as a reference. A series of numerical tests are performed for several element types in order to assess the robustness of the method, with emphasis placed upon polycrystalline domains and multi‐axis loading. Copyright © 2016 John Wiley & Sons, Ltd.

Truster, Timothy J.↗
DOE OSTI2015

Some new function spaces of variable smoothness

A new Besov space of variable smoothness is introduced on which the norm is defined in terms of difference relations. This space is shown to be the trace of a weighted Sobolev space with a weight in the corresponding Muckenhoupt class. Methods of nonlinear spline approximation are applied to derive an atomic decomposition theorem for functions in a Besov space of variable smoothness. A complete description of traces on the hyperplane of a Besov space of variable smoothness and of a weighted Besov space with a weight in the corresponding Muckenhoupt class is given. Bibliography: 27 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Comparison of Spatial Discretization Methods for Solving the S N Equations Using a Three-Dimensional Method of Manufactured Solutions Benchmark Suite with Escalating Order of Nonsmoothness

For this study, a comparison of the accuracy and computational efficiency of spatial discretization methods of the three-dimensional S N equations is conducted, including discontinuous Galerkin finite element methods, the arbitrarily high-order transport method of nodal type (AHOTN), the linear-linear method, the linear-nodal (LN) method, and the higher-order diamond difference method. For this purpose, we have developed a suite of method of manufactured solutions benchmarks that provides an exact solution of the S N equations even in the presence of scattering. Most importantly, our benchmark suite permits the user to set an arbitrary level of smoothness of the exact solution across the singular characteristics. Our study focuses on the computational efficiency of the considered spatial discretization methods. Numerical results indicate that the best-performing method depends on the norm used to measure the discretization error. We employ discrete Lp norms and integral error norms in this work. For configurations with continuous exact angular flux, high-order AHOTNs perform best under L p error norms, while the LN method performs best when measured by integral error norms. Lastly, when the angular flux is discontinuous, a new singular-characteristic tracking method for three-dimensional geometries performs best among the considered methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Comparative and Cost Analysis of a Novel Predictive Power Ramp Rate Control Method: A Case Study in a PV Power Plant in Puerto Rico

One of the most important aspects that need to be addressed to increase solar energy penetration is the power ramp-rate control. In weak grids such as the one found in Puerto Rico, it is important to smooth power fluctuations caused by the intermittence of passing clouds. In this work, a novel power ramp-rate control strategy is proposed. Additionally, a comparison with some of the most common power ramp-rate control methods is performed using a proposed model and real solar radiation data from the Coto Laurel photovoltaic power plant located in Ponce, Puerto Rico. The proposed model was validated using one-year real data from Coto Laurel. The power ramp-rate control methods were compared in real-time simulations using the OP5700 from Opal-RT Technologies considering power ramp rate fluctuations, power ramp-rate violations, fluctuations in the state-of-charge, among other indicators. Moreover, the proposed power ramp-rate control strategy, called predictive dynamic smoothing was explained and compared. Results indicate that the predictive dynamic smoothing produced a considerably reduced Levelized Cost of Storage compared to other power ramp-rate control methods and provided a higher lifetime expectancy for lithium batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Optimal smoothing length scale for actuator line models of wind turbine blades based on Gaussian body force distribution: Wind energy, actuator line model

The actuator line model (ALM) is a commonly used method to represent lifting surfaces such as wind turbine blades within large-eddy simulations (LES). In the ALM, the lift and drag forces are replaced by an imposed body force that is typically smoothed over several grid points using a Gaussian kernel with some prescribed smoothing width e. To date, the choice of e has most often been based on numerical considerations related to the grid spacing used in LES. However, especially for finely resolved LES with grid spacings on the order of or smaller than the chord length of the blade, the best choice of e is not known. In this work, a theoretical approach is followed to determine the most suitable value of e, based on an analytical solution to the linearized inviscid flow response to a Gaussian force. We find that the optimal smoothing width eopt is on the order of 14%-25% of the chord length of the blade, and the center of force is located at about 13%-26% downstream of the leading edge of the blade for the cases considered. These optimal values do not depend on angle of attack and depend only weakly on the type of lifting surface. It is then shown that an even more realistic velocity field can be induced by a 2-D elliptical Gaussian lift-force kernel. Some results are also provided regarding drag force representation.

17 WIND ENERGY↗
DOE OSTI2022

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗
NASA NTRS1989

Analysis of a second-order-accurate finite-volume method for temporally-growing compressible shear layers

A finite-volume method for solving the compressible laminar Navier-Stokes equation in two dimensions is assessed for use in the simulation of transitional flows. The method is second-order-accurate, uses alternating-direction-implicit time integration, and a total-variation-diminishing smoothing operator in the inviscid flux terms. The test problem was a free shear layer with a forced periodicity in x at a specified wavelength and with the parallel mean flow specified as a hyperbolic-tangent profile. A grid refinement investigation verified this method to be second-order-accurate with less than 4 percent error in the growth rate of isolated modes on a 32 by 64 grid. The interaction between modes of similar amplitude was small, less than 1 percent. However, in cases where there was a dominant mode, numerical phase error pumped energy into all other modes. This error results in nonphysical growth rate for modes whose energy content is several orders of magnitude below the dominant mode. In the case when the dominant mode saturates, the pumping action stops and the growth rate of the next largest growing mode regains physical significance.

Atkins, H. L.↗