(315e) Multilevel Monte Carlo Applied for Efficient Estimation of Observables in Multiscale Stochastic Systems
AIChE Annual Meeting
2018
2018 AIChE Annual Meeting
Computing and Systems Technology Division
Advances in Computational Methods and Numerical Analysis
Tuesday, October 30, 2018 - 1:46pm to 2:05pm
In our previous work, we applied MLMC to uncertainty quantification in classic chemical engineering systems [3]. In this work, we extend and adapt MLMC to perform accurate estimation of observables in stochastic multiscale systems. We used as a case study a model of thin film formation by Chemical Vapour Deposition (CVD) [4], which was also previously used to quantify uncertainty and estimation of observables in multiscale systems [5],[6]. In this model, a kinetic Monte Carlo (kMC) solid-on-solid simulation represents the microscale formation of a thin film, while continuum mass, energy and momentum transfer equations are used to represent the macroscale supply of the chemical vapor to the substrate where the film forms. The non-closed-form microscale representation and the deterministic macroscale equations are coupled through a boundary condition and thus form a stochastic multiscale system.
Most of the studies that use MLMC, including our previous work [3], discretize the time domain of continuous equations and use finer discretization for higher levels of approximation accuracy. In the CVD system, the accuracy of approximation increases with larger lattice sizes of the kMC simulation. We established an empirical relationship that relates the kMC time step to the lattice size and used it to control the definition of higher accuracy levels in the MLMC framework. We demonstrate that MLMC can be used to accurately estimate the observables of this multiscale stochastic system in a fraction of the time required when conducting the conventional Monte Carlo sampling of the full multiscale simulation.
References
[1] Eldred, M.S., âRecent Advances in Non-Intrusive Polynomial Chaos and Stochastic Collocation Methods for Uncertainty Analysis and Designâ, In Proceedings of the 50th AIAA Structures, Structural Dynamics and Materials Conference, Palm Springs, California, 2009.
[2] Giles, M. âMultilevel Monte Carlo Path Simulationâ, Operations Research, 56, 3, 2008.
[3] Kimaev, G. and Ricardez-Sandoval, L.A., âMultilevel Monte Carlo Applied to Chemical Engineering Systems Subject to Uncertaintyâ, AIChE Journal, 64, 5, 2018.
[4] Vlachos, D., âMultiscale Integration Hybrid Algorithms for HomogeneousâHeterogeneous Reactorsâ, AIChE Journal, 43, 11, 1997.
[5] Rasoulian, S. and Ricardez-Sandoval, L.A. âUncertainty Analysis and Robust Optimization of Multiscale Process Systems With Application to Epitaxial Thin Film Growthâ, Chemical Engineering Science, 116, 590, 2014.
[6] Rasoulian, S. and Ricardez-Sandoval, L.A. âRobust Multivariable Estimation and Control in an Epitaxial Thin Film Growth Process Under Uncertaintyâ, Journal of Process Control, 34, 70, 2015.