High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used

Ising exponents from the functional renormalisation group

We study the 3d Ising universality class using the functional renormalisation group. With the help of background fields and a derivative expansion up to fourth order we compute the leading index, the subleading symmetric and anti-symmetric corrections to scaling, the anomalous dimension, the scaling solution, and the eigenperturbations at criticality. We also study the cross-correlations of scaling exponents, and their dependence on dimensionality. We find a very good numerical convergence of the derivative expansion, also in comparison with earlier findings. Evaluating the data from all functional renormalisation group studies to date, we estimate the systematic error which is found to be small and in good agreement with findings from Monte Carlo simulations, \epsilon-expansion techniques, and resummed perturbation theory.Comment: 24 pages, 3 figures, 7 table