Abstract
In this work we propose a new numerical method for solving a thermomechanical coupled problem which arises as a mathematical model for the evolution of ice sheets profiles and the corresponding temperatures. From the accumulation–ablation ratio, the atmospheric temperature and the geothermic flux, the ice sheet profile is the solution of a moving boundary problem governed by a nonlinear convection–diffusion equation. Moreover, the mathematical model which governs the temperature distribution involves three nonlinearities: a reaction term due to the viscous dissipation, a Signorini boundary condition associated to the geothermic flux and an enthalpy term issued from the two-phase Stefan formulation of the polythermal regime. The problems are discretized in time with upwind characteristics schemes and in space with finite elements. The nonlinearities are solved either by fixed point methods or by duality techniques. Several numerical simulation examples involving real data sets issued from the Antartic ice sheet are shown.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: Mathematical Models and Methods in Applied Sciences
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
