Researcher profile

A. J. Roberts

A. J. Roberts contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
10topics
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

11 published item(s)

preprint2022arXiv

Backwards theory supports modelling via invariant manifolds for non-autonomous dynamical systems

This article establishes the foundation for a new theory of invariant/integral manifolds for non-autonomous dynamical systems. Current rigorous support for dimensional reduction modelling of slow-fast systems is limited by the rare events in stochastic systems that may cause escape, and limited in many applications by the unbounded nature of PDE operators. To circumvent such limitations, we initiate developing a backward theory of invariant/integral manifolds that complements extant forward theory. Here, for deterministic non-autonomous ODE systems, we construct a conjugacy with a normal form system to establish the existence, emergence and exact construction of center manifolds in a finite domain for systems `arbitrarily close' to that specified. A benefit is that the constructed invariant manifolds are known to be exact for systems `close' to the one specified, and hence the only error is in determining how close over the domain of interest for any specific application. Built on the base developed here, planned future research should develop a theory for stochastic and/or PDE systems that is useful in a wide range of modelling applications.

preprint2022arXiv

Embed to rigorously and accurately homogenise quasi-periodic multi-scale heterogeneous PDEs, with computer algebra

For microscale heterogeneous PDEs, this article further develops novel theory and methodology for their macroscale mathematical/asymptotic homogenization. This article specifically encompasses the case of quasi-periodic heterogeneity with finite scale separation: no scale separation limit is required. Dynamical systems theory frames the homogenization as a slow manifold of the ensemble of all phase-shifts of the heterogeneity. Depending upon any perceived scale separation within the quasi-periodic heterogeneity, the homogenization may be done in either one step, or two sequential steps: the results are equivalent. The theory not only assures us of the existence and emergence of the homogenization, it also provides a practical systematic method to construct the homogenization to any specified order. For a class of heterogeneities, we show that the macroscale homogenization is potentially valid down to lengths which are just twice that of the microscale heterogeneity! This methodology provides a new rigorous and flexible approach to homogenization that potentially also provides correct initial and boundary conditions, treatment of forcing and control, and analysis of uncertainty.

preprint2022arXiv

Learning high-order spatial discretisations of PDEs with symmetry-preserving iterative algorithms

Common techniques for the spatial discretisation of PDEs on a macroscale grid include finite difference, finite elements and finite volume methods. Such methods typically impose assumed microscale structures on the subgrid fields, so without further tailored analysis are not suitable for systems with subgrid-scale heterogeneity or nonlinearities. We provide a new algebraic route to systematically approximate, in principle exactly, the macroscale closure of the spatially-discrete dynamics of a general class of heterogeneous non-autonomous reaction-advection-diffusion PDEs. This holistic discretisation approach, developed through rigorous theory and verified with computer algebra, systematically constructs discrete macroscale models through physics informed by the PDE out-of-equilibrium dynamics, thus relaxing many assumptions regarding the subgrid structure. The construction is analogous to recent gray-box machine learning techniques in that predictions are directed by iterative layers (as in neural networks), but informed by the subgrid physics (or 'data') as expressed in the PDEs. A major development of the holistic methodology, presented herein, is novel inter-element coupling between subgrid fields which preserve self-adjointness of the PDE after macroscale discretisation, thereby maintaining the spectral structure of the original system. This holistic methodology also encompasses homogenisation of microscale heterogeneous systems, as shown here with the canonical examples of heterogeneous 1D waves and diffusion.

preprint2022arXiv

Staggered grids for multidimensional multiscale modelling

Numerical schemes for wave-like systems with small dissipation are often inaccurate and unstable due to truncation errors and numerical roundoff errors. Hence, numerical simulations of wave-like systems lacking proper handling of these numerical issues often fail to represent the physical characteristics of wave phenomena. This challenge gets even more intricate for multiscale modelling, especially in multiple dimensions. When using the usual collocated grid, about two-thirds of the resolved wave modes are incorrect with significant dispersion. But, numerical schemes on staggered grids (with alternating variable arrangement) are significantly less dispersive and preserve much of the wave characteristics. Also, the group velocity of the energy propagation in the numerical waves on a staggered grid is in the correct direction, in contrast to the collocated grid. For high accuracy and to preserve much of the wave characteristics, this article extends the concept of staggered grids in full-domain modelling to multidimensional multiscale modelling. Specifically, this article develops 120 multiscale staggered grids and demonstrates their stability, accuracy, and wave-preserving characteristic for equation-free multiscale modelling of weakly damped linear waves. But most characteristics of the developed multiscale staggered grids must also hold in general for multiscale modelling of many complex spatio-temporal physical phenomena such as the general computational fluid dynamics.

preprint2020arXiv

A multiscale scheme accurately simulates macroscale shocks in an equation-free framework

Scientists and engineers often create accurate, trustworthy, computational simulation schemes - but all too often these are too computationally expensive to execute over the time or spatial domain of interest. The equation-free approach is to marry such trusted simulations to a framework for numerical macroscale reduction - the patch dynamics scheme. This article extends the patch scheme to scenarios in which the trusted simulation resolves abrupt state changes on the microscale that appear as shocks on the macroscale. Accurate simulation for problems in these scenarios requires extending the patch scheme by capturing the shock within a novel patch, and also modifying the patch coupling rules in the vicinity in order to maintain accuracy. With these two extensions to the patch scheme, straightforward arguments derive consistency conditions that match the usual order of accuracy for patch schemes. The new scheme is successfully tested on four archetypal problems. This technique will empower scientists and engineers to accurately and efficiently simulate, over large spatial domains, multiscale multiphysics systems that have rapid transition layers on the microscale.

