Sharif Digital Repository

Recently, a hybrid Lattice Boltzmann Level Set Method (LBLSM) for two-phase incompressible fluids with large density differences, in cases of negligible or a priori known pressure gradients, has been proposed. In the present work, the mentioned LBLSM method is extended to take into account pressure gradient effects. The lattice Boltzmann method is used for calculating velocities, the interface is captured by the level set function, and the surface tension is replaced by an equivalent body force. The method can be applied to simulate two-phase fluid flows with density ratios up to 1000 and viscosity ratios up to 100. In order to validate the method, the evolution and merging of rising bubbles... 

The particle shape is an important factor playing critical role in evaluation of the interactions between particles in high-concentration particle-fluid flows. In this paper, the well-known multisphere (MS) approximation approach and the novel rolling resistance approach are utilized to examine their performance in order to simplify the generalized shaped particle’s interactions within the framework of discrete element method (DEM) and computational fluid dynamics (CFD). The performance of two approaches are compared with the perfect particle’s shape geometry, for the limited cases of cubic-shaped and disk-shaped particle flows in a horizontal well drilling process as a reference scenario.... 

In this study, a coupled numerical approach based on Lattice Boltzmann Method (LBM) and Discrete Element Method (DEM) is employed for two-dimensional simulation of fluid flow in deformable particulate media comprising of movable circular particles. The developed LB-DE code is validated against the results of a bi-axial shear test as well as two well-known benchmark problems including settling of a circular particle under gravity force inside a viscous fluid, and motion of a neutrally buoyant particle released in a Poiseuille flow. The verified code is then utilized for simulation of "Sand Production" phenomenon which is of importance for oil producing wells in weakly cemented sandstone... 

Colloidal Gas Aphron (CGA) are finding increasing application in fields of science and engineering because of their distinctive characteristic. As interest in the application of CGA based fluids grows and in order to select the best procedure for using them in successful petroleum engineering operations, there is a need to gain a better understanding of the factors that affect their properties and behavior. This article discusses the rheological characterization, stability analysis and filtration properties of CGA based fluids for three bio-polymers and two ionic surfactant. The stability and filtration analysis were investigated with the static drain rate technique and API filtration tests,... 

Natural flyers have flexible wings, which deform significantly under the combined inertial and aerodynamic forces. In this study, we focus on the role of chord wise flexibility in 2D pitch and plunge motions. We derive the exact nonlinear 2D equations of motion for a flexible flapping wing with flying support. In achieving the closed-form equations, we use the exact strain field concerning considerable elastic deformations. After numerically solving the novel equations, we validate them in simulations with highly deformable wings. By coupling the derived equations of motion with fluid flow, we study the aerodynamic performance of the geometrically nonlinear flexible flapping wing. Through... 

The paper investigates the effects of fluid flow on the static and dynamic behaviors of electrostatically actuated carbon nanotubes using nonlocal elasticity theory. The influences of various parameters of fluid flow including fluid viscosity, velocity, mass and temperature on the mechanical behaviors of the carbon nanotubes under static and step DC voltages are studied using this theory. The results computed from the nonlocal elasticity theory are compared with those estimated using the classical elasticity theorem and the outcomes demonstrate the applicability of the electrostatically actuated carbon nanotubes as nano sensors and nano actuators in nanofluidic systems. The nanosystem can be... 

A comprehensive study was performed to analyze the unsteady laminar flow characteristics around a porously covered, a fully porous, and a solid squared section cylinder located in the middle of a plane channel. In order to simulate fluid flow inside porous media and porous–fluid interface accurately (minimizing modeling error), the porous region was analyzed in pore scale, using LBM. Additionally, to minimize the LBM-related compressibility error through the porous region, a multi-block multiple relaxation time lattice Boltzmann method (MRT-LBM) was used. Also, to decrease CPU time, a Navier–Stokes flow solver, based on finite volume method and SIMPLE algorithm, was coupled with MRT-LBM to... 

In this paper, a pressure splitting formulation is proposed for Weakly Compressible SPH (WC-SPH) method and its capability in the suppression of the spurious oscillations is studied by conducting a stability analysis. The proposed formulation is implemented within the framework of a consistent SPH method. The predictions from the theoretical analysis are verified by the results of numerical test-cases. This method is applied to the incompressible fluid flow around periodic array of circular cylinders. The accuracy and the convergence of the results are investigated for benchmark problems. The results are also compared with those of the conventional WC-SPH method. In a similar test-case, the... 

Fluid flows generated on soap films by non-uniform alternating electric fields are studied. Two parallel metal rods subjected to an AC voltage are placed perpendicular to the soap film, which is anchored in a dielectric frame. The fluid flow is generated by electrohydrodynamic induction. At very low signal frequencies there is induced surface charge, but there is no tangential electric field at the surface, so there is no force and no flow. Fluid flow is observed increasing the frequency, when there are both surface charge and tangential electric field. The flow velocity increases with decreasing thickness of the soap film  

The current study aims at introducing a 2D and fast-running code for the issues pertinent to design, analysis and safety in modular high temperature reactors. While the porous media approach is only applied to pebble bed type, the analysis in this paper covers both pebble bed and prismatic reactor. A time-dependent mass equation along with energy conservation equation for the cooling gas and a time-dependent energy conservation equation for the solid was solved. Appropriate series of constitutive equations (e.g. heat transfer coefficient, effective heat conductivity of solid, heat transfer coefficient, pressure drop etc.) has been recruited as well. In addition a finite-volume method is... 

