Search NASA⌕ Search

SEARCH · Search NASA

Results for “randomized 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,189 records · Page 66

Rigid body mode pointing accuracy and stability criteria for an orbiting spacecraft

The effect of gravity gradient torque disturbances (which is an order of magnitude greater than the other torque disturbances and is sinusoidal with a frequency of two cycles per orbit) on pointing accuracy can be predicted very precisely. However, since the sensor noise such as star tracker and gyro inherent noise, is random, their effect can only be predicted statistically. The most common method used is computer simulation with appropriate sensor model and adjusting spacecraft control laws for minimum pointing error, as presented by Harris (1972). The general investigation outlined theoretically relates spacecraft controls to pointing accuracy in terms of the sensor noise power spectrum. The results are presented in graphical form.

Sandhu, G. S.↗

Propagation of transients in a random medium

The propagation of transient scalar waves in a three-dimensional random medium is considered. The analysis is based on the smoothing method. An integro-differential equation for the coherent (or average) wave is derived and solved for the case of a statistically homogeneous and isotropic medium and a delta-function source. This yields the coherent Green's function of the medium. It is found that the waveform of the coherent wave depends generally on the distance from the source measured in terms of a certain dimensionless parameter. Based on the magnitude of this parameter, three propagation zones, called the near zone, the far zone, and the intermediate zone, are defined.

Wenzel, A. R.↗

U.S. National Congress of Applied Mechanics, 7th, University of Colorado, Boulder, Colo., June 3-7, 1974, Proceedings

A series of experimental, theoretical, and review papers dealing with various aspects of applied mechanics, including structural analysis, mathematical methods, and fluid dynamics. Some of the topics covered are: the direct formulation of two-dimensional mechanics problems, the mechanical aspects of athletics, Laplace's tidal equations, cochlear mechanics, dynamic constitutive equations, optical techniques for the determination of displacement and strain fields, chemical and biochemical problems involving nonlinear diffusion processes, thermal stresses, the rotating rod viscometer, finite-element methods in plasticity, ships of minimum resistance, and wide-band random vibration of structures. Individual items are announced in this issue.

Source record↗

Compositional inhomogeneities in a single Icelandic tholeiite flow

New trace element analyses of 25 basalt samples from a vertical traverse in an 11 m thick flow of Icelandic tholeiite show compositional differences which exceed analytical uncertainties and are random with respect to height. These variations suggest a short-range segregation model in which the amounts of phenocrysts, groundmass materials, and residual liquid vary randomly in different samples of the flow; a least-squares method is applied to determine if these variations reflect different proportions of crystallizing phases and residual liquid. Monte Carlo computations indicate that short-range segregation results in inhomogeneity of a magnitude of interflow compositional differences; this causes incoherence in elements with different geochemical behaviors while preserving coherence in elements of similar behavior.

Lindstrom, M. M.↗

Simulation and study of small numbers of random events

Random events were simulated by computer and subjected to various statistical methods to extract important parameters. Various forms of curve fitting were explored, such as least squares, least distance from a line, maximum likelihood. Problems considered were dead time, exponential decay, and spectrum extraction from cosmic ray data using binned data and data from individual events. Computer programs, mostly of an iterative nature, were developed to do these simulations and extractions and are partially listed as appendices. The mathematical basis for the compuer programs is given.

Shelton, R. D.↗

Linear and nonlinear dynamic instability of rotating polytropes

A three-dimensional hydrodynamic computer program is used to study the growth of nonaxisymmetric structures in rapidly rotating, self-gravitating polytropes. Models with polytropic index n = 0.8, 1.0, 1.3, 1.5, and 1.8 are studied. The initially axisymmetric equilibria are constructed by the Ostriker-Mark self-consistent-field method. The nonaxisymmetric pattern that develops out of low-amplitude random noise is a two-armed spiral with a well-defined pattern speed and growth rate which closely match properties of the toroidal mode predicted from the linear, second-order tensor-virial equation. A Fourier analysis of each polytrope's azimuthal density distribution shows that, even in the linear amplitude regime, higher-order angular patterns also develop exponentially in time. The higher-order patterns ultimately move in synchronization with the broad two-armed spiral, creating a narrow two-armed spiral. As the polytropic index is decreased, a more open and centrally more barlike pattern develops.

