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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling Turbulence in Galactic Centers

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way’s Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Gravity Field of Venus and Comparison with Earth

The acceleration (gravity) anomaly estimates by spacecraft tracking, determined from Doppler residuals, are components of the gravity field directed along the spacecraft Earth line of sight (LOS). These data constitute a set of vector components of a planet's gravity field, the specific component depending upon where the Earth happened to be at the time of each measurement, and they are at varying altitudes above the planet surface. From this data set the gravity field vector components were derived using the method of harmonic splines which imposes a smoothness criterion to select a gravity model compatible with the LOS data. Given the piecewise model it is now possible to upward and downward continue the field quantities desired with a few parameters unlike some other methods which must return to the full dataset for each desired calculation.

Bowin, C.↗
DOE OSTI2018

A Comparison of Different Methods of Separation of Continuous Overlapping Spectral Lines

Several methods of separation of continuous overlapping spectral lines are compared. In the majority of methods, the profile of each line is modeled by a Gaussian or Lorentzian, and total measured spectrum z is processed. Number of lines N and their parameters are usually estimated by the method of derivatives; however, the differentiation of noisy spectrum z causes large errors. To improve the differentiation accuracy, it is proposed to use smoothing splines. In the Fourier-self-deconvolution method, apodization (artificial truncation of the interferogram) is used to resolve overlapping lines, which makes it possible to resolve lines, but at the expense of a significant decrease in their widths. In this work, reducing the widths of lines is not used for their resolution, but rather, true line profiles are reconstructed by minimizing the residual functional with the modified coordinate-descent method using the decremental-constraint technique, and also, for comparison, with the Nelder–Mead simplex method. In the Manoilov method, the parameters of lines (peaks) are determined from convolutions of spectrum derivatives with individual peaks. In that method, the notion of the degree of overlap has been also introduced. In this work, we introduce a generalized degree of overlap for the case in which the amplitudes, widths, and spacings between neighboring lines are different. Numerical illustrations are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Improved ALE mesh velocities for complex flows

A key choice in the development of arbitrary Lagrangian-Eulerian solution algorithms is how to move the computational mesh. The most common approaches are smoothing and relaxation techniques, or to compute a mesh velocity field that produces smooth mesh displacements. We present a method in which the mesh velocity is specified by the irrotational component of the fluid velocity as computed from a Helmholtz decomposition, and excess compression of mesh cells is treated through a noniterative, local spring-force model. This approach allows distinct and separate control over rotational and translational modes. In conclusion, the utility of the new mesh motion algorithm is demonstrated on a number of 3D test problems, including problems that involve both shocks and significant amounts of vorticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Entropy Stable Conservative Flux Form Neural Networks

We propose an entropy-stable conservative flux form neural network (CFN) to predict the dynamics of unknown governing conservation laws. The design of the network is based on the entropy-stable, second-order, and non-oscillatory Kurganov-Tadmor (KT) scheme. The proposed entropy-stable CFN, hereafter referred to as ESCFN, uses slope limiting as a denoising mechanism, ensuring accurate predictions in both noisy and sparse observation environments, as well as in both smooth and discontinuous regions. Importantly, our method is designed to predict long term dynamics of the unknown conservation law exclusively from a short temporal window of observed data, that is, without oracle knowledge of the PDE or later-time solution profiles. Numerical experiments demonstrate that the ESCFN achieves both stability and conservation while maintaining accuracy over extended time domains, and successfully predicts shock propagation speeds in long-term simulations. Furthermore, it is also robust to both noisy and sparse data environments.

Hyperbolic conservation laws↗
DOE OSTI2019

A higher order approximate static condensation method for multi-material diffusion problems

The paper studies an approximate static condensation method for the diffusion problem with discontinuous diffusion coefficients. The method allows for a general polygonal mesh which is unfitted to the material interfaces. Moreover, the interfaces can be discontinuous across the mesh edges as typical for numerical reconstructions using the volume or moment-of-fluid methods. We apply a mimetic finite difference method to solve local diffusion problems and use P 1 (mortar) edge elements to couple local problems into the global system. The condensation process and the properties of the resulting algebraic system are discussed. It is demonstrated that the method is second order accurate on smooth solutions and performs well for problems with high contrast in diffusion coefficients. Experiments also show the robustness with respect to position of the material interfaces against the underlying mesh.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Evaluation of reconstruction errors and identification of artefacts for JET gamma and neutron tomography

The Joint European Torus (JET) neutron profile monitor ensures 2D coverage of the gamma and neutron emissive region that enables tomographic reconstruction. Due to the availability of only two projection angles and to the coarse sampling, tomographic inversion is a limited data set problem. Several techniques have been developed for tomographic reconstruction of the 2-D gamma and neutron emissivity on JET, but the problem of evaluating the errors associated with the reconstructed emissivity profile is still open. The reconstruction technique based on the maximum likelihood principle, that proved already to be a powerful tool for JET tomography, has been used to develop a method for the numerical evaluation of the statistical properties of the uncertainties in gamma and neutron emissivity reconstructions. The image covariance calculation takes into account the additional techniques introduced in the reconstruction process for tackling with the limited data set (projection resampling, smoothness regularization depending on magnetic field). The method has been validated by numerically simulations and applied to JET data. Different sources of artefacts that may significantly influence the quality of reconstructions and the accuracy of variance calculation have been identified.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

