• Title/Summary/Keyword: Monte-Carlo(MC)

Search Result 216, Processing Time 0.028 seconds

Fast Noise Reduction Approach in Multifocal Multiphoton Microscopy Based on Monte-Carlo Simulation

  • Kim, Dongmok;Shin, Younghoon;Kwon, Hyuk-Sang
    • Current Optics and Photonics
    • /
    • v.5 no.4
    • /
    • pp.421-430
    • /
    • 2021
  • The multifocal multiphoton microscopy (MMM) enables high-speed imaging by the concurrent scanning and detection of multiple foci generated by lenslet array or diffractive optical element. The MMM system mainly suffers from crosstalk generated by scattered emission photons that form ghost images among adjacent channels. The ghost image which is a duplicate of the image acquired in sub-images significantly degrades overall image quality. To eliminate the ghost image, the photon reassignment method was established using maximum likelihood estimation. However, this post-processing method generally takes a longer time than image acquisition. In this regard, we propose a novel strategy for rapid noise reduction in the MMM system based upon Monte-Carlo (MC) simulation. Ballistic signal, scattering signal, and scattering noise of each channel are quantified in terms of photon distribution launched in tissue model based on MC simulation. From the analysis of photon distribution, we successfully eliminated the ghost images in the MMM sub-images. If the priori MC simulation under a certain optical condition is established at once, our simple, but robust post-processing technique will continuously provide the noise-reduced images, while significantly reducing the computational cost.

An Introduction to Kinetic Monte Carlo Methods for Nano-scale Diffusion Process Modeling (나노 스케일 확산 공정 모사를 위한 동력학적 몬테칼로 소개)

  • Hwang, Chi-Ok;Seo, Ji-Hyun;Kwon, Oh-Seob;Kim, Ki-Dong;Won, Tae-Young
    • Journal of the Institute of Electronics Engineers of Korea SD
    • /
    • v.41 no.6
    • /
    • pp.25-31
    • /
    • 2004
  • In this paper, we introduce kinetic Monte Carlo (kMC) methods for simulating diffusion process in nano-scale device fabrication. At first, we review kMC theory and backgrounds and give a simple point defect diffusion process modeling in thermal annealing after ion (electron) implantation into Si crystalline substrate to help understand kinetic Monte Carlo methods. kMC is a kind of Monte Carlo but can simulate time evolution of diffusion process through Poisson probabilistic process. In kMC diffusion process, instead of. solving differential reaction-diffusion equations via conventional finite difference or element methods, it is based on a series of chemical reaction (between atoms and/or defects) or diffusion events according to event rates of all possible events. Every event has its own event rate and time evolution of semiconductor diffusion process is directly simulated. Those event rates can be derived either directly from molecular dynamics (MD) or first-principles (ab-initio) calculations, or from experimental data.

McCARD/MIG stochastic sampling calculations for nuclear cross section sensitivity and uncertainty analysis

  • Ho Jin Park
    • Nuclear Engineering and Technology
    • /
    • v.54 no.11
    • /
    • pp.4272-4279
    • /
    • 2022
  • In this study, a cross section stochastic sampling (S.S.) capability is implemented into both the McCARD continuous energy Monte Carlo code and MIG multiple-correlated data sampling code. The ENDF/B-VII.1 covariance data based 30 group cross section sets and the SCALE6 covariance data based 44 group cross section sets are sampled by the MIG code. Through various uncertainty quantification (UQ) benchmark calculations, the McCARD/MIG results are verified to be consistent with the McCARD stand-alone sensitivity/uncertainty (S/U) results and the XSUSA S.S. results. UQ analyses for Three Mile Island Unit 1, Peach Bottom Unit 2, and Kozloduy-6 fuel pin problems are conducted to provide the uncertainties of keff and microscopic and macroscopic cross sections by the McCARD/MIG code system. Moreover, the SNU S/U formulations for uncertainty propagation in a MC depletion analysis are validated through a comparison with the McCARD/MIG S.S. results for the UAM Exercise I-1b burnup benchmark. It is therefore concluded that the SNU formulation based on the S/U method has the capability to accurately estimate the uncertainty propagation in a MC depletion analysis.

