Author: Riccardo Tosi
Kratos version: 8.1
XMC version: 2.0
PyCOMPSs version: 2.7
Source files: Asynchronous and Synchronous Monte Carlo
Application dependencies: ConvectionDiffusionApplication, ExaquteSandboxApplication, FluidDynamicsApplication, LinearSolversApplications, MappingApplication, MeshingApplication, MultilevelMonteCarloApplication, StatisticsApplication
We solve the fluid dynamics problem of a fluid passing through a bluff body. The problem is characterized by uniform and stochastic wind inlet velocity. The wind inlet velocity follows the following probability density function .
As seen in [Ayoul-Guilmard, Q., Núñez, M., Ganesh, S., Nobile, F., Rossi, R., & Tosi, R. (2020). D5.3 Report on theoretical work to allow the use of MLMC with adaptive mesh refinement.], we could not prove Multilevel Monte Carlo hypotheses for such turbulent and chaotic problem. For this reason, we apply Monte Carlo, exploiting XMC.
The problem can be run with two different algorithms:
- Synchronous Monte Carlo (SMC),
- Asynchronous Monte Carlo (AMC),
and by default AMC is selected. If one is interested in running SMC, it is needed to select asynchronous = false in the solver wrapper settings.
The Quantities of Interest of the problem are the drag force, the base moment and the pressure field. Statistical convergence is assessed for the drag force.
All settings can be observed in the corresponding configuration file of the algorithm, located inside the problem_settings folder.
To run the examples, the user should go inside the source folder and run the run_mc_Kratos.py Python file. In case one wants to use PyCOMPSs, the user should execute run.sh from inside the source folder.
The velocity and pressure fields evolution of the problem are shown next.
The power sums and the h-statistics of both the time averaged and time series drag force, base moment and pressure field can be found here.
In addition,the drag coefficient we estimate from the drag force is consistent with literature [Bruno, L., Salvetti, M. V., & Ricciardelli, F. (2014). Benchmark on the aerodynamics of a rectangular 5:1 cylinder: An overview after the first four years of activity. Journal of Wind Engineering and Industrial Aerodynamics, 126, 87–106. https://doi.org/10.1016/j.jweia.2014.01.005].