diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /534/CH14/EX14.7/14_7_Drug_Medication.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '534/CH14/EX14.7/14_7_Drug_Medication.sce')
| -rw-r--r-- | 534/CH14/EX14.7/14_7_Drug_Medication.sce | 46 |
1 files changed, 46 insertions, 0 deletions
diff --git a/534/CH14/EX14.7/14_7_Drug_Medication.sce b/534/CH14/EX14.7/14_7_Drug_Medication.sce new file mode 100644 index 000000000..150f3e3d5 --- /dev/null +++ b/534/CH14/EX14.7/14_7_Drug_Medication.sce @@ -0,0 +1,46 @@ +clear;
+clc;
+printf('FUNDAMENTALS OF HEAT AND MASS TRANSFER \n Incropera / Dewitt / Bergman / Lavine \n EXAMPLE 14.7 Page 913 \n')// Example 14.7
+
+// Total dosage of medicine delivered to the patient over a one-week time period, sensivity of the dosage to the mass duffusivity of the patch and skin
+
+Dap = .1*10^-12 ;//[m^2/s] Diffusion coefficient of medication with patch
+Das = .2*10^-12 ;//[m^2/s] Diffusion coefficient of medication with skin
+L = .05 ;//[m] patch Length
+rhop = 100 ;//[kg/m^3] Density of medication on patch
+rho2 = 0 ;//[kg/m^3] Density of medication on skin
+K = .5 ;//Partition Coefficient
+t = 3600*24*7 ;//[s] Treatment time
+
+//Applying Conservation of species equation 14.47b
+//By analogy to equation 5.62, 5.26 and 5.58
+D = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das*Dap*t)/(sqrt(Das)+sqrt(Dap)/K);
+
+printf('\n Total dosage of medicine delivered to the patient over a one-week time period is %.1f mg',D*10^6);
+
+//Senstivity of dosage to the patch and skin
+clf();
+//Subplot 1
+Dap1 = .1*10^-12 ;//[m^2/s]
+Das1 = .1*10^-12 ;//[m^2/s]
+Das2 = .2*10^-12 ;//[m^2/s]
+Das3 = .4*10^-12 ;//[m^2/s]
+x = linspace(0,7,50);
+y1 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das1*Dap1*3600*24*x)/(sqrt(Das1)+sqrt(Dap1)/K)*10^6;
+y2 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das2*Dap1*3600*24*x)/(sqrt(Das2)+sqrt(Dap1)/K)*10^6;
+y3 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das3*Dap1*3600*24*x)/(sqrt(Das3)+sqrt(Dap1)/K)*10^6;
+subplot(1,2,1);
+plot(x,y1,x,y2,x,y3);
+xtitle("Dosage vs Time-period at Dap = .1*10^ -12 (m^2/s)", "Day", "Dosage (mg)");
+legend (".1*10^12", ".2*10^12", ".4*10^12");
+
+//Subplot 2
+Dap2 = .01*10^-12 ;//[m^2/s]
+yn1 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das1*Dap2*3600*24*x)/(sqrt(Das1)+sqrt(Dap2)/K)*10^6;
+yn2 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das2*Dap2*3600*24*x)/(sqrt(Das2)+sqrt(Dap2)/K)*10^6;
+yn3 = 2*rhop*L^2/(sqrt(%pi))*sqrt(Das3*Dap2*3600*24*x)/(sqrt(Das3)+sqrt(Dap2)/K)*10^6;
+subplot(1,2,2);
+plot(x,yn1,x,yn2,x,yn3);
+xtitle("Dosage vs Time-period at Dap = .01*10^ -12 (m^2/s)", "Day", "Dosage (mg)");
+legend (".1*10^12", ".2*10^12", ".4*10^12");
+//END
\ No newline at end of file |