Williams, Harold A.↗

Kuhn-Tucker optimization based reliability analysis for probabilistic finite elements

The fusion of probability finite element method (PFEM) and reliability analysis for fracture mechanics is considered. Reliability analysis with specific application to fracture mechanics is presented, and computational procedures are discussed. Explicit expressions for the optimization procedure with regard to fracture mechanics are given. The results show the PFEM is a very powerful tool in determining the second-moment statistics. The method can determine the probability of failure or fracture subject to randomness in load, material properties and crack length, orientation, and location.

Liu, W. K.↗

Typical motions in multiple systems

In very old times, people counted - one, two, many. The author wants to show that they were right. Consider the motions of isolated bodies: (1) N = 1 - simple motion; (2) N = 2 - Keplerian orbits; and (3) N = 3 - this is the difficult problem. In general, this problem can be studied only by computer simulations. The author studied this problem over many years (see, e.g., Agekian and Anosova, 1967; Anosova, 1986, 1989 a,b). The principal result is that two basic types of dynamics take place in triple systems. The first special type is the stable hierarchical systems with two almost Keplerian orbits. The second general type is the unstable triple systems with complicated motions of the bodies. By random choice of the initial conditions, by the Monte-Carlo method, the stable systems comprised about approx. 10% of the examined cases; the unstable systems comprised the other approx. 90% of cases under consideration. In N greater than 3, the studies of dynamics of such systems by computer simulations show that we have in general also the motions roughly as at the cases 1 - 3 with the relative negative or positive energies of the bodies. In the author's picture, the typical trajectories of the bodies in unstable triple systems of the general type of dynamics are seen. Such systems are disrupted always after close triple approaches of the bodies. These approaches play a role like the gravitational slingshot. Often, the velocities of escapers are very large. On the other hand, the movie also shows the dynamical processes of a formation, dynamical evolution and disruption of the temporary wide binaries in triples and a formation of final hard massive binaries in the final evolution of triples.

Anosova, Joanna P.↗

Software testability and its application to avionic software

Randomly generated black-box testing is an established yet controversial method of estimating software reliability. Unfortunately, as software applications have required higher reliabilities, practical difficulties with black-box testing have become increasingly problematic. These practical problems are particularly acute in life-critical avionics software, where requirements of 10 exp -7 failures per hour of system reliability can translate into a probability of failure (POF) of perhaps 10 exp -9 or less for each individual execution of the software. This paper describes the application of one type of testability analysis called 'sensitivity analysis' to B-737 avionics software; one application of sensitivity analysis is to quantify whether software testing is capable of detecting faults in a particular program and thus whether we can be confident that a tested program is not hiding faults. We so 80 by finding the testabilities of the individual statements of the program, and then use those statement testabilities to find the testabilities of the functions and modules. For the B-737 system we analyzed, we were able to isolate those functions that are more prone to hide errors during system/reliability testing.

Voas, Jeffrey M.↗

Forward Monte Carlo Computations of Polarized Microwave Radiation

Microwave radiative transfer computations continue to acquire greater importance as the emphasis in remote sensing shifts towards the understanding of microphysical properties of clouds and with these to better understand the non linear relation between rainfall rates and satellite-observed radiance. A first step toward realistic radiative simulations has been the introduction of techniques capable of treating 3-dimensional geometry being generated by ever more sophisticated cloud resolving models. To date, a series of numerical codes have been developed to treat spherical and randomly oriented axisymmetric particles. Backward and backward-forward Monte Carlo methods are, indeed, efficient in this field. These methods, however, cannot deal properly with oriented particles, which seem to play an important role in polarization signatures over stratiform precipitation. Moreover, beyond the polarization channel, the next generation of fully polarimetric radiometers challenges us to better understand the behavior of the last two Stokes parameters as well. In order to solve the vector radiative transfer equation, one-dimensional numerical models have been developed, These codes, unfortunately, consider the atmosphere as horizontally homogeneous with horizontally infinite plane parallel layers. The next development step for microwave radiative transfer codes must be fully polarized 3-D methods. Recently a 3-D polarized radiative transfer model based on the discrete ordinate method was presented. A forward MC code was developed that treats oriented nonspherical hydrometeors, but only for plane-parallel situations.

