Search results

This paper aims to introduce a three-dimensional smoothed particle hydrodynamics (SPH) framework for simulating supercooled large droplets (SLD) dynamics at aeronautical speeds.

To include the effects of the surrounding air, a multiphase model capable of handling high density-ratio problems is adopted. A diffusive term is incorporated to smooth the density field and avoid numerical instabilities. Additionally, a particle shifting technique is used to eliminate anisotropic particle distributions.

The framework is validated against low-speed droplet impingement experimental results and then applied to the droplet impingement at high speeds typical of SLD conditions. Preliminary parametric studies are conducted to investigate the post-impact splashing. It is observed that a thicker water film can decrease the crown diameter and a smaller impact angle can suppress upward and forward splashing.

A three-dimensional multiphase SPH framework for SLD dynamics at a wide range of impact speed is developed and validated. The effects of particle resolution, water film thickness and impact angle on the post-impact crown evolution are investigated.

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.

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.

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 purpose of this study is to investigate the effects of the shear-thinning viscoelastic behavior of the surrounding matrix on droplet deformation by weakly compressible smoothed particle hydrodynamics (WC-SPH). Also, the effect of the presence of another droplet is examined.

A modified consistent weakly compressible SPH method is proposed. After code verification, a complete parameter study is performed for a drop under the simple shear flow of a Giesekus liquid. The investigated parameters are 0.048≤Ca ≤ 14.4, 0.1≤c ≤ 10, 0.04≤De ≤ 10, 0≤α ≤ 1 and 0.12≤Re ≤ 12.

It is demonstrated that the rheological behavior of the surrounding fluid could dramatically affect the droplet deformation. It is shown that the droplet deformation is increased by increasing Re and Ca. In contrast, the droplet deformation is decreased by increasing a, De and polymer content. Also, it is indicated the presence of another droplet could drastically affect the flow field, and the primary stress difference (N1) is resonated between two droplets.

The main originality of this paper is to introduce a new consistent WC-SPH algorithm. The proposed method is very versatile for tackling the shear-thinning viscoelastic multiphase problems. Furthermore, a complete parameter study is performed for a drop under the simple shear flow of Giesekus liquid. Another novelty of the current paper is studying the effect of the presence of a second droplet. To the best of the authors’ knowledge, this is performed for the first time.

The purpose of this study is to apply the incompressible smoothed particle hydrodynamics (ISPH) method to simulate the natural convection flow from an inner heated Y-fin inside Y-shaped enclosure filled with nanofluid.

The dimensionless governing partial differential equations are described in the Lagrangian form and solved by an implicit scheme of the ISPH method. The embedded Y-fin is kept at a high temperature T_h with variable heights during the simulations. The lower area of Y-shaped enclosure is squared with width L = 1 m and its side-walls are kept at a low temperature T_c. The upper area of the Y-shaped enclosure is V-shaped with width 0.5 L for each side and its walls are adiabatic.

The performed simulations revealed that the height of the inner heated Y-fin plays an important role in the heat transfer and fluid flow inside the Y-shaped enclosure, where it enhances the heat transfer. Rayleigh number augments the buoyancy force inside the Y-shaped enclosure and, consequently, it has a strong impact on temperature distributions and strength of the fluid flow inside Y-shaped enclosure. Adding more concentration of the nanofluid until 10% has a slight effect on the temperature distributions and it reduces the strength of the fluid flow inside Y-shaped enclosure. In addition, the average Nusselt number is measured along the inner heated Y-fin and it grows as the Rayleigh number increases. The average Nusselt number is decreasing by adding more concentrations of the nanofluid.

An improved ISPH method is used to simulate the natural convection flow of Y-fin embedded in the Y-shaped enclosure filled with a nanofluid.

This paper aims to introduce a numerical investigation of aquatic locomotion using the smoothed particle hydrodynamics (SPH) method.

To model this problem, a simple improved SPH algorithm is presented that can handle complex geometries using updatable dummy particles. The computational code is validated by solving the flow over a two-dimensional cylinder and comparing its drag coefficient for two different Reynolds numbers with those in the literature.

Additionally, the drag coefficient and vortices created behind the aquatic swimmer are quantitatively and qualitatively compared with available credential data. Afterward, the flow over an aquatic swimmer is simulated for a wide range of Reynolds and Strouhal numbers, as well as for the amplitude envelope. Moreover, comprehensive discussions on drag coefficient and vorticity patterns behind the aquatic are made.

It is found that by increasing both Reynolds and Strouhal numbers separately, the anguilliform motion approaches the self-propulsion condition; however, the vortices show different pattern with these increments.

The purpose of this paper is to perform the computational fluid dynamics (CFD) simulation with experimental validation to investigate the particle segregation effect in abrupt and smooth shapes circulating fluidized bed (CFB) risers.

The experimental investigations were carried out in lab-scale CFB systems and the CFD simulations were performed by using commercial software BARRACUDA. Special attention was paid to investigate the gas-particle flow behavior at the top of the riser with three different superficial velocities, namely, 4, 6 and 7.7 m/s. Here, a CFD-based noble simulation approach called multi-phase particle in cell (MP-PIC) was used to investigate the effect of traditional drag models (Wen-Yu, Ergun, Wen-Yu-Ergun and Richardson-Davidson-Harrison) on particle flow characteristics in CFB riser.

Findings from the experimentations revealed that the increase in gas velocity leads to decrease the mixing index inside the riser. Moreover, the solid holdup found more in abrupt riser than smooth riser at the constant gas velocity. Despite the more experimental investigations, the findings with CFD simulations revealed that the MP-PIC approach, which was combined with different drag models could be more effective for the practical (industrial) design of CFB riser. Well agreement was found between the simulation and experimental outputs. The simulation work was compared with experimental data, which shows the good agreement (<4%).

The experimental and simulation study performed in this research study constitutes an easy-to-use with different drag coefficient. The proposed MP-PIC model is more effective for large particles fluidized bed, which can be helpful for further research on industrial gas-particle fluidized bed reactors. This study is expected to give throughout the analysis of CFB hydrodynamics with further exploration of overall fluidization.

This paper aims to conduct numerical simulations of unsteady natural/mixed convection in a cavity with fixed and moving rigid bodies and different boundary conditions using the incompressible smoothed particle hydrodynamics (ISPH) method.

In the ISPH method, the pressure evaluation is stabilized by including both of divergence of velocity and density invariance in solving pressure Poisson equation. The authors prevented the particles anisotropic distributions by using the shifting technique.

The proposed ISPH method exhibited good performance in natural/mixed convection in a cavity with fixed, moving and free-falling rigid body. In natural convection, the authors investigated the effects of an inner sloshing baffle as well as fixed and moving circular cylinders on the heat transfer and fluid flow. The heated baffle has higher effects on the heat transfer rate compared to a cooled baffle. In the mixed convection, a free-falling circular cylinder over a free surface cavity and heat transfer in the presence of a circular cylinder in a lid-driven cavity are simulated. Fixed or moving rigid body in a cavity results in considerable effects on the heat transfer rate and fluid flow.

The authors conducted numerical simulations of unsteady natural/mixed convection in a cavity with fixed and moving rigid bodies and different boundary conditions using the ISPH method.

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.