Search NASA⌕ Search

SEARCH · Search NASA

Results for “Randomized methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Advanced reliability method for fatigue analysis

When design factors are considered as random variables and the failure condition cannot be expressed by a closed form algebraic inequality, computations of risk (or probability of failure) may become extremely difficult or very inefficient. This study suggests using a simple and easily constructed second degree polynomial to approximate the complicated limit state in the neighborhood of the design point; a computer analysis relates the design variables at selected points. Then a fast probability integration technique (i.e., the Rackwitz-Fiessler algorithm) can be used to estimate risk. The capability of the proposed method is demonstrated in an example of a low cycle fatigue problem for which a computer analysis is required to perform local strain analysis to relate the design variables. A comparison of the performance of this method is made with a far more costly Monte Carlo solution. Agreement of the proposed method with Monte Carlo is considered to be good.

Wu, Y.-T.↗
DOE OSTI2012

Randomized Control Trial: Evaluating Aluminum-Based Antiperspirant Use, Axilla Skin Toxicity, and Reported Quality of Life in Women Receiving External Beam Radiotherapy for Treatment of Stage 0, I, and II Breast Cancer

Purpose: Standard skin care instructions regarding the use of antiperspirants during radiotherapy to the breast varies across North America. Women have articulated that when instructed to not use antiperspirant, the potential for body odor is distressing. Historical practices and individual opinions have often guided practice in this field. The present study had 2 purposes. To evaluate whether the use of aluminum-based antiperspirant while receiving external beam radiotherapy for stage 0, I, or II breast cancer will increase axilla skin toxicity and to evaluate whether the use of antiperspirant during external beam radiotherapy improves quality of life. Methods: A total of 198 participants were randomized to either the experimental group (antiperspirant) or control group (standard care-wash only). The skin reactions in both groups were measured weekly and 2 weeks after treatment using the National Cancer Institute Common Toxicity Criteria Adverse Events, version 3, toxicity grading criteria. Both groups completed the Functional Assessment for Chronic Illness Therapy's questionnaire for the breast population quality of life assessment tool, with additional questions evaluating the effect of underarm antiperspirant use on quality of life before treatment, immediately after treatment, and 2 weeks after treatment during the study. Results: The skin reaction data were analyzed using the generalized estimating equation. No statistically significant difference was seen in the skin reaction between the 2 groups over time. The quality of life data also revealed no statistically significant difference between the 2 groups over time. Conclusions: Data analysis indicates that using antiperspirant routinely during external beam radiotherapy for Stage 0, I, or II breast cancer does not affect the intensity of the skin reaction or the self-reported quality of life. This evidence supports that in this particular population, there is no purpose to restrict these women from using antiperspirants during their treatment, and the decision to use an antiperspirant or not in this setting should be left to the discretion of the patient.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1978

A comparative study of Conroy and Monte Carlo methods applied to multiple quadratures and multiple scattering

An efficient numerical method of multiple quadratures, the Conroy method, is applied to the problem of computing multiple scattering contributions in the radiative transfer through realistic planetary atmospheres. A brief error analysis of the method is given and comparisons are drawn with the more familiar Monte Carlo method. Both methods are stochastic problem-solving models of a physical or mathematical process and utilize the sampling scheme for points distributed over a definite region. In the Monte Carlo scheme the sample points are distributed randomly over the integration region. In the Conroy method, the sample points are distributed systematically, such that the point distribution forms a unique, closed, symmetrical pattern which effectively fills the region of the multidimensional integration. The methods are illustrated by two simple examples: one, of multidimensional integration involving two independent variables, and the other, of computing the second order scattering contribution to the sky radiance.

Deepak, A.↗
DOE OSTI2012

Role of IVC Filters in Endovenous Therapy for Deep Venous Thrombosis: The FILTER-PEVI (Filter Implantation to Lower Thromboembolic Risk in Percutaneous Endovenous Intervention) Trial

