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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Error analysis of numerical methods for thick diffusive neutron transport problems on Shishkin mesh

A thin layer will develop at the boundary if the incoming angular flux is anisotropic in thick diffusive neutron transport problems. Solving such singularly perturbed problems, which have non-smooth solutions with singularity near the boundary, is computationally challenging. Standard finite difference schemes on a uniform mesh cannot yield ε-uniform convergence, where ε is a small parameter, while it can be achieved on a suitable piecewise-uniform Shishkin mesh. We present a formal error analysis of the diamond difference (DD) method and step difference (SD) method for solving the S{sub N} neutron transport equation. The analysis can be extended to other finite difference methods. Numerical results are presented to confirm the error estimates and the advantages of the Shishkin mesh. (author)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1987

An adaptive grid, unsteady model for two-dimensional magnetohydrodynamic (MHD) flow

An adaptive grid finite difference method for solving multi-dimensional, time-dependent, magnetohydrodynamic (MHD) equations is developed. The method is capable of solving problems that include high gradients due to geometry, propagation of shock waves and unsteady boundary conditions. The grid generation technique is based on variational principles with direct control over grid concentration, smoothness and skewness. An example for a two-dimensional MHD simulation of the propagation of a solar-flare-generated shock wave in solar wind flow in the heliographic equatorial plane is selected for illustration of this method.

Panitchob, Supat↗
DOE OSTI2019

A novel invariant mass method to isolate resonance backgrounds from the chiral magnetic effect

The Chiral Magnetic Effect (CME) refers to charge separation along a strong magnetic field, due to topological charge fluctuations in QCD. Charge correlation (Δγ) signals consistent with CME have been first observed almost a decade ago. It has also been known since then that the Δγ is contaminated by a major background from resonance decays coupled with elliptic flow. In this contribution, we propose differential Δγ measurements as function of the pair invariant mass (m inv ). The Δγ in the high minv region is essentially free of resonance backgrounds. In the low m inv region, the Δγ backgrounds show resonance peaks. The CME signal, presumably smooth in m inv , may thus be extracted from a two-component model fit. We demonstrate the feasibility and effectiveness of this novel method by using the AMPT and toy-model Monte-Carlo simulations. We also discuss an application of the method in data analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2010

Vibration Response Predictions for Heavy Panel Mounted Components from Panel Acreage Environment Specifications

The development of new launch vehicles in the Aerospace industry often relies on response measurements taken from previously developed vehicles during various stages of liftoff and ascent, and from wind tunnel models. These measurements include sound pressure levels, dynamic pressures in turbulent boundary layers and accelerations. Rigorous statistical scaling methods are applied to the data to derive new environments and estimate the performance of new skin panel structures. Scaling methods have proven to be reliable, particularly for designs similar to the vehicles used as the basis for scaling, and especially in regions of smooth acreage without exterior protuberances or heavy components mounted to the panel. To account for response attenuation of a panel-mounted component due to its apparent mass at higher frequencies, the vibroacoustics engineer often reduces the acreage vibration according to a weight ratio first suggested by Barrett. The accuracy of the reduction is reduced with increased weight of the panel-mounted component, and does not account for low-frequency amplification of the component/panel response as a system. A method is proposed that combines acreage vibration from scaling methods with finite element analysis to account for the frequency-dependent dynamics of heavy panel-mounted components. Since the acreage and mass-loaded skins respond to the same dynamic input pressure, such pressure may be eliminated in favor of a frequency-dependent scaling function applied to the acreage vibration to predict the mass-loaded panel response. The scaling function replaces the Barrett weight ratio, and contains all of the dynamic character of the loaded and unloaded skin panels. The solution simplifies for spatially uncorrelated and fully correlated input pressures. Since the prediction uses finite element models of the loaded and unloaded skins, a rich suite of response data are available to the design engineer, including interface forces, stress and strain, as well as acceleration and displacement. An extension of the method is also developed to incorporate the effect of a local protuberance near a heavy component. Acreage environments from traditional scaling methods with and without protuberance effects serve as the basis for the extension. Authors:

Harrison, Phillip↗
NASA NTRS1996

Computational considerations for the simulation of shock-induced sound

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order-accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order-accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks in a quasi-one-dimensional nozzle flow. The design order of accuracy is achieved in the smooth regions of a steady-state test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order-accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustics simulations of flows with shocks that are raised by these results are discussed.

