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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Visual exploration of high-dimensional data through subspace analysis and dynamic projections

Here, we introduce a novel interactive framework for visualizing and exploring high-dimensional datasets based on subspace analysis and dynamic projections. We assume the high-dimensional dataset can be represented by a mixture of low-dimensional linear subspaces with mixed dimensions, and provide a method to reliably estimate the intrinsic dimension and linear basis of each subspace extracted from the subspace clustering. Subsequently, we use these bases to define unique 2D linear projections as viewpoints from which to visualize the data. To understand the relationships among the different projections and to discover hidden patterns, we connect these projections through dynamic projections that create smooth animated transitions between pairs of projections. We introduce the view transition graph, which provides flexible navigation among these projections to facilitate an intuitive exploration. Finally, we provide detailed comparisons with related systems, and use real-world examples to demonstrate the novelty and usability of our proposed framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

TU-F-CAMPUS-J-04: Evaluation of Metal Artifact Reduction Technique for the Radiation Therapy Planning

Purpose: To evaluate the feasibility of using a metal artifact reduction technique in depleting metal artifact and its application in improving dose calculation in External Radiation Therapy Planning. Methods: CIRS electron density phantom was scanned with and without steel drill bits placed in some plug holes. Meta artifact reduction software with Metal Deletion Technique (MDT) was used to remove metal artifacts for scanned image with metal. Hounsfield units of electron density plugs from artifact free reference image and MDT processed images were compared. To test the dose calculation improvement after the MDT processed images, clinically approved head and neck plan with manual dental artifact correction was tested. Patient images were exported and processed with MDT and plan was recalculated with new MDT image without manual correction. Dose profiles near the metal artifacts were compared. Results: The MDT used in this study effectively reduced the metal artifact caused by beam hardening and scatter. The windmill around the metal drill was greatly improved with smooth rounded view. Difference of the mean HU in each density plug between reference and MDT images were less than 10 HU in most of the plugs. Dose difference between original plan and MDT images were minimal. Conclusion: Most metal artifact reduction methods were developed for diagnostic improvement purpose. Hence Hounsfield unit accuracy was not rigorously tested before. In our test, MDT effectively eliminated metal artifacts with good HU reproduciblity. However, it can introduce new mild artifacts so the MDT images should be checked with original images.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Improving solve time of aggregation-based adaptive AMG

In this paper, we propose improving the solve time of a bootstrap algebraic multigrid (AMG) designed previously by the authors. This is achieved by incorporating the information, a set of algebraically smooth vectors, generated by the bootstrap algorithm, in a single hierarchy by using sufficiently large aggregates, and these aggregates are compositions of aggregates already built throughout the bootstrap algorithm. The modified AMG method has good convergence properties and shows significant reduction in both memory and solve time. These savings with respect to the original bootstrap AMG are illustrated on some difficult (for standard AMG) linear systems arising from discretization of scalar and vector function elliptic partial differential equations in both 2D and 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

A special issue on the RCPs

This paper provides an editors' introduction to the special issue of Climatic Change on the RCPs. Scenarios form a crucial element in climate change research. They allow researchers to explore the long-term consequences of decisions today, while taking account of the inertia in both the socio-economic and physical system. Scenarios also form an integrating element among the different research disciplines of those studying climate change, such as economists, technology experts, climate researchers, atmospheric chemists and geologists. In 2007, the IPCC requested the scientific community to develop a new set of scenarios, as the existing scenarios (published in the Special Report on Emissions Scenarios, (Nakicenovic and Swart 2000), and called the 'SRES scenarios') needed to be updated and expanded in scope (see Moss et al. (2010) for a detailed discussion). Researchers from different disciplines worked together to develop a process to craft these new scenarios, as summarized by Moss, et al. (2010). The Integrated Assessment Modeling Consortium (IAMC), founded in response to the IPCC call, played a key role in this process.1 The scenario development process aims to develop a set of new scenarios that facilitate integrated analysis of climate change across the main scientific communities. The process comprises 3 main phases: (1) an initial phase, developing a set of pathways for emissions, concentrations and radiative forcing, (2) a parallel phase, comprising both the development of new socio-economic storylines and climate model projections, and (3) an integration phase, combining the information from the first phases into holistic mitigation, impacts and vulnerability assessments. The pathways developed in the first phase were called 'Representative Concentration Pathways (RCPs)'. They play an important role in providing input for prospective climate model experiments, including both the decadal and long-term projections of climate change. The RCPs also provide an important reference point for new research within the integrated assessment modeling (IAM) community by standardizing on a common set of year-2100 conditions, and exploring alternative pathways and policies that could produce these outcomes. By design, the RCPs, as a set, cover the range of radiative forcing levels examined in the open literature and contain relevant information for climate model runs. This Special Issue documents the main assumptions and characteristics of the RCPs, and, in particular, the various steps that were involved in their development. A number of collaborative activities were initiated and finalized during the last 2-3 years to develop the RCPs. This required the cooperation of researchers from various disciplines involved in climate research, including emission experts, climate modelers, atmospheric chemistry modelers, land use modelers and experts involved in integrated assessment. The four RCPs together reflect the range of year-2100 radiative forcing values found in the literature, i.e. from 2.6 to 8.5 W/m{sup 2}. The papers in this Special Issue describe the individual RCPs, but also the various integrative steps that were necessary within the RCP development process to provide a harmonized set of pathways, that show a smooth transition from the past and extend far into the future for very long-term experiments. Important outcomes of this process included, for instance, the development of new emission inventories, new methods for the harmonization of spatial land use patterns, as well as extensions of the RCP trends beyond 2100. They briefly discuss the content of the individual papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