Objectives: The purpose of this study was to evaluate the necessity of and recommend indications for inferior vena cava (IVC) filter implantation during percutaneous endovenous intervention (PEVI) for deep venous thrombosis (DVT).BackgroundPEVI has emerged as a powerful tool in the management of acute proximal DVT. Instrumentation of extensive fresh thrombus is potentially associated with iatrogenic pulmonary embolism (PE). The true frequency of this complication has not been studied in a randomized fashion. We evaluated IVC filter implantation during PEVI for DVT. Methods: A total of 141 patients with symptomatic proximal DVT undergoing PEVI for symptomatic DVT were randomized to receive an IVC filter (70 patients) or no filter (71 patients; control group). The anticoagulation and PEVI regimen were similar between the two groups. Patients with development of symptoms suggestive of PE underwent objective testing for PE. Results: PE developed in 1 of the 14 symptomatic patients in the filter group and 8 of the 22 patients in the control group (P = 0.048). There was no mortality in any group. Three patients (4.2%) in the control group had transient hemodynamic instability necessitating resuscitory efforts. Predictors of iatrogenic PE were found to be PE at admission; involvement of two or more adjacent venous segments with acute thrombus; inflammatory form of DVT (severe erythema, edema, pain, and induration); and vein diameter of {>=}7 mm with preserved architecture. Conclusions: IVC filter implantation during PEVI reduces the risk of iatrogenic PE by eightfold without a mortality benefit. A selective approach may be exercised in filter implantation during PEVI.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS2003

Evaluation of Treadmill Exercise in a Lower Body Negative Pressure Chamber as a Countermeasure for Weightlessness-Induced Bone Loss: a Bed Rest Study with Identical Twins

Counteracting bone loss is required for future space exploration. We evaluated the ability of treadmill exercise in a LBNP chamber to counteract bone loss in a 30-day bed rest study. Eight pairs of identical twins were randomly assigned to sedentary control or exercise groups. Exercise within LBNP decreased the bone resorption caused by bed rest and may provide a countermeasure for spaceflight. INTRODUCTION: Bone loss is one of the greatest physiological challenges for extended-duration space missions. The ability of exercise to counteract weightlessness-induced bone loss has been studied extensively, but to date, it has proven ineffective. We evaluated the effectiveness of a combination of two countermeasures-treadmill exercise while inside a lower body negative pressure (LBNP) chamber-on bone loss during a 30-day bed rest study. MATERIALS AND METHODS: Eight pairs of identical twins were randomized into sedentary (SED) or exercise/LBNP (EX/LBNP) groups. Blood and urine samples were collected before, several times during, and after the 30-day bed rest period. These samples were analyzed for markers of bone and calcium metabolism. Repeated measures ANOVA was used to determine statistical significance. Because identical twins were used, both time and group were treated as repeated variables. RESULTS: Markers of bone resorption were increased during bed rest in samples from sedentary subjects, including the collagen cross-links and serum and urinary calcium concentrations. For N-telopeptide and deoxypyridinoline, there were significant (p < 0.05) interactions between group (SED versus EX/LBNP) and phase of the study (sample collection point). Pyridinium cross-links were increased above pre-bed rest levels in both groups, but the EX/LBNP group had a smaller increase than the SED group. Markers of bone formation were unchanged by bed rest in both groups. CONCLUSIONS: These data show that this weight-bearing exercise combined with LBNP ameliorates some of the negative effects of simulated weightlessness on bone metabolism. This protocol may pave the way to counteracting bone loss during spaceflight and may provide valuable information about normal and abnormal bone physiology here on Earth.

Twin Study↗
DOE OSTI2016

A generalized multi-resolution expansion for uncertainty propagation with application to cardiovascular modeling

