• Journal of Korea Water Resources Association
• /
• v.36 no.5
• /
• pp.751-760
• /
• 2003

Error analysis of finite difference equation on the Muskingum-Cunge flood routing method with free time and space weighting factor was carried out. The error analysis shows that the numerical solution of the Muskingum-Cunge method becomes diverged with time when the sum of time weighting factor and space weighting factor is greater than 1.0. Numerical diffusion increases when the sum of time weighting factor and space weighting factor decreases. Numerical diffusion and numerical oscillation increase when the grid resolution is coarse. Numerical experiments and field applications show that the Muskingum-Cunge method with free space weighting factor is more effective for simulating the flood routing with great peak diminution than conventional Muskingum-Cunge method with fixed space weighting factor, 0.5.

• Journal of the Korean Society of Industry Convergence
• /
• v.11 no.1
• /
• pp.27-34
• /
• 2008

This study presents techniques for the estimation of parameters in flood routing method of natural channel.. The Muskingum routing method is the most widely used method of hydrologic stream channel routing. In this paper, Genetic Algorithm and Fletcher-Powell method is applied to determine parameters(K and x) of the Muskingum routing method. The results of the approach shows that Genetic Algorithm method can be one of methods to determine parameters of the Muskingum routing method. Based on the analysis for estimated parameters and the comparison with the results from observed data, the applicability of Genetic Algorithm is verified.

• Journal of Korea Water Resources Association
• /
• v.32 no.4
• /
• pp.417-426
• /
• 1999

A series of numerical experiments is performed to compare the characteristics of outflow hydrograph using linear and nonlinear Muskingum-Cunge methods for two cases: (a) sinusoidal inflow hydrographs and (b) rainfall inputs. The nonlinear method shows the steepening of the rising limb, coupled with a corresponding flattening of the receding limb. The linear method conserves mass exactly. In contrast, the nonlinear method is subject to a gain and a loss of mass. The loss of mass and the subsidence of peak outflow increases with a mild slope, a small baseflow $q_b$ and a large peak inflow to baseflow ratio $q_p/q_b$. A shock wave and associated numerical instability results in the increase of mass for a steep slope and a large $q_p/q_b$ ratio. While the linear method depends on the reference flow per unit-width, the nonlinear method depends on a baseflow and the $q_p/q_b$ ratio. It is found that, unlike for the sinusoidal inflow, the outflow for the rainfall inputs conserves mass fairly exactly in the nonlinear method.

• Journal of Korea Water Resources Association
• /
• v.43 no.12
• /
• pp.1051-1061
• /
• 2010

This study hydrologically re-analysed the Muskingum channel routing method to represent it as a linear combination of the linear channel model considering only the translation and the linear reservoir model considering only the storage effect. The resulting model becomes a kind of instantaneous unit hydrograph, whose parameters are identical to those of the Muskingum model. That is, the outflow occurs after the routing interval ${\Delta}t$ or concentration time $T_c$, and among the total amount of inflow, the x portion is translated by the linear channel model and the remaining (1-x) portion is routed by the linear reservoir model with the storage coefficient ��$K_c$. The application result of both the Muskingum channel routing method and its corresponding instantaneous unit hydrograph to an imaginary channel showed that these two models are basically identical. This result was also assured by the application to the channel flood routing to the Kumnam and Gongju Station for the discharge from the Daechung Dam.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.6 no.1
• /
• pp.1-11
• /
• 2003

In this study, TOPMODEL(TOPography based hydrologic MODEL) was tested linked with Muskingum river routing technique for $581.7km^2$ Anseong-cheon watershed. Linear trend surface interpolation was used to give flow direction for flat areas located in downstream watershed. MDF (multiple flow direction) algorithm was adopted to derive the distribution of ln(a/$tan{\beta}$) values of the model. Because the coarser DEM resolution, the greater information loss, the watershed was divided into subwaterhseds to keep DEM resolution, and the simulation result of the upstream watershed was transferred to downstream watershed by Muskingum techniques. Relative error of the simulated result by 500 m DEM resolution showed 27.2 %. On the other hand, the relative error of the simulated result of 300 m DEM resolution by linked 2 subwatersheds with Muskingum method showed 15.8 %.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.21 no.1
• /
• pp.13-23
• /
• 1979

