Surface Tension Prototype V1

.gitignore

@@ -1,5 +1,7 @@
 .vscode
 .DS_Store
+.opencode
+graphify-out
 Manifest.toml
 out
 out_*

NEWS.md

@@ -6,9 +6,51 @@ used in the Julia ecosystem. Notable changes will be documented in this file for
 
 ## Version 0.5.3
 
-### Features
-
-- Added `flush` keyword argument to `InfoCallback` to flush `stdout` after each output,
+### API Changes
+
+- Akinci cohesion and adhesion kernels now use dimensionally consistent, integral-matched
+  normalizations in 2D. To preserve previous pairwise kernel contributions at compact-support
+  radius `h_c`, multiply the surface-tension coefficient by `627 / (790 * h_c)` and the
+  adhesion coefficient by `42 / (65 * h_c)`; migrated kernel coefficients are
+  resolution-independent.
+- `SurfaceTensionAkinci` now defaults to fluid-only colorfield normals, matching Equation 2
+  of Akinci et al. (2013). Pass a finite `boundary_contact_threshold` to
+  `ColorfieldSurfaceNormal` explicitly to include boundary neighbors.
+- Corrected `SurfaceTensionMorris` to apply its local CSF acceleration once per particle and
+  retain the required one-phase surface delta. Previous coefficients compensated implicitly
+  for a dimensionally incomplete force repeated once per fluid neighbor and must be recalibrated.
+- For Morris and CSS models, `ideal_density_threshold` now denotes a fraction of the continuous
+  complete-support kernel moment instead of an integer neighbor-count fraction. The default zero
+  still disables interior filtering; validation configurations migrate the previous explicit
+  value `0.9` to `0.95`.
+
+### Features
+
+- Reworked `SurfaceTensionMomentumMorris` as a balanced one-phase CSS model. The model now
+  retains the physical color-gradient surface delta, applies a conservative scalar reproducing
+  correction without another neighbor pass, and evaluates stress on demand.
+- Added the explicit `WettedAreaContactAngle` model for validated 3D CSS wetting. It applies the
+  complete derivative of Young's corrected wetted-area energy, including equal-and-opposite wall
+  and rigid-body reactions. The no-contact default is unchanged. The rejected, unshipped
+  geometric-normal and contact-line-force candidates were removed.
+- Added C1 interface activation for Morris CSF and CSS. The color-gradient and continuous
+  support-moment indicators taper the physical surface delta without an extra neighbor pass;
+  CSS retains exact pairwise linear-momentum conservation.
+- Added opt-in `CorrectedCSFSurfaceNormal` for the single-fluid free-surface core of the C-CSF
+  method. It provides renormalized eigenvalue-gradient normals, curvature, and a Shepard-corrected
+  surface delta for `SurfaceTensionMorris`. A finite contact angle enables planar boundary-integral
+  geometry on dummy-particle walls with explicit face measures and normal offsets. Hydrodynamic wall
+  coupling continues to use the configured dummy-particle model.
+- Added opt-in one-pass activity-weighted Shepard smoothing for `ColorfieldSurfaceNormal` directions.
+  Interface activation, surface-delta magnitudes, and particle-shifting normals remain unsmoothed.
+- Added opt-in `FreeSurfaceTangentialShifting` for Sun particle shifting with Morris CSF and CSS.
+  It reuses the smooth color-field interface activity, retains full consistent shifting in the
+  interior, and projects shifting onto the local tangent plane at the free surface. Closed-system
+  shifting defaults are unchanged.
+- Added the opt-in 3D `SurfaceTensionAkinciCohesionPhysical` model. It converts a physical
+  surface tension in N/m to the resolution-dependent Akinci cohesion coefficient, supports
+  same-kernel Young-Dupre wall ratios, and contributes a capillary time-step restriction.
+- Added `flush` keyword argument to `InfoCallback` to flush `stdout` after each output,
   useful for monitoring progress in real-time on clusters or batch systems (#1246).
 - Added the computation of boundary normals for `RectangularTank`s and `SphereShape`s.
 

Project.toml

@@ -35,6 +35,7 @@ WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 
 [weakdeps]
 CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
+Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
 OrdinaryDiffEqCore = "bbf590c4-e513-4bbe-9b18-05decba2e5d8"
 OrdinaryDiffEqLowStorageRK = "b0944070-b475-4768-8dec-fb6eb410534d"
 OrdinaryDiffEqSymplecticRK = "fa646aed-7ef9-47eb-84c4-9443fc8cbfa8"
@@ -43,6 +44,7 @@ Static = "aedffcd0-7271-4cad-89d0-dc628f76c6d3"
 [extensions]
 TrixiParticlesOrdinaryDiffEqSymplecticRKExt = ["OrdinaryDiffEqSymplecticRK", "OrdinaryDiffEqCore"]
 TrixiParticlesCUDAExt = "CUDA"
+TrixiParticlesMakieExt = "Makie"
 
 [compat]
 Accessors = "0.1.43"
@@ -60,6 +62,7 @@ GPUArraysCore = "0.2"
 JSON = "1"
 KernelAbstractions = "0.9"
 LinearAlgebra = "1"
+Makie = "0.24"
 OrdinaryDiffEqLowStorageRK = "3"
 OrdinaryDiffEqCore = "4"
 OrdinaryDiffEqSymplecticRK = "2"

README.md

@@ -45,24 +45,24 @@ It offers intuitive configuration, robust pre- and post-processing, and vendor-a
 We provide several example simulation setups in the `examples` folder (which can be accessed from Julia via `examples_dir()`).
 
 <table align="center" border="0">
-  <tr>
-  </tr>
   <tr>
     <td align="center">
-      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/683e9363-5705-49cc-9a5c-3b47d73ea4b8" style="width: 80% !important;"/><br><figcaption>2D Dam Break</figcaption>
+      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/683e9363-5705-49cc-9a5c-3b47d73ea4b8" style="width: 80% !important;"/><br />
+      <figcaption>2D Dam Break</figcaption>
     </td>
     <td align="center">
-      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/c10faddf-0400-47c9-b225-f5d286a8ecb8" style="width: 80% !important;"/><br><figcaption>Moving Wall</figcaption>
+      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/c10faddf-0400-47c9-b225-f5d286a8ecb8" style="width: 80% !important;"/><br />
+      <figcaption>Moving Wall</figcaption>
     </td>
   </tr>
-  <tr>
-  </tr>
   <tr>
     <td align="center">
-      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/e05ace63-e330-441a-a391-eda3d2764074" style="width: 80% !important;"/><br><figcaption>Oscillating Beam</figcaption>
+      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/e05ace63-e330-441a-a391-eda3d2764074" style="width: 80% !important;"/><br />
+      <figcaption>Oscillating Beam</figcaption>
     </td>
     <td align="center">
-      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/ada0d554-e0ba-44ed-923d-2b77ef252258" style="width: 80% !important;"/><br><figcaption>Dam Break with Elastic Plate</figcaption>
+      <img src="https://github.com/trixi-framework/TrixiParticles.jl/assets/10238714/ada0d554-e0ba-44ed-923d-2b77ef252258" style="width: 80% !important;"/><br />
+      <figcaption>Dam Break with Elastic Plate</figcaption>
     </td>
   </tr>
 </table>

compare_akinci/.gitignore

@@ -0,0 +1,2 @@
+/Manifest.toml
+!/*.png