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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A Multilevel, Hierarchical Sampling Technique for Spatially Correlated Random Fields

In this paper, we propose an alternative method to generate samples of a spatially correlated random field with applications to large-scale problems for forward propagation of uncertainty. A classical approach for generating these samples is the Karhunen--Loève (KL) decomposition. However, the KL expansion requires solving a dense eigenvalue problem and is therefore computationally infeasible for large-scale problems. Sampling methods based on stochastic partial differential equations provide a highly scalable way to sample Gaussian fields, but the resulting parametrization is mesh dependent. We propose a multilevel decomposition of the stochastic field to allow for scalable, hierarchical sampling based on solving a mixed finite element formulation of a stochastic reaction-diffusion equation with a random, white noise source function. Lastly, numerical experiments are presented to demonstrate the scalability of the sampling method as well as numerical results of multilevel Monte Carlo simulations for a subsurface porous media flow application using the proposed sampling method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a “powder-averaged” emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-E-J-215: Onboard Split PET Including the Effects of Attenuation, Scatter, and Random Events

Purpose: Positron emission tomography (PET) using split PET geometry was investigated as an on-board system for functional imaging and PET marker tracking, specifically with tomotherapy. The open dual ring PET would allow measurement of both inter and intra-fractional variation, improving the delineation of tumor volume at any stage in the radiation treatment delivery process. We present results from data obtained using Monte Carlo simulations including the effects of attenuation, random events, and scatter. Methods: PET design was accomplished via Monte Carlo simulations with GATE, the Geant4 Application for Emission Tomography. Images were reconstructed using Software for Image Reconstruction (STIR) with a fully 3D ordered subset expectation maximization (OSEM) image reconstruction technique. Results: Monte Carlo simulations of the split PET geometry indicate that including physical factors that degrade image quality such as attenuation, random events, and scatter still prove the feasibility of near real-time PET for both inter and intra-fractional radiation delivery. The image quality of scan times under 1 minute reveals that it is possible to utilize PET scanning and reconstruction during the treatment session intrafractionally. GATE also simulates the depth information but does not correct obliqueness of the path of line of response so the data is much more realistic than the data obtained with ray tracing. Conclusion: Onboard split PET with TomoTherapy can generate quality images under 1 minute scan times without the need to correct for attenuation, scatter, random events, or depth information of the interaction.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Radiation Therapy Oncology Group 0247: A Randomized Phase II Study of Neoadjuvant Capecitabine and Irinotecan or Capecitabine and Oxaliplatin With Concurrent Radiotherapy for Patients With Locally Advanced Rectal Cancer

Purpose: To evaluate the rate of pathologic complete response (pCR) and the toxicity of two neoadjuvant chemoradiotherapy (chemoRT) regimens for Stage T3-T4 rectal cancer in a randomized Phase II study. Methods and Materials: Patients with Stage T3 or T4 rectal cancer of <12 cm from the anal verge were randomized to preoperative RT (50.4 Gy in 1.8-Gy fractions) with concurrent capecitabine (1,200 mg/m{sup 2}/d Mondays through Friday) and irinotecan (50 mg/m{sup 2} weekly in four doses) (Arm 1) or concurrent capecitabine (1,650 mg/m{sup 2}/d Monday through Friday) and oxaliplatin (50 mg/m{sup 2} weekly in five doses) (Arm 2). Surgery was performed 4-8 weeks after chemoRT, and adjuvant chemotherapy 4-6 weeks after surgery. The primary endpoint was the pCR rate, requiring 48 evaluable patients per arm. Results: A total of 146 patients were enrolled. The protocol chemotherapy was modified because of excessive gastrointestinal toxicity after treatment of 35 patients; 96 were assessed for the primary endpoint-the final regimen described above. The patient characteristics were similar for both arms. After chemoRT, the rate of tumor downstaging was 52% and 60% and the rate of nodal downstaging (excluding N0 patients) was 46% and 40%, for Arms 1 and 2, respectively. The pCR rate for Arm 1 was 10% and for Arm 2 was 21%. For Arm 1 and 2, the preoperative chemoRT rate of Grade 3-4 hematologic toxicity was 9% and 4% and the rate of Grade 3-4 nonhematologic toxicity was 26% and 27%, respectively. Conclusions: Preoperative chemoRT with capecitabine plus oxaliplatin for distal rectal cancer has significant clinical activity (10 of 48 pCRs) and acceptable toxicity. This regimen is currently being evaluated in a Phase III randomized trial (National Surgical Adjuvant Breast and Bowel Project R04).

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

