Laboratoire de Recherche Hydrodynamique, Energétique et Environnement Atmosphérique

International audience

This paper presents a brief review of grand challenges of Smoothed Particle Hydrodynamics (SPH) method. As a meshless method, SPH can simulate a large range of applications from astrophysics to free-surface flows, to complex mixing problems in industry and has had notable successes. As a young computational method, the SPH method still requires development to address important elements which prevent more widespread use. This effort has been led by members of the SPH rEsearch and engineeRing International Community (SPHERIC) who have identified SPH Grand Challenges. The SPHERIC SPH Grand Challenges (GCs) have been grouped into 5 categories: (GC1) convergence, consistency and stability, (GC2) boundary conditions, (GC3) adaptivity, (GC4) coupling to other models, and (GC5) applicability to industry. The SPH Grand Challenges have been formulated to focus the attention and activities of researchers, developers, and users around the world. The status of each SPH Grand Challenge is presented in this paper with a discussion on the areas for future development.

In the present work, an in-depth analysis of the theoretical structure of the smoothed-particle hydrodynamics (hereinafter SPH) is provided for an inviscid, weakly compressible, and barotropic flow in the presence of a free surface. The role of the free surface in the SPH scheme is indeed little addressed in literature. In the present analysis, the general continuous formulation of the SPH method is considered. A detailed description of the free-surface influence on the smoothed differential operators is supplied. New and existing forms are analyzed in detail, in terms of convergence and conservation properties. The proposed analysis is based on the principle of virtual works, which permits to exhibit the link with the enforcement of the dynamic free-surface boundary condition. Finally, possible SPH formulations resulting from this analysis are investigated, in terms of consistency, conservation, and dynamic free-surface boundary condition.

The combustion of synthesis gas will play an important role in advanced power systems based on the gasification of fuel feedstocks and combined cycle power production. While the most commonly discussed option is to burn syngas in gas turbine engines, another possibility is to burn the syngas in stationary reciprocating engines. Whether spark ignited or compression ignited, syngas could serve to power large bore stationary engines, such as those presently operated on natural gas. To date, however, there has been little published on the combustion of syngas in reciprocating engines. One area that has received attention is dual-fueled diesel combustion, using a combination of diesel pilot injection and syngas fumigation in the intake air. In this article, we survey some of the relevant published work on the use of synthesis gas in IC engines, highlighting recent work on dual-fuel (syngas + diesel) combustion.

Experiments have been performed in the Shallow Water Wave Basin of DHI (Hørsholm, Denmark), on large arrays of up to 25 heaving point absorber type Wave Energy Converters (WECs), for a range of geometric layout configurations and wave conditions. WEC response and modifications of the wave field are measured to provide data for understanding WEC array interactions and to evaluate array interaction numerical models. Each WEC consists of a buoy with a diameter of 0.315 m and power take-off (PTO) is modeled by realizing friction based energy dissipation through damping of the WEC’s motion. Wave gauges are located within and around the WEC array. Wave conditions studied include regular, polychromatic, long- and short-crested irregular waves. A rectilinear arrangement of WEC support structures is employed such that several array configurations can be studied. In this paper, the experimental arrangement and the obtained database are presented. Also, results for wave height attenuation downwave a rectilinear array of 25 heaving WECs are presented, for the case of irregular waves. Up to 16.3% and 18.1% (long-crested) and 11.2% and 18.1% (short-crested waves) reduction in significant wave height is observed downwave the WEC array, for the radiated wave field only and for the combination of incident-diffracted-radiated (perturbed) wave field, respectively. Using spectra at different locations within and around the array, the wave field modifications are presented and discussed.

