BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Will Trojak

Will Trojak contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phphysics.flu-dynComputational Engineering, Finance, and ScienceMathematical Software

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2023arXiv

Bounds Preserving Temporal Integration Methods for Hyperbolic Conservation Laws

In this work, we present a modification of explicit Runge-Kutta temporal integration schemes that guarantees the preservation of any locally-defined quasiconvex set of bounds for the solution. These schemes operate on the basis of a bijective mapping between an admissible set of solutions and the real domain to strictly enforce bounds. Within this framework, we show that it is possible to recover a wide range of methods independently of the spatial discretization, including positivity preserving, discrete maximum principle satisfying, entropy dissipative, and invariant domain preserving schemes. Furthermore, these schemes are proven to recover the order of accuracy of the underlying Runge-Kutta method upon which they are built. The additional computational cost is the evaluation of two nonlinear mappings which generally have closed-form solutions. We show the utility of this approach in numerical experiments using a pseudospectral spatial discretization without any explicit shock capturing schemes for nonlinear hyperbolic problems with discontinuities.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

On the anti-aliasing properties of entropy filtering for discontinuous spectral element approximations of under-resolved turbulent flows

For large Reynolds number flows, it is typically necessary to perform simulations that are under-resolved with respect to the underlying flow physics. For nodal discontinuous spectral element approximations of these under-resolved flows, the collocation projection of the nonlinear flux can introduce aliasing errors which can result in numerical instabilities. In Dzanic and Witherden (J. Comput. Phys., 468, 2022), an entropy-based adaptive filtering approach was introduced as a robust, parameter-free shock-capturing method for discontinuous spectral element methods. This work explores the ability of entropy filtering for mitigating aliasing-driven instabilities in the simulation of under-resolved turbulent flows through high-order implicit large eddy simulations of a NACA0021 airfoil in deep stall at a Reynolds number of 270,000. It was observed that entropy filtering can adequately mitigate aliasing-driven instabilities without degrading the accuracy of the underlying high-order scheme on par with standard anti-aliasing methods such as over-integration, albeit with marginally worse performance at higher approximation orders.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An extended range of energy stable flux reconstruction methods on triangles

We present an extended range of stable flux reconstruction (FR) methods on triangles through the development and application of the summation-by-parts framework in two-dimensions. This extended range of stable schemes is then shown to contain the single parameter schemes of \citet{Castonguay2011} on triangles, and our definition enables wider stability bounds to be developed for those single parameter families. Stable upwinded spectral difference (SD) schemes on triangular elements have previously been found using Fourier analysis. We used our extended range of FR schemes to investigate the linear stability of SD methods on triangles, and it was found that a only first order SD scheme could be recovered within this set of FR methods.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Utilizing Time-Reversibility for Shock Capturing in Nonlinear Hyperbolic Conservation Laws

In this work, we introduce a novel approach to formulating an artificial viscosity for shock capturing in nonlinear hyperbolic systems by utilizing the property that the solutions of hyperbolic conservation laws are not reversible in time in the vicinity of shocks. The proposed approach does not require any additional governing equations or a priori knowledge of the hyperbolic system in question, is independent of the mesh and approximation order, and requires the use of only one tunable parameter. The primary novelty is that the resulting artificial viscosity is unique for each component of the conservation law which is advantageous for systems in which some components exhibit discontinuities while others do not. The efficacy of the method is shown in numerical experiments of multi-dimensional hyperbolic conservation laws such as nonlinear transport, Euler equations, and ideal magnetohydrodynamics using a high-order discontinuous spectral element method on unstructured grids.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Family of Weighted One-Parameter Flux Reconstruction Schemes

The flux reconstruction (FR) approach offers a flexible framework for describing a range of high-order numerical schemes; including nodal discontinuous Galerkin and spectral difference schemes. This is accomplished through the use of so-called correction functions. In this study we employ a weighted Sobolev norm to define a new extended family of FR correction functions, the stability of which is affirmed through Fourier analysis. Several of the schemes within this family are found to exhibit reduced dissipation and dispersion overshoot. Moreover, many of the new schemes possess higher CFL limits whilst maintaining the expected rate of convergence. Numerical experiments with homogeneous linear convection and Burgers turbulence are undertaken, and the results observed to be in agreement with the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inline Vector Compression for Computational Physics

A novel inline data compression method is presented for single-precision vectors in three dimensions. The primary application of the method is for accelerating computational physics calculations where the throughput is bound by memory bandwidth. The scheme employs spherical polar coordinates, angle quantisation, and a bespoke floating-point representation of the magnitude to achieve a fixed compression ratio of 1.5. The anisotropy of this method is considered, along with companding and fractional splitting techniques to improve the efficiency of the representation. We evaluate the scheme numerically within the context of high-order computational fluid dynamics. For both the isentropic convecting vortex and the Taylor--Green vortex test cases, the results are found to be comparable to those without compression. Performance is evaluated for a vector addition kernel on an NVIDIA Titan V GPU; it is demonstrated that a speedup of 1.5 can be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On Fourier analysis of polynomial multigrid for arbitrary multi-stage cycles

The Fourier analysis of the \emph{p}-multigrid acceleration technique is considered for a dual-time scheme applied to the advection-diffusion equation with various cycle configurations. It is found that improved convergence can be achieved through \emph{V}-cycle asymmetry where additional prolongation smoothing is applied. Experiments conducted on the artificial compressibility formulation of the Navier--Stokes equations found that these analytic findings could be observed numerically in the pressure residual, whereas velocity terms---which are more hyperbolic in character---benefited primarily from increased pseudo-time steps.

0 citations0 reviewsNo code linkedOpen access