Brinkman volume penalization for bathymetry in three-dimensional ocean models

Accurate and stable implementation of bathymetry boundary conditions remains a challenging problem. The dynamics of ocean flow often depend sensitively on satisfying bathymetry boundary conditions and correctly representing their complex geometry. Generalized (e.g. sigma) terrain-following coordinat...

Full description

Bibliographic Details
Main Authors: Debreu, L., Kevlahan, N. K. R., /Marchesiello, Patrick
Format: Text
Language:English
Published: 2020
Subjects:
Online Access:https://www.documentation.ird.fr/hor/fdi:010077471
id ftird:oai:ird.fr:fdi:010077471
record_format openpolar
spelling ftird:oai:ird.fr:fdi:010077471 2024-09-15T18:12:40+00:00 Brinkman volume penalization for bathymetry in three-dimensional ocean models Debreu, L. Kevlahan, N. K. R. /Marchesiello, Patrick 2020 https://www.documentation.ird.fr/hor/fdi:010077471 EN eng https://www.documentation.ird.fr/hor/fdi:010077471 oai:ird.fr:fdi:010077471 Debreu L., Kevlahan N. K. R., Marchesiello Patrick. Brinkman volume penalization for bathymetry in three-dimensional ocean models. 2020, 145, 101530 [13 p.] Algorithms Bathymetry Bottom topography Computation Ocean modelling Penalization text 2020 ftird 2024-08-15T05:57:41Z Accurate and stable implementation of bathymetry boundary conditions remains a challenging problem. The dynamics of ocean flow often depend sensitively on satisfying bathymetry boundary conditions and correctly representing their complex geometry. Generalized (e.g. sigma) terrain-following coordinates are often used in ocean models, but they require smoothing the bathymetry to reduce pressure gradient errors (Mellor a al., 1994). Geopotential z-coordinates are a common alternative that avoid pressure gradient and numerical diapycnal diffusion errors, but they generate spurious flow due to their "staircase" geometry. We introduce a new Brinkman volume penalization to approximate the no-slip boundary condition and complex geometry of bathymetry in ocean models. This approach corrects the staircase effect of z-coordinates, does not introduce any new stability constraints on the geometry of the bathymetry and is easy to implement in an existing ocean model. The porosity parameter allows modelling subgrid scale details of the geometry. We illustrate the penalization and confirm its accuracy by applying it to three standard test flows: upwelling over a sloping bottom, resting state over a seamount and internal fides over highly peaked bathymetry features. In future work we will explore applying the penalization to more realistic bathymetry configurations, and moving boundaries such as melting/freezing ice shelves. Text Ice Shelves IRD (Institute de recherche pour le développement): Horizon
institution Open Polar
collection IRD (Institute de recherche pour le développement): Horizon
op_collection_id ftird
language English
topic Algorithms
Bathymetry
Bottom topography
Computation
Ocean modelling
Penalization
spellingShingle Algorithms
Bathymetry
Bottom topography
Computation
Ocean modelling
Penalization
Debreu, L.
Kevlahan, N. K. R.
/Marchesiello, Patrick
Brinkman volume penalization for bathymetry in three-dimensional ocean models
topic_facet Algorithms
Bathymetry
Bottom topography
Computation
Ocean modelling
Penalization
description Accurate and stable implementation of bathymetry boundary conditions remains a challenging problem. The dynamics of ocean flow often depend sensitively on satisfying bathymetry boundary conditions and correctly representing their complex geometry. Generalized (e.g. sigma) terrain-following coordinates are often used in ocean models, but they require smoothing the bathymetry to reduce pressure gradient errors (Mellor a al., 1994). Geopotential z-coordinates are a common alternative that avoid pressure gradient and numerical diapycnal diffusion errors, but they generate spurious flow due to their "staircase" geometry. We introduce a new Brinkman volume penalization to approximate the no-slip boundary condition and complex geometry of bathymetry in ocean models. This approach corrects the staircase effect of z-coordinates, does not introduce any new stability constraints on the geometry of the bathymetry and is easy to implement in an existing ocean model. The porosity parameter allows modelling subgrid scale details of the geometry. We illustrate the penalization and confirm its accuracy by applying it to three standard test flows: upwelling over a sloping bottom, resting state over a seamount and internal fides over highly peaked bathymetry features. In future work we will explore applying the penalization to more realistic bathymetry configurations, and moving boundaries such as melting/freezing ice shelves.
format Text
author Debreu, L.
Kevlahan, N. K. R.
/Marchesiello, Patrick
author_facet Debreu, L.
Kevlahan, N. K. R.
/Marchesiello, Patrick
author_sort Debreu, L.
title Brinkman volume penalization for bathymetry in three-dimensional ocean models
title_short Brinkman volume penalization for bathymetry in three-dimensional ocean models
title_full Brinkman volume penalization for bathymetry in three-dimensional ocean models
title_fullStr Brinkman volume penalization for bathymetry in three-dimensional ocean models
title_full_unstemmed Brinkman volume penalization for bathymetry in three-dimensional ocean models
title_sort brinkman volume penalization for bathymetry in three-dimensional ocean models
publishDate 2020
url https://www.documentation.ird.fr/hor/fdi:010077471
genre Ice Shelves
genre_facet Ice Shelves
op_relation https://www.documentation.ird.fr/hor/fdi:010077471
oai:ird.fr:fdi:010077471
Debreu L., Kevlahan N. K. R., Marchesiello Patrick. Brinkman volume penalization for bathymetry in three-dimensional ocean models. 2020, 145, 101530 [13 p.]
_version_ 1810450251818991616