Incompressible Smoothed Particle Hydrodynamics Research Articles

Magnetic impacts on double diffusion of a non‐Newtonian NEPCM in a grooved cavity: ANN model with ISPH simulations

AbstractEmploying phase change materials (PCMs) offers the advantage of storing and releasing thermal energy while ensuring temperature stability. This characteristic makes PCMs valuable for reducing energy usage across various industrial applications. To explore the magnetic effects on double diffusion of a non‐Newtonian nano‐encapsulated phase change material (NEPCM) in a grooved cavity, the present study combined the incompressible smoothed particle hydrodynamics (ISPH) approach with an artificial neural network (ANN) model. The grooved shape is made up of three constructed grooves: triangular, curved, and rectangular grooves. In the cavity's walls, three segments of boundaries are considered as , , and . The ANN model correctly predicted the mean Nusselt number and Sherwood number when merged with current ISPH simulations. The study's novelty lies in exploring three distinct thermal and mass scenarios regarding double diffusion of a non‐Newtonian NEPCM within an innovative grooved domain. The relevant parameters include the fractional‐time derivative , power‐law index , Rayleigh number , Hartmann number , Soret–Dufour numbers (Sr and Du), and Lewis number Le. The obtained simulations present the significance of distinct boundary conditions in changing the velocity field, heat capacity ratio, temperature, and concentration in a grooved cavity. The fractional parameter accelerates the shift from unstable to steady condition. The increase in from 1.1 to 1.5 results in a 44.5% drop in the velocity maximum. Because of the Lorentz effect of a magnetic field, increasing from 0 to 50 reduces the maximum velocity by 20.9%.

Incompressible Smoothed Particle Hydrodynamics Simulation of Sediment Erosion around Submarine Pipelines

Sediment erosion around submarine pipelines is a popular topic, widely investigated in both ocean and submarine-pipeline engineering. In this paper, the incompressible smoothed-particle hydrodynamics (ISPH) method is modified for simulation of local scouring process around the submarine pipeline under the action of unidirectional flow. The erosion model is based on the Clear Water Particle–Turbid Water Particle–Critical Shear Stress (CWP-TWP-CSS) concept, and a sand–water two-phase model is proposed to deal with the sediment-entrained flow. The results of the numerical simulation are compared with the experimental data to verify the accuracy and applicability of the numerical model. The scouring process around the pipeline is investigated under different conditions, i.e., pipeline diameters, gap ratios, and flow velocities. The ISPH model is further used to study the flow characteristics of the scour pits around the submarine pipeline and the influence of the vortices on the maximum scour depth, to provide a theoretical basis for the stability design of submarine pipelines.

Extreme wave impacts on a suspended box above free surface

A systematic approach incorporating precision-controlled experiments and high-fidelity numerical simulations has been employed to delve deeper into the physics of extreme wave impact pressures. Using a suite of high-speed imaging and high-resolution pressure sensing techniques that are accurately synchronized with the wave generation, a consistent and well-correlated set of data providing details of the focusing incident wave front, flow velocities, and impact pressure time histories has been successfully obtained. The numerical simulations using a validated incompressible smoothed particle hydrodynamics (ISPH) provided further details that were not tractable in the experiments. Based on the results obtained, there were five broad impact scenarios identified, ranging from focusing crest impacts to plunging jet impacts with varying entrapped air. Among these impact types, the highest impact pressures were those associated with a focusing wave crest. It has been found that high impact pressures were correlated with a convergence of the horizontally impinging overturning crest and the vertical surge of the wall boundary water mass, but with the overturning crest impingement just ahead of the arrival of the wall boundary water mass, which in essence amounted to the impingement of a focusing wave front. Depending on the relative arrival times of the impinging crest and the wall boundary upsurge, the impact scenarios could vary from an up-slosh to a jet impingement with a clear entrapped air pocket, followed by a flip-through of the incident crest. For the scenario with the highest impact pressures, broadly classified as type II impact in this paper, the peak pressures could reach as high as 85ρC2, an order of magnitude higher than impact scenarios without the focusing effect, such as the scenario with a plunging jet impact with entrapped air (classified as types III and IV impact in this study). As the volume of entrapped air for this scenario was relatively small, the one-phase ISPH model used in this study was able to capture the peak pressure characteristics, including the peak pressure amplitudes. However, in scenarios with significantly highly entrapped air during impact, a two-phase model would be necessary.

