A potential energy conserving finite element method for turbulent variable density flow: application to glacier-fjord circulation

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Lundgren, Lukas, Helanow, Christian, Wiskandt, Jonathan, Koszalka, Inga Monika, Ahlkrona, Josefin
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916521441230848
author Lundgren, Lukas
Helanow, Christian
Wiskandt, Jonathan
Koszalka, Inga Monika
Ahlkrona, Josefin
author_facet Lundgren, Lukas
Helanow, Christian
Wiskandt, Jonathan
Koszalka, Inga Monika
Ahlkrona, Josefin
contents We introduce a continuous Galerkin finite element discretization of the non-hydrostatic Boussinesq approximation of the Navier-Stokes equations, suitable for various applications such as coastal ocean dynamics and ice-ocean interactions, among others. In particular, we introduce a consistent modification of the gravity force term which enhances conservation properties for Galerkin methods without strictly enforcing the divergence-free condition. We show that this modification results in a sharp energy estimate, including both kinetic and potential energy. Additionally, we propose a new, symmetric, tensor-based viscosity operator that is especially suitable for modeling turbulence in stratified flow. The viscosity coefficients are constructed using a residual-based shock-capturing method and the method conserves angular momentum and dissipates kinetic energy. We validate our proposed method through numerical tests and use it to model the ocean circulation and basal melting beneath the ice tongue of the Ryder Glacier and the adjacent Sherard Osborn fjord in two dimensions on a fully unstructured mesh. Our results compare favorably with a standard numerical ocean model, showing better resolved turbulent flow features and reduced artificial diffusion.
format Preprint
id arxiv_https___arxiv_org_abs_2409_00972
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A potential energy conserving finite element method for turbulent variable density flow: application to glacier-fjord circulation
Lundgren, Lukas
Helanow, Christian
Wiskandt, Jonathan
Koszalka, Inga Monika
Ahlkrona, Josefin
Numerical Analysis
34A45
G.1.8
We introduce a continuous Galerkin finite element discretization of the non-hydrostatic Boussinesq approximation of the Navier-Stokes equations, suitable for various applications such as coastal ocean dynamics and ice-ocean interactions, among others. In particular, we introduce a consistent modification of the gravity force term which enhances conservation properties for Galerkin methods without strictly enforcing the divergence-free condition. We show that this modification results in a sharp energy estimate, including both kinetic and potential energy. Additionally, we propose a new, symmetric, tensor-based viscosity operator that is especially suitable for modeling turbulence in stratified flow. The viscosity coefficients are constructed using a residual-based shock-capturing method and the method conserves angular momentum and dissipates kinetic energy. We validate our proposed method through numerical tests and use it to model the ocean circulation and basal melting beneath the ice tongue of the Ryder Glacier and the adjacent Sherard Osborn fjord in two dimensions on a fully unstructured mesh. Our results compare favorably with a standard numerical ocean model, showing better resolved turbulent flow features and reduced artificial diffusion.
title A potential energy conserving finite element method for turbulent variable density flow: application to glacier-fjord circulation
topic Numerical Analysis
34A45
G.1.8
url https://arxiv.org/abs/2409.00972