Numerical methods for multiphysics flow problems

Belenli Akbaş, Mine
In this dissertation, efficient and reliable numerical algorithms for approximating solutions of multiphysics flow problems are investigated by using numerical methods. The interaction of multiple physical processes makes the systems complex, and two fundamental difficulties arise when attempting to obtain numerical solutions of these problems: the need for algorithms that reduce the problems into smaller pieces in a stable and accurate way and for large (sometimes intractable) amount of computational resources to resolve all the physical scales. Although these two difficulties are often stated as separate issues, in practice they are quite related. The objective of this thesis is to advance state of the art in algorithms, and their better understanding through analysis, for two types of multiphysics problems: incompressible non-isothermal fluid flow, and magnetohydrodynamic flow. The first component of this thesis is to develop numerical algorithms that decouple the multiphysics systems of equations into smaller, easier to solve sub-problems. However, splitting up problems into components is well known to (sometimes dramatically) reduce accuracy and cause numerical instabilities. It will be rigorously proven that the decoupling algorithms proposed and studied herein are stable and accurate. Numerical tests are used to verify the stability and accuracy. The second component of the thesis is to construct the numerical scheme that use the Variational Multiscale (VMS) method to reduce the computational cost of these problems, by reducing the size of the smallest scale needing to be resolved. At the same time, the algorithm will decouple the VMS modeling/stabilization equations from the multiphysics system, and decouple the multiphysics system into its components. This thesis proposes such an efficient algorithm and rigorously proves it is stable and accurate, as well as giving guidance into picking the stabilization parameters. Numerical experiments verify the theoretical results, and reveal that the algorithm gives accurate solutions on coarse discretizations, i.e. with significantly less computational cost than the requirement to be resolved of the original (unstabilized) physical systems. Lastly, this thesis considers the notion of long-time stability for decoupling algorithms for multiphysics problems. It is quite common for stable numerical methods to be stable only over finite time intervals, and to produce numerical solutions that non-physically grow linearly or even exponentially with time, even when the true solution does not grow. Hence, it is desirable to use algorithms that are stable at all times, so that stability and accuracy can be preserved as long as possible in a numerical simulation. This thesis proves unconditional long time stability results for a particular class of linearized, second order methods for multiphysics problems and also for the usual incompressible Navier-Stokes equations.


Convergence acceleration based on convergence error estimation
Eyi, Sinan (2013-01-01)
New methods are developed for convergence error estimation and convergence acceleration in iteratively solved problems. The convergence error estimation method is based on the eigenvalue analysis of linear systems, but it can also be used for nonlinear systems. Newton's method is used to estimate the magnitude and the phase angle of eigenvalues. The convergence of iterative method is accelerated by subtracting convergence error from the iteratively calculated solutions. The performances of these methods are...
Computation of multi-passage cascade flows with overset and deforming grids
Tuncer, İsmail Hakkı (null; 1997-12-01)
An overset grid method is applied to the solution of single and multi-passage cascade flows with a compressible Navier-Stokes solver. C-type grids around individual blades are overset onto a Cartesian background grid. Overset grids are allowed to move in time relative to each other as prescribed by the oscillatory plunging motion. The overset grid method uses a simple, robust numerical algorithm to localize moving boundary points and to interpolate solution variables across intergrid boundaries. Computation...
Application of fully implicit coupled method for 2D incompressible flows on unstructured grids
Zengin, Şeyda; Tarman, Işık Hakan; Department of Engineering Sciences (2012)
In the subject of Computational Fluid Dynamics (CFD), there seems to be small number of important progress in the pressure-based methods for several decades. Recent studies on the implicit coupled algorithms for pressure-based methods have brought a new insight. This method seems to provide a huge reduction in the solution times over segregated methods. Fully implicit coupled algorithm for pressure-based methods is very new subject with only few papers in literature. One of the most important work in this a...
Approximate factorization using ACDI method on hybrid grids and parallelization of the scheme
Onay, Oğuz Kaan; Kurtuluş, Dilek Funda; Sezer Uzol, Nilay; Department of Aerospace Engineering (2013)
In this thesis study, a fast implicit iteration scheme called Alternating Cell Directions Imp licit method is combined with Approximate Factorization scheme. This application aims to offer a mathematically well defined version of the Alternating Cell Directions Implicit Method and increase the accuracy of the iteration scheme that is being used for the numerical solutions of the partial differential equations. The iteration scheme presented here is tested using unsteady diffusion equation, Laplace equation ...
Numerical method for optimizing stirrer configurations
Schafer, M; Karasözen, Bülent; Uludağ, Yusuf; YAPICI, KEREM; Uğur, Ömür (2005-12-15)
A numerical approach for the numerical optimization of stirrer configurations is presented. The methodology is based on a parametrized grid generator, a flow solver, and a mathematical optimization tool, which are integrated into an automated procedure. The flow solver is based on the discretization of the Navier-Stokes equations by means of the finite-volume method for block-structured, boundary-fitted grids with multi-grid acceleration and parallelization by grid partitioning. The optimization tool is an ...
Citation Formats
M. Belenli Akbaş, “Numerical methods for multiphysics flow problems,” Ph.D. - Doctoral Program, Middle East Technical University, 2016.