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 199 records · Page 11

Computation of Variability in the Average Thermal and Mechanical Properties of a Melt-infiltrated Sic/sic Composite

Thermal conductivity and tensile properties (elastic modulus, proportional limit strength, inplane tensile strength, and strain to failure) of a SiC/SiC composite were experimentally determined at 816 and 1204 C. Tests were performed at room temperature on the same material to obtain interlaminar shear and tensile strength. For each thermal and mechanical property 24 tests were conducted to capture the variation. A random sampling method was used to quantify the variability exhibited by the mean values of the thermal and mechanical properties. The minimum number of tests required to characterize the mean value of each thermal or mechanical property for the composite was determined by varying the group size of the random sample.

Kalluri, Sreeramesh↗

Proof-of-Concept Study for Uncertainty Quantification and Sensitivity Analysis using the BRL Shaped-Charge Example

These are the slides for a graduate presentation at Mississippi State University. It covers the following: the BRL Shaped-Charge Geometry in PAGOSA, mesh refinement study, surrogate modeling using a radial basis function network (RBFN), ruling out parameters using sensitivity analysis (equation of state study), uncertainty quantification (UQ) methodology, and sensitivity analysis (SA) methodology. In summary, a mesh convergence study was used to ensure that solutions were numerically stable by comparing PDV data between simulations. A Design of Experiments (DOE) method was used to reduce the simulation space to study the effects of the Jones-Wilkins-Lee (JWL) Parameters for the Composition B main charge. Uncertainty was quantified by computing the 95% data range about the median of simulation output using a brute force Monte Carlo (MC) random sampling method. Parameter sensitivities were quantified using the Fourier Amplitude Sensitivity Test (FAST) spectral analysis method where it was determined that detonation velocity, initial density, C1, and B1 controlled jet tip velocity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The First Intervention Study in Elder Self-Neglect: A Randomized Clinical Trial to Improve Vitamin D Levels

Despite high mortality rates, elder self-neglect is characterized by refusal of medical and social interventions. To date there have been no tested clinical interventions in elders who self-neglect. Previous research from the TEAM Institute has shown significantly low vitamin D levels in this population. This study aimed to determine the feasibility of a clinical intervention. Replacement of vitamin D was chosen because of its ease of administration and favorable safety profile. Methods: A randomized clinical trial using directly observed therapy of vitamin D was conducted using 50 elders, >65 years of age, with Adult Protective Services (APS) validated self-neglect. A staggered intervention with waiting controls was used to maximize statistical power. One-third (n=17) of the group was administered 50,000 IU vitamin D2 (ergocalciferol) monthly and the remainder (n=33) were administered 400 IU monthly. Serum 25-OH vitamin D was assessed at baseline and 5-months. Results: 69% agreed to participate in the study and of those n=40 (80%) remained at 5-months. At baseline, 12% (n=7) were deficient in vitamin D (<30nmol/L) and approximately 38% (n=22) had inadequate vitamin D levels (<50nmol/L). The baseline 25-OH vitamin D level was 59 nmol/L +25 (mean SD), and increased significantly to 72nmol/L +21 nmol/L at 5-months. Conclusion: These data are the first to provide evidence that clinical interventions are feasible in elders who self-neglect. The increase in vitamin D levels confirmed that the study personnel were able to successfully intervene community-dwelling elders with self-neglect. This study sets the precedent for future intervention and prevention studies

Burnett, Jason↗

Uncertainty Quantification of Activation Due to Cross Section Data in Neutron Shielding Calculation

Quantification of activation in structure materials around the nuclear reactor through calculation is necessary before start of decommissioning. Then a neutron shielding calculation is carried out to estimate activation of structural materials. Neutron shielding calculation is performed out for spatially large regions around a reactor through deep penetration calculation, calculation result would have significant uncertainty. In order to minimize total amount of structure materials treated as radioactive waste, quantitative estimation on uncertainty of structural activation and its reduction is important. The uncertainties due to calculation arises from that of input data and that of calculation method. Examples of uncertainty of input data is those for nuclear reaction cross section and atom density of the nuclide in the structure. In this study, we focus on the uncertainty of cross section, which would be one of the dominant uncertainties, and try to estimate uncertainty of activation due to cross section data in a neutron shielding calculation. The random sampling method (RS) is used as the uncertainty quantification method. With the RS method, uncertainty estimation can be carried out without modification of computer codes for a neutron shielding calculation. (authors)

