Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[GeoMechanicsApplication] Revise water loads conditions #13004

Open
indigocoral opened this issue Jan 13, 2025 · 1 comment
Open

[GeoMechanicsApplication] Revise water loads conditions #13004

indigocoral opened this issue Jan 13, 2025 · 1 comment
Assignees
Labels
GeoMechanics Issues related to the GeoMechanicsApplication Refactor When code is moved or rewrote keeping the same behavior

Comments

@indigocoral
Copy link
Contributor

Currently there are many processes that compute the hydrostatic water pressure:
on a boundary:
2D
ApplyConstantBoundaryHydrostaticPressureProcess
ApplyBoundaryHydrostaticPressureTableProcess
ApplyConstantBoundaryPhreaticLinePressureProcess
ApplyBoundaryPhreaticLinePressureTableProcess
note that the multi line is missing.
3D
ApplyConstantBoundaryPhreaticSurfacePressureProcess
ApplyBoundaryPhreaticSurfacePressureTableProcess

On a domain:
2D
ApplyConstantHydrostaticPressureProcess
ApplyHydrostaticPressureTableProcess
ApplyConstantPhreaticLinePressureProcess
ApplyPhreaticLinePressureTableProcess
ApplyConstantPhreaticMultiLinePressureProcess
ApplyPhreaticMultiLinePressureTableProcess
3D
ApplyConstantPhreaticSurfacePressureProcess
ApplyPhreaticSurfacePressureTableProcess

They all implement projection in gravity of a point onto a phreatic surface and computation of the resulting water pressure.

Proposed unified implementation:
Let the user give a point or series of points on the phreatic surface, gravity direction and water weight ( rho*g ) and possibly a suction cutoff value
2D:

  • sort the points to form a line in horizontal direction
    3D
  • form a triangulated surface ( this is the tricky part, we may need tools for this and I propose to do this later and separately, so handle the 2D first )
    followed by
  • apply phreatic surface coordinate changes due to given tables
  • compute projection and distance of point to the constructed phreatic surface. Let points outside the domain get the nearest given point value. ( don't extrapolate)
  • convert distance to water pressure ( zero or suction cutoff above phreatic surface, hydrostatic below phreatic surface

known hurdle: The json file reader from Kratos does not seem to process a list of point coordinates. In the existing PhreaticMulitLine this is solved in an ugly manner by providing X, Y and Z-coordinates in separate lists.

@indigocoral indigocoral converted this from a draft issue Jan 13, 2025
@indigocoral indigocoral changed the title Revise water loads conditions [GeoMechanicsApplication] Revise water loads conditions Jan 13, 2025
@indigocoral indigocoral added GeoMechanics Issues related to the GeoMechanicsApplication Refactor When code is moved or rewrote keeping the same behavior labels Jan 13, 2025
@indigocoral indigocoral moved this from 👷 In Progress to 📑 Product Backlog in Kratos Product Backlog Jan 13, 2025
@rfaasse
Copy link
Contributor

rfaasse commented Jan 20, 2025

@WPK4FEM will split this parent issue into sub-issues

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
GeoMechanics Issues related to the GeoMechanicsApplication Refactor When code is moved or rewrote keeping the same behavior
Projects
None yet
Development

No branches or pull requests

3 participants