Numerical modeling of general compressible multi-phase flows

Kalpaklı, Bora
In this thesis, some novel methods for solution of compressible, multi-phase flows on unstructured grids were developed. The developed methods are especially advantageous for interface problems, while they are also applicable to multi-phase flows containing mixtures as well as particle suspensions. The first method studied was a multi-dimensional, multi-phase Godunov method for compressible multi-phase flows. This method is based on the solution of a hyperbolic equation system for compressible multi-phase flows. There are several difficulties with this hyperbolic equation system due to non-conservative volume fraction equation and non-conservative terms also known as throttling therms existing in momentum and energy equations. Robust and accurate multi-dimensional discretization of these terms were derived based on Abgrall \cite{Abgrall} criterion. Next a new method based on discrete equations for multi-dimensional and multiphase problems on unstructured grids was developed. This method resolves all the problems associated with the non-conservative equations and terms. The high artificial numerical mixing of phase interfaces associated with available compressible schemes was resolved with a novel volume fraction differencing scheme. The developed differencing scheme used for volume fraction is the only scheme providing comparable resolution of the interfaces with tracking methods on multi-dimensional unstructured grids and very robust compared to other interface capturing methods studied in the related literature. The resulting methods provide ignorable numerical mixing of phase interfaces on ustructured solution grids while giving physically correct results for pressure and energy in contrast to other methods available in the literature. In addition to these solution methods, some special boundary conditions and preconditioning methods for low speed steady flows were applied. For high spatial resolution, combinations of linear reconstruction and Weighted Average Flux (WAF) methods were also applied in some problems.


Numerical simulation of advective Lotka-Volterra systems by discontinuous Galerkin method
Aktaş, Senem; Karasözen, Bülent; Uzunca, Murat; Department of Scientific Computing (2014)
In this thesis, we study numerically advection-diffusion-reaction equations arising from Lotka-Volterra models in river ecosystems characterized by unidirectional flow. We consider two and three species models which include competition, coexistence and extinction depending on the parameters. The one dimensional models are discretized by interior penalty discontinuous Galerkin model in space. For time discretization, fully implicit backward Euler method and semi-implicit IMEX-BDF methods are used. Numerical ...
An Efficient Numerical Solution Method for Reflectarrays of Varying Element Sizes
Ercil, Erdinc; Alatan, Lale; Aydın Çivi, Hatice Özlem (2015-12-01)
A novel method is developed for the efficient solution of reflectarrays of variable size elements. The method relies on using characteristic modes (CM) as macro basis functions and reusing the dominant CM of the resonant element for all elements in the reflectarray. This utilization leads to obtaining a reduced matrix system where the number of unknowns is drastically decreased. As far as the far field is concerned, accurate results are achieved even with a single CM. The accuracy is attained owing to prese...
Incompressible flow simulations using least squares spectral element method on adaptively refined triangular grids
Akdağ, Osman; Sert, Cüneyt; Department of Mechanical Engineering (2012)
The main purpose of this study is to develop a flow solver that employs triangular grids to solve two-dimensional, viscous, laminar, steady, incompressible flows. The flow solver is based on Least Squares Spectral Element Method (LSSEM). It has p-type adaptive mesh refinement/coarsening capability and supports p-type nonconforming element interfaces. To validate the developed flow solver several benchmark problems are studied and successful results are obtained. The performances of two different triangular ...
Design and analysis of a vertical axis water turbine for river applications using computational fluid dynamics
Demircan, Eren; Aksel, Mehmet Haluk; Pınarcıoğlu, Mehmet Melih; Department of Mechanical Engineering (2014)
The main purpose of this study is to design a Darrieus rotor type vertical axis water turbine using Computational Fluid Dynamics (CFD) in order to be used in river currents. The CFD modeling is based on two dimensional numerical solution of the rotor motion using commercial Unsteady Reynolds Averaged Navier-Stokes solvers, Ansys Fluent and CFX. To validate the two dimensional numerical solution, an experimental Darrieus rotor type water turbine from literature is studied and performance of several turbulenc...
Numerical simulation of the Çınarcık dam failure on the Orhaneli River
Bağ, Fırat; Bozkuş, Zafer; Department of Civil Engineering (2005)
This thesis analyzes the probable outcome of the fictitious failure of a dam under a set of pre-defined scenarios, within the framework of a case study, the case subject being the Cinarcik Dam located within Bursa Province of Turkey. The failure of the dam is not analyzed neither structural nor hydraulic-wise but is assumed to be triggered when certain critical criteria are exceeded. Hence, the analyses focus on the aftermath of the failure and strive to anticipate the level of inundation downstream of the ...
Citation Formats
B. Kalpaklı, “Numerical modeling of general compressible multi-phase flows,” Ph.D. - Doctoral Program, Middle East Technical University, 2013.