Sharif Digital Repository

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

Two-dimensional stagnation-point flow of viscoelastic fluids is studied theoretically assuming that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary-layer theory is used to simplify the equations of motion which are further reduced to a single non-linear third-order ODE using the concept of stream function coupled with the technique of the similarity solution. The equation so obtained was solved using Chebyshev pseudo-spectral collocation-point method. Based on the results obtained in the present work, it is concluded that the well-established but controversial prediction that in stagnation-point flows of viscoelastic fluids the velocity inside the boundary layer may exceed... 

The stability of a viscoelastic plate strip, subjected to an axial load with the Kelvin–Voigt fractional order constitutive relationship is studied. Based on the classical plate theory, the structural formulation of the plate is obtained by using the Newton’s second law and the aerodynamic force due to the fluid flow is evaluated by piston theory. The Galerkin method is employed to discretize the equation of motion into a set of ordinary differential equations. To determine the stability margin of plate the obtained set of ordinary differential equations are solved using the Laplace transform method. The effects of variation of the governing parameters such as axial force, retardation time,... 

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

Net fluid flow of soap films induced by (ac or dc) electric fields in asymmetric frames is presented. Previous experiments of controllable soap film flow required the simultaneous use of an electrical current passing through the film and an external electric field or the use of nonuniform ac electric fields. Here a single voltage difference generates both the electrical current going through the film and the electric field that actuates on the charge induced on the film. The film is set into global motion due to the broken symmetry that appears by the use of asymmetric frames. If symmetric frames are used, the film flow is not steady but time dependent and irregular. Finally, we study... 

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  

Several lattice Boltzmann models for multi-phase flow have been developed, but few of them are capable of modeling fluid flows with high density ratio in the order of 1000. Therefore, an advanced chromodynamics, Rothmann-Keller (RK) type model is employed in current study, which can handle liquid-gas density ratio in the order of 1000 and viscosity ratio in the order of 100. Other distinctive characteristics of the proposed model are high stability, and capability of setting parameters such as surface tension independently. In spite of these benefits, the original RK model fails to model wetting tendency of the fluids. As a result, it is impossible to correctly simulate two-fluid phase flow... 

Among numerous methods which have been employed to reinforce the thermal efficiency in many systems, one is the thermal radiation which is a mode of heat transfer. Another way to improve the thermal efficiency is the utilization of the porous media. The present work includes the study of micropolar flow with allowance for thermal radiation through a resistive porous medium between channel walls. The governing coupled partial differential equations representing the flow model are transmuted into ordinary ones by using the suitable dimensionless coordinates, and then, quasi-linearization is employed to solve the set of relevant coupled ODEs. Effects of physical parameters on the flow under... 

Among numerous methods which have been employed to reinforce the thermal efficiency in many systems, one is the thermal radiation which is a mode of heat transfer. Another way to improve the thermal efficiency is the utilization of the porous media. The present work includes the study of micropolar flow with allowance for thermal radiation through a resistive porous medium between channel walls. The governing coupled partial differential equations representing the flow model are transmuted into ordinary ones by using the suitable dimensionless coordinates, and then, quasi-linearization is employed to solve the set of relevant coupled ODEs. Effects of physical parameters on the flow under... 

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

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

Proppants transport is an advanced technique to improve the hydraulic fracture phenomenon, in order to promote the versatility of gas/oil reservoirs. A numerical simulation of proppants transport at both hydraulic fracture (HF) and natural fracture (NF) intersection is performed to provide a better understanding of key factors which cause, or contribute to proppants transport in HF–NF intersection. Computational fluid dynamics (CFD) in association with discrete element method (DEM) is used to model the complex interactions between proppant particles, host fluid medium and fractured walls. The effect of non-spherical geometry of particles is considered in this model, using the multi-sphere... 

In this study, we simulate the motion and reformation of polymer chain in the nanoscale fluid flow motion of the DPD (Dissipative Particle Dynamics) solvent. The behavior of polymer chain through DPD solvent is studied for 2D and 3D considerations. We implement two body forces of Poiseuille flow and electroosmotic flow to the DPD fluid particles. In case of the electroosmotic flow force, we show that the movement of polymer chain via the electroosmotic phenomenon provides less dispersion than that of the Poiseuille flow for the same polymer chain movement  

A hybrid lattice Boltzmann and level set method (LBLSM) for two-phase immiscible fluids with large density differences is proposed. The lattice Boltzmann method is used for calculating the velocities, the interface is captured by the level set function and the surface tension force is replaced by an equivalent force field. The method can be applied to simulate two-phase fluid flows with the density ratio up to 1000. In case of zero or known pressure gradient the method is completely explicit. In order to validate the method, several examples are solved and the results are in agreement with analytical or experimental results. © 2008 Elsevier B.V. All rights reserved  

Numerical study of flow through random packing of non-overlapping spheres in a cylindrical geometry is investigated. Dimensionless pressure drop has been studied for a fluid through the porous media at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), and numerical solution of Navier-Stokes equations in three dimensional porous packed bed illustrated in excellent agreement with those reported by Macdonald [1979] in the range of Reynolds number studied. The results compare to the previous work (Soleymani et al., 2002) show more accurate conclusion because the problem of channeling in a duct geometry. By injection of solute into the system, the... 

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

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

In the present study, the applicability and accuracy of a cell-vertex finite volume method developed are assessed in simulating 2D fluid–structure interaction in inviscid compressible flows where the nonlinear phenomena exist in both the unsteady transonic fluid flows and the large nonlinear deformation of solid structures. The unsteady Euler equations are considered as the governing equations of the fluid flow in the arbitrary Lagrangian–Eulerian form and the large nonlinear deformation of the solid structure is considered to be governed by the Cauchy equations in the total Lagrangian form. Both the domains are discretized by a second-order central-difference cell-vertex finite volume... 

Pushing a fluid with a less viscous one gives rise to the well known Saffman-Taylor instability. This instability is important in a wide variety of applications involving strongly non-Newtonian fluids that often exhibit a yield stress. Here we investigate the Saffmann-Taylor instability in this type of fluid, in longitudinal flows in Hele-Shaw cells. In particular, we study Darcy's law for yield stress fluids. The dispersion equation for the flow is similar to the equations obtained for ordinary viscous fluids but the viscous terms in the dimensionless numbers conditioning the instability now contain the yield stress. This also has repercussions on the wavelength of the instability as it... 

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