Implementation and benchmarking of the local weight window generation function for OpenMC

  • Hu, Yuan;Yan, Sha;Qiu, Yuefeng
    • Nuclear Engineering and Technology
    • /
    • v.54 no.10
    • /
    • pp.3803-3810
    • /
    • 2022
  • OpenMC is a community-driven open-source Monte Carlo neutron and photon transport simulation code. The Weight Window Mesh (WWM) function and an automatic Global Variance Reduction (GVR) method was recently developed and implemented in a developmental branch of OpenMC. This WWM function and GVR method broaden OpenMC's usage in general purposes deep penetration shielding calculations. However, the Local Variance Reduction (LVR) method, which suits the source-detector problem, is still missing in OpenMC. In this work, the Weight Window Generator (WWG) function has been developed and benchmarked for the same branch. This WWG function allows OpenMC to generate the WWM for the source-detector problem on its own. Single-material cases with varying shielding and sources were used to benchmark the WWG function and investigate how to set up the particle histories utilized in WWG-run and WWM-run. Results show that there is a maximum improvement of WWM generated by WWG. Based on the above results, instructions on determining the particle histories utilized in WWG-run and WWM-run for optimal computation efficiency are given and tested with a few multi-material cases. These benchmarks demonstrate the ability of the OpenMC WWG function and the above instructions for the source-detector problem. This developmental branch will be released and merged into the main distribution in the future.

Development and verification of a Monte Carlo two-step method for lead-based fast reactor neutronics analysis

  • Yiwei Wu;Qufei Song;Ruixiang Wang;Yao Xiao;Hanyang Gu;Hui Guo
    • Nuclear Engineering and Technology
    • /
    • v.55 no.6
    • /
    • pp.2112-2124
    • /
    • 2023
  • With the rise of economic and safety standards for nuclear reactors, new concepts of Gen-IV reactors and modular reactors showed more complex designs that challenge current tools for reactor physics analysis. A Monte Carlo (MC) two-step method was proposed in this work. This calculation scheme uses the continuous-energy MC method to generate multi-group cross-sections from heterogeneous models. The multi-group MC method, which can adapt locally-heterogeneous models, is used in the core calculation step. This calculation scheme is verified using a Gen-IV modular lead-based fast reactor (LFR) benchmark case. The influence of homogenized patterns, scatter approximations, flux separable approximation, and local heterogeneity in core calculation on simulation results are investigated. Results showed that the cross-sections generated using the 3D assembly model with a locally heterogeneous representation of control rods lead to an accurate estimation with less than 270 pcm bias in core reactivity, 0.5% bias in control rod worth, and 1.5% bias on power distribution. The study verified the applicability of multi-group cross-sections generated with the MC method for LFR analysis. The study also proved the feasibility of multi-group MC in core calculation with local heterogeneity, which saves 85% time compared to the continuous-energy MC.

Efficient Performance Evaluation Method for IS-95 System (IS-95 시스템 역방향 채널에서의 효율적인 성능평가 기법)

  • 전재춘;고윤진;정미선;황인관
    • The Journal of Korean Institute of Communications and Information Sciences
    • /
    • v.27 no.4B
    • /
    • pp.345-352
    • /
    • 2002
  • In this paper, in order to evaluate the performance of IS-95 system reverse link in white gaussian noise and rayleigh fading environment, we suggest epochal proposal to improve computer run-time and its efficiency is verified in terms of the number of samples. MC(Monte Carlo) simulation is the most popular simulation technique lately, but MC simulation requires a number of samples at low bit error rate. Therefore, MC cannot avoid the limit of computer run-time. To alleviate these problems, we apply the suggested method called central moment technique to the reverse link of the IS-95 system and can obtain discrete probability mass functions from Nth order central moments of the less number of received signal samples than those required in MC. Continuous cumulative probability distribution function can be accurately estimated by using interpolation and the improvement effect for the number of samples is proven.

A hybrid neutronics method with novel fission diffusion synthetic acceleration for criticality calculations

  • Jiahao Chen;Jason Hou;Kostadin Ivanov
    • Nuclear Engineering and Technology
    • /
    • v.55 no.4
    • /
    • pp.1428-1438
    • /
    • 2023
  • A novel Fission Diffusion Synthetic Acceleration (FDSA) method is developed and implemented as a part of a hybrid neutronics method for source convergence acceleration and variance reduction in Monte Carlo (MC) criticality calculations. The acceleration of the MC calculation stems from constructing a synthetic operator and solving a low-order problem using information obtained from previous MC calculations. By applying the P1 approximation, two correction terms, one for the scalar flux and the other for the current, can be solved in the low-order problem and applied to the transport solution. A variety of one-dimensional (1-D) and two-dimensional (2-D) numerical tests are constructed to demonstrate the performance of FDSA in comparison with the standalone MC method and the coupled MC and Coarse Mesh Finite Difference (MC-CMFD) method on both intended purposes. The comparison results show that the acceleration by a factor of 3-10 can be expected for source convergence and the reduction in MC variance is comparable to CMFD in both slab and full core geometries, although the effectiveness of such hybrid methods is limited to systems with small dominance ratios.