36 MATERIALS SCIENCE↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗

Cleanroom software development

The 'cleanroom' software development process is a technical and organizational approach to developing software with certifiable reliability. Key ideas behind the process are well structured software specifications, randomized testing methods and the introduction of statistical controls; but the main point is to deny entry for defects during the development of software. This latter point suggests the use of the term 'cleanroom' in analogy to the defect prevention controls used in the manufacturing of high technology hardware. In the 'cleanroom', the entire software development process is embedded within a formal statistical design, in contrast to executing selected tests and appealing to the randomness of operational settings for drawing statistical inferences. Instead, random testing is introduced as a part of the statistical design itself so that when development and testing are completed, statistical inferences are made about the operation of the system.

Dyer, M.↗

A software‐defined networks ‐based measurement method of network traffic for 6G technologies

Abstract In software‐defined networks (SDN) for 6G technologies, the controller should accurately and efficiently measure the flow traffic in switches for traffic engineering. Fine‐grained flow measurement can more accurately describe the traffic in the network, but it also consumes much more resources. To reduce the overhead incurred in the measurement process and obtain the approximate fine‐grained measurements, we propose a novel lightweight measurement scheme that runs in the controller. The novel lightweight measurement architecture consists of two parts: coarse‐grained measurement and interpolation‐optimization. In the first part, based on the SDN architecture, we use the pull‐based random sampling method to quickly obtain the coarse‐grained measurement of flow traffic through OpenFlow protocol. In the second part, we insert some discrete values into the coarse‐grained measurement with the interpolation theory, then we optimize interpolation results until finding the optimal fine‐grained flow traffic measurement by utilizing the multiconstraint method. We verified the feasibility of the proposed measurement method, and simulation results show that the measurement error of the proposed method is under 25%, but the flow measurement overhead of this method occupies only 3.3% compared with that of the fine‐grained method.

Huo, Liuwei↗

A New Aerodynamic Data Dispersion Method for Launch Vehicle Design

A novel method for implementing aerodynamic data dispersion analysis is herein introduced. A general mathematical approach combined with physical modeling tailored to the aerodynamic quantity of interest enables the generation of more realistically relevant dispersed data and, in turn, more reasonable flight simulation results. The method simultaneously allows for the aerodynamic quantities and their derivatives to be dispersed given a set of non-arbitrary constraints, which stresses the controls model in more ways than with the traditional bias up or down of the nominal data within the uncertainty bounds. The adoption and implementation of this new method within the NASA Ares I Crew Launch Vehicle Project has resulted in significant increases in predicted roll control authority, and lowered the induced risks for flight test operations. One direct impact on launch vehicles is a reduced size for auxiliary control systems, and the possibility of an increased payload. This technique has the potential of being applied to problems in multiple areas where nominal data together with uncertainties are used to produce simulations using Monte Carlo type random sampling methods. It is recommended that a tailored physics-based dispersion model be delivered with any aerodynamic product that includes nominal data and uncertainties, in order to make flight simulations more realistic and allow for leaner spacecraft designs.

Pinier, Jeremy T.↗

Emergence of many-body quantum chaos via spontaneous breaking of unitarity

It is suggested that many-body quantum chaos appears as the spontaneous symmetry breaking of unitarity in interacting quantum many-body systems. It has been shown that many-body level statistics, probed by the spectral form factor (SFF) defined as K(η,t)=$\langle$|Tr exp(-ηH+itH)| 2 $\rangle$, is dominated by a diffuson-type mode in a field theory analysis. The key finding of this Letter is that the “unitary” η = 0 case is different from the η → 0 ± limit, with the latter leading to a finite mass of these modes due to interactions. This mass suppresses a rapid exponential ramp in the SFF, which is responsible for the fast emergence of Poisson statistics in the noninteracting case, and gives rise to a nontrivial random matrix structure of many-body levels. The interaction-induced mass in the SFF shares similarities with the dephasing rate in the theory of weak localization and the Lyapunov exponent of the out-of-time-ordered correlators.