Computational models are used in a variety of fields to improve our understanding of complex physical phenomena. Recently, the realism of model predictions has been greatly enhanced by transitioning from deterministic to stochastic frameworks, where the effects of the intrinsic variability in parameters, loads, constitutive properties, model geometry and other quantities can be more naturally included. A general stochastic system may be characterized by a large number of arbitrarily distributed and correlated random inputs, and a limited support response with sharp gradients or event discontinuities. This motivates continued research into novel adaptive algorithms for uncertainty propagation, particularly those handling high dimensional, arbitrarily distributed random inputs and non-smooth stochastic responses. Here, we generalize a previously proposed multi-resolution approach to uncertainty propagation to develop a method that improves computational efficiency, can handle arbitrarily distributed random inputs and non-smooth stochastic responses, and naturally facilitates adaptivity, i.e., the expansion coefficients encode information on solution refinement. Our approach relies on partitioning the stochastic space into elements that are subdivided along a single dimension, or, in other words, progressive refinements exhibiting a binary tree representation. We also show how these binary refinements are particularly effective in avoiding the exponential increase in the multi-resolution basis cardinality and significantly reduce the regression complexity for moderate to high dimensional random inputs. The performance of the approach is demonstrated through previously proposed uncertainty propagation benchmarks and stochastic multi-scale finite element simulations in cardiovascular flow.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2013

UTILIZATION OF THE WAVEFRONT SENSOR AND SHORT-EXPOSURE IMAGES FOR SIMULTANEOUS ESTIMATION OF QUASI-STATIC ABERRATION AND EXOPLANET INTENSITY

Heretofore, the literature on exoplanet detection with coronagraphic telescope systems has paid little attention to the information content of short exposures and methods of utilizing the measurements of adaptive optics wavefront sensors. This paper provides a framework for the incorporation of the wavefront sensor measurements in the context of observing modes in which the science camera takes millisecond exposures. In this formulation, the wavefront sensor measurements provide a means to jointly estimate the static speckle and the planetary signal. The ability to estimate planetary intensities in as little as a few seconds has the potential to greatly improve the efficiency of exoplanet search surveys. For simplicity, the mathematical development assumes a simple optical system with an idealized Lyot coronagraph. Unlike currently used methods, in which increasing the observation time beyond a certain threshold is useless, this method produces estimates whose error covariances decrease more quickly than inversely proportional to the observation time. This is due to the fact that the estimates of the quasi-static aberrations are informed by a new random (but approximately known) wavefront every millisecond. The method can be extended to include angular (due to diurnal field rotation) and spectral diversity. Numerical experiments are performed with wavefront data from the AEOS Adaptive Optics System sensing at 850 nm. These experiments assume a science camera wavelength {lambda} of 1.1 {mu}, that the measured wavefronts are exact, and a Gaussian approximation of shot-noise. The effects of detector read-out noise and other issues are left to future investigations. A number of static aberrations are introduced, including one with a spatial frequency exactly corresponding the planet location, which was at a distance of Almost-Equal-To 3{lambda}/D from the star. Using only 4 s of simulated observation time, a planetary intensity, of Almost-Equal-To 1 photon ms{sup -1}, and a stellar intensity of Almost-Equal-To 10{sup 5} photons ms{sup -1} (contrast ratio 10{sup 5}), the short-exposure estimation method recovers the amplitudes' static aberrations with 1% accuracy, and the planet brightness with 20% accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Parameter estimation for X-ray scattering analysis with Hamiltonian Markov Chain Monte Carlo

Bayesian-inference-based approaches, in particular the random-walk Markov Chain Monte Carlo (MCMC) method, have received much attention recently for X-ray scattering analysis. Hamiltonian MCMC, a state-of-the-art development in the field of MCMC, has become popular in recent years. It utilizes Hamiltonian dynamics for indirect but much more efficient drawings of the model parameters. We described the principle of the Hamiltonian MCMC for inversion problems in X-ray scattering analysis by estimating high-dimensional models for several motivating scenarios in small-angle X-ray scattering, reflectivity, and X-ray fluorescence holography. Hamiltonian MCMC with appropriate preconditioning can deliver superior performance over the random-walk MCMC, and thus can be used as an efficient tool for the statistical analysis of the parameter distributions, as well as model predictions and confidence analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Semimicroscopic description of basic modes of relaxation of the Gamow-Teller resonance

