Sharif Digital Repository

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.... 

A weakly compressible smoothed particle hydrodynamics (WCSPH) method is used along with a new no-slip boundary condition to simulate movement of rigid bodies in incompressible Newtonian fluid flows. It is shown that the new boundary treatment method helps to efficiently calculate the hydrodynamic interaction forces acting on moving bodies. To compensate the effect of truncated compact support near solid boundaries, the method needs specific consistent renormalized schemes for the first and second-order spatial derivatives. In order to resolve the problem of spurious pressure oscillations in the WCSPH method, a modification to the continuity equation is used which improves the stability of... 

An explicit weakly compressible SPH method is introduced to study movement of suspended solid bodies in Oldroyd-B fluid flows. The proposed formulation does not need further stabilizing treatments and can be efficiently employed to study particulate flows with Deborah to Reynolds number ratios up to around 10. A modified boundary treatment technique is also presented which helps to deal with the movement of solid particles in the flow. The technique is computationally efficient and gives an improved evaluation of fluid-solid interaction forces.A number of test cases are solved to show performance of the proposed method in simulating particulate viscoelastic flows containing circular and... 

In this work, the Immersed Boundary Method (IBM) is adapted and implemented in the context of Smoothed Particle Hydrodynamics (SPH) method to study moving solid bodies in an incompressible fluid flow. The proposed computational algorithm is verified by solving a number of benchmark particulate flow problems. The results are also compared with those obtained using the same SPH scheme along with a direct solid boundary imposition technique  

Numerous studies have investigated driving equations used to predict the slip factor in centrifugal compressors so far. Inevitably, through these studies, the ow field characteristics have been simplified and the effects of the related parameters have been neglected. The present study, experimentally and numerically, investigates the slip phenomenon in a specific centrifugal compressor with complex blade curves and splitter blades, considering the main effective parameters such as the number of blades, exit angle, etc. To this end, a three-dimensional simulation of the viscous flow field of the compressor via an appropriate turbulence method was performed. In addition, an experimental study... 

A numerical model has been developed for studying the flow and heat transfer characteristics of single phase liquid flow through a microchannel. In this work the heat transfer characteristics of pressure driven and electroosmotic flow through microchannels have been studied. The governing equations are the Poisson-Boltzmann and Navier-Stokes equations which have been solved numerically using the standard Galerkin and the Mixed 4-1 finite element methods, respectively. Finally the energy equation is solved numerically using the Stream-wise Upwind Petrov Galerkin (SUPG) method. Two dimensional Poisson-Boltzmann equation was first solved to find the electric potential field and net charge... 

To investigate the factors affecting splashing in Peirce-Smith copper converters of Sarcheshmeh copper complex, South-East Iran, a laboratory scale model with dimensions one eighth of the actual converter was designed and assembled. The model was filled with water. Air was impinged onto the water. The sequence of events was subsequently recorded using high speed cameras. Factors affecting splashing such as blowing angle, volumetric flow rate of air, and the distance of the blowers from water surface were studied. Using a value of 2.65×10-11 for Morton dimensionless number, it was found that splashing may be reduced by 75% if the air speed in the tuyeres is reduced from 155 m/s to 30 m/s... 

A finite element solution procedure is presented for the simulation of transient incompressible fluid flows using triangular meshes. The algorithm is based on the artificial compressibility technique in connection with a dual time-stepping approach. A second-order discretization is employed to achieve the required accuracy in real-time while an explicit multistage Runge-Kutta scheme is used to march in the pseudo-time domain. A standard Galerkin finite element method, stabilized by using an artificial dissipation technique, is used for the spatial discretization. The performance of the proposed algorithm is demonstrated by solving a set of internal and external problems including flows with... 

In this paper, an explicit finite element based numerical procedure is presented for simulating three-dimensional inviscid compressible flow problems. The implementation of the first-order upwind method and a higher-order artificial dissipation technique on unstructured grids, using tetrahedral elements, is described. Both schemes use a multi-stage Runge-Kutta time-stepping method for time integration. The use of an edge-based data structure in the finite element formulation and its computational merits are also elaborated. Furthermore, the performance of the two schemes in solving a benchmark problem involving transonic flow about an ONERA M6 wing is compared and detailed solutions are... 

In this paper, flow of multi-component two-phase fluids in highly heterogeneous anisotropic two-dimensional porous media is studied using computational methods suitable for unstructured triangular and/or quadrilateral grids. The physical model accounts for miscibility and compressibility of fluids while gravity and capillary effects are neglected. The governing equations consist of a pressure equation together with a system of mass conservation equations. For solving pressure equation, a new method called Control Volume Distributed Finite Element Method (CVDFEM) is introduced which uses Control Volume Distributed (CVD) vertex-centered grids. It is shown that the proposed method is able to... 

