Page MenuHomec4science
Files F100691897

FEC_BoundaryLayer.py
No OneTemporary
Actions

File Metadata

Created
Sat, Feb 1, 22:24

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 = [rho*V*d/mu for d in x]
delta = []
for d,Re in zip(x, Re):
if Re>5e5:
delta.append(d*0.37*Re**(-1/5))
else:
if d==0:
delta.append(0)
else:
delta.append(5*d*Re**(-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")

Event Timeline

c4science · Help