Difference between revisions of "Ecfd:ecfd 6th edition"

From Extreme CFD workshop
Jump to: navigation, search
(Turbulence - P. Benard, CORIA & G. Balarac, LEGI)
(Numerics - S. Mendez, IMAG & M. Bernard, LEGI)
Line 52: Line 52:
  
 
* '''Sub-project 1: Multi-level domain decomposition method (DDM) for coupled systems of differential-algebraic equations (A. Quirós Rodrígues, V. Le Chenadec)'''
 
* '''Sub-project 1: Multi-level domain decomposition method (DDM) for coupled systems of differential-algebraic equations (A. Quirós Rodrígues, V. Le Chenadec)'''
 
+
The numerical approximation of multi-physics problems gives rise to complex linear systems, the solution of which leverages preconditioning techniques such as multi-grid or domain decomposition methods. This project aimed at coupling two Julia packages that being actively developed: a two-dimensional Navier-Stokes solver for free-surface and two-phase flows (Flower.jl) on the other, and a Domain Decomposition package for Cartesian grids (DDM.jl). The decomposed matrix-vector product was optimised to reduce the overhead associated with halo exchanges. The implementation of a deflated Conjugate Gradient as well as one- and two-level Additive Schwartz Method were also completed and shown to significant reduce the number of iterations for inverting monolithic systems (i.e. without resorting to operator splitting), shown to be independent of the number of subdomains for constant property flows. Future work will focus on a further optimisation of the implementation for vectorisation and multi-threading, and extension of the deflation to generalised coarse spaces to support highly discontinuous transport properties (GenEO).
  
 
* '''Sub-project 2: Ghost fluid method (GFM) for Electrodeformation (A. Spadotto , S. Mendez)'''
 
* '''Sub-project 2: Ghost fluid method (GFM) for Electrodeformation (A. Spadotto , S. Mendez)'''

Revision as of 20:10, 4 February 2023


Description

ECFD6 workshop logo.
  • Event from 23th of January to 3rd of February 2023
  • Location: Hôtel Club de la Plage, Merville-Franceville, near Caen (14)
  • Two types of sessions:
    • common technical presentations: roadmaps, specific points
    • mini-workshops. Potential workshops are listed below
  • Free of charge
  • More than 60 participants from academics, HPC center/experts and industry.
  • Objectives
    • Bring together experts in high-performance computing, applied mathematics and multi-physics CFDs
    • Identify the technological barriers of exaflopic CFD via numerical experiments
    • Identify industrial needs and challenges in high-performance computing
    • Propose action plans to add to the development roadmaps of the CFD codes


Banniere ECFD6.png


Banniere ECFD6 sponso.png



Agenda

ECFD6 program.png

Thematics / Mini-workshops

These mini-workshops may change and cover more or less topics. This page will be adapted according to your feedback.

Projects

Hackathon - G. Staffelbach, CERFACS & P. Begou, LEGI

Mesh adaptation - R. Letournel, Safran & V. Moureau, CORIA

Numerics - S. Mendez, IMAG & M. Bernard, LEGI

  • Sub-project 1: Multi-level domain decomposition method (DDM) for coupled systems of differential-algebraic equations (A. Quirós Rodrígues, V. Le Chenadec)

The numerical approximation of multi-physics problems gives rise to complex linear systems, the solution of which leverages preconditioning techniques such as multi-grid or domain decomposition methods. This project aimed at coupling two Julia packages that being actively developed: a two-dimensional Navier-Stokes solver for free-surface and two-phase flows (Flower.jl) on the other, and a Domain Decomposition package for Cartesian grids (DDM.jl). The decomposed matrix-vector product was optimised to reduce the overhead associated with halo exchanges. The implementation of a deflated Conjugate Gradient as well as one- and two-level Additive Schwartz Method were also completed and shown to significant reduce the number of iterations for inverting monolithic systems (i.e. without resorting to operator splitting), shown to be independent of the number of subdomains for constant property flows. Future work will focus on a further optimisation of the implementation for vectorisation and multi-threading, and extension of the deflation to generalised coarse spaces to support highly discontinuous transport properties (GenEO).

  • Sub-project 2: Ghost fluid method (GFM) for Electrodeformation (A. Spadotto , S. Mendez)


  • Sub-project 3: Optimization of the high order framework (HOF) for Navier-Stokes incompressible (M. Bernard, P. Bégou, G. Lartigue, G. Balarac)

Over the past years, a framework has been developed to improve the spatial accuracy of numerical schemes on distorted meshes. However, even if the solution is more precise, the computational cost of the overall resolution of Navier-Stokes equations is large. As a consequence, HOF becomes profitable only on thin meshes thanks to a better spatial convergence order. The code has been analized with different analysis tools (MAQAO, Gprof, Scalasca). The main time consuming routines have been identified and improved. Moreover, some algorithms have been refactors such that the resolution of Navier-Stokes equations has been speed-up by a factor 2.5.


  • Sub-project 4: Force coupling method (FCM) for particulate flows (C. Raveleau, S. Mendez)


  • Sub-project 5: Breaking limitations of the linearized implicit time advancement (T. Berthelon, G. Lartigue, G. Balarac)


  • Sub-project 6: Development of a traction open boundary condition (TOBC) in Yales2 (J.B. Lagaert, Guillaume Balarac)


  • Sub-project 7: Development of new spatial differential operators in Yales2 (M. Bernard, G. Lartigue)

