UTBEST3D (University of Texas Bays and Estuaries Simulator — 3D)  

UTBEST3D (University of Texas Bays and Estuaries Simulator — 3D)
The numerical solution algorithm in UTBEST3D considers the system of hydrostatic primitive equations (1) with a free surface. A prismatic mesh (Fig. 1) is obtained by projecting a given triangular mesh in the vertical direction to provide a continuous piecewise linear representations of the topography and of the free surface. The vertical columns are then subdivided into layers. If a bottommost prism is degenerate, it is merged with the one above it. Due to the discontinuous nature of the approximation spaces, no constraints need to be enforced on the element connectivity. Hanging nodes and mismatching elements are allowed and have no adverse effects on stability or conservation properties of the scheme. This flexibility with regard to mesh geometry is exploited in several key parts of the algorithm: vertical mesh construction in areas with varying topography (Fig. 4), local mesh adaptivity (Fig. 2), wetting/drying (Fig. 3).
$\backslash partial\_t\; \{\backslash bf\; u\}\_H\; +\backslash nabla\; \backslash cdot\; \backslash left(\; \{\backslash bf\; u\}\_H\; \{\backslash bf\; u\}\; \; D\_\{\backslash bf\; u\}\; \backslash nabla\{\backslash bf\; u\}\; \backslash right)\; +\; \backslash nabla\_H\; P\; f\_C\; \{\backslash bf\; k\}\; \backslash times\; \{\backslash bf\; u\}\_H\; =0$
$\backslash partial\_t\; T\; +\; \backslash nabla\; \backslash cdot\; \backslash left(\{\backslash bf\; u\}\; T\; \; D\_T\; \backslash nabla\; T\; \backslash right)\; =\; 0$
$\backslash partial\_t\; S\; +\; \backslash nabla\; \backslash cdot\; \backslash left(\{\backslash bf\; u\}\; S\; \; D\_S\; \backslash nabla\; S\; \backslash right)\; =\; 0$
$\backslash partial\_t\; \backslash xi\; +\; \backslash nabla\_H\; \backslash cdot\; \backslash int\_\{z\_b\}^\{\backslash xi\}\; \{\backslash bf\; u\}\; dz\; =0\backslash hspace\{45mm\}(1)$
$\backslash nabla\; \backslash cdot\; \{\backslash bf\; u\}\; =0$
$\backslash partial\_z\; P=\; \; \backslash frac\{\backslash rho\}\{\backslash rho\_0\}g$
$\backslash rho\; =\; \backslash rho(T,S,z)$
$\backslash nabla\_H\; =\; (\backslash partial\_x,\backslash partial\_y)$, $\{\backslash bf\; k\}\; =\; (0,0,1)^T$, $D\_\{X,\; X\; \backslash in\; \backslash \{\{\backslash bf\; u\},S,T\backslash \}\}$ – eddy viscosity/diffusivity, $\backslash xi$ – free surface elevation, $z\_b$ – bathymetry, $\{\backslash bf\; u\}=\; (u,v,w)$ – velocity, $\{\backslash bf\; u\}\_H=\; (u,v)$, T – temperature, S – salinity, P – pressure, $\backslash rho$ – density, $\backslash rho\_0$ – reference density, g – gravity acceleration, $f\_C$ – Coriolis coefficient.
Our DG discretization is based on the local discontinuous Galerkin method (Cockburn, Shu, 1998) that represents a direct generalization of the cellcentered finite volume method, the latter being just the piecewise constant DG discretization. One of the features of this method is a much smaller numerical diffusion exhibited by the linear and higher order DG approximations compared to the finite difference or finite volume discretizations. The method guarantees the elementwise conservation of all primary unknowns including tracers, supports an individual choice of the approximation space for each prognostic and diagnostic variable, demonstrates excellent stability properties, and possesses proven local adaptivity skills (Fig. 2).
UTBEST3D is written in C++ to provide clean interfaces between geometrical, numerical, computational, and communication parts of the code. The objectoriented coding paradigm is designed to enable a labor efficient development lifecycle of the model. The programming techniques were carefully chosen and tested with the view of guaranteeing a smooth portability to different hardware architectures, operating systems, compilers, and software environments; the package relies on no external libraries in the serial configuration and only makes use of METIS (Karypis, Kumar, 1999) grid partitioning library in the parallel configuration.
Bibliography
Vadym Aizinger, A discontinuous Galerkin method for two and threedimensional shallowwater equations, PhD Dissertation, THE UNIVERSITY OF TEXAS AT AUSTIN, 2004.
V. Aizinger, C. Dawson, The local discontinuous Galerkin method for threedimensional shallow water flow, Computer methods in applied mechanics and engineering, 196, 734746, 2007.
B. Cockburn and C.W. Shu, The local discontinuous Galerkin method for timedependent convectiondiffusion systems, SIAM J. Numer. Anal. 35, 24402463, 1998.
 C. Dawson, V. Aizinger, A Discontinuous Galerkin Method for ThreeDimensional Shallow Water Equations, Journal of Scientific Computing, 22, 245267, 2005.
 G. Karypis and V. Kumar, A fast and high quality multilevel scheme for partitioning irregular graphs, SIAM Journal on Scientific Computing, 20(1), 359392, 1999.