A hybrid method combining ISPH with graph neural network for simulating free-surface flows

The incompressible Smoothed Particle Hydrodynamics (ISPH) is a popular Lagrangian Particle method. In the conventional ISPH method for simulating free-surface flows, the pressure-projection phase, which solves the pressure Poisson's equation (PPE), is the most time-consuming. In this paper, we propose a novel hybrid method by combining the graph neural network (GNN) with the ISPH for modelling the free-surface flows. In the new hybrid method, the graph neural network (GNN) is employed to replace solving the PPE for pressure in the conventional ISPH. To the best of knowledge of the authors, this is the first attempt to combine the GNN with ISPH model in a Lagrangian formulation. The performance of the hybrid method will be evaluated by comparing its results with experimental data, analytical solution or numerical results from other methods for three benchmark test cases: dam breaking, sloshing wave and solitary wave propagation. In addition, the potential of generalization of the hybrid method will be studied by applying it with the GNN model trained on data for relatively simple cases to simulate more complex cases. It will be demonstrated that the hybrid method does not only give satisfactory results, but also shows good potential of generalization. In addition, the new method will be demonstrated to require the computation time which can be 80 times less than the conventional ISPH for estimating pressure for cases with a large number of particles that is usually needed in the practical free-surface flows simulation using ISPH.

The ISPH simulation of Rayleigh-Taylor instability and multi-phase flow within porous media

The present study applies the incompressible smoothed particle hydrodynamics (ISPH) method to simulate multi-phase flow instabilities, specifically Rayleigh-Taylor instability (RTI) and mixing in a sloshing tank, within various porous media. This study is the first numerical attempt to simulate RTI in diverse porous structures. Porous media are modeled with the non-Darcy approach based on linear and nonlinear factors. The Smagorinsky model is employed to implement turbulence stress. Two validation tests were performed to verify the effectiveness of the ISPH technique, the first evaluating pressure over still water in a tank, the second RTI in free fluid, and excellent agreement was obtained in both tests. The results indicate that porosity parameters play a crucial role in controlling the RTI phenomenon and mixing in the sloshing tank. Low porosity leads to high porous resistance, which significantly reduces RTI development and the mixing processes between the two fluids by hindering the interaction of the heavier and lighter fluids. This study shows the effectiveness of the ISPH method in handling complex interfaces in two-phase mixing flows within porous media.

Comparative study of WCSPH, EISPH and explicit incompressible-compressible SPH (EICSPH) for multi-phase flow with high density difference

This study presents three Smoothed Particle Hydrodynamics (SPH) methods capable of handling high-density differences in violent incompressible multiphase flows. The conventional Weakly Compressible SPH (WCSPH) is reformulated into a quasi-Lagrangian framework based on Arbitrary Lagrangian–Eulerian (ALE) context. The Explicit Incompressible SPH (EISPH) method is extended to handle multiphase flows and reformulated in the ALE framework. The Explicit Incompressible-Compressible SPH (EICSPH) method, which can handle both compressible and incompressible phases, is developed by combining these two approaches in a fully explicit algorithm. The proposed methods are validated and compared through four benchmark problems including Rayleigh-Taylor instability, hydrostatic, liquid sloshing and dam break problems. Firstly, the stability and accuracy of interface modeling of the three methods were verified through Rayleigh-Taylor instability with small density differences. In the hydrostatic problem with a large density difference, the results from all three methods exhibit good agreement in the pressure distribution. Particularly, EICSPH demonstrates the faster convergence when compared to WCSPH, with EISPH showing fastest convergence overall. In the transient sloshing problem, all three methods quantitatively converged to the experimental results and exhibited good agreement, although slight pressure noise was observed in WCSPH initially. Finally, in the dam break problem, all three methods successfully simulated sharp interfaces without non-physical voids. The temporal variations of pressure and height were well predicted by all three methods, with EISPH showing better conformity to the water-air interface morphology obtained from other incompressible numerical methods compared to WCSPH and EICSPH. Additionally, the impact of air speed of sound was examined in EICSPH, where while the difference in pressure variation at the probe was minimal, differences were observed in energy dissipation and the progression of the free surface due to the air cushion effect.

Adjoined ISPH method and artificial intelligence for thermal radiation on double diffusion inside a porous L-shaped cavity with fins