A Machine Learning-Based Vulnerability Analysis for Cascading Failures of Integrated Power-Gas Systems

This article proposes a cascading failure simulation (CFS) method and a hybrid machine learning method for vulnerability analysis of integrated power-gas systems (IPGSs). The CFS method is designed to study the propagating process of cascading failures between the two systems, generating data for machine learning with initial states randomly sampled. The proposed method considers generator and gas well ramping, transmission line and gas pipeline tripping, island issue handling and load shedding strategies. Then, a hybrid machine learning model with a combined random forest (RF) classification and regression algorithms is proposed to investigate the impact of random initial states on the vulnerability metrics of IPGSs. Extensive case studies are carried out on three test IPGSs to verify the proposed models and algorithms. Simulation results show that the proposed models and algorithms can achieve high accuracy for the vulnerability analysis of IPGSs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2011

Schwinger-Dyson equations in large-N quantum field theories and nonlinear random processes

We propose a stochastic method for solving Schwinger-Dyson equations in large-N quantum field theories. Expectation values of single-trace operators are sampled by stationary probability distributions of the so-called nonlinear random processes. The set of all the histories of such processes corresponds to the set of all planar diagrams in the perturbative expansions of the expectation values of singlet operators. We illustrate the method on examples of the matrix-valued scalar field theory and the Weingarten model of random planar surfaces on the lattice. For theories with compact field variables, such as sigma models or non-Abelian lattice gauge theories, the method does not converge in the physically most interesting weak-coupling limit. In this case one can absorb the divergences into a self-consistent redefinition of expansion parameters. A stochastic solution of the self-consistency conditions can be implemented as a 'memory' of the random process, so that some parameters of the process are estimated from its previous history. We illustrate this idea on the two-dimensional O(N) sigma model. The extension to non-Abelian lattice gauge theories is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Development and application of a standardized flow measurement uncertainty analysis framework to various low-head short-converging intake types across the United States federal hydropower fleet

Turbine discharges at low-head short converging intakes are difficult to measure accurately. The proximity of the measurement section to the intake entrance admits large uncertainties related to asymmetry of the velocity profile, swirl, and turbulence. Existing turbine performance codes [10, 24] do not address this special case and published literature is largely silent on rigorous evaluation of uncertainties associated with this measurement context. The American Society of Mechanical Engineers (ASME) Committee investigated the use of Acoustic transit time (ATT), Acoustic scintillation (AS), and Current meter (CM) in a short converging intake at the Kootenay Canal Generating Station in 2009. Based on their findings, a standardized uncertainty analysis (UA) framework for velocity-area method (specifically for CM measurements) is presented in this paper given the fact that CM is still the most fundamental and common type of measurement system. Typical sources of systematic and random errors associated with CM measurements are investigated, and the major sources of uncertainties associated with turbulence and velocity fluctuations, numerical velocity integration technique (bi-cubic spline), and the number and placement of current meters are being considered for an evaluation. Since the velocity measurements in a short converging intake are associated with complex nonlinear and time varying uncertainties (e.g., Reynolds stress in fluid dynamics), simply applying the law of propagation of uncertainty is known to overestimate the measurement variance while the Monte Carlo method does not. Therefore, a pseudo-Monte Carlo simulation method (random flow generation technique [8]) which was initially developed for the purpose of establishing upstream or initial conditions in the Large-Eddy Simulation (LES) and the Direct Numerical Simulation (DNS) is used to statistically determine uncertainties associated with turbulence and velocity fluctuations. This technique is then combined with a bi-cubic spline interpolation method which converts point velocities into a continuous velocity distribution over the measurement domain. Subsequently the number and placement of current meters are simulated to investigate the accuracy of the estimated flow rates using the numerical velocity-area integration method outlined in ISO 3354 [12]. The authors herein consider that statistics on generated flow rates processed with bi-cubic interpolation and sensor simulations are the combined uncertainties which already accounted for the effects of all those three uncertainty sources. A preliminary analysis based on the current meter data obtained through an upgrade acceptance test of a single unit located in a mainstem plant has been presented.

