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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

An embedded function approach for turbulent flow prediction

In conventional prediction methods for turbulent flow influence on surface properties, very small mesh sizes and considerable computational effort is required to adequately resolve the intense velocity and temperature profile variations that occur in the wall-layer region. In this study, an approach is described wherein an outer region numerical solution is smoothly matched to a set of embedded analytic profile functions in the near-wall region; these wall-layer profile functions having been derived through consideration of the coherent structure of the time-dependent near-wall flow provide analytic expressions for the mean velocity and enthalpy profiles in the wall layer. The concept is demonstrated through example applications to turbulent boundary-layer flows in two dimensions. The technique is shown to be very efficient and it is demonstrated that a reduction of approximately half the mesh points across the 2-D layer may be realized (as compared to conventional methods) with no degradation in accuracy.

Walker, J. D. A.↗
DOE OSTI2021

Manifold Sampling for Optimizing Nonsmooth Nonconvex Compositions

Here we propose a manifold sampling algorithm for minimizing a nonsmooth composition $f= h\circ F$, where we assume $h$ is nonsmooth and may be inexpensively computed in closed form and $F$ is smooth but its Jacobian may not be available. We additionally assume that the composition $h\circ F$ defines a continuous selection. Manifold sampling algorithms can be classified as model-based derivative-free methods, in that models of $F$ are combined with particularly sampled information about $h$ to yield local models for use within a trust-region framework. We demonstrate that cluster points of the sequence of iterates generated by the manifold sampling algorithm are Clarke stationary. We consider the tractability of three particular subproblems generated by the manifold sampling algorithm and the extent to which inexact solutions to these subproblems may be tolerated. Numerical results demonstrate that manifold sampling as a derivative-free algorithm is competitive with state-of-the-art algorithms for nonsmooth optimization that utilize first-order information about $f$.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Extracting the Breakout Distance from the ECOT Trajectories: Gaussian Process Regression Approach

Enhanced Corner Turning (ECOT) experiments provide an important metric of performance of high explosive (HE) formulations. The breakout distance is a single scalar value that characterizes the corner turning efficiency of an HE. Extracting the breakout distance from the raw ECOT results, whether experimental or simulated, is a conceptually straightforward procedure which, however, is non-unique, especially in the presence of noise. More specifically, this procedure involves numerical smoothing and selecting particular values for parameters of this smoothing introduces human bias. In this work, we propose to use the Gaussian process regression to analyze ECOT results. This analysis involves the effective smoothing of the data, thus allowing for accurate extraction of the breakout distance. Most importantly, the parameters of this smoothing can be inferred from the ECOT data itself, rendering the approach effectively parameter-free and thus diminishing the human bias. An additional benefit of the Gaussian process regression, being a statistical inference method, is that not just the value of the breakout distance, but also its confidence interval can be extracted from the data. This report introduces the Gaussian process regression, as applied to ECOT, and demonstrates its usefulness by extracting the breakout distances for a selection of experimental and simulated data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

An Improved Method to Measure the Cosmic Curvature

In this paper, we propose an improved model-independent method to constrain the cosmic curvature by combining the most recent Hubble parameter H ( z ) and supernovae Ia (SNe Ia) data. Based on the H ( z ) data, we first use the model-independent smoothing technique, Gaussian processes, to construct a distance modulus μ {sub H} ( z ), which is susceptible to the cosmic curvature parameter Ω{sub k}. In contrary to previous studies, the light-curve-fitting parameters, which account for the distance estimation of SN (μ {sub SN}( z )), are set free to investigate whether Ω {sub k} has a dependence on them. By comparing μ {sub H} ( z ) to μ {sub SN}(z), we put limits on Ω {sub k}. Our results confirm that Ω {sub k} is independent of the SN light-curve parameters. Moreover, we show that the measured Ω {sub k} is in good agreement with zero cosmic curvature, implying that there is no significant deviation from a flat universe at the current observational data level. We also test the influence of different H(z) samples and different Hubble constant H {sub 0} values, finding that different H(z) samples do not have a significant impact on the constraints. However, different H {sub 0} priors can affect the constraints of Ω {sub k} to some degree. The prior of H {sub 0} = 73.24 ± 1.74 km s{sup −1} Mpc{sup −1} gives a value of Ω {sub k}, a little bit above the 1 σ confidence level away from 0, but H{sub 0} = 69.6 ± 0.7 km s{sup −1} Mpc{sup −1} gives it below 1 σ .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1981

New approach to the solution of large, full matrix equations

A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.

Clark, R. W.↗
NASA NTRS1980

Navier-Stokes calculations with a coupled strongly implicit method. II Spline deferred-corrector solutions