A WELL-POSED KELVIN-HELMHOLTZ INSTABILITY TEST AND COMPARISON

Recently, there has been a significant level of discussion of the correct treatment of Kelvin-Helmholtz instability (KHI) in the astrophysical community. This discussion relies largely on how the KHI test is posed and analyzed. We pose a stringent test of the initial growth of the instability. The goal is to provide a rigorous methodology for verifying a code on two-dimensional KHI. We ran the problem in the Pencil Code, Athena, Enzo, NDSPMHD, and Phurbas. A strict comparison, judgment, or ranking, between codes is beyond the scope of this work, though this work provides the mathematical framework needed for such a study. Nonetheless, how the test is posed circumvents the issues raised by tests starting from a sharp contact discontinuity yet it still shows the poor performance of smoothed particle hydrodynamics (SPH). We then comment on the connection between this behavior to the underlying lack of zeroth-order consistency in SPH interpolation. We comment on the tendency of some methods, particularly those with very low numerical diffusion, to produce secondary Kelvin-Helmholtz billows on similar tests. Though the lack of a fixed, physical diffusive scale in the Euler equations lies at the root of the issue, we suggest that in some methods an extra diffusion operator should be used to damp the growth of instabilities arising from grid noise. This statement applies particularly to moving-mesh tessellation codes, but also to fixed-grid Godunov schemes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

The expanded LaGrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrainted optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗
NASA NTRS1988

The expanded Lagrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrained optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗
NASA NTRS1993

Improved Blending-Function Algebraic Generation Of Grids

Method of algebraic generation of grids used in computing flows incorporates blending-function techniques to obtain high degree of orthogonality of each grid on its bounding surfaces and to blend grid smoothly, making it conform to all boundaries. Demonstrated in two dimensions, and applied in continuing development of software for real-time interactive generation of three-dimensional grids with color graphical displays. Software used routinely in generating grids for computing flows of practical interest.

Ching-Chung Luh, Raymond↗
NASA NTRS2001

Remote Sensing of Non-Aerosol (anomalous) Absorption in Cloud Free Atmosphere

The interaction of sunlight with atmospheric gases, aerosols and clouds is fundamental to the understanding of climate and its variation. Several studies questioned our understanding of atmospheric absorption of sunlight in cloudy or in cloud free atmospheres. Uncertainty in instruments' accuracy and in the analysis methods makes this problem difficult to resolve. Here we use several years of measurements of sky and sun spectral brightness by selected instruments of the Aerosol Robotic Network (AERONET), that have known and high measurement accuracy. The measurements taken in several locations around the world show that in the atmospheric windows 0.44, 0.06, 0.86 and 1.02 microns the only significant absorbers in cloud free atmosphere is aerosol and ozone. This conclusions is reached using a method developed to distinguish between absorption associated with the presence of aerosol and absorption that is not related to the presence of aerosol. Non-aerosol absorption, defined as spectrally independent or smoothly variable, was found to have an optical thickness smaller than 0.002 corresponding to absorption of sunlight less than 1W/sq m, or essentially zero.

Kaufman, Yoram J.↗
NASA NTRS2003

Wind Tunnel Database Development using Modern Experiment Design and Multivariate Orthogonal Functions

