Sharif Digital Repository

Solid fuel ignition transient in dump combustor geometry is numerically studied. The objective of this work is to gain insight into the ignition of the fuel and its controlling parameters in back-step turbulent flow. 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 gas phase convective fluxes. The Spalart-Allmaras turbulence model is used in these simulations. First, the auto-ignition of solid fuel subject to an inlet oxidizer flow at elevated temperatures is considered and distinct intervals in ignition delay time are studied. We have demonstrated that for inlet oxidizer streams with high... 

Multi stage Pulse Tube Cryocoolers (PTCs) have found the interest of researchers due to its reliability, long life and absence of moving part. This thesis represents a comprehensive numerical simulation method to investigate one and two stage PTC. All components of the PTC are modeled employing the nodal analysis technique to discretize the mass, momentum and energy conservation equations. SUTPTC code has been developed in cryocooler laboratory of Sharif University of Technology (SUT) to analyze and optimize one and two stage cryocooler performance. The SUTPTC code has been validated with the existing experimental data. Employing the proposed code, the effect of precooling temperature, the... 

Complicated theories dominated toward interacting of coastal regions resulted from its nature. Besides of that interference of human activities for acquiring beneficial aspects of shorelines causes environmental considerations to be challenged. One of the main factors for resolving this challenge is to obtain adequate information from hydrodynamic process and sediment regime. In spite of scientific researches that carried out along past three decades for understanding the exact mechanism of along shore sediment transport, it was not fully developed. One of the main reasons of uncertainty for forecasting the sediment regime regarded not exact knowledge of how sediment act in different... 

Drop penetration into the capillary tube, as a classical problem, has been studied over 100 years. But there are a few studies that investigate forced impregnation of capillary tubes, which has major application in coating, inkjet printing and rain drop penetration into the soil, have been made. A comprehensive numerical investigation on millimeter droplet impact dynamics on a thick capillary tube with micrometer sized pore is conducted and validated against experimental data. This work is extended to oblique impact of drops into pores and normal impact on nanoscaled pores. We apply the finite volume numerical method to solve the time-dependent governing equations of continuity, momentum in... 

The attention to ice accretion subject and its implementation have recently become very serious in our country due to a substantial need to a further understanding of this phenomenon and improving the design and construction of anti-ice systems for the airplanes in-hand. The subject of this study is the numerical simulation of ice accretion at the surface of airplanes, which is a technical approach in aerospace engineering. In other words, this study simulates the details of ice growth at different sections of an airplane wing considering wide ranges of different affecting parameters and conditions. To achieve this, we numerically solve the governing equations for two-phase flow of air and... 

Designers consider formation of slug flow in pressurized pipes because of hazardous effects of this flow pattern. Therefore, researchers, using different methods, attempt to study formation of slug flow and its effects. For numerical simulation of this phenomenon, In spite of existence of commercial 2D or 3D two-phase models, numerical models used for the slug flow modeling are 1D. In this thesis, a 2D model has been developed to simulate the initiation and formation of slug flow and estimation of the effects in pressurized pipes. Most important assumptions used in this model are existence of turbulence, surface tension and air compressibility. Isothermal condition has been supposed for... 

Nowadays, microfluidic flow appears in many applications, such as medical, biological, and chemical industries. Where as, micromixing, which deals with mixing of microfluidic flow in micro channels, appears to be an important issue to many researchers. In such systems, molecular diffusion plays an important role. On the other hand, lattice Boltzmann method is a relatively new simulation technique for complex fluid systems and has become interesting to many researchers in computational physics. In this study, computational simulation of such mixing process, using LBM is the main objective. Different obstacle layouts inside a microchannel have been investigated. Chaotic advection and jet... 

The objective of this thesis is to use proper orthogonal decomposition (POD) method in numerical simulation of a hydrocarbon reservoir to decrease computational time.The algorithm basically consists of two projection steps: (1) projecting the governing equations onto a subspace spanned by a low number of global empirical basis functions extracted from previous time step solutions, and (2) solving the governing equations in this reduced space and projecting the solution back on the original, high dimensional one. After projection of governing equations to new subspace, because of much less dimensions of this subspace, the computational cost of numerical solution of the equations, which mostly... 

In this research, the numerical simulation of liquid jet ejecting for a series of elliptical orifices into gaseous phase with different aspect ratios, at Rayleigh flow regime and in a range of Weber numbers has been performed.
For simulations, the OpenFoam software was utilized, which its multiphase flow simulation solvers are based on volume of fluid (VOF) method. In order to achieve the most accurate results for axis-switching phenomenon and jet breakup length, the dynamic mesh refinement was used for all the examined cases. The results, which were validated with recent experimental and numerical works, indicate that jet breakup length raises by increasing Weber number from 20 to 300.... 

In this project, a response of a single well in a gas condensate reservoir has been investigated by compositional modeling and implicit numerical method. The most significant goals of this paper are as follows: 1- Examination of the flow produced in a gas condensate reservoir in constant pressure or mass rate. 2- Investigation of variables affecting the productivity of reservoir. 3- The study of effect of injection and gas cycling to the well on the production. A three-zone model has been used in this study as if there is only gas phase in the third and second zones. Condensation occurs since the pressure drops below the dewpoint pressure in the second zone; however, gas phase is only... 

