Applying flux boundary condition on sloping grid

[Paulo Herrera]

Hello,

I am working on a model that has a steep slope on one of the sides (attached screenshot shows grid, colors show liquid saturation, green arrow show position where I would like applying a constant flow rate in the vertical direction). What are the options to do this in PFLOTRAN? Which one would be best? Is it necessary to project the flux onto the sloping plane?

Thanks for any suggestion,

Paulo

[Hammond, Glenn E]

Paulo,

 

The (Neumann) flux is always applied orthogonal of the cell faces, including on the slope. Therefore, you must project the vertical flux onto the normal vector orthogonal to the sloped face. The alternative is projecting the area, but if this is an “implicit” unstructured grid (finite element style defined through elements and nodes), PFLOTRAN calculates the area internally (i.e., you have no ability to project the area, only the flux).

 

Glenn

 

From: pflotr...@googlegroups.com <pflotr...@googlegroups.com> On Behalf Of Paulo Herrera
Sent: Sunday, August 4, 2024 5:41 PM
To: pflotr...@googlegroups.com
Subject: [pflotran-dev: 6278] Applying flux boundary condition on sloping grid

 

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

 

--
You received this message because you are subscribed to the Google Groups "pflotran-dev" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pflotran-dev...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/pflotran-dev/CA%2BGTKiOWu2dM2%2BeCPJiAZeNgSOciYmKHGbhaGBNxcHk%3DSJ0uag%40mail.gmail.com.