Analysis of Relative Error in Perturbation Monte Carlo Simulations of Radiative Transport

SIGNIFICANCE: Perturbation and differential Monte Carlo (pMC/dMC) methods, used in conjunction with nonlinear optimization methods, have been successfully applied to solve inverse problems in diffuse optics. Application of pMC to systems over a large range of optical properties requires optimal placement of baseline conventional Monte Carlo (cMC) simulations to minimize the pMC variance. The inability to predict the growth in pMC solution uncertainty with perturbation size limits the application of pMC, especially for multispectral datasets where the variation of optical properties can be substantial. AIM: We aim to predict the variation of pMC variance with perturbation size without explicit computation of perturbed photon weights. Our proposed method can be used to determine the range of optical properties over which pMC predictions provide sufficient accuracy. This method can be used to specify the optical properties for the reference cMC simulations that pMC utilizes to provide accurate predictions over a desired optical property range. APPROACH: We utilize a conventional error propagation methodology to calculate changes in pMC relative error for Monte Carlo simulations. We demonstrate this methodology for spatially resolved diffuse reflectance measurements with ±20% scattering perturbations. We examine the performance of our method for reference simulations spanning a broad range of optical properties relevant for diffuse optical imaging of biological tissues. Our predictions are computed using the variance, covariance, and skewness of the photon weight, path length, and collision distributions generated by the reference simulation. RESULTS: We find that our methodology performs best when used in conjunction with reference cMC simulations that utilize Russian Roulette (RR) method. Specifically, we demonstrate that for a proximal detector placed immediately adjacent to the source, we can estimate the pMC relative error within 5% of the true value for scattering perturbations in the range of CONCLUSIONS: These findings indicate that reference simulations utilizing continuous absorption weighting (CAW) with the Russian Roulette method and executed using optical properties with a lo

Discutindo a educação ambiental no cotidiano escolar: desenvolvimento de projetos na escola formação inicial e continuada de professores

A presente pesquisa buscou discutir como a Educação Ambiental (EA) vem sendo trabalhada, no Ensino Fundamental e como os docentes desta escola compreendem e vem inserindo a EA no cotidiano escolar., em uma escola estadual do município de Tangará da Serra/MT, Brasil. Para tanto, realizou-se entrevistas com os professores que fazem parte de um projeto interdisciplinar de EA na escola pesquisada. Verificou-se que o projeto da escola não vem conseguindo alcançar os objetivos propostos por: desconhecimento do mesmo, pelos professores; formação deficiente dos professores, não entendimento da EA como processo de ensino-aprendizagem, falta de recursos didáticos, planejamento inadequado das atividades. A partir dessa constatação, procurou-se debater a impossibilidade de tratar do tema fora do trabalho interdisciplinar, bem como, e principalmente, a importância de um estudo mais aprofundado de EA, vinculando teoria e prática, tanto na formação docente, como em projetos escolares, a fim de fugir do tradicional vínculo “EA e ecologia, lixo e horta”.Facultad de Humanidades y Ciencias de la Educació

A new proof of geometric convergence for the adaptive generalized weighted analog sampling (GWAS) method.

Generalized Weighted Analog Sampling is a variance-reducing method for solving radiative transport problems that makes use of a biased (though asymptotically unbiased) estimator. The introduction of bias provides a mechanism for combining the best features of unbiased estimators while avoiding their limitations. In this paper we present a new proof that adaptive GWAS estimation based on combining the variance-reducing power of importance sampling with the sampling simplicity of correlated sampling yields geometrically convergent estimates of radiative transport solutions. The new proof establishes a stronger and more general theory of geometric convergence for GWAS

Hybrid Monte Carlo estimators for multilayer transport problems.

This paper introduces a new family of hybrid estimators aimed at controlling the efficiency of Monte Carlo computations in particle transport problems. In this context, efficiency is usually measured by the figure of merit (FOM) given by the inverse product of the estimator variance Var[ξ] and the run time T: FOM := (Var[ξ] T)-1. Previously, we developed a new family of transport-constrained unbiased radiance estimators (T-CURE) that generalize the conventional collision and track length estimators [1] and provide 1-2 orders of magnitude additional variance reduction. However, these gains in variance reduction are partly offset by increases in overhead time [2], lowering their computational efficiency. Here we show that combining T-CURE estimation with conventional terminal estimation within each individual biography can moderate the efficiency of the resulting "hybrid" estimator without introducing bias in the computation. This is achieved by treating only the refractive interface crossings with the extended next event estimator, and all others by standard terminal estimators. This is because when there are index-mismatched interfaces between the collision location and the detector, the T-CURE computation rapidly becomes intractable due to the large number of refractions and reflections that can arise. We illustrate the gains in efficiency by comparing our hybrid strategy with more conventional estimation methods in a series of multi-layer numerical examples