Search results

The purpose of this paper is to present an efficient smoothed particle hydrodynamics (SPH) method particularly adapted for the geometrically nonlinear analysis of structures.

In order to resolve the inconsistency phenomenon which systematically occurs in the standard SPH method at the domain’s boundaries of the studied structure, the classical kernel function and its spatial derivatives were modified by the use of Taylor series expansion. The well-known tensile instabilities inherent to the Eulerian SPH formulation were attenuated by the use of the Total Lagrangian Formulation (TLF).

In order to demonstrate the effectiveness of the present improved SPH method, several numerical applications involving geometrically nonlinear behaviors were carried out using the explicit dynamics scheme for the time integration of the PDEs. Comparisons of the obtained results using the present SPH model with analytical reference solutions and with those obtained using ABAQUS finite element (FE) commercial software, show its good accuracy and robustness.

An additional application including a multilayered composite structure and involving buckling and delamination was investigated using the present improved SPH model and the results are compared to the FE results, they confirmed both the efficiency and the accuracy of the proposed method.

An efficient 2D-continuum SPH model for the geometrically nonlinear analysis of thin and thick structures is proposed. Contrarily to the classical SPH approaches, here the constitutive material relations are used to link naturally the stresses and strains. The Total Lagrangian approach is investigated to alleviate the tensile instabilities problem, allowing at the same time to avoid the updating procedure of the neighboring particles search and therefore reducing CPU usage. The proposed approach is valid for isotropic and multilayered composites structures undergoing large transformations. CPU time savings and better results with the new 2D-continuum SPH formulation compared to the classical continuum SPH. The explicit dynamic scheme was used for time integration allowing a fast resolution algorithm even for highly nonlinear problems.

This paper aims to propose a new smoothed particle hydrodynamics (SPH)-finite element (FE) algorithm to study fluid–structure interaction (FSI) problems.

The fluid domain is discretized based on the theory of SPH), and solid part is solved through FE method, similar to other SPH-FE methods in the previous studies. Instead of master-slave technique, the interpolating (kernel) functions of immersed boundary method are implemented to couple fluid and solid domains. The procedure of modeling completely follows the classic IB framework where forces and velocities are transferred between interacting parts. Three benchmark FSI problems are simulated and the results are compared with those of similar numerical and experimental works.

The proposed SPH-FE algorithm with promising and acceptable results can be utilized as a reliable method to simulate FSI problems.

Contrary to most SPH-FE algorithms, the calculation of contact force is not required at interacting boundaries and no iterative process is proposed to calculate forces, velocities and positions at new time step.

The purpose of this paper is to investigate the formation process of linear‐shaped charge jet using the smoothed particle hydrodynamics (SPH). Different material yield models are embed to test the performance of SPH method in the simulation of explosive driven metal liner. The effects of different ignition model to the formation of metal jet have also been studied.

The SPH method is used with the correction of artificial viscosity and penalty force to simulate the formation process of linear‐shaped charge jet, which includes the process of explosion and interaction between explosive gas and metal liner. The numerical results which got by SPH method are compared with these obtained by mesh‐based method. Different material yield models are implemented in the numerical examples to show the effect of material model to the formation process of metal jet. The single point and two point ignition models are used to study the effect of ignition models to the process of explosion and formation of metal jet.

Compared with the original mesh‐based method, the SPH method can simulate the physical process of linear‐shaped charge jet naturally, as well as the capturing of explosive wave propagation. The implementation of different material yields models to obtain the same formation tendency of metal jet, but some numerical difference exists. In two‐point ignition model the explosive pressure is superimposed at the location that two detonation waves intersect. Compared with two ignition models, the two point ignition model can form the metal jet faster and get the higher velocity metal jet.

