Efficient Solution of Bimaterial Riemann Problems for Compressible Multi-Material Flow Simulations

by   Wentao Ma, et al.

When solving compressible multi-material flow problems, an unresolved challenge is the computation of advective fluxes across material interfaces that separate drastically different thermodynamic states and relations. A popular idea in this regard is to locally construct bimaterial Riemann problems, and to apply their exact solutions in flux computation. For general equations of state, however, finding the exact solution of a Riemann problem is expensive as it requires nested loops. Multiplied by the large number of Riemann problems constructed during a simulation, the computational cost often becomes prohibitive. The work presented in this paper aims to accelerate the solution of bimaterial Riemann problems without introducing approximations or offline precomputation tasks. The basic idea is to exploit some special properties of the Riemann problem equations, and to recycle previous solutions as much as possible. Following this idea, four acceleration methods are developed, including (1) a change of integration variable through rarefaction fans, (2) storing and reusing integration trajectory data, (3) step size adaptation, and (4) constructing an R-tree on the fly to generate initial guesses. The performance of these acceleration methods are assessed using four example problems in underwater explosion, laser-induced cavitation, and hypervelocity impact. These problems exhibit strong shock waves, large interface deformation, contact of multiple (>2) interfaces, and interaction between gases and condensed matters. In these challenging cases, the solution of bimaterial Riemann problems is accelerated by 37 to 83 times. As a result, the total cost of advective flux computation, which includes the exact Riemann problem solution at material interfaces and the numerical flux calculation over the entire computational domain, is accelerated by 18 to 79 times.


page 30

page 35

page 36

page 38


XIGA: An eXtended IsoGeometric Analysis approach for multi-material problems

Multi-material problems often exhibit complex geometries along with phys...

An energy-splitting high order numerical method for multi-material flows

This chapter deals with multi-material flow problems by a kind of effect...

Numerical simulation of hydraulic fracturing: a hybrid FEM based algorithm

In this paper a problem of numerical simulation of hydraulic fractures i...

Biomechanical modeling and computer simulation of the brain during neurosurgery

Computational biomechanics of the brain for neurosurgery is an emerging ...

A FETI approach to domain decomposition for meshfree discretizations of nonlocal problems

We propose a domain decomposition method for the efficient simulation of...

A Reduced Order Model for Joint Assemblies by Hyper-Reduction and Model-Driven Sampling

The dynamic behavior of jointed assemblies exhibiting friction nonlinear...

Momentum-based Accelerated Mirror Descent Stochastic Approximation for Robust Topology Optimization under Stochastic Loads

Robust topology optimization (RTO) improves the robustness of designs wi...

Please sign up or login with your details

Forgot password? Click here to reset