Sharif Digital Repository

Wellbore stability is a main concern in drilling operation. Troublesome drilling issues are chemically active formations and/or high-pressure–high-temperature environments. These are mainly responsible for most of wellbore instabilities. Wellbore failure is mostly controlled by the interaction between active shales and drilling fluid in shale formations. The factors influencing this interaction consist of fluid pressure, temperature, composition of drilling fluid, and exposure time. In this paper, a non-linear fully coupled chemo-thermo-poroelasticity model is developed. At first, a fully implicit finite difference model is presented to analyze the problem, and then, it is verified through... 

The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

A high-order compact finite-difference total Lagrangian method (CFDTLM) is developed and applied to nonlinear structural dynamic analysis. The two-dimensional simulation of thermo-elastodynamics is numerically performed in generalized curvilinear coordinates by taking into account the geometric and material nonlinearities. The spatial discretization is carried out by a fourth-order compact finite-difference scheme and an implicit second-order accurate dual time-stepping method is applied for the time integration. The accuracy and capability of the proposed solution methodology for the nonlinear structural analysis is investigated through simulating different static and dynamic benchmark... 

The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

This paper investigates implementation of upwind compact implicit and explicit high-order finite difference schemes for solution of the level set equation. The upwind compact implicit and explicit high-order finite difference schemes are well-known techniques to descritize spatial derivatives for convection term in hyperbolic equations. Applying of upwind high-order schemes on the level set equation leads to less error and CPU time reduction compared to essential non-oscillatory (ENO), weighted essential non-oscillatory schemes (WENO), and even different particle level set methods. The results indicate the error based on area loss decreases drastically with applying high-order upwind,... 

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

Utilization of heat storage units in solar energy systems can resolve the challenge of fluctuation and uncertainty of the solar energy. Phase change materials (PCMs) are used as the storage media for solar energy storage systems. In this research, a system including of a solar collector and a PCM-based cascaded energy storage unit was numerically investigated. Air was used as the heat transfer fluid (HTF) and three paraffin-based materials (RT50, RT65, and RT80) were used as PCM for the energy storage unit. The investigated system mainly operates between 15 °C and 90 °C and considering different PCMs, the selected PCMs were appropriate. Paraffin-based PCMs also present acceptable thermal... 

The Adomian Decomposition Method (ADM) for solving the highly non-linear vorticity-stream function formulation of 2D incompressible Navier-Stokes equations has been implemented. The analysis is accompanied by numerical boundary conditions. Also, numerical simulation, using finite difference method (FDM), is performed for comparison purposes. The obtained results only for few terms of the expansion are presented. Because present software such as Mathematica/Maple can not calculate many terms (for example: up to 10 terms) of solution and then ADM approach of this problem is an open problem case  

This paper reports the design of an optimal controller to prevent suppressvertical vibration due to undesired out of plane excitations generated by environment or gripper during manipulation for a CMOS-MEMS Nano-Newton capacitive force sensor applied for biomedical applications. Undesired out of plane excitations generated by environment or gripper during manipulation is the most prevalent source of vertical vibration in this type of sensors. To suppress the vibrational movement a PZT 5A is used as actuation mechanism. Discrete element method DEM model and Modal analysis were used to find dominant natural frequencies and mode shape vectors. To eliminate out of plane excitation an optimal... 

This paper presents the mathematical modeling of single-phase natural circulation of the University of Genoa's rectangular loop (LOOP#1) by a computer program and using RELAP5 system code. The mass, momentum and energy conservation equations in transient form were solved numerically using the finite difference method. One-dimensional linear stability analysis was performed for the single-phase natural circulation loop and the numerical perturbation technique was used in this analysis. The Nyquist criterion was employed to find the stability map of the LOOP#1. The obtained transient results using the first order upwind scheme of the fluid temperatures in various sectors of the LOOP#1 are... 

The present research is devoted to theoretical study of the pull-in performance of double-sided and paddle-type NEMS actuators fabricated from cylindrical nanowire operating in the Casimir regime and in the presence of the centrifugal force. D'Alembert's principle was used to transform the angular velocity into an equivalent static, centrifugal force. Using the couple stress theory, the constitutive equations of the actuators were derived. The equivalent boundary condition technique was applied to obtain the governing equation of the paddle-type actuator. Three distinct approaches, the Duan-Adomian Method (DAM), Finite Difference Method (FDM), and Lumped Parameter Model (LPM), were applied... 

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

