Related papers: Dynamically coupling the non-linear Stokes equatio…
The Shallow Ice Approximation (SIA) model on strong form is commonly used for inferring the flow dynamics of grounded ice sheets. The solution to the SIA model is a closed-form expression for the velocity field. When that velocity field is…
Numerical models for predicting future ice-mass loss of the Antarctic and Greenland ice sheet requires accurately representing their dynamics. Unfortunately, ice-sheet models suffer from a very strict time-step size constraint, which for…
We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past…
A highly accurate method for simulating surfactant-covered droplets in two-dimensional Stokes flow with solid boundaries is presented. The method handles both periodic channel flows of arbitrary shape and stationary solid constrictions. A…
The velocity field and the height at the surface of a dynamic ice sheet are observed. The ice sheets are modeled by the full Stokes equations and shallow shelf/shelfy stream approximations. Time dependence is introduced by a kinematic free…
A key parameter in ice flow modelling is basal slipping at the ice-bed interface as it can have a large effect on the resultant ice thickness. Unfortunately, its contribution to surface observations can be hard to distinguish from that of…
We present the Thickness Stabilization Scheme (TSS), a numerical stabilization scheme suitable for the Shallow Shelf Approximation (SSA), one of the most widely-used models for large-scale Antarctic and Greenland ice sheet simulations. The…
The shallow shelf approximation is a better ``sliding law'' for ice sheet modeling than those sliding laws in which basal velocity is a function of driving stress. The shallow shelf approximation as formulated by \emph{Schoof} [2006a] is…
We present a coupling framework for Stokes-Darcy systems valid for arbitrary flow direction at low Reynolds numbers and for isotropic porous media. The proposed method is based on an overlapping domain decomposition concept to represent the…
Ice sheet models are routinely used to quantify and project an ice sheet's contribution to sea level rise. In order for an ice sheet model to generate realistic projections, its parameters must first be calibrated using observational data;…
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…
The interaction between thin structures and incompressible Newtonian fluids is ubiquitous both in nature and in industrial applications. In this paper we present an isogeometric formulation of such problems which exploits a boundary…
In this work we present an efficient numerical algorithm for the solution of interfacial statistical associating fluid theory (iSAFT) in cylindrical geometry to facilitate the study of inhomogeneous fluids having curvatures. The new…
Sea ice plays a crucial role in the climate system, particularly in the Marginal Ice Zone (MIZ), a transitional area consisting of fragmented ice between the open ocean and consolidated pack ice. As the MIZ expands, understanding its…
In seasonally ice-covered seas and along the margins of perennial ice pack, i.e. in regions with medium ice concentrations, the ice cover typically consists of separate floes interacting with each other by inelastic collisions. In this…
We implement a data assimilation framework for integrating ice surface and terminus position observations into a numerical ice-flow model. The model uses the well-known shallow shelf approximation (SSA) coupled to a level set method to…
Performing highly accurate simulations of droplet systems is a challenging problem. This is primarily due to the interface dynamics which is complicated further by the addition of surfactants. This paper presents a boundary integral method…
(Abridged) We use the statistical tool known as the ``Spectral Correlation Function" [SCF] to intercompare simulations and observations of the atomic interstellar medium. The simulations considered mimic three distinct sets of physical…
State estimation in multi-layer turbulent flow fields with only a single layer of partial observation remains a challenging yet practically important task. Applications include inferring the state of the deep ocean by exploiting surface…
Smoothed Particle Hydrodynamics (SPH) is a popular numerical technique developed for simulating complex fluid flows. Among its key ingredients is the use of nonlocal integral relaxations to local differentiations. Mathematical analysis of…