河流动力学大作业solutions to problem 4.pptVIP

Solution to Problem 4 Problem Consider Run 20 of the flume experiments of Coleman (1986), find the best value of β in the modified suspension index for satisfactory agreement between profiles of suspended sediment concentration calculated from the Rouse formula and measured data. Reference: Coleman, N.L. (1986). “Effects of suspended sediment on the open-channel velocity distribution. ”Water Resources Research, AGU, 22(10), 1377-1384. Contents Problem Abstract 1 Basic Model 2 Improved Model 3 Conclutions 4 Referebces 5 Appendix 6 Use Rouse formula to match the measured data in Run 20. Here the modified suspension index is: Problem abstract where Measured Data z(mm) 6 12 18 24 30 46 c(10-3kg/m3) 23 12 8.2 6.1 4.8 3.3 Z(mm) 69 91 122 137 152 162 c(10-3kg/m3) 2.6 1.6 0.76 0.4 0.2 0.11 D(mm) R u*(m/s) T(℃) g(m/s) h(mm) κ 0.105 1.65 0.041 23.9 9.81 170 0.4 According to Coleman’s experiment (Run 20) in 1986, we can get some basic parameters as below. Table 2 Table 1 So,Here we get,? , , 2.Zhang RJ formula 1. empirical formula. After analysis and calculation, now the problem is simplified into a mathematic problem. Where General model 固定β值 yi compare yi与 差异最小时，β的值最准确 判定yi与 差别的指标： x y C/Cb y C/Cb x x β1=0.5 β2=1 sβ=0.5 sβ=1 β1=0.5 β2=0.75 sβ=0.5 sβ=0.75 β1=0.5 β2=0.625 …… β1=0.5 β=0.51619 Because of the precision limitations in C program, this result is not very accurate. Figure 2 Basic method, c language s=0.0019 A preferable way to calculate with Matlab measured data x y 1 1 0.4817073 0.5217391 0.3089431 0.3565217 0.2225610 0.2652174 0.1707317 0.2086957 0.0986214 0.1434783 0.0535525 0.1130435 0.0317609 0.0695652 0.0143942 0.0330435 0.0088125 0.0173913 0.0043325 0.0086957 0.0018067 0.0047826 Here we set the General model, Then we choose the Curve Fitting Tool in Matlab to find the best value of β . Table 6 The calculation is also based on the least square solution and the result is l