Sharif Digital Repository

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

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

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

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

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

In the present paper, conceptual duct shape design for kinetic energy extraction with hydrokinetic turbines is discussed. The goal is to find a single-passage axisymmetric geometry that holds stable flow with maximum kinetic energy flux at duct throat. For finding the optimum duct shape, the fluid flow was numerically simulated in a wedge shaped space with Flow-Simulation Software. In a multi-stage conceptual design, tabulated configurations were employed to study each geometrical characteristic separately. These include curvature of profile camber, trailing edge shape, profile tip shape, and duct exit cross sectional area. The revolved profile of each duct consists of a well constrained... 

Excess pressure drop induced by inertial effects limits the applicability of Darcy's law for modeling of fluid flow through porous media at high velocities. It is expected such additional pressure drop is influenced by pore/morphology of porous media. This work concerns with fundamental understanding of how throat curvature affects intrinsic properties of porous media at non-Darcy flow conditions using network modeling. Conical, parabolic, hyperbolic, and sinusoidal capillary ducts with three types of imposed anisotropy are used to construct the network in a more realistic manner. Solutions of one dimensional Navier-Stokes equation for incompressible fluid flow through converging/diverging... 

In this paper, the effect of several different parameters on the thermal resistance of a Closed Loop Pulsating Heat Pipe (CLPHP) has been investigated. These parameters include the working fluid, the inclination angle, the filling ratio and the heat influx. Also, the impact of using nanofluids with different nano-particle concentrations has been analyzed. It was observed that a CLPHP can increase the heat transfer up to 11.5 times compared to an empty pipe. Optimum performance for a system with the water-silver nanofluid was achieved at conditions of 50% filling ratio and 0.9 K/W of thermal resistance, and for the water-titanium oxide system, these optimal conditions were found to be 40%... 

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

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In... 

Condensate and water banking around gas condensate wells result in vital well deliverability issues. Wettability alteration of near wellbore region to gas wetting condition is known to be the most novel and the only permanent method, to improve condensate well productivity. In this work, a water based nanofluid is used to change the wettability of sandstone reservoir rocks from strongly liquid wetting to intermediate gas wetting condition. Static contact angle measurements demonstrated significant increase of liquid phase contact angle as a result of chemical treatment with SurfaPore M nanofluid. The characteristics of SurfaPore M adsorption on sandstone rock are quantified through kinetic... 

The semi-empirical Kozeny–Carman (KC) equation is the widely used equation for determining permeability of porous media. Recent studies have shown that KC coefficient (CKC ) is a function of porous media parameters. In this study, the relation between parameters of randomly generated porous media is investigated to improve permeability prediction. In particular, site percolation theory is applied to construct random porous media. The static parameters of porous media, including porosity and specific surface area, are evaluated from porous media structure, and dynamic parameters, tortuosity and permeability, are derived from the results of Lattice Boltzmann fluid flow simulation.... 

In this study, an extrinsically enriched element free Galerkin (EFG) method is proposed for the thermo-hydro-mechanical simulation of saturated porous media. By taking advantage of partition of unity property of moving least square shape functions, weak discontinuities such as material interfaces are modeled using the Ridge enrichment function and impermeable strong discontinuities are simulated using the Heaviside function. Some guidelines are proposed for the selection of EFG numerical parameters to ensure the stability and accuracy of the results. Numerical examples are provided to illustrate the capability of the proposed approach for fully coupled THM analysis of discontinuous porous... 

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In...