Sharif Digital Repository

A novel technique for upscaling of detailed geological reservoir descriptions is presented. The technique aims at reducing both numerical dispersion and homogenization error generated due to incorporating a coarse computational grid and assigning effective permeability to coarse-grid blocks, respectively. In particular, we consider implicit-pressure explicit-saturation scheme where homogenization error impacts the accuracy of the coarse-grid solution of the pressure equation. To reduce the homogenization error, we employ the new vorticity-based gridding that generates a non-uniform coarse grid with high resolution at high vorticity zones. In addition, to control numerical dispersion, we use... 

This paper presents a new procedure for determining the fracture toughness of rock-like specimens using the diametric compression test with the center-cracked horseshoe disk (CCHD) method. Using finite element analysis, a dimensionless stress intensity factor was obtained and a formula was rendered for determining mode I fracture toughness. To evaluate the accuracy of the measurement results produced by the CCHD method, fracture toughness experiments were conducted on the same rock-like material using the notched Brazilian disk (NBD) method. The CCHD tests were simulated using a two-dimensional particle flow code for validation of the experimental results, and a great agreement between the... 

Gas permeation through nanoscale pores is ubiquitous in nature and has an important role in many technologies 1,2 . Because the pore size is typically smaller than the mean free path of gas molecules, the flow of the gas molecules is conventionally described by Knudsen theory, which assumes diffuse reflection (random-angle scattering) at confining walls 3-7 . This assumption holds surprisingly well in experiments, with only a few cases of partially specular (mirror-like) reflection known 5,8-11 . Here we report gas transport through ångström-scale channels with atomically flat walls 12,13 and show that surface scattering can be either diffuse or specular, depending on the fine details of the... 

This paper presents a novel coarse grid generation technique based on using vorticity as a measure of the impact of heterogeneity on flow. Vorticity is a maximum when the total flow is high and perpendicular to a large permeability gradient. Maps of vorticity were generated from single-phase flow simulations and used to generate coarse simulation grids from finely gridded geological models. The resulting grid was more refined in areas of high vorticity and coarser in areas of low vorticity. The method is first demonstrated on a simple five-layered model before being applied to three, 2D models of geologically realistic heterogeneity. The homogenized model generated from vorticity maps shows... 

In the present study, shock-disturbances interaction in high-speed compressible inviscid flow is simulated utilizing the weighted essentially non-oscillatory (WENO) scheme by employing the shock-capturing technique. For this aim, the two-dimensional Euler equations in strong conservative form are discretized by using the explicit third-order TVD Runge–Kutta scheme in time and the fifth-order WENO finite difference scheme in space. The main advantage of using the WENO scheme is its capability for properly solving the discontinuities in the domain without needing any artificial viscosity, limiter function or filter. Hence, this scheme is stable, and thus, it is suitable for simulating very... 

