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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rapid scalable plasma processing of thin-film Li–La–Zr–O solid-state electrolytes

Solid-state electrolytes, such as lithium lanthanum zirconium oxide (LLZO), show promise as technologies for next-generation high-energy-density batteries, but commercial development has been hindered by a lack of scalable processing methods. Current fabrication methods are costly or require long annealing steps to create dense films. We report an atmospheric pressure blown-arc nitrogen plasma jet process to rapidly form sub-micrometer-thick, dense amorphous LLZO (a-LLZO) films from sol-gel precursors. Films are processed in less than 2 min, an order of magnitude faster than what has previously been reported. We demonstrate 500-nm-thick a-LLZO films processed at 350°C with an ionic conductivity of 2 × 10 −6 S/cm at 30°C and 2 × 10 −3 S/cm at 100°C and a conductance of 19 S at 100°C, the highest conductance of any LLZO phase to date. Here, the films exhibit outstanding smooth surface morphology with low defectivity, advancing atmospheric plasma processing as a scalable processing method for solid-state electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2018

Reevaluating growing season length controls on net ecosystem production in evergreen conifer forests

Abstract Growing season length (GSL) is a key unifying concept in ecology that can be estimated from eddy covariance-derived estimates of net ecosystem production (NEP). Previous studies disagree on how increasing GSLs may affect NEP in evergreen coniferous forests, potentially due to the variety of methods used to quantify GSL from NEP. We calculated GSL and GSL-NEP regressions at eleven evergreen conifer sites across a broad climatic gradient in western North America using three common approaches: (1) variable length (3–7 days) regressions of day of year versus NEP, (2) a smoothed threshold approach, and (3) the carbon uptake period, followed by a new approach of a method-averaged ensemble. The GSL and the GSL-NEP relationship differed among methods, resulting in linear relationships with variable sign, slope, and statistical significance. For all combinations of sites and methods, the GSL explained between 6% and 82% of NEP with p -values ranging from 0.45 to < 0.01. These results demonstrate the variability among GSL methods and the importance of selecting an appropriate method to accurately project the ecosystem carbon cycling response to longer growing seasons in the future. To encourage this approach in future studies, we outline a series of best practices for GSL method selection depending on research goals and the annual NEP dynamics of the study site(s). These results contribute to understanding growing season dynamics at ecosystem and continental scales and underscore the potential for methodological variability to influence forecasts of the evergreen conifer forest response to climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Reaction analogy based forcing for incompressible scalar turbulence

Here, we present a novel reaction analogy (RA) based forcing method for generating statistically stationary scalar fields in incompressible turbulence. The new method can produce more general scalar probability density functions (PDFs), for example, quasi-double- δ PDF, than current methods, while ensuring that scalar fields remain bounded, unlike existent forcing methodologies that can potentially violate naturally existing bounds. Such features are useful for generating initial fields in nonpremixed combustion, inlet conditions for spatially developing flows, or for studying non-Gaussian scalar turbulence. The RA method mathematically models hypothetical chemical reactions that convert reactants in a mixed state back into its pure unmixed components. Various types of chemical reactions are formulated and the corresponding mathematical expressions derived such that the reaction term is smooth in scalar space and is consistent with mass conservation. For large values of the scalar forcing rate, the method produces statistically stationary quasi-double- δ scalar PDFs. Quasiuniform, Gaussian, and stretched exponential scalar statistics are recovered for smaller values of the scalar forcing rate. The shape of the scalar PDF can be further controlled by changing the stoichiometric coefficients of the reaction. Finally, the ability of the new method to produce fully developed passive scalar fields with quasi-Gaussian PDFs is also investigated, by exploring the convergence of the scalar variance spectrum to the Obukhov-Corrsin scaling and of the third-order mixed structure function to the “four-thirds” Yaglom's law.

42 ENGINEERING↗
NASA NTRS2004

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗
DOE OSTI2023

Adaptive Online Multivariate Signal Extraction With Locally Weighted Robust Polynomial Regression

