This title appears in the Scientific Report :
2019
Please use the identifier:
http://dx.doi.org/10.1016/j.jcp.2018.07.051 in citations.
Please use the identifier: http://hdl.handle.net/2128/23398 in citations.
On the possibility of track length based Monte-Carlo algorithms for stationary drift-diffusion systems with sources and sinks
On the possibility of track length based Monte-Carlo algorithms for stationary drift-diffusion systems with sources and sinks
The problem of constructing Monte-Carlo solutions of drift-diffusion systems corresponding to Fokker–Planck equations with sources and sinks is revisited. Firstly, a compact formalism is introduced for the specific problem of stationary solutions. This leads to identification of the dwell time as th...
Saved in:
Personal Name(s): | Reiser, D. (Corresponding author) |
---|---|
Romazanov, J. / Linsmeier, Ch. | |
Contributing Institute: |
Plasmaphysik; IEK-4 |
Published in: | Journal of computational physics, 377 (2019) S. 219 - 231 |
Imprint: |
Amsterdam
Elsevier
2019
|
DOI: |
10.1016/j.jcp.2018.07.051 |
Document Type: |
Journal Article |
Research Program: |
Plasma-Wall-Interaction |
Link: |
Published on 2018-07-30. Available in OpenAccess from 2020-07-30. Published on 2018-07-30. Available in OpenAccess from 2020-07-30. |
Publikationsportal JuSER |
Please use the identifier: http://hdl.handle.net/2128/23398 in citations.
The problem of constructing Monte-Carlo solutions of drift-diffusion systems corresponding to Fokker–Planck equations with sources and sinks is revisited. Firstly, a compact formalism is introduced for the specific problem of stationary solutions. This leads to identification of the dwell time as the key quantity to characterize the system and to obtain a proper normalization for statistical analysis of numerical results. Secondly, the question of appropriate track length estimators for drift-diffusion systems is discussed for a 1D model system. It is found that a simple track length estimator can be given only for pure drift motion without diffusion. The stochastic nature of the diffusive part cannot be appropriately described by the path length of simulation particles. Further analysis of the usual situation with inhomogeneous drift and diffusion coefficients leads to an error estimate based on particle trajectories. The result for limits in grid cell size and time step used for the construction of Monte-Carlo trajectories resembles the Courant-Friedrichs-Lewy and von Neumann conditions for explicit methods. |