Smith, Brennan T↗
DOE OSTI2018

On the Stability of a Linear System of Difference Equations with Random Parameters

We study the asymptotic behavior of solutions to a linear system of difference equation whose right-hand side at each time moment depends not only on the value at the previous moment, but also on a random parameter that takes its values in a given set. We obtain conditions of the Lyapunov stability and the asymptotic stability of the equilibrium position that are valid for all values of the random parameter or with probability 1.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Poly-Gaussian model of randomly rough surface in rarefied gas flow

Surface roughness is simulated by the model of non-Gaussian random process. Our results for the scattering of rarefied gas atoms from a rough surface using modified approach to the DSMC calculation of rarefied gas flow near a rough surface are developed and generalized applying the poly-Gaussian model representing probability density as the mixture of Gaussian densities. The transformation of the scattering function due to the roughness is characterized by the roughness operator. Simulating rough surface of the walls by the poly-Gaussian random field expressed as integrated Wiener process, we derive a representation of the roughness operator that can be applied in numerical DSMC methods as well as in analytical investigations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Updated Results and Patterns of Failure in a Randomized Hypofractionation Trial for High-risk Prostate Cancer

Purpose: To report long-term results and patterns of failure after conventional and hypofractionated radiation therapy in high-risk prostate cancer. Methods and Materials: This randomized phase III trial compared conventional fractionation (80 Gy at 2 Gy per fraction in 8 weeks) vs hypofractionation (62 Gy at 3.1 Gy per fraction in 5 weeks) in combination with 9-month androgen deprivation therapy in 168 patients with high-risk prostate cancer. Freedom from biochemical failure (FFBF), freedom from local failure (FFLF), and freedom from distant failure (FFDF) were analyzed. Results: In a median follow-up of 70 months, biochemical failure (BF) occurred in 35 of the 168 patients (21%) in the study. Among these 35 patients, local failure (LF) only was detected in 11 (31%), distant failure (DF) only in 16 (46%), and both LF and DF in 6 (17%). In 2 patients (6%) BF has not yet been clinically detected. The risk reduction by hypofractionation was significant in BF (10.3%) but not in LF and DF. We found that hypofractionation, with respect to conventional fractionation, determined only an insignificant increase in the actuarial FFBF but no difference in FFLF and FFDF, when considering the entire group of patients. However, an increase in the 5-year rates in all 3 endpoints-FFBF, FFLF, and FFDF-was observed in the subgroup of patients with a pretreatment prostate-specific antigen (iPSA) level of 20 ng/mL or less. On multivariate analysis, the type of fractionation, iPSA level, Gleason score of 4+3 or higher, and T stage of 2c or higher have been confirmed as independent prognostic factors for BF. High iPSA levels and Gleason score of 4+3 or higher were also significantly associated with an increased risk of DF, whereas T stage of 2c or higher was the only independent variable for LF. Conclusion: Our results confirm the isoeffectiveness of the 2 fractionation schedules used in this study, although a benefit in favor of hypofractionation cannot be excluded in the subgroup of patients with an iPSA level of 20 ng/mL or less. The {alpha}/{beta} ratio might be more appropriately evaluated by FFLF than FFBF results, at least in high-risk disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2010

