A hybrid stabilization technique for simulating water wave - Structure interaction by incompressible Smoothed Particle Hydrodynamics (ISPH) method

The Smoothed Particle Hydrodynamics (SPH) method is emerging as a potential tool for studying water wave related problems, especially for violent free surface flow and large deformation problems. The incompressible SPH (ISPH) computations have been found not to be able to maintain the stability in certain situations and there exist some spurious oscillations in the pressure time history, which is similar to the weakly compressible SPH (WCSPH). One main cause of this problem is related to the non-uniform and clustered distribution of the moving particles. In order to improve the model performance, the paper proposed an efficient hybrid numerical technique aiming to correct the ill particle distributions. The correction approach is realized through the combination of particle shifting and pressure gradient improvement. The advantages of the proposed hybrid technique in improving ISPH calculations are demonstrated through several applications that include solitary wave impact on a slope or overtopping a seawall, and regular wave slamming on the subface of open-piled structure

A comparative study on violent sloshing with complex baffles using the ISPH method

The Smoothed Particle Hydrodynamics (SPH) method has become one of the most promising methods for violent wave impact simulations. In this paper, the incompressible SPH (ISPH) method will be used to simulate liquid sloshing in a 2D tank with complex baffles. Firstly, the numerical model is validated against the experimental results and the simulations from commercial CFD software STAR-CCM+ for a sloshing tank without any baffle. Then various sloshing tanks are simulated under different conditions to analyze the influence of the excitation frequency and baffle configuration. The results show that the complex baffles can significantly influence the impact pressures on the wall caused by the violent sloshing, and the relevant analysis can help find the engineering solutions to effectively suppress the problem. The main purpose of the paper is to study the practical importance of this effect

An improved solid boundary treatment for wave-float interactions using ISPH method

The Smoothed Particle Hydrodynamics (SPH) method has proved to have great potentials in dealing with the wave-structure interactions. Compared with the Weakly Compressible SPH (WCSPH) method, the ISPH approach solves the pressure by using the pressure Poisson equation rather than the equation of state. This could provide a more stable and accurate pressure field that is important in the study of wave-structure interactions. This paper improves the solid boundary treatment of ISPH by using a high accuracy Simplified Finite Difference Interpolation (SFDI) scheme for the 2D wave-structure coupling problems, especially for free-moving structure. The proposed method is referred as the ISPH_BS. The model improvement is demonstrated by the documented benchmark tests and laboratory experiment covering various wave-structure interaction applications

A review on approaches to solving Poisson’s equation in projection-based meshless methods for modelling strongly nonlinear water waves

Three meshless methods, including incompressible smooth particle hydrodynamic (ISPH), moving particle semi-implicit (MPS) and meshless local Petrov–Galerkin method based on Rankine source solution (MLPG_R) methods, are often employed to model nonlinear or violent water waves and their interaction with marine structures. They are all based on the projection procedure, in which solving Poisson’s equation about pressure at each time step is a major task. There are three different approaches to solving Poisson’s equation, i.e. (1) discretizing Laplacian directly by approximating the second-order derivatives, (2) transferring Poisson’s equation into a weak form containing only gradient of pressure and (3) transferring Poisson’s equation into a weak form that does not contain any derivatives of functions to be solved. The first approach is often adopted in ISPH and MPS, while the third one is implemented by the MLPG_R method. This paper attempts to review the most popular, though not all, approaches available in literature for solving the equation

A Consistent Second Order ISPH for Free Surface Flow

The Incompressible Smoothed Particle Hydrodynamics (ISPH) is now a popular numerical method for modelling free surface flows, in particular the breaking waves and violent wave-structures interaction. The ISPH requires the projection approach, leading to solving a pressure Poisson's equation (PPE). Although the accuracy and convergence of the numerical scheme to discretise the Laplacian operator involved in PPE is critical for securing a satisfactory solution of the PPE, the overall performance of the ISPH is also influenced by other key numerical implementations, including (1) estimation of the viscous terms; (2) calculation of the velocity divergence; (3) discretisation of the boundary conditions for the PPE; and (4) evaluation of the pressure gradient. In our previous paper [29], the quadratic semi-analytical finite difference interpolation scheme (QSFDI), which has a leading truncation error at third order derivatives, has been adopted to discretise the Laplacian operator. In this paper, the QSFDI will be adopted, not only for discretising the Laplacian operator, but also for approximating viscous terms, velocity divergence, boundary conditions and pressure gradient. The performance of the newly formulated consistent second order ISPH is assessed by various cases including the oscillating liquid drop, the wave propagation, and the liquid sloshing. The results do not only demonstrate a second order convergence over a limited range of conditions and a higher computational efficiency, i.e., requiring less computational time to achieve the same accuracy, but also show a better mass/energy conservation property and capacity of reproducing a smooth pressure field, than other ISPH models considered in this study

Time estimation and beta segregation: An EEG study and graph theoretical approach

