FEC_BoundaryLayer.py
01f4d6d80e8bmaster

FEC_BoundaryLayer.py
View Options

			import numpy as np
			import matplotlib.pyplot as plt

			def BoundaryLayer():
			x = np.linspace(0,1,100)
			rho = 1000
			mu = 1.7e-3
			V = 2
			def BoundaryLayer(x):
			Re = [rhoVd/mu for d in x]
			delta = []
			for d,Re in zip(x, Re):
			if Re>5e5:
			delta.append(d0.37Re**(-1/5))
			else:
			if d==0:
			delta.append(0)
			else:
			delta.append(5dRe**(-0.5))
			return delta

			delta = BoundaryLayer(x)

			plt.figure(figsize=(8,5))
			plt.rcParams.update({'font.size': 13})
			plt.plot(x,delta, 'b*');
			plt.xlabel("Downstream distance")
			plt.ylabel("Boundary Layer thickness")