Abstract

The necessity of dealing with uncertainties is growing in many different fields of science and engineering. Due to the constant development of computational capabilities, current solvers must satisfy both statistical accuracy and computational efficiency. The aim of this work is to introduce an asynchronous framework for Monte Carlo and Multilevel Monte Carlo methods to achieve such a result. The proposed approach presents the same reliability of state of the art techniques, and aims at improving the computational efficiency by adding a new level of parallelism with respect to existing algorithms: between batches, where each batch owns its hierarchy and is independent from the others. Two different numerical problems are considered and solved in a supercomputer to show the behavior of the proposed approach.

Highlights

  • The increasing necessity of handling data uncertainties in many different science and engineering scenarios and the development of computational capabilities of supercomputers are leading to the necessity of having a strict integration between Uncertainty Quantification (UQ) techniques and efficient algorithms

  • We propose an asynchronous version of the algorithms belonging to the Monte Carlo (MC) family, together with a reference implementation based on the Kratos Multiphysics software [13,14] and the XMC library [3], exploiting the PyCOMPSs programming model [2,25,34]

  • In the tables presenting the results, B identifies the initial number of batches, defined by the initial hierarchy, it the number of iterations, that is the amount of convergence checks executed, N the total number of realizations computed per level at the end of the of the execution, EMC[Q H ] and EMLMC[Q H ] the expected value estimate, SE the statistical error, time to solution the total time the simulation required to finish, measured in seconds, and C the computational cost of the algorithm run, expressed in CPU hours

Read more

Summary

Introduction

The increasing necessity of handling data uncertainties in many different science and engineering scenarios and the development of computational capabilities of supercomputers are leading to the necessity of having a strict integration between Uncertainty Quantification (UQ) techniques and efficient algorithms. UQ studies how uncertainties propagate in the. The authors thankfully acknowledge the computer resources at MareNostrum and the technical support provided by Barcelona Supercomputing Center (IM-2020-1-0016)

28 Page 2 of 25
Compressible Potential Flow Problem
Wind Engineering Problem
Statistical Overview
28 Page 6 of 25
Monte Carlo
Multilevel Monte Carlo
Adaptivity
28 Page 8 of 25
Computation of Moments
Online Update of Central Moments
Online Update of Power Sums
Stopping Criteria
Algorithms
Asynchronous Framework
Asynchronous Algorithms
Numerical Experiments
Compressible Potential Flow Benchmark
MC Analysis
MLMC Analysis
Traces
Wind Engineering Benchmark
28 Page 18 of 25
28 Page 20 of 25
Scalability of Asynchronous Algorithm
Conclusions
28 Page 24 of 25
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call