Related papers: Weak form Shallow Ice Approximation models with an…
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 propose and implement a new method, called the Ice Sheet Coupled Approximation Levels (ISCAL) method, for simulation of ice sheet flow in large domains under long time-intervals. The method couples the exact, full Stokes (FS) equations…
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…
Numerical ice sheet models compute evolving ice geometry and velocity fields using various stress-balance approximations and boundary conditions. At high spatial resolution, with horizontal mesh/grid resolutions of a few kilometers or…
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…
We adopt the frozen Gaussian approximation (FGA) for modeling seismic waves. The method belongs to the category of ray-based beam methods. It decomposes seismic wavefield into a set of Gaussian functions and propagates these Gaussian…
A nonlinear sea-ice problem is considered in a least-squares finite element setting. The corresponding variational formulation approximating simultaneously the stress tensor and the velocity is analysed. In particular, the least-squares…
The paper presents advancement of the matrix structural analysis technique (MSA) for stiffness modeling of robotic manipulators. In contrast to the classical MSA, it can be applied to both parallel and serial manipulators composed of…
We propose a variant of the $\theta$-scheme for diffuse interface models for two-phase flow, together with three new linearization techniques for the surface tension. These involve either additional stabilizing force terms, or a fully…
Surface integral equation (SIE) methods are of great interest for the efficient electromagnetic modeling of various devices, from integrated circuits to antenna arrays. Existing acceleration algorithms for SIEs, such as the adaptive…
Robust, broadly applicable fluid-structure interaction (FSI) algorithms remain a challenge for computational mechanics. In previous work, we introduced an immersed interface method (IIM) for discrete surfaces and an extension based on an…
An efficient, iterative semi-implicit (SI) numerical method for the time integration of stiff wave systems is presented. Physics-based assumptions are used to derive a convergent iterative formulation of the SI scheme which enables the…
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…
Shallow free surface flows are often characterized by both subdomains that require high modeling complexity and subdomains that can be sufficiently accurately modeled with low modeling complexity. Moreover, these subdomains may change in…
Recent applications (e.g. active gels and self-assembly of elastic sheets) motivate the need to efficiently simulate the dynamics of thin elastic sheets. We present semi-implicit time stepping algorithms to improve the time step constraints…
Two particle irreducible effective actions (2PIEAs) are valuable non-perturbative techniques in quantum field theory; however, finite truncations of them violate the Ward identities (WIs) of theories with spontaneously broken symmetries.…
The Cahn-Hilliard equation is the most common model to describe phase separation processes of a mixture of two components. For a better description of short-range interactions of the material with the solid wall, various dynamic boundary…
The theory of slow invariant manifolds (SIMs) is the foundation of various model-order reduction techniques for dissipative dynamical systems with multiple time-scales, e.g. in chemical kinetic models. The construction of SIMs and many…
A non-hydrostatic depth-averaged model for dry granular flows is proposed, taking into account vertical acceleration. A variable friction coefficient based on the $\mu(I)$ rheology is considered. The model is obtained from an asymptotic…