There are a few references that address the application of SPH to simulate shaped charge explosion process. The feasibility of the SPH method to simulate the formation process of linear shaped charged jet is tested and verified in this paper. From the results which compared with mesh‐based method, it is shown that the SPH method has the advantage in tracking the large deformation of material and capturing the explosive wave propagation. The SPH method can be selected as a good alternative to traditional mesh‐based numerical methods in simulating similar explosively driven metal material problems which can be referenced from this paper.

This paper aims to introduce a two-dimensional smoothed particle hydrodynamics (SPH) framework for simulating the evaporation and combustion process of fuel droplets.

To solve the gas–liquid two-phase flow problem, a multiphase SPH method capable of handling high density-ratio problems is established. Based on the Fourier heat conduction equation and Fick’s law of diffusion, the SPH discrete equations are derived. To effectively characterize the phase transition problem, inspired by volume of fluid method, the concept of liquid phase mass fraction of the SPH particles is proposed. The one-step global reaction model of n-hexane is used for the vapor combustion.

The evaporation and combustion process of single droplet conforms to the law. The framework works out well when the evaporation of multiple droplets involves coalescence process. Three different kinds of flames are observed in succession in the combustion process of a single droplet at different inflow velocity, which agree well with the results of the experiment.

To the best of the authors’ knowledge, this is the first computational framework that has the capability to simulate evaporation and combustion with SPH method. Based on the particle nature of SPH method, the framework has natural advantages in interface tracking.

The paper aims to construct a method to simulate the relationship between the parameters of soil properties and the area of sliding mass of the true slip surface of a landslide.

The smoothed particle hydrodynamics (SPH) algorithm is used to calibrate a response surface function which is adopted to quantify the area of sliding mass of the true slip surface for each failure sample in Monte Carlo simulation. The proposed method is illustrated through a homogeneous and a heterogeneous cohesive soil slope.

The comparison of the results between the proposed method and the traditional method using the slip surface with minimum factor of safety (FS_min) to quantify the failure consequence has shown that the landslide risk tends to be attributed to a variety of risk sources, and that the use of a slip surface with FS_min to quantify the consequence of a landslide underestimates the landslide risk value. The difference of the risk value between the proposed method and the traditional method increases dramatically as the uncertainty of soil properties becomes significant.

A geotechnical engineer could use the proposed method to perform slope failure analysis.

The failure consequence of a landslide can be rationally predicted using the proposed method.

A new coupled model is developed to simulate the interaction between fluid droplet collisions on discrete particles (DPs) by using mathematic function.

In this model, the smoothed particle hydrodynamics (SPH) is used based on the kernel function and the time step which takes into consideration to the fluid domain in accordance with the discrete element method (DEM) with resistance function. The interaction between fluid and DPs consists of three parts, which are repulsive force, viscous shear force and attractive force caused by the capillary action. The numerical simulation of droplet collision on DPs presents the whole process of droplet motion. Otherwise, an experimental data were conducted to record the realistic process for verification.

The comparison result indicated that the numerical simulation is capable of capturing the entire process for droplet collision on DPs.

However, based on the difference of experimental environment, type of the DP and setups, the maximum spreading dimeters of could not fit the experimental data exactly.

In sum, the coupled SPH-DEM method simulation shows that the coupled model of SPH-DEM developed an entire effectiveness process for fluid–solid interaction problem.

Modeling of multi-phase flows for Rayleigh-Taylor instability and natural convection in a square cavity has been investigated using an incompressible smoothed particle hydrodynamics (ISPH) technique. In this technique, incompressibility is enforced by using SPH projection method and a stabilized incompressible SPH method by relaxing the density invariance condition is applied. The paper aims to discuss these issues.

The Rayleigh-Taylor instability is introduced in two and three phases by using ISPH method. The author simulated natural convection in a square/cubic cavity using ISPH method in two and three dimensions. The solutions represented in temperature, vertical velocity and horizontal velocity have been studied with different values of Rayleigh number Ra parameter (103=Ra=105). In addition, characteristic based scheme in Finite Element Method is introduced for modeling the natural convection in a square cavity.

