• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.10 no.4
• /
• pp.174-186
• /
• 1998

A Galerkin's finite element model incorporating infinite elements for modeling of radiation condition at infinity has been developed, which is based on an extended mild-slope equation. To illustrate the validity and applicability of the present model, the example analyses were carried out for a resonance problem in the rectangular harbor of Ippen and Goda (1963) and for wave transformations over circular shoals of Sharp (1968) and Chandrasekera and Cheung (1997). Comparisons with the results obtained by hydraulic experiments and hybrid element method showed that the present model gives very good results in spite of the rapidly varying topography. Numerical experiments were also performed for wave transformations over a circular concave well which may be an alternative to conventional wave barriers.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.15 no.3
• /
• pp.159-166
• /
• 2003

In this study, we develop a finite element model to directly solve the Laplace equation while keeping the same computational efficiency as the models based on the extended mild-slope equation which has been widely used for calculation of wave transformation in shallow water. For this, the computational domain is discretized into finite elements with a single layer in the vertical direction. The velocity potential in the element is then expressed in terms of the potentials at the nodes located at water surface, and the Galerkin method is used to construct the numerical model. A common shape function is adopted in horizontal direction, and the cosine hyperbolic function in vertical direction, which describes the vertical behavior of progressive waves. The model was developed for vertical two-dimensional problems. In order to verify the developed model, it is applied to vertical two-dimensional problems of wave reflection and transmission. It is shown that the present finite element model is comparable to the models based on extended mild-slope equations in both computational efficiency and accuracy.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1741-1745
• /
• 2006

본 연구에서는 하천에 유입되는 오염물질 중 대부분을 차지하는 비보존성 오염물질의 확산거동을 예측하기 위해 2차원 수심평균된 이송분산방정식에 유한요소모형을 적용하였다. 수치모형 구성을 위해 Galerkin법을 이용한 가중잔차법을 사용하였으며, 복잡한 하천경계를 보다 정확히 재현할 수 있도록 삼각 및 사각요소망의 혼용이 가능하도록 하였다. 모의대상 수질인자는 BOD, DO, 질소화합물, 인화합물, 수온, pH 및 대장균군수이며, 이 가운데 BOD와 DO는 상호 쌍을 이루는 방정식을 풀어야 하는 특수한 형태이므로 별도로 취급하였다. 순간주입 및 연속주입에 의한 비보존성 오염물질의 확산거동을 모의하였으며, 시간에 따른 민감도를 분석하기 위해 보존성 오염물질의 확산거동과 비교하였다. 해석해를 이용해 순간주입된 오염물질의 오염운을 구하는 식을 유도하고 해석해와 본 연구에서 개발한 수치모형에 의한 수치해를 비교하였다.

• The Journal of the Korea Contents Association
• /
• v.11 no.1
• /
• pp.404-410
• /
• 2011

For the simulation of one-dimensional unsteady flow, the Taylor-Galerkin finite element method was adopted to the discretization of the Saint Venant equation. The model was applied to the backwater problem in a single channel and the flood routing in dendritic channel networks. The numerical solutions were compared with previously published results of finite difference and finite element methods and good agreement was observed. The model solves the continuity and the momentum equations in a sequential manner and this leads to easy implementation. Since the final system of matrix is tri-diagonal with a few additional entry due to channel junctions, the tri-diagonal matrix solution algorithm can be used with minor modification. So it is fast and economical in terms of memory for storing matrices.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.475-483
• /
• 2008

2D finite element model for analysis of transport of accidentally released pollutants in the flow was developed by SUPG method, and the mass balance of this model was checked though two example problems: line source and point source problem in the straight channel and unidirectional 2D flow field, respectively. All the test cases were simulated with both SUPG and conventional Galerkin method to compare the accuraccy of the numerical mass balance. Test results show that the model with SUPG can adequately conserve the released mass though simulation than the model using Galerkin method, so the developed model verified to be appropriate to solve this accidental mass release problem.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2B
• /
• pp.199-209
• /
• 2010

A finite element model for analyzing surface water flow was developed. Shallow water equation was discretized and solved by Galerkin and Newton-Raphson method. Triangular or rectangular elements can be mixed together to construct meshes. The algebraic equation was solved by frontal method which is very efficient in finite element problem. The developed model was applied to rectangular meandering channel with two bends and transverse velocities and water depth distributions were examined. High velocity was located near the inner bank at the apexes of the bends and velocity distribution was symmetrical about the centerline at the midsection of two bend and super elevation also occurred. Simulation results showed very good agreement with measured data. Another numerical simulation was carried out in mild, steep, adverse and abrupt bottom change slope and channels with weir. 12 water surface profiles of gradually varied flow were correct in terms of hydraulic interpretation.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.3
• /
• pp.133-141
• /
• 2007

A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.201-214
• /
• 1987

하천에서 종, 횡방향이 고려된 확산-이송 현상을 수치기법인 유한요소법을 이용하여 2차원으로 해석하였다. 유한요소법으로는 Galerkin의 가중잔차 방법을 수십에 대해 적분을 취한 연속, 운동량 및 확산-이송방정식에 적용하였고, 선형보간함수와 선형삼각형요소가 이용되었다. 모형의 타당성을 입증하기 위해 단순화된 1차원 수로에서 수차례 검정한 결과 정확해와 거의 일치하는 만족할만한 결과가 도출되었다. 개발된 모형의 실험이 2차원수로에서 행하여져 지류의 유입에 따른 확산-이송현상이 모의되었으며, 실험적용은 개발사업후의 한강본류 9km 구간에 적용되어 탄천과 중량천의 지천 영향을 받는 오염 농도가 2차원적으로 도시되었다.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.439-451
• /
• 1994

A numerical model for harbour oscillation is presented by use of Galerkin finite element method. The governing equation is used by the modified mild slope equation derived from Chen (1986) in which bottom friction is incorporated. Since the existing absorbing boundary condition. however. is shown to be incorrect correct boundary condition and forcing term due to an incident plane wave are rederived. Computation results for a rectangular harbour are shown in comparison with both laboratory data and existing numerical results. After the values of friction factor (f) and reflection coefficient (K$_{r}$) are discussed, the set (K$_{r}$=0, 94, f=0) is found to be best fitted to the laboratory data of the rectangular harbour.