Unstructured grid: assigning coordinates/vertices to regions

[Swamini Khurana]

Hello

I see in the information on the REGION input card that it is possible to assign faces defined by vertex ID/node IDs and face type (triangle, quadrilateral etc.) to a particular region. It is also possible to route this information to pflotran using an external file (*.ss or *.h5).

I was wondering if it is possible to assign coordinates/vertex IDs to a particular region using *.ss or *.h5 files (or any other external file format)? I would like to assign boundary conditions at specific vertices actually, but I am not sure how to do that without assigning these vertices to a particular region. Is there an example problem that I could use for this?

Thank you for your help!

Regards

Swamini

[Hammond, Glenn E]

Vertices have to be assigned to a region. 

 

Glenn

 

From: pflotra...@googlegroups.com <pflotra...@googlegroups.com> On Behalf Of Swamini Khurana
Sent: Wednesday, February 24, 2021 3:32 AM
To: pflotran-users <pflotra...@googlegroups.com>
Subject: [pflotran-users: 6157] Unstructured grid: assigning coordinates/vertices to regions

 

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

 

--
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/0914dc1f-d326-4973-938f-20f7ba86a61en%40googlegroups.com.

[Swamini Khurana]

Hi Glenn,

Yes, I assigned vertices using vertex IDs to a region "bc" in the HDF5 file (feflowtoh5_grid.h5 file here). The error I get is:

"ERROR: Dimension of Regions/bc/Vertex Ids"

Regions/bc/Vertex Ids is a 1D array as I am not assigning the boundary condition by cell/faces but by vertices.

I also then tried to give only the coordinates of the vertices, but I got the error:

"ERROR: Regions/bc is not defined by "Cell Ids" or  "Face Ids" or  "Vertex Ids".

In such a case, what is the dimension that is expected while defining a region by vertices?

Regards

Swamini

[Hammond, Glenn E]

Swamini,

 

From: http://doc-dev.pflotran.org/user_guide/cards/subsurface/region_card.html

 

UNSTRUCTURED Grid (Implicit):

Specify boundary faces through a list of vertex IDs in one of the following formats:

·         ASCII (.ss): each line specifies a face element type (T=triangle and Q=quadrilateral) and associated vertex IDs

·         HDF5 (.h5): 2D array of integers where the vertex IDs for each face are in the short dimension and the long dimension equals the number of faces.

 

 

Take a look at how vertex ids are listed in https://bitbucket.org/pflotran/pflotran/src/master/regression_tests/surfaceflow/inlinesurface/simple.h5

 

The dataset is 2D.  There are 4 vertex ids for each row where the vertex ids define the triangle or quad face.  In the case of a triangle, three ids with a 0 at the fourth.  The number of rows corresponds to the number of faces.  There is no need to define the coordinates as the unstructured grid file defines coordinates.

 

Glenn

 

To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-users/6f8d2fa6-8f32-480c-939b-74e567647276n%40googlegroups.com.

[Swamini Khurana]

Hi Glenn,

Thank you. So I understand that I cannot assign boundary condition to specific points (vertex) or to specific cells, but only to a face of a cell.

Regards

Swamini