The coupled strongly implicit (CSIP) method described previously is combined with a deferred-corrector spline solver for the vorticity-stream function form of the Navier-Stokes equations. Solutions for cavity, channel and cylinder flows are obtained with the fourth-order spline 4 procedure. The strongly coupled spline corrector method converges as rapidly as the finite difference calculations and also allows for arbitrary large time increments for the Reynolds numbers considered (equal to or less than 1000). In some cases fourth-order smoothing or filtering is required in order to suppress high frequency oscillations.

Rubin, S. G.↗
NASA NTRS1990

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗
NASA NTRS1991

Stability and bifurcation analysis of oscillators with piecewise-linear characteristics - A general approach

A general approach is developed for determining the periodic solutions and their stability of nonlinear oscillators with piecewise-smooth characteristics. A modified harmonic balance/Fourier transform procedure is devised for the analysis. The procedure avoids certain numerical differentiation employed previously in determining the periodic solutions, therefore enhancing the reliability and efficiency of the method. Stability of the solutions is determined via perturbations of their state variables. The method is applied to a forced oscillator interacting with a stop of finite stiffness. Flip and fold bifurcations are found to occur. This led to the identification of parameter ranges in which chaotic response occurred.

Noah, S. T.↗
NASA NTRS2005

A Forecast of Reduced Solar Activity and Its Implications for NASA

The "Solar Dynamo" method of solar activity forecasting is reviewed. Known generically as a 'precursor" method, insofar as it uses observations which precede solar activity generation, this method now uses the Solar Dynamo Amplitude (SODA) Index to estimate future long-term solar activity. The peak amplitude of the next solar cycle (#24), is estimated at roughly 124 in terms of smoothed F10.7 Radio Flux and 74 in terms of the older, more traditional smoothed international or Zurich Sunspot number (Ri or Rz). These values are significantly smaller than the amplitudes of recent solar cycles. Levels of activity stay large for about four years near the peak in smoothed activity, which is estimated to occur near the 2012 timeflame. Confidence is added to the prediction of low activity by numerous examinations of the Sun s weakened polar field. Direct measurements are obtained by the Mount Wilson Solar Observatory and the Wilcox Solar Observatory. Further support is obtained by examining the Sun s polar faculae (bright features), the shape of coronal soft X-ray "holes," and the shape of the "source surface" - a calculated coronal feature which maps the large scale structure of the Sun s field. These features do not show the characteristics of well-formed polar coronal holes associated with typical solar minima. They show stunted polar field levels, which are thought to result in stunted levels of solar activity during solar cycle #24. The reduced levels of solar activity would have concomitant effects upon the space environment in which satellites orbit. In particular, the largest influences would affect orbit determination of satellites in LEO (Low Earth Orbit), based upon the altered thermospheric and exospheric densities. A decrease in solar activity would result in smaller satellite decay rates, as well as fewer large solar events that can destroy satellite electronic functions. Other effects of reduced solar activity upon the space environment include enhanced galactic cosmic rays and more space debris at low altitudes (from the decay of old satellite parts, etc.). The reasons are well known: namely, solar activity serves to sweep the inner heliosphere of galactic cosmic rays, and lower exospheric densities result in decreased drag on LEO debris, allowing longer lifetimes.

Schatten, Kenneth↗
NASA NTRS2007

Building Operations Efficiencies into NASA's Ares I Crew Launch Vehicle Design

The U.S. Vision for Space Exploration guides the National Aeronautics and Space Administration's (NASA's) challenging missions that expand humanity's boundaries and open new routes to the space frontier. With the Agency's commitment to complete the International Space Station (ISS) and to retire the venerable Space Shuttle by 2010, the NASA Administrator commissioned the Exploration Systems Architecture Study (ESAS) in 2005 to analyze options for safe, simple, cost-efficient launch solutions that could deliver human-rated space transportation capabilities in a timely manner within fixed budget guidelines. The Exploration Launch Projects (ELP) Office, chartered by the Constellation Program in October 2005, has been conducting systems engineering studies and business planning to successively refine the design configurations and better align vehicle concepts with customer and stakeholder requirements, such as significantly reduced life-cycle costs. As the Agency begins the process of replacing the Shuttle with a new generation of spacecraft destined for missions beyond low-Earth orbit to the Moon and Mars, NASA is designing the follow-on crew and cargo launch systems for maximum operational efficiencies. To sustain the long-term exploration of space, it is imperative to reduce the $4 billion NASA typically spends on space transportation each year. This paper gives toplevel information about how the follow-on Ares I Crew Launch Vehicle (CLV) is being designed for improved safety and reliability, coupled with reduced operations costs. These methods include carefully developing operational requirements; conducting operability design and analysis; using the latest information technology tools to design and simulate the vehicle; and developing a learning culture across the workforce to ensure a smooth transition between Space Shuttle operations and Ares vehicle development.

