FreeConvection.py
57feb8539944master

FreeConvection.py
View Options

			import numpy as np
			import matplotlib.pyplot as plt


			def plotVerticalPLate():

			g = 9.81;
			T_s = 80;
			T_inf = 20;
			DT = T_s - T_inf;
			T_f = (T_s + T_inf)/2;
			N = 1.5;
			beta = 1/T_f;
			mu = 1.957e-5;
			rho = 1/0.9276;

			nu = mu/rho;

			def Gr(t) :
			return gbeta(T_s - T_inf)np.power(t, 3)/(nunu)

			x=np.linspace(0.0000001,40,100)
			y= Nnp.multiply(np.power(0.25Gr(x),-0.25),x)

			plt.figure(figsize=[5,7])
			plt.ylim(0,40)
			plt.xlim(0,0.03)
			plt.xlabel('Boundary Layer Thickness')
			plt.ylabel('Plate')
			plt.fill_between(y,x,40, facecolor = 'red')
			plt.plot(y,x, fillstyle ='left', c = 'red')