Speaker
Description
Buoyant incompressible multiphase flow and transport in porous media typically is simulated by coupling an elliptic flow equation with a hyperbolic transport description. The tight coupling between flow and transport either requires a fully implicit solution algorithm or very small time steps, if solved sequentially. In any case, however, a large linear system has to be solved each time step.
Here a new solution approach is devised, which relies on solving a coupled hyperbolic system of conservation laws with an explicit finite volume method. Consequently, only local operations have to be performed, which is a great advantage in case of massive parallel simulations and if GPUs are employed.
The devised method is based on the isothermal Euler equations with momentum source terms accounting for resistance due to the porous medium and buoyancy. In this system the pressure is proportional to the density and in case of very small Mach numbers and relative density changes, the computed velocity field is divergence free and converges towards the total volumetric flow in the porous domain. To account for saturation transport, the system was augmented by an additional hyperbolic equation. In order to obtain the numerical fluxes, a characteristic based approximative Riemann solver was developed and 2nd order accuracy in space and time is achieved by piecewise linear reconstruction.
Two-phase flow test cases with buoyant plumes and a lock exchange problem demonstrate that the devised hyperbolic approach recovers the
the correct incompressible solutions. Numerical experiments also demonstrate that the method is accurate and efficient.
| Participation | In person |
|---|---|
| Country | Switzerland |
| MDPI Energies Student Poster Award | No, do not submit my presenation for the student posters award. |
| Time Block Preference | Time Block A (09:00-12:00 CET) |
| Acceptance of the Terms & Conditions | Click here to agree |
