## Abstract

A stationary Stokes problem with nonlinear rheology and with mixed no-slip and sliding basal boundary conditions is considered. The model describes the flow of ice in glaciers and ice sheets. A least-squares finite element method is developed and analyzed. The method does not require that the finite element spaces satisfy an inf-sup condition. Moreover, the usage of negative Sobolev norm in the least-squares functional allows for the use of standard piecewise polynomials spaces for both the velocity and pressure approximations. A Picard-type iterative method is used to linearize the Stokes problem. It is shown that the linearized least-squares functional is coercive and continuous in an appropriate solution space so the existence and uniqueness of a weak solution immediately follows as do optimal error estimates for finite element approximations. Numerical tests are provided to illustrate the theory.

Original language | English |
---|---|

Pages (from-to) | 2421-2431 |

Number of pages | 11 |

Journal | Computers and Mathematics with Applications |

Volume | 71 |

Issue number | 11 |

DOIs | |

Publication status | Published - 2016 Jun 1 |

### Bibliographical note

Funding Information:The second author’s research was supported by Republic of Korea NRF grant 2015 R1A5A1009350 . The third author’s research was supported by the United States Department of Energy Office of Science grant DE-SC0008273 .

## All Science Journal Classification (ASJC) codes

- Modelling and Simulation
- Computational Theory and Mathematics
- Computational Mathematics