Dumbacher, Daniel L.↗
DOE OSTI2016

Particle control near reticle and optics using showerhead

A method and an apparatus to protect a reticle against particles and chemicals in an actinic EUV reticle inspection tool are presented. The method and apparatus utilizes a pair of porous metal diffusers in the form of showerheads to provide a continual flow of clean gas. The main showerhead bathes the reticle surface to be inspected in smoothly flowing, low pressure gas, isolating it from particles coming from surrounding volumes. The secondary showerhead faces away from the reticle and toward the EUV illumination and projection optics, supplying them with purge gas while at the same time creating a buffer zone that is kept free of any particle contamination originating from those optics.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Operator split methods in the numerical solution of the finite deformation elastoplastic dynamic problem

The spatial formulation of the elastoplastic dynamic problem for finite deformations is considered. A thermodynamic argument leads to an additive decomposition of the spatial rate of deformation tensor and allows an operator split of the evolutionary equations of the problem into elastic and plastic parts. This operator split is taken as the basis for the definition of a global product algorithm. In the context of finite element discretization the product algorithm entails, for every time step, the solution of a nonlinear elastodynamic problem followed by the application of plastic algorithms that operate on the stresses and internal variables at the integration points and bring in the plastic constitutive equations. Suitable plastic algorithms are discussed for the cases of perfect and hardening plasticity and viscoplasticity. The proposed formalism does not depend on any notion of smoothness of the yield surface and is applicable to arbitrary convex elastic regions, with or without corners. The stabiity properties of the global product algorithm are shown to be identical to those of the algorithm used for the integration of the nonlinear elastodynamic problem. Numerical examples illustrate the accuracy of the method.

Pinsky, P. M.↗
DOE OSTI2024

Using spatio-temporal graph neural networks to estimate fleet-wide photovoltaic performance degradation patterns

Accurate estimation of photovoltaic (PV) system performance is crucial for determining its feasibility as a power generation technology and financial asset. PV-based energy solutions offer a viable alternative to traditional energy resources due to their superior Levelized Cost of Energy (LCOE). A significant challenge in assessing the LCOE of PV systems lies in understanding the Performance Loss Rate (PLR) for large fleets of PV systems. Estimating the PLR of PV systems becomes increasingly important in the rapidly growing PV industry. Precise PLR estimation benefits PV users by providing real-time monitoring of PV module performance, while explainable PLR estimation assists PV manufacturers in studying and enhancing the performance of their products. However, traditional PLR estimation methods based on statistical models have notable drawbacks. Firstly, they require user knowledge and decision-making. Secondly, they fail to leverage spatial coherence for fleet-level analysis. Additionally, these methods inherently assume the linearity of degradation, which is not representative of real world degradation. To overcome these challenges, we propose a novel graph deep learning-based decomposition method called the Spatio-Temporal Graph Neural Network for fleet-level PLR estimation (PV-stGNN-PLR). PV-stGNN-PLR decomposes the power timeseries data into aging and fluctuation components, utilizing the aging component to estimate PLR. PV-stGNN-PLR exploits spatial and temporal coherence to derive PLR estimation for all systems in a fleet and imposes flatness and smoothness regularization in loss function to ensure the successful disentanglement between aging and fluctuation. We have evaluated PV-stGNN-PLR on three simulated PV datasets consisting of 100 inverters from 5 sites. Experimental results show that PV-stGNN-PLR obtains a reduction of 33.9% and 35.1% on average in Mean Absolute Percent Error (MAPE) and Euclidean Distance (ED) in PLR degradation pattern estimation compared to the state-of-the-art PLR estimation methods.

14 SOLAR ENERGY↗
DOE OSTI2021

Combining Spectroscopic Ellipsometry with Transmission Spectroscopy to Derive Accurate Optical Constants for Organic Materials

A method for deriving the optical constants of organic powdered materials in the mid-infrared spectral range is introduced using variable angle spectroscopic ellipsometry and transmission spectroscopy. The approach uses pressed pellets of the powder and is applied to organic solids, which have both strong and weak infrared absorption features. Many powders have significant voids and do not press into smooth, homogenous pellets. To account for pellet non-idealities and to accurately measure both n and k, three different forms of pellets were pressed and measured: A pure analyte pellet, a mixed analyte/KBr pellet with a large analyte percentage, and a KBr transmission pellet with only a small analyte percentage. Using all three pellets in a multi-sample analysis involving both ellipsometric and transmission data, the complex refractive index (n/k) values can be derived for many organic compounds. This method is illustrated to calculate the optical constants for anhydrous lactose from 6000-400 cm-1. The transmission measurements improve the spectral fitting of weak absorption features, and the multi-sample analysis enables a better determination of the significant void space that is present in the pure pellet, leading to lower values for both n and k if not properly accounted for in the multi-oscillator model used to fit the ellipsometric data.

