Changed mu to nu

Author: hkjeldsberg

find_baseflow.py

@@ -12,7 +12,7 @@ def solve_navier_stokes(mesh, boundaries, nu, pin, pout, inflow_marker, outflow_
     Args:
         mesh (Mesh): Mesh of problem domain
         boundaries (MeshFunction): Function determining the boundaries of the mesh
-        nu (float): Dynamic viscosity
+        nu (float): Kinematic viscosity
         pin (float): Predefined pressure value at inlet(s)
         pout (float): Predefined pressure value at outlet(s)
         inflow_marker (list): ID(s) corresponding to inlet(s) of mesh
@@ -161,7 +161,7 @@ def D(u):
 if __name__ == '__main__':
     pin = 2.0
     pout = 1.0
-    mu = 1.0
+    nu = 1.0
     radius = 1.0
     n_elem = 15
 
@@ -185,7 +185,7 @@ def D(u):
         inflow_marker, outflow_marker, no_slip_marker = get_marker_ids(case)
 
     # Solve for NS flow in mesh domain
-    u, p = solve_navier_stokes(mesh, boundaries, mu, pin, pout, inflow_marker, outflow_marker, no_slip_marker)
+    u, p = solve_navier_stokes(mesh, boundaries, nu, pin, pout, inflow_marker, outflow_marker, no_slip_marker)
 
     file = File("Baseflow/u.pvd")
     file << u

find_perturbation.py

@@ -255,7 +255,7 @@ def get_eigenvalue_matrices(mesh, nu, u_init):
 
     Args:
         mesh (Mesh): Mesh of problem domain
-        nu (float): Dynamic viscosity
+        nu (float): Kinematic viscosity
         u_init (Function): Baseflow solution
 
     Returns:
@@ -305,11 +305,11 @@ def all_boundaries(x, on_boundary):
     return A, B, W
 
 
-def EIG_A_cyl(u, p, v, q, mu):
+def EIG_A_cyl(u, p, v, q, nu):
     """
     Set up operators for eigenvalue problem (Left hand side)
     """
-    eig_a = mu * inner(grad(u), grad(v)) * dx - div(v) * p * dx - q * div(u) * dx
+    eig_a = nu * inner(grad(u), grad(v)) * dx - div(v) * p * dx - q * div(u) * dx
     return eig_a
 
 
@@ -332,7 +332,7 @@ def read_command_line():
                                        '1 (C0015_terminal), 2 (C0019), 3 (C0065_healthy), 4 (C0065_healthy)',
                         default='0', type=str, choices={"0", "1", "2", "3", "4", "poise"})
     parser.add_argument('--delta_p', help='Pressure drop in mmHg', default=5, type=float)
-    parser.add_argument('--nu', help='Dynamic viscosity', default=3e-6, type=float)
+    parser.add_argument('--nu', help='Kinematic viscosity', default=3e-6, type=float)
     args = parser.parse_args()
 
     return args