We propose a hybrid algorithm for the time integration of large sets of rate equations coupled by a relatively small number of degrees of freedom. A subset containing fast degrees of freedom evolves deterministically, while the rest of the variables evolves stochastically. The emphasis is put on the coupling between the two subsets, in order to achieve both accuracy and efficiency. The algorithm is tested on the problem of nucleation, growth and coarsening of clusters of defects in iron, treated by the formalism of cluster dynamics. We show that it is possible to obtain results indistinguishable from fully deterministic and fully stochastic calculations, while speeding up significantly the computations with respect to these two cases.
Hybrid deterministic/stochastic algorithm for large sets of rate equations / M. Gherardi, T. Jourdan, S. Le Bourdiec, G. Bencteux. - In: COMPUTER PHYSICS COMMUNICATIONS. - ISSN 0010-4655. - 183:9(2012), pp. 1966-1973.
|Titolo:||Hybrid deterministic/stochastic algorithm for large sets of rate equations|
|Parole Chiave:||Cluster dynamics; Helium; Iron; Irradiation; Monte Carlo; Rate theory; Stiffness; Hardware and Architecture; Physics and Astronomy (all)|
|Settore Scientifico Disciplinare:||Settore FIS/02 - Fisica Teorica, Modelli e Metodi Matematici|
Settore FIS/03 - Fisica della Materia
|Data di pubblicazione:||2012|
|Digital Object Identifier (DOI):||http://dx.doi.org/10.1016/j.cpc.2012.04.020|
|Appare nelle tipologie:||01 - Articolo su periodico|