Elucidation of the neural correlates of time perception constitutes an important research topic in cognitive neuroscience. The focus to date has been on durations in the millisecond to seconds range, but here we used electroencephalography (EEG) to examine brain functional connectivity during much longer durations (i.e., 15 min). For this purpose, we conducted an initial exploratory experiment followed by a confirmatory experiment. Our results showed that those participants who overestimated time exhibited lower activity of beta (18�30 Hz) at several electrode sites. Furthermore, graph theoretical analysis indicated significant differences in the beta range (15�30 Hz) between those that overestimated and underestimated time. Participants who underestimated time showed higher clustering coefficient compared to those that overestimated time. We discuss our results in terms of two aspects. FFT results, as a linear approach, are discussed within localized/dedicated models (i.e., scalar timing model). Second, non-localized properties of psychological interval timing (as emphasized by intrinsic models) are addressed and discussed based on results derived from graph theory. Results suggested that although beta amplitude in central regions (related to activity of BG-thalamocortical pathway as a dedicated module) is important in relation to timing mechanisms, the properties of functional activity of brain networks; such as the segregation of beta network, are also crucial for time perception. These results may suggest subjective time may be created by vector units instead of scalar ticks. © 2018 Ghaderi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Time estimation and beta segregation: An EEG study and graph theoretical approach

Elucidation of the neural correlates of time perception constitutes an important research topic in cognitive neuroscience. The focus to date has been on durations in the millisecond to seconds range, but here we used electroencephalography (EEG) to examine brain functional connectivity during much longer durations (i.e., 15 min). For this purpose, we conducted an initial exploratory experiment followed by a confirmatory experiment. Our results showed that those participants who overestimated time exhibited lower activity of beta (18�30 Hz) at several electrode sites. Furthermore, graph theoretical analysis indicated significant differences in the beta range (15�30 Hz) between those that overestimated and underestimated time. Participants who underestimated time showed higher clustering coefficient compared to those that overestimated time. We discuss our results in terms of two aspects. FFT results, as a linear approach, are discussed within localized/dedicated models (i.e., scalar timing model). Second, non-localized properties of psychological interval timing (as emphasized by intrinsic models) are addressed and discussed based on results derived from graph theory. Results suggested that although beta amplitude in central regions (related to activity of BG-thalamocortical pathway as a dedicated module) is important in relation to timing mechanisms, the properties of functional activity of brain networks; such as the segregation of beta network, are also crucial for time perception. These results may suggest subjective time may be created by vector units instead of scalar ticks. © 2018 Ghaderi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Corrected First-order Derivative ISPH in Water Wave Simulations

The smoothed particle hydrodynamics (SPH) method is a meshless numerical modeling technique. It has been applied in many different research fields in coastal engineering. Due to the drawback of its kernel approximation, however, the accuracy of SPH simulation results still needs to be improved in the prediction of violent wave impact. This paper compares several different forms of correction on the first-order derivative of ISPH formulation aiming to find one optimum kernel approximation. Based on four benchmark case analysis, we explored different kernel corrections and compared their accuracies. Furthermore, we applied them to one solitary wave and two dam-break flows with violent wave impact. The recommended method has been found to achieve much more promising results as compared with experimental data and other numerical approaches

A hybrid method for modelling wake flow of a wind turbine

A fast and accurate prediction of the wake of an upwind turbine is very important to quantify the performance of downwind turbines in offshore wind farms, which become larger and larger. The wake flow and dynamics may be quite accurately simulated by high-fidelity computational fluid dynamics (CFD) software but its computational costs are too high, in particular to simulate a long wake flow often required in engineering practice. Therefore, the wake is often modelled by simplified dynamic wake models in design practices. They are computationally efficient but could not catch all physics, depend on pre-specified empirical parameters, and are not suitable for flow near the turbine. This paper proposes a new hybrid method, in which the near wake flow is simulated by a CFD model based on Navier–Stokes equations with the turbine represented by actuator lines while the far wake flow is modelled by an improved simplified CFD-based dynamic wake model. The two models are two-way coupled at a section downwind the turbine. The newly formulated method is validated by the results of full CFD simulations in the whole domain. Its performances are investigated under different conditions. It will be demonstrated that the new method takes considerably less computational time than the full CFD tool to produce similar results

A set of canonical problems in sloshing. Part 0: Experimental setup and data processing

In a series of attempts to research and document relevant sloshing type phenomena, a series of experiments have been conducted. The aim of this paper is to describe the setup and data processing of such experiments. A sloshing tank is subjected to angular motion. As a result pressure registers are obtained at several locations, together with the motion data, torque and a collection of image and video information. The experimental rig and the data acquisition systems are described. Useful information for experimental sloshing research practitioners is provided. This information is related to the liquids used in the experiments, the dying techniques, tank building processes, synchronization of acquisition systems, etc. A new procedure for reconstructing experimental data, that takes into account experimental uncertainties, is presented. This procedure is based on a least squares spline approximation of the data. Based on a deterministic approach to the first sloshing wave impact event in a sloshing experiment, an uncertainty analysis procedure of the associated first pressure peak value is described