summaryrefslogtreecommitdiff
path: root/3434/CH9/EX9.1.ii/Ex9_1_ii.sce
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /3434/CH9/EX9.1.ii/Ex9_1_ii.sce
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '3434/CH9/EX9.1.ii/Ex9_1_ii.sce')
-rw-r--r--3434/CH9/EX9.1.ii/Ex9_1_ii.sce19
1 files changed, 19 insertions, 0 deletions
diff --git a/3434/CH9/EX9.1.ii/Ex9_1_ii.sce b/3434/CH9/EX9.1.ii/Ex9_1_ii.sce
new file mode 100644
index 000000000..45266d0e5
--- /dev/null
+++ b/3434/CH9/EX9.1.ii/Ex9_1_ii.sce
@@ -0,0 +1,19 @@
+clc
+// given data
+G=39.0 // temperature gradient in K/km.
+h2=10.0 // depth in km
+rhor=2700.0 // kg/m^3
+cr=820.0 // in J/kg-K
+QbyA=0.5 //water flow rate in m^3/sec-km^2
+rhow=1000.0 // density of water in kg/m^3
+cw=4200.0 // specific heat of water in J/kg-K
+h1=120.0/G // T1-T0=120 K is given
+h21=h2-h1 // in km
+t=25 // time in years
+
+thtao=G*h21/2.0 // in degree K
+printf( "Useful initial temp is %.2f degree K",thtao)
+tau=rhor*cr*h21*(1000**3)/(QbyA*rhow*cw) // in seconds
+tau=tau/(2*60*60*24*365) // in years
+thta=thtao*exp(-t/tau) // in degree Kelvin
+printf( " \n Useful average temp after 25 years is %.2f degree K",thta)