Sharif Digital Repository

The gravity currents on the inclined boundaries are formed when the inflow fluid has a density difference with the ambient fluid and a tangential component of gravity becomes the driving force. If the density difference arises from the suspension of particles, the currents are known as particle-driven density currents, in which the local density of the gravity current depends on the concentration of particles. A low Reynolds k-ε turbulence model is used to simulate three dimensional turbidity currents. Also a laboratory apparatus was built to study the 3D flow resulting from the release of particle laden density currents on a sloping surface in a channel of freshwater via a sluice gate and... 

The gravity currents on the inclined boundaries are formed when the inflow fluid has a density difference with the ambient fluid and a tangential component of gravity becomes the driving force. If the density difference arises from the suspended particles, the currents are known as particle-laden density currents, or turbidity currents in which the local density depends on the concentration of particles. A low Reynolds k- Turbulent model is used to simulate three dimensional turbidity currents. Also some laboratory tests were conducted to study the 3D flow resulting from the release of particle laden density currents on a sloping surface in a channel of freshwater via a sluice gate. Kaolin... 

A new characteristic-based method is developed and used for solving the mixed and forced convection problems. The nano-fluid flow with heat transfer is simulated with a novel characteristic-based scheme in closed domains with different aspect ratios. For this purpose, a FORTRAN code has been written and developed. Water as a pure fluid and water-titanium dioxide as a nano-fluid were considered. The governing equations are solved by the finite volume utilizing a characteristic-based scheme for the convective fluxes. The simulation is done at Grashof numbers from 100 to 104, Reynolds numbers from 100 to 1000, and volume fractions of nano-particles from 0% to 10%. Streamlines, isotherms,... 

One of the most important challenges in mobile Ad-hoc networks is simulation. The simulation of these networks based on the accessible simulation techniques, which are based on packet, is a demanding and time-consuming task. Moreover, with the complexity of the network and the increase in the numbers of the nodes, it is likely to take a long time. This is because these simulators, one by one, regard the acts of all moving packets in one part of the network, and process these acts as a series of events. Since the number of the events is high in this method, the simulation takes time. Nowadays, in wired networks, several methods have been suggested to lessen the time of simulation; one example... 

Scale precipitation due to the mixing of incompatible injected water with formation brine and its subsequent deposition in porous media is an unpleasant phenomenon in water injection projects that can lead to severe injectivity and productivity decline. As a result of the complexity of geochemical reactions, modelling scale precipitation and deposition is a challenge. This paper presents a coupled geochemical and fluid flow model to simulate reactive flow in porous media which models pressure difference increase resulting from scale formation during water injection into porous media. To simulate chemical reactions during scale formation and subsequent rock permeability decline, PHREEQC... 

We have investigated the usage of a Dielectric Barrier Discharge (DBD) plasma actuator to improve the aerodynamic performance of an offshore 6 MW wind turbine. By controlling the aerodynamic load combined with pitch angles of 2, 5, and 10 degrees, we studied the plasma actuator effect on the overall harvested power. Actuators were installed in single and tandem configurations in different chord-wise locations to find the optimum design. The improved phenomenological model developed by authors was used in an analysis to simulate the interaction of the electrostatic field, the ionized particles and the fluid flow. A design software was used to estimate the harvested power of the real 3D blade.... 

Despite great advancement in the lattice Boltzmann method and its application in fluid flow problems, there are still major restrictions in treating either the solution domains with complex boundaries or buoyant flow problems. The past experience shows that the heat equation is a source for instabilities which jeopardizes the stable solution of the lattice Boltzmann method in solving fluid flow problems with heat transfer. The instabilities Increase with increasing buoyant force strength. In this work, we suggest a new approach to overcome the restrictions through implementing the advantages of finite volume method in LBM. In this regard, the lattice Boltzznann equation is incorporated with... 

Most available numerical methods face problems, in the presence of variable topographies, due to the imbalance between the source and flux terms. Treatments for this problem generally work well for structured grids, but most of them are not directly applicable for unstructured grids. On the other hand, despite of their good performance for discontinuous flows, most available numerical schemes (such as HLL flux and ENO schemes) induce a high level of numerical diffusion in simulating recirculating flows. A numerical method for simulating shallow recirculating flows over a variable topography on unstructured grids is presented. This mass conservative approach can simulate different flow... 

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

Rock formations are always under in situ stresses due to overburden or tectonic stresses. Drilling a well will lead to stress redistribution around the well. Understanding such a stress redistribution, and adopting a proper failure criterion, play a vital role in predicting any potential wellbore failure. However, most of the published analytical models are based on assumptions that do not satisfy the boundary conditions during production, that is, when the well pressure is less than the pore pressure. This paper is aimed at the modeling of the stress regime around the wellbore through combining the poroelastic model with proper boundary conditions under different flow regimes. As a result,... 

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

