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
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↗
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↗
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↗
DOE OSTI2013

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

Dispersion in porous media flows has been the subject of much experimental, theoretical and numerical study. Here we consider a wavy-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media, where constrictions represent pore throats and expansions pore bodies. A theoretical model for effective (macroscopic) longitudinal dispersion in this system has been developed by volume averaging the microscale velocity field. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier-Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re = 449 for which unsteady flow was observed. Dispersion values were computed using both the volume averaging solution and a random walk particle tracking method, and results of the two methods were shown to be consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for the low-Re, Stokes flow regime. In the steady inertial regime we observe an power-law increase in effective longitudinal dispersion (DL) with Re, consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). For the unsteady case (Re = 449), the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion.

CFD↗
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 OSTI2016

A biorthogonal decomposition for the identification and simulation of non-stationary and non-Gaussian random fields

In this paper, a new method for the identification and simulation of non-Gaussian and non-stationary stochastic fields given a database is proposed. It is based on two successive biorthogonal decompositions aiming at representing spatio–temporal stochastic fields. The proposed double expansion allows to build the model even in the case of large-size problems by separating the time, space and random parts of the field. A Gaussian kernel estimator is used to simulate the high dimensional set of random variables appearing in the decomposition. The capability of the method to reproduce the non-stationary and non-Gaussian features of random phenomena is illustrated by applications to earthquakes (seismic ground motion) and sea states (wave heights).

58 GEOSCIENCES↗
DOE OSTI2013

Solution of stochastic media transport problems using a numerical quadrature-based method

We present a new conceptual framework for analyzing transport problems in random media. We decompose such problems into stratified subproblems according to the number of material pseudo-interfaces within realizations. For a given subproblem we assign pseudo-interface locations in each realization according to product quadrature rules, which allows us to deterministically generate a fixed number of realizations. Quadrature integration of the solutions of these realizations thus approximately solves each subproblem; the weighted superposition of solutions of the subproblems approximately solves the general stochastic media transport problem. We revisit some benchmark problems to determine the accuracy and efficiency of this approach in comparison to randomly generated realizations. We find that this method is very accurate and fast when the number of pseudo-interfaces in a problem is generally low, but that these advantages quickly degrade as the number of pseudo-interfaces increases. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Probing the astrophysical origin of high-energy cosmic-ray electrons with Monte Carlo simulation

High-energy cosmic-ray electrons reveal some remarkable spectral features, the most noteworthy of which is the rise in the positron fraction above 10 GeV. Due to strong energy loss during propagation, these particles can reach Earth only from nearby sources. Yet, the exact nature of these sources, which most likely manifest themselves in the observed anomalies, remains elusive. The many explanations put forward to resolve this case range from standard astrophysics to exotic physics. In this paper, we discuss the possible astrophysical origin of high-energy cosmic-ray electrons through a fully three-dimensional time-dependent Monte Carlo simulation. This approach, which takes advantage of the intrinsic random nature of cosmic-ray diffusive propagation, provides valuable information on the electron-by-electron fluctuations, making it particularly suitable for analyzing in depth the single-source scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Comparing the white dwarf cooling sequences in 47 Tuc and NGC 6397

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci—a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Quantification of Organ Motion During Chemoradiotherapy of Rectal Cancer Using Cone-Beam Computed Tomography

