Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

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,261 records · Page 70

Reducing the number of variational equations in the implementation of multiple shooting

The standard method of multiple shooting for a system of n first-order differential equations with k unknown initial conditions requires the integration of k sets of variational equations on the first shot and n sets of variational equations on every shot thereafter. This paper describes a variant of multiple shooting that requires the solution of k sets of variational equations on every shot. The technique applies to both linear and nonlinear boundary-value problems. Techniques to deal with difficulties unique to the solution of nonlinear problems are suggested.

Krogh, F. T.↗

Strain Gauges Indicate Differential-CTE-Induced Failures

A method of detecting mechanical failure induced by variation in temperature at an adhesive bond between two materials that have different coefficients of thermal expansion (CTEs) involves monitoring of strain-gauge readings. This method can be regarded as an exploitation of the prior observation that the readings of strain gauges commonly used in tensile and compressive testing of material specimens include features indicative of incremental failures in the specimens. In this method, one or more strain gauges are bonded to either or both of the two materials near the bond between the materials. (The adhesive used to bond the strain gauges would not ordinarily be the same as the one used to bond the two materials). Then strain-gauge readings are recorded as the temperature of the materials is varied through a range of interest. Any significant discontinuity in the slope of the resulting strain-versus-temperature curve(s) is taken to be a qualitative indication of a failure of the bond between the two materials and/or a failure within one of the materials in the vicinity of the bond. The method has been demonstrated in experiments on specimens consisting of polyacrylonitrile-fiber/epoxy-matrix laminated composite plates bonded by epoxy to smaller plates made, variously, of aluminum, titanium, and a low-CTE nickel/iron alloy. In preparation for each experiment, strain gauges were bonded, by use of cryogenic-rated adhesives, to the composite plate near the corners of the metal plate (see Figure 1). In each experiment, strain-gauge and temperature readings were taken as the specimen was cooled from room temperature to 20 K. The specimens were then returned to room temperature and ultrasonically inspected for damage in the bond region. No failure events were detectable in the strain-gauge readings from the composite/ titanium and composite/low-thermalexpansion- alloy specimens, and ultrasonic inspection of these specimens revealed no damage. However, failure events were seen in the strain-gauge readings from the composite/aluminum specimens (see Figure 2), and ultrasonic inspection confirmed that there was damage in the bond regions of these specimens.

Harris, Brian↗

Algorithm 1049: The Delaunay Density Diagnostic

Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is sufficient to identify the relevant variations in outputs—i.e., the “geometric scale” of the function—despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this article, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. The algorithm is intended to work on scalar data in moderate dimensions (2–10). Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in L 2 norm) if and only if the data are sampled densely enough to distinguish the feature from noise (assuming sufficiently regular sampling). We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e., static) datasets of input–output pairs. Finally, we include analytical results in support of our numerical findings and have released lightweight code that can be adapted for use in a variety of data science settings.

97 MATHEMATICS AND COMPUTING↗

Influence of Design Variations on Systems Performance

High-risk aerospace components have to meet very stringent quality, performance, and safety requirements. Any source of variation is a concern, as it may result in scrap or rework. poor performance, and potentially unsafe flying conditions. The sources of variation during product development, including design, manufacturing, and assembly, and during operation are shown. Sources of static and dynamic variation during development need to be detected accurately in order to prevent failure when the components are placed in operation. The Systems' Health and Safety (SHAS) research at the NASA Ames Research Center addresses the problem of detecting and evaluating the statistical variation in helicopter transmissions. In this work, we focus on the variations caused by design, manufacturing, and assembly of these components, prior to being placed in operation (DMV). In particular, we aim to understand and represent the failure and variation information, and their correlation to performance and safety and feed this information back into the development cycle at an early stage. The feedback of such critical information will assure the development of more reliable components with less rework and scrap. Variations during design and manufacturing are a common source of concern in the development and production of such components. Accounting for these variations, especially those that have the potential to affect performance, is accomplished in a variety ways, including Taguchi methods, FMEA, quality control, statistical process control, and variation risk management. In this work, we start with the assumption that any of these variations can be represented mathematically, and accounted for by using analytical tools incorporating these mathematical representations. In this paper, we concentrate on variations that are introduced during design. Variations introduced during manufacturing are investigated in parallel work.

Tumer, Irem Y.↗

Completed Beltrami-Michell Formulation in Polar Coordinates

