Sharif Digital Repository

In this paper a smoothed particle hydrodynamics (SPH) method is introduced for simulating two-dimensional incompressible non-Newtonian fluid flows, and the non-Newtonian effects in the flow of a fluid which can be modelled by generalized Newtonian constitutive equations are investigated. Two viscoplastic models including Bingham-plastic and power-law models are considered along with the Newtonian model. The governing equations include the conservation of mass and momentum equations in a pseudo-compressible form. The spatial discretization of these equations is achieved by using the SPH method. The temporal discretization algorithm is a fully explicit two-step predictorcorrector scheme. In... 

This paper deals with numerical modeling of water flow which is generated by the break of a dam. The problem is solved by applying a new Incompressible Smoothed Particle Hydrodynamics (ISPH) algorithm based on projection method. The proposed ISPH model has two steps. In the first step, the incompressibility of fluid is maintained regarding to the changes of intermediate and initial particles densities at the first half-time step (stability step). In the second step, by computing the divergence of the intermediate secondary velocity at the second half-time step (accuracy step), the incompressibility is satisfied completely. In fact, by using this method both stability and accuracy are... 

Several schemes for discretization of first and second derivatives are available in Smoothed Particle Hydrodynamics (SPH). Here, four schemes for approximation of the first derivative and three schemes for the second derivative are examined using a theoretical analysis based on Taylor series expansion both for regular and irregular particle distributions. Estimation of terms in the truncation errors shows that only the renormalized (the first-order consistent) scheme has acceptable convergence properties to approximate the first derivative. None of the second derivative schemes has the first-order consistency. Therefore, they converge only when the particle spacing decreases much faster than... 

This paper deals with numerical modeling of two-phase liquid jet breakup using the smoothed particle hydrodynamics (SPH) method. Simulation of multiphase flows involving fluids with a high-density ratio causes large pressure gradients at the interface and subsequently divergence of numerical solutions. A modified procedure extended by Monaghan and Rafiee is employed to stabilize the sharp interface between the fluids. Various test cases such as Rayleigh-Taylor instability, two-phase still water and air bubble rising in water have been conducted, by which the capability of accurately capturing the physics of multiphase flows is verified. The results of these simulations are in a good... 

The rupturing process of a membrane, located between two fluids at the center of a three-dimensional channel, is numerically investigated. The smoothed particle hydrodynamics (SPH) and the finite element method (FEM) are used, respectively, for modeling the fluid and solid phases. A range of pressure differences and membrane thicknesses are studied and two different rupturing processes are identified. These processes differ in the time scale of the rupture, the location of the rupture initiation, the level of destruction and the driving mechanism. © 2019 by the authors  

Multiphase flow is a challenging area of computational fluid dynamics (CFD) due to their potential large topological change and close coupling between the interface and fluid flow solvers. As such, Lagrangian meshless methods are very well suited for solving such problems. In this paper, we present a new fully explicit incompressible Smoothed Particle Hydrodynamics approach (EISPH) for solving multiphase flow problems. Assuming that the change in pressure between consecutive time-steps is small, due to small time steps in explicit solvers, an approximation of the pressure for following time-steps is derived. To verify the proposed method, several test cases including both single-phase and... 

The Smoothed Particle Hydrodynamics (SPH) method is extended to solve magnetostatic problems involving magnetically interacting solid bodies. In order to deal with the jump in the magnetic permeability at a fluid-solid interface, a consistent SPH scheme is utilized and a modified formulation is proposed to calculate the magnetic force density along the interface. The results of the magnetostatic solver are verified against those of the finite element method. The governing fluid flow equations are discretized using the same SPH scheme, developing an efficient method for simulating the motion of paramagnetic solid bodies in a fluid flow. The proposed algorithm is applied to a benchmark problem... 

Instability of the interface between two immiscible fluids representing the so-called Kelvin-Helmholtz instability problem is studied using smoothed particle hydrodynamics method. Interfacial tension is included, and the fluids are assumed to be inviscid. The time evolution of interfaces is obtained for two low Richardson numbers Ri=0.01 and Ri=0.1 while Bond number varies between zero and infinity. This study focuses on the effect of Bond and Richardson numbers on secondary instability of a two-dimensional shear layer. A brief theoretical discussion is given concerning the linear early time regime followed by numerical investigation of the growth of secondary waves on the main billow....