In this study, a new type of metallic yielding damper called Telescopic damper is proposed and analyzed numerically. The damper comprises of three concentric steel tubes connected together by lead shear keys. The relative axial deformation of steel tubes under axial forces results in shear deformation and yielding in lead shear keys and dissipates energy. The inner steel tube has a lock which limits its movement to some certain extent and beyond this limit the outer tube activates. Therefore, a multi-linear force-deformation behavior with two different yield plateaus is expected for the introduced telescopic damper. The multi-level yield plateaus of damper makes it effective in both low and... 

Safety against fire and its potential destructions is provided by detection, protection and ventilation methods. Knowledge on fire behavior such as ignition and propagation along with heat and smoke release is essential for implementing these methods. Having such deep information, the fire can be detected at the correct time and be controlled or suppressed, with its heat and smoke removed efficiently, so that rescue services will be enhanced and fatalities and damages will be reduced.Numerical simulation of fire requires deep understanding of multiple phenomena including turbulent flows, conductive, convective and radiative heat transfer, large buoyancy forces, chemical reactions and... 

In the present work, the numerical solution of 2D inviscid compressible Magneto-hydrodynamic flow is performed by using the spectral difference (SD) method on quadrilateral grids. In this numerical method, similar to the discontinuous Galerkin (DG) and spectral volume (SV) methods, the concept of the discontinuous and high-order local representations is used to achieve conservation property and high-order accuracy. In the SD method, the test function or the surface integral is not involved and thus it has a simpler formulation than the DG and SV methods. In this numerical method, two sets of structured points, namely unknown points and flux points, are defined in each cell to support the... 

The application of internal combustion engine is vast in different context, such as transportation and power generation. Therefore to investigation of these objects and the importance of flame front development, a combustion chamber with constant volume is designed and manufactured. The chamber designing is based on the experimental knowledge and numerical simulation of chamber and glasses. The combustion chamber made of stainless steel and contained transparent quartz glasses in both side of it to provide the ability of using schlieren photography method and investigation of flame. The great advantage of present chamber is the capability of different type of combustion, including fuel and... 

In the present study, hydrodynamics and particle sedimentation of a rectangular basin is studied numerically. Results are also validated with experimental data. Two major goals are perused in this work: First comparing LES and RNG k-ϵ turbulence models and second, investigating the accuracy of simulation of particle sedimentation by solving the concentration formula.
In order to achieve these goals, a 2D sedimentation tank is simulated using the RNG k-ϵ model. Velocity and concentration profiles are also plotted. Comparing the results with experimental data reveals that simulating sedimentation by the mentioned method has acceptable accuracy and can be used to predict sedimentation. In... 

Investigation of reservoir performance and long-term production prediction are main goal of numerical simulation. Numerical simulation mainly consists of solving sets of non-linear equations developed by combination of continuity equation, equation of state and Darcy equation. Classical methods such as finite difference method (FDM), finite volume method (FVM) and finite element method (FEM) are widely used for solving sets of fluid flow equations that mentioned above. Researchers have used other methods like boundary element method recently. Each of these methods has their own restrictions. But the main source of error in these methods is due to numerical estimation of spatial derivatives... 

Cavitation can occur in many fluid systems such as pumps, nozzles, hydrofoils and submarine vehicles and therefore, numerical modeling of this phenomenon has a significant importance. In this study, the numerical simulation of the cavitating flows through the Euler/Navier-Stokes equations employing the interface capturing method associated with a barotropic state law is performed. The system of governing equations is discretized using a cell-centered finite-volume algorithm and the fluxes are evaluated using a central-difference scheme. To account for density jumps across the cavity interface, the numerical dissipation terms with suitable density and pressure sensors are used. Since... 

Numerical simulations of participating media need careful solution of the radiation transfer equation (RTE) along the fluid flow governing equations. Evidently, the accuracy of achieved solutions highly depends on the accuracy of radiation transfer calculations. Despite numerous efforts performed in the previous researches, the RTE calculation still faces with several challenges from both accuracy and computational cost perspectives. The main objective of this research is to present new contributions in overcoming such challenges. So, a finite-volume (FV) solver is suitably developed to solve the RTE numerically. This extended solver is then coupled with an existing FV flow solver. The... 

In the this study, a numerical simulation of two-dimensional incompressible sloshing flow in a rectangular tank with baffle(s), using Lattice Boltzmann Method (LBM) is presented. Finite difference LBM and two-phase "Lee" model were used. The potential form of intermolecular forces is utilized to guarantee the stability of the numerical scheme and the discretization of the solution domain is performed by a two dimensional structured grid. Two different distribution functions are applied to obtain pressure, momentum, and composition of the particles. Furthermore, the Boltzmann transport equation is discretized, by using standard D2Q9 method. The "dropl" test case is simulated by the present... 

Moving of living organisms appears in many interesting problems, e.g. the growth of bacteria colonies, tumor growth, wound healing, color patterns of animals and etc. There are many ways to model such problems and PDE theory is widely used to investigate these problems. In this thesis, we study two well-known classic models. First, macroscopic “Keller–Segel” model and then kinetic “Othmer–Dunbar–Alt” System. Since these models have a nice behavior in two dimensions that they don’t have in other dimensions, we propose a way to alter them such that they behave in this way in all dimensions. Also none of the known models have the suitable dynamics in one dimension, so our model has the property...