PurposeThe purpose of this study is to adapt the incompressible smoothed particle hydrodynamics (ISPH) method with artificial intelligence to manage the physical problem of double diffusion inside a porous L-shaped cavity including two fins.Design/methodology/approachThe ISPH method solves the nondimensional governing equations of a physical model. The ISPH simulations are attained at different Frank–Kamenetskii number, Darcy number, coupled Soret/Dufour numbers, coupled Cattaneo–Christov heat/mass fluxes, thermal radiation parameter and nanoparticle parameter. An artificial neural network (ANN) is developed using a total of 243 data sets. The data set is optimized as 171 of the data sets were used for training the model, 36 for validation and 36 for the testing phase. The network model was trained using the Levenberg–Marquardt training algorithm.FindingsThe resulting simulations show how thermal radiation declines the temperature distribution and changes the contour of a heat capacity ratio. The temperature distribution is improved, and the velocity field is decreased by 36.77% when the coupled heat Cattaneo–Christov heat/mass fluxes are increased from 0 to 0.8. The temperature distribution is supported, and the concentration distribution is declined by an increase in Soret–Dufour numbers. A rise in Soret–Dufour numbers corresponds to a decreasing velocity field. The Frank–Kamenetskii number is useful for enhancing the velocity field and temperature distribution. A reduction in Darcy number causes a high porous struggle, which reduces nanofluid velocity and improves temperature and concentration distribution. An increase in nanoparticle concentration causes a high fluid suspension viscosity, which reduces the suspension’s velocity. With the help of the ANN, the obtained model accurately predicts the values of the Nusselt and Sherwood numbers.Originality/valueA novel integration between the ISPH method and the ANN is adapted to handle the heat and mass transfer within a new L-shaped geometry with fins in the presence of several physical effects.

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.

Incompressible smoothed particle hydrodynamics solutions of double diffusion of nanoparticle-enhanced phase change materials in vertical mounted circular cylinders

The scheme of incompressible smoothed particle hydrodynamics (ISPH) method is developed to solve the time-conformable systems that control the double diffusion of nanoparticle-enhanced phase change materials (NEPCM) inside two circular cylinders saturated by a porous medium. The use of the time-conformable operator in the ISPH approach for addressing the double diffusion of NEPCM in a new domain is what distinguishes this work. This investigation revealed a comprehensive rate of change that impacts the quickening and slowdown of the monotonicity for the nonlinear systems. The physical parameters are time-conformable operator α , time parameter τ , nanoparticles parameter ϕ , Hartmann number Ha , thermal radiation parameter Rd , Darcy parameter Da , Rayleigh number Ra , and fusion temperature θ f . The principal outcomes of the numerical simulations revealed that the nanofluid flow is slowing down by an augmentation in ϕ , Ha , and Rd . The Rayleigh number is working well in enhancing the nanofluid flow and thermosolutal convection within the two circular cylinders. As a result, variable values of Ra are affecting the phase change material. The fusion temperature parameter controls the position/intensity of the phase change material within the two circular cylinders. Besides, the numerical results indicate the significance of α on the transition process of reaching the steady state. The thermosolutal convection and nanofluid speed are varied with the changes on the time-conformable derivative.

Integrating artificial intelligence with numerical simulations of Cattaneo-Christov heat flux on thermosolutal convection of nano-enhanced phase change materials in Bézier-annulus

The numerical analysis based on incompressible smoothed particle hydrodynamics (ISPH) is introduced to examine the impacts of Cattaneo-Christov (Ca-Ch) heat flux and exothermic chemical reaction on thermosolutal convection of nano-enhanced phase change materials (NEPCM) in Bézier-annulus. The used annulus is formed between inner connected Bézier curves and outer connected spline-Bézier curves. The inner shape of connected Bézier curves is maintained at Th&Ch, left/right walls of spline-Bézier curves are kept at Tc&Cc and other walls are adiabatic. The governing equations, after being converted into non-dimensional form, have been solved by ISPH which is an accurate meshless algorithm the treatment of internal flows inside complex geometries. The simulations are executed for Frank-Kamenetskii number FK, Ca-Ch heat flux δC, buoyancy ratio parameter N, Soret-Dufour numbers SrDu, Rayleigh number Ra, and nanoparticle parameter φ on thermosolutal convection of a suspension fluid. From the numerical simulation values for the average Nusselt and Sherwood numbers (Nu¯, Sh¯) were obtained for some fluid flow scenarios. Then, based on those values an artificial neural network (ANN) model was developed to predict the values of Nu¯, and Sh¯ without the need to perform the ordinary simulations again for the new cases which is a high cost compared to ANN. From the obtained simulations, it was concluded that the ANN model is an accurate tool to be used to predict the needed values. Also, the Frank-Kamenetskii number significantly influences the enhancement process of the temperature distributions and velocity field as well as phase change material in Bézier-annulus.

