Related papers: Simulating linear kinematic features in viscous-pl…
Linear Kinematic Features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. Viscous-plastic sea-ice models start to produce…
For the numerical simulation of earth system models, Arakawa grids are largely employed. A quadrilateral mesh is assumed for their classical definition, and different types of grids are identified depending on the location of the…
We present a stable discretization of sea-ice dynamics on triangular grids that can straightforwardly be coupled to an ocean model on a triangular grid with Arakawa C-type staggering. The approach is based on a nonconforming finite element…
Discretization of the equations of Viscous Plastic and Elastic Viscous Plastic (EVP) sea ice dynamics on triangular meshes can be done by placing discrete velocities at vertices, cells or edges. Since there are more cells and edges than…
Finite-difference (FD) modeling of seismic waves in the vicinity of dipping interfaces gives rise to artifacts. Examples are phase and amplitude errors, as well as staircase diffractions. Such errors can be reduced in two general ways. In…
We present an approach to simulate the 3D isotropic elastic wave propagation using nonuniform finite difference discretization on staggered grids. Specifically, we consider simulation domains composed of layers of uniform grids with…
As global ocean general circulation models are run at eddy-permitting resolutions, reproducing accurate growth rates of baroclinic instabilities is a major concern when choosing a discretization of the equations of motion. From this…
Continuum computational kinetic plasma models evolve the distribution function of a plasma species $f_s$ on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence,…
We propose an extension of the discretization approaches for multilayer shallow water models, aimed at making them more flexible and efficient for realistic applications to coastal flows. A novel discretization approach is proposed, in…
Sea ice mechanical properties are involved in dynamical processes acting from the scale of meters to several hundred kilometers. The current rapid changes in the state of polar sea ice require a better understanding and modeling of these…
Arctic sea ice plays an important role in the global climate. Sea ice models governed by physical equations have been used to simulate the state of the ice including characteristics such as ice thickness, concentration, and motion. More…
In this note, the importance of spectral properties of viscous flux discretization in solving compressible Navier-Stokes equations for turbulent flow simulations is discussed. We studied six different methods, divided into two different…
The problem of natural convection in a laterally heated three-dimensional cubic cavity under the action of an externally imposed magnetic field is revisited. Flows at the Rayleigh number Ra=10^6 and the Hartmann number Ha=100, and three…
We address a numerical instability that arises in the directionally split computation of hydrodynamic flows when shock fronts are parallel to a grid plane. Transverse oscillations in pressure, density and temperature are produced that are…
Flexure of Antarctic ice shelves under excitation from long ocean waves induces mechanical ice shelf stresses that amplify fractures and, hence, contribute to calving events. Here, a solution method is developed for a hydroelastic…
The marginal ice zone is a highly dynamical region where sea ice and ocean waves interact. Large-scale sea ice models only compute domain-averaged responses. As the majority of the marginal ice zone consists of mobile ice floes surrounded…
A significant fraction (4%-13%) of Antarctic sea ice remains stationary as landfast sea-ice ("fast ice"), typically anchored by grounded icebergs. Current global climate models do not represent fast-ice formation due to iceberg grounding,…
We present an efficient matrix-free geometric multigrid method for the elastic Helmholtz equation, and a suitable discretization. Many discretization methods had been considered in the literature for the Helmholtz equations, as well as many…
An ice shelf is a floating extension of a land-based ice sheet into the ocean. It plays a crucial role in slowing down the flow of land ice into the sea, thus stabilizing the ice sheet. However, this stabilizing effect can be weakened by…
We develop a new meshfree geometric multilevel (MGM) method for solving linear systems that arise from discretizing elliptic PDEs on surfaces represented by point clouds. The method uses a Poisson disk sampling-type technique for coarsening…