preprint2020arXiv

A toolbox of Equation-Free functions in Matlab\Octave for efficient system level simulation

The `equation-free toolbox' empowers the computer-assisted analysis of complex, multiscale systems. Its aim is to enable you to immediately use microscopic simulators to perform macro-scale system level tasks and analysis, because micro-scale simulations are often the best available description of a system. The methodology bypasses the derivation of macroscopic evolution equations by computing the micro-scale simulator only over short bursts in time on small patches in space, with bursts and patches well-separated in time and space respectively. We introduce the suite of coded equation-free functions in an accessible way, link to more detailed descriptions, discuss their mathematical support, and introduce a novel and efficient algorithm for Projective Integration. Some facets of toolbox development of equation-free functions are then detailed. Download the toolbox functions (https://github.com/uoa1184615/EquationFreeGit) and use to empower efficient and accurate simulation in a wide range of your science and engineering problems.

preprint2020arXiv

Equation-free patch scheme for efficient computational homogenisation via self-adjoint coupling

Equation-free macroscale modelling is a systematic and rigorous computational methodology for efficiently predicting the dynamics of a microscale system at a desired macroscale system level. In this scheme, the given microscale model is computed in small patches spread across the space-time domain, with patch coupling conditions bridging the unsimulated space. For accurate simulations, care must be taken in designing the patch coupling conditions. Here we construct novel coupling conditions which preserve translational invariance, rotational invariance, and self-adjoint symmetry, thus guaranteeing that conservation laws associated with these symmetries are preserved in the macroscale simulation. Spectral and algebraic analyses of the proposed scheme in both one and two dimensions reveal mechanisms for further improving the accuracy of the simulations. Consistency of the patch scheme's macroscale dynamics with the original microscale model is proved. This new self-adjoint patch scheme provides an efficient, flexible, and accurate computational homogenisation in a wide range of multiscale scenarios of interest to scientists and engineers.

preprint2020arXiv

Rigorous modelling of nonlocal interactions determines a macroscale advection-diffusion PDE

A slowly-varying or thin-layer multiscale assumption empowers macroscale understanding of many physical scenarios from dispersion in pipes and rivers, including beams, shells, and the modulation of nonlinear waves, to homogenisation of micro-structures. Here we begin a new exploration of the scenario where the given physics has non-local microscale interactions. We rigorously analyse the dynamics of a basic example of shear dispersion. Near each cross-section, the dynamics is expressed in the local moments of the microscale non-local effects. Centre manifold theory then supports the local modelling of the system's dynamics with coupling to neighbouring cross-sections as a non-autonomous forcing. The union over all cross-sections then provides powerful new support for the existence and emergence of a macroscale model advection-diffusion PDE global in the large, finite-sized, domain. The approach quantifies the accuracy of macroscale advection-diffusion approximations, and has the potential to open previously intractable multiscale issues to new insights.

preprint2010arXiv

Large deviations for slow-fast stochastic partial differential equations

A large deviation principle is derived for stochastic partial differential equations with slow-fast components. The result shows that the rate function is exactly that of the averaged equation plus the fluctuating deviation which is a stochastic partial differential equation with small Gaussian perturbation. This also confirms the effectiveness of the approximation of the averaged equation plus the fluctuating deviation to the slow-fast stochastic partial differential equations.

preprint2010arXiv

Macroscopic discrete modelling of stochastic reaction-diffusion equations on a periodic domain

Dynamical systems theory provides powerful methods to extract effective macroscopic dynamics from complex systems with slow modes and fast modes. Here we derive and theoretically support a macroscopic, spatially discrete, model for a class of stochastic reaction-diffusion partial differential equations with cubic nonlinearity. Dividing space into overlapping finite elements, a special coupling condition between neighbouring elements preserves the self-adjoint dynamics and controls interelement interactions. When the interelement coupling parameter is small, an averaging method and an asymptotic expansion of the slow modes show that the macroscopic discrete model will be a family of coupled stochastic ordinary differential equations which describe the evolution of the grid values. This modelling shows the importance of subgrid scale interaction between noise and spatial diffusion and provides a new rigourous approach to constructing semi-discrete approximations to stochastic reaction-diffusion partial differential equations.

preprint2010arXiv

Resolve subgrid microscale interactions to discretise stochastic partial differential equations

Constructing discrete models of stochastic partial differential equations is very delicate. Stochastic centre manifold theory provides novel support for coarse grained, macroscale, spatial discretisations of nonlinear stochastic partial differential or difference equations such as the example of the stochastically forced Burgers' equation. Dividing the physical domain into finite length overlapping elements empowers the approach to resolve fully coupled dynamical interactions between neighbouring elements. The crucial aspect of this approach is that the underlying theory organises the resolution of the vast multitude of subgrid microscale noise processes interacting via the nonlinear dynamics within and between neighbouring elements. Noise processes with coarse structure across a finite element are the most significant noises for the discrete model. Their influence also diffuses away to weakly correlate the noise in the spatial discretisation. Nonlinear interactions have two further consequences: additive forcing generates multiplicative noise in the discretisation; and effectively new noise processes appear in the macroscale discretisation. The techniques and theory developed here may be applied to soundly discretise many dissipative stochastic partial differential and difference equations.