Sharif Digital Repository

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a... 

Complicated sedimentary processes control the spatial distribution of geological heterogeneities. This serves to make the nature of the fluid flow in the hydrocarbon reservoirs immensely complex. Proper modeling of these heterogeneities and evaluation of their connectivity are crucial and affects all aspects of fluid flow. Since the natural variability of heterogeneity occurs in a myriad of length scales, accurate modeling of the rock type connectivity requires a very fine scheme, which is computationally very expensive. Hence, this makes other alternative methods such as the percolation approach attractive and necessary. The percolation approach considers the hypothesis that a reservoir can... 

Most of the heavy oil reservoirs contain discontinuous shale which affects fluid flow through porous media as well as recovery efficiency during enhanced oil recovery processes. However, the role of shale geometrical characteristics (including orientation, length, discontinuity, and spacing of the shale) on oil recovery remains a topic of debate in the literature, especially during miscible injection of heavy oils and five-spot systems. Here, a series of hydrocarbon solvent injection tests have been performed on various five-spot glass micromodels containing barriers which are initially saturated with heavy oil under fixed flow rate conditions. Oil recoveries as a function of pore volumes of... 

Body force modeling is studied in the Generalized Newtonian (GN) fluid flow simulation using a single relaxation time lattice Boltzmann (LB) method. First, in a shear thickening Poiseuille flow, the necessity for studying body force modeling in the LB method is explained. Then, a parametric unified framework is constructed for the first time which is composed of a parametric LB model and its associated macroscopic dual equations in both steady state and transient simulations. This unified framework is used to compare the macroscopic behavior of different forcing models. Besides, using this unified framework, a new forcing model for steady state simulations is devised. Finally, by solving a... 

The density jump on an inclined surface is analyzed using an integral method by applying mass and momentum conservation equations. The jump occurs in a two-layered fluid flow in which the upper layer is stagnant and very deep. A relation is derived, which gives the conjugate depth ratio as a function of inlet densimetric Froude number, inlet concentration ratio, bed slope and entrainment. A set of experiments are performed to verify the relation. The theory and the measurements are in good agreement. The analysis reveals that increasing the surface inclination results in a decrease in the conjugate depth ratio. This analysis also shows that the densimetric Froude number just after the jump... 

In this paper, a coupled hydro-mechanical formulation is developed for deformable porous media subjected to crack interfaces in the framework of extended finite element method. Governing equations of the porous medium consist of the momentum balance of the bulk together with the momentum balance and continuity equations of the fluid phase, known as [InlineEquation not available: see fulltext.] formulation. The discontinuity in fractured porous medium is modeled for both opening and closing modes that results in the fluid flow within the fracture, and/or contact behavior at the crack edges. The fluid flow through the fracture is assumed to be viscous and is modeled by employing the Darcy law... 

In this paper, the crack growth simulation is presented in saturated porous media using the extended finite element method. The mass balance equation of fluid phase and the momentum balance of bulk and fluid phases are employed to obtain the fully coupled set of equations in the framework of u - p formulation. The fluid flow within the fracture is modeled using the Darcy law, in which the fracture permeability is assumed according to the well-known cubic law. The spatial discritization is performed using the extended finite element method, the time domain discritization is performed based on the generalized Newmark scheme, and the non-linear system of equations is solved using the... 

The paper presents the effects of fluid flow on the static and dynamic properties of carbon nanotubes that convey a viscous fluid. The mathematical model is based on the modified couple stress theory. The effects of various fluid parameters and boundary conditions on the pull-in voltages are investigated in detail. The applicability of the proposed system as nanovalves or nanosensors in nanoscale fluidic systems is elaborated. The results confirm that the nanoscale system studied in this paper can be properly applied for these purposes  

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a... 

SUMMARY: In this paper, a numerical model is developed for the fully coupled hydro-mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non-wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two-phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid... 

In recent years, some numerical models have been proposed to investigate the effects of the elastic wave such as ultrasonic on fluid flow behavior in porous media. Nevertheless, none of these models are applicable to the fractured reservoirs, especially when the fluid is a Bingham plastic. In this work, the model proposed by P. P. Iassonov and I. A. Beresnev (2003) for flow under exposure of elastic wave in nonfractured porous media is considered and used for development of a new model of steady state flow of a Bingham plastic fluid in fractured reservoirs. The Kazemi's "block and fracture" model assuming negligible vertical permeability in blocks is considered for modeling. In addition,... 

