Add physical 3D Akinci cohesion model - #1260
Conversation
Removed 'Experimental' from GPU support description.
# Conflicts: # docs/src/development.md # docs/src/systems/boundary.md
# Conflicts: # docs/src/systems/fluid.md # src/general/semidiscretization.jl
# Conflicts: # src/callbacks/stepsize.jl # src/schemes/fluid/surface_tension.jl
Codecov Report❌ Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #1260 +/- ##
===========================================
- Coverage 90.13% 79.46% -10.67%
===========================================
Files 138 141 +3
Lines 10797 11758 +961
===========================================
- Hits 9732 9344 -388
- Misses 1065 2414 +1349
Flags with carried forward coverage won't be shown. Click here to find out more. ☔ View full report in Codecov by Harness. 🚀 New features to boost your workflow:
|
|
Closing this draft because the rebuilt Akinci stack intentionally excludes |
Summary
SurfaceTensionAkinciCohesionPhysicalmodel for physical surface tension values in N/mStack
Depends on #1252 and #1024. It also retains the Akinci dependency chain #1259 -> #1258. PR #1024 is not currently contained in #1252, so both heads are explicit ancestors of this branch. The physical model is isolated in commit
3f84349fb; merge commits7a364dbc2and024248c22integrate #1252 and #1024, andbdb842f87tests the physical interface timestep.The logical dependency order is #1252, #1024, #1258, #1259, then this PR. The branch targets
main, so inherited changes remain visible until the dependency PRs merge.Normalization
For physical surface tension
sigma, reference densityrho_0, and compact-support radiush_c, the internal cohesion coefficient isgamma = sigma / ((21 / 7040) * rho_0^2 * h_c^2). The walladhesion_coefficientis a dimensionless multiplier of the same kernel;(1 + cosd(theta)) / 2maps a target contact angle to this ratio.Verification
The complete examples target was not rerun for this stack layer; the affected 3D fluid and wall paths passed through the targeted example.