Example of salinity-based changes in density

[Sarah McKnight]

Hello PFLOTRAN community,

I am reaching out to see if anyone has an example simulation of vertical density-driven flow in PFLOTRAN that is fully saturated and is coupled with salinity. I have been looking into the "1d_soluble_matrix" example, which does have a fully coupled formulation, but it is not fully saturated. I would like to see if I could start from an example that is both fully saturated and coupling density with salinity. Any help would be greatly appreciated.

Thank you,

Sarah

[benjami...@gmail.com]

Sarah,

Here's a not-particularly-realistic, but nevertheless representative, example of a fully saturated, salinity-coupled simulation. I say not-particularly-realistic because I just threw in some liquid pressures as the boundary condition instead of dealing with calculating and adding hydrostatic boundary condition files. Nevertheless, you can see that the salinity is affecting density:

[A Palacios]

Hi Sarah, 

This is another example that is used by developers to train on density driven flow. https://bitbucket.org/pflotran/pflotran/src/master/shortcourse/exercises/density_dependent_flow/

In the users group there are also several conversations about it, and you could find Henry benchmark for example. 

Good luck,
Andrea

[Fukuyama, David Erik]

Hi Sara,

I wrote the implicit salt capability for general mode. If you want to simulate salinity-based changes in density in a non-soluble rock (sandstone, for example), you should remove the SOLUBLE_MATRIX keyword in the general mode options block. The SOLUBLE_MATRIX keyword is for rocks that can dissolve in water, such as halite. For fully saturated flow in a sandstone, you will need to specify liquid pressure, air mole fraction, salt mole fraction, and temperature.  

For example:

FLOW_CONDITION initial

  TYPE

    LIQUID_PRESSURE DIRICHLET

    MOLE_FRACTION DIRICHLET

    SALT_MOLE_FRACTION DIRICHLET

    TEMPERATURE DIRICHLET

  /

  LIQUID_PRESSURE 111325.d0

  MOLE_FRACTION 1.d-10

  SALT_MOLE_FRACTION 1.d-4

  TEMPERATURE 25.d0

/

 

A good starting point is the 1d_solute_mixedbc.in regression test. When the SOLUTE SALT is specified in the GENERAL options block, density is coupled to salinity/salt mole fraction by default with the Batzle and Wang EOS.

 

David

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> on behalf of A Palacios <andr...@gmail.com>
Date: Thursday, October 5, 2023 at 6:43 AM
To: pflotran-users <pflotra...@googlegroups.com>
Subject: [EXTERNAL] [pflotran-users: 7317] Re: Example of salinity-based changes in density

Hi Sarah, 

This is another example that is used by developers to train on density driven flow. https://bitbucket.org/pflotran/pflotran/src/master/shortcourse/exercises/density_dependent_flow/

In the users group there are also several conversations about it, and you could find Henry benchmark for example. 

Good luck,
Andrea

On Thursday, October 5, 2023 at 12:34:26 PM UTC+2 benjami...@gmail.com wrote:

Sarah,

Here's a not-particularly-realistic, but nevertheless representative, example of a fully saturated, salinity-coupled simulation. I say not-particularly-realistic because I just threw in some liquid pressures as the boundary condition instead of dealing with calculating and adding hydrostatic boundary condition files. Nevertheless, you can see that the salinity is affecting density:

 

On Wednesday, 4 October 2023 at 16:11:39 UTC+1 sarah.v....@gmail.com wrote:

Hello PFLOTRAN community,

 

I am reaching out to see if anyone has an example simulation of vertical density-driven flow in PFLOTRAN that is fully saturated and is coupled with salinity. I have been looking into the "1d_soluble_matrix" example, which does have a fully coupled formulation, but it is not fully saturated. I would like to see if I could start from an example that is both fully saturated and coupling density with salinity. Any help would be greatly appreciated.

 

Thank you,

Sarah

--
You received this message because you are subscribed to the Google Groups "pflotran-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-user...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/64818e31-d0df-4d11-9c86-764887ad4cf6n%40googlegroups.com.

[Sarah McKnight]

Thank you so much! This is all very helpful. I will follow up if I have any questions.