Search

We propose a simple boundary condition regularization strategy to reduce error propagation in pressure field reconstruction from corrupted image velocimetry data. The core idea is to replace the canonical Neumann boundary conditions with Dirichlet ones obtained by integrating the tangential part of the pressure gradient along the boundaries. Rigorous analysis and numerical experiments justify the effectiveness of this regularization. [ABSTRACT FROM AUTHOR]

The dynamical core of an atmospheric general circulation model is engineered to satisfy a delicate balance between numerical stability, computational cost, and an accurate representation of the equations of motion. It generally contains either explicitly added or inherent numerical diffusion mechanisms to control the buildup of energy or enstrophy at the smallest scales. The diffusion fosters computational stability and is sometimes also viewed as a substitute for unresolved subgrid-scale processes. A particular form of explicitly added diffusion is horizontal divergence damping. In this paper a von Neumann stability analysis of horizontal divergence damping on a latitude-longitude grid is performed. Stability restrictions are derived for the damping coefficients of both second- and fourth-order divergence damping. The accuracy of the theoretical analysis is verified through the use of idealized dynamical core test cases that include the simulation of gravity waves and a baroclinic wave. The tests are applied to the finite-volume dynamical core of NCAR's Community Atmosphere Model version 5 (CAM5). Investigation of the amplification factor for the divergence damping mechanisms explains how small-scale meridional waves found in a baroclinic wave test case are not eliminated by the damping. [ABSTRACT FROM AUTHOR]