Anisotropic Permeability Field for Hydrology Runs

[Saubhagya Singh Rathore]

Hello ATS community, 

Is it possible to provide a subsurface layer with anisotropic permeability field for an integrated hydrology run? Is there an example available on how this can be achieved?

Thanks

Saubhagya

[Coon, Ethan]

Yes, there is an example of tensor permeability in ats-regression-tests,

https://github.com/amanzi/ats-regression-tests/tree/master/02_richards/seepage_drawdown_full_tensor_mfd.regression.gold

 

We haven’t really tried this too much though.  The discretizations are very ready for this as Konstantin has tested this quite well in Amanzi, but the conceptualization of how to do this is tricky.  For instance, most people would want to align some tensor with the principle directions of some layer, which may be undulating through the terrain.  The tensor we expect is, however, in standard x-y-z coordinates.  So to even do terrain following would require some thought (maybe relating the tensor as a function of aspect?)

 

Konstantin and David may have some advice here – they have done more full tensor work in Amanzi than we have in ATS.

Ethan

 

-- 

-------------------------------------------------------------------------

Ethan Coon

Research Scientist

Oak Ridge National Laboratory

 

865-241-1296

https://www.ornl.gov/staff-profile/ethan-t-coon

-------------------------------------------------------------------------

 

--
You received this message because you are subscribed to the Google Groups "Amanzi-ATS Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ats-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ats-users/CAHh7dX7beaQHN9OsYBoisCybz%2BoF%2Bx_QFXKZFU8xU-dg6X3Z9w%40mail.gmail.com.

[Xueyuan Kang]

Hello all,

I’m wondering whether the current version of ATS (either 1.5 or the master branch) supports anisotropic permeability tensors aligned with the principal directions of geological layers, particularly when these layers are undulating and follow the terrain.

Specifically, can ATS automatically compute the tensor components in the global x–y–z coordinate system based on a user-provided anisotropy ratio (e.g., Kxy​/Kz​), rather than requiring the full tensor to be explicitly defined in the XML input?

If this functionality is supported, is there an example or a recommended workflow to implement it?

Thank you very much!

Best regards,
Xueyuan