scholarly journals Implementation of Open Boundaries within a Two-Way Coupled SPH Model to Simulate Nonlinear Wave–Structure Interactions

Energies ◽  
2019 ◽  
Vol 12 (4) ◽  
pp. 697 ◽  
Author(s):  
Tim Verbrugghe ◽  
Vasiliki Stratigaki ◽  
Corrado Altomare ◽  
J. Domínguez ◽  
Peter Troch ◽  
...  

A two-way coupling between the Smoothed Particle Hydrodynamics (SPH) solver DualSPHysics and the Fully Nonlinear Potential Flow solver OceanWave3D is presented. At the coupling interfaces within the SPH numerical domain, an open boundary formulation is applied. An inlet and outlet zone are filled with buffer particles. At the inlet, horizontal orbital velocities and surface elevations calculated using OceanWave3D are imposed on the buffer particles. At the outlet, horizontal orbital velocities are imposed, but the surface elevation is extrapolated from the fluid domain. Velocity corrections are applied to avoid unwanted reflections in the SPH fluid domain. The SPH surface elevation is coupled back to OceanWave3D, where the originally calculated free surface is overwritten. The coupling methodology is validated using a 2D test case of a floating box. Additionally, a 3D proof of concept is shown where overtopping waves are acting on a heaving cylinder. The two-way coupled model (exchange of information in two directions between the coupled models) has proven to be capable of simulating wave propagation and wave–structure interaction problems with an acceptable accuracy with error values remaining below the smoothing length h S P H .

Author(s):  
Tim Verbrugghe ◽  
José Manuel Dominguez ◽  
Corrado Altomare ◽  
Angelantonio Tafuni ◽  
Peter Troch ◽  
...  

A two-way coupling between the fully non-linear potential flow (FNPF) solver OceanWave3D and the Smoothed Particle Hydrodynamics (SPH) solver DualSPHysics is presented. At the coupling interfaces within the SPH domain, an open boundary formulation is applied. An inlet and outlet zone are filled with bu er particles. At the inlet, horizontal orbital velocities and surface elevations calculated with OceanWave3D are imposed on the bu er particles. At the outlet, horizontal orbital velocities are imposed, but the surface elevation is extrapolated from the fluid domain. Velocity corrections are applied to avoid unwanted reflections in the fluid domain. The SPH surface elevation can be coupled back to OceanWave3D, where the original solution is overwritten. The coupling methodology is validated using a 2-D test case of a floating box. Additionally, a 3-D proof of concept is shown where overtopping waves are acting on a heaving cylinder. The 2-way coupled model proofs to be capable of simulating wave propagation and wave-structure interaction problems with an acceptable accuracy with RMSE values remaining below the smoothing length h.


Author(s):  
Pierre-Yves Wuillaume ◽  
Pierre Ferrant ◽  
Aurélien Babarit ◽  
François Rongère ◽  
Mattias Lynch ◽  
...  

This paper presents validation tests for a new numerical tool for the numerical simulation of marine operations. It involves multibody dynamics modeling, wave-structure interactions with large amplitude body motion and cable’s dynamic modeling. Hydrodynamic loads are computed using the WS_CN weakly nonlinear potential flow solver, based on the weak-scatterer hypothesis. Large deformation of the wetted body surfaces can be taken into account. Firstly the ECN’s WS_CN solver capabilities are extended to multibody simulations. A first validation test is performed by comparing numerical results to the experimental data of [1]. Then, a second validation test is proposed. It consists in the ballasting operation of a spar. The experimental set-up is described.


Author(s):  
Bo Terp Paulsen ◽  
Henrik Bredmose ◽  
Harry B. Bingham ◽  
Signe Schløer

