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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Modeling Alzheimer’s Disease Progression with Fused Laplacian Sparse Group Lasso

Alzheimer’s disease (AD), the most common type of dementia, not only imposes a huge financial burden on the health care system, but also a psychological and emotional burden on patients and their families. There is furthermore an urgent need to infer trajectories of cognitive performance over time and identify biomarkers predictive of the progression. Here, we introduce the multi-task learning with fused Laplacian sparse group lasso model, which can identify biomarkers closely related to cognitive measures due to its sparsity-inducing property, and model the disease progression with a general weighted (undirected) dependency graphs among the tasks. An efficient alternative directions method of multipliers based optimization algorithm is derived to solve the proposed non-smooth objective formulation. The effectiveness of the presented model is demonstrated by its superior prediction performance over multiple state-of-the-art methods and accurate identification of compact sets of cognition-relevant imaging biomarkers that are consistent with prior medical studies.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1982

Transonic flow calculations using Euler's equations

An implicit finite difference method with implicit boundary conditions is employed to solve the steady Euler's equations for flows past arbitrary geometries. The resulting code is used to investigate in a systematic way various aspects of flow past airfoils at transonic speeds such as, method of solution, boundary conditions, grid stretching and generation, shock and sonic point operators, the Kutta condition, and smoothing. Results obtained are in good agreement with results of other codes. Moreover, it appears that the method of solution employed is such that Kutta's condition need not be invoked. This statement appears to be valid for other existing schemes employed in the solution of Euler's equations.

Atkins, H. L.↗
DOE OSTI2024

A comprehensive review of dwell time optimization methods in computer-controlled optical surfacing

Dwell time plays a vital role in determining the accuracy and convergence of the computer-controlled optical surfacing process. However, optimizing dwell time presents a challenge due to its ill-posed nature, resulting in non-unique solutions. To address this issue, several well-known methods have emerged, including the iterative, Bayesian, Fourier transform, and matrix-form methods. Despite their independent development, these methods share common objectives, such as minimizing residual errors, ensuring dwell time's positivity and smoothness, minimizing total processing time, and enabling flexible dwell positions. This paper aims to comprehensively review the existing dwell time optimization methods, explore their interrelationships, provide insights for their effective implementations, evaluate their performances, and ultimately propose a unified dwell time optimization methodology.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Statistical properties of filaments in the cosmic web

ABSTRACT In the context of the cosmological and constrained Exploring the Local Universe with the reConstructed Initial Density field (ELUCID) simulation, this study explores the statistical characteristics of filaments within the cosmic web, focussing on aspects such as the distribution of filament lengths and their radial density profiles. Using the classification of the cosmic web environment through the Hessian matrix of the density field, our primary focus is on how cosmic structures react to the two variables $R_{\rm s}$ and $\lambda _{\rm th}$. The findings show that the volume fractions of knots, filaments, sheets, and voids are highly influenced by the threshold parameter $\lambda _{\rm th}$, with only a slight influence from the smoothing length $R_{\rm s}$. The central axis of the cylindrical filament is pinpointed using the medial-axis thinning algorithm of the COsmic Web Skeleton (COWS) method. It is observed that median filament lengths tend to increase as the smoothing lengths increase. Analysis of filament length functions at different values of $R_{\rm s}$ indicates a reduction in shorter filaments and an increase in longer filaments as $R_{\rm s}$ increases, peaking around $2.5R_{\rm s}$. The study also shows that the radial density profiles of filaments are markedly affected by the parameters $R_{\rm s}$ and $\lambda _{\rm th}$, showing a valley at approximately $2R_{\rm s}$, with increases in the threshold leading to higher amplitudes of the density profile. Moreover, shorter filaments tend to have denser profiles than their longer counterparts.

Zhang, Youcai (ORCID:0000000319674091)↗
NASA NTRS1993

On the prediction of multigrid efficiency through local mode analysis

A single grid local mode analysis is used to predict the smoothing properties of numerical schemes for solving the Navier-Stokes equations with factorization based on Stone's Strongly Implicit Method. Four difference approximations for the convection terms are considered, namely, hybrid, central, second-order upwind, and third-order upwind. Smoothing factors from the analysis are compared with practical convergence factors in a multigrid method for flow over a backward facing step and it is found that the local mode analysis correctly predicts the effects of Reynolds number and higher-order schemes.

Wilson, R. V.↗
DOE OSTI2024

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-F-T-343: Dosimetric Investigation for Fluence Smoothing On VMAT Plans in Monaco for Head - Neck Cancers

Purpose: The study investigates the effect of fluence smoothing parameter on VMAT plans for ten head-neck cancer patients using Monaco5.00.04. Methods: VMAT plans were created using Monaco5.00.04 planning system for 10 head-neck patients. Four plans were generated for each patient using available smoothing parameters i.e. high, medium, low and off. The number of monitor units required to deliver 1 cGy was defined as a modulation degree; and was taken as a measure of plan complexity. Routinely used plan quality parameters Conformity index (CI) and Homogeneity index (HI) were used in the study. As a protocol our center, practices “medium” smoothing for clinical implementation. Plans with medium smoothing were opted as reference plans due to the clinical acceptance and dosimetric verifications made on these plans. Plans were generated by varying the smoothing parameter and re-optimization was done. The PTV was evaluated for D98%, D95%, D50%, D1% and prescription isodose volume (PIV). For critical organs; spine and parotids the parameters recorded were D1cc and Dmean respectively. Results: The cohort had the median prescription as 6000 cGy in the range of 6600 cGy - 4500 cGy. The modulation degree was observed to increase up to 6% from reference to the most complex plan. High smoothing had about 11% increase in segments which marginally (0.5 to 1%) increased the homogeneity index while conformity index remains constant. For spine the maximum D1cc was observed in medium smoothing as 4639.8 cGy, this plan was clinically accepted and dosimetrically verified. Similarly for parotids, the Dmean was 2011.9 cGy and 1817.05 cGy. Conclusion: The sensitivity of plan quality in terms of smoothing options (high, medium, low and off) available in Monaco 5.00.04 was resulted in minimal difference in terms of target coverage, conformity index and homogeneity index. Similarly changing smoothing did not result in any enhanced advantage in sparing of critical organs.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗
NASA NTRS2002