Characterization of electron microscopes with binary pseudo-random multilayer test samples

Verification of the reliability of metrology data from high quality X-ray optics requires that adequate methods for test and calibration of the instruments be developed. For such verification for optical surface profilometers in the spatial frequency domain, a modulation transfer function (MTF) calibration method based on binarypseudo-random (BPR) gratings and arrays has been suggested and and proven to be an effective calibration method for a number of interferometric microscopes, a phase shifting Fizeau interferometer, and a scatterometer. Here we describe the details of development of binarypseudo-random multilayer (BPRML) test samples suitable for characterization of scanning (SEM) and transmission (TEM) electron microscopes. We discuss the results of TEM measurements with the BPRML test samples fabricated from a WiSi{sub 2}/Si multilayer coating with pseudo-randomly distributed layers. In particular, we demonstrate that significant information about the metrological reliability of the TEM measurements can be extracted even when the fundamental frequency of the BPRML sample is smaller than the Nyquist frequency of the measurements. The measurements demonstrate a number of problems related to the interpretation of the SEM and TEM data. Note that similar BPRML testsamples can be used to characterize X-ray microscopes. Corresponding work with X-ray microscopes is in progress.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Formal Estimation of Errors in Computed Absolute Interaction Energies of Protein-Ligand Complexes

A largely unsolved problem in computational biochemistry is the accurate prediction of binding affinities of small ligands to protein receptors. We present a detailed analysis of the systematic and random errors present in computational methods through the use of error probability density functions, specifically for computed interaction energies between chemical fragments comprising a protein-ligand complex. An HIV-II protease crystal structure with a bound ligand (indinavir) was chosen as a model protein-ligand complex. The complex was decomposed into 21 interacting fragment pairs, which were studied using a number of computational methods. The chemically accurate complete basis set coupled cluster theory (CCSD(T)/CBS) interaction energies were used as reference values to generate our error estimates. In our analysis, we observed significant systematic and random errors in most methods, which was surprising, especially for parametrized classical and semiempirical quantum mechanical calculations. After propagating these fragment-based error estimates over the entire protein-ligand complex, our total error estimates for many methods are large compared to the experimentally determined free energy of binding. Thus, we conclude that statistical error analysis is a necessary addition to any scoring function attempting to produce reliable binding affinity predictions.

Chemistry↗
DOE OSTI2015

A Robust and Versatile Method of Combinatorial Chemical Synthesis of Gene Libraries via Hierarchical Assembly of Partially Randomized Modules

A major challenge in gene library generation is to guarantee a large functional size and diversity that significantly increases the chances of selecting different functional protein variants. The use of trinucleotides mixtures for controlled randomization results in superior library diversity and offers the ability to specify the type and distribution of the amino acids at each position. Here we describe the generation of a high diversity gene library using tHisF of the hyperthermophile Thermotoga maritima as a scaffold. Combining various rational criteria with contingency, we targeted 26 selected codons of the thisF gene sequence for randomization at a controlled level. We have developed a novel method of creating full-length gene libraries by combinatorial assembly of smaller sub-libraries. Full-length libraries of high diversity can easily be assembled on demand from smaller and much less diverse sublibraries, which circumvent the notoriously troublesome long-term archivation and repeated proliferation of high diversity ensembles of phages or plasmids. We developed a generally applicable software tool for sequence analysis of mutated gene sequences that provides efficient assistance for analysis of library diversity. Finally, practical utility of the library was demonstrated in principle by assessment of the conformational stability of library members and isolating protein variants with HisF activity from it. Our approach integrates a number of features of nucleic acids synthetic chemistry, biochemistry and molecular genetics to a coherent, flexible and robust method of combinatorial gene synthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Monte Carlo based toy model for fission process