Two-dimensional irregular waves on a sloping bed and their impact on a bottom mounted circular cylinder is modeled by three different numerical methods and the results are validated against laboratory experiments. We here consider the performance of a linear-, a fully nonlinear potential flow solver and a fully nonlinear Navier-Stokes/VOF solver. The validation is carried out in terms of both the free surface elevation and the inline force. Special attention is paid to the ultimate load in case of a single wave event and the general ability of the numerical models to capture the higher harmonic forcing. The test case is representative for monopile foundations at intermediate water depths. The potential flow computations are carried out in a two-dimensional vertical plane and the inline force on the cylinder is evaluated by the Morison equation. The Navier-Stokes/VOF computations are carried out in three-dimensions and the force is obtained by spatial pressure integration over the wettet area of the cylinder. In terms of both the free surface elevation and the inline force, the linear potential flow model is shown to be of limited accuracy and large deviations are generally seen when compared to the experimental measurements. The fully nonlinear Navier-Stokes/VOF computations are accurately predicting both the free surface elevation and the inline force. However, the computational cost is high relative to the potential flow solvers. Despite the fact that the nonlinear potential flow model is carried out in two-dimensions it is shown to perform just as good as the three-dimensional Navier-Stokes/VOF solver. This is observed for both the free surface elevation and the inline force, where both the ultimate load and the higher harmonic forces are accurately predicted. This shows that for moderately steep irregular waves a Morison equation combined with a fully nonlinear two-dimensional potential flow solver can be a good approximation.


Water ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 429
Author(s):  
Gael Verao Fernandez ◽  
Vasiliki Stratigaki ◽  
Nicolas Quartier ◽  
Peter Troch

The study of the potential impact of wave energy converter (WEC) farms on the surrounding wave field at long distances from the WEC farm location (also know as “far field” effects) has been a topic of great interest in the past decade. Typically, “far-field” effects have been studied using phase average or phase resolving numerical models using a parametrization of the WEC power absorption using wave transmission coefficients. Most recent studies have focused on using coupled models between a wave-structure interaction solver and a wave-propagation model, which offer a more complex and accurate representation of the WEC hydrodynamics and PTO behaviour. The difference in the results between the two aforementioned approaches has not been studied yet, nor how different ways of modelling the PTO system can affect wave propagation in the lee of the WEC farm. The Coastal Engineering Research Group of Ghent University has developed both a parameterized model using the sponge layer technique in the mild slope wave propagation model MILDwave and a coupled model MILDwave-NEMOH (NEMOH is a boundary element method-based wave-structure interaction solver), for studying the “far-field” effects of WEC farms. The objective of the present study is to perform a comparison between both numerical approaches in terms of performance for obtaining the “far-field” effects of two WEC farms. Results are given for a series of regular wave conditions, demonstrating a better accuracy of the MILDwave-NEMOH coupled model in obtaining the wave disturbance coefficient (Kd) values around the considered WEC farms. Subsequently, the analysis is extended to study the influence of the PTO system modelling technique on the “far-field” effects by considering: (i) a linear optimal, (ii) a linear sub-optimal and (iii) a non-linear hydraulic PTO system. It is shown that modelling a linear optimal PTO system can lead to an unrealistic overestimation of the WEC motions than can heavily affect the wave height at a large distance in the lee of the WEC farm. On the contrary, modelling of a sub-optimal PTO system and of a hydraulic PTO system leads to a similar, yet reduced impact on the “far-field” effects on wave height. The comparison of the PTO systems’ modelling technique shows that when using coupled models, it is necessary to carefully model the WEC hydrodynamics and PTO behaviour as they can introduce substantial inaccuracies into the WECs’ motions and the WEC farm “far-field” effects.


2021 ◽  
Author(s):  
Vito Zago ◽  
Giuseppe Bilotta ◽  
Annalisa Cappello ◽  
Robert Dalrymple ◽  
Luigi Fortuna ◽  
...  

