Sharif Digital Repository

Numerical study of two-dimensional supersonic hydrogen-air mixing layer is performed to investigate the effects of turbulence and chemical additive on ignition distance. Chemical reaction is treated using detail kinetics. Advection upstream splitting method is used to calculate the fluxes, and one-equation turbulence model is chosen here to simulate the considered problem. Hydrogen peroxide is used as an additive and the results show that inflow turbulence and chemical additive may drastically decrease the ignition delay in supersonic combustion  

In naturally fractured reservoirs, modeling of mass transfer between matrix blocks and fractures is an important subject during gas injection or contaminant transport. This study focuses on developing an exact analytical solution to transient tracer transport problem along a discrete fracture in a porous rock matrix. Using Gauss-Legendre quadrature, an expression was obtained in the form of a double integral which is considered as the general transient solution. This solution has the ability to account the following phenomena: advective transport in fractures and molecular diffusion from the fracture to the matrix block. Certain assumptions are made which allow the problem to be formulated... 

In this paper, compressible and incompressible flows through planar and axisymmetric sudden expansion channels are investigated numerically. Both laminar and turbulent flows are taken into consideration. Proper preconditioning in conjunction with a second-order accurate advection upstream splitting method (AUSM+-up) is employed. General equations for the loss coefficient and pressure ratio as a function of expansion ratio, Reynolds number, and the inlet Mach number are obtained. It is found that the reattachment length increases by increasing the Reynolds number. Changing the flow regime to turbulent results in a decreased reattachment length. Reattachment length increases slightly with a... 

Compressible turbulent flow through the abrupt enlargement in pipes is studied numerically by means of Advection Upstream Splitting Method (AUSM+-up). In low Mach numbers, a pressure correction equation of elliptic type is derived. This equation is compatible with the nature of governing equations and retrieves hyperbolic characteristic at higher Mach numbers. It is shown that the proposed numerical algorithm is computationally more efficient than the preconditioned density-based methods. The flow parameters such as reattachment length, pressure loss coefficient and wall shear stress are predicted. It is found that the loss coefficient of the compressible flow rises drastically with... 

A non-collimated gamma spectroscopy with fewer sampling (counting) and shorter sampling periods is developed and used to estimate the concentration profile of the radiotracers within the column. In this method, the radiotracer activity in the whole column is measured in a few different positions of the detector and then, by using an appropriate advection-dispersion model and the genetic algorithm, the concentration profile within the column is constructed. The concentration profile obtained by the new method fits very well to the experimentally obtained thin-section scanning results with RMSE less than 0.046. The uniqueness of the concentration profile obtained by the method was confirmed by... 

Combining the advection-diffusion equation approach with Monte Carlo simulations we study chaperone driven polymer translocation of a stiff polymer through a nanopore. We demonstrate that the probability density function of first passage times across the pore depends solely on the Péclet number, a dimensionless parameter comparing drift strength and diffusivity. Moreover it is shown that the characteristic exponent in the power-law dependence of the translocation time on the chain length, a function of the chaperone-polymer binding energy, the chaperone concentration, and the chain length, is also effectively determined by the Péclet number. We investigate the effect of the chaperone size on... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases... 

Urmia Lake, located in a closed basin in north-west Iran, is the largest lake (5000-6000 km2) in the Middle East. It is very saline with total dissolved salts reaching 200 g/l compared with a normal seawater salinity of about 35 g/l. The construction of a causeway, which was initiated in 1979 but then abandoned until the early 2000s, is near completion and will provide road access between the western and eastern provinces. The causeway has an opening 1.25 km long and divides Urmia Lake into a northern and southern basin and restricts water exchange. The flow and salinity regimes are affected by the presence of this new causeway, and there are concerns over the well being of the Artemia... 

The hydraulic efficiency of sedimentation basins is reduced by short-circuiting, circulation zones and bottom particleladen jets. Baffles are used to improve the sediment tank performance. In this study, laboratory experiments were used to examine the hydrodynamics of several baffle configurations. An accompanying numerical analysis was performed based on the 2-D Reynolds-averaged Navier-Stokes equations along with the k-ε turbulence closure model. The numerical model was supplemented with the volume-of-fluid technique, and the advection-diffusion equation to simulate the dynamics of particle-laden flow. Model predictions compared well with the experimental data. An empirical function was... 