Integrating ISPH simulations with machine learning for thermal radiation and exothermic chemical reaction on heat and mass transfer in spline/triangle star annulus

This work integrates incompressible smoothed particle hydrodynamics (ISPH) with machine learning (ML) to examine the uses of two different domain shapes, spline star-domain, and triangle star-domain, on double diffusion of nano-encapsulated phase change material (NEPCM). The effects of thermal radiation and exothermic chemical reactions are conducted in this work. The hexagonal-shaped is located inside the center of the star-shape carrying (Th&Ch). The left/right walls of the star shape are kept at (Tc&Cc) and the other walls have an adiabatic condition. It is recommended to examine the smoothness of closed curves in the heat and mass transfer by changing the triangle to spline curves. The present investigations of transmission of heat and mass of NEPCM within a complicated star shape can be employed in cooling systems of electronic/battery, heat exchangers, food processing, and chemical engineering. The ISPH simulations are executed for the pertinent parameters, Darcy number, Frank-Kamenetskii number, Rayleigh number, Dufour number, Soret number, and thermal radiation parameter. The performed simulations showed the significance of the spline star domain compared to the triangle star domain in smoothening the nanofluid flow and double-diffusive convection in an annulus. Frank-Kamenetskii number is a promising factor in the enhancement of heat transfer and controlling the phase change zone within an annulus. At a spline-star cavity, the maximum velocity rises by 32.32 % and 54.95 % when the Frank-Kamenetskii number grows from 0 to 5, and thermal radiation increases from 0 to 5. Soret number works effectively in enhancing the distribution of concentration at large temperature differences. Machine learning plays a vital role in many engineering applications. Hence, ML's powerful capabilities to predict the average Nu‾ and Sh‾. In the proposed ML model, the dimensionless time and thermal radiation parameter are used as input to the ML model, while their corresponding numeral estimated values of Nu‾ and Sh‾ are used as output targets for the ML model. The values of Nu‾ and Sh‾ are predicted by the suggested ML model with a comparatively small error less than 1 %.

A ROTATING WAVY CYLINDER ON BIOCONVECTION FLOW OF NANOENCAPSULATED PHASE CHANGE MATERIALS IN A FINNED CIRCULAR CYLINDER: ISPH SIMULATIONS

The originality of this study is the introduction of numerical investigations on the bioconvection flow of nano-encapsulated phase change materials (NEPCMs) with oxytactic microorganisms in a new configuration of a circular annulus with a rotating wavy inner cylinder. The incompressible smoothed particle hydrodynamics (ISPH) method was applied to solve the governing partial differential equations for the velocity, temperature, concentration, and density of motile microorganisms. Compared with the conventional mesh-based method, this mesh-free, particle-based approach offers strong advantages in the simulation of complex problems with free surfaces and moving boundaries with large displacements. The pertinent parameters are the undulation number (<i>N<sub>und</sub></i> = 2-36), bioconvection Rayleigh number (<i>Ra<sub>b</sub></i> = 1-1000), Darcy parameter (Da = 10<sup>-5</sup>-10<sup>-2</sup>), length of the inner fin (<i>L<sub>Fin</sub></i> = 0.05-0.15), radius of the inner wavy cylinder (<i>R<sub>c</sub></i> = 0.05-0.25), Rayleigh number (Ra = 10<sup>3</sup>-10<sup>5</sup>), undulation amplitude of the inner wavy cylinder surface (<i>A</i> = 0.1-0.4), and frequency parameter (<i>ω </i>= 1-5). The undulation number of the inner wavy cylinder enhanced the flow of the oxytactic microorganisms and isotherms, whereas it had the reverse effect on the velocity, decreasing the maximum velocity by 26.56%. In addition, the comparatively high undulation amplitude and frequency increased the average Nusselt and Sherwood numbers. It was found that the embedded wavy cylinder interacting with fins plays an important role in enhancing heat transfer and the bioconvection flow within a closed domain.

Heat and mass transport of nano-encapsulated phase change materials in a complex cavity: An artificial neural network coupled with incompressible smoothed particle hydrodynamics simulations

<abstract> <p>This work simulates thermo-diffusion and diffusion-thermo on heat, mass transfer, and fluid flow of nano-encapsulated phase change materials (NEPCM) within a complex cavity. It is a novel study in handling the heat/mass transfer inside a highly complicated shape saturated by a partial layer porous medium. In addition, an artificial neural network (ANN) model is used in conjunction with the incompressible smoothed particle hydrodynamics (ISPH) simulation to forecast the mean Nusselt and Sherwood numbers ($ \stackrel{-}{Nu} $ and $ \stackrel{-}{Sh} $). Heat and mass transfer, as well as thermo-diffusion effects, are useful in a variety of applications, including chemical engineering, material processing, and multifunctional heat exchangers. The ISPH method is used to solve the system of governing equations for the heat and mass transfer inside a complex cavity. The scales of pertinent parameters are fusion temperature $ {\theta }_{f} = 0.05-0.95 $, Rayleigh number $ Ra = {10}^{3}-{10}^{6} $, buoyancy ratio parameter $ N = -2-1 $, Darcy number $ Da = {10}^{-2}-{10}^{-5} $, Lewis number $ Le = 1-20 $, Dufour number $ Du = 0-0.25 $, and Soret number $ Sr = 0-0.8 $. Alterations of Rayleigh number are effective in enhancing the intensity of heat and mass transfer and velocity field of NEPCM within a complex cavity. The high complexity of a closed domain reduced the influences of Soret-Dufour numbers on heat and mass transfer especially at the steady state. The fusion temperature works well in adjusting the intensity and location of a heat capacity ratio inside a complex cavity. The presence of a porous layer in a cavity's center decreases the velocity field within a complex cavity at a reduction in Darcy number. The goal values of $ \stackrel{-}{Nu} $ and $ \stackrel{-}{Sh} $ for each data point are compared to those estimated by the ANN model. It is discovered that the ANN model's $ \stackrel{-}{Nu} $ and $ \stackrel{-}{Sh} $ values correspond completely with the target values. The exact harmony of the ANN model prediction values with the target values demonstrates that the developed ANN model can forecast the $ \stackrel{-}{Nu} $ and $ \stackrel{-}{Sh} $ values precisely.</p> </abstract>

Incompressible [formula omitted]-SPH via artificial compressibility

Smoothed particle hydrodynamics using artificial compressibility (ACSPH) is developed, with the inclusion of pressure smoothing terms. Theoretical links between pressure/velocity correction incompressible SPH and artificial compressibility are explored, illustrating that ACSPH may be considered an extension of, or closely related to, the δ-SPH method. An implicit dual-time integration procedure is used to enforce an incompressible solution at every time-step, removing acoustic effects arising from the common assumption of weak compressibility. An established weakly-compressible quasi-Lagrangian δ-SPH method is used for comparison against ACSPH, and a series of test cases show that ACSPH provides a similar solution cost to δ-SPH. However, the residual acoustic effects in δ-SPH are removed entirely in ACSPH, providing improved pressure prediction capabilities across all test cases, including intense fluid impacts. Improved modelling of fluid–structure-interaction cases and coupled energy dissipation are also recorded as a result of correctly capturing incompressible flow.

A Current Loop Model for the Fast Simulation of Ferrofluids.

Ferrofluids are oil-based liquids containing magnetic particles that interact with magnetic fields without solidifying. Leveraging the exploration of new applications of these promising materials (such as in optics, medicine and engineering) requires high fidelity modeling and simulation capabilities in order to accurately explore ferrofluids in silico. While recent work addressed the macroscopic simulation of large-scale ferrofluids using smoothed-particle hydrodynamics (SPH), such simulations are computationally expensive. In their work, the Kelvin force model has been used to calculate interactions between different SPH particles. The application of this model results in a force pointing outwards with respect to the fluid surface causing significant levitation problems. This drawback limits the application of more advanced and efficient SPH frameworks such as divergence-free SPH (DFSPH) or implicit incompressible SPH (IISPH). In this contribution, we propose a current loop magnetic force model which enables the fast macroscopic simulation of ferrofluids. Our new force model results in a force term pointing inwards allowing for more stable and fast simulations of ferrofluids using DFSPH and IISPH.

Circular rotation of different structures on natural convection of nanofluid-mobilized circular cylinder cavity saturated with a heterogeneous porous medium