36 MATERIALS SCIENCE↗

Unmixing hyperspectral images using Markov random fields

This paper proposes a new spectral unmixing strategy based on the normal compositional model that exploits the spatial correlations between the image pixels. The pure materials (referred to as endmembers) contained in the image are assumed to be available (they can be obtained by using an appropriate endmember extraction algorithm), while the corresponding fractions (referred to as abundances) are estimated by the proposed algorithm. Due to physical constraints, the abundances have to satisfy positivity and sum-to-one constraints. The image is divided into homogeneous distinct regions having the same statistical properties for the abundance coefficients. The spatial dependencies within each class are modeled thanks to Potts-Markov random fields. Within a Bayesian framework, prior distributions for the abundances and the associated hyperparameters are introduced. A reparametrization of the abundance coefficients is proposed to handle the physical constraints (positivity and sum-to-one) inherent to hyperspectral imagery. The parameters (abundances), hyperparameters (abundance mean and variance for each class) and the classification map indicating the classes of all pixels in the image are inferred from the resulting joint posterior distribution. To overcome the complexity of the joint posterior distribution, Markov chain Monte Carlo methods are used to generate samples asymptotically distributed according to the joint posterior of interest. Simulations conducted on synthetic and real data are presented to illustrate the performance of the proposed algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying climate–growth relationships at the stand level in a mature mixed‐species conifer forest

Abstract A range of environmental factors regulate tree growth; however, climate is generally thought to most strongly influence year‐to‐year variability in growth. Numerous dendrochronological (tree‐ring) studies have identified climate factors that influence year‐to‐year variability in growth for given tree species and location. However, traditional dendrochronology methods have limitations that prevent them from adequately assessing stand‐level (as opposed to species‐level) growth. We argue that stand‐level growth analyses provide a more meaningful assessment of forest response to climate fluctuations, as well as the management options that may be employed to sustain forest productivity. Working in a mature, mixed‐species stand at the Howland Research Forest of central Maine, USA , we used two alternatives to traditional dendrochronological analyses by (1) selecting trees for coring using a stratified (by size and species), random sampling method that ensures a representative sample of the stand, and (2) converting ring widths to biomass increments, which once summed, produced a representation of stand‐level growth, while maintaining species identities or canopy position if needed. We then tested the relative influence of seasonal climate variables on year‐to‐year variability in the biomass increment using generalized least squares regression, while accounting for temporal autocorrelation. Our results indicate that stand‐level growth responded most strongly to previous summer and current spring climate variables, resulting from a combination of individualistic climate responses occurring at the species‐ and canopy‐position level. Our climate models were better fit to stand‐level biomass increment than to species‐level or canopy‐position summaries. The relative growth responses (i.e., percent change) predicted from the most influential climate variables indicate stand‐level growth varies less from to year‐to‐year than species‐level or canopy‐position growth responses. By assessing stand‐level growth response to climate, we provide an alternative perspective on climate–growth relationships of forests, improving our understanding of forest growth dynamics under a fluctuating climate.

Teets, Aaron↗

Monte Carlo Study of the Effect of Initial States and Structural Defects on Nonequilibrium Critical Behavior of the Three-Dimensional Ising Model

