From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 914/CH14/EX14.1/ex14_1.sce | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) create mode 100755 914/CH14/EX14.1/ex14_1.sce (limited to '914/CH14/EX14.1/ex14_1.sce') diff --git a/914/CH14/EX14.1/ex14_1.sce b/914/CH14/EX14.1/ex14_1.sce new file mode 100755 index 000000000..7b390666b --- /dev/null +++ b/914/CH14/EX14.1/ex14_1.sce @@ -0,0 +1,21 @@ +clc; +warning("off"); +printf("\n\n example14.1 - pg726"); +// given +T=40+273.15; //[K] - temperature +P=1; //[atm] - pressure +sigma=3.711*10^-10; //[m] +etadivkb=78.6; //[K] +A=1.16145; +B=0.14874; +C=0.52487; +D=0.77320; +E=2.16178; +F=2.43787; +Tstar=T/(etadivkb); +// using the formula si=(A/(Tstar^B))+(C/exp(D*Tstar))+(E/exp(F*Tstar) +si=(A/(Tstar^B))+(C/exp(D*Tstar))+(E/exp(F*Tstar)); +M=28.966; //[kg/mole] - molecular weight +// using the formula mu=(2.6693*(10^-26))*(((M*T)^(1/2))/((sigma^2)*si)) +mu=(2.6693*(10^-26))*(((M*T)^(1/2))/((sigma^2)*si)); +printf("\n\n The viscosity of air is \n mu=%eNs/m^2=%fcp",mu,mu*10^3); -- cgit