Finite elements numerical solution of a coupled profile–velocity–temperature shallow ice sheet approximation model

AbstractThis work deals with the numerical solution of a complex mathematical model arising in theoretical glaciology. The global moving boundary problem governs thermomechanical processes jointly with ice sheet hydrodynamics. One major novelty is the inclusion of the ice velocity field computation in the framework of the shallow ice model so that it can be coupled with profile and temperature equations. Moreover, the proposed basal velocity and shear stress laws allow the integration of basal sliding effects in the global model. Both features were not taking into account in a previous paper (Math. Model. Methods Appl. Sci. 12 (2) (2002) 229) and provide more realistic convective terms and more complete Signorini boundary conditions for the thermal problem. In the proposed numerical algorithm, one- and two-dimensional piecewise linear Lagrange finite elements in space and a semi-implicit upwinding scheme in time are combined with duality and Newton's methods for nonlinearities. A simulation example involving real data issued from Antarctic shows the temperature, profile and velocity qualitative behaviour as well as the free boundaries and basal effects

A nonsmooth Newton multigrid method for a hybrid, shallow model of marine ice sheets

The time evolution of ice sheets and ice shelves is model by combining a shallow lubrication approximation for shear deformation with the shallow shelf approximation for basal sliding, along with the mass conservation principle. At each time step two p-Laplace problems and one transport problem are solved. Both p-Laplace problems are formulated as minimisation problems. They are approximated by a finite element truncated nonsmooth Newton multigrid method. As an illustration, we compute the steady state shape of an idealized ice sheet/shelf system

Lessons from the short history of ice sheet model intercomparison

International audienceIntercomparison should include measurement of differences, between model results and observations, among the model results themselves, or between model results and exact solutions. The processes of measuring differences and critically analyzing those differences are vital. Without such measurement as a component of intercomparison, the only expected benefits of an intercomparison project are participation, possibly the discovery of communal confusion, and the establishment of public, non-proprietary data sets

The Fisher-KPP problem with doubly nonlinear "fast" diffusion

The famous Fisher-KPP reaction diffusion model combines linear diffusion with the typical Fisher-KPP reaction term, and appears in a number of relevant applications. It is remarkable as a mathematical model since, in the case of linear diffusion, it possesses a family of travelling waves that describe the asymptotic behaviour of a wide class solutions $0\leq u(x,t)\leq 1$ of the problem posed in the real line. The existence of propagation wave with finite speed has been confirmed in the cases of "slow" and "pseudo-linear" doubly nonlinear diffusion too, see arXiv:1601.05718. We investigate here the corresponding theory with "fast" doubly nonlinear diffusion and we find that general solutions show a non-TW asymptotic behaviour, and exponential propagation in space for large times. Finally, we prove precise bounds for the level sets of general solutions, even when we work in with spacial dimension $N \geq 1$. In particular, we show that location of the level sets is approximately linear for large times, when we take spatial logarithmic scale, finding a strong departure from the linear case, in which appears the famous Bramson logarithmic correction.Comment: 42 pages, 6 figure

Conservation laws for free-boundary fluid layers

Time-dependent models of fluid motion in thin layers, subject to signed source terms, represent important sub-problems within climate dynamics. Examples include ice sheets, sea ice, and even shallow oceans and lakes. We address these problems as discrete-time sequences of continuous-space weak formulations, namely (monotone) variational inequalities or complementarity problems, in which the conserved quantity is the layer thickness. Free boundaries wherein the thickness and mass flux both go to zero at the margin of the fluid layer generically arise in such models. After showing these problems are well-posed in several cases, we consider the limitations to discrete conservation or balance in numerical schemes. A free boundary in a region of negative source -- an ablation-caused margin -- turns out to be a barrier to exact balance for a numerical scheme (in either a continuous- or discrete-space sense). We propose computable \emph{a posteriori} quantities which allow conservation-error accounting in finite volume and element schemes.Comment: 26 pages, 4 figure

Bistable reaction equations with doubly nonlinear diffusion

Reaction-diffusion equations appear in biology and chemistry, and combine linear diffusion with different kind of reaction terms. Some of them are remarkable from the mathematical point of view, since they admit families of travelling waves that describe the asymptotic behaviour of a larger class of solutions $0\leq u(x,t)\leq 1$ of the problem posed in the real line. We investigate here the existence of waves with constant propagation speed, when the linear diffusion is replaced by the "slow" doubly nonlinear diffusion. In the present setting we consider bistable reaction terms, which present interesting differences w.r.t. the Fisher-KPP framework recently studied in \cite{AA-JLV:art}. We find different families of travelling waves that are employed to describe the wave propagation of more general solutions and to study the stability/instability of the steady states, even when we extend the study to several space dimensions. A similar study is performed in the critical case that we call "pseudo-linear", i.e., when the operator is still nonlinear but has homogeneity one. With respect to the classical model and the "pseudo-linear" case, the travelling waves of the "slow" diffusion setting exhibit free boundaries. \\ Finally, as a complement of \cite{AA-JLV:art}, we study the asymptotic behaviour of more general solutions in the presence of a "heterozygote superior" reaction function and doubly nonlinear diffusion ("slow" and "pseudo-linear").Comment: 42 pages, 11 figures. Accepted version on Discrete Contin. Dyn. Sys

A moving point approach to model shallow ice sheets: a study case with radially-symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly and the waiting time behaviours to be modelled efficiently. A finite difference moving point scheme is derived and applied in a simplified context (continental radially-symmetrical shallow ice approximation). The scheme, which is inexpensive, is validated by comparing the results with moving-margin exact solutions and steady states. In both cases the scheme is able to track the position of the ice sheet margin with high precision

A moving-point approach to model shallow ice sheets: a study case with radially symmetrical ice sheets

Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically based moving-point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly. Our approach is also well suited to capture waiting-time behaviour efficiently. A finite-difference moving-point scheme is derived and applied in a simplified context (continental radially symmetrical shallow ice approximation). The scheme, which is inexpensive, is verified by comparing the results with steady states obtained from an analytic solution and with exact moving-margin transient solutions. In both cases the scheme is able to track the position of the ice sheet margin with high accuracy

Mathematical Theory of Water Waves

Water waves, that is waves on the surface of a fluid (or the interface between different fluids) are omnipresent phenomena. However, as Feynman wrote in his lecture, water waves that are easily seen by everyone, and which are usually used as an example of waves in elementary courses, are the worst possible example; they have all the complications that waves can have. These complications make mathematical investigations particularly challenging and the physics particularly rich. Indeed, expertise gained in modelling, mathematical analysis and numerical simulation of water waves can be expected to lead to progress in issues of high societal impact (renewable energies in marine environments, vorticity generation and wave breaking, macro-vortices and coastal erosion, ocean shipping and near-shore navigation, tsunamis and hurricane-generated waves, floating airports, ice-sea interactions, ferrofluids in high-technology applications, ...). The workshop was mostly devoted to rigorous mathematical theory for the exact hydrodynamic equations; numerical simulations, modelling and experimental issues were included insofar as they had an evident synergy effect

Mathematical Theory and Modelling in Atmosphere-Ocean Science

Mathematical theory and modelling in atmosphere-ocean science combines a broad range of advanced mathematical and numerical techniques and research directions. This includes the asymptotic analysis of multiscale systems, the deterministic and stochastic modelling of sub-grid-scale processes, and the numerical analysis of nonlinear PDEs over a broad range of spatial and temporal scales. This workshop brought together applied mathematicians and experts in the disciplinary fields of meteorology and oceanography for a wide-ranging exchange of ideas and results in this area with the aim of fostering fundamental interdisciplinary work in this important science area