This paper presents a novel Gait Pattern Generator ((IPC) developed for the "Alice" social humanoid robot, which up to now lacked an appropriate walking pattern. Due to t he limitations of this robot, the proposed gate pattern generator was formulated based on a nine-mass model to decrease the modeling errors and t he inverse kinematics of the whole lower-body was solved in such a way that the robot remained statically stable during t he movements. The main challenge of this work was to solve t he inverse kinematics of a 7-link chain with 12 degrees of freedom. For this purpose, a new graphical-numerical technique has been provided using the definition of the kinematic equations of the robot... 

In this paper, the general analysis of two-dimensional steady-state thermal and mechanical stresses for a hollow thick infinite cylinder made of Functionally Graded Piezoelectric Materials (2D-FGPMs) is performed and developed. The general form of thermal, mechanical, and electrical boundary conditions is considered on the inside and outside surfaces. A direct method is used to solve the heat conduction equation and the non-homogenous system of partial differential Navier equations, using the complex Fourier series and the power law functions method. The material properties are assumed dependent on the radial and circumferential variables and are expressed as power law functions along the... 

In this paper, the general analysis of two-dimensional steady-state thermal and mechanical stresses for a hollow thick infinite cylinder made of Functionally Graded Piezoelectric Materials (2D-FGPMs) is performed and developed. The general form of thermal, mechanical, and electrical boundary conditions is considered on the inside and outside surfaces. A direct method is used to solve the heat conduction equation and the non-homogenous system of partial differential Navier equations, using the complex Fourier series and the power law functions method. The material properties are assumed dependent on the radial and circumferential variables and are expressed as power law functions along the... 

This paper presents a novel Gait Pattern Generator ((IPC) developed for the "Alice" social humanoid robot, which up to now lacked an appropriate walking pattern. Due to t he limitations of this robot, the proposed gate pattern generator was formulated based on a nine-mass model to decrease the modeling errors and t he inverse kinematics of the whole lower-body was solved in such a way that the robot remained statically stable during t he movements. The main challenge of this work was to solve t he inverse kinematics of a 7-link chain with 12 degrees of freedom. For this purpose, a new graphical-numerical technique has been provided using the definition of the kinematic equations of the robot... 

This article presents a two-dimensional simulation of the coal gasification in tapered fluidized bed gasifier. The effects of tapered angle, gasifier temperature, velocity of gasifying agent, and steam-to-air ratio on the gas compositions, lower heating value (LHV), and higher heating value (HHV) were examined. In order to find the appropriate operating conditions of coal gasification, carbon conversion efficiency (CCE) and cold gas efficiency (CGE) were also explored. It was found that with an increase of the gasifier temperature, CCE and CGE of the tapered gasifier diminishes. Increasing tapered angle results in a decrease of the LHV and HHV of the gas products, whereas the CCE of... 

A rigorous two-dimensional model is developed for simulating the operation of a less-investigated type steam reformer having a considerably lower operating Reynolds number, higher tube diameter, and non-availability of extra steam in the feed compared with conventional steam reformers. Simulation results show that reasonable predictions can only be achieved when certain correlations for wall to fluid heat transfer equations are applied. In all cases, strong radial temperature gradients inside the reformer tubes have been found. Furthermore, the results show how a certain catalyst loading profile will affect the operation of the reformer. © 2008 Canadian Society for Chemical Engineering  

In this study, an integrated 2D numerical model is applied to investigate the mechanical effects of the vegetation and soil type on slope stability. The developed model can assess the mechanical aspects of vegetation in slope stabilization. For this purpose, a case study of the Kheyrud forest located in northern Iran is considered as a real case site. Different scenarios including; three soil types (fine grain with low-, medium-, and high-strength) and two vegetation types (Maple and Common-ash) are assessed in the stability analysis (based on safety factor). The results confirm that the considered vegetation can prevent shallow landslide occasions, but has a limited impact on deep landslide... 

Oil trapping behavior during the pre-flush stage is critically important to evaluate the effectiveness of matrix acidizing for the oil well stimulation. In this study, the visco-capillary behavior of the two-phase flow in the pore-scale is analyzed to investigate the influence of wetting properties for a natural rock sample. A two-dimensional model, based on Cahn-Hilliard phase-field and Navier-Stokes equations, was established and solved using the finite element method. A stability phase diagram for log capillary number (Ca)-log viscosity ratio (M) was constructed and then compared with the reported experimental works. The maximum and minimum ranges of capillary number and viscosity ratio... 

This paper presents a case study and numerical simulations of a corner of a deep excavation in Tehran supported by soldier piles and ground anchors. This study focuses on the differences between 2D and 3D numerical modelling in estimating the wall deflection at the corner locations of the excavation. Furthermore, the performance of modelling with Mohr–Coulomb constitutive law was compared with the result of a hardening soil model. The modelling procedure was calibrated against a full-scale instrumented tieback wall at Texas A&M University and the monitoring data of the excavation project. The results indicated that the hardening soil model yields reasonable predictions of wall deflection in... 

A series of wind tunnel tests are performed to examine the flow field over a swept wing under various conditions. The wing has a laminar flow airfoil section, similar to those of the NACA 6-series. Static pressure distributions over the upper surface of the wing, in both chordwise and spanwise directions, are measured at different angles of attack. The data is employed to predict the transition point at each chordwise section. The skewness parameter of the pressure data shows that this factor drops to zero in the transition region. A comparison of the calculated transition point on the wing surface with that obtained from the 2D computational method shows reasonable agreement over a portion... 

In this paper, a two-dimensional numerical approach is used to study the effect of micro combustor height, 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 is used. The transient gas phase energy and species conservation equations result in an Advection-Diffusion-Reaction system that leads to two stiff systems of PDEs. In the present work, the computational domain is solved through the Strang splitting method, which is suitable for a nonlinear stiff system of PDEs. A revised boundary condition for the velocity equation is applied and its effect on the flow... 

This paper deals with the flexural instability of flexible spinning cylinders partially filled with viscous fluid. Using the linearized Navier-Stokes equations for the incompressible flow, a two-dimensional model is developed for fluid motion. The resultant force exerted on the flexible cylinder wall as the result of the fluid motion is calculated as a function of lateral acceleration of the cylinder axis in the Laplace domain. Applying the Hamilton principle, the governing equations of flexural motion of the rotary flexible cylinder mounted on general viscoelastic supports are derived. Then combining the equations describing the fluid force on the flexible cylinder with the structural... 

In nonlinear static (pushover) methods of analysis as an alternative to time history analysis, the capacity curve of the structure is established with respect to the roof displacement. Disproportionate increases in the roof displacement and even outright reversals of the higher modes can distort the capacity curve of the equivalent single degree of freedom system in these kinds of method, including MPA. To overcome this problem, recently, "Energy-Based" the Modal Pushover Analysis (Energy-Based MPA) method has been introduced. In this method, the absorbed energy and/or the external work in the pushover analysis is considered. Accordingly, the assessment of the Energy-Based MPA method is... 

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid... 

An effect-directed analysis (EDA) of fresh and artificially weathered (evaporated, photooxidized) samples of North Sea crude oil and residual heavy fuel oil is presented. Aliphatic, aromatic, and polar oil fractions were tested for the presence of aryl hydrocarbon receptor (AhR) agonist and androgen receptor (AR) antagonist, demonstrating for the first time the AR antagonist effects in the aromatic and, to a lesser extent, polar fractions. An extension of the typical EDA strategy to include an N-way partial least-squares (N-PLS) model capable of relating the comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry (GC × GC-TOFMS) data set to the bioassay data...