Sharif Digital Repository

A two dimensional rectangular microchannel with circular micropillars was modeled in the presence of an electric field. Continuity and Navier-Stokes equations were solved along with convection-diffusion equation using finite element method. Reaction phenomenon was applied via a partial differential equation on the reaction surfaces and electric force was added as a source term to the transport equations. Velocity, concentration and electric potential distributions were obtained, with the aid of which, capture efficiency and average surface concentration of reaction surfaces were calculated. To ameliorate the reaction rate, different designs of reaction surfaces were investigated; the designs... 

Hydrodynamic instabilities at the interface of stratified shear layers could occur in various modes and have an important role in the mixing process. In this work, the linear stability analysis in the temporal framework is used to study the stability characteristics of a particle-laden stratified two-layer flow for two different background density profiles: smooth (hyperbolic tangent) and piecewise linear. The effect of parameters, such as bed slope, viscosity, and particle size, on the stability is also considered. The pseudospectral collocation method employing Chebyshev polynomials is used to solve two coupled eigenvalue equations. Based on the results, there are some differences in the... 

In this paper, an analytical solution for steady creeping motion of viscoelastic drop falling through a viscous Newtonian fluid is presented. The Oldroyd-B model is used as the constitutive equation. The analytical solutions for both interior and exterior flows are obtained using the perturbation method. Deborah number and capillary numbers are considered as the perturbation parameters. The effect of viscoelastic properties on drop shape and motion are studied in detail. The previous empirical studies indicated that unlike the Newtonian creeping drop in which the drop shape is exactly spherical, a dimpled shape appears in viscoelastic drops. It is shown that the results of the present... 

The density jump on an inclined surface is analyzed using an integral method by applying mass and momentum conservation equations. The jump occurs in a two-layered fluid flow in which the upper layer is stagnant and very deep. A relation is derived, which gives the conjugate depth ratio as a function of inlet densimetric Froude number, inlet concentration ratio, bed slope and entrainment. A set of experiments are performed to verify the relation. The theory and the measurements are in good agreement. The analysis reveals that increasing the surface inclination results in a decrease in the conjugate depth ratio. This analysis also shows that the densimetric Froude number just after the jump... 

The urine formation and excretion system have long been of interest for mathematicians and physiologists to elucidate the obscurities within the process happens in renal tissue. In this study, a novel three-dimensional approach is utilized for modeling the urine concentrating mechanism in rat renal outer medulla which is essentially focused on demonstrating the significance of tubule's architecture revealed in anatomic studies and physiological literature. Since nephrons and vasculatures work interdependently through a highly structured arrangement in outer medulla which is dominated by vascular bundles, a detailed functional unit is proposed based on this specific configuration.... 

The network of actin filaments in the lamellipodium (LP) of stationary and migrating cells flows in a retrograde direction, from the membrane periphery toward the cell nucleus. We have theoretically studied this phenomenon in the circular stationary (fully spread) cells. Adopting a continuum view on the LP actin network, new closed-form solutions are provided for the actin-retrograde-flow (ARF) in a polar coordinate system. Due to discrepancy in the mechanical models of the actin network in the ARF regime, solutions are provided for both assumptions of solid and fluid behavior. Other involved phenomena, including polymerizing machine at the membrane periphery, cytosol drag, adhesion... 

This study describes a multidimensional 3D/lumped parameter (LP) model which contains appropriate inflow/outflow boundary conditions in order to model the entire human arterial trees. A new extensive LP model of the entire arterial network (48 arteries) was developed including the effect of vessel diameter tapering and the parameterization of resistance, conductor and inductor variables. A computer aided-design (CAD) algorithm was proposed to efficiently handle the coupling of two or more 3D models with the LP model, and substantially lessen the coupling processing time. Realistic boundary conditions and Navier-Stokes equations in healthy and stenosed models of carotid artery bifurcation... 

Today, the flow of gas-solid, solid-liquid and liquid-liquid mixtures is broadly used in many industries such as slurry transportation, propulsion, dredging and power generation equipment. In this paper, single solid particle motion through free and forced vortices is analytically studied. The equations are solved for cases in which the drag, pressure gradient, added mass, buoyancy and weight forces are considered individually and simultaneously. Verification has been done for the value of ReP , which confirms the solution for the first t = 0.1 s. The results show that the most important force governing particle motion is the pressure gradient force  

