Commit 4788789a authored by Stefan Hiemer's avatar Stefan Hiemer
Browse files

Introduction of effective viscosity.

parent dfe7155f
Loading
Loading
Loading
Loading
+6 −0
Original line number Diff line number Diff line
@@ -51,6 +51,12 @@ def PS_settings():
    
    return Th,T0

def PP_settings():
    Th =  # K
    T0 = 293.15 # K
    
    return Th,T0

def vft_forces(Us,
               material,
               geometry,
+1 −1
Original line number Diff line number Diff line
@@ -95,7 +95,7 @@ def F_newt_isoth(km, Th, Tm, rho_s, rho_m, h, Cs, T0,
    delta = get_delta(km, Th, Tm, rho_s, Us, h, Cs, T0)
    
    if debug:
        print("alpha",alpha,
        print("alpha",alpha,"\n"
              "U0",U0,
              "delta",delta)