Adaptively detect and accurately resolve macro-scale shocks in an efficient Equation-Free multiscale simulation

by   John Maclean, et al.

The Equation-Free approach to efficient multiscale numerical computation marries trusted micro-scale simulations to a framework for numerical macro-scale reduction – the patch dynamics scheme. A recent novel patch scheme empowered the Equation-Free approach to simulate systems containing shocks on the macro-scale. However, the scheme did not predict the formation of shocks accurately, and it could not simulate moving shocks. This article resolves both issues, as a first step in one spatial dimension, by embedding the Equation-Free, shock-resolving patch scheme within a classic framework for adaptive moving meshes. Our canonical micro-scale problems exhibit heterogeneous nonlinear advection and heterogeneous diffusion. We demonstrate many remarkable benefits from the moving patch scheme, including efficient and accurate macro-scale prediction despite the unknown macro-scale closure. Equation-free methods are here extended to simulate moving, forming and merging shocks without a priori knowledge of the existence or closure of the shocks. Whereas adaptive moving mesh equations are typically stiff, typically requiring small time-steps on the macro-scale, the moving macro-scale mesh of patches is typically not stiff given the context of the micro-scale time-steps required for the sub-patch dynamics.


page 4

page 22

page 24


Large-scale simulation of shallow water waves with computation only on small staggered patches

The multiscale patch scheme is built from given small micro-scale simula...

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

Scientists and engineers often create accurate, trustworthy, computation...

Efficient computational homogenisation of 2D beams of heterogeneous elasticity using the patch scheme

Modern 'smart' materials have complex heterogeneous microscale structure...

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

The `equation-free toolbox' empowers the computer-assisted analysis of c...

A toolbox of Equation-Free functions in Matlabfor efficient system level simulation

The `equation-free toolbox' empowers the computer-assisted analysis of c...

Accurate and efficient multiscale simulation of a heterogeneous elastic beam via computation on small sparse patches

Modern `smart' materials have complex microscale structure, often with u...

A Micro-Macro Parareal Implementation for the Ocean-Circulation Model FESOM2

A micro-macro variant of the parallel-in-time algorithm Parareal has bee...

Please sign up or login with your details

Forgot password? Click here to reset