Kadioglu, Samet Y | Search NASA

DOE OSTI2011

A Gas Dynamics Method Based on The Spectral Deferred Corrections (SDC) Time Integration Technique and The Piecewise Parabolic Method (PPM)

We present a computational gas dynamics method based on the Spectral Deferred Corrections (SDC) time integration technique and the Piecewise Parabolic Method (PPM) finite volume method. The PPM framework is used to define edge averaged quantities which are then used to evaluate numerical flux functions. The SDC technique is used to integrate solution in time. This kind of approach was first taken by Anita et al in [17]. However, [17] is problematic when it is implemented to certain shock problems. Here we propose significant improvements to [17]. The method is fourth order (both in space and time) for smooth flows, and provides highly resolved discontinuous solutions. We tested the method by solving variety of problems. Results indicate that the fourth order of accuracy in both space and time has been achieved when the flow is smooth. Results also demonstrate the shock capturing ability of the method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2010

A Second Order Self-Consistent IMEX Method for Radiation Hydrodynamics

We present a second order non-linearly consistent implicit/explicit (methods that use com- bination of implicit/explicit discretizations are often referred to as IMEX (IMplicit/EXplicit) methods [2,1,3]) time integration technique for solving radiation hydrodynamics problems. The idea is to hybridize an implicit and an explicit discretization in a nonlinearly consis- tent way in order to achieve second order time convergent calculations. In this scope, the hydrodynamics equations are discretized explicitly making use of the capability of well- understood explicit schemes. On the other hand , the radiation diffusion part is solved implicitly. In the view of the non-linearly consistent IMEX method, the hydrodynamics equations are solved as part of the non-linear function evaluation within the framework of the Jacobian-Free Newton Krylov (JFNK) method [16,14]. This is to avoid order reductions in time convergence due to the operator splitting. We present results from several test prob- lems in order to validate the numerical order of the scheme. For each test, we established second order time convergence.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