Implementation and assessment of modern shock-capturing schemes for hypersonic viscous flows

Şahin, Çağatay
Spatially second order accurate Finite Volume Method (FVM) is the most preferred method in Computational Fluid Mechanics (CFD) with its acceptable results in short computation times. FVM's accuracy heavily relies on the particular numerical scheme with which the fluxes are evaluated. Despite the maturity of traditional flux schemes today, simulations of a viscous hypersonic flow are still challenging. Since these cases involve strong shock waves and viscous layers with non-linear gradients, appropriate shock capturing is necessary for robust and accurate solutions. In this thesis, modern shock-capturing schemes based on Liou's AUSM and Von Neumann and Richtmyer's artificial viscosity are studied and implemented in an in-house CFD solver. The accuracy and robustness of the implemented methods for viscous hypersonic flow problems are tested with well-known numerical experiments. The comparison of the obtained results with each other and experimental data in the literature is presented. The advantages and disadvantages of schemes among each other are identified.


Evaluation of Hypersingular Integrals on Non-planar Surfaces
Selcuk, Gokhun; Koç, Seyit Sencer (2014-05-16)
Solving electric field integral equation (EFIE) with Nystrom method requires accurate evaluation of hypersingular surface integrals since this method does not use divergence conforming basis and testing functions. The success of the method also depends on accurate representation of non-planar characteristics of the scattering object. In this study Hadamard finite part interpretation is used to evaluate hypersingular integrals over non-planar surfaces, which are represented by their Taylor series expansions....
Assessment and improvement of elementary force computations for cold forward rod extrusion
Ocal, M; Egemen, N; Tekkaya, AE (2005-06-01)
Two commonly used analytical force computation methods for cold forward rod extrusion are evaluated by means of precise finite element computations. The upperbound model by Avitzur based on the spherical velocity field and the model by Siebel based on a quasi-upper-bound solution are considered. It has been found that the pure deformation forces obtained by summing the ideal force and shear force terms deviate between +25% and -20% from the finite element solutions. Larger deviations, however, occur for the...
Parallel processing of two-dimensional euler equations for compressible flows
Doǧru, K.; Aksel, M.h.; Tuncer, İsmail Hakkı (2008-12-01)
A parallel implementation of a previously developed finite volume algorithm for the solution of two-dimensional, unsteady, compressible Euler equations is given. The conservative form of the Euler equations is discretized with a second order accurate, one-step Lax-Wendroff scheme. Local time stepping is utilized in order to accelerate the convergence. For the parallel implementation of the method, the solution domain is partitioned into a number of subdomains to be distributed to separate processors for par...
Applications of hybrid discrete Fourier transform-moment method to the fast analysis of large rectangular dipole arrays printed on a thin grounded dielectric substrate
Chou, HT; Ko, HK; Aydın Çivi, Hatice Özlem; ERTÜRK, VAKUR BEHÇET (2002-08-05)
Recently a discrete Fourier transform-method of moments (DFT-MoM) scheme was developed for fast analysis of electrically large rectangular planar dipole arrays, which has been shown to be very efficient in terms of number reduction of unknown variables and computational complexity. The applications of this DFT-MoM to treat dipole arrays printed on a grounded dielectric substrate are examined in this Letter. Numerical results are presented to validate its efficiency and accuracy. (C) 2002 Wiley Periodicals, ...
Ciydem, Mehmet; Koç, Seyit Sencer (2014-03-01)
Numerical methods in space-time have long been used to solve Maxwell's partial differential equations (PDEs) accurately. Finite Difference Time Domain (FDTD), one of the most widely used method, solves Maxwell's PDEs directly in computational grid. In FDTD, grid spacings (Delta x, Delta y, Delta z) are selected to properly sample field quantities to avoid aliasing and maximum allowable time-step (Delta t) is determined to ensure numerical stability of algorithm. Due to discretization of PDEs, FDTD inherentl...
Citation Formats
Ç. Şahin, “Implementation and assessment of modern shock-capturing schemes for hypersonic viscous flows,” M.S. - Master of Science, Middle East Technical University, 2022.