Ecfd:ecfd 3rd edition

ECFD workshop, 3rd edition, 2020

Participants

Flyer

Flyer

Presentations

Booklet

Template

Project achievements

Project #1: Hackathon GENCI/ATOS/AMD/CERFACS on AVBP

C. Piechurski (GENCI), S. Jauré (ATOS), B. Pajot (ATOS), P.-A. Harraud (AMD), P. Mohanamuraly (CERFACS), G. Staffelbach (CERFACS), J. Legaux (CERFACS)

We ported the AVBP solver to the AMD Rome system available at GENCI -TGCC ( IRENE Joliot Curie). Characterisation of the application on the architecture showed a 1/3 performance dependency to bandwidth and 2/3 to compute. Strong scaling performance up to 130k cores was measured with openmpi and provided an acceleration of 75% without optimisations. Weak scaling up to 32k MPI ranks suggests that decimation of the processes by a factor 2 improves computational efficiency by up to 30%. This suggests a trade off between mpi imbalance and decimation is possible if imbalance is higher than 30% to improve time to solution.

Currently Openmpi offers the best perfofrmance, intelmpi is still a bit unstable.

During the Hackathon we also introduced colour based cache blocking using ColPack in the code in order to use OpenMP without critical sections. On a 2x18 core Skylake processor the new implementation offered similar speedup using full threading versus full MPI with the best trade off being 4 MPI and 9 threads per MPI. On AMD Rome, Full threading did not offer much acceleration and needs to be inversigated but 8 MPI and 16 threads per MPI seem quite promising.

Final presentation of project #1

Project #2: Hackathon GENCI/ATOS/AMD/CORIA on YALES2

C. Piechurski (GENCI), S. Jauré (ATOS), P.-A. Harraud (AMD), P. Mohanamuraly (CERFACS), G.Lartigue (CORIA), F. Gava (CORIA), P. Begou (LEGI)

Final presentation of project #2

Project #3: Implementation of a secondary atomization model in YALES2

C. G. Guillamon (Safran Tech), L .Voivenel (Safran Tech), R. Mercier (Safran Tech)

In Lagrangian simulations, droplets are transported following a ballistic motion in an eulerian mesh. For non-reactive environments, droplets might undergo secondary atomization due to the aerodynamic interaction. In this work, we implement in YALES2 a breakup model known as Taylor-Analogy Breakup (TAB). This model is based on the analogy between a droplet and a second-order mechanical system, hence making possible to determine the breakup behaviour by means of Newton's second law.

Another model, the stochastic breakup model by Gorokhovski, is also suggested for future work and will be implemented in YALES2.

Final presentation of project #3

Project #4: Application to combustion and lubrication applications

Final presentation of project #4

Project #5: Jet-in-crossflow par une méthode d’interface diffuse

Final presentation of project #5

Project #6: Accurate numerical predicti􏴇on of vorti􏴇cal flows using AMR

Final presentation of project #6 We try to demonstrate that Eulerian method YALES2 using AMR can do a very good job to capture complex vortical flows at moderate Re=10k Here we use an AMR strategy based on vorticity. We investigate the problem of vortex ring collision. We have a gain of 1000 on the numbers of elements compared to a non adaptative approach. We are able to capture the transition from a very simple laminar flow to a complex turbulent flow.

Project #7: Modélisation de parois pour la simulation des grandes échelles

Final presentation of project #7

Project #8: Accurate numerical simulation of contact lines with dynamic mesh adaptation

S. Pertant (LEGI), G. Ghigliotti (LEGI), G. Balarac (LEGI)

The main objective of this project was to develop a methodology to simulate contact lines on unstructured meshes. We especially wanted to get rid of mesh influence on contact line movement when the flow is driven by surface tension and the contact line close to its equilibrium position. A slight modification in the Ghost Fluid Method to apply the pressure jump has been tested and seems promising. The pressure gradient at contact line is indeed less sensitive to mesh elements for high density ratios. Furthermore, dynamic mesh adaptation has been used to simulate a 2D vapour bubble lying on a wall. Due to gravity, the two contact lines are receding until their merging and the bubble departure. The mesh remains fine to capture the contact line dynamics. As a future work, we plan to perform mesh adaptation on 3D contact line cases and to include additional physics such as contact angle imposition (already implemented but not used yet with mesh adaptation).