Complicated sedimentary processes control the spatial distribution of geological heterogeneities. This serves to make the nature of the fluid flow in the hydrocarbon reservoirs immensely complex. Proper modeling of these heterogeneities and evaluation of their connectivity are crucial and affects all aspects of fluid flow. Since the natural variability of heterogeneity occurs in a myriad of length scales, accurate modeling of the rock type connectivity requires a very fine scheme, which is computationally very expensive. Hence, this makes other alternative methods such as the percolation approach attractive and necessary. The percolation approach considers the hypothesis that a reservoir can... 

In this paper, the fluid-structure interaction problem in mechanical systems in which a high frequency vibrating solid structure interacts with the surrounding fluid flow is considered. Such a situation normally appears in many microelectromechanical systems like a wide variety of microfluidic devices. A different implementation of the residual-based variational multiscale flow method is employed within the arbitrary Lagrangian-Eulerian formulation. The combination of the variational multiscale method with appropriate stabilization parameters is used to handle the so-called small time step instability in the finite element analysis of the fluid part in the coupled fluid-structure interaction... 

In this article, the problem of simultaneous heat, mass, and momentum transfer in the developing region of tubular reactors with homogeneous chemical reaction and laminar power-law fluid flow under unsteady-state conditions has been solved. In this regard, the general governing equations were solved using finite-difference method and analyzed carefully. Moreover, the influences of various parameters and dimensionless numbers such as power-law index, heat of reaction, reaction order, and Damköhler number value on the numerical results were investigated. In addition, the numerical results obtained for the fully-developed velocity distribution of Newtonian fluid flow and Sherwood number value... 

Proton exchange membrane (PEM) fuel cell performance is directly related to the bipolar plate design and their channels pattern. Power enhancements can be achieved by optimal design of the type, size, or patterns of the channels. It has been realized that the bipolar plate design has significant role on reactant transport as well as water management in a PEM Fuel cell. Present work concentrates on improvements in the fuel cell performance by optimization of flow-field design and channels configurations. A three-dimensional, multi-component numerical model of flow distribution based on Navier-Stokes equations using individual computer code is presented. The simulation results showed excellent... 

Microfilters are commonly used to block undesirable particles in the fluid flows and to control the flow patterns in MEMS. The main purpose of this study is to understand the effect of gas type on density, pressure, Mach number, and velocity distributions of fluid flows through a microfilter. The Knudsen number is the slip flow regime passing through the microfilter. We use direct simulation Monte Carlo (DSMC) method to simulate the flow of nitrogen, helium, oxygen, air and methane passing through a specific microfilter. The geometry of microfilter is unique in all cases. Our results confirm that every gas performs a different performance passing through a specific microfilter, and that the... 

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

Fluid flow behavior in a porous medium is a function of the geometry and topology of its pore space. The construction of a three dimensional pore space model of a porous medium is therefore an important first step in characterizing the medium and predicting its flow properties. A stochastic technique for reconstruction of the 3D pore structure of unstructured random porous media from a 2D thin section training image is presented. The proposed technique relies on successive 2D multiple point statistics simulations coupled to a multi-scale conditioning data extraction procedure. The Single Normal Equation Simulation Algorithm (SNESIM), originally developed as a tool for reproduction of... 

Automatic control of mesh movement is mandatory in many fluid flow and fluid-solid interaction problems. This paper presents a new strategy, called reduced domain strategy (RDS), which enhances the efficiency of node connectivity-based mesh movement methods and moves the unstructured grid locally and effectively. The strategy dramatically reduces the grid computations by dividing the unstructured grid into two active and inactive zones. After any local boundary movement, the grid movement is performed only within the active zone. To enhance the efficiency of our strategy, we also develop an automatic mesh partitioning scheme. This scheme benefits from a new quasi-structured mesh data... 

Relative permeability curves have practical implications in petroleum reservoir simulations. Study of the effects of reservoir wettability, pore shape geometry, and viscosity ratio of flowing fluids on the relative permeabilities is of great importance in reservoir modeling. In this paper, lattice Boltzmann method (LBM) is employed for analyzing the two-fluid flow in rigid porous media. The developed LBM code proved to be a robust numerical tool for analyzing the factors that influence the relative permeabilities of two immiscible fluids flowing through porous media. The numerically derived relative permeability curves demonstrate that in neutrally wet reservoirs, the effect of viscosity...