Skip to content

Fix higher-order Darcy flux recovery - #3

Merged
mmegally merged 1 commit into
mmegally:darcy-solverfrom
zasexton:review/pr-620-latest
Sep 15, 2026
Merged

mmegally merged 1 commit into
mmegally:darcy-solverfrom
zasexton:review/pr-620-latest

Conversation

@zasexton

Copy link
Copy Markdown

Summary

Corrects singular Darcy flux recovery at TRI6 corner nodes in SimVascular/svMultiPhysics#620. Uses an element-local consistent L2 projection and volume-weighted nodal averaging in post.cpp, while preserving the existing TRI3/TET4 recovery.

Validation

  • 13 manual solver runs passed across serial and three MPI ranks, covering TRI6/TET10, curved and small elements, and centroid TET4 quadrature.
  • Original TRI6 flux error reduced from 5.1e14 to 8.9e-16.

@mmegally
mmegally merged commit 28f48ec into mmegally:darcy-solver Sep 15, 2026
@zasexton
zasexton deleted the review/pr-620-latest branch September 15, 2026 03:59
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants