- 51 Views
- Uploaded on
- Presentation posted in: General

2D finite element modeling of bed elevation change in a curved channel

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

2D finite element modeling of bed elevation change in a curved channel

S.-U. Choi, T.B. Kim, & K.D. Min

Yonsei University

Seoul, KOREA

- Most natural streams are sinuous and meandering.
- In a curved channel, the centrifugal force makes the flow structure and sediment transport mechanism extremely complicated.
- To simulate the flow and morphological change in a curved channel, the secondary currents and the gravity effect due to morphological change should be properly considered (Kassem & Chaudhry, 2002).

- 1D Model
- Impossible to account for sediment transport in the transverse direction

- 3D Model
- Still Expensive
- Not readily applicable to many engineering problems
- Turbulence Closure
- Sediment Transport Model
- Boundary Conditions

- Only applicable to the steady flow condition, constant channel width and constant radius curvature
- Koch & Flokstra (1981), Struiksma et al. (1985), Shimizu & Itakura (1989), Yen & Ho (1990) and so on.

- The coordinate transformed, unsteady FDM & FVM
- Kassem & Chaudhry (2002), Duc et al. (2004), Wu (2004)

- The finite element model for bed elevation change in a curved channel has never been proposed!!
- FEM provides greater flexibility in handling spatial domain.

- To Develop a 2D FEM model
- capable of predicting time-dependent morphological change in a curved channel.

- For flow analysis, the shallow water equations are solved by the SU/PG scheme.
- To assess the be elevation change, Exner’s equation is solved by BG scheme.
- For validation, we applied the model to two laboratory experiments.

- Decoupled modeling approach
- Flow equations and Exner’s equations are solved separately.

- Neglecting armoring or grain sorting effects.

- 2D shallow water equations with the effective stress terms
- Eddy viscosity model

- Exner’s equation
- Total Sediment Load

Engelund & Hansen’s formula

- Flow Equations
- Weighted Residual Equations
- 2D SU/PG Method (Ghanem, 1995)

- Exner’s Equation
- Weighted Residual Equation
- BG Method

- Upstream & Downstream BCs
- Sidewall BC

(Akanbi, 1986)

(a) Under a flat (& fixed) bed condition

- The centrifugal force makes higher flow depth, but lower mean velocity, at the outer bank. This generates the secondary flows satisfying the continuity.

- Observed in Experiments and Numerical simulations.

(b) Under a mobile bed condition

- Secondary flows induces sediment erosion & deposition at the outer & inner banks, respectively.

- The flow depth and mean velocity at the inner bank is lower.

- Observed in natural meandering rivers and Experiment by Yen (1967 & 70)

- Gravity effect on a slope
- Angle of bed shear stress due to the secondary flow effect

(Struiksma et al., 1985)

(Rozovskii, 1957)

1. 180º Curved Channel Experiment

Lab. of Fluid Mech. (LFM) in Delft Univ. of Tech.(Sutmuller & Glerum, 1980)

2. 140º Curved Channel Experiment

Delft Hydraulics Lab. (DHL)

(Struiksma, 1983)

- Experimental Conditions

- 1400 elements, 1551 nodes
- Porosity = 0.4
- 10 times extension of width
- Fr = 0.36
- Fixed bed B.C. for upstream & downstream boundaries

10 min.

150 min.

- At the initial stage, the particles are heading for the inner bank. This induces sediment deposition & erosion at the inner and outer banks.
- After for a while, the gravity effect due to changed bed reduces the secondary flow effect.

10 min.

150 min.

- At the initial stage, the flow depth near the outer bank is higher than that near inner bank.
- A similar pattern at 150 min. But, considering deposition & erosion, the water surface elevation across the width is nearly uniform.

10 min.

150 min.

- At the initial stage, the mean velocity near the inner bank is slightly higher.
- Later, we have an opposite situation after bed deformation.

Measured data by

Sutmuller & Glerum (1980)

Simulated Result

- In the numerical simulation, the bed elevation change became negligible after 150 min.
- A good agreement. But the location of max deposition is slightly different. This may be …

- Overall trend is the same.
- Near the inner bank, the amount of sediment deposition is over-predicted.
- Near the outer bank, the amount of sediment erosion is under-predicted.

- Experimental Conditions
(Struiksma, 1983)

- 945 elements, 804 nodes
- Porosity = 0.4
- 10 & 15 times extension of width
- for US & DS, respectively
- Fixed bed B.C. for both US & DS

- Simulated result is after 10 hr.
- Overall trend is the same.
- Especially good agreement in max deposition & erosion.
- The simulated results fluctuate with distance while...

- Spatial fluctuation increases with time.
- This is due to BG scheme applied to Exner’s eq.

- Development of 2D FEM model for bed elevation change
- SU/PG method for shallow water eqs.
- BG method for Exner’s eq.
- Secondary flow effect and gravity effect on sloping bed

- Applications to 2 curved channel experiments
- The model predicts the flow and bed morphology well.
- Specially, the time-evolution of changing bed morphology from the flat bed.
- Sediment deposition & erosion at the inner & outer banks.

- The model predicts the flow and bed morphology well.
- Necessity of introducing the upwind scheme to Exner’s eq.
- Spatial fluctuations in the simulated bed profiles increase with time.
- Weighting is required in the upwind direction along the trajectory of sediment particles.