EN FR
EN FR


Section: Application Domains

Energy conversion

Fluid-structure interaction

We apply the methods developed in our team to the domain of wind engineering and sea-wave converters. In Figure 1, we show results of a numerical model for a sea-wave energy converter. We here rely on a monolithic model to describe the interaction between the rigid floater, air and water; material properties such as densities, viscosities and rigidity vary across the domain. The appropriate boundary conditions are imposed at interfaces that arbitrarily cross the grid using adapted schemes built thanks to geometrical information computed via level set functions [46]. The background method for fluid-structure interface is the volume penalization method [27] where the level set functions is used to improve the degree of accuracy of the method [3] and also to follow the object. The underlined mathematical model is unsteady, and three dimensional; numerical simulations based on a grid with 𝒪(108) degrees of freedom are executed in parallel using 512 CPUs .

Figure 1. numerical modeling of a sea-wave converter by a monolithic model and Cartesian meshes.
IMG/energy_conversion.png

In the context of the Aerogust (Aeroelastic gust modelling) European project, together with Valorem, we investigated the behavior of wind turbine blades under gust loading. The aim of the project was to optimize the design of wind turbine blades to maximize the power extracted. A meteorological mast (Figure 2(a)) has been installed in March 2017 in Brittany to measure wind on-site: data provided by the mast have been exploited to initialize the mathematical model. Due to the large cost of the full-order mathematical model, we relied on a simplified model [38] to optimize the global twist. Then, we validated the optimal configuration using the full-order Cartesian model based on the NaSCar solver. Figure 2(b) shows the flow around the optimized optimized wind turbine rotor.

Figure 2. Aerogust project. Left: met mast after its installation. Right: flow around the optimized wind turbine rotor (as predicted by NaSCar).
IMG/met_mast.png IMG/wind_turbine.png

Schemes for turbulent flow simulations using Octrees

We have initially developed and tested a 3D first-order Octree code for unsteady incompressible Navier-Stokes equations for full windmill simulations with an LES model and wall laws. We have validated this code on Occigen for complex flows at increasing Reynolds numbers. This step implied identifying stable and feasible schemes compatible with the parallel linear Octree structure. The validation has been conducted with respect to the results of a fully Cartesian code (NaSCAR) that we run on Turing (with significantly more degrees of freedom) and with respect to experimental results.

Subsequently, we have developed a second-order Octree scheme that has been validated on Occigen for a sphere at a moderate Reynolds number ( Re =500), see Table 1. Then, for a cylinder at ( Re =140000) (Figures 3(a) and 3(b)), close to real applications, we have preliminary validation results for the second-order scheme with respect to experimental drag coefficient (Table 2). Additional resources will be asked on Occigen to complete the study.

Table 1. Flow past a sphere at Re =500. Results in the literature are spread between C D =0.48 and C D =0.52.
Mesh Δ x min number of cells CD (1st-order scheme) CD (2nd-order scheme)
1 0 . 094 0 . 72 · 10 5 N.A. 0 . 526
2 0 . 047 4 . 9 · 10 5 0 . 595 0 . 522
3 0 . 023 4 . 7 · 10 6 0 . 546 0 . 492
4 0 . 012 37 . 6 · 10 6 0 . 555 0 . 496
Table 2. Flow past a sphere at Re =14000.
Case C D
Octree, 1st-order scheme 1 . 007
Octree, 2nd-order scheme 1 . 157
Cartesian 1 . 188
Experimental estimate [34] 1 . 237
Figure 3. flow past a cylinder at Re =140000. Left: vorticity contour lines. Right: streamwise velocity section and grid for the second-order Octree scheme.
IMG/contour.png IMG/streamline_grid.png