It exists different philosophies for computing differential operators on distorted meshes. In a HPC context, the 2 main approaches are the Green-Gauss operators and the Least-Squares operators. During ECFD#6, 2 new types of "non-compact" Hessian operators have been implemented by computing successively the gradient operator, eather with Green-Gauss gradient, or with Least-Squares gradient. Those operators lead to good convergence order, even on distorted mehes. However, their application on low-resolution signals lead to large error magnitude due to their extended stencil.\\ Another pair of gradient \& hessian Least-Squares operators have been implemented, leading to 2nd and 1st order accuracy for the gradient and hessian respectively. Those operators have very interesting characteristics as their stencil is restricted to the direct neighbors only and their computational cost remains low.


  • Sub-project 8: DOROTHY optimization (M. Roperch, H. Mulakaloori, G. Pinon, P. Bénard)


  • Sub-project 9: Anamika, a tool to improve programming productivity (K. Mohana Muraly, G. Staffelbach)

Turbulence - P. Benard, CORIA & G. Balarac, LEGI

  • Sub-project 1: Explore hybrid RANS/LES strategies (T. Berthelon, G. Balarac)

For complex industrial applications, LES can still lead to a too long restitution time. In other hand, statistical approaches can lead too a lack of accuracy. In this project, the potentiality of hybrid approaches combining both have been explored. Conventional hybrid RANS/LES approaches consider a unique solution field, with an unique transport equation and a clusre terme modeled using RANS or LES models depending of the regions. The main idea is to evaluate a strategy based on a separation between mean fields and fluctuations with distinct coupled transport equations. First elements of validation using YALES2 code are shown that it was possible to correct the prediction of a RANS models, by performing LES of the fluctuations. Next steps should be to consider disctinct meshes, or even computational domains for RANS and LES with this strategy.


  • Sub-project 2: Flow Instabilities over Rotating curved Surfaces (S. Sawaf, M. Shadloo, A. Hadjadj, S. Moreau, S. Poncet)

For evaluating the effect of the clearance between the blade tip and the casing of axial ducted fans on noise emissions, LES offers excellent tool to capture the consitricted flow around the blade tip especially for small clearances where RANS fails because of unsteady flow conditions. LES simulation of the aerodynamics is the first step toward extracting accoustics data that helps to improve the design of axial ducted fans so they comply with the noise emission regulations in admistrative buildings. noise emmisions are estimated using analytical aeroacoustic models informed by data that are extracted from the LES simulations.


  • Sub-project 3: Automatic statistical convergence metric (C. Papagiannis, G. Balarac, O. Le Maitre, P. Congedo)

Statistics accumulation can be an important part of the restitution time in unsteady simulations (DNS/LES). In this project, the goal was to estimate uncertainties on the "finite time statistics". For time correlated data, it can be shown that the variance of the mean estimator (i.e. the fluctuation of the estimation of the mean) is dependent of the correlation time. Modeling this correlation time based on the integral time scale of the turbulence appears as a first way to define a practical metric to evaluate the statistic convergence on-fly during simulations. Next step should be to explore procedures to accelerate the statistics accumulation step.


  • Sub-project 4: Wall law for Immersed Boundaries & Rough surfaces (M. Cailler, A. Cuffaro, P. Benez, S. Meynet)

Conservative Lagrangian Immersed Boundaries (CLIB) are now a useful way to take into account complex geometries in YALES2. During the workshop, a brand-new data-structure for modular and generic immersed-body has been developed. This data-structure paves the way for various new capabilities for IB methods: penalization mask shape optimization for improved velocity imposition, better control of near wall discretization based on a reliable evaluation of wall units, wall-modeling, etc... For this purpose the periodic hill test case has been considered. Simulations of this configuration has been performed by using body-fitted meshes, and CLIB for both smooth and rough surfaces. This will allow to assess the accuracy of the IB methods, and will constitute a database for IB models improvement, and the development of wall-modeling strategies.


  • Sub-project 8: Atmospheric flow (U. Vigny, L. Voivenel, P. Benard, S. Zeoli)

Atmospheric flow such as Atmospheric Boundary Layer (ABL) and thermal stratification have an impact on wind turbines aerodynamic and wakes. Mostly at a wind farm scale, the change of wind turbine wake size and recovery can modify the global power production. During the workshop, the Coriolis force implementation has been validated through neutral case (where no thermal stratification i.e. no temperature gradient). It also allowed to validate the pressure forcing term, needed to drive the flow in a periodic box. YALES2 results showed a good agreement with other numerical and experimental results. Afterwards, the stable case (i.e. temperature gradient downwards) has been studied. A surface temperature as boundary condition has been developed. Yet, results are not as expected and further investigation is needed.