Journal of Korean Society of Coastal and Ocean Engineers (한국해안·해양공학회논문집)

Korean Society of Coastal and Ocean Engineers (한국해안·해양공학회)
권호 표지

Wiggle-free Finite Element Model for extended Boussinesq equations

확장형 Boussinesq FEM model의 수치진동오차 개선

  • Woo, Seung-Buhm (Department of Ocean Science, Inha University) ;
  • Choi, Young-Kwang (Department of Coastal Oceanography and Engineering, GeoSystem Research Corporation) ;
  • Gonzalez-Ondina, Jose M. (School of Civil & Environmental Engineering, Cornell University)
  • 우승범 (인하대학교 자연과학계열 해양과학과) ;
  • 최영광 ((주)지오시스템리서치 연안해양.공학부) ;
  • Published : 2010.02.28
Download PDF
⟨ Previous Next ⟩

Abstract

Subgrid scale stabilization method is applied to Woo and Liu(2004)'s extended Boussinesq FEM numerical model to eliminate the 2dx wiggles. In order to optimize the computational efficiency, Hessian operator is introduced and the matrix of velocity vector is combined to one matrix for solving matrix equations. The mass lumping technique is also applied to the matrix equations of auxiliary variables. The newly developed code is applied to simulate Vincent and Briggs(1989)' wave transformation experiments and the results show that the numerical solution is almost wiggle-free and it matches very well with experimental data. Due to improvement of computational efficiency and wiggle reduction, it is plausible to apply this model to a realistic problem such as harbor oscillation problems.

Woo and Liu (2004)의 확장형 Boussinesq FEM 수치모형에서 한계점으로 지적되었던 수치진동현상과 계산 효율성이 크게 개선되었다. 수치진동을 해결하기 위해 subgrid scale stabilization method를 사용하였고, 계산효율성을 높이기 위해서 Hessian 연산자를 도입하였으며, 유속벡터에 대한 행렬 구성을 하나의 행렬로 구성하였다. 또한 추가변수에 대한 행렬은 mass lumping technique을 사용하여 대각행렬로 구성하였다. Vincent and Briggs(1989)의 파랑 굴절 및 회절에 대한 수치실험 결과 수치진동현상이 확연히 줄어 들은 것을 관찰할 수 있었으며, 수리실험 결과와도 상당히 일치하는 결과를 보였다. 이전 모형에 비해 약 10배의 계산소요시간이 줄어 향후 항만부진동이나 퇴적물 이동과 같은 현실적인 문제에 적용이 가능할 것으로 기대된다.

Keywords

References

  1. Chen, Q., Kirby, J.T., Dalrymple, R.A., Kennedy, A.B. and Chawla, A. (2000). Boussinesq modeling of wave transformation, breaking, and runup. II:2D. Journal of Waterway, Port, Coastal, and Ocean Engineering, 126(1), 48-56. https://doi.org/10.1061/(ASCE)0733-950X(2000)126:1(48)
  2. Codina, R. (2001). A stabilized finite element method for generalized stationary incompressible flows. Computer methods in applied mechanics and engineering, 190, 2681-2706. https://doi.org/10.1016/S0045-7825(00)00260-7
  3. Goodbody, A.M. (1982). Cartesian tensors : With applications to mechanics, fluid mechanics and elasticity. JOHN WILEY &SONS, New York, N.Y.
  4. Johnson, C. (1987). Numerical solution of partial differential equations by the finite element method. Cambridge University Press.
  5. Kennedy, A.B., Chen, Q., Kirby, J.T. and Dalrymple, R.A. (2000). Boussinesq modeling of wave transformation, breaking, and runup. I: 1D. Journal of Waterway, Port, Coastal, and Ocean Engineering, 126(1), 39-47. https://doi.org/10.1061/(ASCE)0733-950X(2000)126:1(39)
  6. Kirby, J.T., Wei, G., Chen, Q., Kennedy, A.B. and Dalrymple, R.A. (1998). Fully nonlinear Boussinesq wave model, Documentation and User's Manual. Technical Report CACR-98-06, University of Delaware.
  7. Kirby, J.T., Long, W. and Shi, F. (2003). Fully nonlinear Boussinesq wave model on curvilinear coordinates, Documentation and User's Manual. Technical Report CACR-02-xx, University of Delaware.
  8. Losada, I.J., Gonzalez-Ondina, J.M., Diaz, G. and Gonzalez, E.M. (2008). Numerical modeling of nonlinear resonance of semienclosed water bodies: Description and experimental validation. Coastal Engineering, 55, 21-34. https://doi.org/10.1016/j.coastaleng.2007.06.002
  9. Nwogu, O. (1993). Alternative form of Boussinesq equations for nearshore wave propagation. Journal of Waterway, Port, Coastal, and Ocean Engineering, 119(6), 618-638. https://doi.org/10.1061/(ASCE)0733-950X(1993)119:6(618)
  10. Toro, E.F. (1999). Riemann solvers and numerical methods for fluid dynamics – A practical introduction. Springer, Berlin.
  11. Vincent, C.L. and Briggs, M.J. (1989). Refraction-diffraction of irregular waves over a mound. Journal of Waterway, Port, Coastal,and Ocean Engineering, 115(2), 269-284. https://doi.org/10.1061/(ASCE)0733-950X(1989)115:2(269)
  12. Walkley, M. and Berzins, M. (2002). A finite element method for the two-dimensional extended Boussinesq equations. International Journal for Numerical Methods in Fluids, 39, 865-885 https://doi.org/10.1002/fld.349
  13. Wei, G. and Kirby, J.T. (1995). Time dependent numerical code for extended Boussinesq equations. Journal of Waterway, Port, Coastal, and Ocean Engineering, 121(5), 251-261. https://doi.org/10.1061/(ASCE)0733-950X(1995)121:5(251)
  14. Woo, S.-B. and Liu, P.-H. (2004). Finite element model for modified Boussinesq equations. I: Model development. Journal of Waterway, Port, Coastal, and Ocean Engineering, 130(1), 1-16. https://doi.org/10.1061/(ASCE)0733-950X(2004)130:1(1)
  15. Yoon, S.B., Cho, Y.S. and Lee, C. (2004). Effects of breakinginduced currents on refraction-diffraction of irregular waves over submerged shoal. Ocean Engineering, 31, 633-652. https://doi.org/10.1016/j.oceaneng.2003.07.008
  16. Zhong, Z and Wang, K.H. (2007). Time-accurate stabilized finiteelement model for weakly nonlinear and weakly dispersive water waves. International Journal for Numerical methods in fluids, accepted.