Is it possible to use Perple_X to compute H20 and CO2 EOS?

ASPECT
1

Dear Developers,

ASPECT should be able to handle any material within the Steinberger model. I was wondering if it is possible to compute a tabulated EOS fully composed of H2O and CO2, in order to simulate a free-fluid system within magma.

Additionally, I would like to ask whether ASPECT is the appropriate simulator to model the injection of CO2 or H2O into the subsurface, and the subsequent cracking and hydrofracturing that might occur as the fluid advances upward within the simulation box.

Francesco

2

Dear Francesco,

I was wondering if it is possible to compute a tabulated EOS fully composed of H2O and CO2

You can certainly create an table with density, heat capacity etc., yes. This is not a PerpleX messaging board though, so if you want to do this (see below) I’d recommend you ask for help on the PerpleX boards.

in order to simulate a free-fluid system within magma.

You need to check for yourself which governing equations are needed for your intended simulations, and the values of the parameters required. Then check ASPECT’s manual to see whether it solves those equations.

I would like to ask whether ASPECT is the appropriate simulator to model the injection of CO2 or H2O into the subsurface, and the subsequent cracking and hydrofracturing that might occur as the fluid advances upward within the simulation box.

Again, check for yourself which governing equations are needed for such simulations, and whether ASPECT solves those equations.

Best wishes,
Bob

3

Hi Francesco,

As a follow-up to @bobmyhill message above, I will note that that hydrofracturing (fluid pressure coupled to failure) is currently only an experimental feature in ASPECT and not implemented within the main material models.

There is a general framework in place to combine different material models with reactive fluid transport, but I think what you have in mind would require a few years of very technical code development and testing.

Cheers,
John

4

Surely for my project may require years, but an early version is also appreciated! Where can I find the extension?

Francesco

5

@Francyrad:

The very preliminary hydrofracturing implementation and tests (physics is correct, not tested for time-dependent behavior) are as follows:
melt_visco_plastic.cc
melt_visco_plastic_mohr_circle_compaction.cc
melt_visco_plastic_mohr_circle_shear.cc
melt_visco_plastic_mohr_circle_compaction.prm
melt_visco_plastic_mohr_circle_shear.prm

For the reaction part of the problem, I would start by looking at cookbooks and tests that use the reactive fluid transport material model.