A diffuse interface model of reactive-fluids and solid-dynamics. (1st October 2021)

• Record Type: Journal Article
• Authors: Wallis, Tim; Barton, Philip T.; Nikiforakis, Nikolaos
• Is Part Of: Computers & structures. Volume 254(2021)
• Publication Date: 2021-10-01
• Subjects: Multi-physics -- Multi-phase -- Reactive fluids -- Elastoplastic solids -- Interface sharpening -- Diffuse interface; Structural engineering -- Data processing -- Periodicals; Electronic data processing -- Structures, Theory of -- Periodicals; 624.171
• DOI: 10.1016/j.compstruc.2021.106578
• Languages: English
• ISSNs: 0045-7949

Abstract

Highlights: A novel method to simulate elastoplastic solids, reactive mixtures and inert fluids. The diffuse interface method allows full physical interaction between all phases. Interface reconstruction maintains sharp interfaces while preserving conservation. The multi-physics approach shows excellent agreement with demanding experiments. Abstract: This article presents a multi-physics methodology for the numerical simulation of physical systems that involve the non-linear interaction of multi-phase reactive fluids and elastoplastic solids, inducing high strain-rates and high deformations. Each state of matter is governed by a single system of non-linear, inhomogeneous partial differential equations, which are solved simultaneously on the same computational grid, and do not require special treatment of immersed boundaries. To this end, the governing equations for solid and reactive multiphase fluid mechanics are written in the same mathematical form and are discretised on a regular Cartesian mesh. All phase and material boundaries are treated as diffuse interfaces. An interface-steepening technique is employed at material boundaries to keep interfaces sharp whilst maintaining the conservation properties of the system. These algorithms are implemented in a highly-parallelised hierarchical adaptive mesh refinement platform, and are verified and validated using numerical and experimental benchmarks. Results indicate very good agreement with experiment and an improvement of numerical performance compared to certain existing Eulerian methods, without loss of conservation.