Multi-scale waves in sound-proof global simulations with EULAG

EULAG is a computational model for simulating flows across a wide range of scales and physical scenarios. A standard option employs an anelastic approximation to capture nonhydrostatic effects and simultaneously filter sound waves from the solution. In this study, we examine a localized gravity wave packet generated by instabilities in Held-Suarez climates. Although still simplified versus the Earth’s atmosphere, a rich set of planetary wave instabilities and ensuing radiated gravity waves can arise. Wave packets are observed that have lifetimes ≤ 2 days, are negligibly impacted by Coriolis force, and do not show the rotational effects of differential jet advection typical of inertia-gravity waves. Linear modal analysis shows that wavelength, period, and phase speed fit the dispersion equation to within a mean difference of ∼ 4%, suggesting an excellent fit. However, the group velocities match poorly even though a propagation of uncertainty analysis indicates that they should be predicted as well as the phase velocities. Theoretical arguments suggest the discrepancy is due to nonlinearity — a strong southerly flow leads to a critical surface forming to the southwest of the wave packet that prevents the expected propagation

A new pressure relaxation closure model for one-dimensional two-material Lagrangian hydrodynamics

We present a new model for closing a system of Lagrangian hydrodynamics equations for a two-material cell with a single velocity model. We describe a new approach that is motivated by earlier work of Delov and Sadchikov and of Goncharov and Yanilkin. Using a linearized Riemann problem to initialize volume fraction changes, we require that each material satisfy its own p dV equation, which breaks the overall energy balance in the mixed cell. To enforce this balance, we redistribute the energy discrepancy by assuming that the corresponding pressure change in each material is equal. This multiple-material model is packaged as part of a two-step time integration scheme. We compare results of our approach with other models and with corresponding pure-material calculations, on two-material test problems with ideal-gas or stiffened-gas equations of state

Volume Tracking of Interfaces Having Surface Tension in Two and Three Dimensions

. Solution algorithms are presented for tracking interfaces with piecewise linear (PLIC) volume-of-fluid (VOF) methods on Eulerian grids (structured and unstructured) in two and three dimensions. We review the theory of volume tracking methods, derive appropriate volume evolution equations, identify and present solutions to the basic geometric functions needed for interface reconstruction and volume fluxing, and provide algorithm templates for modern 2-D and 3-D PLIC VOF interface tracking methods. We discuss some key issues for PLIC VOF methods, namely the method used for time integration of fluid volumes (operator splitting, unsplit, Runge-Kutta, etc.) and the estimation of interface normals. We also present our latest developments in the continuum surface force (CSF) model for surface tension, namely extension to 3-D and variable surface tension effects. We identify and focus on CSF model issues that become especially critical on fine meshes with high density ratio interfacial flows..