Sharif Digital Repository

We study numerically the heat transfer of steady laminar flow in a square cavity filled with electrically conducting fluids, in the presence of an external uniform magnetic field. Imposing a large temperature gradient between two opposite vertical walls, there are substantial temperature and density variations in the domain. The fluid is treated as an ideal gas. Indeed, high temperature gradient thermo-buoyant cavity flows result in natural convection flow domains with high Rayleigh number. To implement the temperature variation effect, the fluid properties, including the conductivity and viscosity coefficients, are considered to vary with temperature in accordance to the Sutherland's law.... 

We investigate the efficiency of a parallel direct simulation Monte Carlo (PDSMC) algorithm in solving the rarefied subsonic flow through a nanochannel. We use MPI library to transfer data between the processors. It is observed that PDSMC solver shows ideal speed up if sufficient workload is provided for each of processors. Additionally, this study shows that the computational time and speed up of the extended PDSMC solver do not depend (or slightly depend) on the number of cells. In contrary, increasing the total number of particles would result in a better efficiency of the PDSMC  

Benefitting from an analogy between compressible and incompressible governing equations, a novel dual-purpose, pressure-based finite-volume algorithm is suitably extended to simulate laminar mixing and reacting flows in low-Mach-number regimes. In our test cases, the Mach number is as high as 0.00326. Definitely, such low-Mach-number flows cannot be readily solved by either regular density-based solvers or most of their extensions. To examine the accuracy and performance of the extended formulation and algorithm, we simulate two benchmark cases including the mixing natural-convection flow in a square cavity with strong temperature gradients and the premixed reacting flow through annuli with... 

Mass-transfer rates to and from drops in liquid-liquid extraction processes are often reduced by the presence of contaminants. To design an industrial extractor, it is essential to consider this contamination effect in a quantitative manner. To achieve this goal, an experimental investigation was conducted on the mass transfer into single drops for n-butanol-succinic acid-water, as the recommended test system by the European Federation of Chemical Engineering (EFCE). The effects of anionic (sodium dodecyl sulfate, SDS), cationic (dodecyl trimethyl ammonium chloride, DTMAC), and nonionic (octylphenol decaethylene glycol ether, Triton X-100) surfactants on the hydrodynamic and mass-transfer... 

This paper presents an advanced grid adaptation strategy to be used by unstructured grid users. The idea behind this strategy originates from the need for automatic control of computational grids during iterative procedures utilized by fluid flow solvers. This strategy eliminates unnecessary grid computations by dividing the unstructured grid into active and inactive zones automatically. The active zones are extended automatically in order to capture the propagation of disturbances in solution domain. In this work, we focus to solve the grid deformations which are imposed in some portions of the main body and are propagated into computational domain during the iterative solutions. To achieve... 

A very simple and computationally low cost numerical algorithm is developed to generate a quasi-structured data structure for an unstructured grid. To achieve this purpose, the data structure in the matrices of an unstructured grid is classified to address the element layers and node lines in the computational domain. In this regard, elements and nodes of the unstructured grid are renumbered in a directional ordering-based strategy. The elements and nodes arrangement in each layer and line is accomplished in a unique direction around an interior object of the grid either clockwise or counterclockwise. Furthermore, a new searching scheme is introduced which guarantees a quick search inside... 

This paper presents the results of validation of an open source Direct Simulation Monte Carlo (DSMC) code for general application to rarefied gas flows. The new DSMC code, called dsmcFoam, has been written within the framework of the open source C++ CFD toolbox OpenFOAM. The main features of dsmcFoam code include the capability to perform both steady and transient solutions, to model arbitrary 2D/3D geometries, and unlimited parallel processing. Test cases have been selected to cover a wide range of benchmark examples from 1D to 3D. These include relaxation to equilibrium, 2D flow over a flat plate and a cylinder, and 3D supersonic flows over complex geometries. In all cases, dsmcFoam shows... 

Entropy generation due to heat transfer and fluid friction irreversibility has been investigated in a square cavity subjected to different side wall temperatures for compressible and incompressible natural convection flows. Based on the obtained velocity and temperature values, the distributions of local entropy generation, average entropy generation and average Bejan number are determined and compared for compressible and incompressible regimes. It is found that the entropy generated for compressible flow always is more than incompressible flow. The study is performed for Ra=10 4-10 8, j{cyrillic, ukrainian}=0.01(incompressible regime) and 0.6 (compressible regime), Ge=10 -5 and Pr=0.7  

A new approach has been presented to solve gas-liquid flow numerically in vertical pipes of air-lift pumps. To improve modelling, a new strategy has been employed with the capability of coupling the continuity and momentum equations and enforcing the role of pressure directly in the continuity equation. This is achieved via applying a novel scheme called the Physical Influence Scheme (PIS). The current finite volume solution is compared with other available numerical solutions. Indeed, they are in fair agreement. However, the present predictions are far superior to those obtained from an existing simple method, which is widely used in airlift pump modelling  

In this study, an arbitrary Lagrangian-Eulerian (ALE) approach is incorporated with a mixed finite-volume-element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non-stationary meshes. The method collects the advantages of both finite-volume and finite-element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical-influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving... 