This paper presents a model to analyze contact phenomenon in microsystems actuated by ramp voltages, which has applications in frequency sweeping. First-order shear deformation theory is used to model dynamical system using finite element method, while finite difference method is applied to model squeeze film damping. The model is validated by static pull-in results. The presented hybrid FEMFDM model is utilized to compute values of contact time and dynamic behavior. Considering this model, effects of different geometrical and mechanical parameters on contact time are studied. The influence of imposing the additional reverse voltage on dynamic characteristics of the system is also... 

A hexagonal-structured reactor core (e.g. VVER-type) is mostly modeled by structured triangular and hexagonal mesh zones. Although both the triangular and hexagonal models give good approximations over the neutronic calculation of the core, there are some differences between them that seem necessary to be clarified. For this purpose, the neutronic calculations of a hexagonal-structured reactor core have to be performed using the structured triangular and hexagonal meshes based on box method of discretisation and then the results of two models should be benchmarked in different cases. In this paper, the box method of discretisation is derived for triangular and hexagonal meshes. Then, two 2-D... 

Neutron noise induced by propagating disturbances in VVER-type reactor core is addressed in this paper. The spatial discretization of the governing equations is based on the box-scheme finite difference method for triangular-z geometry. Using the derived equations, a 3-D 2-group neutron noise simulator (called TRIDYN-3) is developed for hexagonal-structured reactor core, by which the discrete form of both the forward and adjoint reactor dynamic transfer functions (in the frequency domain) can be calculated. In addition, both types of noise sources, namely point-like and traveling perturbations, can be modeled by TRIDYN-3. The results are then benchmarked in different cases. Considering the... 

Neutron noise induced by propagating thermal-hydraulic disturbances (propagation noise for short) in pressurized water reactors is investigated in this paper. A closed-loop neutron-kinetic/thermal-hydraulic noise simulator (named NOISIM) has been developed, with the capability of modeling the propagation noise in both Western-type and VVER-type pressurized water reactors. The neutron-kinetic/thermal-hydraulic noise equations are on the basis of the first-order perturbation theory. The spatial discretization among the neutron-kinetic noise equations is based on the box-scheme finite difference method (BSFDM) for rectangular-z, triangular-z and hexagonal-z geometries. Furthermore, the finite... 

In this paper, the development of a neutron noise simulator for hexagonal-structured reactor cores using both the forward and the adjoint methods is reported. The spatial discretisation of both 2-D 2-group static and dynamic equations is based on a developed box-scheme finite difference method for hexagonal mesh boxes. Using the power iteration method for the static calculations, the 2-group neutron flux and its adjoint with the corresponding eigenvalues are obtained by the developed static simulator. The results are then benchmarked against the well-known CITATION computer code. The dynamic calculations are performed in the frequency domain which leads to discarding of the time... 

The deformation and breakup dynamics of a compound ferrofluid droplet under shear flow and uniform magnetic field are numerically studied in this paper. Utilizing magnetic field provides the possibility to obtain better control over the compound droplet morphology and breakup in a simple shear flow. To solve the governing equations for interfaces motion and hydrodynamics, the conservative phase field lattice Boltzmann model is employed, and a finite difference approach is applied for calculating the magnetic field. To verify the accuracy of present simulations, the results are validated with those of four relevant benchmarks including liquid lens between two stratified fluids, three-phase... 

Fast fault detection (FD) and reconfiguration is necessary for fault tolerant power electronic converters in safety critical applications to prevent further damage and to make the continuity of service possible. The aim of this study is to minimize the number of the used additional voltage sensors in a fault tolerant three-phase converter. In this paper, first a practical implementation of a very fast FD scheme with reduced sensor number is discussed. Then, an optimization in this scheme is also presented to decrease the detection time. For FD, special time and voltage criterion are applied to observe the error in the estimated phase-to-phase voltages for a specific period of time. The... 

This paper gives details of the refinement and application of a two-dimensional horizontal model for rivers. An explicit finite-difference algorithm was used for solving the governing differential equations, which includes the conservation of mass and momentum to predict hydrodynamic parameters. The model includes different turbulence closure models - that is, constant eddy viscosity, Prandtl simple mixing length and Smagorinsky methods. An experimental programme was designed and carried out in a laboratory flume to measure the length of eddy produced at the entrance of the intake. Model predictions have been compared with experimental results for a lateral intake. The effect of different...