The proper understanding of rheological characteristics of CGA based fluids is of crucial importance in determining the performance of the fluid, in order to maintain the most effective fluid properties for safe, efficient, and economical drilling operation. This paper presents a concise investigation on the effect of concentration of the three main components of a novel environmentally friendly lightweight CGA based drilling fluid, i.e., xanthan gum biopolymer, starch, and biosurfactant, to the Herschel-Bulkley rheological model parameters. The three parameters of Herschel-Bulkley model, i.e., yield stress, fluid consistency, and fluid flow index were calculated by fitting the experimental... 

A precise thermal-hydraulics model is of great importance for developing more effective designs of High Temperature Gas Cooled Reactors (HTGR). Recently, several advancements have been made in the methods of analysis of porous media which could be of significant value in the development of more precise and robust codes. The objective of this research is to incorporate some of the most recent improvements in the development of a new 2D program for thermal-hydraulics analysis of modular high temperature reactors. The program is mainly based on the solution of a coupled set of mass, energy and momentum conservation equations for the gas flow, along with the energy conservation equation in the... 

Production from gas condensate reservoirs requires precise determination of reservoir fluid properties along with their positive impact on real reservoir performance evaluation and fluid in place volume calculation. This fact is particularly important because liquid drop out phenomena occurs as a result of pressure drop due to fluid production and condensate remains in reservoir [Mohebzadeh, 2005. In this study the challenges and problems of fluid phase behavior simulation in southern Iranian retrograde gas condensate reservoirs is discussed. Most of the problems and challenges are inaccuracy of equation of states near the critical point, special conditions and produced fluid flow... 

Field evidence indicates that the thermal and chemical regimes in wellbore considerably affect the wellbore stability. This study presents a general coupled model for transport of solute, solvent and heat including their combined effects on the wellbore stability. Optimization of drilling fluid parameters is crucial for wellbore stability analysis particularly in high pressure-high temperature environments. The coupled effects of chemical potential and temperature gradients on fluid flow significantly change the pore pressure and stress around a borehole. The effects of various parameters such as mud weight, solute concentration gradient, shale properties, and temperature gradient on... 

During the well drilling process, particles are produced in different shapes. The shape of particles can influence the characteristics of particles transport process. The aim of this work is to analyze the effects of particle shape on the transportation mechanism. For this purpose, a three-dimensional model is prepared for simulation of particle transportation with spherical and non-spherical shapes, during deviated well drilling. The motion of particles and the non-Newtonian fluid flow are simulated via discrete element method and CFD, respectively. The two-way coupling scheme is used to incorporate the effects of fluid-particle interactions. Three different samples of non-spherical shapes... 

The selection of interface boundary conditions between porous-medium and clear-fluid regions is very important for the wide range of engineering applications. In this paper, the difference between two common types of fluid flow interfacial conditions between clear fluid and porous medium is analyzed in detail. These two types of fluid flow interfacial condition are stress-jump and stress-continuity conditions. The effects of porosity on these types of interface condition are studied. The results are presented for different Reynolds numbers in the range 1–40, porosity equal to 0.4 and 0.8 and Darcy number Da=5×10-4. In this study, the Darcy–Brinkmann–Forchheimer model is used to model the... 

This paper presents an advanced grid adaptation strategy to be used by unstructured grid users. The idea behind this strategy originates from the need for automatic control of computational grids during iterative procedures utilized by fluid flow solvers. This strategy eliminates unnecessary grid computations by dividing the unstructured grid into active and inactive zones automatically. The active zones are extended automatically in order to capture the propagation of disturbances in solution domain. In this work, we focus to solve the grid deformations which are imposed in some portions of the main body and are propagated into computational domain during the iterative solutions. To achieve... 

The Poiseuille flow through slit-like nanochannels is investigated using the nonequilibrium molecular dynamics simulations. To drive a dense flow through the channel, we use two self-adjusting vertical plates strategy. These plates force the liquid to flow through the nanochannel under adjustable inlet and outlet boundary conditions. Comparing with the dual-control-volume grand-canonical molecular dynamics method, the current strategy provides many advantages. The current strategy does not need particle insertion and deletion, therefore, the system dynamics would not be affected at all. Moreover, the number of particles in the simulation system is fixed due to inserting the two... 

A numerical study has been carried out to investigate the fluid flow structure and convective heat transfer due to a circular jet impinging on a rotating disk. The temperature distribution and convection heat transfer coefficient on the disk are calculated. Flow is considered to be steady, incompressible and turbulent. k-e RNG model is used to model the turbulent flow. Two new criteria are introduced and used to evaluate the performance of cooling process which are maximum temperature difference on the disk and the average temperature of the disk. The first parameter shows the uniformity of temperature distribution in the disk and the second shows the effect of both thermo physical... 

This work deals with the problem of controlling the outlet temperature of a tubular heat exchanger system by means of flow pressure. The usual industrial case is to try to control the outlet temperature by either the temperature or the flow of the fluid, which flows through the shell tube. But, in some situations, this is not possible, due to the fact that the whole of system coefficients variation cannot quite be covered by control action. In this case, the system behavior must precisely be modeled and appropriate control action needs to be obtained based on novel techniques. A new multiple models control strategy using the well-known linear generalized predictive control (LGPC) scheme has...