Search for: All records

Numerical weather prediction is pushing the envelope of grid resolution at local and global scales alike. Aiming to model topography with higher precision, a handful of articles introduced unstructured vertical grids and tested them for dry atmospheres. The next step toward effective high‐resolution unstructured grids for atmospheric modeling requires that also microphysics is independent of any vertical columns, in contrast to what is ubiquitous across operational and research models. In this paper, we present a non‐column based continuous and discontinuous spectral element implementation of Kessler's microphysics with warm rain. We test the proposed algorithm against standard three‐dimensional benchmarks for precipitating clouds and show that the results are comparable with those presented in the literature across all of the tested effective resolutions. While presented for both continuous and discontinuous spectral elements in this paper, the method that we propose can be adapted to any numerical method used in other codes, as long as the code can already handle vertically unstructured grids.

 

The paper discusses a reuse of matrix factorization as a building block in the Augmented Lagrangian (AL) and modified AL preconditioners for nonsymmetric saddle point linear algebraic systems. The strategy is applied to solve two‐dimensional incompressible fluid problems with efficiency rates independent of the Reynolds number. The solver is then tested to simulate motion of a surface fluid, an example of a two‐dimensional flow motivated by an interest in lateral fluidity of inextensible viscous membranes. Numerical examples include the Kelvin–Helmholtz instability problem posed on the sphere and on the torus. Some new eigenvalue estimates for the AL preconditioner are derived.

 

Within OpenFOAM, we develop a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe non-hydrostatic atmospheric flow. For the stabilization of the Euler equations and to capture sub-grid processes, we consider two Large Eddy Simulation models: the classical Smagorinsky model and the one equation eddy-viscosity model. To achieve high computational efficiency, our solver uses a splitting scheme that decouples the computation of each variable. The numerical results obtained with our solver are validated against numerical data available in the literature for two classical benchmarks: the rising thermal bubble and the density current. Through qualitative and quantitative comparisons, we show that our approach is accurate. This paper is meant to lay the foundation for a new open-source package specifically created for the quick assessment of new computational approaches for the simulation of atmospheric flows at the mesoscale level. 

Abstract This work introduces a novel approach for data-driven model reduction of time-dependent parametric partial differential equations. Using a multi-step procedure consisting of proper orthogonal decomposition, dynamic mode decomposition, and manifold interpolation, the proposed approach allows to accurately recover field solutions from a few large-scale simulations. Numerical experiments for the Rayleigh-Bénard cavity problem show the effectiveness of such multi-step procedure in two parametric regimes, i.e., medium and high Grashof number. The latter regime is particularly challenging as it nears the onset of turbulent and chaotic behavior. A major advantage of the proposed method in the context of time-periodic solutions is the ability to recover frequencies that are not present in the sampled data. 

We construct and analyze a CutFEM discretization for the Stokes problem based on the Scott–Vogelius pair. The discrete piecewise polynomial spaces are defined on macro-element triangulations which are not fitted to the smooth physical domain. Boundary conditions are imposed via penalization through the help of a Nitsche-type discretization, whereas stability with respect to small and anisotropic cuts of the bulk elements is ensured by adding local ghost penalty stabilization terms. We show stability of the scheme as well as a divergence–free property of the discrete velocity outside an O ( h ) neighborhood of the boundary. To mitigate the error caused by the violation of the divergence–free condition, we introduce local grad–div stabilization. The error analysis shows that the grad–div parameter can scale like O ( h −1 ), allowing a rather heavy penalty for the violation of mass conservation, while still ensuring optimal order error estimates. 

The paper considers a system of equations that models a lateral flow of a Boussinesq–Scriven fluid on a passively evolving surface embedded in [Formula: see text]. For the resulting Navier–Stokes type system, posed on a smooth closed time-dependent surface, we introduce a weak formulation in terms of functional spaces on a space-time manifold defined by the surface evolution. The weak formulation is shown to be well-posed for any finite final time and without smallness conditions on data. We further extend an unfitted finite element method, known as TraceFEM, to compute solutions to the fluid system. Convergence of the method is demonstrated numerically. In another series of experiments we visualize lateral flows induced by smooth deformations of a material surface.