USE OF GEOMETRICAL OPTIC METHODS FOR DISADVANTAGES OF FDTD METHOD

2014-03-01
Ciydem, Mehmet
Koç, Seyit Sencer
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 inherently suffers from numerical dispersion, which results in numerical velocity errors and anisotropy in the grid. Anisotropy and different velocities result in numerical phase errors in the solution and it accumulates within the grid. Moreover, some modes in the grid propagate faster than light. In this study, contrary to FDTD, Geometrical Optic methods have been utilized and a new computational method called as Ray-Based Finite Difference (RBFD) method has been proposed for computational electromagnetics. Discontinuities in the fields and their successive time derivatives can only exist on the wavefronts and propagate along the rays. They are transported in computational domain by transport equations that are ordinary differential equations (ODEs). In isotropic media, energy flows in ray direction, which is perpendicular to the wavefronts. RBFD mainly utilizes directional energy flow property for grid generation and ODE nature of transport equations for numerical computations. Simulation results show that RBFD can be exploited to eliminate disadvantges of FDTD.
JOURNAL OF THE FACULTY OF ENGINEERING AND ARCHITECTURE OF GAZI UNIVERSITY

Suggestions

Implementation of different flux evaluation schemes into a two-dimensional Euler solver
Eraslan, Elvan; Aksel, Mehmet Haluk; Department of Mechanical Engineering (2006)
This study investigates the accuracy and efficiency of several flux splitting methods for the compressible, two-dimensional Euler equations. Steger-Warming flux vector splitting method, Van Leer flux vector splitting method, The Advection Upstream Splitting Method (AUSM), Artificially Upstream Flux Vector Splitting Scheme (AUFS) and Roe’s flux difference splitting schemes were implemented using the first- and second-order reconstruction methods. Limiter functions were embedded to the second-order reconstruc...
Development of an incompressible navier-stokes solver with alternating cell direction implicit method on structured and unstructured quadrilateral grids
Baş, Onur; Tuncer, İsmail Hakkı; Department of Aerospace Engineering (2007)
In this research, the Alternating Cell Direction Implicit method is used in temporal discretisation of the incompressible Navier-Stokes equations and compared with the well known and widely used Point Gauss Seidel scheme on structured and quadrilateral unstructured meshes. A two dimensional, laminar and incompressible Navier-Stokes solver is developed for this purpose using the artificial compressibility formulation. The developed solver is used to obtain steady-state solutions with implicit time stepping m...
Two dimensional finite volume weighted essentially non-oscillatory euler schemes with uniform and non-uniform grid coefficients
Elfarra, Monier Ali; Akmandor, İbrahim Sinan; Department of Aerospace Engineering (2005)
In this thesis, Finite Volume Weighted Essentially Non-Oscillatory (FV-WENO) codes for one and two-dimensional discretised Euler equations are developed. The construction and application of the FV-WENO scheme and codes will be described. Also the effects of the grid coefficients as well as the effect of the Gaussian Quadrature on the solution have been tested and discussed. WENO schemes are high order accurate schemes designed for problems with piecewise smooth solutions containing discontinuities. The key ...
Aerodynamic design optimization of three dimensional rocket nozzles using adjoint method
Eyi, Sinan (null; 2013-09-13)
A design optimization method based on three dimensional Euler equations is developed. A finite volume method is implemented to discretize the Euler equations. Newton's method is used to solve the discretized form of Euler equations. Newton's method requires the calculation of the Jacobian matrix which is the derivative of the residual vector with respect to the flux vector. Different upwind methods are used in the calculation of flux vectors. Numerical and analytical methods are utilized in the evaluation o...
Implementation and assessment of modern shock-capturing schemes for hypersonic viscous flows
Şahin, Çağatay; Baran, Özgür Uğraş; Atik, Hediye; Department of Mechanical Engineering (2022-8-22)
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...
Citation Formats
M. Ciydem and S. S. Koç, “USE OF GEOMETRICAL OPTIC METHODS FOR DISADVANTAGES OF FDTD METHOD,” JOURNAL OF THE FACULTY OF ENGINEERING AND ARCHITECTURE OF GAZI UNIVERSITY, pp. 121–129, 2014, Accessed: 00, 2020. [Online]. Available: https://hdl.handle.net/11511/52528.