Journal of Hydrology
|
|
- Alexander Fox
- 5 years ago
- Views:
Transcription
1 Journal of Hydrology 59 (24) Contents lsts avalable at ScenceDrect Journal of Hydrology journal homepage: Hgh resoluton numercal schemes for solvng knematc wave equaton Chunshu Yu, Jennfer G. Duan,2 Department of Cvl Engneerng and Engneerng Mechancs, Unversty of Arzona, USA artcle nfo summary Artcle hstory: Receved 3 July 23 Receved n revsed form 3 June 24 Accepted 2 August 24 Avalable onlne 3 August 24 Ths manuscrpt was handled by Konstantne P. Georgakakos, Edtor-n-Chef, wth the assstance of Ioanns K Tsans, Assocate Edtor Keywords: Knematc wave equaton Godunov-type scheme MacCormack scheme Shock wave Rarefacton wave Ranfall-runoff overland flow Ths paper compares the stablty, accuracy, and computatonal cost of several numercal methods for solvng the knematc wave equaton. The numercal methods nclude the second-order MacCormack fnte dfference scheme, the MacCormack scheme wth a dsspatve nterface, the second-order fnte volume scheme, and the ffth-order fnte volume scheme. These numercal schemes are tested aganst several synthetc cases and an overland flow experment, whch nclude shock wave, rarefacton wave, wave steepenng, unform/non-unform ranfall generated overland flows, and flow over a channel of varyng bed slope. The results show that the MacCormack scheme s not a Total Varaton Dmnshng (TVD) scheme because oscllatory solutons occurred at the presence of shock wave, rarefacton wave, and overland flow over rapdly varyng bed slopes. The MacCormack scheme wth a dsspatve nterface s free of oscllaton but wth consderable dffusons. The Godunov-type schemes are accurate and stable when dealng wth dscontnuous waves. Furthermore the Godunov-type schemes, lke and scheme, are needed for smulatng surface flow from spatally non-unformly dstrbuted ranfalls over rregular terrans usng moderate computng resources on current personal computers. Ó 24 Elsever B.V. All rghts reserved.. Introducton The knematc wave equaton was frst developed by Lghthll and Whtham (955). The equaton s based on the assumptons that the acceleraton term and the pressure gradent term n the momentum equaton are neglgble, so that the energy slope s equal to the bottom slope. The knematc wave model s commonly used to smulate the overland flow (Ponce, 99; Sngh, 2). Henderson (966) showed that natural flood waves behave nearly the same as the knematc wave n steep slopes (S > :2). Vera (93) concluded that the knematc wave equaton can be used on natural slopes wth the knematc wave number, k 5. Ponce (99) compared the knematc wave equaton wth the unt hydrograph as a practcal method of overland flow routng. Sngh (2) concluded that the knematc wave equaton s applcable to surface water routng, vadose zone hydrology, rverne and costal processes, eroson and sedment transport, etc. The knematc wave equaton s a frst-order hyperbolc partal dfferental equaton (PDE). For a hyperbolc equaton, the Correspondng author. E-mal addresses: chunshu@emal.arzona.edu (C. Yu), gduan@emal.arzona.edu (J.G. Duan). Post-doc Research Assocate. 2 Assocate Professor. dsturbance wll travel along the characterstcs of the equaton n a fnte propagaton speed. Ths feature dstngushes the hyperbolc equatons from ellptc and parabolc equatons. On the other hand, the knematc wave equaton also belongs to a knd of equatons called conservaton laws (LeVeque, 22; Toro, 29). Snce the flux term s a nonlnear functon of conservatve varables, the soluton does not propagate unformly but deforms as t evolves. Even the ntal condtons are contnuous and smooth, the hyperbolc conservaton laws can develop dscontnutes n the soluton, for example, shock waves. Both shock and rarefacton waves are the ntrnsc features of hyperbolc equatons. Lghthll and Whtham (955) dscussed the formatons of shock wave and rarefacton wave. Kbler and Woolhser (97) nvestgated the structure and general propertes of shock waves and developed a numercal procedure for shock fttng. Eagleson (97) found that usng non-unform flow depth as ntal condton, non-unform ranfall n the source term, or ncreasng nflows as the boundary condton may cause the formaton of knematc shock wave. Borah et al. (9) presented the propagatng shock-fttng scheme (PSF) to smulate overland flow wth shock waves. Sngh (2) found three factors that affect the shock wave formaton: () ntal and boundary condtons; (2) lateral nflow and outflow, and (3) watershed geometrc characterstcs. Due to the complex geometry, non-unform roughness /Ó 24 Elsever B.V. All rghts reserved.
2 24 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Nomenclature c wave celerty (m/s) f nfltraton rate (mm/h) h flow depth (m) h flow depth at the center of cell (m) h +/2 reconstructed depth at the nterface + /2 of cell (m) h L, h R flow depth at the left and rght of wave front (m) ranfall excess (mm/h) ntensty of ranfall (mm/h) L(h) operator L(h) defned n Eq. (2) ( ) L channel length (m) m exponental n Eq. (9) ( ) n Mannng s roughness coeffcent (sm /3 ) q dscharge per unt wdth (m 2 /s) q L outflow dscharge (m 2 /s) r +/2 dstance from the cell center to the nterface + /2 (m) S shock wave speed (m/s) S bed slope ( ) TV(h) total varaton of flow depth ( ) t smulaton tme (s) t r duraton of ranfall (s) tme of concentraton (s) t c t p v x Subscrpt Superscrpt tme defned n Eq. (3) (s) flow velocty (m/s) downslope dstance (m) spatal ndex n, n þ, and n + and n +, temporal ndces Greek a coeffcent n Eq. (9) ( ) c j lnear weghts ( ) h coeffcent for the dsspatve nterface ( ) e = 6 truncaton error for actual calculatons ( ) Dt tme step (s) Dx space step (m). rh lmted depth gradent of cell ( ) rh L, rh R lmted depth gradent at the left and rght of cell nterface ( ) and non-unform ranfall pattern, t s mpossble to derve a general analytcal soluton for the knematc wave equaton. Sngh (2) summarzed three numercal technques for solvng the knematc wave equaton: () method of characterstc, (2) Lax Wendroff fnte dfference method, and (3) fnte element method. Numercal dffuson and numercal dsperson were observed when usng the fnte dfference schemes (Ponce, 99). Kazezyılmaz- Alhan et al. (25) evaluated several fnte dfference schemes for solvng knematc wave equaton: the lnear explct scheme, the four-pont Pressmann mplct scheme, and the MacCormack scheme. The study (Kazezyılmaz-Alhan et al., 25) found the Mac- Cormack scheme s better than the four-pont mplct fnte dfference scheme for shock capture. However, Kazezyılmaz-Alhan et al. (25) dd not explctly examne the dsperson occurred at the shock and rarefacton waves from non-unformly dstrbuted ranfall. The stablty of the classcal MacCormack scheme at the presence of shock and rarefacton wave remans unknown. Recently, the Godunov-type fnte volume method has been wdely used n solvng shallow water equatons (LeVeque, 22; Toro, 29) because of ts wde applcablty, strong stablty, and hgh accuracy. One of the most popular Godunov-type methods s a second-order, TVD (Total Varaton Dmnshng) scheme, namely the (Monotone Upstream-centered Schemes for Conservaton Laws) scheme (van Leer, 979). The scheme s a hghresoluton scheme because () the spatal accuracy of the scheme s equal to or hgher than second order; (2) the scheme s free from numercal oscllatons or wggles; (3) hgh-resoluton s produced around dscontnutes. In general, the hgh-resoluton schemes are consdered as tradeoffs between computatonal cost and desred accuracy (Harten, 93; Toro, 29). Another popular but relatvely new method s the hgh-order (Weghted Essentally Non-Oscllatory) fnte volume scheme (Shu, 999). Hghorder means the order of accuracy s equal to or hgher than the thrd-order. Accordng to Shu (29), the scheme s sutable for the complcated problems, such as flow havng both shocks and complcated smooth structures (e.g., small perturbaton). Although the computatonal cost of hgh-order scheme can be three to ten tmes than a second-order hgh-resoluton scheme, t s stll preferable because of ts hgh-order accuracy n both tme and space. The applcatons of those two hgh resoluton schemes to solve the knematc wave equaton have not been studed. Whether or not these fnte volume schemes have advantages over the commonly used fnte dfference schemes are examned n ths paper. Ths study compares the Godunov-type fnte volume method usng scheme and scheme wth the fnte dfference method usng MacCormack scheme. The paper s organzed as follows: Secton 2 ntroduces the knematc wave equaton and ts analytcal solutons; Secton 3 dscusses the numercal schemes: the MacCormack scheme, the scheme and the scheme; Secton 4 shows the results of typcal test cases. Fnally, several concludng remarks are gven n Secton Governng equatons The one-dmensonal knematc wave equaton for flows over a slope s gven by (Eagleson, 97; Lghthll and @x ¼ where h s the depth of flow; q s the dscharge per unt wdth; = f s the ran excess; s the ntensty of ranfall; f s the nfltraton rate; t s the tme; x s the downslope dstance. For the overland flow, the dscharge q s defned as: q ¼ ah m where m s the exponental, and a s the coeffcent. For fully turbulent flow, the coeffcents are gven by Ponce (99): a ¼ p ffffffffff 5 S ; m ¼ ð3þ n 3 where n s the Mannng s roughness coeffcent; S s the bottom slope. It s obvous that the flux functon q(h) s a convex functon (Jacovks and Tabak, 996; Toro, 29) because the second order dervatve s postve: d 2 q dh ¼ amðm 2 Þhm 2 > ; for h > ð4þ The analytcal soluton for Eq. () has been found by (Eagleson, 97) n whch the outflow hydrograph s a functon of ranfall ntensty and the tme of concentraton. ðþ ð2þ
3 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Numercal schemes Snce the knematc wave equaton s a nonlnear hyperbolc partal dfferental equaton, dfferent numercal schemes exhbt dfferent amounts of numercal dffuson and dsperson dependng on the nature of schemes. Numercal dffuson often presents tself as the attenuaton of the knematc wave, whle numercal dsperson s responsble for oscllatons or negatve outflows near large surface gradents. To provde comparsons wth classcal fnte dfference schemes, a study of two hgh-resoluton Godunov-type fnte volume schemes, the fnte volume scheme and the fnte volume scheme, s presented n ths paper. All the selected schemes are explct, but dffer n the order of accuracy. The scheme s a second-order scheme whle the scheme s ffth-order. 3.. MacCormack fnte dfference scheme The MacCormack scheme (MacCormack, 23) s a commonly used fnte dfference scheme to solve hyperbolc PDEs. Ths twostep scheme s second-order accurate n both tme and space. Compared to the frst-order scheme, the MacCormack scheme does not ntroduce numercal dffuson n the soluton. However, the numercal dsperson can be ntroduced n the regon of large surface gradents. The MacCormack scheme s also a varaton of the two-step Lax Wendroff scheme. It ncludes two steps: a predctor step followed by a corrector step. The predctor step uses forward dfference approxmatons whle the corrector step uses backward dfference approxmatons for spatal dervatves. The order of dfferencng can be reversed from tme step to tme step (.e., forward/ backward dfferencng followed by backward/forward dfferencng). The tme step used n the predctor step s Dt n contrast to Dt/2 used n the corrector step. The dscretzed knematc wave equaton usng the MacCormack scheme s below: h nþ Predctor step: ¼ h n Dt Dx ðqn þ qn Þþ Dt Corrector step: h nþ ¼ 2 hn þ h nþ 2 Dt Dx qnþ q nþ þ Dt ð5aþ ð5bþ where subscrpt s the spatal ndex; superscrpt n, n þ, and n + are the temporal ndces; Dt s tme step; Dx s the space step. The classcal MacCormack scheme creates spurous oscllatons at the fronts of shock and rarefacton waves (Garca-Navarro et al., 992; Macchone and Morell, 23). A remedy s to smooth out the oscllatons usng the dsspatve nterface (DI) after the corrector step (Abbott and Mnns, 99). Mathematcal formulaton for the DI s as follows: >< >: h nþ; h nþ; h nþ; ¼ hh nþ þ ¼ hh nþ þ ¼ð hþh nþ þð 2hÞhnþ þ hh nþ ; for nternal nodes þð hþhnþ ; for left boundary node þ hh nþ ; for rght boundary node where h nþ; s the averaged soluton; h s the coeffcent for the dsspatve nterface fnte volume scheme The scheme was ntroduced by van Leer (979). It s the frst second-order TVD Godunov-type fnte volume scheme. MUS- CL uses pecewse lnear approxmaton to reconstruct the depths at the nterfaces of cells: h þ=2 ¼ h þ rh r þ=2 ð6þ ð7þ where h s flow depth at the center of cell ; h +/2 s the reconstructed depth at the nterface + /2 of cell; rh s the lmted depth gradent of cell ; r +/2 s the dstance from the cell center to the nterface + /2. The lmted gradents can be calculated by several slope lmters (e.g., Mnmod, Superbee, etc.). The van Leer lmter (Causon et al., 2; van Leer, 974) s used n the study: rh ¼ rh Ljrh R jþjrh L jrh R jrh L jþjrh R j where rh L ¼ h h Dx ðþ and rh R ¼ h þ h. The local Lax Fredrchs (LLF) Dx method (LeVeque, 22) s used to calculate the fluxes across the cell nterfaces: Q =2 ¼ 2 ½q þ q a =2 ðh h ÞŠ a =2 ¼ maxðjc j; jc jþ ð9þ ðþ where c s the celerty and defned n Eq. (9). Based on the calculated numercal fluxes across the nterfaces, the soluton can be updated by a tme marchng method. Instead of usng the predctor corrector method of -Hancock scheme, the second-order TVD Runge Kutta method (Gottleb and Shu, 99) s employed n the study to dscretze the tme dervatve. The second-order TVD Runge Kutta method (RK2) can be wrtten as: h ðþ ¼ h n þ Dt L ðh n Þ h nþ ¼ 2 hn þ 2 hðþ þ Dt 2 LðhðÞ Þ where the operator L(h) s defned as: ðþ LðhÞ ¼ Dx X QðhÞþ ð2þ For a TVD Runge Kutta method, t s guaranteed that each ntermedate soluton also satsfes the TVD crtera to avod spurous oscllatons n the soluton. Because of ths, t s a better choce to ncorporate a TVD Runge Kutta method for solvng hyperbolc equaton fnte volume scheme The frst scheme was provded by Lu et al. (994). Jang and Shu (996) presented a general framework to construct arbtrary hgh order schemes. Detals about schemes can be found n Shu (999, 29). Instead of usng the pecewse lnear reconstructon procedure n the scheme, the fnte volume scheme uses the reconstructon procedure to obtan an approxmated functon value at the cell nterface. The reconstructon procedure conssts of four steps: Step. Calculate the smoothness ndcators: >< >: b ¼ 3 2 ðh 2 2h þ h Þ 2 þ 4 ðh 2 4h þ 3h Þ 2 b 2 ¼ 3 2 ðh 2h þ h þ Þ 2 þ 4 ðh h þ Þ 2 b 3 ¼ 3 2 ðh 2h þ þ h þ2 Þ 2 þ 4 ð3h 4h þ þ h þ2 Þ 2 ð3þ Step 2. Calculate the thrd-order approxmatons at cell nterfaces: h ðþ þ=2 ¼ h h 6 þ h 6 >< h ð2þ þ=2 ¼ h 6 þ 5 h 6 þ h 3 þ ð4þ >: h ð3þ þ=2 ¼ h 3 þ 5 h 6 þ h 6 þ2 Step 3. Calculate the nonlnear weghts: ~w j w j ¼ ; wth ~w j ¼ ; j ¼ ; 2; 3 ð5þ ~w þ ~w 2 þ ~w 2 3 ðe þ b j Þ c j
4 26 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) where e ¼ 6 for actual calculatons; c j s the lnear weghts, and s gven by: c ¼ 6 ; c 2 ¼ 5 ; c 3 ¼ 5 6 ð6þ Step 4. Calculate the ffth-order approxmaton as a convex combnaton of the three thrd-order approxmatons: h þ=2 ¼ w h ðþ þ=2 þ w h ð2þ þ=2 þ w h ð3þ þ=2 ð7þ The dervatve n tme s dscretzed by the thrd-order TVD Runge Kutta method (Gottleb and Shu, 99), whch s a threestep method: h ðþ ¼ h n þ Dt Lðh n Þ h ð2þ ¼ 3 4 hn þ 4 hðþ þ Dt 4 LðhðÞ Þ h nþ ¼ 3 hn þ 2 3 hð2þ þ 2 Dt 3 Lðhð2Þ Þ where the operator L(h) s gven by Eq. (2). 4. Results ðþ Ths secton evaluates the selected schemes n a varety of flow condtons because a gven scheme does well for a test case does not guarantee t wll do well for another test case. In order to evaluate the performances of these schemes, t s necessary to choose test cases wth exact solutons. For ths reason, the one-dmensonal propagatons of a shock wave and a rarefacton wave are selected. Then the behavors of these schemes are tested by a wave steepenng example. In addton, two synthetc ranfall-runoff cases and one overland flow over a vared slope are selected to demonstrate the applcablty of those schemes. Fnally, a real expermental case s smulated to test the performance of these schemes. The followng abbrevatons are used n the fgures: MUS- CL stands for the results by usng the scheme; for the scheme; for the MacCormack scheme; for the MacCormack scheme wth a dsspatve nterface. h =.25 s the coeffcent for the dsspatve nterface used n all the test cases. 4.. Shock wave The shock wave and rarefacton wave are the ntrnsc features of the hyperbolc equatons, and thus the knematc wave equaton. The celerty of knematc wave s defned as (Lghthll and Whtham, 955): c ¼ amh m ¼ mv ð9þ where v s flow velocty. Consder the followng ntal-value problem for the knematc wave equaton: hðx; Þ ¼ h L; f x < x ð2þ h R ; f x > x If we assume h L > h R, we wll have c L > c R. Ths means that a shock wave wll arse. The exact shock wave soluton for the knematc wave equaton s: hðx; tþ ¼ h L; forðx x Þ=t < S ð2þ h R ; forðx x Þ=t > S where S s the shock wave speed. For the knematc wave equaton, accordng to the Rankne-Hugonot jump condton (Toro, 29), S s equal to: S ¼ a hm R hm L h R h L ð22þ In the study, the followng parameters are used: the channel length L =. m; the bed slope S =.6; the Mannng s coeffcent n =.25 s/m /3 ; the smulaton tme t = 3.5 s; the tme step Dt =. s; the grd spacng Dx =. m. The ntal condton s: h L =. m and h R =.5 m. The boundary condtons are zero-depth gradent at both the entrance and the ext of the channel. The flow depth profles calculated by the numercal schemes and the exact soluton are plotted n Fg.. The results show that the numercal dffuson for flow depth near the front of shock wave s small for the and schemes. But the results of flow depth from the MacCormack scheme are oscllatory. The wggled surface oscllatons are largest at the dscontnuous wave front, and the solutons eventually dverge. No oscllatons are present usng the MacCormack scheme wth the dsspatve nterface. But, the results are not accurate because the dsspatve nterface has smoothed out the sharp front of the shock wave. To determne f a numercal scheme wll ntroduce oscllatons n the soluton, the total varaton defned n Eq. (23) can be used. TVðhÞ ¼ X ¼ jh h j ð23þ where TV(h) s the total varaton of flow depth, s the ndex of computatonal cells. If the total varaton s ncreasng wth tme, the numercal scheme wll ntroduce oscllatons; otherwse, t preserves the monotoncty of ntal monotone functons (LeVeque, 22). To avod oscllatory solutons, the total varaton needs to decrease or reman constant wth tme. The total varatons of flow depth from those schemes are shown n Fg. 2, whch ndcate that the total varaton ncreases wth tme usng the MacCormack scheme, whle t remans constant for, and schemes. Ths mples that the MacCormack scheme wll ntroduce oscllatons at the presence of shock waves, but, and schemes can preserve the functonal propertes of the shock wave. Therefore, the solutons usng the MacCormack scheme are.2..6 EXACT Channel Dstance (m) EXACT Channel Dstance (m) Fg.. Results of shock wave test case: global vew; local vew.
5 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Total Varaton (m) EXACT.5.4 oscllatory, and those usng the scheme are naccurate for capturng the shock wave Rarefacton wave Reconsder the ntal-value problem descrbed n Eq. (2), and assume h L < h R, we have c L < c R. Ths tme, nstead of generatng a shock wave, a rarefacton wave s generated near the dscontnuty snce the celerty at the head of the dscontnuty s greater than that at the tal and, consequently, the dscontnuty contnually expands as t propagates. For the knematc wave equaton, the exact soluton of rarefacton wave s: >< h ¼ >: Tme (s) Fg. 2. Calculated total varatons for the shock wave test case. h L ; for x x t < c L h L þð x x c t L Þ h R h L c R c L ; for c L 6 x x 6 c t R h R ; for x x t > c R ð24þ The parameters used n the rarefacton test case are summarzed as: the channel length L =. m; the bed slope S =.6; the Mannng s coeffcent n =.25 s/m /3 ; the smulaton tme t = 3. s; the tme step Dt =. s; and the grd spacng Dx =. m. The ntal condton s: h L =.5 m and h R =. m. The boundary condtons are zero-depth gradent at both ends of the channel. The calculated flow depth profles are plotted n Fg. 3. The results are smlar to the shock wave case: the MacCormack scheme generated notable oscllatons at the tal of the rarefacton wave. Taylor et al. (972) tested the MacCormack scheme for solvng the Burgers equaton and found that the MacCormack scheme s unstable for rarefacton wave under some condtons. Ths study fnds that the same phenomenon occurred to the knematc wave equaton. Snce the oscllatons occur at the tal of the rarefacton wave, where flow depth s smaller, unrealstc negatve flow depths can be nduced at the tals of rarefacton waves. Ths can also be proved by calculatng the total varaton of flow depth as shown n Fg. 4. The total varatons from, and schemes are constant. But the values from the MacCormack scheme ncrease suddenly to a peak value and then gradually decrease to a constant value (about.9) greater than the ntal total varaton (=.5). Although the total varaton shows a decreasng trend, the scheme cannot keep the total varaton not to ncrease wth tme for the entre smulaton tme. Therefore, oscllatons occur at the dscontnuous front, the tal of rarefacton wave. Therefore, for the rarefacton wave, the MacCormack scheme s not as stable as the, or scheme. Both the shock and rarefacton wave test cases suggest that the MacCormack scheme s not a Total Varaton Dmnshng (TVD) scheme, and the oscllatons wll be generated near the dscontnuous wave Channel Dstance (m) EXACT fronts. Although the scheme has elmnated spurous oscllatons n the soluton, t also smears out the dscontnuous wave fronts by ntroducng redundant numercal dffuson to the soluton. Therefore, the scheme and schemes should be used to smulate unsteady flows wth shock and rarefacton waves Wave steepenng Channel Dstance (m) Fg. 3. Results of rarefacton wave test case: global vew; local vew. Total Varaton (m) Tme (s) Fg. 4. Calculated total varaton for the rarefacton wave test case. One of the promnent features of the hyperbolc equaton s that dscontnutes wll be generated even the ntal water surface s smooth. So t s mportant for a scheme to preserve the stablty
6 2 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) and sharpness of dscontnuous fronts n a smulaton. Ths test case s to test ths behavor of these schemes, especally the ablty of ant-dffuson. Here are the parameters n the test case: the channel length L =. m; the bed slope S =.6; the Mannng s coeffcent n =.25 s/m /3 ; the smulaton tme t = 3.5 s; the tme step Dt =. s; and the grd spacng Dx =. m. The ntal condton s: h h ¼ max :5; : ðx :Þ 2 ð25þ Ths ntal condton wll create a parabolc perturbaton n the channel as shown n Fg. 5. The boundary condtons are perodc flow depths at both ends of the channel. The parabolc perturbaton wll produce a wave steepenng effect n the doman (Henderson, 966; Toro, 29). The reason s that the celerty of the perturbaton s faster than the ambent flud. The head of the perturbaton s a compressve regon whle the tal of the perturbaton s an expansve regon. Ths perturbaton s a combnaton of a shock wave (head) and a rarefacton wave (tal). Ths scenaro s detrmental to the MacCormack scheme because the oscllatons occurred behnd the shock wave wll drown out the parabolc perturbaton mmedately (Fg. 6). Consequently, the MacCormack scheme cannot converge to a stable soluton, but amplfed oscllatons. Ths further proves that the MacCormack scheme s an unstable scheme for flows of dscontnuous waves. The results of, and schemes are plotted n Fg. 7. After 3.5 s, the fronts of the perturbaton reman sharp usng the and schemes, whereas the results from the scheme cannot preserve the sharp front because of the numercal dffuson Unform ranfall-runoff overland flow The unform ranfall-runoff overland flow can be solved by the method of characterstcs (Eagleson, 97; Kazezyılmaz-Alhan et al., 25). At a gven ranfall excess n a specfed duraton, the outflow hydrograph (q = q L, x = L) can be solved analytcally. Assumng a constant ranfall excess ( ) and an ntal zero flow depth, h ¼ ðt ¼ and 6 x 6 LÞ and the boundary condton at the channel entrance, h ¼ ðt > and x ¼ Þ ð26þ ð27þ The solutons of two possble outflow hydrographs are summarzed below: Case (t r P t c ): Channel Dstance (m) Fg. 6. Unstable result of the wave steepenng test case Channel Dstance (m) Fg. 7. Stable results of the wave steepenng test case. >< h L ¼ t; for t 6 t c q L ¼ ah m L ;to solve h L use : h L ¼ t c ; for t c < t 6 t r >: L ¼ ah m L ½h L þ mðt t rþš; for t > t r ð2þ Case 2 (t r < t c ): >< h L ¼ t; for t 6 t r q L ¼ ah m L ;to solve h L use : h L ¼ t r ; for t r < t 6 t p >: L ¼ ah m L ½h L þ mðt t rþš; for t > t p ð29þ where t r s the duraton of ranfall, and t c s the tme of concentraton: 2 t c ¼ L m a! =m ð3þ Channel Dstance (m) Fg. 5. Intal condton of wave steepenng case. and t p s defned as: t p ¼ t r þ t c t r m ; t c ¼ L ah m Lr ; h Lr ¼ t r ð3þ In Kazezyılmaz-Alhan et al. (25) s hypothetcal experment, the duraton of ranfall s longer than the tme of concentraton (t r P t c ). The expermental parameters are summarzed here: the channel length L = 2. m; the bed slope S =.6; the Mannng s coeffcent n =.25 s/m /3 ; the duraton of ranfall t r =.5 h; the ranfall excess = 5. mm/h; smulaton tme t = h; the tme step Dt =. s; and the grd spacng Dx =.3 m. The ntal condton s h =. m. The boundary condtons are zero-depth gradent at the outlet and zero depth at the nlet.
7 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Fg. plots the outflow hydrographs calculated by the numercal schemes and the exact soluton calculated by Eq. (2). Snce there s no dscontnuty/perturbaton n the doman, the results obtaned by the MacCormack, and schemes are very close to the exact soluton. The results usng the scheme showed remarkable numercal dffusons than those usng other schemes. The comparsons of the peak of the hydrographs show that the scheme yelds the closest results to the exact soluton wthout any oscllaton. When the outflow dscharge reaches the peak, flow n the channel reaches steady state. To check f the mass s conserved n the smulaton, the dfferences between the total outflow volume and the total ranfall volume are calculated, and then non-dmensonalzed by the total ranfall volume. The percentage flow dfferences relatve to the total ranfall are.65% for the MacCormack scheme,.59% for the scheme,.35% for the scheme, and.3% for the scheme, respectvely. All schemes preserve mass conservaton very well at the absence of dscontnuous waves. The scheme and scheme are slghtly more conservatve than the MacCormack scheme Steady non-unform ranfall-runoff overland flow To compare the stablty of those schemes, a steady non-unform ranfall-runoff case, Case I n Moramarco and Sngh (22), was smulated. The smulated plane s m long wth a slope of.5 3, and the Mannng s roughness coeffcent s.3 s/m /3. The followng steady non-unform ranfall excess s used n ths case (Fg. 9): ð 5 4:43; x 6 5 m m=sþ ¼ ð32þ :; x > 5 m Ranfall Intensty (m/s) 6 x Channel Dstance (m) Fg. 9. Non-unform ranfall dstrbuton. The smulated results are showed n Fgs. and. The depth of runoff generated by the non-unform ranfall s not smooth and a shock wave s generated near the front of runoff. Oscllatons occurred n the soluton of the MacCormack scheme. Solutons by the, and schemes are stable, and free of oscllatons. But, only the results usng the and schemes can capture the sharp front of shock wave generated by the runoff. The computaton s run wth a DELL XPS notebook (Intel 5 M GHz CPU wth 4 GB memory). The programs are developed usng Matlab 2b runnng on Mcrosoft Wndows 7. OS. For ths smulaton, the CPU tmes usng the MacCormack,,, and schemes are 3.4 s, 3.7 s, 4. s, and 27.3 s, Dscharge (m 3 /s) Dscharge (m 3 /s) EXACT Tme (s) EXACT x Channel Dstance (m) x Tme (s) Channel Dstance (m) Fg.. Hydrograph of the unform ranfall-runoff case: global vew; local vew. Fg.. Results of non-unform ranfall-runoff test case: global vew; local vew.
8 3 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Dscharge (m 3 /s) Flow Dscharge (m 2 /mn) Dscharge (m 3 /s) Tme (s) Tme (s) Flow Dscharge (m 2 /mn) Tme (s) Tme (s) Fg.. Hydrographs of non-unform ranfall-runoff test case: global vew; local vew. Fg. 2. Hydrographs of the non-unform overland slope test case: global vew; local vew. respectvely. As expected, the MacCormack scheme runs the fastest, whle the scheme s two tmes slower than the scheme Ranfall-runoff over non-unform overland slope Snce overland flow always travels n channels of varyng slopes, ths requres the numercal scheme for solvng the knematc wave equaton to be vald for changng bed slopes. Pror test cases have shown that the MacCormack scheme s unstable for flows wth dscontnuous waves (e.g., shock, rarefacton). Therefore, the MacCormack scheme s not applcable to non-unform ranfall on a unform slope. To verfy ts applcablty on a changng slope, ths study hypotheszes an experment of ranfall and overland flow n a watershed havng the same expermental condtons as n Kazezyılmaz-Alhan et al. (25) except the upstream half of the channel has a slope dfferent from the downstream half (S =.6). Ths study adopted S =.6 as the bed slope for the upstream half of the channel. The smulated water surfaces are shown n Fg. 2. As the slope n the upstream half channel ncreases to.6, ten tmes of the downstream slope, the MacCormack scheme generated severe oscllatons. Ths test verfes that the MacCormack scheme s not only unstable for shock and rarefacton waves but also for flows on rapdly varyng slopes Overland flow experment (Schreber, 97) Ths laboratory experment was conducted n the Eroson Laboratory, Palouse Conservaton Feld Staton, USDA-ARS-SWC (Schreber, 97). The experment was conducted on a sloped plane to mmc one dmensonal overland flow. Durng the experment, the ranfall s generated by an array of nozzles. The runoff from the plane was draned nto a tank and the dscharge was measured by a stran gauge. The length of the plane s 4. m. The plane slope s.465, and the Mannng s roughness coeffcent s.25 s/m /3. The ntensty of ranfall s 27 mm/h wth duraton of 4. mn. For the model, the smulaton tme s. mn. The upstream boundary condton s zero flux, the same as the sold wall boundary condton. The downstream boundary condton s defned by the characterstc equaton (Eq. 9). The smulated hydrographs by usng the above mentoned schemes are compared wth the measured ones n Fg. 3. Fg. 3 shows the smulated results are nearly the same usng dfferent schemes. A detaled plot n Fg. 3 shows the MacCormack scheme over-predcted the dscharge. However, none of the smulated results matched the expermental observatons. Ths attrbutes to the assumpton of knematc wave equaton, whch has a smaller celerty than the dynamc wave, and unable to capture the shock wave front. Therefore, the proposed and solver for the knematc wave equaton s vald only for overland flow when the knematc wave assumpton s vald. To evaluate the accuracy of those numercal schemes, the rootmean-square error (RMSE) of three test cases,.e. the shock wave case, the rarefacton wave case, and the ranfall-runoff case, are calculated by comparng to the exact solutons n Table. The RMSE showed that the MacCormack scheme wth the dsspatve nterface has the maxmum error, although ts results are free of oscllaton. Ths mples that the dsspatve nterface brought n consderable numercal dffusons to the solutons. On the other hand, the and the scheme are not only free of oscllaton, but also yelded more accurate results than the MacCormack scheme. In summary, the testng cases lsted n Sectons demonstrated the lmtatons of the MacCormack scheme for solvng the
9 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Dscharge (m 3 /s) x Measured scheme, the oscllatons serously affect the stablty of the soluton, and make the soluton dverge. The dsspatve nterface can elmnate spurous oscllatons seen n the MacCormack scheme, but t also brngs n redundant numercal dffusons to soluton. On the other hand, the and schemes consstently perform better than the MacCormack scheme. There s no oscllaton n the solutons usng the scheme or the scheme. The hgh-order scheme shows the best resoluton power n all test cases. Although the computatonal costs are hgher than that of the MacCormack scheme, to ensure numercal stablty, the and schemes are the stable and accurate technques for solvng the knematc wave equaton. Dscharge (m 3 /s) Tme (mn) x -5 knematc wave equaton and proved the valdty of hgh resoluton schemes, the and schemes. Although the MacCormack scheme has been wdely used to solve the knematc wave equaton, t should be avoded at the presence of shock/rarefacton waves, or flow over rapdly varyng sloped land surfaces. The dsspatve nterface for the MacCormack scheme resulted n excessve numercal dffusons, and s not accurate for solvng unsteady flows wth dscontnuous waves. Instead, the and the scheme should be used n spte of the slght ncrease of computaton cost. 5. Conclusons Tme (mn) Fg. 3. Hydrographs of the Schreber (97) overland flow experment: global vew; local vew. Table Calculated root-mean-square errors (RMSE). Test case Shock wave Rarefacton wave Ranfall-runoff In ths study, the applcablty of four numercal schemes (Mac- Cormack scheme, MacCormack scheme wth the dsspatve nterface, scheme and scheme) to solve the knematc wave equaton s nvestgated. The schemes are examned by usng several test cases: shock wave, rarefacton wave, wave steepenng, overland flow from unform or non-unform ranfall, overland flow over varyng bed slopes. The results show that oscllatons appeared n the solutons usng the MacCormack scheme when the shock and rarefacton waves are present, or bed slopes changng rapdly. Snce the MacCormack scheme s not a TVD Acknowledgements The authors are grateful for research fundng provded by Unted States Natonal Scence Foundaton Award EAR to the Unversty of Arzona. The authors wsh to acknowledge n partcular Chunyan Gao for valuable dscussons on the knematc wave equaton. The authors also would lke to thank the revewers and the edtors for ther constructve comments. References Abbott, M.B., Mnns, A.W., 99. Computatonal Hydraulcs. Ashgate, Aldershot, Brookfeld, USA, vol. xv, p Causon, D.M., Ingram, D.M., Mngham, C.G., Yang, G., Pearson, R.V., 2. Calculaton of shallow water flows usng a Cartesan cut cell approach. Adv. Water Resour. 23 (5), Eagleson, P.S., 97. Dynamc Hydrology. McGraw-Hll, New York, vol. xv, p Garca-Navarro, P., Alcrudo, F., Savron, J.M., D Open-channel flow smulaton usng Tvd-Mccormack scheme. J. Hydraul. Eng. ASCE (), Gottleb, S., Shu, C.W., 99. Total varaton dmnshng Runge Kutta schemes. Math Comput. 67 (22), Harten, A., 93. Hgh-resoluton schemes for hyperbolc conservaton-laws. J. Comput. Phys. 49 (3), Henderson, F.M., 966. Open Channel Flow. Macmllan seres n cvl engneerng, Macmllan, New York, vol. xx, p Jacovks, P.M., Tabak, E.G., 996. A knematc wave model for rvers wth flood plans and other rregular geometres. Math Comput. Model 24 (), 2. Jang, G.S., Shu, C.W., 996. Effcent mplementaton of weghted ENO schemes. J. Comput. Phys. 26 (), Kazezyılmaz-Alhan, C.M., Medna Jr, M.A., Rao, P., 25. On numercal modelng of overland flow. Appl. Math. Comput. 66 (3), Kbler, D.F., Woolhser, D.A., 97. The knematc cascade as a hydrologc model. Colorado State Unversty. LeVeque, R.J., 22. Fnte volume methods for hyperbolc problems. Cambrdge texts n appled mathematcs. Cambrdge Unversty Press, Cambrdge; New York, xx, p. 55. Lghthll, M.J., Whtham, G.B., 955. On knematc waves. I. Flood movement n long rvers. Proc. R. Soc. Lond. A 229 (7), Lu, X.-D., Osher, S., Chan, T., 994. Weghted essentally non-oscllatory schemes. J. Comput. Phys. 5 (), MacCormack, R.W., 23. The effect of vscosty n hypervelocty mpact craterng (Reprnted from AlAA Paper , 969). J. Spacecraft Rockets 4 (5), Macchone, F., Morell, M.A., 23. Practcal aspects n comparng shock-capturng schemes for dam break problems. J Hydraul. Eng. ASCE 29 (3), Moramarco, T., Sngh, V.P., 22. Accuracy of knematc wave and dffuson wave for spatal varyng ranfall excess over a plane. Hydrol. Process. 6 (7), Borah, D.K., Prasad, S.N., Alonso, C.V., 9. Knematc wave routng ncorporatng shock fttng. Water Resour. Res. 6 (3), Ponce, V.M., 99. Engneerng hydrology: prncples and practces. Prentce Hall, Englewood Clffs, N.J., vol. xv, p. 64. Ponce, V.M., 99. The knematc wave controversy. J. Hydraul. Eng. ASCE 7 (4), Schreber, D.L., 97. Overland flow smulaton by a nonlnear dstrbuted parameter model. Thess Thess. Washngton State Unversty., vol. xv, 25 leaves. pp. Shu, C.W., 999. Hgh order ENO and schemes for computatonal flud dynamcs. In: Deconnck (Ed.), Hgh-order methods for computatonal physcs. Sprnger, pp Shu, C.W., 29. Hgh order weghted essentally nonoscllatory schemes for convecton domnated problems. SIAM Rev. 5 (), Sngh, V.P., 2. Knematc wave modellng n water resources: a hstorcal perspectve. Hydrol. Process. 5 (4), Taylor, T.D., Ndefo, E., Masson, B.S., 972. A study of numercal methods for solvng vscous and nvscd flow problems. J. Comput. Phys. 9 (), 99 9.
10 32 C. Yu, J.G. Duan / Journal of Hydrology 59 (24) Toro, E.F., 29. Remann Solvers and Numercal Methods for Flud Dynamcs: A Practcal Introducton. New York, Sprnger, Dordrecht, vol. xxv, p van Leer, B., 974. Towards the ultmate conservatve dfference scheme. II. Monotoncty and conservaton combned n a second-order scheme. J. Comput. Phys. 4 (4), van Leer, B., 979. Towards the ultmate conservatve dfference scheme. 5. 2nd-Order sequel to Godunovs method. J. Comput. Phys. 32 (), 36. Vera, J.H.D., 93. Condtons governng the use of approxmatons for the Sant-Venant equatons for shallow surface-water flow. J. Hydrol. 6 ( 4), 43 5.
Two Dimensional Modeling of Dam-Break Flows
Rver Flow 1 - Dttrch, Koll, Aberle & Gesenhaner (eds) - 1 Bundesanstalt für Wasserbau ISBN 978-3-9393--7 wo Dmensonal Modelng of Dam-Break Flows Ahmad ahershams & Mahnaz Ghaen Hessaroeyeh Department of
More informationS1 Note. Basis functions.
S1 Note. Bass functons. Contents Types of bass functons...1 The Fourer bass...2 B-splne bass...3 Power and type I error rates wth dfferent numbers of bass functons...4 Table S1. Smulaton results of type
More informationMathematics 256 a course in differential equations for engineering students
Mathematcs 56 a course n dfferental equatons for engneerng students Chapter 5. More effcent methods of numercal soluton Euler s method s qute neffcent. Because the error s essentally proportonal to the
More informationDiscontinuous Galerkin methods for flow and transport problems in porous media
T COMMUNICATIONS IN NUMERICA METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2; :1 6 [Verson: 2/3/22 v1.] Dscontnuous Galerkn methods for flow and transport problems n porous meda Béatrve Rvère and Mary
More informationSolitary and Traveling Wave Solutions to a Model. of Long Range Diffusion Involving Flux with. Stability Analysis
Internatonal Mathematcal Forum, Vol. 6,, no. 7, 8 Soltary and Travelng Wave Solutons to a Model of Long Range ffuson Involvng Flux wth Stablty Analyss Manar A. Al-Qudah Math epartment, Rabgh Faculty of
More informationS.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION?
S.P.H. : A SOLUTION TO AVOID USING EROSION CRITERION? Célne GALLET ENSICA 1 place Emle Bloun 31056 TOULOUSE CEDEX e-mal :cgallet@ensca.fr Jean Luc LACOME DYNALIS Immeuble AEROPOLE - Bat 1 5, Avenue Albert
More informationNUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS
ARPN Journal of Engneerng and Appled Scences 006-017 Asan Research Publshng Network (ARPN). All rghts reserved. NUMERICAL SOLVING OPTIMAL CONTROL PROBLEMS BY THE METHOD OF VARIATIONS Igor Grgoryev, Svetlana
More informationSimulation of a Ship with Partially Filled Tanks Rolling in Waves by Applying Moving Particle Semi-Implicit Method
Smulaton of a Shp wth Partally Flled Tanks Rollng n Waves by Applyng Movng Partcle Sem-Implct Method Jen-Shang Kouh Department of Engneerng Scence and Ocean Engneerng, Natonal Tawan Unversty, Tape, Tawan,
More informationA HIGH-ORDER SPECTRAL (FINITE) VOLUME METHOD FOR CONSERVATION LAWS ON UNSTRUCTURED GRIDS
AIAA-00-058 A HIGH-ORDER SPECTRAL (FIITE) VOLUME METHOD FOR COSERVATIO LAWS O USTRUCTURED GRIDS Z.J. Wang Department of Mechancal Engneerng Mchgan State Unversty, East Lansng, MI 88 Yen Lu * MS T7B-, ASA
More informationR s s f. m y s. SPH3UW Unit 7.3 Spherical Concave Mirrors Page 1 of 12. Notes
SPH3UW Unt 7.3 Sphercal Concave Mrrors Page 1 of 1 Notes Physcs Tool box Concave Mrror If the reflectng surface takes place on the nner surface of the sphercal shape so that the centre of the mrror bulges
More informationX- Chart Using ANOM Approach
ISSN 1684-8403 Journal of Statstcs Volume 17, 010, pp. 3-3 Abstract X- Chart Usng ANOM Approach Gullapall Chakravarth 1 and Chaluvad Venkateswara Rao Control lmts for ndvdual measurements (X) chart are
More informationAn Optimal Algorithm for Prufer Codes *
J. Software Engneerng & Applcatons, 2009, 2: 111-115 do:10.4236/jsea.2009.22016 Publshed Onlne July 2009 (www.scrp.org/journal/jsea) An Optmal Algorthm for Prufer Codes * Xaodong Wang 1, 2, Le Wang 3,
More informationParallelism for Nested Loops with Non-uniform and Flow Dependences
Parallelsm for Nested Loops wth Non-unform and Flow Dependences Sam-Jn Jeong Dept. of Informaton & Communcaton Engneerng, Cheonan Unversty, 5, Anseo-dong, Cheonan, Chungnam, 330-80, Korea. seong@cheonan.ac.kr
More informationHermite Splines in Lie Groups as Products of Geodesics
Hermte Splnes n Le Groups as Products of Geodescs Ethan Eade Updated May 28, 2017 1 Introducton 1.1 Goal Ths document defnes a curve n the Le group G parametrzed by tme and by structural parameters n the
More informationNORMALE. A modied structured central scheme for. 2D hyperbolic conservation laws. Theodoros KATSAOUNIS. Doron LEVY
E COLE NORMALE SUPERIEURE A moded structured central scheme for 2D hyperbolc conservaton laws Theodoros KATSAOUNIS Doron LEVY LMENS - 98-30 Département de Mathématques et Informatque CNRS URA 762 A moded
More informationDynamic wetting property investigation of AFM tips in micro/nanoscale
Dynamc wettng property nvestgaton of AFM tps n mcro/nanoscale The wettng propertes of AFM probe tps are of concern n AFM tp related force measurement, fabrcaton, and manpulaton technques, such as dp-pen
More information3D vector computer graphics
3D vector computer graphcs Paolo Varagnolo: freelance engneer Padova Aprl 2016 Prvate Practce ----------------------------------- 1. Introducton Vector 3D model representaton n computer graphcs requres
More informationNumerical Investigation on 1D and 2D Embankment Dams Failure Due to Overtopping Flow
Internatonal Journal of Hydraulc Engneerng 2016, 5(1): 9-18 DOI: 10.5923/.he.20160501.02 Numercal Investgaton on 1D and 2D Embankment Dams Falure Due to Overtoppng Flow Omd Saber 1,*, Gerald Zenz 2 1 PhD
More informationA Binarization Algorithm specialized on Document Images and Photos
A Bnarzaton Algorthm specalzed on Document mages and Photos Ergna Kavalleratou Dept. of nformaton and Communcaton Systems Engneerng Unversty of the Aegean kavalleratou@aegean.gr Abstract n ths paper, a
More informationVery simple computational domains can be discretized using boundary-fitted structured meshes (also called grids)
Structured meshes Very smple computatonal domans can be dscretzed usng boundary-ftted structured meshes (also called grds) The grd lnes of a Cartesan mesh are parallel to one another Structured meshes
More informationNAG Fortran Library Chapter Introduction. G10 Smoothing in Statistics
Introducton G10 NAG Fortran Lbrary Chapter Introducton G10 Smoothng n Statstcs Contents 1 Scope of the Chapter... 2 2 Background to the Problems... 2 2.1 Smoothng Methods... 2 2.2 Smoothng Splnes and Regresson
More informationContent Based Image Retrieval Using 2-D Discrete Wavelet with Texture Feature with Different Classifiers
IOSR Journal of Electroncs and Communcaton Engneerng (IOSR-JECE) e-issn: 78-834,p- ISSN: 78-8735.Volume 9, Issue, Ver. IV (Mar - Apr. 04), PP 0-07 Content Based Image Retreval Usng -D Dscrete Wavelet wth
More informationAn Accurate Evaluation of Integrals in Convex and Non convex Polygonal Domain by Twelve Node Quadrilateral Finite Element Method
Internatonal Journal of Computatonal and Appled Mathematcs. ISSN 89-4966 Volume, Number (07), pp. 33-4 Research Inda Publcatons http://www.rpublcaton.com An Accurate Evaluaton of Integrals n Convex and
More informationElectrical analysis of light-weight, triangular weave reflector antennas
Electrcal analyss of lght-weght, trangular weave reflector antennas Knud Pontoppdan TICRA Laederstraede 34 DK-121 Copenhagen K Denmark Emal: kp@tcra.com INTRODUCTION The new lght-weght reflector antenna
More informationSimulation Based Analysis of FAST TCP using OMNET++
Smulaton Based Analyss of FAST TCP usng OMNET++ Umar ul Hassan 04030038@lums.edu.pk Md Term Report CS678 Topcs n Internet Research Sprng, 2006 Introducton Internet traffc s doublng roughly every 3 months
More informationAn inverse problem solution for post-processing of PIV data
An nverse problem soluton for post-processng of PIV data Wt Strycznewcz 1,* 1 Appled Aerodynamcs Laboratory, Insttute of Avaton, Warsaw, Poland *correspondng author: wt.strycznewcz@lot.edu.pl Abstract
More informationA MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS
Proceedngs of the Wnter Smulaton Conference M E Kuhl, N M Steger, F B Armstrong, and J A Jones, eds A MOVING MESH APPROACH FOR SIMULATION BUDGET ALLOCATION ON CONTINUOUS DOMAINS Mark W Brantley Chun-Hung
More informationJournal of Computational Physics
Journal of Computatonal Physcs 8 (9) 447 465 Contents lsts avalable at ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp A shock-capturng methodology based on adaptatve
More informationA Five-Point Subdivision Scheme with Two Parameters and a Four-Point Shape-Preserving Scheme
Mathematcal and Computatonal Applcatons Artcle A Fve-Pont Subdvson Scheme wth Two Parameters and a Four-Pont Shape-Preservng Scheme Jeqng Tan,2, Bo Wang, * and Jun Sh School of Mathematcs, Hefe Unversty
More informationAVO Modeling of Monochromatic Spherical Waves: Comparison to Band-Limited Waves
AVO Modelng of Monochromatc Sphercal Waves: Comparson to Band-Lmted Waves Charles Ursenbach* Unversty of Calgary, Calgary, AB, Canada ursenbach@crewes.org and Arnm Haase Unversty of Calgary, Calgary, AB,
More informationAn improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws
Avalable onlne at www.scencedrect.com Journal of Computatonal Physcs 7 (8) 9 www.elsever.com/locate/jcp An mproved weghted essentally non-oscllatory scheme for hyperbolc conservaton laws Rafael Borges
More informationUser Authentication Based On Behavioral Mouse Dynamics Biometrics
User Authentcaton Based On Behavoral Mouse Dynamcs Bometrcs Chee-Hyung Yoon Danel Donghyun Km Department of Computer Scence Department of Computer Scence Stanford Unversty Stanford Unversty Stanford, CA
More informationReview of approximation techniques
CHAPTER 2 Revew of appromaton technques 2. Introducton Optmzaton problems n engneerng desgn are characterzed by the followng assocated features: the objectve functon and constrants are mplct functons evaluated
More informationLecture 5: Multilayer Perceptrons
Lecture 5: Multlayer Perceptrons Roger Grosse 1 Introducton So far, we ve only talked about lnear models: lnear regresson and lnear bnary classfers. We noted that there are functons that can t be represented
More informationHelsinki University Of Technology, Systems Analysis Laboratory Mat Independent research projects in applied mathematics (3 cr)
Helsnk Unversty Of Technology, Systems Analyss Laboratory Mat-2.08 Independent research projects n appled mathematcs (3 cr) "! #$&% Antt Laukkanen 506 R ajlaukka@cc.hut.f 2 Introducton...3 2 Multattrbute
More informationFEATURE EXTRACTION. Dr. K.Vijayarekha. Associate Dean School of Electrical and Electronics Engineering SASTRA University, Thanjavur
FEATURE EXTRACTION Dr. K.Vjayarekha Assocate Dean School of Electrcal and Electroncs Engneerng SASTRA Unversty, Thanjavur613 41 Jont Intatve of IITs and IISc Funded by MHRD Page 1 of 8 Table of Contents
More informationCluster Analysis of Electrical Behavior
Journal of Computer and Communcatons, 205, 3, 88-93 Publshed Onlne May 205 n ScRes. http://www.scrp.org/ournal/cc http://dx.do.org/0.4236/cc.205.350 Cluster Analyss of Electrcal Behavor Ln Lu Ln Lu, School
More informationFinite Element Analysis of Rubber Sealing Ring Resilience Behavior Qu Jia 1,a, Chen Geng 1,b and Yang Yuwei 2,c
Advanced Materals Research Onlne: 03-06-3 ISSN: 66-8985, Vol. 705, pp 40-44 do:0.408/www.scentfc.net/amr.705.40 03 Trans Tech Publcatons, Swtzerland Fnte Element Analyss of Rubber Sealng Rng Reslence Behavor
More informationFRICTION SLOPE FORMULAE FOR THE TWO-DIMENSIONAL SHALLOW WATER MODEL
Annual Journal of Hydraulc Engneerng, JSCE, Vol.58, 2014, February FRICTION SLOPE FORMULAE FOR THE TWO-DIMENSIONAL SHALLOW WATER MODEL Hdekazu Yoshoka 1, Koch Unam 2 and Masayuk Fuhara 3 1 Member of JSCE,
More informationVISCOELASTIC SIMULATION OF BI-LAYER COEXTRUSION IN A SQUARE DIE: AN ANALYSIS OF VISCOUS ENCAPSULATION
VISCOELASTIC SIMULATION OF BI-LAYER COEXTRUSION IN A SQUARE DIE: AN ANALYSIS OF VISCOUS ENCAPSULATION Mahesh Gupta Mchgan Technologcal Unversty Plastc Flow, LLC Houghton, MI 49931 Hancock, MI 49930 Abstract
More informationJournal of Computational Physics
Journal of Computatonal Physcs 230 (2011) 2540 2561 Contents lsts avalable at ScenceDrect Journal of Computatonal Physcs journal homepage: www.elsever.com/locate/jcp A grd based partcle method for solvng
More informationSolving two-person zero-sum game by Matlab
Appled Mechancs and Materals Onlne: 2011-02-02 ISSN: 1662-7482, Vols. 50-51, pp 262-265 do:10.4028/www.scentfc.net/amm.50-51.262 2011 Trans Tech Publcatons, Swtzerland Solvng two-person zero-sum game by
More informationA gradient smoothing method (GSM) for fluid dynamics problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluds 2008; 58:1101 1133 Publshed onlne 27 March 2008 n Wley InterScence (www.nterscence.wley.com)..1788 A gradent smoothng method
More informationWishing you all a Total Quality New Year!
Total Qualty Management and Sx Sgma Post Graduate Program 214-15 Sesson 4 Vnay Kumar Kalakband Assstant Professor Operatons & Systems Area 1 Wshng you all a Total Qualty New Year! Hope you acheve Sx sgma
More informationSupport Vector Machines
/9/207 MIST.6060 Busness Intellgence and Data Mnng What are Support Vector Machnes? Support Vector Machnes Support Vector Machnes (SVMs) are supervsed learnng technques that analyze data and recognze patterns.
More informationProper Choice of Data Used for the Estimation of Datum Transformation Parameters
Proper Choce of Data Used for the Estmaton of Datum Transformaton Parameters Hakan S. KUTOGLU, Turkey Key words: Coordnate systems; transformaton; estmaton, relablty. SUMMARY Advances n technologes and
More informationAnalysis of Continuous Beams in General
Analyss of Contnuous Beams n General Contnuous beams consdered here are prsmatc, rgdly connected to each beam segment and supported at varous ponts along the beam. onts are selected at ponts of support,
More informationImprovement of Spatial Resolution Using BlockMatching Based Motion Estimation and Frame. Integration
Improvement of Spatal Resoluton Usng BlockMatchng Based Moton Estmaton and Frame Integraton Danya Suga and Takayuk Hamamoto Graduate School of Engneerng, Tokyo Unversty of Scence, 6-3-1, Nuku, Katsuska-ku,
More informationBackpropagation: In Search of Performance Parameters
Bacpropagaton: In Search of Performance Parameters ANIL KUMAR ENUMULAPALLY, LINGGUO BU, and KHOSROW KAIKHAH, Ph.D. Computer Scence Department Texas State Unversty-San Marcos San Marcos, TX-78666 USA ae049@txstate.edu,
More information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(Ths s a sample cover mage for ths ssue. The actual cover s not yet avalable at ths tme.) Ths artcle appeared n a journal publshed by Elsever. The attached copy s furnshed to the author for nternal non-commercal
More informationAn Image Fusion Approach Based on Segmentation Region
Rong Wang, L-Qun Gao, Shu Yang, Yu-Hua Cha, and Yan-Chun Lu An Image Fuson Approach Based On Segmentaton Regon An Image Fuson Approach Based on Segmentaton Regon Rong Wang, L-Qun Gao, Shu Yang 3, Yu-Hua
More informationWavefront Reconstructor
A Dstrbuted Smplex B-Splne Based Wavefront Reconstructor Coen de Vsser and Mchel Verhaegen 14-12-201212 2012 Delft Unversty of Technology Contents Introducton Wavefront reconstructon usng Smplex B-Splnes
More informationSixth-Order Difference Scheme for Sigma Coordinate Ocean Models
064 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 7 Sxth-Order Dfference Scheme for Sgma Coordnate Ocean Models PETER C. CHU AND CHENWU FAN Department of Oceanography, Naval Postgraduate School, Monterey, Calforna
More informationLearning the Kernel Parameters in Kernel Minimum Distance Classifier
Learnng the Kernel Parameters n Kernel Mnmum Dstance Classfer Daoqang Zhang 1,, Songcan Chen and Zh-Hua Zhou 1* 1 Natonal Laboratory for Novel Software Technology Nanjng Unversty, Nanjng 193, Chna Department
More informationA Cartesian grid finite-difference method for 2D incompressible viscous flows in irregular geometries
Journal of Computatonal Physcs 204 (2005) 302 31 www.elsever.com/locate/jcp A Cartesan grd fnte-dfference method for 2D ncompressble vscous flows n rregular geometres E. Sanmguel-Rojas a, J. Ortega-Casanova
More informationType-2 Fuzzy Non-uniform Rational B-spline Model with Type-2 Fuzzy Data
Malaysan Journal of Mathematcal Scences 11(S) Aprl : 35 46 (2017) Specal Issue: The 2nd Internatonal Conference and Workshop on Mathematcal Analyss (ICWOMA 2016) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES
More informationModeling of Airfoil Trailing Edge Flap with Immersed Boundary Method
Downloaded from orbt.dtu.dk on: Sep 27, 2018 Modelng of Arfol Tralng Edge Flap wth Immersed Boundary Method Zhu, We Jun; Shen, Wen Zhong; Sørensen, Jens Nørkær Publshed n: ICOWEOE-2011 Publcaton date:
More informationPositivity-Preserving Well-Balanced Discontinuous Galerkin Methods for the Shallow Water Equations on Unstructured Triangular Meshes
Postvty-Preservng Well-Balanced Dscontnuous Galerkn Methods for the Shallow Water Equatons on Unstructured Trangular Meshes The MIT Faculty has made ths artcle openly avalable. Please share how ths access
More informationPHANTOM VORTICITY IN EULER SOLUTIONS ON HIGHLY STRETCHED GRIDS
ICAS 000 CONGRESS PHANTOM VORTICITY IN EULER SOLUTIONS ON HIGHLY STRETCHED GRIDS S. A. Prnce, D. K. Ludlow, N. Qn Cranfeld College of Aeronautcs, Bedfordshre, UK Currently DERA Bedford, UK Keywords: Phantom
More informationSLAM Summer School 2006 Practical 2: SLAM using Monocular Vision
SLAM Summer School 2006 Practcal 2: SLAM usng Monocular Vson Javer Cvera, Unversty of Zaragoza Andrew J. Davson, Imperal College London J.M.M Montel, Unversty of Zaragoza. josemar@unzar.es, jcvera@unzar.es,
More informationHigh-Boost Mesh Filtering for 3-D Shape Enhancement
Hgh-Boost Mesh Flterng for 3-D Shape Enhancement Hrokazu Yagou Λ Alexander Belyaev y Damng We z Λ y z ; ; Shape Modelng Laboratory, Unversty of Azu, Azu-Wakamatsu 965-8580 Japan y Computer Graphcs Group,
More informationUrbaWind, a Computational Fluid Dynamics tool to predict wind resource in urban area
UrbaWnd, a Computatonal Flud Dynamcs tool to predct wnd resource n urban area Karm FAHSSIS a, Gullaume DUPONT a, Perre LEYRONNAS a a Meteodyn, Nantes, France Presentng Author: Karm.fahsss@meteodyn.com,
More informationAn Iterative Solution Approach to Process Plant Layout using Mixed Integer Optimisation
17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 An Iteratve Soluton Approach to Process Plant Layout usng Mxed
More informationA mathematical programming approach to the analysis, design and scheduling of offshore oilfields
17 th European Symposum on Computer Aded Process Engneerng ESCAPE17 V. Plesu and P.S. Agach (Edtors) 2007 Elsever B.V. All rghts reserved. 1 A mathematcal programmng approach to the analyss, desgn and
More informationDetermining the Optimal Bandwidth Based on Multi-criterion Fusion
Proceedngs of 01 4th Internatonal Conference on Machne Learnng and Computng IPCSIT vol. 5 (01) (01) IACSIT Press, Sngapore Determnng the Optmal Bandwdth Based on Mult-crteron Fuson Ha-L Lang 1+, Xan-Mn
More informationThe Research of Ellipse Parameter Fitting Algorithm of Ultrasonic Imaging Logging in the Casing Hole
Appled Mathematcs, 04, 5, 37-3 Publshed Onlne May 04 n ScRes. http://www.scrp.org/journal/am http://dx.do.org/0.436/am.04.584 The Research of Ellpse Parameter Fttng Algorthm of Ultrasonc Imagng Loggng
More informationGSLM Operations Research II Fall 13/14
GSLM 58 Operatons Research II Fall /4 6. Separable Programmng Consder a general NLP mn f(x) s.t. g j (x) b j j =. m. Defnton 6.. The NLP s a separable program f ts objectve functon and all constrants are
More informationSolutions to Programming Assignment Five Interpolation and Numerical Differentiation
College of Engneerng and Coputer Scence Mechancal Engneerng Departent Mechancal Engneerng 309 Nuercal Analyss of Engneerng Systes Sprng 04 Nuber: 537 Instructor: Larry Caretto Solutons to Prograng Assgnent
More informationA Fast Visual Tracking Algorithm Based on Circle Pixels Matching
A Fast Vsual Trackng Algorthm Based on Crcle Pxels Matchng Zhqang Hou hou_zhq@sohu.com Chongzhao Han czhan@mal.xjtu.edu.cn Ln Zheng Abstract: A fast vsual trackng algorthm based on crcle pxels matchng
More informationUSING GRAPHING SKILLS
Name: BOLOGY: Date: _ Class: USNG GRAPHNG SKLLS NTRODUCTON: Recorded data can be plotted on a graph. A graph s a pctoral representaton of nformaton recorded n a data table. t s used to show a relatonshp
More informationTsinghua University at TAC 2009: Summarizing Multi-documents by Information Distance
Tsnghua Unversty at TAC 2009: Summarzng Mult-documents by Informaton Dstance Chong Long, Mnle Huang, Xaoyan Zhu State Key Laboratory of Intellgent Technology and Systems, Tsnghua Natonal Laboratory for
More informationSteps for Computing the Dissimilarity, Entropy, Herfindahl-Hirschman and. Accessibility (Gravity with Competition) Indices
Steps for Computng the Dssmlarty, Entropy, Herfndahl-Hrschman and Accessblty (Gravty wth Competton) Indces I. Dssmlarty Index Measurement: The followng formula can be used to measure the evenness between
More informationA New Approach For the Ranking of Fuzzy Sets With Different Heights
New pproach For the ankng of Fuzzy Sets Wth Dfferent Heghts Pushpnder Sngh School of Mathematcs Computer pplcatons Thapar Unversty, Patala-7 00 Inda pushpndersnl@gmalcom STCT ankng of fuzzy sets plays
More informationProblem Definitions and Evaluation Criteria for Computational Expensive Optimization
Problem efntons and Evaluaton Crtera for Computatonal Expensve Optmzaton B. Lu 1, Q. Chen and Q. Zhang 3, J. J. Lang 4, P. N. Suganthan, B. Y. Qu 6 1 epartment of Computng, Glyndwr Unversty, UK Faclty
More informationSubspace clustering. Clustering. Fundamental to all clustering techniques is the choice of distance measure between data points;
Subspace clusterng Clusterng Fundamental to all clusterng technques s the choce of dstance measure between data ponts; D q ( ) ( ) 2 x x = x x, j k = 1 k jk Squared Eucldean dstance Assumpton: All features
More informationThe Codesign Challenge
ECE 4530 Codesgn Challenge Fall 2007 Hardware/Software Codesgn The Codesgn Challenge Objectves In the codesgn challenge, your task s to accelerate a gven software reference mplementaton as fast as possble.
More informationy and the total sum of
Lnear regresson Testng for non-lnearty In analytcal chemstry, lnear regresson s commonly used n the constructon of calbraton functons requred for analytcal technques such as gas chromatography, atomc absorpton
More informationOrder of Accuracy Study of Unstructured Grid Finite Volume Upwind Schemes
João Luz F. Azevedo et al. João Luz F. Azevedo joaoluz.azevedo@gmal.com Comando-Geral de Tecnologa Aeroespacal Insttuto de Aeronáutca e Espaço IAE 12228-903 São José dos Campos, SP, Brazl Luís F. Fguera
More informationDesign of Structure Optimization with APDL
Desgn of Structure Optmzaton wth APDL Yanyun School of Cvl Engneerng and Archtecture, East Chna Jaotong Unversty Nanchang 330013 Chna Abstract In ths paper, the desgn process of structure optmzaton wth
More informationHigh resolution 3D Tau-p transform by matching pursuit Weiping Cao* and Warren S. Ross, Shearwater GeoServices
Hgh resoluton 3D Tau-p transform by matchng pursut Wepng Cao* and Warren S. Ross, Shearwater GeoServces Summary The 3D Tau-p transform s of vtal sgnfcance for processng sesmc data acqured wth modern wde
More informationModule 6: FEM for Plates and Shells Lecture 6: Finite Element Analysis of Shell
Module 6: FEM for Plates and Shells Lecture 6: Fnte Element Analyss of Shell 3 6.6. Introducton A shell s a curved surface, whch by vrtue of ther shape can wthstand both membrane and bendng forces. A shell
More information6.854 Advanced Algorithms Petar Maymounkov Problem Set 11 (November 23, 2005) With: Benjamin Rossman, Oren Weimann, and Pouya Kheradpour
6.854 Advanced Algorthms Petar Maymounkov Problem Set 11 (November 23, 2005) Wth: Benjamn Rossman, Oren Wemann, and Pouya Kheradpour Problem 1. We reduce vertex cover to MAX-SAT wth weghts, such that the
More informationModeling of Fillets in Thin-walled Beams Using Shell/Plate and Beam Finite. Elements
Modelng of Fllets n Thn-walled Beams Usng Shell/Plate and Beam Fnte Elements K. He Graduate Research Assstant W.D. Zhu Professor and Correspondng Author Emal: wzhu@umbc.edu; Tel: 410-455-3394 Department
More informationOverview. Basic Setup [9] Motivation and Tasks. Modularization 2008/2/20 IMPROVED COVERAGE CONTROL USING ONLY LOCAL INFORMATION
Overvew 2 IMPROVED COVERAGE CONTROL USING ONLY LOCAL INFORMATION Introducton Mult- Smulator MASIM Theoretcal Work and Smulaton Results Concluson Jay Wagenpfel, Adran Trachte Motvaton and Tasks Basc Setup
More informationAPPLICATION OF A COMPUTATIONALLY EFFICIENT GEOSTATISTICAL APPROACH TO CHARACTERIZING VARIABLY SPACED WATER-TABLE DATA
RFr"W/FZD JAN 2 4 1995 OST control # 1385 John J Q U ~ M Argonne Natonal Laboratory Argonne, L 60439 Tel: 708-252-5357, Fax: 708-252-3 611 APPLCATON OF A COMPUTATONALLY EFFCENT GEOSTATSTCAL APPROACH TO
More informationChapter 9. Model Calibration. John Hourdakis Center for Transportation Studies, U of Mn
Chapter 9 Model Calbraton John Hourdaks Center for Transportaton Studes, U of Mn Hourd00@tc.umn.edu Why Calbrate? Computers Cannot Magcally Replcate Realty! Smulaton Models Are Desgned to be General Drver
More informationREFRACTION. a. To study the refraction of light from plane surfaces. b. To determine the index of refraction for Acrylic and Water.
Purpose Theory REFRACTION a. To study the refracton of lght from plane surfaces. b. To determne the ndex of refracton for Acrylc and Water. When a ray of lght passes from one medum nto another one of dfferent
More informationLearning-Based Top-N Selection Query Evaluation over Relational Databases
Learnng-Based Top-N Selecton Query Evaluaton over Relatonal Databases Lang Zhu *, Wey Meng ** * School of Mathematcs and Computer Scence, Hebe Unversty, Baodng, Hebe 071002, Chna, zhu@mal.hbu.edu.cn **
More informationQuality Improvement Algorithm for Tetrahedral Mesh Based on Optimal Delaunay Triangulation
Intellgent Informaton Management, 013, 5, 191-195 Publshed Onlne November 013 (http://www.scrp.org/journal/m) http://dx.do.org/10.36/m.013.5601 Qualty Improvement Algorthm for Tetrahedral Mesh Based on
More informationUniversität Stuttgart Direkte numerische Simulation von Strömungslärm in komplexen Geometrien
Unverstät Stuttgart Drekte numersche Smulaton von Strömungslärm n komplexen Geometren Claus-Deter Munz Gregor Gassner, Floran Hndenlang, Andreas Brkefeld, Andrea Beck, Marc Staudenmaer, Thomas Bolemann,
More informationAn implicit finite difference method for drainage basin evolution
WATER RESOURCES RESEARCH, VOL. 38, NO. 7, 10.1029/2001WR000721, 2002 An mplct fnte dfference method for dranage basn evoluton Sergo Fagherazz, 1 Alan D. Howard, and Patrca L. Wberg Department of Envronmental
More informationTN348: Openlab Module - Colocalization
TN348: Openlab Module - Colocalzaton Topc The Colocalzaton module provdes the faclty to vsualze and quantfy colocalzaton between pars of mages. The Colocalzaton wndow contans a prevew of the two mages
More information2x x l. Module 3: Element Properties Lecture 4: Lagrange and Serendipity Elements
Module 3: Element Propertes Lecture : Lagrange and Serendpty Elements 5 In last lecture note, the nterpolaton functons are derved on the bass of assumed polynomal from Pascal s trangle for the fled varable.
More informationConstruction of second-order accurate monotone and stable residual distribution schemes for steady problems
Journal of Computatonal Physcs 195 (2004) 474 507 www.elsever.com/locate/cp Constructon of second-order accurate monotone and stable resdual dstrbuton schemes for steady problems Rem Abgrall *,1, Mohamed
More informationBioTechnology. An Indian Journal FULL PAPER. Trade Science Inc.
[Type text] [Type text] [Type text] ISSN : 0974-74 Volume 0 Issue BoTechnology 04 An Indan Journal FULL PAPER BTAIJ 0() 04 [684-689] Revew on Chna s sports ndustry fnancng market based on market -orented
More informationIP mobility support is becoming very important as the
706 IEEE TRANSACTIONS ON COMPUTERS, VOL. 52, NO. 6, JUNE 2003 A New Scheme for Reducng Lnk and Sgnalng Costs n Moble IP Young J. Lee and Ian F. Akyldz, Fellow, IEEE Abstract IP moblty support s provded
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
Ths artcle appeared n a journal publshed by Elsever. The attached copy s furnshed to the author for nternal non-commercal research and educaton use, ncludng for nstructon at the authors nsttuton and sharng
More informationApplication of Residual Distribution (RD) schemes to the geographical part of the Wave Action Equation
Applcaton of Resdual Dstrbuton (RD) schemes to the geographcal part of the Wave Acton Equaton Aron Roland Insttute for Hydraulc and Water Resources Engneerng, Technsche Unverstät Darmstadt, Germany Abstract
More informationMachine Learning 9. week
Machne Learnng 9. week Mappng Concept Radal Bass Functons (RBF) RBF Networks 1 Mappng It s probably the best scenaro for the classfcaton of two dataset s to separate them lnearly. As you see n the below
More informationComputer models of motion: Iterative calculations
Computer models o moton: Iteratve calculatons OBJECTIVES In ths actvty you wll learn how to: Create 3D box objects Update the poston o an object teratvely (repeatedly) to anmate ts moton Update the momentum
More information