Sharif Digital Repository

A fast and efficient reduced order formulation is presented for the first time to study dynamics and stability of conical/cylindrical shells with internal fluid flows. The structural and fluid formulations are developed based on general assumptions to avoid any deficiency due to modeling. Their respective solutions and the final solution to the coupled field problem are also developed in a way to be capable of capturing any desirable set of boundary conditions. In addition to the flexibility provided by the solution methodology and generalization provided by the formulation, current solution proposes an additional advantage over others which is the minimal computational cost due to the... 

In this paper, a coupled hydro-mechanical formulation is developed for the simulation of interaction between hydraulic fractures and natural discontinuities within saturated porous media. The momentum balance equation of the bulk together with the momentum balance and continuity equations of the fluid phase are employed to obtain the fully coupled set of governing equations. The hydro-fracture is modeled for fluid flow using the Darcy law. The natural discontinuity on the other hand is modeled for both opening and closing modes which induce fluid flow and/or contact behavior along the crack edges. The discontinuity in the displacement field is incorporated by using the Heaviside enrichment... 

A sequential implicit numerical method based on discrete-fracture model and the Galerkin Finite Element method, for time-dependent coupled fluid flow and geomechanics problems in fractured subsurface formations is presented. Discrete-fracture model has been used to explicitly represent the fracture network inside porous media. The Galerkin Finite Element method with adaptive unstructured gridding is implemented to numerically solve the spatially three-dimensional and time-dependent problem. The presented method is validated with previously obtained solutions. Two problems are numerically solved by applying the presented methodology in a three-dimensional fractured petroleum reservoir under... 

Purpose - A numerical method is developed to capture the interaction of solid object with two-phase flow with high density ratios. The current computational tool would be the first step of accurate modeling of wave energy converters in which the immense energy of the ocean can be extracted at low cost. Design/methodology/approach - The full two-dimensional Navier-Stokes equations are discretized on a regular structured grid, and the two-step projection method along with multi-processing (OpenMP) is used to efficiently solve the flow equations. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. The full... 

In this study, an extrinsically enriched element free Galerkin (EFG) method is proposed for the thermo-hydro-mechanical simulation of saturated porous media. By taking advantage of partition of unity property of moving least square shape functions, weak discontinuities such as material interfaces are modeled using the Ridge enrichment function and impermeable strong discontinuities are simulated using the Heaviside function. Some guidelines are proposed for the selection of EFG numerical parameters to ensure the stability and accuracy of the results. Numerical examples are provided to illustrate the capability of the proposed approach for fully coupled THM analysis of discontinuous porous... 

In this study, a novel type of LHP is innovated, fabricated, and its performance is assessed by numerous experiments. Compared to conventional LHPs, this novel design has some modifications in its evaporator and reservoir configuration. This particular type has a simpler and less costly fabrication procedure compared to other LHPs and it yields acceptable performance. Additionally, a novel method is introduced to control the temperature distribution in the system. A steel ball is placed in the evaporator to actively control the operating temperature. The ball is moved by two magnets installed outside of the evaporator. Moreover, a steady-state one-dimensional mathematical model of the... 

COMSOL Multiphysics is a comprehensive simulation software environment for a wide range of applications. COMSOL has an interactive interface that facilitates the modeling procedure and allows an easy coupling of different physical processes. The Subsurface Flow module extends the COMSOL modeling environment to applications related to fluid flow in saturated and variably saturated porous media. COMSOL is increasingly used in the investigation of geophysical, hydrogeological and environmental phenomena. The main goal of this work is to explore the ability of COMSOL for simulating seawater intrusion (SWI) in fractured coastal aquifers. Numerical modeling of such a problem is of high interest as... 

In this paper, a three-dimensional simulation of fully coupled multiphase fluid flow and heat transfer through deforming porous media is presented in the context of EFG mesh-less method. Spatial discretization of the system of governing equations is performed using EFG and a fully implicit finite difference scheme is employed for temporal discretization. Penalty method is used for imposition of essential boundary conditions. The developed numerical tool is employed to simulate two problems of nuclear waste disposal and CO2 sequestration in deep underground strata. The obtained results demonstrate the capability and robustness of the developed EFG code. © 2018 Elsevier Ltd  

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In... 

Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a... 

The three-dimensional power law fluid flow through rough microchannels has been studied numerically to determine the effects of the topographic structures on the thermal and hydrodynamic characteristics of the system. Rectangular, triangular and sinusoidal element shapes have been considered in order to investigate the effects of roughness height, width, pitch and channel separation on the pressure drop and heat transfer. Uniform wall heat flux boundary condition has been applied for all the peripheral walls. The results indicate that the global heat transfer performance can be improved or reduced by the roughness elements at the expense of pressure head when compared with the smooth... 