The objective of this paper is to investigate cavitating flows around a pitching hydrofoil via combined physical and numerical studies. The aims are to (1) improve the understanding of the interplay between unsteady cavitating flow, hydrofoil motion, and hydrodynamic performance, (2) quantify the influence of pitching rate on subcavitating and cavitating responses, and (3) quantify the influence of cavitation on the hydrodynamic load coefficients and surrounding flow structures. Results are presented for a NACA66 hydrofoil undergoing controlled, slow $(\dot \alpha = 6^\circ /{\rm s})$(α̇=6∘/s) and fast $(\dot \alpha = 63^\circ /{\rm s})$(α̇=63∘/s) pitching motions from α = 0° to α = 15° and back to α = 0° for both subcavitating and cavitating conditions at a moderate Reynolds number of Re = 750 000. The experimental studies were conducted in a cavitation tunnel at the French Naval Academy, France. The numerical simulations are performed by solving the incompressible, multiphase Unsteady Reynolds-Averaged Navier-Stokes Equations via the commercial code CFX using a transport equation-based cavitation model; a modified k-ω SST turbulence model is used to account for the effect of local compressibility on the turbulent eddy viscosity. The results showed that increases in the pitching rate suppressed laminar to turbulent transition, delayed stall, and significantly modified post-stall behavior. Cavitation inception at the leading edge modified the pressure distribution, which in turn significantly changed the interaction between leading edge and trailing edge vortices, and hence the magnitude as well as the frequency of the load fluctuations. For a fixed cavitation number, increases in pitching rate lead to increase in cavitation volume, which in turn changed the cavity shedding frequencies and significantly modified the hydrodynamic loads. Inversely, the leading edge cavitation observed for the low pitching velocity case tends to stabilize the stall because of the decrease of the pressure gradient due to the formation of the cavity. The results showed strong correlation between the cavity and vorticity structures, which suggest that the inception, growth, collapse and shedding of sheet/cloud cavities are important mechanisms for vorticity production and modification.

The ability of linear stochastic response analysis to estimate coherent motions is investigated in turbulent channel flow at the friction Reynolds number $\text{Re}_{\unicode[STIX]{x1D70F}}=1007$ . The analysis is performed for spatial scales characteristic of buffer-layer and large-scale motions by separating the contributions of different temporal frequencies. Good agreement between the measured spatio-temporal power spectral densities and those estimated by means of the resolvent is found when the effect of turbulent Reynolds stresses, modelled with an eddy-viscosity associated with the turbulent mean flow, is included in the resolvent operator. The agreement is further improved when the flat forcing power spectrum (white noise) is replaced with a power spectrum matching the measures. Such a good agreement is not observed when the eddy-viscosity terms are not included in the resolvent operator. In this case, the estimation based on the resolvent is unable to select the right peak frequency and wall-normal location of buffer-layer motions. Similar results are found when comparing truncated expansions of measured streamwise velocity power spectral densities based on a spectral proper orthogonal decomposition to those obtained with optimal resolvent modes.

During the process of wave slamming on a structure with sharp corners, the wave receding after wave impingement can induce strong negative pressure (relative to the atmospheric pressure) at the bottom of the structure, which is called the suction effect. From the practical point of view, the suction force induced by the negative pressure, coinciding with the gravity force, pulls the structure down and hence increases the risk of structural damage. In this work, the smoothed particle hydrodynamics (SPH) method, more specifically the δ+SPH model, is adopted to simulate the freak wave slamming on a fixed platform with the consideration of the suction effect, i.e., negative pressure, which is a challenging issue because it can cause the so-called tensile instability in SPH simulations. The key to overcome the numerical issue is to use a numerical technique named tensile instability control (TIC). Comparative studies using SPH models with and without TIC will show the importance of this technique in capturing the negative pressure. It is also found that using a two-phase simulation that takes the air phase into account is essential for an SPH model to accurately predict the impact pressure during the initial slamming stage. The freak wave impacts with different water depths are studied. All the multiphase SPH results are validated by our experimental data. The wave kinematics/dynamics and wave impact features in the wave-structure interacting process are discussed, and the mechanism of the suction effect characterized by the negative pressure is carefully analyzed.

The Smoothed Particle Hydrodynamics (SPH) method, often used for the modelling of the Navier–Stokes equations by a meshless Lagrangian approach, is revisited from the point of view of Large Eddy Simulation (LES). To this aim, the LES filtering procedure is recast in a Lagrangian framework by defining a filter that moves with the positions of the fluid particles at the filtered velocity. It is shown that the SPH smoothing procedure can be reinterpreted as a sort of LES Lagrangian filtering, and that, besides the terms coming from the LES convolution, additional contributions (never accounted for in the SPH literature) appear in the equations when formulated in a filtered fashion. Appropriate closure formulas are derived for the additional terms and a preliminary numerical test is provided to show the main features of the proposed LES-SPH model.

Soliton gases represent large random soliton ensembles in physical systems that exhibit integrable dynamics at the leading order. Despite significant theoretical developments and observational evidence of ubiquity of soliton gases in fluids and optical media, their controlled experimental realization has been missing. We report a controlled synthesis of a dense soliton gas in deep-water surface gravity waves using the tools of nonlinear spectral theory [inverse scattering transform (IST)] for the one-dimensional focusing nonlinear Schrödinger equation. The soliton gas is experimentally generated in a one-dimensional water tank where we demonstrate that we can control and measure the density of states, i.e., the probability density function parametrizing the soliton gas in the IST spectral phase space. Nonlinear spectral analysis of the generated hydrodynamic soliton gas reveals that the density of states slowly changes under the influence of perturbative higher-order effects that break the integrability of the wave dynamics.

