Herbert Owen

Alya: Multiphysics engineering simulation toward exascale

Alya is a multi-physics simulation code developed at Barcelona Supercomputing Center (BSC). From its inception Alya code is designed using advanced High Performance Computing programming techniques to solve coupled problems on supercomputers efficiently. The target domain is engineering, with all its particular features: complex geometries and unstructured meshes, coupled multi-physics with exotic coupling schemes and physical models, ill-posed problems, flexibility needs for rapidly including new models, etc. Since its beginnings in 2004, Alya has scaled well in an increasing number of processors when solving single-physics problems such as fluid mechanics, solid mechanics, acoustics, etc. Over time, we have made a concerted effort to maintain and even improve scalability for multi-physics problems. This poses challenges on multiple fronts, including: numerical models, parallel implementation, physical coupling models, algorithms and solution schemes, meshing process, etc. In this paper, we introduce Alyas main features and focus particularly on its solvers. We present Alyas performance up to 100.000 processors in Blue Waters, the NCSA supercomputer with selected multi-physics tests that are representative of the engineering world. The tests are incompressible flow in a human respiratory system, low Mach combustion problem in a kiln furnace, and coupled electro-mechanical contraction of the heart. We show scalability plots for all cases and discuss all aspects of such simulations, including solver convergence.

Mesh generation, sizing and convergence for onshore and offshore wind farm Atmospheric Boundary Layer flow simulation with actuator discs

Abstract A new mesh generation process for wind farm modeling is presented together with a mesh convergence and sizing analysis for wind farm flow simulations. The generated meshes are tailored to simulate Atmospheric Boundary Layer (ABL) flows on complex terrains modeling the wind turbines as actuator discs. The wind farm mesher is fully automatic and, given the topography and the turbine characteristics (location, diameter and hub height), it generates a hybrid mesh conformal with the actuator discs and refined upwind and downstream. Moreover, it presents smooth element size transitions across scales and avoids extending high-resolution areas to all the domain. We take advantage of our automatic and robust mesher to study the mesh convergence of our RANS solver with linear elements, obtaining quadratic mesh convergence for a quantity of interest in all the tested cases. In addition, we quantify the mesh resolution at the terrain surface and at the actuator discs required to achieve a given numerical error in simulations in onshore and offshore frameworks. Finally, we present the generated meshes and the simulation results for an offshore and an onshore wind farm. We analyze in detail one particular wind direction for both cases, and for the onshore wind farm we use our automatic framework to estimate the yearly production of energy and measuring the error against the actual produced one.