There are many models and calculation techniques to obtain visible image of fission yield process. In particular, fission yield can be calculated by using two calculations approach, namely macroscopic approach and microscopic approach. This work proposes another calculation approach in which the nucleus is treated as a toy model. Hence, the fission process does not represent real fission process in nature completely. The toy model is formed by Gaussian distribution of random number that randomizes distance like the distance between particle and central point. The scission process is started by smashing compound nucleus central point into two parts that are left central and right central points. These three points have different Gaussian distribution parameters such as mean (μ{sub CN}, μ{sub L}, μ{sub R}), and standard deviation (σ{sub CN}, σ{sub L}, σ{sub R}). By overlaying of three distributions, the number of particles (N{sub L}, N{sub R}) that are trapped by central points can be obtained. This process is iterated until (N{sub L}, N{sub R}) become constant numbers. Smashing process is repeated by changing σ{sub L} and σ{sub R}, randomly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Calibration of the modulation transfer function of surface profilometers with binary pseudo-random test standards: expanding the application range

A modulation transfer function (MTF) calibration method based on binary pseudo-random (BPR) gratings and arrays [Proc. SPIE 7077-7 (2007), Opt. Eng. 47, 073602 (2008)] has been proven to be an effective MTF calibration method for a number of interferometric microscopes and a scatterometer [Nucl. Instr. and Meth. A616, 172 (2010)]. Here we report on a further expansion of the application range of the method. We describe the MTF calibration of a 6 inch phase shifting Fizeau interferometer. Beyond providing a direct measurement of the interferometer's MTF, tests with a BPR array surface have revealed an asymmetry in the instrument's data processing algorithm that fundamentally limits its bandwidth. Moreover, the tests have illustrated the effects of the instrument's detrending and filtering procedures on power spectral density measurements. The details of the development of a BPR test sample suitable for calibration of scanning and transmission electron microscopes are also presented. Such a test sample is realized as a multilayer structure with the layer thicknesses of two materials corresponding to BPR sequence. The investigations confirm the universal character of the method that makes it applicable to a large variety of metrology instrumentation with spatial wavelength bandwidths from a few nanometers to hundreds of millimeters.

42 ENGINEERING↗
DOE OSTI2011

Multilevel Method for Modeling Large-Scale Networks

Understanding the behavior of real complex networks is of great theoretical and practical significance. It includes developing accurate artificial models whose topological properties are similar to the real networks, generating the artificial networks at different scales under special conditions, investigating a network dynamics, reconstructing missing data, predicting network response, detecting anomalies and other tasks. Network generation, reconstruction, and prediction of its future topology are central issues of this field. In this project, we address the questions related to the understanding of the network modeling, investigating its structure and properties, and generating artificial networks. Most of the modern network generation methods are based either on various random graph models (reinforced by a set of properties such as power law distribution of node degrees, graph diameter, and number of triangles) or on the principle of replicating an existing model with elements of randomization such as R-MAT generator and Kronecker product modeling. Hierarchical models operate at different levels of network hierarchy but with the same finest elements of the network. However, in many cases the methods that include randomization and replication elements on the finest relationships between network nodes and modeling that addresses the problem of preserving a set of simplified properties do not fit accurately enough the real networks. Among the unsatisfactory features are numerically inadequate results, non-stability of algorithms on real (artificial) data, that have been tested on artificial (real) data, and incorrect behavior at different scales. One reason is that randomization and replication of existing structures can create conflicts between fine and coarse scales of the real network geometry. Moreover, the randomization and satisfying of some attribute at the same time can abolish those topological attributes that have been undefined or hidden from researchers. We propose to develop multilevel methods to model complex networks. The key point of the proposed strategy is that it will help to preserve part of the unknown structural attributes by guaranteeing the similar behavior of the real and artificial model on different scales.

99 GENERAL AND MISCELLANEOUS↗