Net fluid flow of soap films induced by (ac or dc) electric fields in asymmetric frames is presented. Previous experiments of controllable soap film flow required the simultaneous use of an electrical current passing through the film and an external electric field or the use of nonuniform ac electric fields. Here a single voltage difference generates both the electrical current going through the film and the electric field that actuates on the charge induced on the film. The film is set into global motion due to the broken symmetry that appears by the use of asymmetric frames. If symmetric frames are used, the film flow is not steady but time dependent and irregular. Finally, we study... 

Microfluidic devices have been widely used for biological and cellular studies. Microbioreactors for three-dimensional (3D) multicellular spheroid culture are now considered as the next generation in in vitro diagnostic tools. The feasibility of using 3D cell aggregates to form multicellular spheroids in a microbioreactor with U-shaped barriers has been demonstrated experimentally. A barrier array is an alternative to commonly used microwell traps. The present study investigates oxygen and glucose concentration distributions as key parameters in a U-shaped array microbioreactor using finite element simulation. The effect of spheroid diameter, inlet concentration and flow rate of the medium... 

Optical and thermal analysis of the most well-known solar concentrator system; parabolic trough collector (PTC) are investigated and analyzed. To evaluate performance of the PTC, four cities of Iran with different weather conditions are chosen as case studies. Effective parameters such as concentration ratio, incident angle correction factor, collector mass flow rate are considered. The main objective of this work is evaluation of the solar energy potential using PTC in under consideration cities with different climates. Numerical modeling of the analysis is done using MATLAB software. Simulation results shows that Shiraz, with an average annual thermal efficiency of 13.91% and annual useful... 

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

In this paper, an extended finite element method is presented for simulation of interaction between hydraulic fracturing and natural fractures in saturated porous media. The well-known u−p formulation is employed in order to obtain the fully coupled set of governing equations. Natural faults are modeled for both opening and closure modes where the fluid inflow and contact conditions are considered at the interface, respectively. The Darcy law is employed in conjunction with an aperture dependent permeability for the fracture channel to describe the interfacial inflow. The contact constraints of both the solid and fluid phases are imposed using the Penalty method. The Heaviside and modified... 

Successful designs for an injection process to improve reservoir recovery seek more accurate values of saturation dependent functions, including relative permeability curves. The ignorance on underlying heterogeneity (e.g. layering) might lead to inaccurate predicting of flow behavior in layered heterogeneous porous media. Analysis of experimental core flood results in composite core systems (the system) along with numerical simulation at the core scale can provide an insight into this problem, especially when the capillary effects are present. In this study unsteady state two phase displacements have been conducted on the system with different internal arrangements. The main contribution is... 

In this work, the electrokinetic properties of asphaltene particles have been investigated. Micro-electrophoresis method by applying DC electric field, was utilized to different mixtures containing asphaltene to determine its electric charge. It was observed that in the case of using n-heptane and its mixture with toluene (heptol), the asphaltene particles were showed to be positively charged however for toluene itself, they expressed no tendency toward the electrodes. While it is expected that larger asphaltene aggregates carry higher electric charge, the results contradictorily showed that they are mainly governed by gravity rather than electro-static force and that “aggregation” reduces... 

The algebraic multigrid (AMG) algorithm as a preconditioner to the Krylov subspace methods has drawn the attention of many researchers in solving fluid flow and heat transfer problems. However, the efficient employment of this solver needs experience, because users have to quantify several important parameters. In this work, we choose a hybrid finite-volume element method and quantify the optimum magnitudes for those parameters. To generalize our results, two sets of fluid flow governing equations, the thermobuoyant flow and confined diffusion flame, are studied and the optimum values are determined. The results indicate that the AMG can be very effective if a proper storage method is chosen... 

The dynamics of a two dimensional plane jet injected at the base of a step, parallel to the wall, in backward facing step flow geometry is numerically studied. The objective of this work is to gain insight into the dynamics of the igniter flow field in solid fuel ramjet motors. Solid fuel ramjets operate by ingestion of air and subsequent combustion with a solid fuel grain such as polyethylene. The system of governing equations is solved with a finite volume approach using a structured grid in which the AUSM+ scheme is used to calculate the convective fluxes. The Spalart and Allmaras turbulence model is used in these simulations. Experimental data have been used to validate the flow solver...