<p>Numerical simulation is a fundamental aspect of modern volcanology, providing tools for the forecasting of lava flows behavior, so as to assist in the design of mitigation actions for volcanic risk. In addition to the prediction of the emplacement topology, numerical simulation can be useful to study the possible outcomes of the interaction between a lava flow and a building. This kind of information can help to estimate the vulnerability of buildings so as to produce more accurate risk evaluations. Smoothed Particle Hydrodynamics (SPH) is a particle-based numerical method, particularly suited for the simulation of fluids with a high level of complexity, that can intrinsically deal with all of the physical properties of lava. GPUSPH is a simulation engine based on the SPH method that has been developed in order to take into account the challenging aspects of lava simulations and has been successfully applied to the simulation of lava-related benchmark tests. Here we use the SPH method, coupled within the framework of GPUSPH with a rigid body mechanics solver provided by the Project Chrono engine, for the realistic study of lava-buildings interaction. The resulting coupled model is able to simulate masonry with a brick-level accurate description, providing insights on any damages happening to the structure. We will show the simulation of a lava flow interacting with an elementary masonry piece, where a total collapse of the structure is induced by the action of the lava.</p>


2012 ◽  
Vol 11 (1-2) ◽  
pp. 68 ◽  
Author(s):  
E. Didier ◽  
D. R. C. B. Neves ◽  
R. Martins ◽  
M. G. Neves

This work presents the new developments and the validation of a Smoothed Particle Hydrodynamics (SPH) numerical model used in the National Laboratory of Civil Engineering (Laboratório Nacional de Engenharia Civil - LNEC) for studies in coastal engineering processes. Although the model requires a high CPU time, it proved to be very promising in the simulation of complex flows, such as the wave-structure interaction and the wave breaking phenomenon. For the SPH model validation, physical modeling tests were performed in one LNEC’s flume to study the interaction between an impermeable structure and an incident regular wave. The comparison between numerical and experimental results, i.e. free surface elevation, overtopping volume and pressure, shows the good accuracy of the SPH model to reproduce the various phenomena involving on the wave propagation and interaction with the structure, namely the wave breaking, the wave overtopping and the pressure field on the structure.


Author(s):  
Steven J. Lind ◽  
Benedict D. Rogers ◽  
Peter K. Stansby

This paper presents a review of the progress of smoothed particle hydrodynamics (SPH) towards high-order converged simulations. As a mesh-free Lagrangian method suitable for complex flows with interfaces and multiple phases, SPH has developed considerably in the past decade. While original applications were in astrophysics, early engineering applications showed the versatility and robustness of the method without emphasis on accuracy and convergence. The early method was of weakly compressible form resulting in noisy pressures due to spurious pressure waves. This was effectively removed in the incompressible (divergence-free) form which followed; since then the weakly compressible form has been advanced, reducing pressure noise. Now numerical convergence studies are standard. While the method is computationally demanding on conventional processors, it is well suited to parallel processing on massively parallel computing and graphics processing units. Applications are diverse and encompass wave–structure interaction, geophysical flows due to landslides, nuclear sludge flows, welding, gearbox flows and many others. In the state of the art, convergence is typically between the first- and second-order theoretical limits. Recent advances are improving convergence to fourth order (and higher) and these will also be outlined. This can be necessary to resolve multi-scale aspects of turbulent flow.


Energies ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3745
Author(s):  
Tristan Revaz ◽  
Fernando Porté-Agel

