From cf0250247958e96f4ab7c8f4f26cae756ba3fe88 Mon Sep 17 00:00:00 2001
From: Martin Bauer <martin.bauer@fau.de>
Date: Sat, 8 Sep 2018 20:18:26 +0200
Subject: [PATCH] Phasefield 2D angle testcase
---
phasefield/contact_angle_circle_fitting.py | 8 ++++++++
1 file changed, 8 insertions(+)
diff --git a/phasefield/contact_angle_circle_fitting.py b/phasefield/contact_angle_circle_fitting.py
index 606c5f12..d5bae464 100644
--- a/phasefield/contact_angle_circle_fitting.py
+++ b/phasefield/contact_angle_circle_fitting.py
@@ -94,6 +94,14 @@ def neumann_angles_from_surface_tensions(surface_tensions):
return [np.rad2deg(a) for a in angles]
+def surface_tension_from_kappas(kappas, surface_width):
+ def surface_tensions(i, j):
+ if i == j:
+ return 0
+ return (kappas[i] + kappas[j]) / 6 * surface_width
+ return surface_tensions
+
+
def liquid_lens_neumann_angles(concentration, drop_phase_idx=2, enclosing_phase1=0, enclosing_phase2=1):
"""Assumes a liquid lens setup, where a drop is enclosed between two other phases in the middle of the domain.
--
GitLab