blob: d7be0dab33a0ce971c312b687bfce0d68568abf4 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
|
// A Textbook of Fluid Mecahnics and Hydraulic Machines - By R K Bansal
// Chapter 1-Properties of Fluid
// Problem 1.3
y=poly(0,"y")
u=2/3*y-y^2 // Defining the Velocity Function
a=derivat(u) //Taking Derivative of the Velocity
visc=8.63/10 //Converting Dynamic Viscosity from poise to N s/m^2
ss1=visc*horner(a,0) //Shear stress=(Dynamic viscosity *Velocity Gradient) at y=0
mprintf("The shear stress at y=0 is %f N/m^2 \n",ss1)
ss2=visc*horner(a,0.15) //Shear stress=(Dynamic viscosity *Velocity Gradient) at y=0.15
mprintf("The shear stress at y=0.15 is %f ",ss2)
|