Battaglia, A.↗

Cloud-Computing and Machine Learning in Support of Country-Level Land Cover and Ecosystem Extent Mapping in Liberia and Gabon

Liberia and Gabon joined the Gaborone Declaration for Sustainability in Africa (GDSA), established in 2012, with the goal of incorporating the value of nature intonational decision making by estimating the multiple services obtained from ecosystems using the natural capital accounting framework. In this study, we produced 30-m resolution 10 classes land cover maps for the 2015 epoch for Liberia and Gabon using the Google Earth Engine (GEE) cloud platform to support the ongoing natural capital accounting efforts in these nations. We pro-pose an integrated method of pixel-based classification using Landsat 8 data, the Random Forest(RF) classifier and ancillary data to produce high quality land cover products to fit abroad range of applications, including natural capital accounting. Our approach focuses on a pre-classification filtering (Masking Phase) based on spectral signature and ancillary data to reduce the number of pixels prone to be misclassified; therefore, increasing the quality of the final product. The proposed approach yields an overall accuracy of 83% and 81% for Liberia and Gabon, respectively, out performing prior land cover products for these countries in both thematic content and accuracy. Our approach, while relatively simple and highly replicable, was able to produce high quality land cover products to fill an observational gap in up to date land cover data at national scale for Liberia and Gabon.

Celio de Sousa↗

Use of Uncertainty in Modeling DFLAW Treatment at Hanford - 17115

Average or mean is understood by most; standard deviation (SD) is understood by few. Certainty is easy to model; uncertainty is not. The flowsheet simulation, described in this paper, placed uncertainties in the G2 model of the Hanford Tank Waste Treatment and Immobilization Plant's (WTP) Direct Feed Low-Activity Waste (DFLAW) process for performance and sensitivity evaluations. G2 is a dynamic material balance model that starts and proceeds with operating logic until all the feed is processed. Dynamic material balance is achieved by tracking the flows, storage, and change of all materials within the plant as time increments. Processing equipment and limitations of the flowsheet are modeled (e.g., equipment types, piping, volumes, flowrates, efficiencies, and physical and chemical environments that impact separations, phase equilibriums, and chemical reactions). Operating logic represents the rules and strategies of a running plant. The G2 model simulated both certainty and uncertainty scenarios for comparison purposes. For the uncertainty scenarios, process decisions were made using presumed values to simulate various operational measurement errors and process variabilities. Presumed means that the known G2 values are adjusted to replicate measurement errors and variabilities. These include: 1) volume measurement, 2) sample analysis, 3) melter decontamination factor (DF), 4) glass former weighing, and 5) glass former composition variations. The errors were randomly generated in accordance with SD values of measurement devices, laboratory methods, and glass former composition ranges. The G2 model runs give insight on the DFLAW process' controllability and the expected range of glass properties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Hypergraph Random Walks, Laplacians, and Clustering

We propose a flexible framework for clustering hypergraph-structured data based on recently proposed random walks utilizing edge-dependent vertex weights. When incorporating edge-dependent vertex weights (EDVW), a weight is associated with each vertex-hyperedge pair, yielding a weighted incidence matrix of the hypergraph. Such weightings have been utilized in term-document representations of text data sets. We explain how random walks with EDVW serve to construct different hypergraph Laplacian matrices, and then develop a suite of clustering methods that use these incidence matrices and Laplacians for hypergraph clustering. Using 20Newsgroup, U.S. patent, Reuters' Corpus Volume 1, and genetics data sets, we compare the performance of these clustering algorithms experimentally against a variety of existing hypergraph clustering methods. We show that the proposed methods produce higher-quality clusters.

hypergraphs, clustering, laplacian, random walk, M↗

