Sharif Digital Repository

The formation damage (FD) caused by the invasion of drilling fluid severely affects reservoir performance during production. Most of the published research studies which address this type of FD have been carried out at the core or field scale. Thus, the main aim of the paper is to investigate the pore-scale mechanisms of FD induced by drilling fluids and their control with silica nanoparticles (NPs) using a microfluidic approach. The proper identification of the mechanisms of FD can lead to the proper selection of NP type and concentration as well as a suitable method to remediate FD. The micromodel was designed in a way to closely simulate the cross-flow at the wellbore surface. A... 

We consider a cooperative wireless communication network comprising two full-duplex (FD) nodes transmitting to a common destination based on slotted Aloha protocol. Each node has exogenous arrivals and also may relay some of the unsuccessfully transmitted packets of the other node. In this article, we find the optimal static policies of nodes in order to minimize the sum of the average transmission delays, while the average transmission delay of each node is constrained. The static policy of each node specifies the probability of accepting an unsuccessfully transmitted packet of the other node and how the node prioritizes its transmissions. We show that in the optimal policies, just the node... 

This note considers the effect of different Darcy numbers on the onset of natural convection in a horizontal, fluid-saturated porous layer with uniform internal heating. It is assumed that the two bounding surfaces are maintained at constant but equal temperatures and that the fluid and porous matrix are in local thermal equilibrium. Linear stability theory is applied to the problem, and numerical solutions obtained using compact fourth order finite differences are presented for all Darcy numbers between Da = 0 (Darcian porous medium) and Da → ∞ (the clear fluid limit). The numerical work is supplemented by an asymptotic analysis for small values Da. © 2007 Elsevier Masson SAS. All rights... 

The optimal design of a casting feeding system is considered. The problem is formulated as the volume constrained topology optimization and is solved with the finite element analysis, explicit design sensitivity, and numerical optimization. In contrast to the traditional topology optimization where the objective function is defined on the design space, in the presented method, the design space is a subset of the complement of the objective function space. To accelerate optimization procedure, the nonlinear unsteady heat transfer equation is approximated with a Poisson-like equation. The feasibility of the presented method is supported with illustrative examples. © 2007 Springer-Verlag  

To describe immobilized cells in porous microcapsule membranes with straight pores, a novel model called corrugated parallel bundle model (CPBM) was utilized. In this model, a network was developed with 10 main pores each composing 10 pore elements. Cell growth kinetic in the network was examined using non-structural models. Effectiveness factor and pore plugging time were calculated by solving reaction-diffusion equation set via finite difference method. The findings revealed that diffusion coefficient for lower order reactions will create a lesser impact on the reduction of effectiveness factor. These findings also indicated that the use of such supporting carrier for cell immobilization... 

The numerical solution to the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock-fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously... 

The main objective of this work is to develop a novel moving-mesh finite-volume method capable of solving the seepage problem in domains with arbitrary geometries. One major difficulty in analysing the seepage problem is the position of phreatic boundary which is unknown at the beginning of solution. In the current algorithm, we first choose an arbitrary solution domain with a hypothetical phreatic boundary and distribute the finite volumes therein. Then, we derive the conservative statement on a curvilinear co-ordinate system for each cell and implement the known boundary conditions all over the solution domain. Defining a consistency factor, the inconsistency between the hypothesis... 

A two-dimensional Boussinesq-type model is presented accurate to O(μ)6 , μ = h0/l0, in dispersion and all consequential order for non-linearity with arbitrary bottom boundary, where h0 is the water depth and l0 is the characteristic wave length. The mathematical formulation is an extension of (4,4) the Padé approximant to include varying bottom boundary in two horizontal dimensions. A higher order perturbation method is used for mathematical derivation of the presented model. A two horizontal dimension numerical model is developed based on derived equations using the Finite Difference Method in higher-order scheme for time and space. The numerical wave model is verified successfully in... 

Simulation of brittle regime machining of materials (such as ceramics) is often difficult because of the complex material removal mechanisms involved. In this study, the discrete element method is used to simulate the dynamic process for machining of slip-cast fused silica ceramics. Flat-joint contact model is exploited to model contacts between particles in synthetic discrete element method models. This contact model is suitable for modeling of brittle materials with high ratios (higher than 10) of unconfined compressive strength to tensile strength. The discrete element method has the ability to simulate initiation, propagation, and coalescence of cracks leading to chip formation in the... 

High-order accurate solutions of parabolized Navier-Stokes (PNS) schemes are used as basic flow models for stability analysis of hypersonic axisymmetric flows over blunt and sharp cones at Mach 8. Both the PNS and the globally iterated PNS (IPNS) schemes are utilized. The IPNS scheme can provide the basic flow field and stability results comparable with those of the thin-layer Navier-Stokes (TLNS) scheme. As a result, using the fourth-order compact IPNS scheme, a high-order accurate basic flow model suitable for stability analysis and transition prediction can be efficiently provided. The numerical solution of the PNS equations is based on an implicit algorithm with a shock fitting procedure... 

