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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Probabilistic structural analysis computer code (NESSUS)

Probabilistic structural analysis has been developed to analyze the effects of fluctuating loads, variable material properties, and uncertain analytical models especially for high performance structures such as SSME turbopump blades. The computer code NESSUS (Numerical Evaluation of Stochastic Structure Under Stress) was developed to serve as a primary computation tool for the characterization of the probabilistic structural response due to the stochastic environments by statistical description. The code consists of three major modules NESSUS/PRE, NESSUS/FEM, and NESSUS/FPI. NESSUS/PRE is a preprocessor which decomposes the spatially correlated random variables into a set of uncorrelated random variables using a modal analysis method. NESSUS/FEM is a finite element module which provides structural sensitivities to all the random variables considered. NESSUS/FPI is Fast Probability Integration method by which a cumulative distribution function or a probability density function is calculated.

Shiao, Michael C.↗
DOE OSTI2020

Continued Discussion of Failure Mode Modeling and Overall Component Reliability: Is the Data Missing or Censored?

This paper is a continuation of a paper presented at the 13th Probabilistic Safety Assessment and Management Conference, in which a discussion of modeling failure modes of complex components was presented; see Paulos and Smith (2016). The paper demonstrated how the typical method of treating failure modes as being exponential in nature may yield optimistic predictions when estimating how improvements to components will perform in the real world. Instead of relying on textbook methods, it is more accurate to model the failure modes as a race in time; unfortunately, this does not give a closed-form solution. This work uses a simulation to solve for the model of the world, and is compared to the standard methodology of treating the failure modes as exponential random failures. The standard method is shown to have optimistic predictions, which will lead to prediction errors when failure modes were removed or “fixed.” The failure mode methodology presented in the first paper treated the data as being censored when the test stopped. In this paper, we will compare the results from treating the data as both censored and missing, and comparisons will be shown showing which method is more accurate and under what conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗
DOE OSTI2019

Two-Year Mortality After Angioplasty of the Femoro-Popliteal Artery with Uncoated Balloons and Paclitaxel-Coated Balloons—A Pooled Analysis of Four Randomized Controlled Multicenter Trials

Purpose: In view of a recent meta-analysis reporting increased mortality following angioplasty with paclitaxel-coated devices in peripheral arteries, we performed a patient-level 2-year mortality analysis based on pooled original data of four randomized controlled trials (THUNDER, FEMPAC, PACIFIER and CONSEQUENT). Methods and Results: Clinical data of four randomized controlled trial were pooled to assess 2-year mortality following paclitaxel-coated balloon (PCB) angioplasty compared to angioplasty without paclitaxel (control group). A logistic regression model was applied to identify potential predictors of mortality. At two years, 13 of 185 (7.0%) patients had died in the control group and 16/184 (8.7%) in the PBC group, p = 0.55. Kaplan–Meier analysis revealed no significant difference from all-cause death at 2 years (log rank p = 0.54). Causes of death were well balanced between the groups with no pattern or trend in favour of any specific causes in the PBC group. Logistic regression revealed that treatment groups (controls or PBC) were not a predictor of 2-year mortality. The only predictor for mortality was patient age ≥ 75 years. The delivered paclitaxel doses per patient were not significantly different in patients that died and those who did not die during the 24-month follow-up (5.300 ± 4.224 μg vs. 6.248 ± 4.629 μg, p = 0.433). Conclusions: Based on original patient-level data of four pooled randomized controlled trials, we found no increase in 2-year mortality in patients treated with PCB compared to control patients treated with uncoated balloons. Causes of death were well balanced between PCB and control patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2020

W$^+$W$^-$ boson pair production in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A measurement of the W+W- boson pair production cross section in proton-proton collisions at s=13 TeV is presented. The data used in this study are collected with the CMS detector at the CERN LHC and correspond to an integrated luminosity of 35.9 fb-1. The W+W- candidate events are selected by requiring two oppositely charged leptons (electrons or muons). Two methods for reducing background contributions are employed. In the first one, a sequence of requirements on kinematic quantities is applied allowing a measurement of the total production cross section, 117.6±6.8 pb, which agrees well with the theoretical prediction. Fiducial cross sections are also reported for events with zero or one jet, and the change in the zero-jet fiducial cross section with the jet transverse momentum threshold is measured. Normalized differential cross sections are reported within the fiducial region. A second method for suppressing background contributions employs two random forest classifiers. The analysis based on this method includes a measurement of the total production cross section and also a measurement of the normalized jet multiplicity distribution in W+W- events. Finally, a dilepton invariant mass distribution is used to probe for physics beyond the standard model in the context of an effective field theory, and constraints on the presence of dimension-6 operators are derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1975

A primer for structural response to random pressure fluctuations

A review was made of power spectral methods for determining linear response of structures to random pressure fluctuations. Various simplifying assumptions are made for the purpose of obtaining useful formula for structural response. The transmission of sound through a flexible structure into an interior cavity was also treated.