A wind tunnel experiment for characterizing the aerodynamic and propulsion forces and moments acting on a research model airplane is described. The model airplane called the Free-flying Airplane for Sub-scale Experimental Research (FASER), is a modified off-the-shelf radio-controlled model airplane, with 7 ft wingspan, a tractor propeller driven by an electric motor, and aerobatic capability. FASER was tested in the NASA Langley 12-foot Low-Speed Wind Tunnel, using a combination of traditional sweeps and modern experiment design. Power level was included as an independent variable in the wind tunnel test, to allow characterization of power effects on aerodynamic forces and moments. A modeling technique that employs multivariate orthogonal functions was used to develop accurate analytic models for the aerodynamic and propulsion force and moment coefficient dependencies from the wind tunnel data. Efficient methods for generating orthogonal modeling functions, expanding the orthogonal modeling functions in terms of ordinary polynomial functions, and analytical orthogonal blocking were developed and discussed. The resulting models comprise a set of smooth, differentiable functions for the non-dimensional aerodynamic force and moment coefficients in terms of ordinary polynomials in the independent variables, suitable for nonlinear aircraft simulation.

Morelli, Eugene A.↗
NASA NTRS1998

Method and apparatus for reducing the drag of flows over surfaces

An apparatus, and its accompanying method, for reducing the drag of flows over a surface includes arrays of small disks and sensors. The arrays are embedded in the surface and may extend above, or be depressed below, the surface, provided they remain hydraulically smooth either when operating or when inactive. The disks are arranged in arrays of various shapes, and spaced according to the cruising speed of the vehicle on which the arrays are installed. For drag reduction at speeds of the order of 30 meters/second, preferred embodiments include disks that are 0.2 millimeter in diameter and spaced 0.4 millimeter apart. For drag reduction at speeds of the order of 300 meters/second, preferred embodiments include disks that are 0.045 millimeter in diameter and spaced 0.09 millimeter apart. Smaller and larger dimensions for diameter and spacing are also possible. The disks rotate in the plane of the surface, with their rotation axis substantially perpendicular to the surface. The rotating disks produce velocity perturbations parallel to the surface in the overlying boundary layer. The sensors sense the flow at the surface and connect to control circuitry that adjusts the rotation rates and duty cycles of the disks accordingly. Suction and blowing holes can be interspersed among, or made coaxial with, the disks for creating general three-component velocity perturbations in the near-surface region. The surface can be a flat, planar surface or a nonplanar surface, such as a triangular riblet surface. The present apparatus and method have potential applications in the field of aeronautics for improving performance and efficiency of commercial and military aircraft, and in other industries where drag is an obstacle, including gas and oil delivery through long-haul pipelines.

Keefe, Laurence R.↗
NASA NTRS1986

Proposed method of producing large optical mirrors Single-point diamond crushing followed by polishing with a small-area tool

Faster production of large optical mirrors may result from combining single-point diamond crushing of the glass with polishing using a small area tool to smooth the surface and remove the damaged layer. Diamond crushing allows a surface contour accurate to 0.5 microns to be generated, and the small area computer-controlled polishing tool allows the surface roughness to be removed without destroying the initial contour. Final contours with an accuracy of 0.04 microns have been achieved.

Wright, G.↗
NASA NTRS1989

Application of a self-adaptive grid method to complex flows

A directional-split modular user-friendly grid-point distribution code is applied to several test problems. Grid-point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flowfield gradients and/or geometric features, are used to determine the point distribution.

Deiwert, G. S.↗
NASA NTRS2017

LPV Modeling of a Flexible Wing Aircraft Using Modal Alignment and Adaptive Gridding Methods

One of the earliest approaches in gain-scheduling control is the gridding based approach, in which a set of local linear time-invariant models are obtained at various gridded points corresponding to the varying parameters within the flight envelop. In order to ensure smooth and effective Linear Parameter-Varying control, aligning all the flexible modes within each local model and maintaining small number of representative local models over the gridded parameter space are crucial. In addition, since the flexible structural models tend to have large dimensions, a tractable model reduction process is necessary. In this paper, the notion of s-shifted H2- and H Infinity-norm are introduced and used as a metric to measure the model mismatch. A new modal alignment algorithm is developed which utilizes the defined metric for aligning all the local models over the entire gridded parameter space. Furthermore, an Adaptive Grid Step Size Determination algorithm is developed to minimize the number of local models required to represent the gridded parameter space. For model reduction, we propose to utilize the concept of Composite Modal Cost Analysis, through which the collective contribution of each flexible mode is computed and ranked. Therefore, a reduced-order model is constructed by retaining only those modes with significant contribution. The NASA Generic Transport Model operating at various flight speeds is studied for verification purpose, and the analysis and simulation results demonstrate the effectiveness of the proposed modeling approach.

LPV Modeling↗
DOE OSTI2019

Estimating PV power from aggregate power measurements within the distribution grid