High-Order Spectral Volume Method for 2D Euler Equations

The Spectral Volume (SV) method is extended to the 2D Euler equations. The focus of this paper is to study the performance of the SV method on multidimensional non-linear systems. Implementation details including total variation diminishing (TVD) and total variation bounded (TVB) limiters are presented. Solutions with both smooth features and discontinuities are utilized to demonstrate the overall capability of the SV method.

Wang, Z. J.↗
DOE OSTI2020

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Soft chemical synthesis of silicon nanosheets and their applications

Two-dimensional silicon nanomaterials are expected to show different properties from those of bulk silicon materials by virtue of surface functionalization and quantum size effects. Since facile fabrication processes of large area silicon nanosheets (SiNSs) are required for practical applications, a development of soft chemical synthesis route without using conventional vacuum processes is a challenging issue. We have recently succeeded to prepare SiNSs with sub-nanometer thicknesses by exfoliating layered silicon compounds, and they are found to be composed of crystalline single-atom-thick silicon layers. In this review, we present the synthesis and modification methods of SiNSs. These SiNSs have atomically flat and smooth surfaces due to dense coverage of organic moieties, and they are easily self-assembled in a concentrated state to form a regularly stacked structure. We have also characterized the electron transport properties and the electronic structures of SiNSs. Finally, the potential applications of these SiNSs and organic modified SiNSs are also reviewed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A tesselation-based model for intensity estimation and laser plasma interactions calculations in three dimensions

Here, a three-dimensional laser propagation model for computation of laser-plasma interactions is presented. It is focused on indirect drive geometries in inertial confinement fusion and formulated for use at large temporal and spatial scales. A modified tesselation-based estimator and a relaxation scheme are used to estimate the intensity distribution in plasma from geometrical optics rays. Comparisons with reference solutions show that this approach is well-suited to reproduce realistic 3D intensity field distributions of beams smoothed by phase plates. It is shown that the method requires a reduced number of rays compared to traditional rigid-scale intensity estimation. Using this field estimator, we have implemented laser refraction, inverse-bremsstrahlung absorption, and steady-state crossed-beam energy transfer with a linear kinetic model in the numerical code VAMPIRE. Probe beam amplification and laser spot shapes are compared with experimental results and pf3d paraxial simulations. These results are promising for the efficient and accurate computation of laser intensity distributions in holhraums, which is of importance for determining the capsule implosion shape and risks of laser-plasma instabilities such as hot electron generation and backscatter in multi-beam configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Parametric interaction of optical waves in metamaterials under low-frequency pumping

The influence of phase effects under three-wave parametric interaction and low-frequency pumping in metamaterials is studied in the case of a negative refractive index at a signal-wave frequency. It is found that the efficiency of the backward signal-wave amplification is the higher, the greater the ratio of the intensities of the idler and signal waves at the input to the metamaterial. An increase in the idler wave intensity at the input by five times, as compared to the signal-wave intensity, leads to a nonlinear increase in the signal-wave amplification by almost 20 times. According to the analytic expressions obtained in the constant-intensity approximation, the choice of the optimal parameters for the pump intensity, total length of the metamaterial and phase detuning will facilitate the implementation of regimes of effective amplification and generation of the signal wave. A comparison is made with the results obtained in the constant-field approximation, and a numerical estimate of the expected efficiency of the frequency conversion is presented. Control of frequency and pump power is shown to make possible the smooth tuning of the parametric converter frequency. The developed method can be used to design frequency converters based on nonlinear metamaterials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Attribute-Aware RBFs: Interactive Visualization of Time Series Particle Volumes Using RT Core Range Queries

Smoothed-particle hydrodynamics (SPH) is a mesh-free method used to simulate volumetric media in fluids, astrophysics, and solid mechanics. Visualizing these simulations is problematic because these datasets often contain millions, if not billions of particles carrying physical attributes and moving over time. Radial basis functions (RBFs) are used to model particles, and overlapping particles are interpolated to reconstruct a high-quality volumetric field; however, this interpolation process is expensive and makes interactive visualization difficult. Existing RBF interpolation schemes do not account for color-mapped attributes and are instead constrained to visualizing just the density field. To address these challenges, we exploit ray tracing cores in modern GPU architectures to accelerate scalar field reconstruction. We use a novel RBF interpolation scheme to integrate per-particle colors and densities, and leverage GPU-parallel tree construction and refitting to quickly update the tree as the simulation animates over time or when the user manipulates particle radii. We also propose a Hilbert reordering scheme to cluster particles together at the leaves of the tree to reduce tree memory consumption. Finally, we reduce the noise of volumetric shadows by adopting a spatially temporal blue noise sampling scheme. Our method can provide a more detailed and interactive view of these large, volumetric, time-series particle datasets than traditional methods, leading to new insights into these physics simulations.

Particle Volumes↗