Dowell, E. H.↗
NASA NTRS1989

Evaluation of a multi-point method for determining acoustic impedance

A multipoint method for determining acoustic impedance was evaluated in comparison with the traditional standing wave and two-microphone methods using 30 test samples covering the reflection factor magnitude range 0.004-0.999. The multipoint method is shown to combine the strengths of the standing wave and two-microphone methods while avoiding some of their inherent weaknesses. In particular, the results obtained suggest that the multipoint method will be less subject to flow induced random error than the two-microphone method in the presence of significant broadband noise levels associated with mean flow.

Jones, Michael G.↗
NASA NTRS1990

Limits to derived flare properties using estimates for the background fluxes - Examples from GOES

A method has been presented for restricting the random derived flare temperatures and emission measures. The method is based on three tests and is applicable not only to disk-integrated observations, but also for any high-resolution observations that cannot distinguish the emissions from flaring and nonflaring plasma. The first two tests are based on the assumption that both the temperature and the emission measure should increase as the flare flux rises. The third assumes that the preflare temperature of the plasma that flares should exceed the temperature averaged over all other emitting sources. These tests were applied to two solar flares to illustrate the improvements in the derived flare properties. It is shown that flare properties do not always lie between those derived by the use or subtraction of the total preflare flux. The flare properties can, however, be bounded by subtracting and using the long- and short-wavelength preflare fluxes, respectively.

Bornmann, P. L.↗
DOE OSTI2020

Data-Driven Preemptive Voltage Monitoring and Control Using Probabilistic Voltage Sensitivities

Increased penetration levels of distributed variable renewable generation can cause random voltage fluctuations and violations at multiple nodes. Traditional methods of voltage control typically involve reactionary responses of capacitor banks, tap changers, and recently even smart inverters. But because of the lack of foresight in voltage violations, these controls are ineffective to completely mitigate the issue. Therefore, new methods of predicting voltage violations subject to random power injection changes in the distribution network are needed, which can be used to guide optimal and dynamic methods of voltage control. This work lays the foundation for such preemptive voltage monitoring and control by proposing an analytical and sensor data-driven voltage sensitivity analysis method. Driven by stochastic data and forecasts, the method can be used to develop probabilistic voltage sensitivities and consequently to predict system nodes with high likelihood of voltage limit violations. The effectiveness of this method is tested on IEEE 69-node distribution system integrated with distributed solar. The results demonstrate the proposed method's ability to successfully predict nodes with high probability of voltage violations for a specific time-series simulation. The results also demonstrate the ability to guide timely power injection control actions to mitigate future voltage violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Data-Driven Preemptive Voltage Monitoring and Control Using Probabilistic Voltage Sensitivities: Preprint

Increased penetration of distributed variable re-newable generation may cause random voltage fluctuations and violations at multiple nodes. Traditional methods of voltage control typically involve reactionary responses of capacitor banks, tap changers, and recently even smart inverters. But due to the lack of foresight in voltage violations, these controls are ineffective to completely mitigate the issue. Therefore, new methods of predicting voltage violations subject to random power injection changes in the distribution network is needed, which can be used to guide optimal and dynamic methods of voltage control. This work lays the foundation for such preemptive voltage monitoring and control by proposing an analytical and sensor data-driven voltage sensitivity analysis method. Driven by stochastic data and forecasts, the method can be used to develop probabilistic voltage sensitivities and consequently to predict system nodes with high likelihood of voltage limit violations. The effectiveness of this method is tested on IEEE 69 node distribution system integrated with distributed solar. The results demonstrate the proposed method's ability to successfully predict nodes with high probability of voltage violations for a specific time series simulation. The results also demonstrate the ability to guide timely power injection control actions to mitigate future voltage violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1972

Geometry of quantum signal detection.

Consideration of a binary quantum signal detection problem in a two-dimensional Hilbert space. The optimum detection problem is reduced to the problem of finding the locus of a vector which has a maximum projection along the fixed a priori probability vector of hypotheses. It is shown that the desired locus can be determined by a geometrical method involving the use of a randomized decision strategy. It is further shown that this geometrical method can be applied to the optimum solution of a binary detection problem described in a product Hilbert space.

Harger, R. O.↗
NASA NTRS1991

Estimating The SNR Of AVIRIS Data

Combination of established statistical and filtering techniques yields improved estimates of signal-to-noise ratio (SNR) of multispectral imagery produced by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Method yields more accurate estimate of random variability in AVIRIS imagery. Combination of techniques called "geostatistical" method. In first step of method, periodic noise removed by notch filtering in frequency domain. In subsequent steps, sensor noise and variability within picture elements isolated by use of semi-variance of picture elements.

Curran, Paul J.↗
NASA NTRS2004

Reduced Order Methods for Prediction of Thermal-Acoustic Fatigue

