def compute_tbl_properties(y, uMean, nu, flip):
"""Compute various parameters of a TBL."""
y = y[np.nonzero(y)]
uMean = uMean[np.nonzero(uMean)]
if flip:
y = np.flipud(y)
uMean = np.flipud(uMean)
theta = momentum_thickness(y, uMean)
delta = delta_99(y, uMean)
deltaStar = delta_star(y, uMean)
uTau = np.sqrt(nu*uMean[1]/y[1])
u0 = uMean[-1]
yPlus1 = y[1]*uTau/nu
return theta, deltaStar, delta, uTau, u0, yPlus1
评论列表
文章目录