The predictien of a design flood hydrograph at a particular site on a river may be based on the derivation of a discharge or stage hydrograph at an upstream section, together with a method to route this hydrograph along the rest of the river. In order to limit this investigation to cases where the assumption like uniform rainfall may be reasonably valid, the derivation of unit hydrographs has been limited to catchment with an area less than 500 km2. Consequently, flood routing methods provide a useful tool for the analysis of flooding in all but the smaller catchment, particularly where the shape of the hydrograph as well as the peak value is required. The author, therefore, will introduce here a flood routing method on the open channel with a peak discharge of the catchment area concerned. The importance of being able to route floods accurately is also reflected in the large number of flood routing method which have been developed since the year 1900. There are the modified puls method, Steinberg method, Goodrich method, Ekdahl method, Tatum's mean continuously Equation, wisler-Brater method, Muskingum, chung, and Muskingum-cunge (M-C) method and so on. The author will try to introduce a flood routing method which is revised Muskingum-cunge method. In calculating flood routing by the M-C method, whole variable parameters on the river were assumed to almost uniform values from the upstream to the downstream. In the results, the controlled flood rates at the 40km downstream on the river is appeared to decrease 22m$^3$/sec or 12 percent of the peak flood 170m$^3$/sec.

• Water for future
• /
• v.18 no.3
• /
• pp.265-270
• /
• 1985

The prediction of a design-flood hydrograph at a particular site on a river may be based on the derivation of discharge or stage hydrograph at an upstream section, togeater with a method to route this hydrograph along the rest of river. On the other hand, flood routing methods provide a useful tool for the analysis of flooding in all but the smaller catchment, and these methods are largely stored into hydrological method and hydraulic method. Although the Muskingum Method as a hydrological method ignores dynamic effects on the flood wave, Muskingum-Cunge Method based on hydraulic method is possible to improve the method so that it gives a good approximation to the solution of the linear convective-diffusion equation. This is made on the basis of the finite diffeience equation for the Muskingum Method. In the study, the outflows predicted by Muskingum-Cunge Method are campared with the observed outflows of the Pyung Chang River.

• Water for future
• /
• v.29 no.5
• /
• pp.139-150
• /
• 1996

This paper presents the effect of numerical parameters, such as grid size and grid ratio, on the outflow hydrograph of a unit-width plane in the linear Muskingum-Cunge method. The numerical results depend on Courant number C and cell Reynolds number D, two physically and numerically meaningful parameters. As C approache 1 and D increases, the numerical dispersion-relating oscillations are difficult to occur. The numerical oscillations occur in the front of a propagating wave for C < 1, while smaller oscillations occur behind the wave for C > 1 due to the numerical diffusion effect. For a plane with a small value of characteristic reach length L (e.g., a steep plane), the numerical solution of the Muskingum-Cunge method is similar to that of the kinematic wave method, which shows no wave attenuation. However, for a plane with a large value of L (e.g., a mild plane), the Muskingum-Cunge method leads to the diffusion waves which are essentially independent of the Courant number. Accordingly, the Muskingum-Cunge method will be suited for the routing of the catchment with relatively mild slopes.

• Journal of Korea Water Resources Association
• /
• v.43 no.2
• /
• pp.233-243
• /
• 2010

A method is proposed of estimating Muskingum-Cunge parameters for natural streams using cross-sectional and longitudinal channel geometry and roughness coefficient data. Firstly, for various water-surface levels at a cross section cross-sectional areas and hydraulic radii are calculated. Corresponding discharges are then calculated using Manning's equation. This procedure is repeated for all cross-sections in the reach. Finally, routing parameters are estimated from the calculated cross-sectional area and discharge value pairs by regression analysis. The procedures for estimating Muskingum-Cunge parameters are applied to the South Han River. Flows calculated by Muskingum-Cunge model with estimated parameters showed much better agreement with those by dynamic wave model in peak discharge, time to peak discharge, and normalized RMS errors than those calculated by the HEC-1 Muskingum-Cunge model.

• Journal of Korea Water Resources Association
• /
• v.46 no.5
• /
• pp.449-463
• /
• 2013

This study proposes the method for determining the Muskingum channel routing model parameters based on the assumption of linear system. The proposed method was applied to the Chungju dam basin for the evaluation. Additionally, the rainfall-runoff was repeated for the Yeongchun-Chungju dam reach using seven rainfall events observed. Summarizing the results is as follows. First, the concentration time and storage coefficient of a channel reach formed by the subdivision can be expressed as the difference between the concentration times and storage coefficients of upstream and downstream basins. The storage coefficients of the channel reach estimated is equal to the storage coefficient of the Muskingum channel routing model and the weight factor can be simply estimated using the ratio between the concentration time and storage coefficient. Second, the weight factor of the Muskingum model is in inverse proportion to the Russel coefficient, which is in between 0.4166 and 0.625 when considering the Russel coefficients generally applied. Finally the application to the Yeongchun-Chungju dam reach showed that the proposed method is still valid regardless of the limitations such as the uncertainty of the observed data.