Search results

This paper aims to examine the accuracy of several higher-order incompressible smoothed particle hydrodynamics (ISPH) Laplacian models and compared with the classic model (Shao and Lo, 2003).

The numerical errors in solving two-dimensional elliptic partial differential equations using the Laplacian models are investigated for regular and highly irregular node distributions over a unit square computational domain.

The numerical results show that one of the Laplacian models, which is newly developed by one of the authors (Shobeyri, 2019) can get the smallest errors for various used node distributions.

The newly proposed model is formulated by the hybrid of the standard ISPH Laplacian model combined with Taylor expansion and moving least squares method. The superiority of the proposed model is significant when multi-resolution irregular node distributions commonly seen in adaptive refinement strategies used to save computational cost are applied.

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.

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 develop an efficient and accurate mesh-less method to simulate free flows with continuous deformation in boundary positions.

A two-step pressure projection method in a Lagrangian form is used to solve the governing equations of mass and momentum conservation. In the first step, velocity field is calculated in which incompressibility is not enforced. In the second step, a pressure Poisson equation is applied to satisfy incompressibility conditions. The numerical proposed method is used for spatial discretization of the governing equations. Three benchmark-free surface problems, namely, dam break, solitary wave propagation and evolution of an elliptical bubble with available experimental results and analytical solutions, are used to test the accuracy of the proposed method. The results prove the accuracy of the method in simulating free surface problems.

The Voronoi diagram instead of kernel function summation can be used to estimate the particle or nodal volume concept in particle-based (mesh-less) methods for function approximation. This idea probably works well especially for highly irregular node distributions.

The continuous moving least squares shape functions are applied for function approximation, and the Voronoi diagram concept is also used to estimate region influence of computational nodal points or particle volumes. Combinations of these two concepts and finite differences formulation for first derivatives gives an accurate numerical model for Laplacian operator in the proposed method.

The purpose of this paper is to show how a surface tension model and an eddy viscosity based on the Smagorinsky sub‐grid scale model, which belongs to the Large‐Eddy Simulation (LES) theory for turbulent flow, have been introduced into ISPH (Incompressible smoothed particle hydrodynamics) method. In addition, a small modification in the source term of pressure Poisson equation has been introduced as a stabilizer for robust simulations. This stabilization generates a smoothed pressure distribution and keeps the total volume of fluid, and it is analogous to the recent modification in MPS.

The surface tension force in free surface flow is evaluated without a direct modeling of surrounding air for decreasing computational costs. The proposed model was validated by calculating the surface tension force in the free surface interface for a cubic‐droplet under null‐gravity and the milk crown problem with different resolution models. Finally, effects of the eddy viscosity have been discussed with a fluid‐fluid interaction simulation.

From the numerical tests, the surface tension model can handle free surface tension problems including high curvature without special treatments. The eddy viscosity has clear effects in adjusting the splashes and reduces the deformation of free surface in the interaction. Finally, the proposed stabilization appeared in the source term of pressure Poisson equation has an important role in the simulation to keep the total volume of fluid.

An incompressible smoothed particle hydrodynamics is developed to simulate milk crown problem using a surface tension model and the eddy viscosity.

The jet impingement usually accompanying large interface movement is studied by the in-house solver MLParticle-SJTU based on the modified moving particle semi-implicit (MPS) method, which can provide more accurate pressure fields and deformed interface shape. The comparisons of the pressure distribution and the shape of free surface between the presented numerical results and the analytical solution are investigated. The paper aims to discuss these issues.

To avoid the instability in traditional MPS, a modified MPS method is employed, which include mixed source term for Poisson pressure equation (PPE), kernel function without singularity, momentum conservative gradient model and highly precise free surface detection approach. Detailed analysis on improved schemes in the modified MPS is carried out. In particular, three kinds of source term in PPE are considered, including: particle number density (PND) method, mixed source term method and divergence-free method. Two typical kernel functions containing original kernel function with singularity and modified kernel function without singularity are analyzed. Three kinds of pressure gradient are considered: original pressure gradient (OPG), conservative pressure gradient (CPG) and modified pressure gradient (MPG). In addition, particle convergence is performed by running the simulation with various spatial resolutions. Finally, the comparison of the pressure fields by the modified MPS and by SPH is presented.