A set of conditions had not been formulated on the boundary of an elastic continuum since the time of Saint-Venant. This limitation prevented the formulation of a direct stress calculation method in elasticity for a continuum with a displacement boundary condition. The missed condition, referred to as the boundary compatibility condition, is now formulated in polar coordinates. The augmentation of the new condition completes the Beltrami-Michell formulation in polar coordinates. The completed formulation that includes equilibrium equations and a compatibility condition in the field as well as the traction and boundary compatibility condition is derived from the stationary condition of the variational functional of the integrated force method. The new method is illustrated by solving an example of a mixed boundary value problem for mechanical as well as thermal loads.

Patnaik, Surya N.↗

Borehole muography of subsurface reservoirs

Imaging subsurface rock formations or geological objects like oil and gas reservoirs, mineral deposits, cavities or even magmatic plumbing systems under active volcanoes has been for many years a major quest of geoscientists. Since these subsurface objects cannot be observed directly, different indirect methods have been developed. These methods are all based on variations of certain physical properties of the subsurface materials that can be detected from the ground surface or from boreholes. To determine the density distribution, a new imaging technique using cosmic-ray muon detectors deployed in a borehole has been developed and a first prototype of a borehole muon detector successfully tested. In addition to providing a static image of the subsurface density in three dimensions (or three-dimensional tomography), borehole muography can also inform on the variations of density with time, which recently became of major importance with the injection of large volumes of fluids, mainly water and CO2, in porous subsurface reservoirs (e.g. aquifer storage and recovery, wastewater disposal, enhanced oil recovery and carbon sequestration). This raises several concerns about the risk of leakage and the mechanical integrity of the reservoirs. Determining the field scale induced displacement of fluids by geophysical methods like muography is thus a priority.This article is part of the Theo Murphy meeting issue 'Cosmic-ray muography'.

58 GEOSCIENCES↗

Matrix methods for the design of cascades to prescribed surface velocity distributions and for fully compressible flow

This paper describes matrix methods that have been developed for calculating compressible flow on a blade-to-blade surface of revolution. The methods have been fully tested to date only for the design of plane cascades to prescribed blade surface distributions; the methods will be illustrated here for that problem only. Similar methods are presently being applied to both the direct and indirect problems and for flow on arbitrary surfaces of revolution in annular cascades with stream sheet thickness variations. It is believed that by such methods, both the direct and indirect calculations can be reduced to about 60 to 90 seconds of computing.

Silvester, M. E.↗

Enhanced Separation of Sm and La from Sm{sub 2}O{sub 3}-La{sub 2}O{sub 3} Mixtures Based on a Reduction–Oxidation Method

The Sm{sub 2}O{sub 3}-La{sub 2}O{sub 3} mixture was converted to the component chlorides, and Sm and La were separated from the chlorides using the reduction–oxidation method based on the chemical valence variation of Sm. In the reduction stage, samarium (III) chloride was reduced to samarium (II) chloride using hydrogen. In the oxidation stage, samarium (II) chloride was oxidized by water to SmCl{sub 3} and water-insoluble SmOCl. After filter separation of SmOCl, the Sm{sub 2}O{sub 3} content in the initial La{sub 2}O{sub 3}-Sm{sub 2}O{sub 3} mixture was reduced by ~ 50 pct. These results demonstrate an innovative method for pyrometallurgical separation of Sm and La from Sm{sub 2}O{sub 3}-La{sub 2}O{sub 3} mixtures.

36 MATERIALS SCIENCE↗

Increased signal-to-noise ratios within experimental field trials by regressing spatially distributed soil properties as principal components

Environmental variability poses a major challenge to any field study. Researchers attempt to mitigate this challenge through replication. Thus, the ability to detect experimental signals is determined by the degree of replication and the amount of environmental variation, noise, within the experimental system. A major source of noise in field studies comes from the natural heterogeneity of soil properties which create microtreatments throughout the field. In addition, the variation within different soil properties is often nonrandomly distributed across a field. We explore this challenge through a sorghum field trial dataset with accompanying plant, microbiome, and soil property data. Diverse sorghum genotypes and two watering regimes were applied in a split-plot design. We describe a process of identifying, estimating, and controlling for the effects of spatially distributed soil properties on plant traits and microbial communities using minimal degrees of freedom. Importantly, this process provides a method with which sources of environmental variation in field data can be identified and adjusted, improving our ability to resolve effects of interest and to quantify subtle phenotypes.

60 APPLIED LIFE SCIENCES↗

Errant Beam Prognostics with Machine Leaning at SNS Accelerator

