A method for post-processing the velocity after a pressure projection is developed that helps to maintain stability in an under-resolved, inviscid, discontinuous element-based simulation for use in environmental fluid mechanics process studies. The post-processing method is needed because of spurious divergence growth at element interfaces due to the discontinuous nature of the discretization used. This spurious divergence eventually leads to a numerical instability. Previous work has shown that a discontinuous element-local projection onto the space of divergence-free basis functions is capable of stabilizing the projection method, but the discontinuity inherent in this technique may lead to instability in under-resolved simulations. By enforcing inter-element discontinuity and requiring a divergence-free result in the weak sense only, a new post-processing technique is developed that simultaneously improves smoothness and reduces divergence in the pressure-projected velocity field at the same time. When compared against a non-post-processed velocity field, the post-processed velocity field remains stable far longer and exhibits better smoothness and conservation properties.

On the New Jersey continental shelf ambient sound levels were recorded during tropical storm Ernesto that produced wind speeds up to 40 knots in early September 2006. The seabed at the position of the acoustic measurements can be approximately described as coarse sand. Differences between the ambient noise levels for the New Jersey shelf measurements and deep water reference measurements are modeled using both normal mode and ray methods. The analysis is consistent with a nonlinear frequency dependent seabed attenuation for the New Jersey site.

No abstract

A combination of block-Jacobi and deflation preconditioning is used to solve a high-order discontinuous collocation-based discretization of the Schur complement of the Poisson-Neumann system as arises in the operator splitting of the incompressible Navier-Stokes equations. The preconditioners and deflation vectors are chosen to mitigate the effects of ill-conditioning due to highly-elongated domains and to achieve GMRES convergence independent of the size of the grid. The ill-posedness of the PoissonNeumann system manifests as an inconsistency of the Schur complement problem, but it is shown that this can be accounted for with appropriate projections out of the null space of the Schur matrix without affecting the accuracy of the solution. The combined deflation/block-Jacobi preconditioning is compared with two-level non-overlapping additive Schwarz preconditioning of the Schur problem, and while both methods achieve convergence independent of the grid size, deflation is shown to require half as many GMRES iterations and 25% less wall-clock time for a variety of grid sizes and domain aspect ratios. The deflation methods shown to be effective for the two-dimensional Poisson-Neumann problem are extensible to the three-dimensional problem assuming a Fourier discretization in the third dimension. A Fourier discretization results in a two-dimensional Helmholtz problem for each Fourier component that is solved using deflation/block-Jacobi preconditioning on its Schur complement. Here again deflation is shown to be superior to two-level non-overlapping additive Schwarz preconditioning, requiring about half as many GMRES iterations and 15% less time. While the methods here are demonstrated on a spectral multidomain penalty method discretization, they are readily extensible to any discontinuous elementbased discretization of an elliptic problem, and are particularly well-suited for high-order methods.