Optical constants, refractive index, extinction co↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗
DOE OSTI2011

DENSITY ESTIMATION FOR PROJECTED EXOPLANET QUANTITIES

Exoplanet searches using radial velocity (RV) and microlensing (ML) produce samples of 'projected' mass and orbital radius, respectively. We present a new method for estimating the probability density distribution (density) of the unprojected quantity from such samples. For a sample of n data values, the method involves solving n simultaneous linear equations to determine the weights of delta functions for the raw, unsmoothed density of the unprojected quantity that cause the associated cumulative distribution function (CDF) of the projected quantity to exactly reproduce the empirical CDF of the sample at the locations of the n data values. We smooth the raw density using nonparametric kernel density estimation with a normal kernel of bandwidth {sigma}. We calibrate the dependence of {sigma} on n by Monte Carlo experiments performed on samples drawn from a theoretical density, in which the integrated square error is minimized. We scale this calibration to the ranges of real RV samples using the Normal Reference Rule. The resolution and amplitude accuracy of the estimated density improve with n. For typical RV and ML samples, we expect the fractional noise at the PDF peak to be approximately 80 n{sup -log2}. For illustrations, we apply the new method to 67 RV values given a similar treatment by Jorissen et al. in 2001, and to the 308 RV values listed at exoplanets.org on 2010 October 20. In addition to analyzing observational results, our methods can be used to develop measurement requirements-particularly on the minimum sample size n-for future programs, such as the microlensing survey of Earth-like exoplanets recommended by the Astro 2010 committee.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Inhibition of cell proliferation by a selective inhibitor of the Ca{sup 2+}-activated Cl{sup -} channel, Ano1

Highlights: Black-Right-Pointing-Pointer T16A{sub inh}-A01 blocked Ano1 currents in HEK cells expressing Ano1. Black-Right-Pointing-Pointer T16A{sub inh}-A01 reduced proliferation in ICC primary cultures and CFPAC-1 cell line. Black-Right-Pointing-Pointer T16A{sub inh}-A01 reduced proliferation of ICC in intact smooth muscle strips. -- Abstract: Background: Ion channels play important roles in regulation of cellular proliferation. Ano1 (TMEM16A) is a Ca{sup 2+}-activated Cl{sup -} channel expressed in several tumors and cell types. In the muscle layers of the gastrointestinal tract Ano1 is selectively expressed in interstitial cells of Cajal (ICC) and appears to be required for normal gastrointestinal slow wave electrical activity. However, Ano1 is expressed in all classes of ICC, including those that do not generate slow waves suggesting that Ano1 may have other functions. Indeed, a role for Ano1 in regulating proliferation of tumors and ICC has been recently suggested. Recently, a high-throughput screen identified a small molecule, T16A{sub inh}-A01 as a specific inhibitor of Ano1. Aim: To investigate the effect of the T16A{sub inh}-A01 inhibitor on proliferation in ICC and in the Ano1-expressing human pancreatic cancer cell line CFPAC-1. Methods: Inhibition of Ano1 was demonstrated by whole cell voltage clamp recordings of currents in cells transfected with full-length human Ano1. The effect of T16A{sub inh}-A01 on ICC proliferation was examined in situ in organotypic cultures of intact mouse small intestinal smooth muscle strips and in primary cell cultures prepared from these tissues. ICC were identified by Kit immunoreactivity. Proliferating ICC and CFPAC-1 cells were identified by immunoreactivity for the nuclear antigen Ki67 or EdU incorporation, respectively. Results: T16A{sub inh}-A01 inhibited Ca{sup 2+}-activated Cl{sup -} currents by 60% at 10 {mu}M in a voltage-independent fashion. Proliferation of ICC was significantly reduced in primary cultures from BALB/c mice following treatment with T16A{sub inh}-A01. Proliferation of the CFPAC-1 human cell-line was also reduced by T16A{sub inh}-A01. In organotypic cultures of smooth muscle strips from mouse jejunum, the proliferation of ICC was reduced but the total number of proliferating cells/confocal stack was not affected, suggesting that the inhibitory effect was specific for ICC. Conclusions: The selective Ano1 inhibitor T16A{sub inh}-A01 inhibited Ca{sup 2+}-activated Cl{sup -} currents, reduced the number of proliferating ICC in culture and inhibited proliferation in the pancreatic cancer cell line CFPAC-1. These data support the notion that chloride channels in general and Ano1 in particular are involved in the regulation of proliferation.

60 APPLIED LIFE SCIENCES↗