Sharif Digital Repository

In common catamaran vessels, demihulls are parallel to each other. In this paper, the total resistance of a catamaran vessel with non-parallel demihulls is investigated experimentally and numerically. Experiments are carried out at different Separation Ratios (S.R.), that is the ratio of fore to aft separation of the catamaran demihullsand also in two ratios of length to separation in amidships (L/Sm). The FLUENT solver, based on the Finite Volume Method (FVM), was used for numerical solution. Applying the VOF model, the free surface around the catamaran vessel and total resistance are calculated and compared with experimental results. Finally, the frictional resistance of the catamaran from... 

In this paper, the ball valve performance is numerically simulated using an unstructured CFD (Computational Fluid Dynamics) code based on the finite volume method. Navier-Stokes equations in addition to a transport equation for the vapor volume fraction were coupled in the RANS solver. Separation is modeled very well with a modification of turbulent viscosity. The results of CFD calculations of flow through a ball valve, based on the concept of experimental data, are described and analyzed. Comparison of the flow pattern at several opening angles is investigated. Pressure drop behind the ball valve and formation of the vortex flow downstream the valve section are also discussed. As the... 

We study numerically the heat transfer of steady laminar flow in a square cavity filled with electrically conducting fluids, in the presence of an external uniform magnetic field. Imposing a large temperature gradient between two opposite vertical walls, there are substantial temperature and density variations in the domain. The fluid is treated as an ideal gas. Indeed, high temperature gradient thermo-buoyant cavity flows result in natural convection flow domains with high Rayleigh number. To implement the temperature variation effect, the fluid properties, including the conductivity and viscosity coefficients, are considered to vary with temperature in accordance to the Sutherland's law.... 

A computational fluid dynamics code was developed to compute the flow inside and around a supersonic external compression axisymmetric intake. The code solves the Reynolds-averaged Navier-Stokes equations using an explicit finite volume method in a structured grid and uses the Baldwin-Lomax algebraic model to compute the turbulent viscosity coefficient. Experiments were performed to validate the predicted results and good agreements are achieved. In the next part of the research, a parametric study was undertaken using the designed base case at a constant Mach number of 2 and at 0° angle of attack. The effects of various important parameters such as free stream Mach number, spike deflection... 

The brain is one of the vital organs in the body. The main cerebral distribution center of blood flow in the brain is the circle of Willis (CoW). In more than 50% of healthy brains and in more than 80% of dysfunctional ones, at least one artery of the circle of Willis is absent or underdeveloped. These variations reduce the collateral flow availability and increase the risk of stroke and transient ischemic attack in patients with atherosclerosis. Thus it is essential to simulate the circle of Willis and investigate the effects of stenosis. In this work the systemic arteries along with the circle of Willis are simulated using the finite volume method and one-dimensional equations of... 

A vast amount of death in the world has been attributed to atherosclerosis. This disease causes plaque formation and finally lack of blood supply to an organ. The prominent aim of this study is proposing an accurate and simple model to investigate the process of arterial wall thickening. In order to investigate LDL (low density lipoprotein) accumulation in arterial wall which is considered the first stage of atherosclerosis, a four layer model for arterial wall consisting of endothelium, intima, IEL and media is presented. All layers are treated as homogenous porous media. The four-layer arterial wall model is the most powerful and reliable tool for modeling LDL transport within arterial... 

In this study the extrudate swell of polymer solutions is estimated using the microstructure of polymer molecules. When a flexible polymer chain goes through a narrow die shear stress exerting on the chain will cause the polymer chain to be stretched along the flow direction. After emerging from die all external stresses vanish immediately and the chains tend to recover their previous state due to elastic recovery. This phenomenon will results in a gradual increase in extrudate diameter and this is used as the key idea for estimating swell ratio. A Giesekus based conformational model was used in order to predict polymer chains microstructure everywhere in the domain. The resulting PDE set... 

This paper presents a new landslide-generated wave (LGW) model based on incompressible Euler equations with Savage-Hutter assumptions. A two-layer model is developed including a layer of granular-type flow beneath a layer of an inviscid fluid. Landslide is modeled as a two-phase Coulomb mixture. A well-balanced second-order finite volume formulation is applied to solve the model equations. Wet/dry transitions are treated properly using a modified non-linear method. The numerical model is validated using two sets of experimental data on subaerial and submarine LGWs. Impulsive wave characteristics and landslide deformations are estimated with a computational error less than 5 %. Then, the... 

Nature presents a variety of propulsion, maneuvering, and stabilization mechanisms which can be inspired to design and construction of manmade vehicles and the devices involved in them, such as stabilizers or control surfaces. This study aims to elucidate and compare the propulsive vortical signature and performance of a foil in two important natural mechanisms: flapping and undulation. Navier–Stokes equations are solved in an ALE framework domain containing a 2D NACA 0012 foil moving with prescribed kinematics. All simulations are carried out using a pressure-based finite volume method solver. The results of time-averaged inline force versus Strouhal number (St) show that in a given... 

Design of an appropriate ventilation system is one of the main concerns in the construction of subways. In this paper, thermal comforts of side-platform stations are investigated by numerical methods. For the numerical simulation, grids are generated by structured methods and governing equations are discretized by finite volume methods. In the numerical simulation, second order upwind and second order central difference scheme are used to consider the convection and diffusion terms of momentum and energy equations, respectively. Results of this study show that in term of thermal comfort, region next to the hallway and the middle part of the platform is the worst zones of the station. In... 

