Skip to content

Use Darcy Velocity in entropy_viscosity.cc#6937

Open
danieldouglas92 wants to merge 4 commits intogeodynamics:mainfrom
danieldouglas92:use_darcy_velocity_for_diffusion
Open

Use Darcy Velocity in entropy_viscosity.cc#6937
danieldouglas92 wants to merge 4 commits intogeodynamics:mainfrom
danieldouglas92:use_darcy_velocity_for_diffusion

Conversation

@danieldouglas92
Copy link
Copy Markdown
Contributor

When a field is advected with the darcy field advection method, I realized when doing some tests where I changed the stabilization parameters cR and beta that they didn't seem to have much of an influence on the diffusion of the porosity composition, especially when compared to the same test with the porosity advected with the McKenzie equations. It turns out that the Darcy velocity was never computed when calculating the entropy viscosity, and so the solid velocity was always used.

This PR modifies entropy_viscosity.cc to account for the Darcy velocity when advecting a field with the darcy field advection method.

@danieldouglas92 danieldouglas92 force-pushed the use_darcy_velocity_for_diffusion branch from 29e2baf to 631fe38 Compare April 21, 2026 18:56
@danieldouglas92
Copy link
Copy Markdown
Contributor Author

Tests that are failing all use the darcy field advection method, so this is expected. Once the implementation gets approved I'll update the tests.

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.

1 participant