Purpose: This paper aims to numerically study the compositional flow of two- and three-phase fluids in one-dimensional porous media and to make a comparison between several upwind and central numerical schemes. Design/methodology/approach: Implicit pressure explicit composition (IMPEC) procedure is used for discretization of governing equations. The pressure equation is solved implicitly, whereas the mass conservation equations are solved explicitly using different upwind (UPW) and central (CEN) numerical schemes. These include classical upwind (UPW-CLS), flux-based decomposition upwind (UPW-FLX), variable-based decomposition upwind (UPW-VAR), Roe’s upwind (UPW-ROE), local Lax–Friedrichs... 

This paper presents an extension of the multiscale finite volume (MsFV) method to multi-resolution coarse grid solvers for single phase incompressible flows. To achieve this, a grid one level coarser than the coarse grids used in the MsFV method is constructed and the local problems are redefined to compute the basis and correction functions associated with this new grid. To construct the coarse-scale pressure equations, the coarse-scale transmissibility coefficients are calculated using a new multi-point flux approximation (MPFA) method. The estimated coarse-scale pressures are utilized to compute the multiscale pressure solution. Finally a reconstruction step is performed to produce a... 

In this paper, a modified space-time method is presented for the simulation of convection-diffusion equations. The new method differs from the original space-time method in the sense that the weight functions for space and time are different. The performance of the proposed algorithm is studied for numerical simulation of incompressible immiscible two-phase flow in porous media. The governing equations consist of one conservation of mass equation for each phase, the Darcy law and one capillary-saturation correlation for the flow. By defining a global pressure, the governing equations lead to a system of nonlinear equations in terms of this global pressure, the velocity components and the... 

Discretization of spatial derivatives is an important issue in meshfree methods especially when the derivative terms contain non-linear coefficients. In this paper, various methods used for discretization of second-order spatial derivatives are investigated in the context of Smoothed Particle Hydrodynamics. Three popular forms (i.e."double summation","second-order kernel derivation", and"difference scheme") are studied using one-dimensional unsteady heat conduction equation. To assess these schemes, transient response to a step function initial condition is considered. Due to parabolic nature of the heat equation, one can expect smooth and monotone solutions. It is shown, however in this... 

Adirect numerical simulation approach is used to investigate the effective non-linear viscoelastic stress response of non-gap-spanning magnetic chains suspended in a Newtonian fluid. The suspension is confined in a channel and the suspended clusters are formed under the influence of a constant external magnetic field. Large amplitude oscillatory shear (LAOS) tests are conducted to study the non-linear rheology of the system. The effect of inertia on the intensity of non-linearities is discussed for both magnetic and non-magnetic cases. By conducting magnetic sweep tests, the intensity and quality of the non-linear stress response are studied as a function of the strength of the external... 

An improved incompressible smoothed particle hydrodynamics (ISPH) method is presented, which employs first-order consistent discretization schemes both for the first-order and second-order spatial derivatives. A recently introduced wall boundary condition is implemented in the context of ISPH method, which does not rely on using dummy particles and, as a result, can be applied more efficiently and with less computational complexity. To assess the accuracy and computational efficiency of this improved ISPH method, a number of two-dimensional incompressible laminar internal flow benchmark problems are solved and the results are compared with available analytical solutions and numerical data.... 

In this paper, a modified finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) extended for convection-diffusion problems with a non-linear flux function is introduced. Tracking schemes are discussed using viscous Burgers' equation. It is shown that in order to have smooth results, only the new time level values should be used in tracking process. Then, the proposed method is employed to study immiscible incompressible two-phase flows in porous media. Various one- and two-dimensional test cases involving internal sources and sinks are solved and accuracy of solution and performance of the method are investigated by comparing the results obtained using FVELLAM with those of... 

An analytical solution is presented to study the heat transfer characteristics of the combined pressure - electroosmotically driven flow in planar microchannels. The physical model includes the Joule heating effect to predict the convective heat transfer coefficient in two dimensional microchannels. The velocity field, which is a function of external electrical field, electroosmotic mobility, fluid viscosity and the pressure gradient, is obtained by solving the hydrodynamically fully-developed laminar Navier-Stokes equations considering the electrokinetic body force for low wall zeta potentials. Then, assuming a thermally fully-developed flow, the temperature distribution and the Nusselt... 

The present study employs Discrete Element Method (DEM) to establish a rheological model that relates the apparent viscosity of a granular material to shear rate, normal stress, and water saturation. In addition, a theoretical model was developed to determine water distribution and water-induced forces between particles for different saturations. The resulting forces were embedded in a 3D shear cell as a numerical rheometer, and a wet specimen was sheared between two walls. A power law rheological model was then obtained as a function of inertia number and saturation. It was found that up to a critical saturation, the apparent viscosity increased with saturation that was higher than that of... 

Purpose - The purpose of this paper is to investigate the performance of a specific class of high-resolution central schemes in conjunction with the black oil models for hydrocarbon reservoir simulation. Design/methodology/approach - A generalized black oil model is adopted, in which the solubility of gas in both oil and water and evaporation of oil are considered, leading to a system of equations prone to degeneracy. A computer code is generated and three test cases are solved to evaluate the performance of various schemes in terms of accuracy and discontinuity handling. Findings - It is shown that, although some of the central schemes are highly sensitive to the choice of...