Semimicroscopic approach based on the continuum version of the random-phase approximation and on a phenomenological method for taking into account the fragmentation effect is used to describe quantitatively basic relaxation parameters of the Gamow-Teller resonance in magic and nearmagic nuclei. For the {sup 208}Bi nucleus, the results obtained by calculating the resonance parameters in question are compared with respective experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Girsanov's transformation based variance reduced Monte Carlo simulation schemes for reliability estimation in nonlinear stochastic dynamics

The study considers the problem of simulation based time variant reliability analysis of nonlinear randomly excited dynamical systems. Attention is focused on importance sampling strategies based on the application of Girsanov's transformation method. Controls which minimize the distance function, as in the first order reliability method (FORM), are shown to minimize a bound on the sampling variance of the estimator for the probability of failure. Two schemes based on the application of calculus of variations for selecting control signals are proposed: the first obtains the control force as the solution of a two-point nonlinear boundary value problem, and, the second explores the application of the Volterra series in characterizing the controls. The relative merits of these schemes, vis-à-vis the method based on ideas from the FORM, are discussed. Illustrative examples, involving archetypal single degree of freedom (dof) nonlinear oscillators, and a multi-degree of freedom nonlinear dynamical system, are presented. The credentials of the proposed procedures are established by comparing the solutions with pertinent results from direct Monte Carlo simulations. - Highlights: • The distance minimizing control forces minimize a bound on the sampling variance. • Establishing Girsanov controls via solution of a two-point boundary value problem. • Girsanov controls via Volterra's series representation for the transfer functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Graphite Characterization: Baseline Variability Analysis Report

The Idaho National INL’s Baseline Material Property program provides many measurements of different material properties on a variety of nuclear graphite grades being investigated for irradiation-induced changes in material properties. The baseline material properties provide an unirradiated reference for comparison with measurements of those properties after irradiation. This study describes how uncertainty, or variance, in estimates of important material properties, like compressive strength, may depend on sample size, sample distribution and parameter estimation methods. Compressive strength data from two large collections of specimens from a billet of PCEA graphite (sample size = 230) and a billet of IG-110 graphite (sample size = 48) are used to illustrate these dependencies. The grades of graphite used represent approximate end members in terms of relative variability in strength measurements. Monte Carlo simulation and bootstrapping methods were used to examine variance in Weibull-parameter estimates as a function of sample size. Results were consistent with what previous statistical analyses have shown about the relationship between sample size and width of the confidence interval for maximum likelihood (ML) estimated parameters when samples are repeatedly drawn from a Weibull distribution. As a function of sample size, the precision of the modulus estimate increases faster than that of the characteristic strength. The modulus maximum likelihood estimates (MLEs) are biased for low sample sizes. The American Society for Testing and Materials document ASTM D7846-16 provides equations that describe the dependence of the variance of parameter MLEs on sample size. Using the same methodology, we show that the polynomials provided for variance dependence on sample size are inaccurate for sample sizes greater than ~120. We provide approximations to that relationship that illustrate relative constancy for sample size > ~10, in the logspace relationship of the modulus variance and the log transformed characteristic strength variance. Bootstrapping provided a means of comparing Weibull parameters obtained using ASME’s guidance for specimen collection to those obtained using different weighting in a case where strength is related to position in a billet and to grain orientation and a single distribution is used to characterize the billet. ASME specifies that samples be collected with equal representation in all locations and grain orientations. Parameter values obtained using ASME’s method were different from those obtained using random sampling of the complete data sets, and this was particularly true for the PCEA data. In addition, while the ASME approach yielded more conservative parameters for the PCEA data, the opposite was true for the IG-110 data. The difference observed for the PCEA data set likely reflects the larger number of against-grain specimens collected, which had generally higher strength than the with-grain specimens. This emphasizes the importance of the dependence of strength on location or grain orientation and suggests that it may be prudent to weight specimen collection toward the weaker subgroups if a conservative estimate as desired, as the equal weighting method recommended by ASME would will not always provide the most conservative Weibull parameters. Because spatial correlations commonly exist in certain grades of graphite, future work may involve developing a methodology for parameter inference that relaxes the assumption of independent specimens. ASTM and ASME codes recommend different parameter estimation methods and models of the Weibull distribution in different situations. Descriptions of the effects on parameter estimates for the example billets used in this study provide an example of how the method may affect the parameter estimate, and how—in this case—that compares to differences due to other factors, such as sample size and dependence on grain orientation or sample location.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Resolution of 100 photons and quantum generation of unbiased random numbers

