multicomponent single-phase flow

[leofor...@gmail.com]

Dear all,

Is it possible to simulate multicomponent single-phase flow by the compositional module? We want the flow assumed to be compressible.

Best

YT

[Olav Møyner]

Hello,

Yes, this is very possible. Your options are either to use equation of state under pressure / temperature ranges where the components are single-phase, or define equilibrium constants and properties using the K-value EOS model.

The upcoming release of MRST will have some improvements to the K-value model, which was fairly experimental at the time of the previous release.

Best regards,

Olav

---

Olav Møyner, PhD

SINTEF Digital

Department of Mathematics & Cybernetics

+47 99 02 49 49

---

From: sinte...@googlegroups.com <sinte...@googlegroups.com> on behalf of leofor...@gmail.com <leofor...@gmail.com>
Sent: Sunday, July 1, 2018 3:13:04 PM
To: MRST-users: The Matlab Reservoir Simulation Toolbox User Group
Subject: [MRST Users] multicomponent single-phase flow

 

--
You received this message because you are subscribed to the Google Groups "MRST-users: The Matlab Reservoir Simulation Toolbox User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sintef-mrst...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/sintef-mrst/b9853059-30e4-4ad7-96cc-9dafed96d632%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[李伊腾]

Dear Olav,

Thank you for your answer. The equation of state is preferable. 

Actually, I have done a simple test by modifying one of your example compositionalBoundaryConditionsExample. However, I'm not quite sure if the simulation result is correct or not, because I have some questions about setting the initial condition and boundary condition. Here I attached the code block of initial and boundary condition in the original script so that you can understand my question.

The first thing I'm confused about is why you define phase properties in the above code, since all phase properties can be computed by the equation of state. In addition, phase properties depend on the pressure and composition, and why you define these properties as constant?  It seems the equation of state is not used in this example. If I want to use the equation of state, how can I call it? 

In my test, I consider an example that a pure component (C1) is injected to displace another pure component (C3) in the 1D domain. At the low enough pressure and high enough temperature, the injected component and the mixture of both components are always in the vapor phase. For single-phase flow, there is no relative permeability. In this case, how to define the flowfluid if it is necessary? Initially, I define the flowfluid as follows

flowfluid = initSimpleADIFluid('phases','G','rho', 0.3477, 'mu', 3e-2*centi*poise, 'c', 1/barsa).

This gives me an error, since relative permeability is still considered but I only define one phase. To make it work, the flowfluid is defined as

flowfluid = initSimpleADIFluid('phases','OG','rho', [500, 0.3477], ...

                      'mu', [1, 3e-2]*centi*poise, 'n', [1, 1], 'c', [0, 1]/barsa).

Although the code can be run, the simulation results depend on the exponents of the relative permeability, which does't make sense.

The above questions trouble me for a couple of days. Hope you can give me some idea or help. Look forward to hearing from you as soon as possible, thank you very much.

Best Regards

Yiteng

2018-07-02 15:11 GMT+03:00 Olav Møyner <Olav....@sintef.no>:

Hello,

Yes, this is very possible. Your options are either to use equation of state under pressure / temperature ranges where the components are single-phase, or define equilibrium constants and properties using the K-value EOS model.

The upcoming release of MRST will have some improvements to the K-value model, which was fairly experimental at the time of the previous release.

Best regards,

Olav

---

Olav Møyner, PhD

SINTEF Digital

Department of Mathematics & Cybernetics

+47 99 02 49 49

---

From: sinte...@googlegroups.com <sinte...@googlegroups.com> on behalf of leofor...@gmail.com <leofor...@gmail.com>
Sent: Sunday, July 1, 2018 3:13:04 PM
To: MRST-users: The Matlab Reservoir Simulation Toolbox User Group
Subject: [MRST Users] multicomponent single-phase flow

 

Dear all,

Is it possible to simulate multicomponent single-phase flow by the compositional module? We want the flow assumed to be compressible.

Best

YT

--
You received this message because you are subscribed to the Google Groups "MRST-users: The Matlab Reservoir Simulation Toolbox User Group" group.

To unsubscribe from this group and stop receiving emails from it, send an email to sintef-mrst+unsubscribe@googlegroups.com.

[Olav Møyner]

Dear Yiteng,

I think I understand the issue you are having, and I can hopefully clear this up.

The boundary conditions you have set up seem correct to me. The main reason why you have to specify a two-phase fluid with rel.perms. is that, while the simulation is single-phase at your pressure and temperature conditions, the solver cannot know this: It is still a two-phase solver. In principle, the components could form two phases and for this reason it always evaluates the rel.perms. (in this case, always at the endpoints, resulting in krO = 0 and krG = 1).

The only parts of the two-phase fluid which is used in this case is the surface density (used to define surface rates) and the rel.perms. (which should always evaluate to 1/0).

I hope this clarifies things.

Best regards,

Olav

---

Olav Møyner, PhD

SINTEF Digital

Department of Mathematics & Cybernetics

+47 99 02 49 49

---

From: sinte...@googlegroups.com <sinte...@googlegroups.com> on behalf of 李伊腾 <leofor...@gmail.com>
Sent: Tuesday, July 3, 2018 2:23:48 PM
To: Olav Møyner
Cc: MRST-users: The Matlab Reservoir Simulation Toolbox User Group
Subject: Re: [MRST Users] multicomponent single-phase flow

 

To unsubscribe from this group and stop receiving emails from it, send an email to sintef-mrst...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/sintef-mrst/CAE%3DEQjfVCk5z5zanmgyxUSkUH0tQ%2B9rr6acE6F5SJ4N3mp99JA%40mail.gmail.com.