Density currents (DCs) or gravity currents are driven by gravity in a fluid environment with density variation. Smoothed Particle Hydrodynamics (SPH) has been proved to have capabilities such as free surface modeling and accurate tracking of the immiscible-fluids interface that can be useful in the context of gravity currents. However, SPH applications to gravity currents have been limited to often-coarse simulations of high density-ratio currents. In this work, the SPH projection method is tried to solve currents with very low density-ratios (close to one), at a resolution, that captures the Kelvin-Helmholtz instabilities at the fluids interface. Existing implementations of the SPH... 

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

Turbidity current is produced when a particle-laden fluid flows under lighter ambient fresh fluid. The streaming of particle-laden fluid is called a turbidity current and this kind of current is an important mechanism for sediment transportation in lakes and oceans. In the present research, the main concentration is on the effect of obstacle with an isosceles right triangular cross section on the behavior of turbidity current. A series of laboratory experiments were carried out with various obstacle heights and different inlet densimetric Froude numbers. In each experiment, velocity profiles upstream and downstream of the obstacle were measured, using an acoustic Doppler velocimeter. Kaolin... 

A one dimensional model of the pulmonary arterial network along with systemic circulation was introduced to quantify the human pulmonary artery hemodynamics. A lumped parameter model of the right ventricle was used as the inlet boundary condition. A time varying elastance related the volume and pressure of the right ventricle was imposed in the model. By using a four element Windkessel model, blood flow rate throughout the right ventricle was in hand. An asymmetric structured tree was chosen as the outlet boundary conditions. Simulation was based on solving one-dimensional equations of conservation of mass and momentum by using finite volume method. Blood flow rate, pressure and velocity... 

The Fontan surgery is performed on patients with a single ventricle heart defect to prevent the combination of highlyoxygenated and poorly-oxygenated blood. Blood flow in total cavopulmonary connection (TCPC) which culminates an ordinary Fontan operation is practically steady-state but this flow is not appropriate for respiratory systems. This article investigates an approach in Fontan surgery that has been recently proposed in order to make the pulmonary blood flow pulsating. Moreover, for investigating the compliance of vessels and its effects on blood flow in TCPC, we have used the FSI (Fluid Structure Interaction) method as well as rigid wall assumption for comparison purposes. Our TCPC... 

The front part of a cell is divided to two regions called lamellum and lamellipodium (lamellipodial). This part plays an essential role for cell migration. Indeed, there are many protein filaments called actin in lamellum and lamellipodium, which induce the cell motion with polymerization in the leading edge of the cell. The actin filaments adhere to the extracellular matrix (ECM) by means of focal adhesions and they have contact by myosin motor proteins. The myosin motor proteins cause actin retrograde and anterograde flow exerted contractile stress on them. The focal adhesions exert frictional stress on the actin filaments. In this work, we developed a two-dimensional continuum model of... 

In this study, at first the cardiovascular system is modeled based on the 1D method and then the simulation of the reactive hyperemia experiment has been applied on the model. In this simulation, by applying a cuff at the brachial artery, the flow downstream of the cuff is occluded. Then with releasing the cuff immediately, a large amount of shear stress, about 4 times the basal amount, is applied downstream of the cuff and consequently to the endothelial cells in a very short moment. Considering a reported experimental transfer function between shear stress and vasodilation, the increase of the artery diameter due to the sudden increase of the shear stress is obtained. Finally, the... 

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

This study considers blood fl ow in Total Cavopulmonary Connection (TCPC) morphology, created in Fontan surgical procedure in patients with single ventricle heart disease. Ordinary process of TCPC operation reduces pulmonary blood fl ow pulsatility; because of right ventricle being bypassed. This reduction may limit the long term outcome of Fontan circulation. There is an idea stating that keeping Main Pulmonary Artery (MPA) partially open, would increase pulmonary fl ow pulsations. MPA gets closed in ordinary TCPC operation. The purpose of current study is to verify effects of Antegrade Flow (AF) coming through stenosed MPA on pulmonary fl ow pulsations, by means of Computational Fluid... 

In the present study, the motion of Newtonian and non-Newtonian liquid drops has been investigated experimentally. In order to investigate the effect of bulk fluid on drops, we have used water and air, as two fluids with different properties, and various industrial and biological applications. Image processing is utilized to analyze the images obtained by a high speed camera. The research has been separated into two parts. The first part has been devoted to the experiments in which air is the bulk fluid, and the second is related to the experiment carried out in water. The range of Reynolds number is, approximately, 50