Macroscopic quantum phenomena, such as observed in superfluids and superconductors, have led to promising technological advancements and some of the most important tests of fundamental physics. At present, quantum detection of light is mostly relegated to the microscale, where avalanche photodiodes are very sensitive to distinguishing single-photon events from vacuum but cannot differentiate between larger photon-number events. Beyond this, the ability to perform measurements to resolve photon numbers is highly desirable for a variety of quantum information applications including computation, sensing, and cryptography. True photon-number resolving detectors do exist, but they are currently limited to the ability to resolve on the order of 10 photons, which is too small for certain proposals. In this work, we extend photon measurement into the mesoscopic regime by implementing a detection scheme based on multiplexing highly quantum-efficient transition-edge sensors to accurately resolve photon numbers between zero and 100. Further, we then demonstrate the use of our system by implementing a quantum random number generator with no inherent bias. This method is based on sampling a coherent state in the photon-number basis and is robust against environmental noise, phase and amplitude fluctuations in the laser, loss and detector inefficiency as well as eavesdropping. Beyond true random number generation, our detection scheme serves as a means to implement quantum measurement and engineering techniques valuable for photonic quantum information processing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Comparison of Prophylactic Naftopidil, Tamsulosin, and Silodosin for {sup 125}I Brachytherapy-Induced Lower Urinary Tract Symptoms in Patients With Prostate Cancer: Randomized Controlled Trial

Purpose: To compare the efficacy of three {alpha}{sub 1A}/{alpha}{sub 1D}-adrenoceptor (AR) antagonists-naftopidil, tamsulosin, and silodosin-that have differing affinities for the {alpha}{sub 1}-AR subtypes in treating urinary morbidities in Japanese men with {sup 125}I prostate implantation (PI) for prostate cancer. Methods and Materials: This single-institution prospective randomized controlled trial compared naftopidil, tamsulosin, and silodosin in patients undergoing PI. Patients were randomized and received either naftopidil, tamsulosin, or silodosin. Treatment began 1 day after PI and continued for 1 year. The primary efficacy variables were the changes in total International Prostate Symptom Score (IPSS) and postvoid residual urine (PVR). The secondary efficacy variables were changes in IPSS storage score and IPSS voiding score from baseline to set points during the study (1, 3, 6, and 12 months). Results: Two hundred twelve patients were evaluated in this study between June 2006 and February 2009: 71, 70, and 71 patients in the naftopidil, tamsulosin, and silodosin groups, respectively. With respect to the primary efficacy variables, the mean changes in the total IPSS at 1 month after PI in the naftopidil, tamsulosin, and silodosin groups were +10.3, +8.9, and +7.5, respectively. There were significantly greater decreases with silodosin than naftopidil at 1 month in the total IPSS. The mean changes in the PVR at 6 months were +14.6, +23.7, and +5.7 mL in the naftopidil, tamsulosin, and silodosin groups, respectively; silodosin showed a significant improvement in the PVR at 6 months vs. tamsulosin. With respect to the secondary efficacy variables, the mean changes in the IPSS voiding score at 1 month in the naftopidil, tamsulosin, and silodosin groups were +6.5, +5.6, and +4.5, respectively; silodosin showed a significant improvement in the IPSS voiding score at 1 month vs. naftopidil. Conclusions: Silodosin has a greater impact on improving PI-induced lower urinary tract symptoms than the other two agents.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Multi-element flow-driven spectral chaos (ME-FSC) method for uncertainty quantification of dynamical systems