The effect of various initial magnetizations m{sub 0} and structural defects the nonequilibrium critical behavior of the three-dimensional Ising model is numerically studied. Based on an analysis of the time dependence of the magnetization and the two-time dependence of the autocorrelation function and dynamic susceptibility, the significant effect of initial states on relaxation magnetizations and aging effects characterized by anomalous relaxation inhibition and correlation in the system with increasing waiting time was revealed. The fluctuation–dissipation theorem violation was studied, and the values of the limit fluctuation–dissipation ratio (FDR) are calculated. It is shown that two universality subclasses can be distinguished in the nonequilibrium critical behavior of the three-dimensional Ising model with random initial magnetization m{sub 0} These subclasses correspond to the system evolution from the high-temperature (m{sub 0} = 0) and low-temperature (m{sub 0} = 1) initial states with limit FDRs characteristic of these states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Extended Main-sequence Turn-offs in Intermediate-age Star Clusters: Stellar Rotation Diminishes, but Does Not Eliminate, Age Spreads

Extended main-sequence turn-off (eMSTO) regions are a common feature in color–magnitude diagrams of young- and intermediate-age star clusters in the Magellanic Clouds. The nature of eMSTOs remains debated in the literature. The currently most popular scenarios are extended star formation activity and ranges of stellar rotation rates. Here we study details of differences in main-sequence turn-off (MSTO) morphology expected from spreads in age versus spreads in rotation rates, using Monte Carlo simulations with the Geneva syclist isochrone models that include the effects of stellar rotation. We confirm a recent finding of Niederhofer et al. that a distribution of stellar rotation velocities yields an MSTO extent that is proportional to the cluster age, as observed. However, we find that stellar rotation yields MSTO crosscut widths that are generally smaller than observed ones at a given age. We compare the simulations with high-quality Hubble Space Telescope data of NGC 1987 and NGC 2249, which are the two only relatively massive star clusters with an age of ∼1 Gyr for which such data is available. We find that the distribution of stars across the eMSTOs of these clusters cannot be explained solely by a distribution of stellar rotation velocities, unless the orientations of rapidly rotating stars are heavily biased toward an equator-on configuration. Under the assumption of random viewing angles, stellar rotation can account for ∼60% and ∼40% of the observed FWHM widths of the eMSTOs of NGC 1987 and NGC 2249, respectively. In contrast, a combination of distributions of stellar rotation velocities and stellar ages fits the observed eMSTO morphologies very well.

79 ASTRONOMY AND ASTROPHYSICS↗

Signatures of multiple charge excitations in resonant inelastic x-ray scattering spectra of metals

In this work, we study how multiple charge excitations appear in the resonant inelastic x-ray scattering (RIXS) spectra of metals. Typically, the focus is upon single excitations in the problem, the plasmons and electron-hole pairs, while the multiexcitation processes are neglected. However, at small momentum transfer the multiexcitation contributions may dominate the signal and one needs to understand how to interpret the data in such cases. In particular, we demonstrate how to “decode” the total multiexcitation intensity and extract the plasmon dispersion. Furthermore, while our calculations are based on the random phase approximation for jellium, which does not allow us to obtain quantitatively precise results in the entire region of parameters, we expect them to capture semiqualitatively all features expected for charged Fermi-liquid states, including universal and singular properties of the RIXS spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Randomized algorithms for generalized singular value decomposition with application to sensitivity analysis

The generalized singular value decomposition (GSVD) is a valuable tool that has many applications in computational science. However, computing the GSVD for large-scale problems is challenging. Motivated by applications in hyper-differential sensitivity analysis (HDSA), in this work we propose new randomized algorithms for computing the GSVD which use randomized subspace iteration and weighted QR factorization. Detailed error analysis is given which provides insight into the accuracy of the algorithms and the choice of the algorithmic parameters. We demonstrate the performance of our algorithms on test matrices and a large-scale model problem where HDSA is used to study subsurface flow.

97 MATHEMATICS AND COMPUTING↗

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations

In this work, we present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the butterfly algorithm and its hierarchical matrix extension for compressing and factorizing large frontal matrices via graph-distance guided entry evaluation or randomized matrix-vector multiplication-based schemes. Complexity analysis and numerical experiments demonstrate $\mathcal{O}(N\log^2 N)$ computation and $\mathcal{O}(N)$ memory complexity when applied to an $N\times N$ sparse system arising from 3D high-frequency Helmholtz and Maxwell problems.

97 MATHEMATICS AND COMPUTING↗