A new proof of geometric convergence for general transport problems based on sequential correlated sampling methods
نویسندگان
چکیده
In [J. Halton, Sequential Monte Carlo, Proc. Comb. Phil. Soc. 58 (1962), J. Halton, Sequential Monte Carlo Techniques for the Solution of Linear Systems, J. Sci. Comp. 9 (1994) 213-257] Halton introduced a strategy to be used in Monte Carlo algorithms for the efficient solution of certain matrix problems. We showed in [R. Kong, J. Spanier, Sequential correlated sampling methods for some transport problems, in: Harold Niederreiter, Jerome Spanier (Eds.), Monte-Carlo and Quasi Monte-Carlo Methods 1998: Proceedings of a Conference at the Claremont Graduate University, Springer-Verlag, New York, 2000, R. Kong, J. Spanier, Error analysis of sequential Monte Carlo methods for transport problems, in: Harold Niederreiter, Jerome Spanier (Eds.), Monte-Carlo and Quasi Monte-Carlo Methods 1998: Proceedings of a Conference at the Claremont Graduate University, Springer-Verlag, New York, 2000] how Halton's method based on correlated sampling can be extended to continuous transport problems and established their geometric convergence for a family of transport problems in slab geometry. In our algorithm, random walks are processed in batches, called stages, each stage producing a small correction that is added to the approximate solution developed from the previous stages. In this paper, we demonstrate that strict error reduction from stage to stage can be assured under rather general conditions and we illustrate this rapid convergence numerically for a simple family of two dimensional transport problems.
منابع مشابه
Comparison of Monte Carlo Algorithms for Obtaining Geometric Convergence for Model Transport Problems
Two quite diierent methods for accelerating the convergence of global Monte Carlo solutions of continuous transport problems have been developed recently in Claremont. One of these is based on a sequential form of correlated sampling , rst proposed for matrix problems by Haltonn1]. The second method makes use of importance sampling transformations applied adaptivelyy2]. These two methods are co...
متن کاملGeometrically Convergent Learning Algorithms for Global Solutions of Transport Problems
In 1996 Los Alamos National Laboratory initiated an ambitious ve year research program aimed at achieving geometric convergence for Monte Carlo solutions of diicult neutron and photon transport problems. Claremont students, working with the author in Mathematics Clinic projects that same year and subsequently , have been partners in this undertaking. This paper summarizes progress made at Clare...
متن کاملAdaptive Monte Carlo Algorithms for General Transport Problems
Recently there has been a concerted effort to develop adaptively modified Monte Carlo algorithms that converge geometrically to solutions of the radiative transport equation. We have concentrated on algorithms that extend to integral equations methods first proposed for matrix equations by Halton in 1962 [Halton, J., Proc. Camb. Phil. Soc., 58, 57–78 (1962)]. Geometric convergence has been rigo...
متن کاملA TRUST-REGION SEQUENTIAL QUADRATIC PROGRAMMING WITH NEW SIMPLE FILTER AS AN EFFICIENT AND ROBUST FIRST-ORDER RELIABILITY METHOD
The real-world applications addressing the nonlinear functions of multiple variables could be implicitly assessed through structural reliability analysis. This study establishes an efficient algorithm for resolving highly nonlinear structural reliability problems. To this end, first a numerical nonlinear optimization algorithm with a new simple filter is defined to locate and estimate the most ...
متن کاملSequential Implicit Numerical Scheme for Pollutant and Heat Transport in a Plane-Poiseuille Flow
A sequential implicit numerical scheme is proposed for a system of partial differential equations defining the transport of heat and mass in the channel flow of a variable-viscosity fluid. By adopting the backward difference scheme for time derivative and the central difference scheme for the spatial derivatives, an implicit finite difference scheme is formulated. The variable-coefficient diffu...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Journal of computational physics
دوره 227 23 شماره
صفحات -
تاریخ انتشار 2008