Identification and Preservation of Arm Lymphatic System in Axillary Dissection for Breast Cancer to Reduce Arm Lymphedema Events: A Randomized Clinical Trial

Background: Controversy in axillary reverse mapping in axillary lymph node dissection (ALND) possibly results from incomplete recognition of the arm lymphatic system (ALS) and its compromise to oncological safety. The iDEntification and Preservation of ARm lymphaTic system (DEPART) technique facilitates complete identification of ALS; therefore, its use may decrease the occurrence of arm lymphedema. This study aimed to examine the arm lymphedema rate, locoregional recurrence, and feasibility to perform DEPART in ALND. Methods: Patients from February 2013 to October 2017 from two tertiary referral centers were randomly assigned to two groups. In the study group, indocyanine green and methylene blue (MB) were utilized to identify arm sentinel nodes, and 0.1 ml MB was injected into the arm sentinel nodes to reveal the subsequent-echelon nodes and lymphatics. Gross arm lymph nodes were examined by intraoperative partial frozen section and were removed if positive. Arm lymphedema, local recurrence, regional recurrence, and distant metastasis were recorded at different follow-up examinations. Results: Arm sentinel nodes were identified in 573 (83.2%) patients. Subsequent-echelon nodes and lymphatics were visualized in 558 (97.4%) patients. Metastatic arm nodes were identified in 38 (6.8%) patients. The arm lymphedema rate was 3.3% (18/543) in the study group versus 15.3% (99/648) in the control group (p < 0.001) after 37-month median follow-up. Regional recurrence showed no difference between the two groups (1.4% and 1.2%, respectively) (p = 0.392). Conclusions: DEPART can benefit breast cancer patients who undergo ALND, reducing the arm lymphedema rate without adversely affecting the morbidity of regional recurrence.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Prediction of atomization energy using graph kernel and active learning

Data-driven prediction of molecular properties presents unique challenges to the design of machine learning methods concerning data structure/dimensionality, symmetry adaption, and confidence management. In this paper, we present a kernel-based pipeline that can learn and predict the atomization energy of molecules with high accuracy. The framework employs Gaussian process regression to perform predictions based on the similarity between molecules, which is computed using the marginalized graph kernel. To apply the marginalized graph kernel, a spatial adjacency rule is first employed to convert molecules into graphs whose vertices and edges are labeled by elements and interatomic distances, respectively. We then derive formulas for the efficient evaluation of the kernel. Specific functional components for the marginalized graph kernel are proposed, while the effects of the associated hyperparameters on accuracy and predictive confidence are examined. We show that the graph kernel is particularly suitable for predicting extensive properties because its convolutional structure coincides with that of the covariance formula between sums of random variables. Using an active learning procedure, we demonstrate that the proposed method can achieve a mean absolute error of 0.62 ± 0.01 kcal/mol using as few as 2000 training samples on the QM7 dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Active Suppression Of Vibrations On Aircraft Structures

Method of active suppression of nonlinear and nonstationary vibrations developed to reduce sonic fatigue and interior noise in high-speed aircraft. Structure of aircraft exhibits periodic, chaotic, and random vibrations when forced by high-intensity sound from jet engines, shock waves, turbulence, and separated flows. Method of suppressing vibrations involves feedback control: Strain gauges or other sensors mounted in paths of propagation of vibrations on structure sense vibrations; outputs of sensors processed into control signal applied to actuator mounted on structure, inducing compensatory forces.

Maestrello, Lucio↗

Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.

58 GEOSCIENCES↗

Formation and structural phase transition in Co atomic chains on a Cu(775) surface

The formation of Co atomic chains on a Cu(775) surface is investigated by the kinetic Monte Carlo method. It is found that the length of Co atomic chains formed as a result of self-organization during epitaxial growth is a random quantity and its mean value depends on the parameters of the experiment. The existence of two structural phases in atomic chains is detected using the density functional theory. In the first phase, the separations between an atom and its two nearest neighbors in a chain are 0.230 and 0.280 nm. In the second phase, an atomic chain has identical atomic spacings of 0.255 nm. It is shown that the temperature of the structural phase transition depends on the length of the atomic chain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