Synchrotron x-ray study of a low roughness and high efficiency K 2 CsSb photocathode during film growth

Reduction of roughness to the nm level is critical of achieving the ultimate performance from photocathodes used in high gradient fields. The thrust of this paper is to explore the evolution of roughness during sequential growth, and to show that deposition of multilayer structures consisting of very thin reacted layers results in an nm level smooth photocathode. Synchrotron x-ray methods were applied to study the multi-step growth process of a high efficiency K 2 CsSb photocathode. We observed a transition point of the Sb film grown on Si at the film thickness of similar to 40 angstrom with the substrate temperature at 100 degrees C and the growth rate at 0.1 Å s -1 . The final K 2 CsSb photocathode exhibits a thickness of around five times that of the total deposited Sb film regardless of how the Sb film was grown. The film surface roughening process occurs first at the step when K diffuses into the crystalline Sb. Furthermore, the photocathode we obtained from the multi-step growth exhibits roughness in an order of magnitude lower than the normal sequential process. X-ray diffraction measurements show that the material goes through two structural changes of the crystalline phase during formation, from crystalline Sb to K 3 Sb and finally to K 2 CsSb.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Revealing Fine Structure in the Antineutrino Spectra from a Nuclear Reactor

Here, we calculate the Inverse Beta Decay (IBD) antineutrino spectrum generated by nuclear reactors using the summation method to understand deviations from the smooth Huber-Mueller model due to the decay of individual fission products, showing that plotting the ratio of two adjacent spectra points can effectively reveal these deviations. We obtain that for binning energies of 0.1 MeV or lower, abrupt changes in the spectra due to the jagged nature of the individual antineutrino spectra could be observed for highly precise experiments. Surprisingly, our calculations also reveal a peak-like feature in the adjacent points ratio plot at 4.5 MeV even with a 0.25 MeV binning interval, which we find is present in the IBD spectrum published by Daya Bay in 2016. We show that this 4.5 MeV feature is caused by the contributions of just four fission products, 95 Y, 98,101 Nb and 102 T c . This would be the first evidence of the decay of a few fission products in the IBD antineutrino spectrum from a nuclear reactor. This result is supported by applying the same numerical technique to measured aggregate electron spectra.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Giant impacts stochastically change the internal pressures of terrestrial planets

Pressure is a key parameter in the physics and chemistry of planet formation and evolution. Previous studies have erroneously assumed that internal pressures monotonically increase with the mass of a body. Using smoothed particle hydrodynamics and potential field method calculations, we demonstrate that the hot, rapidly rotating bodies produced by giant impacts can have much lower internal pressures than cool, slowly rotating planets of the same mass. Pressures subsequently increase because of thermal and rotational evolution of the body. Using the Moon-forming impact as an example, we show that the internal pressures after the collision could have been less than half that in present-day Earth. The current pressure profile was not established until Earth cooled and the Moon receded, a process that may take up to tens of millions of years. Our work defines a new paradigm for pressure evolution during accretion of terrestrial planets: stochastic changes driven by impacts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Extending Dynamic Range for Untargeted Ion Mobility-MS Workflows

Preprocessing and identification algorithms can boost the quality and reproducibility of results in untargeted ion mobility (IM)-MS workflows. For this CRADA we developed two methods for pre-processing to improve the quality of Agilent metabolomics data: i) a new saturation repair preprocessing method and ii) multi-dimensional smoothing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Modeling the Radio Foreground for Detection of CMB Spectral Distortions from the Cosmic Dawn and the Epoch of Reionization

Cosmic baryon evolution during the Cosmic Dawn and Reionization results in redshifted 21-cm spectral distortions in the cosmic microwave background (CMB). These encode information about the nature and timing of first sources over redshifts 30–6 and appear at meter wavelengths as a tiny CMB distortion along with the Galactic and extragalactic radio sky, which is orders of magnitude brighter. Therefore, detection requires precise methods to model foregrounds. We present a method of foreground fitting using maximally smooth (MS) functions. We demonstrate the usefulness of MS functions over traditionally used polynomials to separate foregrounds from the Epoch of Reionization (EoR) signal. We also examine the level of spectral complexity in plausible foregrounds using GMOSS, a physically motivated model of the radio sky, and find that they are indeed smooth and can be modeled by MS functions to levels sufficient to discern the vanilla model of the EoR signal. We show that MS functions are loss resistant and robustly preserve EoR signal strength and turning points in the residuals. Finally, we demonstrate that in using a well-calibrated spectral radiometer and modeling foregrounds with MS functions, the global EoR signal can be detected with a Bayesian approach with 90% confidence in 10 minutes’ integration.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

