Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 3.654 IF 3.654
  • IF 5-year<br/> value: 5.212 IF 5-year
  • SNIP value: 1.228 SNIP 1.228
  • IPP value: 4.021 IPP 4.021
  • SJR value: 2.365 SJR 2.365
  • h5-index value: 35 h5-index 35
Geosci. Model Dev., 6, 2023-2047, 2013
© Author(s) 2013. This work is distributed
under the Creative Commons Attribution 3.0 License.
22 Nov 2013
A hybrid Eulerian–Lagrangian numerical scheme for solving prognostic equations in fluid dynamics
E. Kaas1, B. Sørensen1, P. H. Lauritzen2, and A. B. Hansen3,4 1Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
2Climate and Global Dynamics Division, Boulder, Colorado, USA
3Aarhus University, Department of Environmental Science, Roskilde, Denmark
4National Institute of Water and Atmospheric Research, Lauder, New Zealand
Abstract. A new hybrid Eulerian–Lagrangian numerical scheme (HEL) for solving prognostic equations in fluid dynamics is proposed. The basic idea is to use an Eulerian as well as a fully Lagrangian representation of all prognostic variables.

The time step in Lagrangian space is obtained as a translation of irregularly spaced Lagrangian parcels along downstream trajectories. Tendencies due to other physical processes than advection are calculated in Eulerian space, interpolated, and added to the Lagrangian parcel values. A directionally biased mixing amongst neighboring Lagrangian parcels is introduced. The rate of mixing is proportional to the local deformation rate of the flow.

The time stepping in Eulerian representation is achieved in two steps: first a mass-conserving Eulerian or semi-Lagrangian scheme is used to obtain a provisional forecast. This forecast is then nudged towards target values defined from the irregularly spaced Lagrangian parcel values. The nudging procedure is defined in such a way that mass conservation and shape preservation is ensured in Eulerian space.

The HEL scheme has been designed to be accurate, multi-tracer efficient, mass conserving, and shape preserving. In Lagrangian space only physically based mixing takes place; i.e., the problem of artificial numerical mixing is avoided. This property is desirable in atmospheric chemical transport models since spurious numerical mixing can impact chemical concentrations severely.

The properties of HEL are here verified in two-dimensional tests. These include deformational passive transport on the sphere, and simulations with a semi-implicit shallow water model including topography.

Citation: Kaas, E., Sørensen, B., Lauritzen, P. H., and Hansen, A. B.: A hybrid Eulerian–Lagrangian numerical scheme for solving prognostic equations in fluid dynamics, Geosci. Model Dev., 6, 2023-2047, doi:10.5194/gmd-6-2023-2013, 2013.
Publications Copernicus