The goal of this investigation is to assess the quality of high-cycle-fatigue life estimation via a reduced order method, for structures undergoing random nonlinear vibrations in a presence of thermal loading. Modal reduction is performed with several different suites of basis functions. After numerically solving the reduced order system equations of motion, the physical displacement time history is obtained by an inverse transformation and stresses are recovered. Stress ranges obtained through the rainflow counting procedure are used in a linear damage accumulation method to yield fatigue estimates. Fatigue life estimates obtained using various basis functions in the reduced order method are compared with those obtained from numerical simulation in physical degrees-of-freedom.

Przekop, A.↗
DOE OSTI2018

Transportation of Low-Level Radioactive Waste from KAERI to the Final Repository - 18179

The Korea Atomic Energy Research Institute (KAERI) collects the radioactive wastes generated during the operation of nuclear facilities and R and D processes. Radioactive wastes are treated and repackaged according to the national radioactive waste management policy, and then transferred to the final disposal site at the Korea Radioactive Waste Agency (KORAD). After 2015, the wastes delivered to the final disposal site have consisted of spent filters and combustible wastes such as clothes and gloves generated during R and D or facility operation processes at KAERI. The wastes generated by each facility are collected at the radioactive waste treatment facility, then treated and repackaged in 200 L drums according to their generation characteristics, i.e. each drum contains waste of the same type from the same facility, and with the same date of generation. Then, a radionuclide inventory is obtained for the repackaged drum by using a direct measurement method such as the drum assay system, radiochemical analysis through random sampling or an indirect measurement method such as applying a scaling factor. When performing sampling of the radioactive waste drums, the waste samples with a high contamination are selected in a conservative approach as representative. The radioactive waste drums are inspected and checked for acceptability by the Quality Management Division of KAERI, and a preliminary inspection is carried out on-site by KORAD. At this time, the acceptability of the radioactive waste drums is assessed through visual inspection, the appropriateness of their nuclide analysis value is confirmed, and a survey or smear test of sample drums is performed. Drums that meet the acceptance criteria are transferred to KORAD after a radiation environmental impact assessment of the transportation. Since 2015, 1,084 low-level operational radioactive wastes have been transferred from KAERI to the final disposal site, and further preparation work is now being made for the transportation of additional drums. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Large scale Brownian dynamics of confined suspensions of rigid particles

We introduce methods for large-scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method [F. Balboa Usabiaga et al., Commun. Appl. Math. Comput. Sci. 11(2), 217–296 (2016)] at a cost comparable to the cost of deterministic simulations. We demonstrate that we can efficiently generate deterministic and random displacements for many particles using preconditioned Krylov iterative methods, if kernel methods to efficiently compute the action of the Rotne-Prager-Yamakawa (RPY) mobility matrix and its “square” root are available for the given boundary conditions. These kernel operations can be computed with near linear scaling for periodic domains using the positively split Ewald method. Here we study particles partially confined by gravity above a no-slip bottom wall using a graphical processing unit implementation of the mobility matrix-vector product, combined with a preconditioned Lanczos iteration for generating Brownian displacements. We address a major challenge in large-scale BD simulations, capturing the stochastic drift term that arises because of the configuration-dependent mobility. Unlike the widely used Fixman midpoint scheme, our methods utilize random finite differences and do not require the solution of resistance problems or the computation of the action of the inverse square root of the RPY mobility matrix. We construct two temporal schemes which are viable for large-scale simulations, an Euler-Maruyama traction scheme and a trapezoidal slip scheme, which minimize the number of mobility problems to be solved per time step while capturing the required stochastic drift terms. We validate and compare these schemes numerically by modeling suspensions of boomerang-shaped particles sedimented near a bottom wall. Using the trapezoidal scheme, we investigate the steady-state active motion in dense suspensions of confined microrollers, whose height above the wall is set by a combination of thermal noise and active flows. We find the existence of two populations of active particles, slower ones closer to the bottom and faster ones above them, and demonstrate that our method provides quantitative accuracy even with relatively coarse resolutions of the particle geometry.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Numerical solution of random singular integral equation appearing in crack problems

The solution of several elasticity problems, and particularly crack problems, can be reduced to the solution of one-dimensional singular integral equations with a Cauchy-type kernel or to a system of uncoupled singular integral equations. Here a method for the numerical solution of random singular integral equations of Cauchy type is presented. The solution technique involves a Chebyshev series approximation, the coefficients of which are the solutions of a system of random linear equations. This method is applied to the problem of periodic array of straight cracks inside an infinite isotropic elastic medium and subjected to a nonuniform pressure distribution along the crack edges. The statistical properties of the random solution are evaluated numerically, and the random solution is used to determine the values of the stress-intensity factors at the crack tips. The error, expressed as the difference between the mean of the random solution and the deterministic solution, is established. Values of stress-intensity factors at the crack tip for different random input functions are presented.

Sambandham, M.↗
DOE OSTI2025

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton↗
DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