The increased integration of photovoltaic (PV) systems in distribution grids reduces visibility and situational awareness for utilities because the PV systems' power production is usually not monitored by them. To address this problem, a method called Contextually Supervised Source Separation (CSSS) has been recently adapted for real-time estimation of aggregate PV active power generation from aggregate net active and reactive power measurements at a point in a radially configured distribution grid (e.g., substation). In its original version, PV disaggregation is formulated as an optimization problem that fits linear regression models for the aggregate PV active power generation and true substation active power load. This paper extends the previous work by adding regularization terms in the objective function to capture additional contextual information such as smoothness, by adding new constraints, by introducing new regressors such as ambient temperature, and by investigating the use of time-varying regressors. Furthermore, we perform extensive parametric analysis to study tuning of the objective function weighting factors in a way that maximizes performance and robustness. The proposed PV disaggregation method can be applied to networks with either a single PV system (e.g., megawatt scale) or many distributed ones (e.g., residential scale) connected downstream of the substation. Simulation studies with real field recorded data show that the enhancements of the proposed method reduce disaggregation error by 58% in winter and 35% in summer compared to previous CSSS-based work. When compared against a commonly used transposition model based approach, the reduction in disaggregation error is more pronounced (78% reduction in winter and 45% in summer). Additional simulations indicate that the proposed algorithm is also applicable for PV systems with time-varying power factors. In conclusion, overall, our results show that—with appropriate modeling and tuning—it is possible to accurately estimate the aggregated PV active power generation of a distribution feeder with minimal or no additional sensor deployment.

14 SOLAR ENERGY↗
DOE OSTI2012

The second-order differential phase contrast and its retrieval for imaging with x-ray Talbot interferometry

Purpose: The x-ray differential phase contrast imaging implemented with the Talbot interferometry has recently been reported to be capable of providing tomographic images corresponding to attenuation-contrast, phase-contrast, and dark-field contrast, simultaneously, from a single set of projection data. The authors believe that, along with small-angle x-ray scattering, the second-order phase derivative {Phi}{sup Double-Prime }{sub s}(x) plays a role in the generation of dark-field contrast. In this paper, the authors derive the analytic formulae to characterize the contribution made by the second-order phase derivative to the dark-field contrast (namely, second-order differential phase contrast) and validate them via computer simulation study. By proposing a practical retrieval method, the authors investigate the potential of second-order differential phase contrast imaging for extensive applications. Methods: The theoretical derivation starts at assuming that the refractive index decrement of an object can be decomposed into {delta}={delta}{sub s}+{delta}{sub f}, where {delta}{sub f} corresponds to the object's fine structures and manifests itself in the dark-field contrast via small-angle scattering. Based on the paraxial Fresnel-Kirchhoff theory, the analytic formulae to characterize the contribution made by {delta}{sub s}, which corresponds to the object's smooth structures, to the dark-field contrast are derived. Through computer simulation with specially designed numerical phantoms, an x-ray differential phase contrast imaging system implemented with the Talbot interferometry is utilized to evaluate and validate the derived formulae. The same imaging system is also utilized to evaluate and verify the capability of the proposed method to retrieve the second-order differential phase contrast for imaging, as well as its robustness over the dimension of detector cell and the number of steps in grating shifting. Results: Both analytic formulae and computer simulations show that, in addition to small-angle scattering, the contrast generated by the second-order derivative is magnified substantially by the ratio of detector cell dimension over grating period, which plays a significant role in dark-field imaging implemented with the Talbot interferometry. Conclusions: The analytic formulae derived in this work to characterize the second-order differential phase contrast in the dark-field imaging implemented with the Talbot interferometry are of significance, which may initiate more activities in the research and development of x-ray differential phase contrast imaging for extensive preclinical and eventually clinical applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1997

Exponential Approximations Using Fourier Series Partial Sums

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Banerjee, Nana S.↗
NASA NTRS1973

Experimental and data analysis techniques for deducing collision-induced forces from photographic histories of engine rotor fragment impact/interaction with a containment ring

An analysis method termed TEJ-JET is described whereby measured transient elastic and inelastic deformations of an engine-rotor fragment-impacted structural ring are analyzed to deduce the transient external forces experienced by that ring as a result of fragment impact and interaction with the ring. Although the theoretical feasibility of the TEJ-JET concept was established, its practical feasibility when utilizing experimental measurements of limited precision and accuracy remains to be established. The experimental equipment and the techniques (high-speed motion photography) employed to measure the transient deformations of fragment-impacted rings are described. Sources of error and data uncertainties are identified. Techniques employed to reduce data reading uncertainties and to correct the data for optical-distortion effects are discussed. These procedures, including spatial smoothing of the deformed ring shape by Fourier series and timewise smoothing by Gram polynomials, are applied illustratively to recent measurements involving the impact of a single T58 turbine rotor blade against an aluminum containment ring. Plausible predictions of the fragment-ring impact/interaction forces are obtained by one branch of this TEJ-JET method; however, a second branch of this method, which provides an independent estimate of these forces, remains to be evaluated.

Yeghiayan, R. P.↗