Extreme CFD workshop - User contributions [en]

Ecfd:ecfd 7th edition

2024-02-05T09:55:38Z

Corentin.raveleau: /* Projects */

{{DISPLAYTITLE: ECFD workshop, 7th edition, 2024}}

== Description ==

* Event from '''22th of January to 2nd of February 2024'''
* Location: [https://www.hotelclubdelaplage.com 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 70 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

[[File:ecfd7.png|600px|link=https://ecfd.coria-cfd.fr/index.php/Ecfd:ecfd_6th_edition]]

[[File:sponsor_ecfd7.png|text-bottom|600px]]



== Agenda ==

[[File:agenda_ecfd7.png|text-bottom|600px]]

== Thematics / Mini-workshops ==

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

To come...

== Projects ==

=== Hackathon GENCI - P. Begou, LEGI ===
The '''GENCI Hackathon''' will be devoted to porting two CFD codes to the Mi250 GPUs of the Adastra supercomputer deployed by GENCI at CINES.

For the '''YALES2''' code the goal is to obtain a first reference version giving the expected results then, if possible, to start its optimization to gain performance. The approach is OpenACC based with the objective of an implementation as least intrusive as possible in the existing code and which remains portable with the work done on the Nvidia GPUs of the Jean-Zay supercomputer at IDRIS.

The porting of the '''AVBP''' code is more advanced with a prototype already functional on Adastra but "hard-coded". The objective is to rationalize this first implementation, to integrate the latest developments in the code, to centralize memory management (host and device), to work on porting the Lagrangian part of the code and, of course, to improve the global performance.

This Hackathon is supported by GENCI, HPE, AMD and CINES with the presence on site of several development experts on AMD GPUS.

=== Mesh adaptation - R. Letournel, Safran ===

==== M1: ASMR for reheat chamber applications - Paul Pouech (CERFACS), Thibault Duranton & Luis Carbajal Carrasco (Safran) ====

Combustion in reheat chambers feature a wide range of lenght scales. Mesh refinement is thus mandatory to capture the flow characteristics within a reasonnable CPU cost for LES computations using the AVBP code. The purpose of this project is to consolidate mesh refinement criteria and strategy in an academic reference case. The retained workflow is supported by the [https://lemmings.readthedocs.io/en/latest/readme_copy.html Lemmings] code that calls the Tékigô wrapper for the mesh adaptations. During the ECFD7, the convergence time needed to have significant distribution of quantities of interest was analysed. An optimum runtime, based on a characteristic flow time-scale, was thus identified and led to a reduced running time for each adaptation step. As a second step, discussions with the ECFD7 participants led to the identification of interesting refinement criteria, namely the flame sensor or the mach rms for instance. Parametric analysis showed the robustness of the workflow based on a ponderation of different criteria. Finally, in order to facilitate the use of the workflow, efforts were made to improve the user experience by making it more human readable.

=== Numerics - S. Mendez, IMAG & G. Balarac, LEGI ===

==== WIP N1: High-Order schemes for Navier-Stokes - M. Bernard & G. Balarac (LEGI) & G. Lartigue (Total Eneriges) ====

Finite Volumes method.
Distinction between integrate and pointwise quantities.
Spacial accuracy depend on ability to evaluate accurately fluxes, starting from integrated quantities.
Accurate evaluation of fluxes through interface of each control volume (CV).
OK inside the domain, but chalenging in case of boundary conditions : by definition non uniformity of the stencil.
BC on faces instead of nodes.
Enforcement of the BC by modifying the normal component of the wall gradient in order to evaluate accurate diffusif flux.
Treatment of inlet BC by integrating convective flux through each boundary-facelet.

==== N5: Optimization of the RBC solver - F. Rojas & S. Mendez (IMAG) ====

==== N6: Electrodeformation of red blood cells, extension to 3D and improved accuracy at membrane - A. Spadotto & S. Mendez (IMAG), M. Bernard (LEGI) ====
The Leaky Dielectric Model is a popular framework to describe electric stresses over micro-scale membranes. We have adopted it to simulate the effect of a DC electric field on a red blood cell using the YALES2BIO solver. The goal of the project is to reproduce the electric charging process of the membrane, as well as the resulting stresses, which may yield to electrodeformation of the cell. From the point of view of the implementation, the grid is represented by a 2D surface mesh embedded in a 3D eulerian grid. The need to make variables stored on the surface interact with quantities stored on the Eulerian grid calls for a proper bidirectional 2D-membrane/3D-grid dynamic connectivity. The advancement of theis task during this ECFD has led to the first 3D simulation of a charging fixed spherical shell. Moreover, the estimation of grid variables on elements cut by the membrane has been improved thanks to a High-Order extrapolation. The latter has been successfully tested on 2D configurations. The project opens the way for a series of validation tests. In particular, future work will demand treatment of instabilities emerging in symmetrical configurations.

=== Turbulence - P. Benard, CORIA & L. Bricteux, UMONS ===

=== T4: Atmospheric solver ===
Wind turbines, bigger and bigger, are now influenced by atmospheric flows. An atmospheric solver has already been developed in YALES2 to represents some of its effects (Coriolis, veer, thermal stratification). In this continuum, the project has been divided into two work-packages.
- Work-package 1: The use of the Variable density solver (VDS).
Before ECFD7, thermal stratification was taken into account using the Boussinesq buoyancy approximation within the incompressible solver framework. Now, VDS can be used, taking into account all thermal effect. Results are promissing.
- Work-package 2: Wall law velocity filtering.
Wall law are using velocity at the first grid node to compute wall shear stress. Before ECFD7, atmospheric wall law were using the local velocity, leading sometimes to convergence errors. Now a gather-scatter filter can be used to average velocity (and temperature) at first grid node.

=== Two Phase Flow - M. Cailler, Safran Tech & V. Moureau, CORIA ===

==== P3: Blood platelets adhesion model - C. Raveleau, S. Mendez, F. Nicoud (IMAG) ====

Medical devices in contact with blood (e.g. artificial valves) are used to treat various cardiovascular diseases, but their thrombogenicity remains the main unresolved issue in their development. A numerical model of blood platelets is being constructed to help to understand the effect of microstructuration on the thrombogenicity of artificial surface. The Force Coupling Method (FCM) was previously implemented and allows the modelisation of ellipsoidal particle and their interaction with the surrounding fluid. During the workshop, the particle model was extended to include adhesive and repulsive interactions with walls or with other particles. The adhesive bonds are modeled with springs forming when the distance between a node of a particle surface and a node of the wall or another particle is smaller than a given threshold. The stiffness of the bond is increased after a given formation time to mimic the 2-step adhesion process of platelets to von Willebrand Factor. A Lennard-Jones potential was used to model the collision of particles. Future work will aim at generalizing these implementations for an arbitrary number of particles (currently only working for 2 particles) and ensuring the interactions are unaltered by the crossing of a periodic boundary.

=== Combustion - K. Bioche, CORIA & R. Mercier, Safran ===

=== User Experience & Data - L. Korzeczek, GDTech ===

Ecfd:ecfd 6th edition

2023-02-06T18:37:57Z

Corentin.raveleau: /* Numerics - S. Mendez, IMAG & M. Bernard, LEGI */

{{DISPLAYTITLE: ECFD workshop, 6th edition, 2022}}

== Description ==
{| align="right" style="text-align:center;" cellpadding="2"
| [[File:Logo_ECFD6.png | center | thumb | 350px | ECFD6 workshop logo.]]
|}
* Event from '''23th of January to 3rd of February 2023'''
* Location: [https://www.hotelclubdelaplage.com 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

[[File:Banniere_ECFD6.png|600px|link=https://ecfd.coria-cfd.fr/index.php/Ecfd:ecfd_6th_edition]]

[[File:Banniere_ECFD6_sponso.png|text-bottom|600px]]



== Agenda ==

[[File:ECFD6_program.png|text-bottom|600px]]

== 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)'''
According to the Leaky Dielectric Model, red blood cells (RBCs) are subject to a force which is proportional to the jump of the Maxwell tensor. This latter is a quantity scaling as the square of the electric field, which under the quasi-static hypothesis is defined as the gradient of the electrostatic potential. To work out the potential, an elliptic interface problem must be solved, taking into account the presence of the RBC membrane. The aim of the project was implementing the Ghost Fluid Method (GFM) to face the interface problem. Good results were obtained on unstructured meshes. Secondly, a gradient calculation was performed applying the Green-Gauss scheme, modified in the style of the GFM. Future work will focus on interpolation of the gradient field onto the membrane to get an estimation of the effort. Possibly, high-order schemes for the gradient calculation will be explored. In a second time, the effort calculation will be merged into an Immersed Boundary solver for the RBC dynamics.

* '''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)'''
The Force Coupling Method (FCM) allows the representation of particles in flows from 0 to finite Reynolds number based on a regularization of the multipole expansion for Stokes flows, providing enough particle resolution to match the exact Stokes flow away from the particle. This is done by using a Gaussian function to spread the singular forces over a domain whose size is derived from physical quantities of interest such as the settling velocity of a particle under gravity in Stokes flow and the particle radius. During the workshop, the first term of the multipole expansion (monopole) and the antisymetric part of the second term (antisymetric dipole) were implemented and validated against test cases from the litterature. The monopole was tested in the case of a particle held fixed against an incoming flow and the dipole was tested by applying a torque on a particle in a still fluid. In both cases, the velocity profiles matched the results from the litterature and approximated well the Stokes solution at a distance of about 1.5 radii from the particle center. It was also compared to the conservative immersed boundary method (CLIB solver) implemented in YALES2 for the monopole test case. Both methods give good results, although the FCM is able to predict the expected solution with a coarser mesh than CLIB. Cases where the particle moves are under investigation.

* '''Sub-project 5: Breaking limitations of the linearized implicit time advancement (T. Berthelon, G. Lartigue, G. Balarac)'''
The explicit time advancement classically used in ics solver is limited by CFL constraint. In order to get ride of this constraint, an implicit time advancement method, based on the linearization of the convective term, has been recently developed.
However, the method is limited by difficulties to solve linear system, with the BiCGSTAB2 algorithm, during the prediction step. The objective of this project was to understand these limitations. The correction of a bug on the boundary conditions (viscosity imposed at zero) was identified. In addition, the spatial scheme and the presence of a buffer zone at the end of the domain showed a great influence on the convergence of the prediction. The perspectives for a more robust and efficient use of this temporal integration consist in working on the spatial schemes and on the pre-conditioning.

* '''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: Aerodynamics of floating wind turbines (R. Amaral, F. Houtin-Mongrolle, E. Muller)'''

* '''Sub-project 5: Implementation of the Risoe dynamic stall model for YALES2 (V. Maronnier, E. Muller, B. Duboc)'''

* '''Sub-project 6: YALES2-OpenFAST coupling (A. Parinam, A. Viré, D. Von Terzi, B. Duboc, F. Houtin-Mongrolle, P. Bénard)'''

* '''Sub-project 7: 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.

=== Two Phase Flow - C. Merlin, Ariane Group & M. Cailler, Safran Tech ===

* '''Sub-project 1: Convergent computation of interface curvature (G. Ghigliotti, M. Benard, G. Balarac, J. Carmona, R. Mercier, G. Lartigue)'''

Though Level-set distance evaluation through GPMM (Janodet et al., 2022) converges at order 2, the interface curvature convergence is as best 0 using the non-compact Goldman formulation.
Following progresses obtained during ECFD5, a strategy based on parabolic fit of the interface has been explored during the workshhop. This method aims at fitting a parabola through least squares using the interface markers stored in the interface vicinity. First the method was applied on a 2-D perfectly spherical droplet with exact projection of the marker on the circle. This results in a first order convergent curvature. Without projection of the markers, the fiting strategy allows a slight decrease of the error but no improve on the curvature convergence order in comparison with the standard non-compact formulation. As a persective, these results will be validated on dynamic and 3-D cases (MMG3D meshes). Also, the sensitivity on the number of markers and their redundancy will be investigated.

=== Combustion - K. Bioche, CORIA & R. Mercier, Safran Tech ===

* '''Sub-project 1: Multi-regime F-TACLES (S. Dillon, R. Mercier, B. Fiorina)'''

Filtered tabulated chemistry for large eddy simulations is currently a common tool to model premixed flames or diffusion flames. Tabulation using 1D counterflow flames, as a function of the mixture fraction and progress variable, was previously tested on laminar and turbulent cases. It resulted in difficulties to describe the outer mixing zone and yield a very stiff evolution of SGS source terms in the phase space. The model was modified to include the mixture fraction scalar dissipation rate as a table dimension. This solves previous limitations, but using 1D counterflow flames yields empty table zones, making the method numerically infeasible. Tabulation using both 1D counterflow flames and 1D partially premixed flames gives well-built tables, and was tested on 1D flames for various strain rates.

* '''Sub-project 2: Limiter model for turbulence combustion interaction in MILD combustion (E. Stendardo, L. Bricteux, K. Bioche)'''

MILD combustion yields intense turbulence and widespread reaction zones, requiring expensive mesh refinement over large areas. Practical mesh won’t be fine enough, leading to sub-grid heterogeneousness and effects of sub-grid turbulent fluctuations. A generic limiter type combustion model was implemented to solve for turbulence combustion interaction. This family of models includes Partially Stirred Reactor, Quasi Laminar Finite Rate and Laminar Finite Rate models. In these models, the source term is multiplied by a limiter factor and the residence time in inner cell reactive structure can be modelled. This implementation will permit testing the various limiter formulations in the future.

* '''Sub-project 3: Evaluate spatial discretization schemes on scalar transport (K. Bioche, Y. Bechane, R. Mercier, G. Lartigue, V. Moureau, J. Carmona, M. Bernard, L. Voivenel)'''

Common practice in combustion solvers is to use centred spatial schemes. Such low-dissipation schemes can prove unstable when applied to under-resolved scalar transport in presence of strong gradients. This is typically the case for H2/air combustion. Initial low-resolution simulations require thus adapted numerical schemes. Various spatial schemes were evaluated on the scalar transport problem, including: 4th order, 3rd order, 2nd order, WENO3, high order schemes, MUSCL schemes with various limiters (overbee, superbee, sweby, van leer, minmod). Their application to various configurations was discussed to emphasize on their robustness and accuracy. Tests cases include: 1D scalar convection Jiang Shu test case, 2D scalar bump convection for convergence analysis, a 2D reactive Bunsen burner and finally the 3D Preccinsta burner.

* '''Sub-project 4: Phenomenological plasma model for reacting systems (S. Wang, Y. Bechane, B. Fiorina)'''

Plasma assisted combustion consists in stabilizing flames in near extinction conditions thanks to electric discharges. Stabilization of lean premixed flames with Nanosecond Repetitively Pulsed electric Discharges is a strategy to reduce NOx emissions. Full 3D simulations of plasma assisted combustion are extremely expensive, so that the use of a semi-empirical strategy to model NRPD is preferred in CFD solvers. During the workshop, Castela’s model was implemented in a variable density solver. This model was extended to an explicit compressible solver. The model of Blanchard was also implemented in both frameworks. A 2D pin-to-pin configuration was successfully simulated with both models and frameworks.

* '''Sub-project 5: Development and assessment of combustion in an explicit compressible solver (Y. Bechane, L. Voivenel, R. Mercier, K. Bioche)'''

The implementation of reactive physics in the Explicit Compressible Solver (ECS) of the YALES2 platform was undertaken. To this aim, reactive gases thermochemical functions were implemented. Specific schemes were developed to increase the temperature and species diffusion schemes from 2nd to 4th order. Finally, a 2D methane-air Bunsen flame was simulated with low order numerical schemes (RK1 and SLAU2).

* '''Sub-project 6: Clustering for finite rate chemistry using PCA (R. Mercier, A. Stock)'''
To reduce the cost of species source terms computation, a clustering method was adopted. It consists in detecting nodes with similar properties and compute chemical source terms only once for these. Still, considering each species in this process creates a high dimensional cluster, while replacing species by a user-set progress variable may not well describe species. The strategy adopted here relies on the application of a PCA on species. It can be viewed as an automated “progress variable” creation. The use of such strategy was shown to reduce the simulation cost of source term computation by a factor 6 on a simple 2D flame ball case.

=== User Experience - J. Leparoux, Safran Tech & A. Pushkarev, GE Renewable Energy===

* '''Sub-project 1: External coupling with CWIPI (R. Letournel, V. Moureau, C. Merlin, M. Cailler, P. Bégou, S. Meynet)'''

* '''Sub-project 2: Automated Grid Convergence refactoring (J. Leparoux, M. Cailler, R. Letournel)'''

* '''Sub-project 3: Advanced Liquid spray post-processing (J. Carmona, J. Leparoux, N. Gasnier, C. Brunet, I. El Yamani)'''

* '''Sub-project 4: YALES2 as industrial solver for GE design optimization tools (A. Pushkarev, H. Lam, G. Balarac)'''

* '''Sub-project 5: YALES2 History and Geography (T. Marzlin, A. Dauptain, P. Bénard)'''

* '''Sub-project 6: Improve the HT solver: refactoring of linear solver operators & Robin BC (C. Merlin, V. Moureau, R. Letournel)'''