Journal Metrics

  • IF value: 6.086 IF 6.086
  • IF 5-year<br/> value: 6.174 IF 5-year
    6.174
  • SNIP value: 1.812 SNIP 1.812
  • IPP value: 5.140 IPP 5.140
  • SJR value: 3.969 SJR 3.969
  • h5-index value: 29 h5-index 29
Geosci. Model Dev., 6, 2023-2047, 2013
www.geosci-model-dev.net/6/2023/2013/
doi:10.5194/gmd-6-2023-2013
© Author(s) 2013. This work is distributed
under the Creative Commons Attribution 3.0 License.
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.
 
Search GMD
Final Revised Paper
PDF XML
Citation
Discussion Paper
Share