Samet Y. Kadioglu

A fourth-order auxiliary variable projection method for zero-Mach number gas dynamics

A fourth-order numerical method for the zero-Mach-number limit of the equations for compressible flow is presented. The method is formed by discretizing a new auxiliary variable formulation of the conservation equations, which is a variable density analog to the impulse or gauge formulation of the incompressible Euler equations. An auxiliary variable projection method is applied to this formulation, and accuracy is achieved by combining a fourth-order finite-volume spatial discretization with a fourth-order temporal scheme based on spectral deferred corrections. Numerical results are included which demonstrate fourth-order spatial and temporal accuracy for non-trivial flows in simple geometries.

Adaptive solution techniques for simulating underwater explosions and implosions

Adaptive solution techniques are presented for simulating underwater explosions and implosions. The liquid is assumed to be an adiabatic fluid and the solution in the gas is assumed to be uniform in space. The solution in water is integrated in time using a semi-implicit time discretization of the adiabatic Euler equations. Results are presented either using a non-conservative semi-implicit algorithm or a conservative semi-implicit algorithm. A semi-implicit algorithm allows one to compute with relatively large time steps compared to an explicit method. The interface solver is based on the coupled level set and volume-of-fluid method (CLSVOF) M. Sussman, A second order coupled level set and volume-of-fluid method for computing growth and collapse of vapor bubbles, J. Comput. Phys. 187 (2003) 110-136; M. Sussman, E.G. Puckett, A coupled level set and volume-of-fluid method for computing 3D and axisymmetric incompressible two-phase flows, J. Comput. Phys. 162 (2000) 301-337]. Several underwater explosion and implosion test cases are presented to show the performances of our proposed techniques.

A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems

We present a fully second order implicit/explicit time integration technique for solving hydrodynamics coupled with nonlinear heat conduction problems. The idea is to hybridize an implicit and an explicit discretization in such a way 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 nonlinear heat conduction is solved implicitly. Such methods are often referred to as IMEX methods [2,1,3]. The Jacobian-Free Newton Krylov (JFNK) method (e.g. [10,9]) is applied to the problem in such a way as to render a nonlinearly iterated IMEX method. We solve three test problems in order to validate the numerical order of the scheme. For each test, we established second order time convergence. We support these numerical results with a modified equation analysis (MEA) [21,20]. The set of equations studied here constitute a base model for radiation hydrodynamics.

A second order self-consistent IMEX method for radiation hydrodynamics

We present a second order self-consistent implicit/explicit (methods that use the combination of implicit and explicit discretizations are often referred to as IMEX (implicit/explicit) methods [2,1,3]) time integration technique for solving radiation hydrodynamics problems. The operators of the radiation hydrodynamics are splitted as such that the hydrodynamics equations are solved explicitly making use of the capability of well-understood explicit schemes. On the other hand, the radiation diffusion part is solved implicitly. The idea of the self-consistent IMEX method is to hybridize the implicit and explicit time discretizations in a nonlinearly consistent way to achieve second order time convergent calculations. In our self-consistent IMEX method, we solve the hydrodynamics equations inside the implicit block as part of the nonlinear function evaluation making use of the Jacobian-free Newton Krylov (JFNK) method [5,20,17]. This is done to avoid order reductions in time convergence due to the operator splitting. We present results from several test calculations in order to validate the numerical order of our scheme. For each test, we have established second order time convergence.

Marker Redistancing/Level Set Method for High-Fidelity Implicit Interface Tracking

A hybrid of the front tracking (FT) and the level set (LS) methods is introduced, combining advantages and removing drawbacks of both methods. The kinematics of the interface is treated in a Lagrangian (FT) manner, by tracking markers placed at the interface. The markers are not connected—instead, the interface topology is resolved in an Eulerian (LS) framework, by wrapping a signed distance function around Lagrangian markers each time the markers move. For accuracy and efficiency, we have developed a high-order “anchoring” algorithm and an implicit PDE-based redistancing. We have demonstrated that the method is 3rd-order accurate in space, near the markers, and therefore 1st-order convergent in curvature; this is in contrast to traditional PDE-based reinitialization algorithms, which tend to slightly relocate the zero level set and can be shown to be nonconvergent in curvature. The implicit pseudo-time discretization of the redistancing equation is implemented within the Jacobian-free Newton-Krylov (JFNK) framework combined with ILU(k) preconditioning. Due to the LS localization, the bandwidth of the Jacobian matrix is nearly constant, and the ILU preconditioning scales as

A Second Order JFNK-Based IMEX Method for Single and Multi-Phase Flows

\sim N\log(\sqrt{N})

A Comparative Study of the Harmonic and Arithmetic Averaging of Diffusion Coefficients for Non-linear Heat Conduction Problems

in two dimensions, which implies efficiency and good scalability of the overall algorithm. We have demonstrated that the steady-state solutions in pseudo-time can be achieved very efficiently, with

Multiphysics Analysis of Spherical Fast Burst Reactors

\approx10

A Novel Hyperbolization Procedure for The Two-Phase Six-Equation Flow Model

iterations (

An IMEX Method for the Euler Equations that Posses Strong Non-Linear Heat Conduction and Stiff Source Terms (Radiation Hydrodynamics)

\mathrm{CFL}\approx10^4