Search

Numerical simulation of fluid structure interaction (FSI) problems is one of the most challenging topics in computational fluid dynamics. In this paper, coupling edge-based smoothed finite element method (ES-FEM) and smoothed particle hydrodynamics (SPH) method (ES-FEM-SPH) is proposed for solving FSI problems, where the edge-based smoothed finite element method is used to model the movement and deformation of structures, and the smoothed particle hydrodynamics is used to model the fluid flow. In ES-FEM, the gradient smoothing technique is applied over the smoothing domain and it can effectively overcome the "overly-stiff" effect in conventional FEM model. Some correction algorithms including density correction, kernel gradient correction and particle shift technique are integrated into the SPH method to improve computational stability and accuracy. A virtual particle coupling scheme is used to implement the coupling of ES-FEM and SPH with complex geometry interface. As ES-FEM is more accurate than conventional FEM, and it is expected that this ES-FEM-SPH coupling approach should be superior than existing FEM-SPH coupling approaches. A number of test examples with FSI are investigated with the presented ES-FEM-SPH, and compared with results from other approaches including FEM-SPH. From the obtained numerical results, we can conclude that the ES-FEM-SPH coupling approach is effective to simulate FSI problems. • Coupling edge-based smoothed finite element method (ES-FEM) with SPH is proposed for solving FSI problems. • The coupling strategy is based on the automatically generated virtual particles. • The gradient smoothing technique in ES-FEM overcomes the "overly-stiff" effect in conventional FEM. • The SPH method with integrated algorithms is effective and robust in modeling fluid flows. • ES-FEM-SPH is demonstrated to be effective in modeling diversified FSI problems. [ABSTRACT FROM AUTHOR]

• Four integrated SPH models are presented and compared for simulating natural convection. • The most suitable SPH model for natural convection is provided. • Natural convection for R a = 10 9 and Pr = 0.71 is successfully modeled for the first time by using the SPH method. • Mechanisms of the natural convection in a square cavity at different Rayleigh numbers are discussed. This paper conducts the simulation of natural convection in a differentially heated square cavity at high Rayleigh numbers by using the smoothed particle hydrodynamics (SPH) method. Due to the decrease of the accuracy and stability, it is challenging for the SPH method to simulate natural convection at high Rayleigh numbers, and there are few reported SPH literatures of natural convection at R a > 10 6 for air (Pr = 0.71). In this study, four integrated SPH models are presented to simulate the natural convection and their accuracy and stability are assessed. These four SPH models are associated with Kernel Gradient Correction (KGC) to improve approximation accuracy and Particle Shifting Technology (PST) to regularize particle distribution while they are different in treating density diffusion and calculating the pressure term. The numerical results show that SPH model_4 (KGC, PST, δ -SPH and asymmetric pressure approximation) is the most suitable for simulating the closed natural convection problems, especially at high Rayleigh numbers. Good agreements with reference solutions are obtained by SPH model_4 for the natural convection at 10 4 ≤ R a ≤ 10 8. Furthermore, the simulation of natural convection at R a = 10 9 is conducted by SPH model_4. The evolutions of thermal convection are described in detail. It is found that dynamics characteristic reveals that the dominant force is the pressure gradient, rather than the buoyancy force before the quasi-steady state. In addition, the chaotic motion at R a = 10 9 has significant influence to the heat transfer characteristic in the vertical boundary layers. [ABSTRACT FROM AUTHOR]

In order to understand the protective mechanism of helmet against shock wave, a head model under far-field blast is investigated experimentally and numerically. The distribution of overpressure on the surface of the head model are obtained and compared. The numerical simulations are carried out using the Multi-Material Arbitrary Lagrangian Eulerian (MM-ALE) technique with suitable material models for a Kevlar helmet and soft pads. The differences in the distributions of overpressure on the head models with and without the protection of the helmet are analyzed. The interaction of shock waves around the head and the helmet is investigated. Our experimental observations and numerical results indicate that, if protected by a helmet in the case of a shock wave acting from the front, the peaks of overpressure on the head are reduced at the top area, but increased at the rear area. The numerical simulations also show that the enhancement of the overpressure at the rear is mainly due to the collision of the shock wave. The underwash phenomenon also occurs when the shock wave acts from the side or the rear. [ABSTRACT FROM AUTHOR]

Predicting the propagation and transformation of the solitary wave is very important in coastal engineering. This paper presents a numerical investigation of the solitary wave breaking over a slope by using a finite particle method (FPM), which is an enhanced smoothed particle hydrodynamics (SPH). We firstly conduct a comparative study on SPH and FPM in modelling solitary wave and it is demonstrated that FPM performs better than SPH qualitatively and quantitatively. A modified particle shifting technique (PST) is integrated into FPM to avoid possible ill corrective matrix due to extremely disordered particle distribution. We find that the artificial viscous coefficient can greatly influence the wave run-up and propose an empirical equation to quickly determine the optimal value of the artificial viscous coefficient used in the FPM simulations. The solitary wave breaking is then modelled for scenarios with relative wave heights and slopes, and three typical breaking types, including surging breaker (SU), plunging breaker (PL) and spilling breaker (SP), are analyzed. The result indicates that SP affects the pressure field more extensive than PL and SU. When the breaker type transits from SP to SU, the breaker tends to break at a narrower region close to the shoreline, the breaking depth decreases, the breaking index gradually turns into an unstable value, and the relative run-up height of breaking waves becomes larger. Although the wave celerity at breaking has no clear relationship with breaker types, the wave celerity at breaking increases with the increase of the relative wave height and slightly decreases with a steeper slope. • For simulations of the solitary wave breaking over slopes, FPM performs better than SPH in terms of stability and accuracy. • An empirical equation is proposed to decide the artificial viscous coefficient for simulating wave transformation by FPM. • The spilling breaker affects the pressure field more extensive than plunging breaker and surging breaker. • The wave celerity at breaking has no clear relationship with breaker types. • The wave dissipation turns into lower when wave breaking types change from spilling breaker to surging breaker. [ABSTRACT FROM AUTHOR]