The modified MPS method can provide a reliable pressure distribution and the shape of the free surface compared to the analytical solution in a steady state after the water jet impinging on the wall. Specifically, mixed source term in PPE can give a reasonable profile of the shape of free surface and pressure distribution, while PND method adopted in the traditional MPS is not stable in simulation, and divergence-free method cannot produce rational pressure field near the wall. Two kernel functions show similar pressure field, however, the kernel function without singularity is preferred in this case to predict the profile of free surface and pressure on the wall. The shape of free surface by CPG and MPG is agreement with the analytical solution, while a great discrepancy can be observed by OPG. The pressure peak by MPG is closer to the analytical solution than that by CPG, while the pressure distribution on the right hand side of the pressure peak by latter is better match with the analytical solution than that by former. Besides, fine spatial resolution is necessary to achieve a good agreement with analytical results. In addition, the pressure field by the modified MPS is also quite similar to that by SPH, and this can further validate the reliable of current modified MPS.

The present modified MPS appears to be a stable and reliable tool to deal with the impinging jet flow problems involving large interface movement. Mixed source term in PPE is superior to PND adopted in the traditional MPS and divergence-free method. The kernel function without singularity is preferred to improve the computational accuracy in this case. CPG is a good choice to obtain the shape of free surface and the pressure distribution by jet impingement.

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.

The purpose of this paper is to propose an adaptive refinement strategy based on a posteriori error estimate for the efficient simulation of free surface flows using discrete least squares meshless (DLSM) method.

A pressure projection method is employed to discretize the governing equations of mass and momentum conservation in a Lagrangian form. The semi‐discretized equations are then discretized in space using the DLSM method, in which the sum of squared residual of the governing equations and their boundary conditions are minimized with respect to the unknown nodal parameters.

Since the position of the free surface is of great significant in free surface problems, a posteriori error estimator which automatically associates higher error to the nodes near the free surface is proposed and used along with a node moving refinement strategy to simulate the free surface problems more efficiently. To test the ability and efficiency of the proposed adaptive simulation method, two test problems, namely dam break and evolution of a water bubble, are solved and the results are presented and compared to those of analytical and experimental results.

Error estimate and adaptive refinement have been mostly used in confined and steady‐state flow. Here in this paper, a new attempt has been made to use these concepts in moving boundary problem.

The purpose of this study is to simulate the thermo-solutal convection resulting from a circular cylinder hanging in a rod inside a ∧-shaped cavity.

The two dimensional ∧-shaped cavity is filled by Al₂O₃-water nanofluid and saturated by three different levels of heterogeneous porous media. An incompressible smoothed particle hydrodynamics (ISPH) method is adopted to solve the governing equations of the present problem. The present simulations have been performed for the alteration of buoyancy ratio (−2≤N≤2), radius of a circular cylinder (0.05≤Rc≤0.3), a height of a rod (0.1≤Lh≤0.4), Darcy parameter (10−3≤Da≤10−5), Lewis number (1≤Le≤40), solid volume fraction (0≤ϕ≤0.06), porous levels (0≤η1=η2≤1.5)and various boundary-wall conditions.

The performed numerical simulations indicated the importance of embedded shapes on the distributions of temperature, concentration and velocity fields inside ∧-shaped cavity. Increasing buoyancy ratio parameter enhances thermo-solutal convection and nanofluid velocity. Adiabatic conditions of the vertical-walls of ∧-shaped cavity augment the distributions of the temperature and concentration. Regardless the Darcy parameter, a homogeneous porous medium gives the lowest values of a nanofluid velocity.

ISPH method is used to simulate thermo-solutal convection of a nanofluid inside a novel ∧-shaped cavity containing a novel embedded shape and heterogeneous porous media.

The purpose of this paper is to highlight the possibilities of a novel Lagrangian formulation in dealing with the solution of the incompressible Navier‐Stokes equations with very large time steps.

The design of the paper is based on introducing the origin of this novel numerical method, originally inspired on the Particle Finite Element Method (PFEM), summarizing the previously published theory in its moving mesh version. Afterwards its extension to fixed mesh version is introduced, showing some details about the implementation.

The authors have found that even though this method was originally designed to deal with heterogeneous or free‐surface flows, it can be competitive with Eulerian alternatives, even in their range of optimal application in terms of accuracy, with an interesting robustness allowing to use large time steps in a stable way.

With this objective in mind, the authors have chosen a number of benchmark examples and have proved that the proposed algorithm provides results which compare favourably, both in terms of solution time and accuracy achieved, with alternative approaches, implemented in in‐house and commercial codes.