SUMMARY In this paper, an in‐depth study of SPH method, in its original weakly compressible version, is achieved on dedicated 2D and 3D free‐surface flow test cases. These rather critical prototype problems shall constitute suitable test cases to get through when building a free‐surface SPH model. The present work aims at investigating various numerical aspects of this method, often little mentioned in literature. In particular, a great care is paid to the dynamic part of the solution, which is critical to the local hydrodynamic load prediction. The role of numerical errors in the development of acoustic frequencies in the pressure signals is discussed, as well as the influence of the choice of the sound velocity. On the shown test problems, it is also evidenced that some numerical tools are crucial to ensure the robustness and accuracy of the standard SPH method. The convergence of our model is heuristically proved on these nonlinear prototype tests, showing at the same time the very satisfactory level of accuracy reached. Through these tests, some other numerical specificities of the SPH method are discussed, such as the self‐redistribution of the particles occurring during the Lagrangian evolution. A higher order model is also proposed, and its advantages and drawbacks are discussed. Copyright © 2013 John Wiley & Sons, Ltd.

The monitoring of the water cycle at the Earth surface which tightly interacts with the climate change processes as well as a number of practical applications (agriculture, soil and water quality assessment, irrigation and water resource management, etc…) requires surface temperature measurements at local scale. Such is the goal of the Indian-French high spatio-temporal TRISHNA mission (Thermal infraRed Imaging Satellite for High-resolution Natural resource Assessment). The scientific objectives of the mission and research work conducted to consolidate the mission specifications are presented. Progress in modelling of surface fluxes is then discussed. The main specifications of the mission such as the revisit, the spatial resolution, the overpass time, the spectral bands and the orbit are analyzed and justified. The resulting baseline of the mission is given.

Abstract A series of wave basin experiments is reported that investigates the flexural response of one or two floating thin elastic discs to monochromatic waves. The work is motivated by numerical model validation. Innovative techniques are used to ensure the experimental configuration is consistent with the model. This demands linear motions, time-harmonic conditions, homogeneity of the plate and the restriction of horizontal motions of the disc or discs. An optical remote sensing device is employed to record the deflection of the discs accurately. Tests involving a single disc and two discs are conducted for a range of disc thicknesses, incident wave steepnesses, frequencies and, in the case of two discs, geometrical arrangements. A data processing technique is used to decompose the raw data into its spectral harmonics and filter the higher-order components. Pointwise comparisons of the linear first-order component of the experimental deflection with numerical predictions are presented. Satisfying agreement is found, although the model consistently over predicts the deflection. Disc–disc interactions are observed in the two-disc tests. A brief discussion of the shortcomings of the pointwise analysis, with associated possible sources of discrepancy, provides a link to the study reported in Part 2 (Montiel et al. J. Fluid Mech. , vol. 723, 2013, pp. 629–652).

Non-reflecting boundary conditions (NRBCs) play an important role in computational fluid dynamics (CFD). A novel NRBC based on the method of characteristics using timeline interpolations is proposed for fluid dynamics solved by smoothed particle hydrodynamics (SPH). It is performed by four layers of particles whose pressures and velocities are obtained through the Lagrange interpolation in the time domain which is derived from the propagation of characteristic waves between particles. The proposed NRBC can allow outward travelling pressure and velocity messages to pass through the boundary without obvious reflection. That is, with the implementation of the NRBC, the solution in a finite computational domain of interest is close to that in an infinite domain. Several numerical tests show that this NRBC is robust and applicable for a broad variety of hydrodynamics ranging from low to high speed.

The paper presents a study on four adaptive sampling methods of a multi-fidelity (MF) metamodel, based on stochastic radial basis functions (RBF), for global design optimisation based on expensive CFD computer simulations and adaptive grid refinement. The MF metamodel is built as the sum of a low-fidelity-trained metamodel and an error metamodel, based on the difference between high- and low-fidelity simulations. The MF metamodel is adaptively refined using dynamic sampling criteria, based on the prediction uncertainty in combination with the objective optimum and the computational cost of high- and low-fidelity evaluations. The adaptive sampling methods are demonstrated by four analytical benchmark and two design optimisation problems, pertaining to the resistance reduction of a NACA hydrofoil and a destroyer-type vessel. The performance of the adaptive sampling methods is assessed via objective function convergence.

