No abstract

SUMMARYWe consider a Leray model with a deconvolution-based indicator function for the simulation of incompressible fluid flow at moderately large Reynolds number (in the range of few thousand) with underresolved meshes. For the implementation of the model, we adopt a three-step algorithm called evolve-filterrelax (EFR) that requires (i) the solution of a Navier-Stokes problem, (ii) the solution of a Stokes-like problem to filter the Navier-Stokes velocity field, and (iii) a final relaxation step. We take advantage of a reformulation of the EFR algorithm as an operator splitting method to analyze the impact of the filter on the final solution vs a direct simulation of the Navier-Stokes equations. In addition, we provide some direction for tuning the parameters involved in the model based on physical and numerical arguments. Our approach is validated against experimental data for fluid flow in an idealized medical device (consisting of a conical convergent, a narrow throat, and a sudden expansion, as recommended by the U.S. Food and Drug Administration). Numerical results are in good quantitative agreement with the measured axial components of the velocity and pressures for two different flow rates corresponding to turbulent regimes, even for meshes with a mesh size more than 40 times larger than the smallest turbulent scale. After several numerical experiments, we perform a preliminary sensitivity analysis of the computed solution to the parameters involved in the model.

Abstract. We introduce MPAS-Albany Land Ice (MALI) v6.0, a new variable-resolution land ice model that uses unstructured Voronoi grids on a plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework for developing variable-resolution Earth system model components and the Albany multi-physics code base for the solution of coupled systems of partial differential equations, which itself makes use of Trilinos solver libraries. MALI includes a three-dimensional first-order momentum balance solver (Blatter–Pattyn) by linking to the Albany-LI ice sheet velocity solver and an explicit shallow ice velocity solver. The evolution of ice geometry and tracers is handled through an explicit first-order horizontal advection scheme with vertical remapping. The evolution of ice temperature is treated using operator splitting of vertical diffusion and horizontal advection and can be configured to use either a temperature or enthalpy formulation. MALI includes a mass-conserving subglacial hydrology model that supports distributed and/or channelized drainage and can optionally be coupled to ice dynamics. Options for calving include “eigencalving”, which assumes that the calving rate is proportional to extensional strain rates. MALI is evaluated against commonly used exact solutions and community benchmark experiments and shows the expected accuracy. Results for the MISMIP3d benchmark experiments with MALI's Blatter–Pattyn solver fall between published results from Stokes and L1L2 models as expected. We use the model to simulate a semi-realistic Antarctic ice sheet problem following the initMIP protocol and using 2 km resolution in marine ice sheet regions. MALI is the glacier component of the Energy Exascale Earth System Model (E3SM) version 1, and we describe current and planned coupling to other E3SM components.

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40‐day prescribed‐sea‐surface‐temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non‐hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long‐standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.