Particle Accelerators are complex machine with many pieces of equipment running in synchronization to deliver required beam. However, faults in particle accelerators reduce the availability of the beam for experiments affecting the overall science output. To avoid these faults, we apply anomaly detection techniques to predict any unusual behavior and perform preemptive actions to improve the total availability. Many researchers have adopted semi-supervised Machine Learning (ML) methods such as auto-encoders and variational auto-encoders for such tasks. However, supervised ML techniques designed for similarity learning such as Siamese Neural Network (SNN) can outperform semi-supervised or unsupervised methods for anomaly prediction. One of the challenges associated with application of ML models to particle accelerators is the variability in observed data over time due to system configuration changes. We employ conditional models such as Conditional Siamese Neural Networks (CSNN), and Conditional-VAE (CVAE) to learn the variability in the data by using beam configuration parameters as conditional input. We apply these models for errant beam prediction at Spallation Neutron Source accelerator under different system configurations and compare their performance. We demonstrate that CSNN outperforms CVAE in our application. This talk will present the data source, collection, analysis, data-preparation, model development, hyper-parameter studies and the results.

Rajput, Kishansingh↗

Tracking Excited States in Wave Function Optimization Using Density Matrices and Variational Principles

We introduce an approach for finding individual excited states' energy stationary points in complete active space self-consistent field theory that is compatible with standard optimization methods and highly effective at overcoming difficulties due to root flipping and near-degeneracies. Inspired by both the maximum overlap method and recent progress in excited-state variational principles, our approach combines these ideas in order to track individual excited states throughout the orbital optimization process. In a series of tests involving root flipping, near-degeneracies, charge transfers, and double excitations, we show that this approach is more effective for state-specific optimization than either the naive selection of roots on the basis of energy ordering or a more direct generalization of the maximum overlap method. We offer evidence that this state-specific approach improves the performance of complete active space perturbation theory for vertical excitation energies. Furthermore, we find that the state-specific optimization can help avoid state-averaging-induced discontinuities on potential energy surfaces. With a simple implementation, a low cost, and compatibility with large active space methods, the approach is designed to be useful in a wide range of excited-state investigations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Completed Beltrami-Michell Formulation for Analyzing Radially Symmetrical Bodies

A force method formulation, the completed Beltrami-Michell formulation (CBMF), has been developed for analyzing boundary value problems in elastic continua. The CBMF is obtained by augmenting the classical Beltrami-Michell formulation with novel boundary compatibility conditions. It can analyze general elastic continua with stress, displacement, or mixed boundary conditions. The CBMF alleviates the limitations of the classical formulation, which can solve stress boundary value problems only. In this report, the CBMF is specialized for plates and shells. All equations of the CBMF, including the boundary compatibility conditions, are derived from the variational formulation of the integrated force method (IFM). These equations are defined only in terms of stresses. Their solution for kinematically stable elastic continua provides stress fields without any reference to displacements. In addition, a stress function formulation for plates and shells is developed by augmenting the classical Airy's formulation with boundary compatibility conditions expressed in terms of the stress function. The versatility of the CBMF and the augmented stress function formulation is demonstrated through analytical solutions of several mixed boundary value problems. The example problems include a composite circular plate and a composite circular cylindrical shell under the simultaneous actions of mechanical and thermal loads.

Kaljevic, Igor↗

Variational Approach to Monte Carlo Renormalization Group

We present a Monte Carlo method for computing the renormalized coupling constants and the critical exponents within renormalization theory. The scheme, which derives from a variational principle, overcomes critical slowing down, by means of a bias potential that renders the coarse grained variables uncorrelated. Here, the two-dimensional Ising model is used to illustrate the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nitrogen Fertilization Restructured Spatial Patterns of Soil Organic Carbon and Total Nitrogen in Switchgrass and Gamagrass Croplands in Tennessee USA

Nitrogen (N) fertilizers can potentially alter spatial distribution of soil organic carbon (SOC) and total nitrogen (TN) concentrations in croplands such as switchgrass (SG: Panicum virgatum L.) and gamagrass (GG: Tripsacum dactyloides L.), but it remains unclear whether these effects are the same between crops and under different rates of fertilization. 13 C and 15 N are two important proxy measures of soil biogeochemistry, but they were rarely examined as to their spatial distributions in soil. Based on a three-year long fertilization experiment in Middle Tennessee, USA, the top mineral horizon soils (0–15 cm) were collected using a spatially explicit design within two 15-m 2 plots under three fertilization treatments in SG and GG croplands. A total of 288 samples were collected based on 12 plots and 24 samples in each plot. The fertilization treatments were no N input (NN), low N input (LN: 84 kg N ha -1 in urea) and high N input (HN: 168 kg N ha -1 in urea). The SOC, TN, SOC/TN (C: N), δ 13 C and δ 15 N were quantified and their within-plot variations and spatial distributions were achieved via descriptive and geostatistical methods. Results showed that SG generally displayed 10~120% higher plot-level variations in all variables than GG, and the plot-level variations were 20~77% higher in NN plots than LN and HN plots in SG but they were comparable in unfertilized and fertilized plots in GG. Relative to NN, LN and HN showed more significant surface trends and spatial structures in SOC and TN in both croplands, and the fertilization effect appeared more pronounced in SG. Spatial patterns in C: N, δ 13 C and δ 15 N were comparable among different fertilization treatments in both croplands. The descending within-plot variations were also identified among variables (SOC > TN > δ 15 N > C: N > δ 13 C). This study demonstrated that N fertilizations generally reduced the plot-level variance and simultaneously re-established spatial structures of SOC and TN in bioenergy croplands, which little varied with fertilization rate but was more responsive in switchgrass cropland.