The present paper summarizes analyses of a two-body floating wave energy converter (WEC) to determine the mooring tension and the effect of the mooring system on energy capture. Also, the effect of the power take-off (PTO) is assessed. An axisymmetric Wavebob-type WEC is chosen as the object of investigation. However, the PTO system is modeled in a simplified manner as ideal linear damping and spring terms that couple the motions of the two bodies. The analysis is performed using SIMO, which is a time domain simulation tool that accommodates the simulation of multibody systems with hydrodynamic interactions. In SIMO, docking cone features between the two bodies allow movement as per actual operation, and fenders are applied to represent end stops. Six alternative mooring configurations are applied to investigate the effect of mooring on power capture. Mooring analysis is performed to determine the necessary capacity of mooring lines for each configuration to carry the tension due to the WEC motion in extreme conditions. Hydrodynamic loads are determined using WAMIT. We assumed that the WEC will be operated to capture wave power at the Yeu site in France. The analysis is performed for several regular and irregular wave conditions according to wave data available for that site. Simulations are performed to study the effect of the PTO system, end stops settings and several mooring configurations on power capture.

Abstract Validation of a linear numerical model of wave interactions with floating compliant discs is sought using data obtained from the wave basin experiments reported in Part 1 (Montiel et al. J. Fluid Mech. , vol. 723, 2013, pp. 604–628). Comparisons are made for both single-disc tests and the two-disc tests in which wave interactions between discs are observed. The deflection of the disc or discs is separated into the natural modes of vibration in vacuo . The decomposition allows the rigid-body motions and flexural motions to be analysed separately. Rigid-body motions are accurately replicated by the numerical model but, although passable agreement is found, the amplitudes of flexural modes are consistently overestimated. Extensions of the numerical model are used to discount the experimental configuration as a source of the discrepancies. An enhanced viscoelastic model for the discs is also proposed, which results in improved model/data agreement for the flexural motions but cannot account for all of the disagreement.

International audience

The present work deals with some recent developments regarding the inclusion of the Large-Eddy Simulation (LES) in the weakly compressible Smoothed Particle Hydrodynamics (SPH) framework. Previously {see the work of Di Mascio et al. [Phys. Fluids 29, 4 (2017)]}, this goal was achieved by applying a Lagrangian filter to the Navier–Stokes equations for compressible fluids and, then, approximating the differential operators in a SPH fashion. Since the Lagrangian nature of the derived scheme turned out to be an obstacle for accurate simulations of high Reynolds number problems, the above approach is here modified to obtain a quasi-Lagrangian LES-SPH model. This relies on the addition of a small velocity deviation to the actual Lagrangian velocity based on the particle shifting technique and on the inclusion of the tensile instability control technique for eliminating the onset of the tensile instability in the fluid regions characterized by large vorticity and negative pressure. The proposed model is successfully tested in both two-dimensional and three-dimensional frameworks by simulating the evolution of freely decaying turbulence problems and comparing the outputs with the available theoretical results and solutions from other numerical models.

In the present work, a boundary layer developing over a rough-wall consisting of staggered cubes with a plan area packing density λp = 25% is studied within the wind tunnel using Particle Image Velocimetry (PIV) to investigate the Turbulent Kinetic Energy (TKE) budget. To access the full TKE budget, an estimation of the dissipation (ε) using both the transport equation of the resolved-scale kinetic energy and Large-Eddy (LE) PIV models based on the use of a subgrid-scale model following the methodology used in large-eddy simulations is employed. A low-pass filter, larger than the Taylor microscale, is applied to the data prior to the computation of the velocity gradients ensuring a clear cutoff in the inertial range where the models are valid. The presence of the cube roughness elements has a significant influence on the TKE budget due to the region of strong shear that develops over the cubes. The shear layer is shown to produce and dissipate energy, as well as transport energy through advection, turbulent transport, and pressure transport. The recirculation region that forms through the interaction of the shear layer and the canopy layer, which is the region below the height of the cube roughness, creates rapid longitudinal evolution of the mean flow thereby inducing weak production. Finally, through stochastic estimation of the conditional average, it is shown that localized regions of backscatter (energy transfer from unresolved to resolved scales) and forward scatter (energy transfer from resolved to unresolved scales) occur as a result of coherent vortical structures.