High-frequency, multivariate data collected in real-time and used to control or make decisions regarding a process’ operation often contain some noise and outliers. Thus, a method to extract the signal is needed in order to reduce the number and magnitude of control-based adjustments that are implemented. Such a method must be (i) online, depending only on past and current observations; (ii) fast, producing a smooth value more quickly than the measurement frequency; (iii) robust, ignoring brief bursts of erroneously measured values; (iv) multivariate, ignoring observations that are jointly unusual; (v) adaptive, adjusting to periods of rapid fluctuation in the signal versus periods of stability; and (vi) purely data-driven, not incorporating any information about the process from which the data are collected. Most existing methods are only able to address a subset of these six features. Furthermore, we also require the method to be nonlinear, providing a local nonlinear estimate of the signal. In this work, we propose a novel, real-time signal extraction method based on a local, robust polynomial fit. We demonstrate the performance of our method compared to a state-of-the-art competitor through simulation. For illustration, the methodology is applied to data collected from a reverse osmosis water treatment process.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2013

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cross Inference of Throughput Profiles Using Micro Kernel Network Method

Dedicated network connections are being increasingly deployed in cloud, centralized and edge computing and data infrastructures, whose throughput profiles are critical indicators of the underlying data transfer performance. Due to the cost and disruptions to physical infrastructures, network emulators, such as Mininet, are often used to generate measurements needed to estimate throughput profiles, typically expressed as a function of the connection round trip time. The profiles estimated using measurements from such emulated networks are usually inaccurate for high bandwidth and high latency connections, since they do not accurately reflect the critical network transport dynamics mainly due to computing and memory constraints of the host. We present a machine learning (ML) method to estimate the throughput profiles using emulation measurements to closely match the testbed and production network profiles. In particular, we propose a micro Kernel Network (mKN) that provides baseline throughput measurements on the host running Mininet emulations, which are used to learn a regression map that converts them to the corresponding testbed measurement estimates. Once initially learned, this map is applied to measurements from subsequent network emulations on the same host. We present experimental measurements to illustrate this approach, and derive generalization equations for the proposed mKN-ML method. Using a four-site scenario emulation, we show the effectiveness of this method in providing accurate concave throughput profiles from inaccurate convex or non-smooth ones indicated by Mininet emulation.

Rao, Nageswara↗
DOE OSTI2022

Search for heavy particles in the b -tagged dijet mass distribution with additional b -tagged jets in proton-proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS experiment

A search optimized for new heavy particles decaying to two b-quarks and produced in association with additional b-quarks is reported. The sensitivity is improved by b-tagging at least one lower-p T jet in addition to the two highest-p T jets. The data used in this search correspond to an integrated luminosity of 103 fb –1 collected with a dedicated trijet trigger during the 2017 and 2018 √s = 13 TeV proton–proton collision runs with the ATLAS detector at the LHC. The search looks for resonant peaks in the b-tagged dijet invariant mass spectrum over a smoothly falling background. The background is estimated with an innovative data-driven method based on orthonormal functions. The observed b-tagged dijet invariant mass spectrum is compatible with the background-only hypothesis. Upper limits at 95% confidence level on a heavy vector-boson production cross section times branching ratio to a pair of b-quarks are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Determination of Refractive Indices of a Layered Medium under Pulsed Irradiation

Problems related to determination of refractive indices of layered strongly scattering biological media by means of an optical tomography are investigated. Methods of reconstruction of refractive indices of layers, which are based on the smoothness properties of the solution to the initial-boundary problem, are proposed within the framework of the nonstationary radiative-transfer model. The influence of the probe-pulse duration on the accuracy of recovery of the sought-for characteristics of the medium is investigated numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Neutron Electric Dipole Moments with Clover Fermions

We present preliminary results for the contributions to the neutron EDM arising from the QCD θ-term, the Weinberg three-gluon and the quark chromo-EDM operators from our ongoing lattice calculations using clover valence quarks on the MILC HISQ lattices. We use the gradient-flow technique to smooth the lattices and renormalize the gluonic operators, and use the Schwinger source method to incorporate the quark chromo-EDM interactions in the quark propagator. For the QCD θ-term and the Weinberg three-gluon operator, we report results in the gradient-flow scheme from 8 ensembles at four lattice spacings and three pion masses, including 2 physical pion mass ensembles. For the quark chromo-EDM, unrenormalized results are presented at two lattice spacings, a=0.12 and 0.09 fm, and two pion masses, M π =310 MeV and 220 MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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↗