Casper, Jay↗
DOE OSTI2018

Applying a difference ratio method in spin-polarized scanning tunneling microscopy to determine crystalline anisotropies and antiferromagnetic spin alignment in Cr(0 0 1) c(2 $\times$ 2)

Here, we present spin-polarized scanning tunneling microscopy results for 8 micrometer thick, chromium (0 0 1)-oriented thin films deposited using molecular beam epitaxy on MgO(0 0 1) substrates, achieving a smooth, atomically-stepped staircase morphology. The c(2 × 2) surface structure of the Cr(0 0 1) surface was found via atomically-resolved scanning tunneling microscopy images and dI/dV spectroscopy. A practical method of interpreting spin-polarized scanning tunneling microscopy dI/dV images was developed and applied to the c(2 × 2) Cr(0 0 1) surface structure. A 180° spin reversal between atomic layers, characteristic of a topological, layer-wise antiferromagnetic spin structure, was quantitatively verified. In addition, the spin anisotropy directions for 3 different antiferromagnetic, atomic staircase domains occurring on the c(2 × 2) surface were quantitatively determined to coincide with the [1 0 0], [0 1 0], and [1 1 0] crystallographic axes of chromium

36 MATERIALS SCIENCE↗
DOE OSTI2022

FSISPH: An SPH formulation for impacts between dissimilar materials

Here, we present an SPH formulation with several new features designed to better model the fully-compressible interaction of dissimilar materials. We developed the new method to simulate the atmospheric entry and break-up of small celestial bodies in planetary atmospheres. The formulation uses a unity-based, density-energy discretization of the hydrodynamic conservation laws with linear-corrected kernel gradients. To account for variations in compressibility, we use an HLLC approximate Riemann solver to adjust the velocity gradient at material interfaces. To handle large transverse velocity discontinuities, we introduce a simple slip interface model that limits the artificial viscosity at material interfaces. Diffusion is optionally applied through the velocity gradient and this allows the density and specific thermal energy to evolve in a manner more consistent with the first law of thermodynamics in comparison to other more direct diffusion schemes. We also introduce a material-local second-order artificial conduction scheme used to smooth the specific thermal energy field. Material damage fits neatly under this framework by treating the damage front as a material interface. The method has been implemented as a solver, FSISPH, within the code, Spheral++, and is publicly available on github. We test our new solver on a number of classic shock, mixing, and multi-material problem. The components we outline can significantly improve accuracy of SPH for problems with sharp contact discontinuities.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Using High Resolution Design Spaces for Aerodynamic Shape Optimization Under Uncertainty

This paper explains why high resolution design spaces encourage traditional airfoil optimization algorithms to generate noisy shape modifications, which lead to inaccurate linear predictions of aerodynamic coefficients and potential failure of descent methods. By using auxiliary drag constraints for a simultaneous drag reduction at all design points and the least shape distortion to achieve the targeted drag reduction, an improved algorithm generates relatively smooth optimal airfoils with no severe off-design performance degradation over a range of flight conditions, in high resolution design spaces parameterized by cubic B-spline functions. Simulation results using FUN2D in Euler flows are included to show the capability of the robust aerodynamic shape optimization method over a range of flight conditions.

Li, Wu↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
NASA NTRS2003

Robust Airfoil Optimization in High Resolution Design Space

The robust airfoil shape optimization is a direct method for drag reduction over a given range of operating conditions and has three advantages: (1) it prevents severe degradation in the off-design performance by using a smart descent direction in each optimization iteration, (2) it uses a large number of B-spline control points as design variables yet the resulting airfoil shape is fairly smooth, and (3) it allows the user to make a trade-off between the level of optimization and the amount of computing time consumed. The robust optimization method is demonstrated by solving a lift-constrained drag minimization problem for a two-dimensional airfoil in viscous flow with a large number of geometric design variables. Our experience with robust optimization indicates that our strategy produces reasonable airfoil shapes that are similar to the original airfoils, but these new shapes provide drag reduction over the specified range of Mach numbers. We have tested this strategy on a number of advanced airfoil models produced by knowledgeable aerodynamic design team members and found that our strategy produces airfoils better or equal to any designs produced by traditional design methods.

Li, Wu↗
DOE OSTI2016