59 BASIC BIOLOGICAL SCIENCES↗

Spatial variation of the low-temperature micro-photoluminescence of THM-grown CdZnTeSe

The quaternary compound CdZnTeSe (CZTS) has emerged as a next-generation detector material. Unlike CdZnTe (CZT), CZTS has the distinct advantages of little or no sub-grain boundary networks and a much lower concentration of Te inclusions, plus better compositional homogeneity. Thus, the material is expected to offer better spatial charge-transport homogeneity compared to CZT. However, an inhomogeneous distribution of point defects in both CZT and CZTS, such as Cd vacancies and residual impurities, can impose additional spatial inhomogeneity of the resistivity, carrier mobilities, and carrier lifetimes of the material. To investigate the spatial distribution of such defects in CZTS, low-temperature photoluminescence (PL) spectroscopic studies were performed at different positions along a single crystalline CZTS sample surface, which was grown by the traveling heater method (THM). In conclusion, the intensity variation of the PL emission of excitons bound to a neutral acceptor defect (A 0 , X), which is likely dominated by an acceptor like Cu-related level, with respect to the PL emission line from a neutral donor-bound exciton (D 0 , X) exhibited significant spatial variation, while the peak energy positions were approximately uniform due to the high compositional and bandgap homogeneity of the material.

II-VI semiconductor↗

Parameter Calibration for Johnson Cook and Preston-Tonks-Wallace Material Strength Models with Uncertainty Quantification

In this study we perform a Bayesian calibration of the parameters in the Johnson Cook (JC) material strength model, with uncertainty, using experiments with for a range of low and medium strain rates. For the parameter calibration, we used a variational Bayesian approach with experimental data from Hopkinson bar and quasi-static tests done on Oxygen Free High Conductivity (OFHC) copper. The estimated parameter values matched well with experimental data for a modest range of strain rates and temperatures. Through this method, we also recovered uncertainty and correlation information for the estimated parameter values. We also compared our results to a calibration of parameters in the Preston-Tonks-Wallace (PTW) material strength model, using the same variational Bayesian method. The parameters estimated for both models provided good agreement with the experimental data.

36 MATERIALS SCIENCE↗

Evaluating the Methodologies of Assessing Long-Term Variability of Wind Speed

Wind resources vary from year to year, and hence the long-term variability of wind speed is a key component of the overall uncertainty in the wind resource assessment process. Because wind-turbine power generation is highly sensitive to wind speed, the variability of wind resources has important implications for renewable power production. As many approaches of calculating the spread of distribution exist, we aim to identify the methods most appropriate for all types of wind-speed distributions. We review and compare various methods, such as the Coefficient of Variation (CoV) and the Robust Coefficient of Variation (RCoV). Because the length of data available also affects the resultant spread, we determine the minimum length of wind-speed records to quantify variability accurately. We use the net monthly energy production of 607 wind farms in the United States from the Energy Information Administration and the wind speeds of the MERRA-2 reanalysis dataset from NASA. After filtering the energy data via linear regression, we assess the accuracy of the variability calculations by comparing the correlations between the wind-speed variability estimates and the variabilities of wind-farm power productions. Results from statistically robust and resistant methods yield higher correlations between wind-speed and energy-production variabilities, as robust methods do not assume Gaussian wind-speed distributions and resistant methods are insensitive to wind speed extremes. We recommend RCoV to systematically estimate long-term variability, with 10 years of monthly mean data to achieve 90% statistical confidence. We will also provide generalized guidance on how to estimate long-term variability including case studies where we identify regions with high wind speeds and low wind-speed RCoVs, such as the Great Plains and the Upper Midwest in central United States.

17 WIND ENERGY↗