A study of the diurnal variation in the thermosphere as derived by satellite drag.

Results of an analysis performed on 29,574 densities derived from the drag of ten satellites to determine the amplitude and shape of the diurnal-variation curve for the thermosphere. Although a considerable amount of smoothing is inherent in the drag method, it seemed useful to determine whether any change in the shape of the diurnal-variation curve with height, latitude, and solar activity could be detected. None was detected: the curve remains remarkably stable, symmetric, with a maximum at 1420 hr LST and a minimum at 0220 hr LST. A systematic variation of the temperature range with height is observed when static models are used to derive it.

Jacchia, L. G.↗
NASA NTRS1973

A study of the diurnal variation in the thermosphere as derived by satellite drag

An analysis was performed on 29,574 densities derived from the drag of 10 satellites to determine simultaneously the parameters of the solar-activity effect in the thermosphere on the one hand, and the amplitude and shape of the diurnal-variation curve on the other. This paper reports on the study of the diurnal variation only. Although a considerable amount of smoothing is inherent in the drag method, it seemed useful to see whether any change could be detected in the shape of the diurnal-variation curve with height, latitude and solar activity. None was detected; the curve remains remarkably stable, symmetric, with a maximum at 14 hr 20 min L.S.T. and a minimum at 2 hr 20 min L.S.T. A systematic variation of the temperature range with height is observed when static models are used to derive it.

Jacchia, L. G.↗
NASA NTRS1975

An experimental investigation of heat transfer to reusable surface insulation tile array gaps in a turbulent boundary layer with pressure gradient

An experimental investigation was performed to determine the effect of pressure gradient on the heat transfer to space shuttle reusable surface insulation (RSI) tile array gaps under thick, turbulent boundary layer conditions. Heat transfer and pressure measurements were obtained on a curved array of full-scale simulated RSI tiles in a tunnel wall boundary layer at a nominal freestream Mach number of 10.3 and freestream unit Reynolds numbers of 1.6, 3.3, and and 6.1 million per meter. Transverse pressure gradients were induced over the model surface by rotating the curved array with respect to the flow. Definition of the tunnel wall boundary layer flow was obtained by measurement of boundary layer pitot pressure profiles, and flat plate wall pressure and heat transfer. Flat plate wall heat transfer data were correlated and a method was derived for prediction of smooth, curved array heat transfer in the highly three-dimensional tunnel wall boundary layer flow and simulation of full-scale space shuttle vehicle pressure gradient levels was assessed.

Throckmorton, D. A.↗
NASA NTRS1977

Landsat-D thematic mapper simulation using aircraft multispectral scanner data

A simulation of imagery from the upcoming Landsat-D Thematic Mapper was accomplished by using selected channels of aircraft 24-channel multispectral scanner data. The purpose was to simulate Thematic Mapper 30-meter resolution imagery, to compare its spectral quality with the original aircraft MSS data, and to determine changes in thematic classification accuracy for the simulated imagery. The original resolution of approximately 7.5 meters IFOV and simulated resolution of 15, 30, and 60 meters were used to indicate the trend of spectral quality and classification accuracy. The study was based in a 6.5 square kilometer area of urban Los Angeles having a diversity of land use. The original imagery was reduced in resolution by two related methods: pixel matrix averaging, and matrix smoothing with a unity box filter, followed by matrix averaging. Thematic land use classification using training sites and a Bayesian maximum-likelihood algorithm was performed at three levels of standard deviation - 1.0, 2.0, and 3.0 sigma. Plots of relative standard deviation showed that for larger training sites with a normal distribution of data, as the resolution decreased, the distribution range of density values decreased. Also, the classification accuracies for three levels of standard deviation increased as resolution decreased. However, the indication is that a point of diminishing returns had been reached, and 30 meters IFOV should be the best for multispectral classification of urban scenes.

Clark, J.↗
NASA NTRS1984

Measurement of local connective heat transfer coefficients of four ice accretion shapes

In the analytical study of ice accretions that form on aerodynamic surfaces (airfoils, engine inlets, etc.) it is often necessary to be able to calculate convective heat transfer rates. In order to do this, local convective heat transfer coefficients for the ice accretion shapes must be known. In the past, coefficients obtained for circular cylinders were used as an approximation to the actual coefficients since no better information existed. The purpose of this experimental study was to provide local convective heat transfer coefficients for four shapes that represent ice accretions. The shapes were tested with smooth and rough surfaces. The experimental method chosen was the thin-skin heat rate technique. Using this method local Nusselt numbers were determined for the ice shapes. In general it was found that the convective heat transfer was higher in regions where the model's surfaces were convex and lower in regions where the model's surfaces were concave. The effect of roughness was to increase the heat transfer in the high heat transfer regions by approximately 100% while little change was apparent in the low heat transfer regions.

Smith, M. E.↗