Recently various versions of alternating group explicit or alternating group explicit-implicit methods were proposed for solution of diffusion equation. The main benefits of these methods are: good stability, accuracy and parallelizing. But these methods were developed for 1D case and stability and accuracy were investigated for 1D case too. In the present study we extend the new group explicit method [R. Tavakoli, P. Davami, New stable group explicit finite difference method for solution of diffusion equation, Appl. Math. Comput. 181 (2006) 1379-1386] to 2D with operator splitting method. The implementation of method is discussed in details. Our numerical experiment shows that such 2D... 

A new parallel Gauss-Seidel method is presented for solution of system of linear equations related to finite difference discretization of partial differential equations. This method is based on domain decomposition method and local coupling between interfaces of neighbor sub-domains, same as alternating group explicit method. This method is convergent and number of iterations for achieving convergence criteria is near the original Gauss-Seidel method (sometimes better and sometimes worse but difference is very small). The convergence theory is discussed in details. Numerical results are given to justify the convergence and performance of the proposed iterative method. © 2006 Elsevier Inc.... 

A two-dimensional depth-integrated numerical model is developed using a fourth-order Boussinesq approximation for an arbitrary time-variable bottom boundary and is applied for submarine-landslide-generated waves. The mathematical formulation of model is an extension of (4,4) Padé approximant for moving bottom boundary. The mathematical formulations are derived based on a higher-order perturbation analysis using the expanded form of velocity components. A sixth-order multi-step finite difference method is applied for spatial discretization and a sixth-order Runge-Kutta method is applied for temporal discretization of the higher-order depth-integrated governing equations and boundary... 

An unconditionally stable fully explicit finite difference method for solution of conduction dominated phase-change problems is presented. This method is based on an asymmetric stable finite difference scheme for approximation of diffusion terms and application of the temperature recovery method as a phase-change modeling method. The computational cost of the presented method is the same as the explicit method per time-step, while it is free from time-step limitation due to stability criteria. It robustly handles isothermal and nonisothermal phase-change problems and is very efficient when the global temperature field is desirable (not accurate front position). The robustness, stability,... 

A new group explicit method for solution of diffusion equation is presented. This method is based on domain decomposition concept and using asymmetric Saul'yev schemes for internal nodes of each sub-domain and alternating group explicit method for sub-domain's interfacial nodes. This new method has several advantages such as: good parallelism, unconditional stability, fully explicit nature and better accuracy than original Saul'yev schemes. The details of implementation and proving stability are briefly discussed. Numerical experiments on stability and accuracy are also presented. © 2006 Elsevier Inc. All rights reserved  

In the present study a finite difference method has been developed to model the transient fluid flow and heat transfer. A single fluid has been selected for modeling of mold filling and The SOLA-VOF 3D technique was modified to increase the accuracy of simulation of filling phenomena for shape castings. The model was then evaluated with the experimental methods. Refereeing to the experimental and simulation results a good consistency and the accuracy of the suggested model are confirmed. © 2005 Published by Elsevier B.V  

This paper presents a model to analyze pull-in phenomenon and dynamics of multi layer microplates using coupled finite element and finite difference methods. Firstorder shear deformation theory is used to model dynamical system using finite element method, while Finite difference method is applied to solve the nonlinear Reynolds equation of squeeze film damping. Using this model, Pull-in analysis of single layer and multi layer microplates are studied. The results of pull-in analysis are in good agreement with literature. Validating our model by pull-in results, an algorithm is presented to study dynamics of microplates. These simulations have many applications in designing multi layer... 

The numerical solution of the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock fitting procedure is utilized in both the compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are... 

Large-stencil schemes which their spectral properties are acceptable in the vicinity of ω = π are analyzed for the first time. A machine independent model for evaluating the efficiency of generalized time-marching finite-difference algorithms over periodic domains is developed. This model which is based on operation count reveals that for small values of Total Computational Cost(TCC), the previous low-order small-stencil schemes are more efficient while for moderate TCC, the efficiency of optimized large-stencil schemes abruptly increases. This important result is the motivation for developing optimized large-stencil schemes. The current schemes are successfully implemented in a full... 

This study is concerned with numerical modeling of viscous surface wave motion using boundary element method (BEM). The equations of motion for thin boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to obtain the viscosity-related quantities such as wave damping rate, shear stress, and velocity distribution inside the boundary layer. The technique is presented for standing surface wave motion. An excellent agreement is obtained between the numerical predictions and the previous results. The extension to other free surface problems is straightforward. © 2005 Elsevier Ltd. All rights reserved