[media:ecfd3_final_project8.pdf | Final presentation of project #8]]

Project #9: Remeshed particle method at high Schmidt and Reynolds number

S. Santoso (LJK), J.-B. Lagaert (Math Orsay), G.Balarac (LEGI)

We study the advection of a scalar function in turbulent flows with a multimesh method. The finite volume method is used to solve Navier-Stokes equations on an unstructured mesh (YALES2). The advection equation is solved with remeshed particle method on a cartesian mesh. In the context of parallel computing, we face a very unbalanced problem since a large number of particles are created in a very fine meshed zone. Our strategy to load-balance the problem is to give a weight to every element group which is equal to the density of particle.

Final presentation of project #9

Project #10: Adaptive mesh refinement for turbulent premixed combustion

W. Agostinelli, O. Dounia, , T. Jaravel, O. Vermorel

The objective of the project was to evaluate the potential of adaptive mesh refinement (AMR) for premixed combustion in unsteady systems. Three target cases were identified: a semi-vented deflagration with laminar to turbulent transition, a planar detonation wave, and a bluff-body stabilized burner subjected to thermoacoustic oscillations. The simulations were performed with AVBP and coupled to the AMR implementation of YALES2. Several metrics and remeshing criterions were developed to identify and correctly resolve both the combustion wave front and the turbulent flow. The comparison of numerical results with reference simulations showed that the main features of the physics could be recovered with a significant speed-up in term of computational cost.

Final presentation of project #10

Project #11: Multiphysics coupling for wind turbine wake modeling

F.Houtin-Mongrolle (CORIA), B. Duboc (SGRE), P. Benard (CORIA)

The goal of this project was to evaluate the coupling of YALES2 (flow solver) and BHawC(Aero-Servo-Elastic solver).

Final presentation of project #11

Project #12: Stability of a semi-implicit compressible cavitation solver

The compressible cavitation solver is used to simulate cavitation inception in an initially liquid flow behind an obstacle. This solver is based on the implicit compressible solver, that has been modified to include a « barotropic » pressure-density relationship playing the role of an equation of state independent from the temperature. While this strategy has proven to be effective for DNS simulations of the implosion of vapour bubbles, the simulation of cavitation inception in an initially liquid flow was leading to strong instabilities in the simulation shortly after the appearence of vapour. The test case chosen is a flow behind a 2D cylinder. The analysis of the results has shown that instabilities were correlated with very low (and even unphysically negative) values of the pressure, that were triggering negative density values leading to code instability. Using limiters to ensure a positive pressure and a density within the range of the equation of state improved the stability and allowed to perform a preliminary simulation of a cavitating flow behind an obstacle. Ultimately instabilities appear anyways, so that the spatial discretisation is now under study, notably through upwinding.

Final presentation of project #12

Project #13: Validations and comparisons of Diffuse / Sharp interface methods in a structured DNS solver (Titan)

V. Boniou (EM2C), J.M. Dupays (EM2C), M. Pelletier (EM2C), T. Schmitt (EM2C), A. Vié (EM2C)

The project aimed at using academic test cases to compare the sharp (incompressible) and diffuse (compressible) models. In particular, the test case of an inviscid initially elliptical oscillating droplet has been carried out. The solvers features are the following:

- incompressible VOF solver (sharp): Numerical Method: Projection Method, Interface reconstruction: VOF, Surface tension: CSF

- compressible multifluid solver (diffuse): Advection scheme: MUSCL + RK2 + minmod limiter, Surface tension: CSF.

The source term is integrated with operator-splitting, and the curvature computation relies on a 2nd-order differentiation of the liquid volume fraction, which is previously smooth by filtering. This test case showed good agreement on the oscillation period, while exhibiting a slight numerical diffusion in the incompressible case and a strong numerical diffusion in the compressible case. In the compressible case, the use of higher-order splitting (Strang [SIAM Num. An. 1968]) has been tested, yielding no noticeable improvement. Reduction of the number of filtering iterations on the liquid volume fraction provides a slight improvement, which may indicate that a better curvature computation could participate to reduce the numerical diffusion.

Final presentation of project #13

Project #14: Méthode d'ordre élevé

M. Bernard (LEGI), G. Lartigue (CORIA), G. Balarac (LEGI), V. Moureau (CORIA)

Final presentation of project #14

Project #15: Utilisation d’éléments finis du second ordre dans le SMS

T. Fabbri (LEGI), G. Lartigue (CORIA), G. Balarac (LEGI), V. Moureau (CORIA)

Final presentation of project #15

Project #16: Development of a RANS solver in YALES2

G. Sahut (LEGI), G. Balarac (LEGI), V. Moureau (CORIA), G. Lartigue (CORIA), P. Bénard (CORIA), A. Grenouilloux (CORIA)

While the accuracy of LES usually approaches the one of DNS, LES are still too time-consuming for daily use in industrial applications. In this context, we started the development of a RANS solver in YALES2. We are first only interested in the steady state of the solution. In order to remove the CFL constraint, we developed, implemented and validated an implicit projection method for the resolution of the Navier-Stokes equations without turbulence models. The method is based on the implicitation of the velocity predictor ; the Poisson equation and the correction step of the velocity are then solved and applied as in the explicit incompressible solver. We validated the method on a stationary 2D Poiseuille flow with periodic boundary conditions: the simulation runs fine for CFL and Fourier numbers which are inaccessible with the explicit incompressible solver. The advection-diffusion equation for scalars has also been implicited and will be used to add turbulence models to the new implicit incompressible solver developped during this Workshop. More complex boundary conditions will also be addressed in a near future.

Final presentation of project #16

Project #17: COUPLING OF A FLUID PLASMA SOLVER WITH A LAGRANGIAN SOLVER FOR THE MODELING OF DUSTY

Final presentation of project #17

Project #18: L’Evaporo O Maıtre

Final presentation of project #18

Project #19: The Clone Wars

Final presentation of project #19

Project #20: Stiff complex fluid simulation with YALES2

Sam Whitmore, Yves Dubief, M2CE, University of Vermont

The objective was to simulate (1) ionized gases and (2) polymer solutions in flows using YALES2. Both problems are challenging owing to their stiff thermodynamics (1) or polymer dynamics (2). Significant gains were achieved in the implementation of the respective models thanks to the stiff integrator library CVODE. The plasma flow demonstrated an increase in time step of two orders of magnitude compared to previous implementation of the plasma chemistry in the variable density solver. Polymer models are notoriously prone to numerical instability. Again the use of CVODE showed equivalent if not superior stability of the solution at a fraction of the cost of commonly employed algorithms designed to address the stiffness of the problem.

Final presentation of project #20

Project #21: AVBP Dense Gases

Final presentation of project #21

Project #22: Numerical prediction of wind turbine wakes using AMR

Final presentation of project #22

Extreme CFD workshop