Large-eddy simulation (LES) with actuator models has become the state-of-the-art numerical tool to study the complex interaction between the atmospheric boundary layer (ABL) and wind turbines. In this paper, a new evaluation of actuator disk models (ADMs) for LES of wind turbine flows is presented. Several details of the implementation of such models are evaluated based on a test case studied experimentally. In contrast to other test cases used in previous similar studies, the present test case consists of a wind turbine immersed in a realistic turbulent boundary-layer flow, for which accurate data for the turbine, the flow, the thrust and the power are available. It is found that the projection of the forces generated by the turbine into the flow solver grid is crucial for rotor predictions, especially for the power, and less important for the wake flow prediction. In this context, the projection of the forces into the flow solver grid should be as accurate as possible, in order to conserve the consistency between the computed axial velocity and the projected axial force. Also, the projection of the force is found to be much more important in the rotor plane directions than in the streamwise direction. It is found that for the case of a wind turbine immersed in a realistic turbulent boundary-layer flow, the potential spurious numerical oscillations originating from sharp force projections are not harmful to the results. By comparing an advanced model which computes the non-uniform distribution of the turbine forces over the rotor with a simple model which assumes uniform effects of the turbine forces, it is found that both can lead to accurate results for the far wake flow and the thrust and power predictions. However, the comparison shows that the advanced model leads to better results for the near wake flow. In addition, it is found that the simple model overestimates the rotor velocity prediction in comparison to the advanced model. These elements are explained by the lack of local feedback between the axial velocity and the axial force in the simple model. By comparing simulations with and without including the effects of the nacelle and tower, it is found that the consideration of the nacelle and tower is relatively important both for the near wake and the power prediction, due to the shadow effects. The grid resolution is not found to be critical once a reasonable resolution is used, i.e. in the order of 10 grid points along each direction across the rotor. The comparison with the experimental data shows that an accurate prediction of the flow, thrust, and power is possible with a very reasonable computational cost. Overall, the results give important guidelines for the implementation of ADMs for LES.


Water ◽  
2019 ◽  
Vol 11 (11) ◽  
pp. 2314 ◽  
Author(s):  
Shu Wang ◽  
Anping Shu ◽  
Matteo Rubinato ◽  
Mengyao Wang ◽  
Jiping Qin

Non-homogeneous viscous debris flows are characterized by high density, impact force and destructiveness, and the complexity of the materials they are made of. This has always made these flows challenging to simulate numerically, and to reproduce experimentally debris flow processes. In this study, the formation-movement process of non-homogeneous debris flow under three different soil configurations was simulated numerically by modifying the formulation of collision, friction, and yield stresses for the existing Smoothed Particle Hydrodynamics (SPH) method. The results obtained by applying this modification to the SPH model clearly demonstrated that the configuration where fine and coarse particles are fully mixed, with no specific layering, produces more fluctuations and instability of the debris flow. The kinetic and potential energies of the fluctuating particles calculated for each scenario have been shown to be affected by the water content by focusing on small local areas. Therefore, this study provides a better understanding and new insights regarding intermittent debris flows, and explains the impact of the water content on their formation and movement processes.


2005 ◽  
Vol 18 (13) ◽  
pp. 2441-2459 ◽  
Author(s):  
J. Zavala-Garay ◽  
C. Zhang ◽  
A. M. Moore ◽  
R. Kleeman

Abstract The possibility that the tropical Pacific coupled system linearly amplifies perturbations produced by the Madden–Julian oscillation (MJO) is explored. This requires an estimate of the low-frequency tail of the MJO. Using 23 yr of NCEP–NCAR reanalyses of surface wind and Reynolds SST, we show that the spatial structure that dominates the intraseasonal band (i.e., the MJO) also dominates the low-frequency band once the anomalies directly related to ENSO have been removed. This low-frequency contribution of the intraseasonal variability is not included in most ENSO coupled models used to date. Its effect in a coupled model of intermediate complexity has, therefore, been studied. It is found that this “MJO forcing” (τMJO) can explain a large fraction of the interannual variability in an asymptotically stable version of the model. This interaction is achieved via linear dynamics. That is, it is the cumulative effect of individual events that maintains ENSOs in this model. The largest coupled wind anomalies are initiated after a sequence of several downwelling Kelvin waves of the same sign have been forced by τMJO. The cumulative effect of the forced Kelvin waves is to persist the (small) SST anomalies in the eastern Pacific just enough for the coupled ocean–atmosphere dynamics to amplify the anomalies into a mature ENSO event. Even though τMJO explains just a small fraction of the energy contained in the stress not associated with ENSO, a large fraction of the modeled ENSO variability is excited by this forcing. The characteristics that make τMJO an optimal stochastic forcing for the model are discussed. The large zonal extent is an important factor that differentiates the MJO from other sources of stochastic forcing.


Sign in / Sign up

Export Citation Format

Share Document