Search

The effects of kinetic-energy preservation errors due to Runge-Kutta (RK) temporal integrators have been analyzed for the case of large-eddy simulations of incompressible turbulent channel flow. Simulations have been run using the open-source solver Xcompact3D with an implicit spectral vanishing viscosity model and a variety of temporal Runge-Kutta integrators. Explicit pseudo-symplectic schemes, with improved energy preservation properties, have been compared to standard RK methods. The results show a marked decrease in the temporal error for higher-order pseudo-symplectic methods; on the other hand, an analysis of the energy spectra indicates that the dissipation introduced by the commonly used three-stage RK scheme can lead to significant distortion of the energy distribution within the inertial range. A cost-vs-accuracy analysis suggests that pseudo-symplectic schemes could be used to attain results comparable to traditional methods at a reduced computational cost., Comment: 13 pages, 11 figures. F.C. and G.C. contributed equally to this work

This study proposes a novel spatial discretization procedure for the compressible Euler equations which guarantees entropy conservation at a discrete level when an arbitrary equation of state is assumed. The proposed method, based on a locally-conservative discretization, guarantees also the spatial conservation of mass, momentum, and total energy and is kinetic energy-preserving. In order to achieve the entropy-conservation property for an arbitrary non-ideal gas, a general strategy is adopted based on the manipulation of discrete balance equations through the imposition of global entropy conservation and the use of a summation by parts rule. The procedure, which is extended to an arbitrary order of accuracy, conducts to a general form of the internal-energy numerical flux which results in a nonlinear function of thermodynamic and dynamic variables and still admits the mass flux as a residual degree of freedom. The effectiveness of the novel entropy-conservative formulation is demonstrated through numerical tests making use of some of the most popular cubic equations of state., Comment: 16 pages, 3 figures

Employing physically-consistent numerical methods is an important step towards attaining robust and accurate numerical simulations. When addressing compressible flows, in addition to preserving kinetic energy at a discrete level, as done in the incompressible case, additional properties are sought after, such as the ability to preserve the equilibrium of pressure that can be found at contact interfaces. This paper investigates the general conditions of the spatial numerical discretizations to achieve the pressure equilibrium preserving property (PEP). Schemes from the literature are analyzed in this respect, and procedures to impart the PEP property to existing discretizations are proposed. Additionally, new PEP numerical schemes are introduced through minor modifications of classical ones. Numerical tests confirmed the theory hereby presented and showed that the modifications, beyond the enforcement of the PEP property, have a generally positive impact on the performances of the original schemes., Comment: 15 pages, 2 figures

This paper proposes a hierarchy of numerical fluxes for the compressible flow equations which are kinetic-energy and pressure equilibrium preserving and asymptotically entropy conservative, i.e., they are able to arbitrarily reduce the numerical error on entropy production due to the spatial discretization. The fluxes are based on the use of the harmonic mean for internal energy and only use algebraic operations, making them less computationally expensive than the entropy-conserving fluxes based on the logarithmic mean. The use of the geometric mean is also explored and identified to be well-suited to reduce errors on entropy evolution. Results of numerical tests confirmed the theoretical predictions and the entropy-conserving capabilities of a selection of schemes have been compared., Comment: 9 pages, 4 figures

Extremal dependence describes the strength of correlation between the largest observations of two variables. It is usually measured with symmetric dependence coefficients that do not depend on the order of the variables. In many cases, there is a natural asymmetry between extreme observations that can not be captured by such coefficients. An example for such asymmetry are large discharges at an upstream and a downstream stations on a river network: an extreme discharge at the upstream station will directly influence the discharge at the downstream station, but not vice versa. Simple measures for asymmetric dependence in extreme events have not yet been investigated. We propose the asymmetric tail Kendall's $\tau$ as a measure for extremal dependence that is sensitive to asymmetric behaviour in the largest observations. It essentially computes the classical Kendall's $\tau$ but conditioned on the extreme observations of one of the two variables. We show theoretical properties of this new coefficient and derive a formula to compute it for existing copula models. We further study its effectiveness and connections to causality in simulation experiments. We apply our methodology to a case study on river networks in the United Kingdom to illustrate the importance of measuring asymmetric extremal dependence in hydrology. Our results show that there is important structural information in the asymmetry that would have been missed by a symmetric measure. Our methodology is an easy but effective tool that can be applied in exploratory analysis for understanding the connections among variables and to detect possible asymmetric dependencies.

We analyze the conservation properties of various discretizations of the system of compressible Euler equations for shock-free flows, with special focus on the treatment of the energy equation and on the induced discrete equations for other thermodynamic quantities. The analysis is conducted both theoretically and numerically and considers two important factors characterizing the various formulations, namely the choice of the energy equation and the splitting used in the discretization of the convective terms. The energy equations analyzed are total and internal energy, total enthalpy, pressure, speed of sound and entropy. In all the cases examined the discretization of the convective terms is made with locally conservative and kinetic-energy preserving schemes. Some important relations between the various formulations are highlighted and the performances of the various schemes are assessed by considering two widely used test cases. Together with some popular formulations from the literature, also new and potentially useful ones are analyzed.

Several environmental phenomena can be described by different correlated variables that must be considered jointly in order to be more representative of the nature of these phenomena. For such events, identification of extremes is inappropriate if it is based on marginal analysis. Extremes have usually been linked to the notion of quantile, which is an important tool to analyze risk in the univariate setting. We propose to identify multivariate extremes and analyze environmental phenomena in terms of the directional multivariate quantile, which allows us to analyze the data considering all the variables implied in the phenomena, as well as look at the data in interesting directions that can better describe an environmental catastrophe. Since there are many references in the literature that propose extremes detection based on copula models, we also generalize the copula method by introducing the directional approach. Advantages and disadvantages of the non-parametric proposal that we introduce and the copula methods are provided in the paper. We show with simulated and real data sets how by considering the first principal component direction we can improve the visualization of extremes. Finally, two cases of study are analyzed: a synthetic case of flood risk at a dam (a 3-variable case), and a real case study of sea storms (a 5-variable case)., Comment: Article with supplementary material in the appendix

We show that the Gamma distribution is not an adequate fit for the probability density function of drop diameters using the Kolmogorov-Smirnov goodness of fit test. We propose a different parametrization of drop size distributions, which not depending by any particular functional form, is based on the adoption of standardized central moments. The first three standardized central moments are sufficient to characterize the distribution of drop diamters at the ground. These parameters together with the drop count form a 4-tuple which fully describe the variability of the drop size distributions. The Cartesian product of this 4-tuple of parameters is the rainfall phase space. Using disdrometer data from 10 different locations we identify invariant, not depending on location, properties of the rainfall phenomenon.

A long standing issue in Hydrology is the arbitrariness of the rain "event" definition. In this manuscript, we show that 1) the event definition resting on the occurrence of a minimum rainless period and the one resting a sequence of consecutive wet intervals are statistically equivalent. 2) In the case of stratiform rain, a non arbitrary definition of rain event is possible. The dynamical properties of stratiform rain indicate the range [1.5,4] h as the proper one for the choice of a minimum rainless period for Chilbolton, UK. 3) The intra event dynamical variability is "described" by an alternate sequence of quiescent and active phases., Comment: On the statistical definition of rain events