The results for Rayleigh-Taylor instability and natural convection flow had been compared with the previous researches.

Modeling of multi-phase flows for Rayleigh-Taylor instability and natural convection in a square cavity has been investigated using an ISPH technique. In ISPH method, incompressibility is enforced by using SPH projection method and a stabilized incompressible SPH method by relaxing the density invariance condition is introduced. The Rayleigh-Taylor instability is introduced in two and three phases by using ISPH method. The author simulated natural convection in a square/cubic cavity using ISPH method in two and three dimensions.

The purpose of this paper is to confirm that the axisymmetric finite element and smoothed particle hydrodynamics (FE-SPH) adaptive coupling method is effective to solve explosion problem in concrete based on the experiments.

Axisymmetric FE-SPH adaptive coupling method is first presented to simulate dynamic deformation process of concrete under internal blast loading. Using calculation codes of FE-SPH coupling method, numerical model of explosion is approximated initially by finite element method (FEM), and distorted finite elements are automatically converted into meshless particles to simulate damage, splash of concrete by SPH method, when equivalent plastic strain of elements reaches a specified value.

In this paper, damage process and pressure curve of concrete around explosive are analyzed and buried depth of explosive in concrete influence on damage effect under internal blast loading are obtained. Numerical analyses show that FE-SPH coupling method integrates high computational efficiency of FEM and advantages of SPH method, such as natural simulation to damage, splash and other characteristics of explosion in concrete.

This work shows that FE-SPH coupling method has good performance to solve the explosion problem.

This paper sets out to present a fully explicit smoothed particle hydrodynamics (SPH) method to solve non‐Newtonian fluid flow problems.

The governing equations are momentum equations along with the continuity equation which are described in a Lagrangian framework. A new treatment similar to that used in Eulerian formulations is applied to viscous terms, which facilitates the implementation of various inelastic non‐Newtonian models. This approach utilizes the exact forms of the shear strain rate tensor and its second principal invariant to calculate the shear stress tensor. Three constitutive laws including power‐law, Bingham‐plastic and Herschel‐Bulkley models are studied in this work. The imposition of the incompressibility is fulfilled using a penalty‐like formulation which creates a trade‐off between the pressure and density variations. Solid walls are simulated by the boundary particles whose positions are fixed but contribute to the field variables in the same way as the fluid particles in flow field.

The performance of the proposed algorithm is assessed by solving three test cases including a non‐Newtonian dam‐break problem, flow in an annular viscometer using the aforementioned models and a mud fluid flow on a sloping bed under an overlying water. The results obtained by the proposed SPH algorithm are in close agreement with the available experimental and/or numerical data.

In this work, only inelastic non‐Newtonian models are studied. This paper deals with 2D problems, although extension of the proposed scheme to 3D is straightforward.

This study shows that various types of flow problems involving fluid‐solid and fluid‐fluid interfaces can be solved using the proposed SPH method.

Using the proposed numerical treatment of viscous terms, a unified and consistent approach was devised to study various non‐Newtonian flow models.

The purpose of this paper is to improve the 2D incompressible smoothed particle hydrodynamics (ISPH) method by working on the wall boundary conditions in ISPH method. Here, two different wall boundary conditions in ISPH method including dummy wall particles and analytical kernel renormalization wall boundary conditions have been discussed in details.

The ISPH algorithm based on the projection method with a divergence velocity condition with improved boundary conditions has been adapted.

The authors tested the current ISPH method with the improved boundary conditions by a lid-driven cavity for different Reynolds number 100 ≤ Re ≤ 1,000. The results are well validated with the benchmark problems.

In the case of dummy wall boundary particles, the homogeneous Newman boundary condition was applied in solving the linear systems of pressure Poisson equation. In the case of renormalization wall boundary conditions, the authors analytically computed the renormalization factor and its gradient based on a quintic kernel function.