feat: ✨ compute Leverett J-function with CapPressure class

ImperialCollegeLondon · Oct 29, 2024 · 312d54d · 312d54d
1 parent 08654cc
commit 312d54d
Showing 1 changed file with 16 additions and 6 deletions.
diff --git a/src/CapPressure.m b/src/CapPressure.m
@@ -12,7 +12,7 @@
             if nargin==0
                 return;
             end
-            
+
             obj.contact_angle = contact_angle;
             obj.surface_tension = surface_tension;
             obj.leverett_j = leverett_j;
@@ -38,20 +38,30 @@
             arguments
                 obj (1,1) CapPressure
                 pc  (1,1) double
-                poro double {mustBePositive}
+                poro double {mustBeNonnegative}
                 perm double {mustBeNonnegative}
             end
 
-            perm = obj.transform_perm(poro,perm);
-
-            lj = pc ./ obj.mult .* sqrt(perm./poro);
-
+            lj = obj.inv_lj(pc,poro,perm);
             lj = min(lj,obj.leverett_j.data(1,2));
             lj = max(lj,obj.leverett_j.data(end,2));
 
             sw = obj.leverett_j.inv(lj);
         end
 
+        function lj = inv_lj(obj,pc,poro,perm)
+            arguments
+                obj (1,1) CapPressure
+                pc  double
+                poro double {mustBeNonnegative}
+                perm double {mustBeNonnegative}
+            end
+
+            perm = obj.transform_perm(poro,perm);
+
+            lj = pc ./ obj.mult .* sqrt(perm./poro);
+        end
+
         function dpc_dsw = deriv(obj,sw,poro,perm)
             arguments
                 obj (1,1) CapPressure