This paper introduces a well-balanced second-order finite volume scheme, based on the Q-scheme of Roe, for simulating granular type flows. The proposed method is applied to solve the incompressible Euler equations under Savage-Hutter assumptions. The model is derived in a local coordinate system along a non-erodible bed to take its curvature into account. Moreover, simultaneous appearance of flowing/static regions is simulated by considering a basal friction resistance which keeps the granular flow from moving when the angle of granular flow is less than the angle of repose. The proposed scheme preserves stationary solutions up to second order and deals with different situations of wet/dry... 

Analysis of thermal-hydraulics fluctuations in pressurized water reactors (e.g., local and global temperature or density fluctuations, as well as primary and charging pumps fluctuations) has various applications in calculation or measurement of the core dynamical parameters (temperature or density reactivity coefficients) in addition to thermal-hydraulics surveillance and diagnostics. In this paper, the thermal-hydraulics fluctuations in PWRs are investigated. At first, the single-phase thermal-hydraulics noise equations (in the frequency domain) are originally derived, without any simplifying assumptions. The fluctuations of all the coolant parameters, as well as the radial distribution of... 

In this paper, two adaptive refinement techniques are presented for enhancing the capability of the new kind of finite volume method called meshless finite volume (MFV) method in which the moving least squares (MLS) approximation technique is employed. The proposed algorithms perform by inserting new nodes in large error regions identified using the Zienkiewicz-Zhu (Z-Z) error estimator and the T-Belytschko (TB) stress recovery scheme. In the first method referred as CV-based adaptive refinement, the new nodes are inserted at the vertices of control volumes with higher errors. The second method, referred as GA-based adaptive refinement, contains two schemes where an adaptive refinement... 

This paper simulates Al2O3-water nanofluid flow and forced convection around a rotating circular cylinder. The governing parameters are Reynolds number (1 ≤ Re ≤ 100), solid volume fraction of nanoparticles (0 ≤ φ ≤ 0.05) and non-dimensional rotation rate (0 ≤ α ≤ 3). The simulations are performed to study the effects of mentioned parameters on the heat transfer rate and fluid flow characteristics. The governing equations including the continuity, momentum, and energy equations are solved with a finite volume method. It is observed that the reduction of heat transfer with increase in rotation rate is in the vicinity of 6.9% and 32% for Re = 5 and 100, respectively at φ = 0.05. Furthermore,... 

The World Health Organization (WHO) in 2013 reported that more than seven million unexpected losses every year are credited to air contamination. Because of incredible adaptability and expense viability of fibrous filters, they are broadly used for removing particulates from gasses. The influence of appropriate parameters, e.g., the fiber arrangement, solid volume fraction (SVF or α), fluid flow face velocity (mean inlet velocity), and filter thickness (Ix), on pressure drop and deposition efficiency are researched. Furthermore, to study the effects of variation of the laminar flow regime and fiber’s cross-sectional shape on the deposition of particles, only a single square fiber has been... 

An unstructured grid, finite-volume, and three-dimensional (3D) primitive equation ocean model has been developed to predict oceanic pollutant dispersions in depth and surface of the Caspian Sea for non-chemical reactions of dissolved constituents. The model consists of momentum, continuity, temperature, salinity, and density equations. Physical and mathematical closure has been achieved using Mellor and Yamada turbulent closure sub-models. Since determining a practical definition of salinity that enjoys acceptable accuracy is difficult; therefore, various definitions have been used in this work. A recent definition of salinity stated in the UNESCO Practical Salinity Scale of 1978, PSS78,... 

Entropy generation analysis for the Cu–water nanofluid flow through a heat exchanger tube equipped with perforated conical rings is numerically investigated. Frictional and thermal entropy generation rates are defined as functions of velocity and temperature gradients. Governing equations are solved by using finite volume method, and Reynolds number is in the range of 5000–15,000. The effects of geometrical and physical parameters such as Reynolds number, number of holes and nanoparticles volume fraction on the thermal and viscous entropy generation rates and Bejan number are investigated. The results indicate that the thermal irreversibility is dominant in most part of the tube. But it... 

The development of density or turbidity currents causes serious problems for environmental hydraulics in reservoirs. The stream entered to a reservoir can carry sediments, nutrients and chemicals as density or turbidity currents. The fate of sediment and other substances transported by the current depends on the characteristics of the turbidity current itself, i.e. the velocity of fluid, the amount of mixing with reservoir water and the rates of sediment deposition and resuspension. These are important factors for water quality in reservoirs. A two-dimensional, depth-averaged, finite-volume numerical model is developed to study density currents, driven by non-cohesive sediments. The model... 

Multiple square cross section jets into a cross flow at three different velocity ratios, namely 0.5, 1.0 and 1.5, have been computationally simulated, using the Large Eddy Simulation (LES) approach. The finite volume method is applied in the computational methodologies, using an unsteady SIMPLE algorithm and employing a non-uniform staggered grid. All spatial and temporal terms in the Navier-Stokes equations have been discretized using the Power-Law and Crank-Nicolson schemes, respectively. Mean velocity profiles at different X-locations are compared with the existing experimental and Reynolds Averaged Navier-Stokes (RANS) computational results. Although the RANS computations require much... 

In resin injection/compression molding (RI/CM), a preform often comprises layers of different fiber reinforcements. Each fiber reinforcement has unique through thickness and in-plane permeabilities as well as compressibility, creating a heterogeneous porous medium in the mold cavity. In the present article, numerical simulation is utilized to investigate the filling process of RI/CM in such a heterogeneous porous medium. The filling stage is simulated in a full three-dimensional space by using control volume/finite element method and based upon an appropriate filling algorithm. The flow in the open gap which may be present in the mold cavity is modeled by Darcy's law using an equivalent...