A new approach has been presented to solve gas-liquid flow numerically in vertical pipes of air-lift pumps. To improve modelling, a new strategy has been employed with the capability of coupling the continuity and momentum equations and enforcing the role of pressure directly in the continuity equation. This is achieved via applying a novel scheme called the Physical Influence Scheme (PIS). The current finite volume solution is compared with other available numerical solutions. Indeed, they are in fair agreement. However, the present predictions are far superior to those obtained from an existing simple method, which is widely used in airlift pump modelling  

We study the natural convection heat transfer in a tilted square cavity with different tilt angles. The cavity is subject to a high gradient temperature resulting in high Rayleigh number flows. The fluid is air and is treated as an ideal gas. The flow is laminar. The fluid properties change with temperature variation using Sutherland's law. Because of imposing large temperature gradients to the two cavity opposite walls, there is substantial density variation in the domain. We use a novel non-Boussinesq algorithm to model the density variation fully. Therefore, the current results are considerably different from those obtained using the classical Boussinesq-based methods, which replace the... 

Because of its unique properties, graphene has attracted the attentions of many academic research groups and recently, the industry. One of the promising applications of the graphene is in micro/nano-sensors, e.g. using it as a pressure sensor. To use it in mechanical-based nano-sensors, it is very important to investigate the mechanical behavior of the nano-sized graphene sheet and its sensitivity to the medium changes applied on its faces. In this work, we use the molecular dynamics MD method and simulate the behavior of graphene sheet under differential water pressure influences. In this regard, a square straight monolayer graphene sheet is placed as a separator diaphragm between two... 

The current research uses an unstructured direct simulation Monte Carlo (DSMC) method to numerically investigate supersonic and subsonic flow behavior in micro convergent-divergent nozzle over a wide range of rarefied regimes. The current unstructured DSMC solver has been suitably modified via using uniform distribution of particles, employing proper subcell geometry, and benefiting from an advanced molecular tracking algorithm. Using this solver, we study the effects of back pressure, gas/surface interactions (diffuse/specular reflections), and Knudsen number, on the flow field in micronozzles. We show that high viscous force manifesting in boundary layers prevents supersonic flow formation... 

Rarefied gas flow in micro/nano electro mechanical systems (MEMS/NEMS) does not perform exactly as that in macro-scale devices. The main goal in this study is to investigate mixed subsonic-supersonic flows in micro/nano channels and nozzles and to provide physical descriptions on their behaviors. We use DSMC method as a reliable numerical tool to extend our simulation. It is because the DSMC provides accurate solution for the Boltzmann equations over the entire range of rarefied flow regime or Knudsen numbers. As is known, the appearance of oblique/normal shocks at the inlet of a channel or a nozzle adds to the complexity of internal flow field analyses. We found some very unique physical... 

Kinetic theory based numerical scheme such as direct simulation Monte Carlo (DSMC) and information preservation (IP) schemes properly solve micro-nano flow problems in transition and free molecular regimes. However, the high computational cost of these methods encourages the researchers toward extending the applicability of the continuumbased equations beyond the slip flow regime. In addition to correct velocity profile, the continuum-based equations should predict accurate mass flow rate magnitude. The secondorder velocity slip models derived from the kinetic theory provide accurate velocity profiles up to Kn=0.5; however, they yield erroneous mass flow rate magnitudes because the basic... 

This study is an extension of our previous work that was conducted on the aerodynamic characteristics of an ultra light weight airplane, which could be a human-powered airplane or an unmanned-aerial vehicle. One of the major aerodynamic tasks which should be done in this type of airplane is to reduce drag coefficient. The mainly effectual strategy for drag reduction is to contemplate on the components that make up the largest percentage of the overall drag and take into account that small improvements on large quantities can become significant aerodynamic improvements. Our experience showed that the use of light material in constructing human-powered airplanes and unmanned-air-vehicles has... 

We use the dissipative particle dynamics (DPD) method to simulate the non-Newtonian electroosmotic flow (EOF) through nanochannels. Contrary to a large amount of past computational efforts dedicated to the study of EOF profile, this work pays attention to the EOF of non-Newtonian fluids, which has been rarely touched in past publications. Practically, there are many MEMS/NEMS devices, in which the EOF behaviour should be treated assuming both non-continuum and non-Newtonian conditions. Therefore, our concern in this work is to simulate the EOF through nanochannels considering both non-Newtonian fluid properties and non-continuum flow conditions. We have chosen DPD as our working tool because... 

In this study, a jet propulsion JP combustor is studied numerically to investigate the combustor wall temperature influences on the soot characteristics emitted at its exhaust. There are a number of ways to control the combustor wall temperatures benefiting from different wall cooling technologies. Irrespective of using different high technology cooling systems, it is important to recognize how the wall temperature can affect the soot emission from one specific JP engine. Before examining the main combustor, it is important to assess the accuracy of the computational fluids dynamics (CFD) tool via solving a benchmark problem. In this regard, the predicted flame structure for the benchmark... 

A wind tunnel study was performed on the FFA-W-3-270 airfoil, which form a segment of a 1.25 MW wind turbine blade, to examine the effect of fixed roughness height and position using a zigzag tape boundary layer trip strip. Tests were conducted at a Reynolds number of 1×106 over a wide range of angles of attack. The zigzag tape, as an artificial roughness device, not only triggers a premature transition in the flow whereby laminar flow regimes change to turbulent, but also increases the momentum thickness of the turbulent boundary layer and change the airfoil camber. The 60° zigzag tape of 0.5 mm and 1 mm height was placed on the suction side of the airfoil at different chord wise locations....