×

Semi-implicit, semi-Lagrangian modelling for environmental problems on staggered Cartesian grids with cut cells. (English) Zbl 1143.76484

Summary: We propose an extension of semi-implicit models for the shallow water or Euler equations on staggered Cartesian meshes to handle the presence of cut cells at the boundaries of the computational domain. Semi-Lagrangian advection algorithms are also extended to staggered Cartesian meshes with cut boundary cells. A linear reconstruction algorithm and radial basis function interpolation are used, respectively, to compute approximate trajectories and to recover an accurate velocity field in cells intersected by the boundary. The effective accuracy of radial basis function interpolators for semi-Lagrangian algorithms is assessed by simple numerical experiments. Timestep control techniques in the computation of the semi-Lagrangian trajectories are also discussed. The accuracy of the resulting models for environmental flows is demonstrated by a number of test cases simulating nonlinear open channel flows and baroclinic flows over an isolated obstacle.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
Full Text: DOI

References:

[1] Adcroft, A.; Hill, C.; Marshall, J., Representation of topography by shaved cells in a height coordinate ocean model, Monthly Weather Review, 125, 2293-2315 (1997)
[2] Bartello, P.; Thomas, S. J., The cost effectiveness of semi-Lagrangian advection, Monthly Weather Review, 124, 2883-2895 (1996)
[3] Behrens, J.; Iske, A., Grid-Free Adaptive semi-Lagrangian advection using radial basis functions, Computers and Mathematics with Applications, 43, 319-327, (2002) · Zbl 0999.65104
[4] Bonaventura, L., A semi-implicit, semi-Lagrangian scheme using the height coordinate for a nonhydrostatic and fully elastic model of atmospheric flows, Journal of Computational Physics, 158, 186-213 (2000) · Zbl 0963.76058
[5] Bonaventura, L., A second order, semi-Lagrangian scheme with accurate approximation of trajectories, (Taylor, C.; Cross, J., Proceedings of the 10th International Conference on Numerical Methods in Fluids (1997), Pineridge Press: Pineridge Press Swansea)
[6] Buhmann, M. D., Radial basis functions, Acta Numerica, 1-38 (2000) · Zbl 1004.65015
[7] Calhoun, D.; LeVeque, R. J., A Cartesian grid finite-volume method for the advection diffusion equation in irregular geometry, Journal of Computational Physics, 157, 143-180 (2000) · Zbl 0952.65075
[8] Casulli, V., Semi-implicit finite difference methods for the two dimensional shallow water equations, Journal of Computational Physics, 86, 56-74 (1990) · Zbl 0681.76022
[9] Casulli, V.; Cattani, E., Stability, accuracy and efficiency of a semi-implicit method for three-dimensional shallow water flow, Computers and Mathematics with Applications, 27, 99-112 (1994) · Zbl 0796.76052
[10] Casulli, V.; Stelling, G. S., Numerical simulation of 3D quasi-hydrostatic free-surface flows, Journal of Hydraulic Engineering ASCE, 124, 678-686 (1998)
[11] Casulli, V., A semi-implicit finite difference method for non-hydrostatic, free-surface flows, International Journal of Numerical Methods in Fluids, 30, 425-440 (1999) · Zbl 0944.76050
[12] Causon, D. M.; Ingram, D. M.; Mingham, C. G., Calculation of shallow water flows using a Cartesian cut cell approach, Advances in Water Resources Research, 23, 545-562 (2000)
[13] Causon, D. M.; Ingram, D. M.; Mingham, C. G., A Cartesian cut cell method for shallow water flows with moving boundaries, Advances in Water Resources Research, 24, 899-911 (2001)
[14] Chow, V. T., Open Channel Hydraulics (1959), McGraw-Hill: McGraw-Hill New York
[15] Dudhia, J., A nonhydrostatic version of the Penn State - NCAR mesoscale model: validation, tests and simulation of an Atlantic cyclone and cold front, Monthly Weather Review, 121, 1493-1513 (1993)
[16] Gallus, W.; Klemp, J., Behaviour of flow over step orography, Monthly Weather Review, 128, 1153-1164 (2000)
[17] Giraldo, F. X., Trajectory computations for spherical geodesic grids in Cartesian space, Monthly Weather Review, 127, 1651-1662 (1999)
[18] Giraldo, F. X., The Lagrange-Galerkin method for the two-dimensional shallow water equations on adaptive grids, International Journal of Numerical Methods in Fluids, 33, 789-832 (2000) · Zbl 0989.76047
[19] Gross, E. S.; Koseff, J. R.; Monismith, S. G., Three-dimensional salinity simulations of South San Francisco Bay, ASCE Journal of Hydraulic Engineering, 125, 1199-1209 (1999)
[20] Harlow, F. H.; Welch, J. E., Numerical calculation of time dependent viscous incompressible flow, Physics of Fluids, 8, 2182-2189 (1965) · Zbl 1180.76043
[21] Henderson, F. M., Open Channel Flow (1966), MacMillan: MacMillan New York · Zbl 0193.55002
[22] The Hamburg Ocean Primitive Equation Model HOPE, DKRZ Report No. 13; The Hamburg Ocean Primitive Equation Model HOPE, DKRZ Report No. 13
[23] Iske, A., On the approximation order and numerical stability of local Lagrange interpolation by polyharmonic splines, (Haussmann, W.; Jetter, K.; Reimer, M.; Stöckler, J., Modern Developments in Multivariate Approximation. Modern Developments in Multivariate Approximation, International Series of Numerical Mathematics (2003), Birkhäuser Verlag: Birkhäuser Verlag Basel), 153-165 · Zbl 1040.41002
[24] Johansen, H.; Colella, P., A Cartesian grid embedded boundary method for Poisson’s equation on irregular domains, Journal of Computational Physics, 147, 60-85 (1998) · Zbl 0923.65079
[25] Kasahara, A.; Washington, W. W., General circulation experiment with a six-layer NCAR model, Monthly Weather Review, 28, 657-701 (1971)
[26] Le Roux, D. Y.; Lin, C. C.; Staniforth, A., An accurate interpolating scheme for semi-Lagrangian advection on an unstructured grid for ocean modelling, Tellus A, 49, 119-138 (1997)
[27] Le Roux, D. Y.; Lin, C. C.; Staniforth, A., A semi-implicit, semi-Lagrangian finite element shallow water ocean model, Monthly Weather Review, 128, 1384-1401 (1998)
[28] Li, X.; Micchelli, C. A., Approximation by radial bases and neural networks, Numerical Algorithms, 25, 242-262 (2000) · Zbl 0976.65014
[29] Miglio, E.; Quarteroni, A.; Saleri, F., Finite element approximation of quasi-3d shallow water equations, Comparative Methods in Applied Mechanics and Engineering, 174, 355-369 (1999) · Zbl 0958.76046
[30] Pinty, J. P.; Benoit, R.; Richard, E.; Laprise, R., Simple tests of a semi-implicit semi-Lagrangian model on 2D mountain wave problems, Monthly Weather Review, 123, 3042-3058 (1995)
[31] Pironneau, O., On the transport-diffusion algorithm and its applications to the Navier-Stokes equations, Numerische Mathematik, 38, 309-332 (1982) · Zbl 0505.76100
[32] Schaback, R.; Wu, Z., Local error estimates for radial basis function interpolation of scattered data, IMA Journal of Numerical Analysis, 13, 13-27 (1993) · Zbl 0762.41006
[33] Schär, C.; Leuenberger, O.; Fuhrer, D.; Lüthi, D.; Girard, C., A new vertical coordinate formulation for atmospheric prediction models, Monthly Weather Review, 130, 2459-2480 (2002)
[34] Smolarkiewicz, P. K.; Pudykiewicz, J., A class of semi-Lagrangian approximations for fluids, Journal of Atmospheric Sciences, 49, 2082-2096 (1992)
[35] Sonar, T., Optimal recovery using thin plate splines in finite volume methods for the numerical solution of hyperbolic conservation laws, IMA Journal of Numerical Analysis, 16, 549-581 (1996) · Zbl 0860.65096
[36] Staniforth, A.; Coté, J., Semi-Lagrangian integration schemes for atmospheric models - a review, Monthly Weather Review, 119, 2206-2223 (1991)
[37] Steppeler, J.; Bitzer, H. W.; Minotte, M.; Bonaventura, L., Nonhydrostatic atmospheric modelling using a \(Z\) coordinate representation, Monthly Weather Review, 130, 2143-2149 (2002)
[38] Steppeler, J.; Doms, G.; Schättler, U., Meso-gamma scale forecasts using the nonhydrostatic model LM, Meteorology and Atmospheric Physics, 82, 75-96 (2003)
[39] Wang, P. F.; Cheng, R. T.; Richter, K.; Gross, E. S.; Sutton, D.; Gartner, J. W., Modeling tidal hydrodynamics of San Diego Bay California, Journal of the American Water Resources Association, 34, 1123-1140 (1998)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.