Purpose: There has been no previously published data related to the quantification of rectal motion using cone-beam computed tomography (CBCT) during standard conformal long-course chemoradiotherapy. The purpose of the present study was to quantify the interfractional changes in rectal movement and dimensions and rectal and bladder volume using CBCT and to quantify the bony anatomy displacements to calculate the margins required to account for systematic ({Sigma}) and random ({sigma}) setup errors. Methods and Materials: CBCT images were acquired from 16 patients on the first 3 days of treatment and weekly thereafter. The rectum and bladder were outlined on all CBCT images. The interfraction movement was measured using fixed bony landmarks as references to define the rectal location (upper, mid, and low), The maximal rectal diameter at the three rectal locations was also measured. The bony anatomy displacements were quantified, allowing the calculation of systematic ({Sigma}) and random ({sigma}) setup errors. Results: A total of 123 CBCT data sets were analyzed. Analysis of variance for standard deviation from planning scans showed that rectal anterior and lateral wall movement differed significantly by rectal location. Anterior and lateral rectal wall movements were larger in the mid and upper rectum compared with the low rectum. The posterior rectal wall movement did not change significantly with the rectal location. The rectal diameter changed more in the mid and upper than in the low rectum. No consistent relationship was found between the rectal and bladder volume and time, nor was a significant relationship found between the rectal volume and bladder volume. Conclusions: In the present study, the anterior and lateral rectal movement and rectal diameter were found to change most in the upper rectum, followed by the mid rectum, with the smallest changes seen in the low rectum. Asymmetric margins are warranted to ensure phase 2 coverage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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↗
DOE OSTI2023

Designing an Intrusion Detection for an Adjustable Speed Drive System Controlling a Critical Process

In this article, we address the cyber-security problem of industrial control systems (ICSs) when their sensor measurements may be compromised due to an attacker who has intercepted those measurements via a network. We introduce a general-purpose method “Dynamic Watermarking (DW)” to detect potential cyber-intrusions on speed sensor measurements within industrial control systems, which deploy an adjustable speed drive (ASD) to control a critical process. The DW method is injecting a random private low-amplitude signal with a zero mean Gaussian distribution, “watermark”, into one of the input phase voltages powering the ASD system. The watermark signal propagates through the system including pulse width modulation (PWM) power conversion stage and motor, then ultimately appears in the speed sensor measurements. By deploying two statistical DW tests with two proper thresholds, the system can detect potential cyber-intrusions or unobservable cyber-attacks such as replay attacks and false data injection attacks (FDIA). The DW method tested on a laboratory-scale ASD system experimentally to protect the system against cyber-intrusions. This system, powered by a commercial PWM drive operating at 208 V, 3-phase, and 3.7 kW, served as our experimental platform.

42 ENGINEERING↗
DOE OSTI2013

Cubature on Wiener Space: Pathwise Convergence

Cubature on Wiener space (Lyons and Victoir in Proc. R. Soc. Lond. A 460(2041):169-198, 2004) provides a powerful alternative to Monte Carlo simulation for the integration of certain functionals on Wiener space. More specifically, and in the language of mathematical finance, cubature allows for fast computation of European option prices in generic diffusion models.We give a random walk interpretation of cubature and similar (e.g. the Ninomiya-Victoir) weak approximation schemes. By using rough path analysis, we are able to establish weak convergence for general path-dependent option prices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

An atomistic methodology of energy release rate for graphene at nanoscale

Graphene is a single layer of carbon atoms packed into a honeycomb architecture, serving as a fundamental building block for electric devices. Understanding the fracture mechanism of graphene under various conditions is crucial for tailoring the electrical and mechanical properties of graphene-based devices at atomic scale. Although most of the fracture mechanics concepts, such as stress intensity factors, are not applicable in molecular dynamics simulation, energy release rate still remains to be a feasible and crucial physical quantity to characterize the fracture mechanical property of materials at nanoscale. This work introduces an atomistic simulation methodology, based on the energy release rate, as a tool to unveil the fracture mechanism of graphene at nanoscale. This methodology can be easily extended to any atomistic material system. We have investigated both opening mode and mixed mode at different temperatures. Simulation results show that the critical energy release rate of graphene is independent of initial crack length at low temperature. Graphene with inclined pre-crack possesses higher fracture strength and fracture deformation but smaller critical energy release rate compared with the graphene with vertical pre-crack. Owing to its anisotropy, graphene with armchair chirality always has greater critical energy release rate than graphene with zigzag chirality. The increase of temperature leads to the reduction of fracture strength, fracture deformation, and the critical energy release rate of graphene. Also, higher temperature brings higher randomness of energy release rate of graphene under a variety of predefined crack lengths. The energy release rate is independent of the strain rate as long as the strain rate is small enough.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