A numerical method (SIMPLE DIRK Method) for unsteady incompressible viscous flow simulation is presented. The proposed method can be used to achieve arbitrarily high order of accuracy in time-discretization which is otherwise limited to second order in majority of the currently used simulation techniques. A special class of implicit Runge-Kutta methods is used for time discretization in conjunction with finite volume based SIMPLE algorithm. The algorithm was tested by solving for velocity field in a lid-driven square cavity. In the test case calculations, power law scheme was used in spatial discretization and time discretization was performed using a second-order implicit Runge-Kutta method. Time evolution of velocity profile along the cavity centerline was obtained from the proposed method and compared with that obtained from a commercial computational fluid dynamics software program, FLUENT 6.2.16. Also, steady state solution from the present method was compared with the numerical solution of Ghia, Ghia, and Shin and that of Erturk, Corke, and Goökçöl. Good agreement of the solution of the proposed method with the solutions of FLUENT; Ghia, Ghia, and Shin; and Erturk, Corke, and Goökçöl establishes the feasibility of the proposed method.

A numerical method (SIMPLE DIRK Method) for unsteady incompressible viscous flow simulation is presented. The proposed method can be used to achieve arbitrarily high order of accuracy in time-discretization which is otherwise limited to second order in majority of the currently used simulation techniques. A special class of implicit Runge-Kutta methods is used for time discretization in conjunction with finite volume based SIMPLE algorithm. The algorithm was tested by solving for velocity field in a lid-driven square cavity. In the test case calculations, power law scheme was used in spatial discretization and time discretization was performed using a second-order implicit Runge-Kutta method. Time evolution of velocity profile along the cavity centerline was obtained from the proposed method and compared with that obtained from a commercial computational fluid dynamics software program, FLUENT 6.2.16. Also, steady state solution from the present method was compared with the numerical solution of Ghia, Ghia, and Shin and that of Erturk, Corke, and Goökçöl. Good agreement of the solution of the proposed method with the solutions of FLUENT; Ghia, Ghia, and Shin; and Erturk, Corke, and Goökçöl establishes the feasibility of the proposed method.

This thesis is concerned with analysis and implementation of Time discontinuous Galerkin method. Important part of it is constructing of algorithm for solving nonlinear convection-diffusion equations, which combines Discontinuous Galerkin method in space (DGFEM) with Time discontinuous Galerkin method (TDG). Nonlinearity of the problem is overcome by damped Newton-like method. This approach provides easy adaptivity manipulation as well as high order approximation with respect to both space and time variables. The second part of the thesis is focused on Time discontinuous Galerkin method, applied to ordinary differential equations. It is shown that the solution of Time discontinuous Galerkin equals the solution obtained by Radau IIA implicit Runge-Kutta method in the roots of right Radau Quadrature. By virtue of this relation, error estimates of the order higher by one than the standard order can be obtained in these points. Furthermore, almost two times higher order can be achieved in the endpoints of the intervals of time discretization. Finally, the thesis deals with the phenomenon of stiffness, which may dramatically decrease the order of the applied method. The theoretical results are verified by numerical experiments. Powered by TCPDF (www.tcpdf.org)