A comparison study of approximate and Monte Carlo radiative transfer methods for late type galaxy models

  • Lee, Dukhang;Baes, Maarten;Seon, Kwang-il
    • The Bulletin of The Korean Astronomical Society
    • /
    • v.41 no.1
    • /
    • pp.49.3-50
    • /
    • 2016
  • Two major radiative transfer (RT) techniques have been developted to model late-type galaxies: approximate RT and Monte Carlo (MC) RT. In the approximate RT, first proposed by Kylafis & Bahcall, only two terms of unscattered (direct) and single-scattered intensities are computed and higher-order multiple scattering components are approximated, saving computing time and cost compared to MC RT. However, the approximate RT can yield errors in regions where multiple scattering effect is significant. In order to examine how significant the errors of the approximate RT are, we compare results of the approximate RT with those of SKIRT, a state-of-the-art MC RT code, which is basically free from the approximation errors by fully incorporating all the multiple scattered intensities. In this study, we present quantitative errors in the approximate RT for late type galaxy models with various optical depths and inclination angles. We report that the approximate RT is not reliable if the central face-on optical depth is intermediate or high (${\tau}_V$ > 3).

  • PDF

Random Balance between Monte Carlo and Temporal Difference in off-policy Reinforcement Learning for Less Sample-Complexity (오프 폴리시 강화학습에서 몬테 칼로와 시간차 학습의 균형을 사용한 적은 샘플 복잡도)

  • Kim, Chayoung;Park, Seohee;Lee, Woosik
    • Journal of Internet Computing and Services
    • /
    • v.21 no.5
    • /
    • pp.1-7
    • /
    • 2020
  • Deep neural networks(DNN), which are used as approximation functions in reinforcement learning (RN), theoretically can be attributed to realistic results. In empirical benchmark works, time difference learning (TD) shows better results than Monte-Carlo learning (MC). However, among some previous works show that MC is better than TD when the reward is very rare or delayed. Also, another recent research shows when the information observed by the agent from the environment is partial on complex control works, it indicates that the MC prediction is superior to the TD-based methods. Most of these environments can be regarded as 5-step Q-learning or 20-step Q-learning, where the experiment continues without long roll-outs for alleviating reduce performance degradation. In other words, for networks with a noise, a representative network that is regardless of the controlled roll-outs, it is better to learn MC, which is robust to noisy rewards than TD, or almost identical to MC. These studies provide a break with that TD is better than MC. These recent research results show that the way combining MC and TD is better than the theoretical one. Therefore, in this study, based on the results shown in previous studies, we attempt to exploit a random balance with a mixture of TD and MC in RL without any complicated formulas by rewards used in those studies do. Compared to the DQN using the MC and TD random mixture and the well-known DQN using only the TD-based learning, we demonstrate that a well-performed TD learning are also granted special favor of the mixture of TD and MC through an experiments in OpenAI Gym.

Development of Electron Beam Monte Carlo Simulation and Analysis of SEM Imaging Characteristics (전자빔 몬테 카를로 시물레이션 프로그램 개발 및 전자현미경 이미징 특성 분석)

  • Kim, Heung-Bae
    • Journal of the Korean Society for Precision Engineering
    • /
    • v.29 no.5
    • /
    • pp.554-562
    • /
    • 2012
  • Processing of Scanning electron microscope imaging has been analyzed in both secondary electron (SE) imaging and backscattered electron (BSE) image. Because of unique characteristics of both secondary electron and backscattered electron image, mechanism of imaging process and image quality are quite different each other. For the sake of characterize imaging process, Monte Carlo simulation code have been developed. It simulates electron penetration and depth profile in certain material. In addition, secondary electron and backscattered electron generation process as well as their spatial distribution and energy characteristics can be simulated. Geometries that has fundamental feature have been imaged using the developed Monte Carlo code. Two, SE and BSE images generation process will be discussed. BSE imaging process can be readily used to discriminate in both material and geometry by simply changing position and direction of BSE detector. The developed MC code could be useful to design BSE detector and their position. Furthermore, surface reconstruction technique is possibly developed at the further research efforts. Basics of Monte Carlo simulation method will be discussed as well as characteristics of SE and BSE images.