The incompressible smoothed particle hydrodynamics (ISPH) method is utilized for studying the circular rotations of three different structures, circular cylinder, rectangle and triangle centered in a circular cylinder cavity occupied by Al2O3 nanofluid. The novelty of this work is appearing in simulating the circular rotations of different solid structures on natural convection of a nanofluid-occupied a circular cylinder. The circular cylinder cavity is suspended by heterogeneous/homogeneous porous media. The embedded structures are taken as a circular cylinder, rectangle and triangle with equal areas. The first thermal condition considers the whole structure is heated, the second thermal condition considers the half of the structure is heated and the other is cooled and the third thermal condition considers the quarter of the structure is heated and the others are cooled. The outer boundary of cylinder cavity is cooled. Due to the small angular velocity ω=3.15 (low rotational speeds), then the natural convection case will be considered only. The results are representing the temperature, velocity fields. The simulations revealed that the presence of the inner hot/cold structures affects on the velocity distributions and temperature field inside a circular cylinder cavity. The triangle shape has introduced the highest temperature distributions and maximum values of the velocity fields compare to other shapes inside a circular cylinder cavity. The homogeneous porous level reduces the maximum values of velocity field by 25% compared to the heterogeneous porous level.

A fully explicit incompressible smoothed particle hydrodynamics approach for modeling transient heat transfer and thermo-capillary flows

This work presents a Lagrangian meshless method for modeling transient heat transfer and thermo-capillary effects at the interface of two-phase incompressible fluids. To simulate the thermo-capillary effect, we add the Marangoni force to the continuum surface force (CSF) model, by considering the forces caused by the surface tension gradient which acts tangentially to the interface position. In the current study, we extend our previously proposed model which uses a fully explicit incompressible smoothed particle hydrodynamics (EISPH) approach along with corrected SPH and VKF kernel function, to obtain stable and accurate results in non-isothermal single and multi-phase problems.The model is validated for several test cases, including transient heat conduction, natural convection heat transfer, and thermo-capillary droplet migration, against conventional numerical methods. The validated model is used to study the effects of several dimensionless parameters on thermo-capillary droplet migration. The results show that the proposed EISPH method can model accurately complex heat transfer problems such as thermo-capillary induced motion.

Simulation of granular surface flows using incompressible non-Newtonian SPH (INNSPH) method

This study presents a comprehensive continuum approach to modeling granular surface flows using the Incompressible Non-Newtonian Smoothed Particle Hydrodynamics (INNSPH) method. The model integrates an incompressible SPH method, where pressure is computed by solving a pressure-Poisson equation with a non-Newtonian rheology. Projection-based SPH has a strong and rigorous theoretical basis can lead to stable and accurate simulations. This paper aims to develop a projection-based SPH method for modeling complex granular flows. INNSPH method combines incompressible SPH and μ(I) rheology models to simulate non-Newtonian viscoplastic fluid behavior in granular flows. Surface deposition, and velocity and pressure fields are estimated using a rheology model. Validation is performed by comparing simulation results, including surface profiles and run-out distances with experimental data, showing good agreement. INNSPH method is then employed to model a 2D column collapse and a dam-break scenario on horizontal and inclined planes for sands and glass beads. The method captures three distinct final deposition morphologies; truncated cone, conic, and Mexican hat shapes in collapsing columns. INNSPH computational efficiency enables simulations of large-scale granular flows, and its open-source nature facilitates collaboration and further modifications by researchers. The method provides a good understanding of granular flow dynamics in various industrial applications including powder technology.

Boundary Treatment with Additive Boundary Particles for Incompressible Smoothed Particle Hydrodynamics Method

A new boundary treatment considering uniform particle distribution using additive boundary particles for the incompressible smoothed particle hydrodynamics (ISPH) method is proposed. With this new boundary treatment, the solid boundary can be easily modeled just following the geometry configuration and consequently the simulations of fluid–structure interaction problems can be simplified. The efficiency of this new boundary treatment is analyzed, and the implementations are compared with other exiting boundary treatments such as repulsive force and ghost particles. Simulations of the 2D dam-breaking case, static tank, and 2D u-tube test are carried out as examples to demonstrate the performance of this new boundary treatment. It is shown that better predictions can be gained with the new boundary treatment for the pressure distribution in these cases than those obtained from other boundary treatments.

An incompressible SPH numerical model for simulating wave and non-Newtonian mud interaction

An incompressible SPH numerical model for simulating wave and non-Newtonian mud interaction