Adaptively detect and accurately resolve macro-scale shocks in an efficient equation-free multiscale simulation

DOI10.1137/21M1437172MaRDI QIDQ5097608zbMATH OpenOpenAlexFDO

Authors John MacLean, J. E. Bunder, I. G. Kevrekidis, A. J. Roberts

Publication date 25 August 2022

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2108.11568

zbMATH Keywords

equation-freematerials modelingmultiscale simulationshock-resolving

Mathematics Subject Classification ID

Homogenization in context of PDEs; PDEs in media with periodic structure (35B27) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12)

Abstract: 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.

Recommendations

• Equation-free computation: An overview of patch dynamics
• Equation-free, coarse-grained multiscale computation: enabling microscopic simulators to perform system-level analysis
• Finite difference patch dynamics for advection homogenization problems
• Equation-Free, Multiscale Computation for Unsteady Random Diffusion
• Patch dynamics for macroscale modelling in one dimension

Cites work

• A dynamical systems approach to simulating macroscale spatial dynamics in multiple dimensions
• A holistic finite difference approach models linear dynamics consistently
• A note on implementations of the boosting algorithm and heterogeneous multiscale methods
• A toolbox of equation-free functions in Matlab/Octave for efficient system level simulation
• Adaptivity with moving grids
• Advanced methods of continuum mechanics for materials and structures
• An analysis of stability and convergence of a finite-difference discretization of a model parabolic PDE in 1D using a moving mesh
• Couple microscale periodic patches to simulate macroscale emergent dynamics
• DeepXDE: a deep learning library for solving differential equations
• Equation-free computation: An overview of patch dynamics
• Equation-free patch scheme for efficient computational homogenisation via self-adjoint coupling
• Equation-free, coarse-grained multiscale computation: enabling microscopic simulators to perform system-level analysis
• General Tooth Boundary Conditions for Equation Free Modeling
• Good coupling for the multiscale patch scheme on systems with microscale heterogeneity
• Heterogeneous multiscale methods: a review
• Model emergent dynamics in complex systems
• Model reduction for transport-dominated problems via online adaptive bases and adaptive sampling
• Multiscale modelling couples patches of non-linear wave-like simulations
• Numerical homogenization of non-linear parabolic problems on adaptive meshes
• Numerical modeling of concrete cracking. Papers based on the presentations at the advanced school, Udine, Italy, May 2009.
• On convergence of higher order schemes for the projective integration method for stiff ordinary differential equations
• On convergence of the projective integration method for stiff ordinary differential equations
• Patch dynamics with buffers for homogenization problems
• Resolution of subgrid microscale interactions enhances the discretisation of nonautonomous partial differential equations
• Stability of Moving Mesh Systems of Partial Differential Equations
• The Exact Order of Convergence for Finite Difference Approximations to Ordinary Boundary Value Problems
• The gap-tooth method in particle simulations
• The Gap-Tooth Scheme for Homogenization Problems
• The heterogeneous multiscale method

Cited in

(5)

• Equation free projective integration: a multiscale method applied to a plasma ion acoustic wave
• Two novel families of multiscale staggered patch schemes efficiently simulate large-scale, weakly damped, linear waves
• Numerical homogenization by an adaptive Fourier spectral method on non-uniform grids using optimal transport
• Generalised patch dynamics schemes in equation-free multiscale modelling
• Efficient prediction of static and dynamical responses of functional graded beams using sparse multiscale patches

This page was built for publication: Adaptively detect and accurately resolve macro-scale shocks in an efficient equation-free multiscale simulation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5097608)