Sharif Digital Repository

In the present paper, a numerical model is developed based on a combination of the extended finite element method and an equivalent continuum model to simulate the two-phase fluid flow through fractured porous media containing fractures with multiple length scales. The governing equations involve the linear momentum balance equation and the flow continuity equation for each fluid phase. The extended finite element method allows for an explicit and accurate representation of cracks by enriching the standard finite element approximation of the field variables with appropriate enrichment functions, and captures the mass transfer between the fracture and the matrix. Due to the high computational... 

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

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

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

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

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

Thermopneumatic micropump is one type of positive displacement micropump, which has many applications due to its relatively large stroke volume, low working voltage, and simple fabrication in microscale. In this paper, a numerical study of heat transfer and fluid flow in a valveless thermopneumatically driven micropump is presented. For rectifying the bidirectional flow, a nozzle and a diffuser are used as the inlet and outlet channels of the chamber. Since the fluid flow is induced by the motion of a diaphragm, the numerical simulation includes fluid structure interaction, which requires applying a dynamic mesh. The domain of solution is divided into two sections; the actuator unit, which... 

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

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 comprehensive set of experiments was performed on weirs placed obliquely in a rectangular open channel. The main objectives were to investigate the effect of the weir angle relative to the approach flow direction on the behaviour and the hydraulic flow characteristics and to examine methods such as using upstream guide vanes and inclined aprons to increase the effective oblique weir length. Results indicate that by increasing the oblique angle, the effective length of the oblique weir increases significantly. Relationships are provided for the threshold submergence and the effective length of the oblique weir. Design recommendations are also presented for the use of guide vanes and... 

In this paper, a fluid-structure interaction model for stability analysis of shells conveying fluid is developed. This model is developed for shells of arbitrary geometry and structure and is based on incompressible potential flow. The boundary element method is applied to model the potential flow. The fluid dynamics model is derived by using an inflow/outflow model along with the impermeability condition at the fluid-shell interface. This model is applied to obtain the flow modes and eigenvalues, which are used for the modal representation of the flow field in the shell. Based on the mode shapes and natural frequencies of the shell obtained from an FEM model, the modal analysis technique is... 

Fluid flow in microchannels has some characteristics, which one of them is rarefaction effect related with gas flow. In the present work, hydrodynamically and thermally fully developed laminar forced convection heat transfer of a rarefied gas flow in two microgeometries is studied, namely, microannulus and parallel plate microchannel. The rarefaction effects are taken into consideration using first-order slip velocity and temperature jump boundary conditions. Viscous heating is also included for either the wall heating or the wall cooling case. Closed form expressions are obtained for dimensionless temperature distribution and Nusselt number. The results demonstrate that for both geometries,... 

This study concerns with the microscopic and macroscopic fluid distribution and flow behavior during water alternating solvent (WAS) injection process to heavy oil using micromodel generated from thin section of a real rock which has rarely attended in the available literature. In this study, a one-quarter five-spot glass micromodel was deployed to examine the effect of flow media topology on microscopic displacements as well as macroscopic efficiency of WAS process. The micromodel was initially saturated with the heavy oil, and then the hydrocarbon solvent and water were injected alternately into it. The observations confirmed that WAS injection scheme is an effective method for the... 

Intracranial pressure (ICP) is derived from cerebral blood pressure and cerebrospinal fluid (CSF) circulatory dynamics and can be affected in the course of many diseases. Computer analysis of the ICP time pattern plays a crucial role in the diagnosis and treatment of those diseases. This study proposes the application of Linninger et al.s [IEEE Trans. Biomed. Eng. , vol. 52, no. 4, pp. 557565, Apr. 2005] fluidsolid interaction model of CSF hydrodynamic in ventricular system based on our clinical data from a group of patients with brain parenchyma tumor. The clinical experiments include the arterial blood pressure (ABP), venous blood pressure, and ICP in the subarachnoid space (SAS). These... 

Thermal conductivity, thermal dispersion, and structural effects in resin transfer molding (RTM) process are studied numerically. The injection part of the RTM process is modeled as a transport of resin flow through a fibrous porous medium in a long rectangular channel. The fluid flow is modeled using the Darcy-Brinkman-Forchheimer model and the heat transfer process using the energy equation based on local thermal equilibrium assumption. Both isotropic and anisotropic heat transfer in porous media are investigated. The governing equations are solved numerically for the isotropic heat transfer case and analytically for the anisotropic case. The numerical results are fitted to the available... 

The dynamics of a pressure regulator valve have been studied using the through Bondgraph simulation technique. This valve consists of several elements that can transmit, transform, store, and consume hydraulic energy. The governing equations of the system have been derived from the dynamic model. In solving system equations numerically, various pressure-flow characteristics across the regulator ports and orifices have been taken into consideration. This simulation study identifies some critical parameters that have significant effects on the transient response of the system. The results have been obtained using the MATLAB-SIMULINK environment. The main advantage of the proposed methodology... 

In the present research, weld pool geometry, thermal cycle, temperature and velocity fields during gas tungsten arc welding of aluminum alloys were predicted by solving three-dimensional equations of conservation of mass, energy and momentum under steady-state conditions. Welding experiments were then conducted on several samples with different thicknesses and chemical compositions. The geometries of the weld pools as well as the thermal cycles were measured. It is found that the calculated geometry of the weld fusion zone and the weld thermal cycles are in good agreement with the corresponding experimental results. In addition, the magnitude of the maximum velocities under different... 

In this paper, a fully coupled numerical model is developed for the modeling of the hydraulic fracture propagation in porous media using the extended finite element method in conjunction with the cohesive crack model. The governing equations, which account for the coupling between various physical phenomena, are derived within the framework of the generalized Biot theory. 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. By taking the advantage of the cohesive crack model, the nonlinear fracture processes developing along the fracture process zone are simulated. The spatial discretization...