Saturation index

[Ibrahim HARB]

Hello all,

I'm just a bit confused regarding the output of "MINERAL_SATURATION_INDEX". Are the output values the saturation index log(IAP/Ksp) or the saturation (IAP/Ksp)? 

Sincerely,

[Hammond, Glenn E]

Ibrahim,

 

It should be IAP/KSP:

 

https://bitbucket.org/pflotran/pflotran/src/4d65d2e5d507ed5677be33cfb1e9e10a2006da67/src/pflotran/patch.F90#lines-6577

https://bitbucket.org/pflotran/pflotran/src/4d65d2e5d507ed5677be33cfb1e9e10a2006da67/src/pflotran/reaction_mineral.F90#lines-1299

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> On Behalf Of Ibrahim HARB
Sent: Thursday, August 10, 2023 4:02 AM
To: pflotran-users <pflotra...@googlegroups.com>
Subject: [pflotran-users: 7220] Saturation index

 

Check twice before you click! This email originated from outside PNNL.

 

Hello all,

 

I'm just a bit confused regarding the output of "MINERAL_SATURATION_INDEX". Are the output values the saturation index log(IAP/Ksp) or the saturation (IAP/Ksp)? 

 

Sincerely,

--
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/6ab02cc5-af92-4b8a-8f03-351ca997e575n%40googlegroups.com.

[Peter Alt-Epping]

Hi all,

Is it possible to use anisotropic permeabilities in unstructured grids?

Thanks,

Peter

To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/PH0PR09MB78367FAA7CA89CB6C285A3909A13A%40PH0PR09MB7836.namprd09.prod.outlook.com.

-- 
Peter Alt-Epping, PD PhD
RWI Group, Institute of Geological Sciences
University of Bern
Baltzerstrasse 3
3012 Bern
Switzerland

Tel.: 0041 31 684 8765

[Moise Rousseau]

Hi Peter,

Yes this is possible. You can assign permeability on the whole material with PERM_X, PERM_Y or PERM_Z keywords or setting it on a cell by by cell basis with HDF5 dataset such as:

DATASET PERM_X
  FILENAME perm_dataset_aniso.h5
END
DATASET PERM_Y
  FILENAME perm_dataset_aniso.h5
END
DATASET PERM_Z
  FILENAME perm_dataset_aniso.h5
END

and

MATERIAL_PROPERTY Mat
  ID 1
  POROSITY 0.3
  TORTUOSITY 0.1
  SATURATION_FUNCTION default
  PERMEABILITY
    ANISOTROPIC
    DATASET PERM_
  /
END

I already do that for full tensor simulations in https://cdnsciencepub.com/doi/10.1139/cgj-2021-0462 with supplementary material at https://zenodo.org/record/5197189.

See also the documentation: https://www.pflotran.org/documentation/user_guide/cards/subsurface/material_property_card.html

Moise

[Hammond, Glenn E]

Peter,

 

If your question is about how we project a diagonal anisotropic permeability tensor onto the unstructured face that is not aligned with the coordinate system, we have several algorithms to do so.  Please see

 

https://documentation.pflotran.org/user_guide/cards/subsurface/grid_card.html?highlight=perm_tensor_to_scalar_model

 

The source code for these calculations is located at

 

https://bitbucket.org/pflotran/pflotran/src/4d65d2e5d507ed5677be33cfb1e9e10a2006da67/src/pflotran/material_aux.F90#lines-378

 

Glenn

[peter.al...@unibe.ch]

Thanks Glenn! Yes, this was my question. 

So I define perm_x, perm_y and perm_z in the MATERIAL_PROPERTY card (as described by Moise) and in addition I specify a PERM_TENSOR_TO_SCALAR_MODEL under GRIDS that defines how the permeability is weighted at the cell faces? Is POTENTIAL a suitable weighting scheme for tet meshes?

Thanks,

Peter

---

From: 'Hammond, Glenn E' via pflotran-users <pflotra...@googlegroups.com>
Sent: Tuesday, August 15, 2023 5:38 PM
To: pflotra...@googlegroups.com <pflotra...@googlegroups.com>
Subject: RE: [pflotran-users: 7224] Re: anisotropic permeability in unstructured grids

 

--

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/PH0PR09MB783668B15823D7757DBBB1E59A14A%40PH0PR09MB7836.namprd09.prod.outlook.com.

[Hammond, Glenn E]

Note that PERM_X, PERM_Y, and PERM_Z can be defined by material (i.e., no need for a dataset).

 

Both POTENTIAL and FLOW are suitable for tets; NOT LINEAR.

 

Actually nothing is suitable for tets as their accuracy with the two-point flux approximation is poor (though one could argue that we cannot accurately characterize the subsurface in the first place).

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> On Behalf Of peter.al...@unibe.ch
Sent: Tuesday, August 15, 2023 10:16 AM
To: pflotra...@googlegroups.com
Subject: Re: [pflotran-users: 7226] Re: anisotropic permeability in unstructured grids

 

Check twice before you click! This email originated from outside PNNL.

 

Thanks Glenn! Yes, this was my question. 

To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/ZR2P278MB10511774E743B9F605626888A314A%40ZR2P278MB1051.CHEP278.PROD.OUTLOOK.COM.