Flux-form Eulerian-Lagrangian method for solving advective transport of scalars in free-surface flows

被引：0
作者：
Hu D. [1 ,2 ]
Yao S. [2 ]
Qu G. [2 ]
Zhong D. [3 ]
机构：
[1] School of Hydropower and Information Engineering, Huazhong Univ. of Science and Technology, Wuhan
[2] Dept. of River Engineering, Yangtze River Scientific Research Institute, 23 Huangpu St., Wuhan
[3] State Key Laboratory of Hydroscience and Engineering, Tsinghua Univ., Beijing
来源：
Journal of Hydraulic Engineering | 2019年 / 145卷 / 03期
关键词：
Conservative Eulerian-Lagrangian method; Free-surface flow; Multiscalar transport; Parallel computing; Scalar transport; Unstructured grid;
D O I：
10.1061/(asce)hy.1943-7900.0001578
中图分类号：
学科分类号：
摘要：
A two-dimensional (2D) flux-form Eulerian-Lagrangian method (FFELM) on unstructured grid is proposed for solving the advection equation in free-surface scalar transport models. The scalar concentrations of backtracking points are combined with timeinterpolated cell-face velocities to evaluate cell-face advective fluxes. A G-correction is defined as an additional mechanism to eliminate potential nonphysical oscillations by correcting the cell-face advective fluxes. A flux-form cell update is finally carried out to obtain new cell concentrations. The role of the G-correction in the FFELM is clarified using a test of scalar transport in unsteady open-channel flows. A solidbody rotation test, a laboratory bend-flume test, and a real river test (using a 600-km river reach of the upper Yangtze River) are used to demonstrate the FFELM. The FFELM is revealed in tests to achieve almost the same accuracy as a pure Eulerian-type method [the subcycling finite-volume method (SCFVM)] and a conservative ELM [the finite-volume ELM (FVELM)]. Relative to explicit Eulerian methods, the FFELM uses the information of backtracking points over an extended upwind dependence domain in evaluating cell-face advective fluxes, and allows larger time steps for which the Courant-Friedrichs-Lewy number (CFL) is greater than 1. In the real river test, stable and accurate FFELM simulations can be achieved at a time step for which the CFL is as large as 5. Efficiency issues of the FFELM are clarified using the bend-flume test (193,536 cells) and the real river test (213,363 cells). In solving a transport problem (using 1-32 kinds of scalars and 16 cores), a parallel run using the FFELM is 1.0-3.3 times faster than a parallel run using the SCFVM. The FFELM has a computational cost slightly less (15%-17%) than that of the FVELM. Moreover, the implementation of the FFELM is much easier than that of the FVELM, and extension of the 2D FFELM to its one-dimensional (1D) and three-dimensional (3D) versions is straightforward. © 2019 American Society of Civil Engineers.
引用
收藏
相关论文
共 26 条
  • [1] Budgell W.P., Oliveira A., Skogen M.D., Scalar advection schemes for ocean modeling on unstructured triangular grids., Ocean Dyn, 57, 4, pp. 339-361, (2007)
  • [2] Cao Z.X., Hu P.H., Hu K.H., Pender G., Liu Q., Modelling roll waves with shallow water equations and turbulent closure., J. Hydraul. Res, 53, 2, pp. 161-177, (2015)
  • [3] Casulli V., Zanolli P., Semi-implicit numerical modeling of nonhydrostatic free-surface flows for environmental problems., Math. Comput. Modell, 36, 9-10, pp. 1131-1149, (2002)
  • [4] Chang Y.C., "Lateral mixing in meandering channels.", (1971)
  • [5] Dimou K., "3-D hybrid Eulerian-Lagrangian/particle tracking model for simulating mass transport in coastal water bodies.", (1992)
  • [6] Dukowicz J.K., Baumgardner J.R., Incremental remapping as a transport/advection algorithm., J. Comput. Phys, 160, 1, pp. 318-335, (2000)
  • [7] Falconer R.A., Wolanski E., Mardapitta-Hadjipandeli L., Modeling tidal circulation in an island's wake., J. Waterw. Port Coastal Ocean Eng, 112, 2, pp. 234-254, (1986)
  • [8] Fischer H.B., Longitudinal dispersion and turbulent mixing in open channel flow., Ann. Rev. Fluid Mech, 5, 1, pp. 59-78, (1973)
  • [9] Gross E.S., Bonaventura L., Giorgio R., Consistency with continuity in conservative advection schemes for free-surface models., Int. J. Numer. Methods Fluids, 38, 4, pp. 307-327, (2002)
  • [10] Harris L.M., Lauritzen P.H., Mittal R., A flux-form version of the conservative semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed sphere grid., J. Comput. Phys, 230, 4, pp. 1215-1237, (2011)
123