ePrints@IIScePrints@IISc Home | About | Browse | Latest Additions | Advanced Search | Contact | Help

Evaluation of finite difference based asynchronous partial differential equations solver for reacting flows

Kumari, K and Cleary, E and Desai, S and Donzis, DA and Chen, JH and Aditya, K (2023) Evaluation of finite difference based asynchronous partial differential equations solver for reacting flows. In: Journal of Computational Physics, 477 .

jou_com_phy_477_2023.pdf - Published Version

Download (1MB) | Preview
Official URL: https://doi.org/10.1016/j.jcp.2023.111906


Next-generation exascale machines with extreme levels of parallelism will provide massive computing resources for large scale numerical simulations of complex physical systems at unprecedented parameter ranges. However, novel numerical methods, scalable algorithms and re-design of current state-of-the art numerical solvers are required for scaling to these machines with minimal overheads. One such approach for partial differential equations based solvers involves computation of spatial derivatives with possibly delayed or asynchronous data using high-order asynchrony-tolerant (AT) schemes to facilitate mitigation of communication and synchronization bottlenecks without affecting the numerical accuracy. In the present study, an effective methodology of implementing temporal discretization using a multi-stage Runge-Kutta method with AT schemes is presented. Together these schemes are used to perform asynchronous simulations of canonical reacting flow problems, demonstrated in one-dimension including auto-ignition of a premixture, premixed flame propagation and non-premixed autoignition. Simulation results show that the AT schemes incur very small numerical errors in all key quantities of interest including stiff intermediate species despite delayed data at processing element (PE) boundaries. For simulations of supersonic flows, the degraded numerical accuracy of well-known shock-resolving WENO (weighted essentially non-oscillatory) schemes when used with relaxed synchronization is also discussed. To overcome this loss of accuracy, high-order AT-WENO schemes are derived and tested on linear and non-linear equations. Finally the novel AT-WENO schemes are demonstrated in the propagation of a detonation wave with delays at PE boundaries. © 2023 Elsevier Inc.

Item Type: Journal Article
Publication: Journal of Computational Physics
Publisher: Academic Press Inc.
Additional Information: The copyright for this article belongs to the Authors.
Keywords: Data handling; Detonation; Numerical methods; Partial differential equations, Asynchronous computing; Asynchrony; DNS; Essentially non-oscillatory; Essentially non-oscillatory schemes; Exascale; Finite-difference; High-order; Reacting flows; Weighted essentially non-oscillatory, Runge Kutta methods
Department/Centre: Division of Interdisciplinary Sciences > Computational and Data Sciences
Date Deposited: 16 Feb 2023 03:52
Last Modified: 16 Feb 2023 03:52
URI: https://eprints.iisc.ac.in/id/eprint/80285

Actions (login required)

View Item View Item