Electrohydrodynamic (EHD) printing is a novel technology used for fabricating high-resolution part features from a wide range of materials. Due to the multiphysics dynamics and the multiphase nature of the microdroplet formation in the EHD printers, modeling of this phenomenon is complicated. In this paper, the formation of a droplet in an EHD printer—under a pulsed electrical field—is simulated using a new numerical model which couples the fluid flow, the electric field distribution and the movement of the electric charges under dynamic and transient conditions. The level-set method is applied to the entire multiphysics domain in order to study the formation of the droplet. The presented... 

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

The purpose of this work is to develop a physical simulation model for the purpose of studying the effect of various design parameters on the performance of a packaged liquid chiller. A computer model which simulates the steady-state cyclic performance of a vapor compression chiller is developed for the purpose of performing detailed physical design analysis of actual industrial chillers. The model can be used for optimizing design and for detailed energy efficiency analysis of packaged liquid chillers. The simulation model takes into account presence of all chiller components such as compressor, shell-and-tube condenser and evaporator heat exchangers, thermostatic expansion valve and... 

A numerical model has been developed for studying the flow and heat transfer characteristics of single phase liquid flow through a microchannel. In this work the heat transfer characteristics of pressure driven and electroosmotic flow through microchannels have been studied. The governing equations are the Poisson-Boltzmann and Navier-Stokes equations which have been solved numerically using the standard Galerkin and the Mixed 4-1 finite element methods, respectively. Finally the energy equation is solved numerically using the Stream-wise Upwind Petrov Galerkin (SUPG) method. Two dimensional Poisson-Boltzmann equation was first solved to find the electric potential field and net charge... 

In the present research a numerical investigation has been made to evaluate the results of analytical approach in the analysis of micro channel heat sink (MCHS) which is increasingly used in the cooling of high heat dissipating electronic devices. In this regard the effects of geometrical and flow parameters affecting the absolute thermal resistance of MCHS have been analyzed. The effect of using different fluids in the value of the thermal resistance is considered as well. The accuracy of analytical method in predicting heat transfer and fluid flow regimes in micro channel is checked in comparison with numerical method. The effects of hydraulic and thermal entrance lengths of these systems... 

Present investigation analyzes the issue of entropy generation in a uniformly heated microchannel heat sink (MCHS). Analytical approach used to solve forced convection problem across MCHS, is porous medium model based on modified Darcy equation for fluid flow and two-equation model for heat transfer between solid and fluid phases. Furthermore, closed form solution of velocity distribution is employed to capture z-direction velocity gradient of flow, which plays a salient role on entropy generation through fluid flow. Analytical expressions for total and thermal entropy generation number (stems from heat transfer), and Bejan number are derived and cast into dimensionless form using velocity... 

Diffusivity equation which can provide us with the pressure distribution, is a Partial Differential Equation (PDE) describing fluid flow in porous media. The quadratic pressure gradient term in the diffusivity equation is nearly neglected in hydrology and petroleum engineering problems such as well test analysis. When a compressible liquid is injected into a well at high pressure gradient or when the reservoir possess a small permeability value, the effect of ignoring this term increases. In such cases, neglecting this parameter can result in high errors. Previous models basically focused on numerical and semi-analytical methods for semi-infinite domain. To the best of our knowledge, no... 

In this paper, an extended finite element method is presented for simulation of interaction between hydraulic fracturing and natural fractures in saturated porous media. The well-known u−p formulation is employed in order to obtain the fully coupled set of governing equations. Natural faults are modeled for both opening and closure modes where the fluid inflow and contact conditions are considered at the interface, respectively. The Darcy law is employed in conjunction with an aperture dependent permeability for the fracture channel to describe the interfacial inflow. The contact constraints of both the solid and fluid phases are imposed using the Penalty method. The Heaviside and modified... 

A lot of laboratory experiments have been carried out to investigate each aspect of ultrasonic wave s role on fluid flow behavior through porous media. Despite all experimental works, little attention has been paid for modeling the ultrasonic wave influence on capillary imbibition, which is the main mechanism of production in fractured reservoirs. At this work the process of imbibition with and without applying ultrasonic waves is mathematically modeled by modification of piston-like model. In contrast to this model s assumption, in which mobility ratio is assumed to be constant, here permeability variation due to increase in water saturation is considered in numerical solution. To evaluate...