MODELING THE Ly α FOREST IN COLLISIONLESS SIMULATIONS

Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining the high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present “Iteratively Matched Statistics” (IMS), a novel method to accurately model the Ly α forest with collisionless N -body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Ly α forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N -body simulation, which we construct from the matter density. We demonstrate that our method can reproduce the PDF, line of sight and 3D power spectra of the Ly α forest with good accuracy (7%, 4%, and 7% respectively). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20%–80%), especially for N -body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic “mock” skies for Ly α forest surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

VALIDATION AND DETERMINATION OF SIGNIFICANT SIMULATION PARAMETERS USING THE SMOOTHED PARTICLE HYDRODYNAMIC CODE NEUTRINO

In recent years there has been an increase in flooding concern to nuclear power plants; thus, tools to reduce the cost of mitigation efforts are valuable to industry. Simulating flooding events is achieved through different methods and at varied levels of detail. However, these simulations must be able to match application conditions, which requires code validation as well as knowledge of which parameters affect the accuracy of the simulation. To validate and determine the significant parameters for these tools, a sloshing experiment was conducted at the George Washington University and then simulated with the smoothed particle hydrodynamic code Neutrino. The simulation setup was constructed to match the experiment as closely as possible, and the end-wall pressure results were analysed. The 90% and 50% end-wall pressure bounds were compared along with the pressure impulse. The simulation bounds fell mostly within the experimental bounds, with exceptions in insignificant low-pressure areas, such as at the beginning (~40% difference) and the end of the cycle (~40% difference). Another exception was seen at a point in the peak pressure (~75% difference). The simulated pressure impulse matched the experiments well and was within 10%. Additionally, seven parameters were qualitatively assessed and then investigated to determine their significance. The parameters’ significance was determined by sampling five values across the parameters’ range using the Risk Analysis Virtual Environment, then comparing the average pressure for each run. All seven parameters were identified as having some significance. However, the particle size, interaction-radius to particle-size ratio, and fluid settling seem to have a greater effect due to the larger fluctuation of the results. Now that the significant parameters are identified, further research is being done to quantify the significance of each parameter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Laccase-mediated synthesis of lignin-core hyperbranched copolymers

Lignin, one of the major chemical constituents of woody biomass, is the second most abundant biopolymer found in nature. The pulp and paper industry has long produced lignin on the scale of millions of tons annually as a by-product of the pulping process, and the dawn of cellulosic ethanol production has further contributed to this amount. Historically, lignin has been perceived as a waste material and burned as a fuel for the pulping process. But, recent research has been geared toward developing cost-effective technologies to convert lignin into valuable commodities. Attributing to the polyphenolic structure of lignin, enzymatic modification of its surface using laccases (benzenediol:oxygen oxidoreductases, EC 1.10.3.2) has demonstrated to be highly successful. The current study aims at developing lignin-core hyperbranched copolymers via the laccase-assisted copolymerization of kraft lignin with methylhydroquinone and a trithiol. Based on the physical properties of the resulting material, it is likely that crosslinking reactions have taken place during the drying process to produce a copolymeric network rather than discrete hyperbranched copolymers, with NMR data providing evidence of covalent bonding between monomers. A preliminary thermal analysis data reveals that the copolymeric material possesses a moderate glass transition temperature and exhibits good thermostability, thus may have potential application as a lignin-based thermoplastic. Scanning electron microscopy images confirm the smooth, glossy surface of the material and that it is densely packed. Our results are a sustainable, ecofriendly, economic method to create an exciting novel biomaterial from a renewable feedstock while further enhancing lignin valorization.

09 BIOMASS FUELS↗
DOE OSTI2021

Six-dimensional matching of intense beam with linear accelerating structure

Beam matching is a common technique that is routinely employed in accelerator design with the aim of minimizing beam losses and preservation of beam brightness. Despite being widely used, a full theoretical understanding of beam matching in 6D remains elusive. In this work, we present an analytical treatment of 6D beam matching of a high-intensity beam onto an RF structure. We begin our analysis within the framework of a linear model, and apply the averaging method to a set of 3D beam envelope equations. Accordingly, we obtain a matched solution that is comprised of smoothed envelopes and periodic terms, describing envelope oscillations with the period of the focusing structure. We then consider the nonlinear regime, where the beam size is comparable with the separatrix size. Stating with a Hamiltonian analysis in 6D phase space, we attain a self-consistent beam profile and show that it is significantly different from the commonly used ellipsoidal shape. Subsequently, we analyze the special case of an equilibrium with equal space charge depression between all degrees of freedom. Comparison of beam dynamics for equipartitioned, equal space charge depression, and equal emittances beams is given. Finally, we present experimental results on beam matching in the LANSCE linac.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Impact of Short-Range Forces on Defect Production from High-Energy Collisions

Primary radiation damage formation in solid materials typically involves collisions between atoms that have up to a few hundred keV of kinetic energy. The distance between two colliding atoms can approach 0.05 nm during these collisions. At such small atomic separations, force fields fitted to equilibrium properties tend to significantly underestimate the potential energy of the colliding dimer. To enable molecular dynamics simulations of high-energy collisions, it is common practice to use a screened Coulomb force field to describe the interactions and to smoothly join this to the equilibrium force field at a suitable interatomic spacing. But, there is no accepted standard method for choosing the parameters used in the joining process, and our results prove that defect production is sensitive to how the force fields are linked. A new procedure is presented that involves the use of ab initio calculations to determine the magnitude and spatial dependence of the pair interactions at intermediate distances, along with systematic criteria for choosing the joining parameters. Results are presented for the case of nickel, which demonstrate the use and validity of the procedure.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

On the determination of the crystal-vapor surface free energy, and why a Gaussian expression can be accurate for a system far from Gaussian

The crystal-vapor surface free energy γ is an important physical parameter governing physical processes, such as wetting and adhesion. We explore exact and approximate routes to calculate γ based on cleaving an intact crystal into non-interacting sub-systems with crystal-vapor interfaces. We do this by turning off the interactions, ΔV, between the sub-systems. Using the soft-core scheme for turning off ΔV, we find that the free energy varies smoothly with the coupling parameter λ, and a single thermodynamic integration yields the exact γ. We generate another exact method, and a cumulant expansion for γ by expressing the surface free energy in terms of an average of e{sup −βΔV} in the intact crystal. The second cumulant, or Gaussian approximation for γ is surprisingly accurate in most situations, even though we find that the underlying probability distribution for ΔV is clearly not Gaussian. We account for this fact by developing a non-Gaussian theory for γ and find that the difference between the non-Gaussian and Gaussian expressions for γ consist of terms that are negligible in many situations. Exact and approximate methods are applied to the (111) surface of a Lennard-Jones crystal and are also tested for more complex molecular solids, the surface of octane and nonadecane. Alkane surfaces were chosen for study because their crystal-vapor surface free energy has been of particular interest for understanding surface freezing in these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Enhanced emissive and lasing characteristics of nano-crystalline MAPbBr 3 films grown via anti-solvent precipitation

Here, control of the nano-scale crystal size distribution in films of organic-inorganic lead-tri-bromide perovskites is achieved through a recently reported solution-based, anti-solvent treatment method [H. Cho et al., Science 350, 1222 (2015)]. The treated films are shown to be high quality, optically smooth with excellent emissive and optical gain properties including higher photoluminescence yield and reduced threshold for laser action. The improved lasing is shown to originate from a synergistic combination of a reduction in waveguide losses and a decrease in the non-radiative decay rate of the excited state population, compared to untreated films. The net gain is increased by a factor of two upon anti-solvent treatment and reaches a respectable value of ≈300 cm ~1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Detecting hidden transient events in noisy nonlinear time-series

The information impulse function (IIF), running Variance, and local Hölder Exponent are three conceptually different time-series evaluation techniques. These techniques examine time-series for local changes in information content, statistical variation, and point-wise smoothness, respectively. Using simulated data emulating a randomly excited nonlinear dynamical system, this study interrogates the utility of each method to correctly differentiate a transient event from the background while simultaneously locating it in time. Computational experiments are designed and conducted to evaluate the efficacy of each technique by varying pulse size, time location, and noise level in time-series. Our findings reveal that, in most cases, the first instance of a transient event is more easily observed with the information-based approach of IIF than with the Variance and local Hölder Exponent methods. While our study highlights the unique strengths of each technique, the results suggest that very robust and reliable event detection for nonlinear systems producing noisy time-series data can be obtained by incorporating the IIF into the analysis.

42 ENGINEERING↗