Source author record

Rob Watson

Rob Watson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis Mathematical Software physics.comp-ph

Catalog footprint

What is connected

3works

4topics

3close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

An extended range of stable flux reconstruction schemes on quadrilaterals for various polynomial bases

An extended range of energy stable flux reconstruction schemes, developed using a summation-by-parts approach, is presented on quadrilateral elements for various sets of polynomial bases. For the maximal order bases, a new set of correction functions which result in stable schemes is found. However, for a range of orders it is shown that only a single correction function can be cast as a tensor-product. Subsequently, correction functions are identified using a generalised analytic framework that results in stable schemes for total order and approximate Euclidean order polynomial bases on quadrilaterals -- which have not previously been explored in the context of flux reconstruction. It is shown that the approximate Euclidean order basis can provide similar numerical accuracy as the maximal order basis but with fewer points per element, and thus lower cost.

preprint2022arXiv

On the Use of RBF Interpolation for Flux Reconstruction

Flux reconstruction provides a framework for solving partial differential equations in which functions are discontinuously approximated within elements. Typically, this is done by using polynomials. Here, the use of radial basis functions as a methods for underlying functional approximation is explored in one dimension, using both analytical and numerical methods. At some mesh densities, RBF flux reconstruction is found to outperform polynomial flux reconstruction, and this range of mesh densities becomes finer as the width of the RBF interpolator is increased. A method which avoids the poor conditioning of flat RBFs is used to test a wide range of basis shapes, and at very small values, the polynomial behaviour is recovered. Changing the location of the solution points is found to have an effect similar to that in polynomial FR, with the Gauss--Legendre points being the most effective. Altering the location of the functional centres is found to have only a very small effect on performance. Similar behaviours are determined for the non-linear Burgers' equation.

preprint2021arXiv

Hyperbolic Diffusion in Flux Reconstruction: Optimisation through Kernel Fusion within Tensor-Product Elements

Novel methods are presented in this initial study for the fusion of GPU kernels in the artificial compressibility method (ACM), using tensor product elements with constant Jacobians and flux reconstruction. This is made possible through the hyperbolisation of the diffusion terms, which eliminates the expensive algorithmic steps needed to form the viscous stresses. Two fusion approaches are presented, which offer differing levels of parallelism. This is found to be necessary for the change in workload as the order of accuracy of the elements is increased. Several further optimisations of these approaches are demonstrated, including a generation time memory manager which maximises resource usage. The fused kernels are able to achieve 3-4 times speedup, which compares favourably with a theoretical maximum speedup of 4. In three dimensional test cases, the generated fused kernels are found to reduce total runtime by ${\sim}25\%$, and, when compared to the standard ACM formulation, simulations demonstrate that a speedup of $2.3$ times can be achieved.