Search results

A new fractional step method in conjunction with the finite element method is proposed for the analysis of the thermal convection and conduction in a fluid region expressed by the momentum equations, the equation of continuity and the energy equation. This paper focuses on the features of the present finite element method which gives a simple way of treating the Neumann boundary condition for the pressure Poisson equation. The applicability and effectiveness of the proposed scheme are illustrated through the numerical examples of the two‐dimensional natural convection flow in enclosures with several Rayleigh numbers.

As the penalized vortex-in-cell (pVIC) method is based on the vorticity-velocity form of the Navier–Stokes equation, the pressure variable is not incorporated in its solution procedure. This is one of the advantages of vorticity-based methods such as pVIC. However, dynamic pressure is an essential flow property in engineering problems. In pVIC, the pressure field can be explicitly evaluated by a pressure Poisson equation (PPE) from the velocity and vorticity solutions. How to specify far-field boundary conditions is then an important numerical issue. Therefore, this paper aims to robustly and accurately determine the boundary conditions for solving the PPE.

This paper introduces a novel non-iterative method for specifying Dirichlet far-field boundary conditions to solve the PPE in a bounded domain. The pressure field is computed using the velocity and vorticity fields obtained from pVIC, and the solid boundary conditions for pressure are also imposed by a penalization term within the framework of pVIC. The basic idea of our approach is that the pressure at any position can be evaluated from its gradient field in a closed contour because the contour integration for conservative vector fields is path-independent. The proposed approach is validated and assessed by a comparative study.

This non-iterative method is successfully implemented to the pressure calculation of the benchmark problems in both 2D and 3D. The method is much faster than all the other methods tested without compromising accuracy and enables one to obtain reasonable pressure field even for small computation domains that are used regardless of a source distribution (the right-hand side in the Poisson equation).

The strategy introduced in this paper provides an effective means of specifying Dirichlet boundary conditions at the exterior domain boundaries for the pressure Poisson problems. It is very efficient and robust compared with the conventional methods. The proposed idea can also be adopted in other fields dealing with infinite-domain Poisson problems.

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.

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 find the best solver for parallelizing particle methods based on solving Pressure Poisson Equation (PPE) by taking Moving Particle Semi-Implicit (MPS) method as an example because the solution for PPE is usually the most time-consuming part difficult to parallelize.

To find the best solver, the authors compare six Krylov solvers, namely, Conjugate Gradient method (CG), Scaled Conjugate Gradient method (SCG), Bi-Conjugate Gradient Stabilized (BiCGStab) method, Conjugate Gradient Squared (CGS) method with Symmetric Lanczos Algorithm (SLA) method and Incomplete Cholesky Conjugate Gradient method (ICCG) in terms of convergence, time consumption, parallel efficiency and memory consumption for the semi-implicit particle method. The MPS method is parallelized by the hybrid Open Multi-Processing (OpenMP)/Message Passing Interface (MPI) model. The dam-break flow and channel flow simulations are used to evaluate the performance of different solvers.

It is found that CG converges stably, runs fastest in the serial way, uses the least memory and has highest OpenMP parallel efficiency, but its MPI parallel efficiency is lower than SLA because SLA requires less synchronization than CG.

With all these criteria considered and weighed, the recommended parallel solver for the MPS method is CG.

The purpose of this study is to simulate the natural convection of a heated square shape embedded in a circular enclosure filled with nanofluid using an incompressible smoothed particle hydrodynamics (ISPH) method.

In the ISPH method, the evaluated pressure was stabilized by using a modified source term in solving the pressure Poisson equation. The divergence of the velocity was corrected, and the dummy particles were used to treat the rigid boundary. Dummy wall particles were initially settled in outer layers of the circular enclosure for preventing particle penetration and reducing the error of truncated kernel. The circular enclosure was partially filled with a porous medium near to the outer region. The single-phase model was used for the nanofluid, and the Brinkman–Forchheimer-extended Darcy model was used for the porous medium. Dummy wall particles were initially settled in outer layers of circular enclosure for preventing particle penetration and reducing error from the truncated kernel on the boundary.

The length of the inner square shape plays an important role in enhancing the heat transfer and reducing the fluid flow inside a circular enclosure. The porous layer represents a resistance force for the fluid flow and heat transfer, and, consequently, the velocity field and temperature distributions are reduced at the outer region of the circular cylinder. Then, the radius of the inner square shape, Darcy parameter and radius of the porous layer were considered the main factors for controlling the fluid flow and heat transfer inside a circular enclosure. The average Nusselt number decreases as the inner square length, radius of the porous layer and solid volume fraction increase.

The stabilized ISPH method is corrected for simulating the natural convection from an inner hot square inside a nanofluid-filled circular enclosure saturated with a partial layer of a porous medium.

To study the steady viscous incompressible electrically conducting fluid flow past a circular cylinder under the influence of an external magnetic field at high Reynolds numbers (Re).

The finite difference method is applied to solve the governing non‐linear Navier‐Stokes equations. First order upwind difference scheme is applied to the convective terms. The multigrid method with coarse grid correction is used to enhance the convergence rate. The defect correction technique is employed to achieve the second order accuracy.

A non‐monotonic behavior in separation angle when N≥5 and separation length when N≥3 is found with the increase of external magnetic field. The drag coefficient is found to increase with increase of N. The pressure drag coefficient, total drag coefficient and rear pressure are found to exhibit a linear dependence with N^0.5. The pressure Poisson equation is solved to find pressure fields in the flow region. It is found that the upstream base pressure increases with increase of external magnetic field while the downstream base pressure decreases with the increase of the external magnetic field.

The non‐monotonic behaviors in the separation angle and separation length at high Re are explained through pressure fields which are found first time for this problem. The linear dependence of the pressure drag coefficient, total drag coefficient and the pressure at rear stagnation point with N^0.5 is in agreement with experimental findings.

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.

A finite element method based on the Eulerian velocity correction method has been used to analyse the laminar natural convection in an annular cavity. Unsteady, incompressible, axisymmetric Navier‐Stokes equations have been made use of. Different radius ratios of the annular cavity have been considered to investigate the effect of the radius of curvature on the heat transfer coefficient.

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.