The flow-driven spectral chaos (FSC) is a recently developed method for tracking and quantifying uncertainties in the long-time response of stochastic dynamical systems using the spectral approach. The method uses a novel concept called enriched stochastic flow maps as a means to construct an evolving finite-dimensional random function space that is both accurate and computationally efficient in time. In this paper, we present a multi-element version of the FSC method (the ME-FSC method for short) to tackle (mainly) those dynamical systems that are inherently discontinuous over the probability space. In ME-FSC, the random domain is partitioned into several elements, and then the problem is solved separately on each random element using the FSC method. Subsequently, results are aggregated to compute the probability moments of interest using the law of total probability. To demonstrate the effectiveness of the ME-FSC method in dealing with discontinuities and long-time integration of stochastic dynamical systems, four representative numerical examples are presented in this paper, including the Van-der-Pol oscillator problem and the Kraichnan-Orszag three-mode problem. Results show that the ME-FSC method is capable of solving problems that have strong nonlinear dependencies over the probability space, both reliably and at low computational cost.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Scattering of electromagnetic waves from random media with strong permittivity fluctuations

By taking into account the singularity of the dyadic Green's function in the renormalization method, a theory is derived for vector electromagnetic wave propagation in a random medium with large permittivity fluctuations and with anisotropic correlation function. The strong fluctuation theory is then applied to a discrete scatterer problem in which the permittivity can assume only two values. The results are found to be consistent with those derived from discrete scatterer theory for all values of dielectric constants of the scatterers.

Tsang, L.↗
DOE OSTI2011

A dynamically adaptive wavelet approach to stochastic computations based on polynomial chaos - capturing all scales of random modes on independent grids

Highlights: {yields} New approach for stochastic computations based on polynomial chaos. {yields} Development of dynamically adaptive wavelet multiscale solver using space refinement. {yields} Accurate capture of steep gradients and multiscale features in stochastic problems. {yields} All scales of each random mode are captured on independent grids. {yields} Numerical examples demonstrate the need for different space resolutions per mode. - Abstract: In stochastic computations, or uncertainty quantification methods, the spectral approach based on the polynomial chaos expansion in random space leads to a coupled system of deterministic equations for the coefficients of the expansion. The size of this system increases drastically when the number of independent random variables and/or order of polynomial chaos expansions increases. This is invariably the case for large scale simulations and/or problems involving steep gradients and other multiscale features; such features are variously reflected on each solution component or random/uncertainty mode requiring the development of adaptive methods for their accurate resolution. In this paper we propose a new approach for treating such problems based on a dynamically adaptive wavelet methodology involving space-refinement on physical space that allows all scales of each solution component to be refined independently of the rest. We exemplify this using the convection-diffusion model with random input data and present three numerical examples demonstrating the salient features of the proposed method. Thus we establish a new, elegant and flexible approach for stochastic problems with steep gradients and multiscale features based on polynomial chaos expansions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Brownian dynamics of confined rigid bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Metropolis–Hastings Monte Carlo Method for Neutron Emissivity Tomographic Inversion in Tokamak Plasma

The article presents a new approach to plasma neutron emissivity reconstruction based on the Metropolis–Hastings Monte Carlo algorithm. The algorithm is based on a biased random walk. A dedicated computer code generates pseudo-random samples within the domain that contains the solution. The properly chosen objective function ensures convergence to the desired solution. The method allows to overcome the problem of trapping in local minima by accepting, with some finite probability, changes in the direction opposite to that indicated by the objective function. The developed method was tested using a synthetic data set. The results show that the method provides reliable reconstructions and thus can be used as an alternative or complementary method to the commonly applied ones. The method capabilities were also tested against noise present in the line-integrated projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