DynamicsPARTA-Ted

BULTI ABEBE TADESSE ID.No. D2015042 Submitted to Assoc. Prof. Xi Li HOHAI UNIVERSITY Nanjing china EXRCISE ON SEDIMENT DYNAMICS Sediment dynamics Exercise application on Mat lab 1. Write MATLAB CODE and Plot the graph of settling velocity versus sediment concentration based on the following given condition; () 2m 1.5n 4b 0.15a ion,concentratinitial0.1,C Given tcoefficienempiricaldependentsedimentm tcoefficienempiricaldependentsedimentisn tcoefficienempiricaldependentsedimentisb tcoefficienempiricaldependentsedimentisa (g/l)literpergramionconcentratsedimentisC (m/s)secondpermeterinvelocitysettlingisWWhere 0 s 22 = = = = = + = is bC aC W m n s 2.00 0.00 5.00 15.00 30.00 60.00 120.00 180.00 1.55 4.60 0.72 0.40 0.30 0.00 0.00 0.00 1.30 4.60 4.00 0.40 0.30 0.24 0.22 0.20 1.05 4.60 4.60 0.40 0.30 0.24 0.22 0.20 0.80 4.60 4.60 2.00 0.30 0.24 0.22 0.20 0.55 4.60 4.60 4.60 0.30 0.24 0.22 0.20 0.30 4.60 4.60 4.60 0.31 0.24 0.22 0.20 0.15 4.60 4.60 13.00 23.00 14.00 22.00 0.30 0.05 4.60 19.00 30.00 32.00 38.00 43.00 50.00 1.1. Tutorial 2.4 for data input Running the mat lab file “set1.m”: Settling Velocity and Concentration Profile Simulations Number of elevation for settling column concentration sampling, nele, is: nele = 8 Number of times for concentration sampling, ntime, is: Sediment dynamics Exercise application on Mat lab ntime = 7 Input file name is: setm2.txt Estimated coefficient, m, is: m = 2.0 Estimated coefficient, n, is: n = 1.4 Estimate suspension concentration for maximum settling velocity, C2 (kg/m3), is: C2 = 6.0 Calculate coefficient, b, based on equation (4.18) page 4-28 is: b = 8.176622 Estimate maximum settling velocity, Ws2 (m/s), is: W Ws2 = 0.02 Calculate coefficient, a, based on equation (4.17) page 4-26 is: a = 0.043003 Do you want to replot for y-axis range? Y/N Y: Y Minimum settling velocity for velocity-concentration plot: y_min = 1e-6 Maximum settling velocity for velocity-concentration plot: y_max = 1e-2 Sediment dynamics Exercise application on Mat lab Do you want to estimation of the four parameters for settling velocity formula? Y/N Y: Y Estimated coefficient, a, is: a = 0.15 Estimated coefficient, b, is: b = 4.0 Estimated coefficient, a, is: m = 2.0 Estimated coefficient, n, is: n = 1.5 Limiting free settling concentration,C1 (kg/m3), is: C1 = 0.2 % MAT profile %* * a=impor b=a(2:9 c=a(2:9 liter (g/l x=logsp semilog axis(1e axis hold on title(Cu xlabel(S ylabel(E hold off TLAB code * * rt data(setm ,1); % Thes ,2:8); %The l) pace(-2,3); % gx(c,b,-o); % e-2,1e3,0,2) urves of Tim Sediment Co Elevation (m f to import da * * m2.txt); % im e are array o ese are array % This is x-a % plots seve );% The lim me evolution oncentration m); ata and plot * mports file s of elevation y of sedimen axis in log sc en curves for it of sedime of sediment n (g/l); Sed curves of th * set2.txt from data as y co t concentrat cale r elevation a nt concentra t concentrati iment dynam he Time-evo * m other folde oordinate val tion data as x against sedim ation on x-ax ion profile) mics Exercise olution of sed * er location lues in mete x coordinate ment concen xis and the l ; e application diment conc * ers (m) in lin e values in g ntration data limit of elev n on Mat lab centration * near scale gram per vation on y- b Sediment dynamics Exercise application on Mat lab % MATLAB CODE to plot the curve for sediment size against Percentage finer by weight % It is also called size distribution curve of sediments or Grading curve %* * x=0.1 0.05 0.02 0.01 0.005 0.001; % Size of sediment , d in microns y=100 94.68 72.10 41.82 27.13 11.69; % Percentage finer by weight in (%) x1=log10(x); % vectors containing values of size of sediment in microns y1=y; % vectors containing values of percent finer by weight xx=-3.0:0.04:-1;% points in the vector xx pp=spline(x1, y1); % spline function for cubic interpolation y2=ppval(pp,xx); % ppvalt function returns the coefficients of a polynomial that best fits the data, x2=10.xx; semilogx(x,y,-,x2,y2) hold on title(Size distribution curves of sediments initially in suspension and at steady state ); xlabel(Size, d (microns) ylabel(Percentage finer by weight (%) hold off Sediment dynamics Exercise application on Mat lab %suspended sediment concentration variation. clear clc disp(total time); dd=input(dd=); disp(Time Step); dt=input(dt=); ff1=(0.1-y)*1.0*0.5027*sin(0.5027*t)+y*0.15*y.1.5/(y.2+4.2).2)/(5.0- 1.0*cos(0.5027*t); f1=inline(ff1,t,y); n=dd*24/dt; cc=0.0