A numerical code for supersonic intake design with a proper simulation of the normal and/or oblique shocks, boundary layer development, interaction of the shock and the boundary layer, as well as prediction of the flow separation is of great help to the designers. In this research, a numerical code is developed to solve the inner and outer flow fields of the intake and validated with various experimental tests. The intake is an axisymmetric external compression one. Roe scheme and new schemes, AUSM+-up (for all speed) and Advection Upstream Splitting Method with Pressure-Based Weight function (AUSMPW), are used to compute the convective fluxes. The original version of the AUSMPW scheme has... 

In this paper, a one-dimensional numerical approach is used to study the effect of various parameters such as micro combustor diameter, mass flow rate and external convection heat transfer coefficient on the temperature and species mass fraction profiles. A premixed mixture of H2-Air with a multi-step chemistry (9 species and 19 reactions) is used and thermal conductivity of the mixture is considered as a function of species thermal conductivity and temperature by using a set of new relations. The transient gas phase energy and species conservation equations result in an Advection-Diffusion-Reaction system (A-D-R) that leads to two stiff systems of PDEs, which can not be solved by... 

The duct inverse design in fluid flow problems usually involves finding the wall shape associated with a prescribed distribution of wall pressure or velocity. In this investigation, an iterative inverse design method for 2-D subsonic ducts is presented. In the proposed method, the duct walls shape is changed under a novel algorithm based on the deformation of a virtual flexible string in flow. The deformation of the string due to the local flow conditions resulting from changes in wall geometry is observed until the target shape satisfying the prescribed wall's pressure distribution is reached. The flow field at each step is analysed using Euler equations and the advection upstream splitting... 

In this work, an efficient bi-implicit strategy is suitably developed within the context of a hybrid finite volume element method to solve axisymmetric turbulent reactive flow in a model gas turbine combustor. Based on the essence of a control-volume-based finite-element method, the formulation benefits from the geometrical flexibility of the finite element methods while the discrete algebraic governing equations are derived through applying the conservation laws to discrete cells distributed in the solution domain. To enhance the efficiency of method, we extend the physical influence upwinding scheme to cylindrical coordinates. This extension helps to improve the advection flux... 

In this paper, transient depressurization of high pressure pipelines containing initially subcooled liquid is simulated numerically by using thermodynamic non-equilibrium and choking condition model. The numerical method relies on finite volume. The convective terms of cell boundaries are discretized by Advection Upstream Splitting Method (AUSM+ - up) with a proposal of partially implicit approach for source terms. Different void fraction correlations are applied to simulate two phase shock tubes as well as the depressurization process. By comparison between the present results and previous experimental data, the best void fraction correlation is introduced. The results indicate that the... 

Alumina nanofluids are one of the most useful nanofluids, especially for increasing the thermal conductivity. Due to importance of porous media in the improvement of heat transfer, this study investigates the transport and retention of gamma alumina/water nanofluid in the water-saturated porous media. For this purpose, alumina nanofluids were introduced to the porous media consisting of water-saturated glass beads possessing various pH values (4, 7 and 10) and different ionic strengths (0.001 M of KCl, CaCl2, AlCl3, K2SO4, CaSO4, Al2(SO4)3, K2CO3 and CaCO3). Then the break through curve of each experiment was drawn and modeled by combining classical filtration theory with... 

The performance of a gas centrifuge that is used for isotopes separation is dependent on the gas flow inside it. In this study, for modeling the UF6 gas flow, an Implicit Coupled Density-Based (ICDB) solver, was developed in OpenFOAM. To validate the ICDB solver, the gas flow within the rotor in total reflux state was compared with the analytical solution obtained by Onsager model and the numerical solution obtained by the Fluent software. The results showed that the ICDB solver had acceptable accuracy and validity. Also the computational efficiency of Roe, AUSM (Advection Upstream Splitting Method) and AUSM+ up schemes were compared. The results showed AUSM+ up scheme is efficient. Then,... 

In this paper the development of a quasi-three-dimensional numerical model that can be used for quantifying groundwater inputs and associated contaminant discharged from coastal aquifers into the coastal zone at a regional scale is presented. The present model is called MODSharp. In order to handle problems at a regional scale, the sharp interface approach which is used for conceptualising seawater intrusion, is applied to this model. This model can be used for the simulation of groundwater flow and contaminant transport in layered coastal aquifers at a regional scale. The method of characteristics is used to solve the advection-dispersion equation, which governs contaminant transport in... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration...