Spatiotemporal variations in soil cultivable mycobiota at the Arava desert (Israel) along latitudinal and elevational gradients

Isabella Grishkan; Isabella Grishkan

doi:10.3934/microbiol.2018.3.502

AIMS Microbiology

2018, Volume 4, Issue 3: 502-521. doi: 10.3934/microbiol.2018.3.502

Previous Article Next Article

Research article

Spatiotemporal variations in soil cultivable mycobiota at the Arava desert (Israel) along latitudinal and elevational gradients

Isabella Grishkan ^,

Institute of Evolution, University of Haifa, 199 Aba Khoushy Ave., Mount Carmel, Haifa 3498838, Israel

Received: 23 May 2018 Accepted: 21 June 2018 Published: 28 June 2018

Regional, local, and seasonal distribution of soil culturable microfungi in the Arava Valley, Israel, was examined along altitudinal and latitudinal gradients. A total of 198 species from 86 genera were isolated using the soil dilution plate method. Melanin-containing species with large multi-cellular spores dominated the majority of microfungal communities, while species with picnidial fruit bodies mostly prevailed in the northern part of the Arava Valley located at 190 m below sea level. Aspergilli (mainly Aspergillus fumigatus) and teleomorphic ascomycetes comprised the basic part of thermotolerant mycobiota obtained at 37 °C. The soil at the northern part of the desert held the highest number of microfungal isolates and, at the same time, was characterized by significantly lower species richness. The open sun-exposed localities harbored a significantly higher number of species than the localities under shrub canopies. Isolate density displayed the opposite trend and was significantly lower in the open than in shrub localities. The mycobiota characteristics such as species composition, contribution of major groupings to mycobiota structure, diversity level, and isolate density showed significant correlations with measured edaphic parameters—organic matter content, water content, pH, and especially, with electrical conductivity. Among the environmental aspects, locality position along altitudinal and latitudinal gradients accompanied by locality type (open sun-exposed or under shrubs), strongly influenced the community’s characteristics, thus demonstrating the effect of the unique altitudinal position of the northern part of the Arava Valley as well as the ability of microfungal communities to be sensitive to the microscale environmental variability.

Keywords:

Citation: Isabella Grishkan. Spatiotemporal variations in soil cultivable mycobiota at the Arava desert (Israel) along latitudinal and elevational gradients[J]. AIMS Microbiology, 2018, 4(3): 502-521. doi: 10.3934/microbiol.2018.3.502

Related Papers:

[1]	Stefania Fresca, Federico Fatone, Andrea Manzoni . Long-time prediction of nonlinear parametrized dynamical systems by deep learning-based reduced order models. Mathematics in Engineering, 2023, 5(6): 1-36. doi: 10.3934/mine.2023096
[2]	Michele Dolce, Ricardo Grande . On the convergence rates of discrete solutions to the Wave Kinetic Equation. Mathematics in Engineering, 2024, 6(4): 536-558. doi: 10.3934/mine.2024022
[3]	Jason Howell, Katelynn Huneycutt, Justin T. Webster, Spencer Wilder . A thorough look at the (in)stability of piston-theoretic beams. Mathematics in Engineering, 2019, 1(3): 614-647. doi: 10.3934/mine.2019.3.614
[4]	Boubacar Fall, Filippo Santambrogio, Diaraf Seck . Shape derivative for obstacles in crowd motion. Mathematics in Engineering, 2022, 4(2): 1-16. doi: 10.3934/mine.2022012
[5]	Ludovica Cicci, Stefania Fresca, Stefano Pagani, Andrea Manzoni, Alfio Quarteroni . Projection-based reduced order models for parameterized nonlinear time-dependent problems arising in cardiac mechanics. Mathematics in Engineering, 2023, 5(2): 1-38. doi: 10.3934/mine.2023026
[6]	Simon Lemaire, Julien Moatti . Structure preservation in high-order hybrid discretisations of potential-driven advection-diffusion: linear and nonlinear approaches. Mathematics in Engineering, 2024, 6(1): 100-136. doi: 10.3934/mine.2024005
[7]	Francesca Tedeschi, Giulio G. Giusteri, Leonid Yelash, Mária Lukáčová-Medvid'ová . A multi-scale method for complex flows of non-Newtonian fluids. Mathematics in Engineering, 2022, 4(6): 1-22. doi: 10.3934/mine.2022050
[8]	Alessia Nota, Juan J. L. Velázquez . Homoenergetic solutions of the Boltzmann equation: the case of simple-shear deformations. Mathematics in Engineering, 2023, 5(1): 1-25. doi: 10.3934/mine.2023019
[9]	François Alouges, Aline Lefebvre-Lepot, Jessie Levillain . A limiting model for a low Reynolds number swimmer with $ N $ passive elastic arms. Mathematics in Engineering, 2023, 5(5): 1-20. doi: 10.3934/mine.2023087
[10]	Yangyang Qiao, Qing Li, Steinar Evje . On the numerical discretization of a tumor progression model driven by competing migration mechanisms. Mathematics in Engineering, 2022, 4(6): 1-24. doi: 10.3934/mine.2022046

Abstract

1. Introduction

In this note we combine some results on the long-time averaging of fluid equations with the recent techniques for reduced order model (ROM) development. In this preliminary work we start proving some analytical results that characterize the time-averaged effect of the exchange of energy between various modes, both in the case of the computable decomposition made with proper orthogonal decomposition (POD) type basis functions and with the abstract basis made with eigenfunctions. We will show that the results obtainable with a generic (but computable) basis are less precise than those obtainable with the abstract spectral basis, the difference coming from the lack of orthogonality of the gradients of the POD basis functions.

We then provide a few numerical examples. Concerning the analytical results we will prove partial results for the energy exchange between large and small scales, showing the difference between the use of a spectral type basis and a POD one. In particular, we are interested in results connected to the statistical equilibrium problem, which can be deduced in a computable way by a long-time averaging of the solutions. We want to investigate the possible forward and backward average transfer of energy. The properties of a turbulent flow are computable (and relevant) only in an average sense. In this respect, we want to follow the classical approach dating back to Stokes, Reynolds, and Prandtl of considering long-time averages of the solution as the key quantity to be computed or observed. Therefore, we do not need to consider statistical averaging and link it with time averaging by means of (unproved) ergodic hypotheses.

To introduce the problem that we will consider, we recall that a Newtonian incompressible flow (with constant density) can be described by the Navier-Stokes equations (NSE in the sequel)

$\begin{equation} \left\{\begin{aligned} \partial_t {\bf u} - \nu\Delta {\bf u} +({\bf u} \cdot \nabla)\, {\bf u} + \nabla p & = {\bf f} \qquad \mbox{in }\Omega\times (0, T), \\ \nabla \cdot {\bf u}& = 0 \qquad \mbox{in }\Omega\times (0, T), \\ {\bf u}& = \mathbf{0} \qquad \mbox{on }\Gamma\times (0, T) , \\ {\bf u}(\cdot, 0)& = {\bf u}(0) \qquad \mbox{in }\Omega, \end{aligned} \right. \end{equation}$ $

(1.1)

with Dirichlet conditions when the motion takes place in a smooth and bounded domain $ \Omega\subset \mathbb{R}^3 $ with solid walls $ \Gamma: = \partial\Omega $. The unknowns are the velocity field $ {\bf u} $ and the scalar pressure $ p $, while the positive constant $ \nu > 0 $ is the kinematic viscosity. The key parameter to detect the nature of the problem is the non-dimensional Reynolds number, which is defined as

$\begin{equation*} Re = \frac{U L}{\nu} , \end{equation*}$ $

where $ U $ and $ L $ are a characteristic velocity and length of the problem. In realistic problems, the Reynolds number can be extremely large (in many cases of the order of $ 10^{8} $, but up to the order of $ 10^{20} $ in certain geophysical problems). For simplicity in the notation, we use the viscosity as a control parameter and assume that it is very small. Hence, the effect of the regularization (similar to the diffusion in heat transfer) due to the Laplacian is negligible and the behavior of solutions is really turbulent and rather close to the motion of ideal fluids. Due to the well-known difficulties in performing direct numerical simulations (DNS), it is nowadays a well-established technique that of trying to reduce the computational efforts by simulating only the largest scales, which, for limited computational and experimental resources, are the only ones computable and observable. In this framework, the large eddy simulation (LES) methods, which emerged in the last 30 years, are among the most popular, and they found a very relevant role within both theorists' and practitioners' communities. For recent LES reviews, see, for instance, the monographs ^[2,5,27,39].

The LES methods are, in many cases, very well developed and both theoretically and computationally appealing, especially for problems without boundaries, but many difficulties and open problems arise when facing solid boundaries. In most cases the design of efficient LES methods is guided by deep properties of the solutions, as emerging from deep theorems of mathematical analysis. Furthermore, the ultimate goal of having a golden standard is far from being obtained, and large families of methods (wave-number asymptotics, differential filters, $ \alpha $-models) attracted the interest of different communities, spanning from the pure mathematicians, to the applied geophysicist and mechanical engineers, to the computational practitioners. For these reasons, we believe that it is important to have some well-defined and clearly stated guidelines, that can be adapted to different problems. In this way the methods can be improved with insight not only from experts in modeling, but also from mathematicians, physicists, and computational scientists.

In this respect, we point out that very recently the use of other (more flexible and computationally simpler) ways of finding approximate systems has become very popular. The LES methods itself can be specialized or even glued with other ways of determining much smaller approximate systems, which are computable in a very efficient way. For instance, reduced order modeling is increasingly becoming an accepted paradigm, in which applications to fluids are still being developed ^{[15,25,35,36,38]}.

The basic ansatz at the basis of the use of these models is the approximation of the velocity by a truncation of the series

$\begin{equation*} {\bf u} = \sum\limits_{k = 1}^\infty {\bf u}_k {\bf w}_k, \end{equation*}$ $

where $ \{{\bf w}_k\}_{k\in \mathbb{N}} $ is a basis constructed by using the POD, not necessarily made with eigenfunctions of the Stokes operator, and the coefficients of the $ L^2 $-projections are evaluated as follows

$\begin{equation*} {\bf u}_k = \frac{\int_{\Omega}{\bf u}\cdot{\bf w}_k\, dx}{\int_\Omega|{\bf w}_k|^2\, dx}. \end{equation*}$ $

The appealing property of this approach is that the choice of the basis is adapted to and determined by the problem itself. Generally the basis is chosen after a preliminary numerical computation, hence it contains at least the basic features of the solution and geometry of the problem to be studied. The other basic fact is that the kinetic energy of the problem is the key quantity under consideration; in fact the $ L^2 $-projection is generally used to determine the approximate velocity and also the energy content in the basis construction. To determine the number $ r\in \mathbb{N} $ such that the solution is projected over the space generated by the (orthogonal) functions $ V_r: = \text{span}\{{\bf w}_1, \dots, {\bf w}_r\} $, generally it is assumed that the projection of the solution over $ V_r $ contains a large percentage (say 80$ \% $) of the total kinetic energy of the underlying system.

It turns out that a basis associated with the problem at hand can greatly improve the effectiveness of the ROM. Its proper choice can be of great interest in applications to fluid flows ^[41,42,44]. The present paper combines mathematical results on the long-time behavior of fluid flows (especially in the case of statistical equilibrium) with reduced order modeling. The main goal of this approach is that it provides a mathematical description of both the long-time averages of ROMs and the energy exchange between ROM modes. Furthermore, preliminary numerical results that support the theoretical developments are presented. Specifically, we are extending to the POD setting the results for statistical solutions by Foias et al. ^[11,12] and those more recently obtained for time-averages in ^[3,28]. In this respect, the main theoretical results of this paper, stated in Theorems 4.2 and 4.3 below, can be viewed as a spectral version of the results of ^[3,28]. These results show that low frequency modes yield a Reynolds stress that is dissipative in the mean, its total spatial mean work being larger than the long time average of the dissipation of the fluctuations, which is consistent with observations and results in ^[3,11]. However, the analysis shows that the triad interaction between high and low frequency modes yields an additional non positive term in the budget between the Reynolds stress of high modes and the corresponding mean dissipation. This term may be non dissipative and may permit an inverse energy cascade, which is not in contradiction with the fact that the total Reynolds stress is dissipative in mean.

The paper is organized as follows: In Section 2, we outline the general framework for ROMs of fluid flows, and we display the exchange of energy between large scales and small scales for two ROM bases: the POD and the Stokes eigenfunctions. In Section 3, we present some preliminaries on long-time averages. In Section 4, we prove the main theoretical results for the average transfer of energy for ROMs constructed with the POD and the Stokes eigenfunctions. In Section 5, we perform a preliminary numerical study in which we investigate the theoretical results in the numerical simulation of the one-dimensional Burgers equation. Finally, in Section 6, we draw conclusions and outline future research directions.

2. Reduced order modeling

As outlined in the introduction, one key quantity in the pure and applied analysis of the Navier-Stokes equations is the kinetic energy, since it is a meaningful physical quantity and the analysis of its budget is at the basis of abstract existence results for weak solutions (cf. Constantin and Foias ^[6]) and also of the conventional turbulence theories of Kolmogorov ^[19]. It is well-known that after testing the NSE (1.1) by $ {\bf u} $ and integrating over the space-time, one (formally) obtains the global energy balance

$\begin{equation*} \frac{1}{2}\|{\bf u}(t)\|^2+\nu\int_0^t\|\nabla {\bf u}(s)\|^2\, ds = \frac{1}{2}\|{\bf u}(0)\|^2+\int_0^t\int_\Omega {\bf f}\cdot{\bf u}\, dx ds. \end{equation*}$ $

At present, it is possible to prove that the above balance holds true with the sign of "less or equal" for the class of weak (or turbulent) solutions for which there are results of existence globally in time, without restrictions on the viscosity and on the size of square summable initial velocity $ {\bf u}(0) $ and external force $ {\bf f} $. It is of fundamental importance in many problems in pure mathematics to understand under which hypotheses the equality holds true. We are now focusing on the "global energy" which is an averaged quantity, since it is the integral of the square modulus of the velocity over the entire domain. We also point out that at the other extreme one can deduce, without the integration over the domain, the point-wise relation

$\begin{equation*} \partial_t\frac{|{\bf u}|^2}{2}+|\nabla {\bf u}|^2-\Delta\frac{|{\bf u}|^2}{2}+\text{div }\left(\frac{{\bf u}|{\bf u}|^2}{2}+p\, {\bf u}\right) = {\bf f}\cdot{\bf u}. \end{equation*}$ $

In between there is the so called "local energy" which can be obtained by multiplying the NSE by $ {\bf u}\, \phi $, where $ \phi $ is a bump function, before integrating in the space-time variables. The goal is to show that

$\begin{equation*} \int_0^T\int_{\Omega}|\nabla {\bf u}|^2\phi\, dx dt = \int_{0}^{T}\int_{\Omega}\left[\frac{|{\bf u}|^2}{2}\left(\partial_t\phi+\Delta\phi\right) +\left(\frac{|{\bf u}|^2}{2}+p\right){\bf u}\cdot\nabla \phi+{\bf f}\cdot {\bf u}\, \phi\right]\, dx dt \end{equation*}$ $

holds (at least with the inequality sign) for all smooth scalar functions $ \phi\in C^\infty_0((0, T)\times\Omega) $ such that $ \phi\geq0 $. The validity of such an inequality is one of the requirements to prove partial regularity results, but it is also one of the conditions to be satisfied by weak solutions constructed by numerical or LES methods. In this respect, see Guermond, Oden and Prudhomme ^[14] and ^[4].

In this paper we study the global energy in the perspective that it can be reconstructed in a computable way or it can be well approximated by the POD basis functions $ \{{\bf w}_k\} $.

The fact that the functions $ \{{\bf w}_k\} $ can be constructed to be orthonormal with respect to the scalar product in $ (L^2(\Omega)^{3}, \, \|\cdot\|) $ allows us to evaluate the kinetic energy easily by the following numerical series

$\begin{equation*} \mathrm{E}({\bf u}) = \frac{1}{2}\sum\limits_{k = 1}^\infty\|{\bf u}_k\|^2. \end{equation*}$ $

Since we are going to use only a reduced number of ROM modes, it is relevant to consider the energy contained in functions described by a restricted set of indices. Hence, following the notation in ^[12], if we define

$\begin{equation*} {\bf u}_{m', m''}: = \sum\limits_{k = m'}^{m''}{\bf u}_k{\bf w}_k, \end{equation*}$ $

then the kinetic energy content of $ {\bf u}_{m', m''} $ is simply evaluated as follows

$\begin{equation*} \mathrm{E}({\bf u}_{m', m''}) = \frac{1}{2}\sum\limits_{k = m'}^{m''}\|{\bf u}_k\|^2. \end{equation*}$ $

We want to investigate the energy transfer between the various modes, together with averaged long-time behavior associated with this splitting.

We are going to adapt well-known studies on the decomposition in small and large eddies. This would be the case if the functions $ {\bf w}_k $ are chosen to be the eigenfunctions of the Stokes operator, hence associated with large and small frequencies. In our setting the basis is determined by the solution of a simplified problem, which can be treated computationally, and the basis functions are orthonormal in $ L^2(\Omega) $, but we cannot expect that their gradients are also orthogonal.

For the NSE, the standard ROM is constructed as follows:

(ⅰ) choose modes $ \{ {\bf w}_1, \ldots, {\bf w}_d\} $, which represent the recurrent spatial structures of the given flow;

(ⅱ) choose the dominant modes $ \{{\bf w}_1, \ldots, {\bf w}_m\} $, with $ m \leq d $, as basis functions for the ROM;

(ⅲ) use a Galerkin truncation $ {\bf u}_m = \sum_{j = 1}^{m} a_j \, {\bf w}_j $;

(ⅳ) replace $ {\bf u} $ with $ {\bf u}_m $ in the NSE;

(ⅴ) use a Galerkin projection of NSE ($ {\bf u}_m $) onto the ROM space $ V_m: = \mbox{span} \{{\bf w}_1, \ldots, {\bf w}_m\} $ to obtain a low-dimensional dynamical system, which represents the ROM:

$\begin{equation*} \dot{a} = A\, a+a^{\top}\, B\, a\, , \end{equation*}$ $

where $ a $ is the vector of unknown ROM coefficients and $ A, B $ are ROM operators;

(ⅵ) in an offline stage, compute the ROM operators;

(ⅶ) in an online stage, repeatedly use the ROM (for various parameter settings and/or longer time intervals).

Hence, there is a very natural splitting of the velocity field $ {\bf u} $ into two components, the part coherent with the basis expansion associated with the more energetic modes ($ {\bf y} $), and the remainder ($ {\bf z} $). This can be formalized as follows:

$\begin{equation*} {\bf u} = {\bf y}+{\bf z}, \end{equation*}$ $

where

$\begin{equation} {\bf y} = \sum\limits_{k = 1}^{m}{\bf u}_k{\bf w}_k = {P}_m{\bf u}\qquad\text{and}\qquad {\bf z} = \sum\limits_{k = m+1}^{+\infty}{\bf u}_k{\bf w}_k = ({I}-{P}_m)\, {\bf u} = :{Q}_m\, {\bf u} . \end{equation}$ $

(2.1)

In (2.1), $ m\in \mathbb{N} $ can be selected computationally (based, e.g., on the relative kinetic energy content in the first $ m $ POD-modes, but other choices relative to the enstrophy are possible) in order to have a significant amount of the energy content of the flow in $ {\bf y} $. Furthermore, $ {P}_m $ is the projection operator over the subspace $ V_m $ spanned by the first $ m $-functions $ \{{\bf w}_k\}_{k = 1, \dots, m} $.

We observe that we are considering the functions $ \{{\bf w}_k\}_{k} $ as divergence-free, at least in a weak and/or approximate sense. Even if generally they are not "exactly divergence-free", numerically we can consider that they have vanishing divergence, by assuming that the solution of the problems used to construct the basis is accurate enough to have negligible divergence. Hence, in the computations that will follow, we will drop the pressure terms by a standard Leray projection. It will be nevertheless interesting to extend our study to bases that are not divergence-free, e.g., those derived by the combination of ROM with artificial compressibility methods ^[8,9,13].

In addition, we consider the external force as stationary, that is $ {\bf f} = {\bf f}(\mathbf{x})\in L^2(\Omega)^{3} $ and we look for conditions holding at statistical equilibrium. Our purpose is to determine –if possible– the long-time behavior of $ {\bf y} $ and to analyze the energy budget between low and high modes in the orthogonal decomposition determined by the functions $ {\bf w}_k $.

As usual in many problems fluid mechanics, we use the Hilbert space functional setting with

$\begin{equation*} \begin{aligned} & {\cal V} = \{ \boldsymbol{\varphi} \in {\cal D}(\Omega)^3, \, \, \nabla \cdot \boldsymbol{\varphi} = 0 \}, \\ &H = \left\{{\bf u}\in L^2(\Omega)^3, \ \nabla\cdot {\bf u} = 0, \ {\bf u}\cdot\mathbf{n} = 0 \text{ on }\Gamma\right\}, \\ &V = \left\{{\bf u}\in H^1_{0}(\Omega)^3, \, \nabla\cdot {\bf u} = 0\right\}, \end{aligned} \end{equation*}$ $

where $ \mathbf{n} $ denotes the outward normal unit vector. Moreover, $ V' $ is the topological dual space of $ V $. We will also denote by $ < \cdot, \cdot > $ the duality pairing between $ V $ and $ V' $. We recall that $ {\cal V} $ is dense in $ H $ and $ V $ for their respective topologies ^[10,29].

Once we project $ L^2(\Omega)^3 $ over the subspace $ H $ of divergence-free and tangential vector fields by means of the Leray projection operator $ \mathrm{P} $, we have the following abstract (functional) equation in $ H $

$\begin{equation*} \frac{d {\bf u}}{dt}+\nu A{\bf u}+B({\bf u}, {\bf u}) = \mathrm{P}\, {\bf f}, \end{equation*}$ $

where $ A: = -\mathrm{P}\, \Delta $, while $ B({\bf u}, {\bf u}): = \mathrm{P} \, (({\bf u}\cdot\nabla)\, {\bf u}) $. As usual in this analysis (see for instance ^[12]), we can start by assuming that the input force can be decomposed within a finite sum of basis functions (or that it belongs to $ V_{m} $, which will be clarified in the next section, in particular by Theorem 2.1), hence

$\begin{equation*} \mathbf{P}_m \, {\bf f} = {\bf f}. \end{equation*}$ $

We split the Navier-Stokes equations into a coupled system for $ {\bf y}\in {P}_m H $ and $ {\bf z}\in ({P}_m H)^\perp = {Q}_m H $ as follows

$\begin{equation} \begin{aligned} &\frac{d {\bf y}}{dt}-\nu\, {P}_m (\Delta{\bf u})+{P}_m B({\bf y}+{\bf z}, {\bf y}+{\bf z}) = {P}_m \, {\bf f}, \\ & \frac{d {\bf z}}{dt}-\nu\, {Q}_m(\Delta{\bf u})+{Q}_m B({\bf y}+{\bf z}, {\bf y}+{\bf z}) = \mathbf{0}, \end{aligned} \end{equation}$ $

(2.2)

where we used that both $ {P}_m $ and $ {Q}_m $ commute with the time derivative.

Once we evaluate the kinetic energy, since $ {P}_m{\bf y} = {\bf y} $ and $ {Q}_m\, {\bf z} = {\bf z} $ we get (by integrating by parts and by using the fact that functions vanish at the boundary) that

$\begin{equation*} \begin{aligned} & -\nu\int_\Omega {P}_m (\Delta{\bf u})\cdot {\bf y}\, dx = -\nu\int_\Omega (\Delta{\bf u})\cdot {\bf y}\, dx = -\nu\int_\Omega (\Delta{\bf y}+\Delta{\bf z})\cdot {\bf y}\, dx = \nu\, \|\nabla{\bf y}\|^2+\nu\int_\Omega\nabla {\bf y}:\nabla{\bf z}\, dx, \\ & -\nu\int_\Omega {Q}_m (\Delta{\bf u})\cdot {\bf z}\, dx = -\nu\int_\Omega (\Delta{\bf u})\cdot {\bf z}\, dx = -\nu\int_\Omega (\Delta{\bf y}+\Delta{\bf z})\cdot {\bf z}\, dx = \nu\, \|\nabla{\bf z}\|^2+\nu\int_\Omega\nabla {\bf y}:\nabla{\bf z}\, dx. \end{aligned} \end{equation*}$ $

In this way we can obtain the following equality and inequality

$\begin{equation} \begin{aligned} &\frac{1}{2}\frac{d }{dt}\|{\bf y}\|^2+\nu \|\nabla {\bf y}\|^2+\nu(\nabla{\bf y}, \nabla{\bf z})+(B({\bf y}+{\bf z}, {\bf y}+{\bf z}), {\bf y}) = ({\bf f}, {\bf y}), \\ & \frac{1}{2}\frac{d }{dt}\|{\bf z}\|^2+\nu \|\nabla {\bf z}\|^2+\nu(\nabla{\bf y}, \nabla{\bf z})+(B({\bf y}+{\bf z}, {\bf y}+{\bf z}), {\bf z})\leq\mathbf{0}, \end{aligned} \end{equation}$ $

(2.3)

These are the two basic balance equations that we will use to infer the behavior and transfer of the kinetic energy between $ {\bf y} $ and $ {\bf z} $. Notice that the balance relation for $ {\bf y} $, involving just a finite combination of rather smooth functions is an equality, while the second one is an inequality. In fact, the second one can be derived by a limiting argument and in the limit the lower semi-continuity of the norm will produce the inequality.

Since the tri-linear term $ (B({\bf u}, {\bf u}), {\bf w}) $ is skew-symmetric with respect to the last two variables, we obtain from (2.3)

$\begin{equation} \begin{aligned} &\frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf y})+\nu \|\nabla {\bf y}\|^2+\nu(\nabla{\bf y}, \nabla{\bf z}) = (B({\bf y}, {\bf y}), {\bf z})-(B({\bf z}, {\bf z}), {\bf y})+({\bf f}, {\bf y}), \\ & \frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf z})+\nu \|\nabla {\bf z}\|^2+\nu(\nabla{\bf y}, \nabla{\bf z}) \leq-(B({\bf y}, {\bf y}), {\bf z})+(B({\bf z}, {\bf z}), {\bf y}). \end{aligned} \end{equation}$ $

(2.4)

This is a formal setting, which is obviously true for strong solutions of the NSE, where the inequality in (2.4) is an equality. When considering weak solutions, the integral $ (B({\bf z}, {\bf z}), {\bf z}) $ might be not defined in $ L^1(0, T) $ for regularity issues. However, one can still rigorously derive (2.4) by a double frequency truncation or a regularization of the operator $ B $ by considering $ (B({\bf z}\star\rho_\varepsilon, {\bf z}), {\bf z}) $ for a given standard mollifier $ \rho_ \varepsilon $ and passing to the limit when $ \varepsilon \to 0 $. Details are standard and out of the scope of the present paper.

We observe that $ -(B({\bf y}, {\bf y}), {\bf z}) $ is the energy flux induced in the more energetic terms by the inertial forces associated to less energetic modes, while $ -(B({\bf z}, {\bf z}), {\bf y}) $ is the energy flux induced in the less energetic terms by the inertial forces associated to more energetic modes. In a schematic way we can decompose the rate of transfer of kinetic energy $ \mathrm{e}_m({\bf u}) $ into two terms as follows

$\begin{equation} \mathrm{e}_m({\bf u}): = e^{\uparrow}({\bf u})- e^{\downarrow}({\bf u})\qquad \text{ with}\qquad e^{\uparrow}({\bf u}): = -(B({\bf y}, {\bf y}), {\bf z}), \qquad e^{\downarrow}({\bf u}): = (B({\bf z}, {\bf z}), {\bf y}). \end{equation}$ $

(2.5)

We also use the following notation:

$\begin{equation} \mathcal{E}_{m}({\bf u}) : = - \nu(\nabla{\bf y}, \nabla{\bf z}). \end{equation}$ $

(2.6)

Hence, we can rewrite (2.4) as follows

$\begin{equation} \begin{aligned} & \frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf y})+\nu \|\nabla {\bf y}\|^2 = \mathcal{E}_{m}({\bf u})-\mathrm{e}_m({\bf u})+({\bf f}, {\bf y}), \\ & \frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf z})+\nu \|\nabla {\bf z}\|^2 \leq \mathcal{E}_{m}({\bf u})+\mathrm{e}_m({\bf u}). \end{aligned} \end{equation}$ $

(2.7)

Remark 2.1. We recall that apart from extremely simple geometries and provided one is willing to use in a systematic way special functions as the Bessel ones or the spherical harmonics (which are nevertheless time consuming in their evaluation), the explicit calculations in numerical tests will not be so easy to be obtained in a precise and efficient way. Hence, the solution of (2.2) and the long-time integration of its solution pose hard numerical problems.

We point out for the reader that we have a first fundamental difference with respect to the classical splitting based on the use of a spectral basis (which will be recalled in Section 2.1), where the latter integral vanishes exactly. For this reason, in the next section we will show the derivation of the corresponding system of equations, which holds when the eigenfunctions are used.

2.1. On the spectral decomposition

In this section, we compare the results of the previous section with the well-established ones that can be proved if the spectral decomposition (i.e., that made with eigenfunctions of the Stokes operator $ \{\mathcal{W}_k\} $) is used instead of utilizing a generic POD basis. We recall that, by classical results about compact operators in Hilbert spaces there exists a sequence of smooth functions $ \{\mathcal{W}_{k}\} $ (and their regularity is depending on the smoothness of the bounded domain $ \Omega $) and an increasing sequence of positive numbers $ \{\lambda_{k}\} $ such that

$\begin{equation*} A\mathcal{W}_k = \lambda_k\mathcal{W}_k\qquad \text{and}\qquad \int_\Omega\mathcal{W}_k\cdot \mathcal{W}_j\, dx = \delta_{k j}. \end{equation*}$ $

Since each function $ \mathcal{W}_k $ solves the following Stokes system $ A\mathcal{W}_k = \lambda_k \mathcal{W}_k $, it follows by an integration by parts that

$\begin{equation*} \int_\Omega\nabla\mathcal{W}_k: \nabla \mathcal{W}_j\, dx = 0\qquad\text{for }k\not = j, \end{equation*}$ $

hence also the $ V $-orthogonality of the family $ \{\mathcal{W}_k\}_{k\in \mathbb{N}} $.

We consider now the usual decomposition by eigenfunctions associated with low and high frequencies

$\begin{equation*} {\bf u} = {\bf y}+{\bf z}: = \sum\limits_{k = 1}^m c_k \mathcal{W}_k+\sum\limits_{k = m+1}^\infty\ c_k \mathcal{W}_k = \mathbf{P}_m\, {\bf u}+\mathbf{Q}_m\, {\bf u}, \end{equation*}$ $

where $ \mathbf{P}_m $ is the projection over the subspace generated by $ \{\mathcal{W}_k\}_{k = 1, \dots, m} $. Our main result is based on a standard result about the projector $ \mathbf{P}_m $, that can be found in [30,App. A.4,Thm. 4.11]:

Theorem 2.1. The projector $ \mathbf{P}_m $ can be defined as a continuous endomorphism over $ V $, $ H $ and $ V' $, and one has

$\begin{equation*} \| \mathbf{P}_m \|_{ {\cal L} (V, V)} \le 1, \qquad \| \mathbf{P}_m \|_{ {\cal L} (H, H)} \le 1, \qquad \| \mathbf{P}_m \|_{ {\cal L} (V', V')} \le 1. \end{equation*}$ $

The result is mainly based on the regularity of solutions of elliptic equations, and thanks to this fact, it is possible to decompose the equations for the velocity, which yields,

$\begin{equation*} \begin{aligned} & -\nu\int_\Omega \mathbf{P}_m (\Delta{\bf u})\cdot {\bf y}\, dx = -\nu\int_\Omega \Delta{\bf y}\cdot {\bf y}\, dx = \nu\, \|\nabla{\bf y}\|^2, \\ &-\nu\int_\Omega \mathbf{Q}_m\, (\Delta{\bf u})\cdot {\bf z}\, dx = -\nu\int_\Omega \Delta{\bf z}\cdot {\bf z}\, dx = \nu\, \|\nabla{\bf z}\|^2, \end{aligned} \end{equation*}$ $

since $ \mathbf{P}_m\Delta{\bf u} = \Delta \mathbf{P}_m\, {\bf u} = \Delta{\bf y} $ and also $ \mathbf{Q}_m (\Delta{\bf u}) = \Delta\mathbf{Q}_m\, {\bf u} = \Delta{\bf z} $.

Thus, we directly obtain the system

$\begin{equation} \begin{aligned} &\frac{1}{2}\frac{d }{dt}\|{\bf y}\|^2+\nu \|\nabla {\bf y}\|^2+(B({\bf y}+{\bf z}, {\bf y}+{\bf z}), {\bf y}) = ({\bf f}, {\bf y}), \\ & \frac{1}{2}\frac{d }{dt}\|{\bf z}\|^2+\nu \|\nabla {\bf z}\|^2+(B({\bf y}+{\bf z}, {\bf y}+{\bf z}), {\bf z})\leq\mathbf{0}, \end{aligned} \end{equation}$ $

(2.8)

which can be rewritten also as

$\begin{equation} \begin{aligned} & \frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf y})+\nu \|\nabla {\bf y}\|^2 = -\mathrm{e}_m({\bf u})+({\bf f}, {\bf y}), \\ & \frac{1}{2}\frac{d }{dt} \mathrm{E}({\bf z})+\nu \|\nabla {\bf z}\|^2 \leq \mathrm{e}_m({\bf u}). \end{aligned} \end{equation}$ $

(2.9)

We note that the equations in (2.9) do not contain the term $ \mathcal{E}_{m}({\bf u}) $, whereas the equations in (2.7) contain $ \mathcal{E}_{m}({\bf u}) $. This will have important consequences in the analysis in Section 4.

3. Preliminaries on long-time averages

Since we consider long-time averages for the NSE, we must consider solutions which are global-in-time (defined for all positive times). Due to the well-known open problems related to the NSE, this forces us to restrict ourselves to Leray-Hopf weak solutions ^[6,29]. By using a natural setting, we take the initial datum $ {\bf u}(0) $ in $ H $. The classical Leray-Hopf result of existence (but not uniqueness) of a global weak solution $ {\bf u} $ to the NSE holds when $ {\bf f}\in V' $, and the velocity $ {\bf u} $ satisfies

$\begin{equation*} {\bf u} \in L^2( \mathbb{R}_+; V) \cap L^\infty( \mathbb{R}_+; H). \end{equation*}$ $

Notice that we consider in this paper the case where $ {\bf f} $ is time-independent, for simplicity. However, the following results can be extended to the case where $ {\bf f} = {\bf f}(t) $ is time dependent, for $ {\bf f} $ belonging to a suitable class (see ^[3]).

In order to properly set what we mean by long-time-averaging, let $ \psi:\ \mathbb{R}^+\times\Omega\to \mathbb{R}^N $ be any tensor field related to a given turbulent flow ($ N $ being its order). The time-average of $ \psi $ over a time interval $ [0, t] $ is defined by

$\begin{equation} M_{t}(\psi )(x) : = \frac{1}{t}\int_0^t \psi (s, x)\, ds\qquad \text{for }t \gt 0. \end{equation}$ $

(3.1)

According to the standard turbulence modeling process, we then apply the averaging operator $ M_t $ to NSE (1.1) and also to (2.2), to study the limits when $ t\to+\infty $. We recall that the long-time averages represent one of the few observable and computable quantities associated to a highly variable turbulent flow. We will adopt the following standard notation for the long-time average of any field $ \psi $

$\begin{equation*} \overline{{\psi}}(x) : = \lim\limits_{t\to+\infty} M_t( \psi) (x), \end{equation*}$ $

whenever the limit exists. (Without too much restrictions, we can suppose that the limits we write do exist, at least after extracting sub-sequences leaving the mathematical difficulties, which can be treated with generalized Banach limits, for a more general and abstract framework.) Within this theory we can decompose the velocity as follows

$\begin{equation*} {\bf u} = \overline{{\bf u}}+{\bf u}', \end{equation*}$ $

where $ {\bf u}': = {\bf u}-\overline{{{\bf u}}} $ represents the so-called turbulent fluctuations. We recall that time-averaging has been introduced by O. Reynolds ^[37], at least for large values of $ t $, and the ideas have been widely developed by L. Prandtl ^[32] in the case of turbulent channel flows. The same ideas have been also later considered in the case of fully developed homogeneous and isotropic turbulence, such as grid-generated turbulence. In this case the velocity field is postulated as oscillating around a mean smoother steady state, see for instance G.-K. Batchelor ^[1]. For further details on the role of time averaging in turbulence, after the work of Stokes and Reynolds, we can recall a few recent papers and books ^{[2,3,5,11,18,26,28]}, where aspects of computation and modeling are studied.

We now observe that, by taking the time-averages of the NSE we have the following estimates, see [28,Prop. 2.1]

$\begin{equation} \begin{aligned} \|{\bf u}(t)\|^2&\leq \|{\bf u}(0)\|^2\, \mathrm{e}^{-\nu\, C_P t}+\frac{\|{\bf f}\|^2}{\nu^2\, C_{P}}\big(1- \mathrm{e}^{-\nu\, C_P t}\big), \qquad \forall\, t \gt 0 \\ \frac{1}{t}\int_0^t\|\nabla {\bf u}(s)\|^2\, ds&\leq\frac{\|{\bf f}\|^2}{\nu^2}+\frac{\|{\bf u}(0)\|^2}{\nu t}, \qquad \forall\, t \gt 0, \end{aligned} \end{equation}$ $

(3.2)

where $ C_P $ is the best constant in the Poincaré inequality

$\begin{equation*} C_P\|{\bf u}\|^2\leq\|\nabla {\bf u}\|^2\qquad \forall\, {\bf u}\in H^1_0(\Omega). \end{equation*}$ $

The above inequalities show that both $ \|{\bf u}(t)\|^2 $ and $ \frac{1}{t}\int_0^t\|\nabla{\bf u}(s)\|^2\, ds $ are uniformly bounded for all $ t\geq 1 $ (any other positive time will be enough), hence we have the following result:

Theorem 3.1 (cf. ^[3,28]). Let $ {\bf u}(0) \in H $, $ {\bf f}\in V' $, and let $ {\bf u} $ be a global-in-time weak solution to the NSE (1.1). Then, there exist

1. a sequence $ \{t_n\}_{n\in \mathbb{N}} $ such that $ \lim\limits_{n \to \infty }t_n = +\infty $;

2. a vector field $ \overline{{{\bf u}}}\in V $;

3. a vector field $ {\bf B}\in L^{3/2}(\Omega)^3 $;

4. a second order tensor field $ \boldsymbol {\sigma}^{(\text{R})} \in L^3(\Omega)^9 $;

such that it holds:

i) When $ n \to \infty $,

$\begin{equation*} \begin{aligned} &M_{t_n}({\bf u})\rightharpoonup \overline{{{\bf u}}}\qquad\mathit{\text{weakly in }}V, \\ & M_{t_{n}}\big (({\bf u}\cdot\nabla)\, {\bf u}\big )\rightharpoonup{\bf B} \qquad \mathit{\text{weakly in }} L^{3/2}(\Omega)^3, \\ &M_{t_n}({\bf u}{'} \otimes {\bf u}{'} )\rightharpoonup \boldsymbol {\sigma}^{(\text{R})}\qquad\mathit{\text{weakly in }}L^3(\Omega)^9; \end{aligned} \end{equation*}$ $

ii) The Reynolds averaged equations:

$\begin{equation} \left\{ \begin{aligned} (\overline{{{\bf u}}}\cdot\nabla)\, \overline{{{\bf u}}} -\nu\Delta \overline{{{\bf u}}}+\nabla \overline{{p}}+\nabla\cdot\boldsymbol {\sigma}^{(\text{R})}& = \overline{{{\bf f}}}\qquad &\mathit{\text{in }}\Omega, \\ \nabla\cdot \overline{{{\bf u}}}& = 0\qquad &\mathit{\text{in }}\Omega, \\ \overline{{{\bf u}}}& = \mathbf{0}\qquad &\mathit{\text{on }}\Gamma, \end{aligned} \right. \end{equation}$ $

(3.3)

hold true in the weak sense;

iii) The equality $ {\bf B} - (\overline{{{\bf u}}} \cdot \nabla)\, \overline{{{\bf u}}} = \nabla \cdot \boldsymbol {\sigma}^{(\text{R})} $ is valid in $ \mathcal{D}'(\Omega) $;

iv) The following energy balance (equality) holds true

$\begin{equation*} \nu\, \|\nabla\overline{{{\bf u}}}\|^2+(\nabla \cdot \boldsymbol {\sigma}^{(\text{R})}, \overline{{{\bf u}}} ) = \, \lt \overline{{{\bf f}}}, \overline{{{\bf u}}} \gt ; \end{equation*}$ $

v) The tensor $ \boldsymbol {\sigma}^{(\text{R})} $ is dissipative on the average or, more precisely, the following inequality

$\begin{equation} \epsilon: = \nu\, \overline{{\|\nabla{{\bf u}}'\|^2}}\leq\int_\Omega(\nabla\cdot\boldsymbol {\sigma}^{(\text{R})})\cdot\overline{{{\bf u}}}\, dx, \end{equation}$ $

(3.4)

holds true.

It is important to observe that the long-time limit is characterized by the solution of the system (3.3), which is similar to the Navier-Stokes equations, but which contains an extra term, coming from the effect of fluctuations, which has the mean effect of increasing the dissipation.

We observe that this is related to the long-time behavior of solutions close to statistical equilibrium. The study of the long-time behavior dates back to pioneering works of Foias and Prodi on deterministic statistical solutions, see, for instance ^[12]. Their interest is mainly devoted to finding measure in the space of initial data to be connected with the long-time limits. Here, we follow a slightly different path, as in ^[3,28], in order to characterize in a less technical way the long-time behavior, without resorting to any ergodic-type result and also with the perspective that long time averages are computable or at least can be approximated in a clear way.

4. Average transfer of energy at equilibrium

Our goal in this section is to characterize in some sense the energy transfer between the two functions $ {\bf y} $ and $ {\bf z} $ of the expansion and to determine –if possible– the sign of $ \mathrm{e}_m({\bf u}) $, at least in an average sense. We will consider both the POD case and the spectral one.

4.1. The POD case

The point concerning the exchange of energy between low and high modes is in the same spirit as the results recalled in Foias, Manley, Rosa, and Temam [12,Chap. 5] and follows from results obtained in a more heuristic way by Kolmogorov ^[19].

We first observe that the $ L^2 $-orthogonality of the POD decomposition implies that

$\begin{equation*} \|{\bf u}\|^2 = \|{\bf y}+{\bf z}\|^2 = \|{\bf y}\|^2+\|{\bf z}\|^2. \end{equation*}$ $

Hence, from the uniform $ L^2 $-bound on $ {\bf u} $ it follows that both $ {\bf y} $ and $ {\bf z} $ are uniformly bounded in time. From this observation we can deduce the following result, reminding that $ \mathrm{e}_m $ and $ \mathcal{E}_{m} $ are defined by equations (2.5) and (2.5), and $ M_t $ is defined by equation (3.1).

Theorem 4.1. Let $ {\bf z} $ be the projection onto the less energetic POD modes of a weak solution $ {\bf u} $ to the Navier-Stokes equations. Then, there exists a sequence $ \{t_n\} $ such that $ t_n\to+\infty $ and a field $ \overline{{{\bf z}}}\in H $ such that

$\begin{equation} {\bf Z}_{t_n} = M_{t_n}({\bf z})\rightharpoonup \overline{{{\bf z}}}\qquad\mathit{\text{weakly in }}H, \end{equation}$ $

(4.1)

and

$\begin{equation} \liminf\limits_{n\to+\infty}M_{t_{n}}( \mathrm{e}_m({\bf u})+\mathcal{E}_{m}({\bf u})) \geq0. \end{equation}$ $

(4.2)

Proof. Let us observe first that by the energy inequality (3.2), we easily deduce that $ (M_t ({\bf z}))_{t > 0} $ is bounded in $ H $, hence the first assertion of the statement and (4.1). We next prove (4.2). To do so, we average with respect to time with the operator $ M_t $ the balance equation (2.7) for $ \mathrm{E}({\bf z}) $, which yields

$\begin{equation} \frac{1}{2t}\|{\bf z}(t)\|^2-\frac{1}{2t}\|{\bf z}(0)\|^2+\nu M_t (\|\nabla{\bf z}\|^2)\leq M_t( \mathrm{e}_m({\bf u})+\mathcal{E}_{m}({\bf u})). \end{equation}$ $

(4.3)

By using the energy inequality (3.2) once again, we see that the first two terms vanish as $ t\to+\infty $ and also that $ M_t (\|\nabla {\bf z}\|^2) $ is bounded. Therefore, (4.3) yields

$\begin{equation*} 0 \le \nu \liminf\limits_{n \to +\infty} M_{t_n} (\|\nabla {\bf z}\|^2) \le \liminf\limits_{n\to+\infty}M_{t_n} (\mathrm{e}_m({\bf u})+\mathcal{E}_{m}({\bf u})), \end{equation*}$ $

hence (4.2). We observe that in this case we do not have any direct estimation on the behavior of the $ H^1 $-norm.

In the case we can assume that the limit exists, we also have the following result.

Corollary 4.1. Let us assume the limit of $ M_{t_n}(\mathrm{e}_m({\bf u})) $ for $ n\to+\infty $ exists, and that

$\begin{equation*} \liminf\limits_{T\to+\infty}\frac{\nu}{T}\int_0^T(\nabla{\bf z}(s), \nabla{\bf y}(s))\, ds = \liminf\limits_{t \to \infty} \mathcal{E}_{m}({\bf u}) \geq0. \end{equation*}$ $

Then, it follows

$\begin{equation*} \overline{\mathrm{e}_m({\bf u})} = \lim\limits_{n\to+\infty}\frac{1}{t_n}\int_0^{t_n} \mathrm{e}_m({\bf u}(s))\, ds\geq0. \end{equation*}$ $

This result can be interpreted as that, beyond the range of injection of energy, the average net transfer of energy occurs only into the small scales. This occurs if the term of interaction between gradients of large and small scales is negligible, in the limit of long time. This latter assumption is not proved rigorously, but we will see it is satisfied in the numerical tests, with a good degree of approximation (see Section 5). However, when one uses the eigenfunctions of the Stokes operator as POD basis, this is automatically satisfied since this basis is also orthogonal for the $ H^1 $-scalar product, so that in this case $ \mathcal{E}_{m}({\bf u}) = 0 $.

4.2. The spectral case

The results of the previous section can be made much more precise in the case of decomposition made by a spectral basis of eigenfunctions of the Stokes operator. We present the results, which are in some sense new and not fully completely included in ^[12], in the sense of time-averaging. This procedure is applied to $ {\bf u} $, which is a weak solution of the Navier-Stokes equations, satisfying the uniform estimates (3.2). In this way, the orthogonality (in both $ H $ and $ V $) of the basis implies that

$\begin{equation*} \|{\bf u}\|^2 = \|{\bf y}\|^2+\|{\bf z}\|^2\qquad\text{and}\qquad \|\nabla {\bf u}\|^2 = \|\nabla {\bf y}\|^2+\|\nabla {\bf z}\|^2. \end{equation*}$ $

The results in this case are more precise than those from Theorem 4.1, since we have at disposal a larger set of a priori estimates and also the set of equations (2.8) has a better structure than (2.3).

We now prove the following results in the case of a decomposition of the velocity into small and large frequencies. The first one aims at taking the time average and then let $ t $ go to infinity in the equations (2.3) satisfied by $ {\bf y} $ and $ {\bf z} $. The second one aims at comparing the amount of turbulent dissipation of small and large frequencies with respect to the total work of the corresponding Reynolds stresses $ \boldsymbol {\sigma}^{(\text{R})}_{\bf y} $ and $ \boldsymbol {\sigma}^{(\text{R})}_{\bf z} $.

Theorem 4.2. Let $ {\bf u}(0) \in H $, $ {\bf f}\in \mathbf{P}_m\, H $, and let $ {\bf u} $ be a global-in-time weak solution to the NSE (1.1). Then, there exist

1. a sequence $ \{t_n\}_{n\in \mathbb{N}} $ such that $ \lim\limits_{n \to \infty }t_n = +\infty $;

2. vector fields $ \overline{{{\bf y}}}, \overline{{{\bf z}}}\in V $;

3. vector fields $ {\bf B}_1, {\bf B}_2\in V' $;

such that it holds:

i) When $ n \to \infty $,

$\begin{equation*} \begin{aligned} \begin{aligned} &M_{t_n}({\bf y})\rightharpoonup \overline{{{\bf y}}}\qquad\mathit{\text{weakly in }}V, \\ & M_{t_{n}}\big (({\bf y}\cdot\nabla)\, {\bf y}\big )\rightharpoonup{\bf B}_1 \qquad \mathit{\text{weakly in }}V', \end{aligned}& \quad \begin{aligned} &M_{t_n}({\bf z})\rightharpoonup \overline{{{\bf z}}}\qquad\mathit{\text{weakly in }}V, \\ &M_{t_{n}}\big (({\bf z}\cdot\nabla)\, {\bf z}\big )\rightharpoonup{\bf B}_2 \qquad \mathit{\text{weakly in }} V', \end{aligned} \end{aligned} \end{equation*}$ $

ii) The Reynolds averaged equations:

$\begin{equation} \left\{ \begin{aligned} -\nu\Delta \overline{{{\bf y}}}+\nabla \overline{{p_{\bf y}}}+{\bf B}_1& = \mathbf{P}_m{\bf f}\qquad &\mathit{\text{in }}\Omega, \\ \nabla\cdot \overline{{{\bf y}}}& = 0\qquad &\mathit{\text{in }}\Omega, \\ \overline{{{\bf y}}}& = \mathbf{0}\qquad &\mathit{\text{on }}\Gamma, \end{aligned} \right. \end{equation}$ $

(4.4)

and

$\begin{equation} \left\{ \begin{aligned} -\nu\Delta \overline{{{\bf z}}}+\nabla \overline{{p_{\bf z}}}+{\bf B}_2& = \mathbf{0}\qquad &\mathit{\text{in }}\Omega, \\ \nabla\cdot \overline{{{\bf z}}}& = 0\qquad &\mathit{\text{in }}\Omega, \\ \overline{{{\bf z}}}& = \mathbf{0}\qquad &\mathit{\text{on }}\Gamma, \end{aligned} \right. \end{equation}$ $

(4.5)

holds true in $ V' $.

Arguing as in ^[3,28], using the relations $ (\overline{{{\bf z}}}\cdot\nabla)\, \overline{{{\bf z}}} = \nabla\cdot (\overline{{{\bf z}}} \otimes \overline{{{\bf z}}}) $ and $ (\overline{{{\bf y}}}\cdot\nabla)\, \overline{{{\bf y}}} = \nabla\cdot (\overline{{{\bf y}}} \otimes \overline{{{\bf y}}}) $, we get the existence of "small frequencies" and "large frequencies" Reynolds stresses $ \boldsymbol {\sigma}^{(\text{R})}_{\bf y} $ and $ \boldsymbol {\sigma}^{(\text{R})}_{\bf z} $ in $ V' $, such that

$\begin{equation*} B_1 = \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y} + (\overline{{{\bf y}}} \cdot \nabla)\, \overline{{{\bf y}}}\qquad\text{and}\qquad B_2 = \nabla\cdot\boldsymbol {\sigma}^{(\text{R})}_{\bf z}+(\overline{{{\bf z}}} \cdot \nabla)\, \overline{{{\bf z}}}, \end{equation*}$ $

or equivalently, if we write the Reynolds decomposition as

$\begin{equation*} {\bf y} = \overline{{{\bf y}}} + {\bf y}' \quad\text{and} \quad {\bf z} = \overline{{{\bf z}}} + {\bf z}', \end{equation*}$ $

then

$\begin{equation*} \boldsymbol {\sigma}^{(\text{R})}_{\bf y} = \overline { {\bf y}' \otimes {\bf y}'}\qquad\text{and} \qquad \boldsymbol {\sigma}^{(\text{R})}_{\bf z} = \overline { {\bf z}' \otimes {\bf z}'}, \end{equation*}$ $

where the bar operator denotes the limit of the $ M_{t_n} $'s in $ V' $ as $ n \to \infty $ (eventually after having extracted another sub-sequence).

According to the budget (3.4), we aim to compare the turbulent dissipation of small and large scales, denoted as $ \varepsilon^\downarrow $ and $ \varepsilon^\uparrow $ respectively, to the total work of the Reynolds stresses, namely $ (\nabla \cdot\boldsymbol {\sigma}^{(\text{R})}_{\bf y}, \overline{{{\bf y}}}) $ and $ (\nabla \cdot\boldsymbol {\sigma}^{(\text{R})}_{\bf z}, \overline{{{\bf z}}}) $, where

$\begin{equation*} \epsilon^\downarrow: = \nu\, \overline{{\|\nabla{{\bf y}}'\|^2}}\quad\text{and} \quad \epsilon^\uparrow: = \nu\, \overline{{\|\nabla{{\bf z}}'\|^2}}. \end{equation*}$ $

When we compare to (3.4), we observe that triad nonlinear effect between small and large frequencies will be felt, that means the nonlinear interactions due to the convection will be provided by the term

$\begin{equation} \Phi_{{\bf z}}({\bf y}): = (\mathbf{Q}_m [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})], {\bf z}) = -(\mathbf{P}_m [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})], {\bf y}) = -\Phi_{{\bf y}}({\bf z}). \end{equation}$ $

(4.6)

Notice that due to the regularity of $ {\bf y} $, it is easily checked that the following energy balance holds true (this property will be shown with more details in the proof of Theorem 4.2 below)

$\begin{equation*} \nu\, \|\nabla\overline{{{\bf y}}}\|^2+ ( \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y}, \overline{{{\bf y}}}) = \, \lt {{\bf f}}, \overline{{{\bf y}}} \gt . \end{equation*}$ $

Finally, to prove Theorem 4.3 we will use the following orthogonality relation (see e.g., [3,Lemma 4.4]), formally written as follows

$\begin{equation} \overline {\| \nabla \psi \|^2} = \| \nabla \overline{{\psi }} \|^2 + \overline { \| \nabla \psi' \|^2 }. \end{equation}$ $

(4.7)

which is valid for any field $ \psi:\, \mathbb{R}^+\to V $, such that $ \overline{{\psi}} $ is well-defined and the fluctuations are defined as $ \psi' = \psi-\overline{{\psi}} $.

Theorem 4.3. The families $ (M_t(\Phi_{{\bf z}}({\bf y})))_{t > 0} $ and $ (M_t(\Phi_{{\bf y}}({\bf z})))_{t > 0} $ converge (along certain subsequences) as $ t \to \infty $. Let $ \overline{{\Phi_{{\bf z}}({\bf y})}} $ and $ \overline{{\Phi_{{\bf y}}({\bf z})}} $ denote the corresponding limits. One has

$\begin{equation} \overline{{\Phi_{{\bf z}}({\bf y})}} = - \overline{{\Phi_{{\bf y}}({\bf z})}}\le 0, \end{equation}$ $

(4.8)

and the following dissipation balances hold true

$\begin{eqnarray} && \epsilon^\downarrow + \overline{{\Phi_{{\bf y}}({\bf z})}} = ( \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y}, \overline{{{\bf y}}}), \end{eqnarray}$ $

(4.9)

$\begin{eqnarray} && \epsilon^\uparrow + \overline{{\Phi_{{\bf z}}({\bf y})}} \le ( \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf z}, \overline{{{\bf z}}}). \end{eqnarray}$ $

(4.10)

Remark 4.1. Notice that by equations (4.8) and (4.9) we see that $ \boldsymbol {\sigma}^{(\text{R})}_{\bf y} $ is dissipative in mean, and follows the same law (3.4) as the complete Reynolds stress, namely

$\begin{equation*} \epsilon^\downarrow \le ( \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y}, \overline{{{\bf y}}}). \end{equation*}$ $

However, nothing similar can be concluded from (4.9) about $ \boldsymbol {\sigma}^{(\text{R})}_{\bf z} $, that might be at this stage non dissipative in mean, which permits an inverse energy cascade to occur.

The results of Theorems 4.2 and 4.3 are original, even if similar results have been already obtained in ^[11] and reported also in ^[12]. In that case, the results are based on the notion of deterministic statistical solutions and on a sort of ergodic hypothesis. Even if statements could look very similar to ours, the main difference is that we do not average over the set $ H $ of initial data, and we do not introduce probability measures on $ H $, as suggested by the work by Prodi ^[33,34]. Our approach is based on a more elementary functional setting and also amenable to include treatment of sets of external forces, as those in several numerical or practical experiments. The main point is an extension of the results in ^[3].

Proof of Theorem 4.2. We know, from the results in ^[3,28] that $ \mathbf{U}_t = M_t ({\bf u}) $ is such that

$\begin{equation*} \begin{aligned} & \mathbf{U}_t\rightharpoonup \overline{{\bf u}}\qquad \text{weakly in }V, \\ &M_t (({\bf u}\cdot\nabla)\, {\bf u})\rightharpoonup \mathbf{B}\qquad \text{ in }L^{3/2}(\Omega)\subset V' , \end{aligned} \end{equation*}$ $

hence, if we define $ \mathbf{F}: = \mathbf{B}-(\overline{{\bf u}}\cdot\nabla)\, \overline{{\bf u}} $, we get

$\begin{equation*} \nu\, (\nabla \overline{{\bf u}}, \nabla {\pmb\phi})+\big((\overline{{\bf u}}\cdot\nabla)\, \overline{{\bf u}}, {\pmb\phi}\big) + \lt \mathbf{F}, {\pmb\phi} \gt \, = \, \lt {\bf f}, {\pmb\phi} \gt , \end{equation*}$ $

and using $ \overline{{\bf u}}\in V $ as test function we obtain

$\begin{equation*} \nu\, \|\nabla \overline{{\bf u}}\|^2+ \lt \mathbf{F}, \overline{{\bf u}} \gt \, = \, \lt {\bf f}, \overline{{\bf u}} \gt . \end{equation*}$ $

We assume now that $ \mathbf{P}_m {\bf f} = {\bf f} $, and we consider the equations satisfied by $ {\bf y} = \mathbf{P}_m{\bf u} $ and $ {\bf z} = \mathbf{Q_m}{\bf u} $. In particular, the equation for $ {\bf y} $ reads, as an abstract equation in $ V_m = \mathbf{P}_m V $, as follows:

$\begin{equation*} \frac{d{\bf y}}{dt} +\nu A{\bf y}+\mathbf{P}_m[({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})] = \mathbf{P}_m{\bf f}. \end{equation*}$ $

The uniform estimates on $ {\bf u} $ from Theorem 3.1 combined with Theorem 2.1, about the properties of the projection operator $ \mathbf{P}_m $ as a continuous operator over $ V' $, yield

$\begin{equation*} \begin{aligned} & M_{t}({\bf y}) = \mathbf{Y}_t\rightharpoonup \overline{{\bf y}} \qquad\text{weakly in }V, \\ &\mathbf{P}_m M_t (({\bf u}\cdot\nabla)\, {\bf u}) = \mathbf{P}_m M_t [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})]\rightharpoonup \mathbf{P}_m \mathbf{B} = \mathbf{B}_1\qquad \text{weakly in }V', \end{aligned} \end{equation*}$ $

in such a way that $ \overline{{{\bf y}}} $ satisfies, for the spectral basis $ \mathcal{W}_k $ introduced in Section 2.1,

$\begin{equation*} \nu\, (\nabla\overline{{\bf y}}, \nabla\mathcal{W}_k)+ \lt (\overline{{\bf y}}\cdot\nabla)\, \overline{{\bf y}}, \mathcal{W}_k \gt + \lt \mathbf{F}_{{\bf y}}, \mathcal{W}_k \gt \, = \, \lt {\bf f}, \mathcal{W}_k \gt \qquad \text{for all }1\leq k\leq m, \end{equation*}$ $

where $ \mathbf{F}_{{\bf y}}: = \mathbf{B}_1-(\overline{{\bf y}}\cdot\nabla)\, \overline{{\bf y}} $, which leads to (4.4) by De Rham Theorem, if written in a strong formulation. Arguing as in ^[3] (which was already mentioned above), it is easily checked that there exists $ \boldsymbol {\sigma}^{(\text{R})}_{\bf y} $ such that $ \mathbf{F}_{{\bf y}} = \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y} $. Hence, being $ \overline{{\bf y}}\in V_m\subset V $ a legitimate test function, we get

$\begin{equation} \nu\, \|\nabla\overline{{\bf y}}\|^2+ \lt \mathbf{F}_{{\bf y}}, \overline{{\bf y}} \gt \, = \, \nu\, \|\nabla\overline{{{\bf y}}}\|^2+ ( \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf y}, \overline{{{\bf y}}}) = \, \lt {\bf f}, \overline{{\bf y}} \gt . \end{equation}$ $

(4.11)

The other term $ {\bf z} $ of the decomposition satisfies

$\begin{equation*} \frac{d}{dt} {\bf z}+\nu A{\bf z}+\mathbf{Q}_m[({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})] = 0. \end{equation*}$ $

The uniform estimates on $ {\bf u} $ and the boundedness of the linear operator $ \mathbf{P}_m $ imply the following convergence (up to a sub-sequence), as already shown in Theorem 4.1

$\begin{equation*} \begin{aligned} &M_{t}({\bf z}) = \mathbf{Z}_t\rightharpoonup \overline{{\bf z}} \qquad \text{weakly in }V, \\ &\mathbf{Q}_m M_t (({\bf u}\cdot\nabla)\, {\bf u}) = \mathbf{Q}_m M_t [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})]\rightharpoonup \mathbf{Q}_m\mathbf{B} = \mathbf{B}_2 \qquad \text{weakly in }V'. \end{aligned} \end{equation*}$ $

By using that $ \mathbf{B} = \mathbf{P}_m \mathbf{B}+\mathbf{Q}_m\, \mathbf{B} $, we get

$\begin{equation} \nu\, (\nabla\overline{{\bf z}}, \nabla\mathcal{W}_j)+ \lt (\overline{{\bf z}}\cdot\nabla)\, \overline{{\bf z}}, \mathcal{W}_j \gt + \lt \mathbf{F}_{{\bf z}}, \mathcal{W}_j \gt \, = 0\qquad \text{for all } j\geq m+1, \end{equation}$ $

(4.12)

for $ \mathbf{F}_{{\bf z}} = \mathbf{B}_2-(\overline{{\bf z}}\cdot\nabla)\, \overline{{\bf z}} = \nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf z} $. Hence, (4.5) follows again by De Rham Theorem. Notice that, $ \overline{{{\bf z}}}\in V_m^\perp\subset V $ and in particular $ \overline{{{\bf z}}}\in \text{span}\{\mathcal{W}_j\}_{j\geq m+1} $, so by linearity it can be used as test function in (4.12). Next, since $ \nabla\cdot\overline{{{\bf z}}} = 0 $, it follows that $ < (\overline{{\bf z}}\cdot\nabla)\, \overline{{\bf z}}, \overline{{{\bf z}}} > = 0 $ and we have the following energy equality:

$\begin{equation} \nu\, \|\nabla \overline{{\bf z}}\|^2+(\nabla \cdot \boldsymbol {\sigma}^{(\text{R})}_{\bf z}, \overline{{\bf z}}) \, = 0, \end{equation}$ $

(4.13)

which concludes the proof.

Proof of Theorem 4.3. We now write the energy inequality for $ {\bf z} $, obtaining

$\begin{equation*} \frac{1}{2} \frac{d}{dt}\|{\bf z}\|^2+\nu\, \|\nabla {\bf z}\|^2+(\mathbf{Q}_m [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})], {\bf z})\leq0, \end{equation*}$ $

and hence, by using the orthogonality of the basis, we have that $ \mathbf{Q}_m{\bf z} = {\bf z} $ and

$\begin{equation*} \frac{1}{2}\frac{d}{dt}\|{\bf z}\|^2+\nu\, \|\nabla {\bf z}\|^2+\big( [({\bf y}+{\bf z})\cdot\nabla({\bf y}+{\bf z})], {\bf z}\big) = \frac{1}{2} \frac{d}{dt}\|{\bf z}\|^2+\|\nabla {\bf z}\|^2+\Phi_{{\bf z}}({\bf y})\leq0, \end{equation*}$ $

recalling the definition of $ \Phi_{{\bf z}}({\bf y}) $ in (4.6).

Averaging the above equation over a fixed time interval $ (0, t) $, we get

$\begin{equation*} \frac{1}{2t}\|{\bf z}(t)\|^2- \frac{1}{2t}\|{\bf z}(0)\|^2+ \nu M_t( \|\nabla{\bf z}\|^2) +M_t (\Phi_{{\bf z}}({\bf y})) \leq0. \end{equation*}$ $

The $ L^2 $-uniform bounds on $ {\bf z} $ imply that $ \frac{1}{2t}\|{\bf z}(t)\|^2\to0 $, hence, possibly after having extracted another sub-sequence to ensure the convergence of the term $ M_t(\|\nabla{\bf z}\|^2) $ (that is known to be bounded by the energy inequality (3.2)) we get

$\begin{equation} \limsup\limits_{n \to \infty} M_{t_n} (\Phi_{{\bf z}}({\bf y})) \le - \nu\, \overline{{\|\nabla{\bf z}\|^2}}\le 0. \end{equation}$ $

(4.14)

We now combine the orthogonality relation (4.7) with the energy balance (4.13), so that (4.14) yields

$\begin{equation*} \epsilon^\uparrow+\limsup\limits_{n \to \infty} M_{t_n}(\Phi_{{\bf z}}({\bf y})) \le ( \nabla \cdot\boldsymbol {\sigma}^{(\text{R})}_{\bf z}, \overline{{{\bf z}}}), \end{equation*}$ $

which is almost inequality (4.10), up to the convergence of $ (M_t(\Phi_{{\bf z}}({\bf y})))_{t > 0} $ that remains to be proved. To prove this, we deal with the budget for $ {\bf y} $, recall the definition of $ \Phi_{{\bf z}}({\bf y}) $ and $ \Phi_{{\bf y}}({\bf z}) $ from (4.6).

Then, averaging the energy equality (in this case we have equality since $ {\bf y} $ solves a finite dimensional system of ordinary differential equations) which is satisfied for $ {\bf y} $, we get

$\begin{equation} \frac{1}{2{t}}\|{\bf y}({t})\|^2- \frac{1}{2{t}}\|{\bf y}(0)\|^2+ \nu M_t(\|\nabla{\bf y}\|^2)+M_t(\Phi_{{\bf y}}({\bf z})) = \, \lt {\bf f}, M_t({\bf y}) \gt . \end{equation}$ $

(4.15)

By the same argument, eventually after having extracted a further sub-sequence, $ (M_t(\|\nabla{\bf y}\|^2))_{t > 0} $ is convergent as $ t_n \to \infty $, as well as $ (<{\bf f}, M_t({\bf y}) > )_{t > 0} $. Therefore, $ \{M_{t_n}(\Phi_{{\bf y}}({\bf z}))\} $ is also convergent by (4.15). Let $ \overline{{\Phi_{{\bf y}}({\bf z})}} $ denote its limit. In particular, by (4.6), $ \{M_{t_n} (\Phi_{{\bf z}}({\bf y})))\} $ is also convergent, with limit $ \overline{{\Phi_{{\bf z}}({\bf y})}} = - \overline{{\Phi_{{\bf y}}({\bf z})}} $. We are done with (4.10).

It remains to check (4.9). Taking the limit as $ n \to \infty $ in (4.15) gives the equality

$\begin{equation*} \nu\, \overline{ \|\nabla{\bf y}\|^2 } + \overline{{\Phi_{{\bf z}}({\bf y})}} = \, \lt {\bf f}, \overline{{\bf y}} \gt , \end{equation*}$ $

which, combined with the energy balance (4.11) and the orthogonality relation (4.7), yields (4.9), ending the proof.

5. Numerical results

In Theorem 4.1, we showed that

$\begin{equation} \liminf\limits_{T\to+\infty}\frac{1}{T}\int_0^{T} \mathrm{e}_m({\bf u}(s)) + \mathcal{E}_{m}({\bf u}(s)) \, ds \geq0 . \end{equation}$ $

(5.1)

In this section, we investigate numerically whether the inequality (5.1) holds. To this end, we consider the one-dimensional Burgers equation with homogeneous Dirichlet boundary conditions as a simplified, yet relevant test case ^{[16,17,20,21,22,23,24,31,40,43]}:

$\begin{equation} \begin{cases} u_t -\nu u_{xx} + u u_x = f \; \; \; \; \; \; \; \; \; \; \; \; (x, t) \in \Omega\times [0, 1], \\ \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; u = 0 \; \; \; \; \; \; \; \; \; \; (x, t) \in \partial\Omega\times [0, 1]. \end{cases} \end{equation}$ $

(5.2)

To calculate the long-time average of $ \mathrm{e}_m(u) $ in (5.1), we use the composite trapezoidal rule:

$\begin{eqnarray} \frac{1}{T} \int_{0}^{T} \mathrm{e}_m(u(s)) \, ds \approx \frac{1}{2n} \sum\limits_{i = 1}^{n} \bigl( \mathrm{e}_m(u(t_i))+\mathrm{e}_m(u(t_{i+1})) \bigr) , \end{eqnarray}$ $

(5.3)

where $ t_i = (i-1)*\frac{T}{n}, \; i = 1, \dots, n+1. $ We also use the composite trapezoidal rule to calculate the long-time average of $ \mathcal{E}_{m}(u) $.

5.1. Numerical results with step function initial condition

Our numerical results are obtained by using the one-dimensional Burgers equation (5.2) with a step function initial condition ^[16,43]:

$\begin{equation} u_0(x) = \begin{cases} \; 1, & x \in (0, 1/2], \\ \; 0, & x \in (1/2, 1]. \end{cases} \end{equation}$ $

(5.4)

We use the following parameters in the finite element discretization of the Burgers equation (5.2): $ \Omega = [0, 1] $, $ \nu = 10^{-2} $, $ f = 0 $, mesh size $ h = 1/128 $, piecewise linear finite element spatial discretization, and backward Euler time discretization.

5.1.1. Case 1:

For this test case, we consider the time interval $ [0, T] = [0, 1] $ and the time step $ \Delta t = 10^{-2} $. We utilize all the snapshots to build the POD basis, whose dimension is $ d = 37 $. In the composite trapezoidal rule, we use $ n = 100 $. In Figure 1, we plot the DNS results (which are used to generate the snapshots). In Table 1, we list the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for different $ m $ values. We note that the time-average of $ \mathrm{e}_m(u) $ is positive for all $ m $ values. The time-average of $ \mathcal{E}_{m}(u) $ is positive for the low $ m $ values and negative for the largest $ m $ values. Furthermore, the magnitude of the time-average of $ \mathcal{E}_m(u) $ is generally lower than the magnitude of the time-average of $ \mathrm{e}_m(u) $. Thus, we conclude that the time average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ in (5.1) is positive for all $ m $ values.

Figure 1. DNS solution obtained by using a piecewise linear finite element spatial discretization and the backward Euler time discretization.

DownLoad: Full-Size Img PowerPoint

Table 1. Case 1: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ d = 37 $ and different $ m $ values.

$ m $	$ \int_{0}^{1} \mathrm{e}_m(u(s)) \, ds $	$ \int_{0}^{1} \mathcal{E}_m(u(s)) \, ds $
3	2.6170e-02	1.0451e-03
6	7.3208e-03	2.2524e-03
9	1.4934e-03	1.5977e-03
15	3.6181e-05	3.6287e-05
20	1.2776e-06	7.0356e-07
25	4.6207e-08	9.5638e-09
30	2.0257e-09	-1.1127e-10
35	8.2806e-11	-2.2595e-11

| Show Table

DownLoad: CSV

In Case 1, we showed that the time average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ in (5.1) is positive. In the remainder of this section, we investigate whether this time average remains positive if we make the following changes in our computational setting: (ⅰ) we increase/decrease the time-interval; and (ⅱ) we use more quadrature points (i.e., subintervals) in the composite trapezoidal rule (5.3).

5.1.2. Case 2:

In this case, we use a longer time interval, i.e., $ [0, T] = [0, 10] $ (instead of $ [0, T] = [0, 1] $, as we used in Case 1). We also use different time step ($ \Delta t $) values to generate the snapshots and a different number of quadrature points to evaluate the time average $ \frac{1}{T}\int_0^{T}\mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $.

In Tables 2–5, we list the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for different time steps ($ \Delta t $) values, different number of equally spaced quadrature points ($ n $), and different $ m $ values. We note that the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ are positive for all $ \Delta t $, $ n $, and $ m $ values. Moreover, the magnitude of the time-average of $ \mathcal{E}_m(u) $ is generally lower than the magnitude of the time-average of $ \mathrm{e}_m(u) $. Thus, we conclude that the time average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ in (5.1) is positive for all $ \Delta t $, $ n $, and $ m $ values. Furthermore, we note that decreasing the time step while keeping the same number of snapshots (i.e., $ n = 10000 $) does not change the time average $ \frac{1}{T}\int_0^{T}\mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ significantly (see Tables 3–5).

Table 2. Case 2: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ d = 38 $, $ \Delta t = 10^{-2} $, $ 1000 $ equally spaced quadrature points, and different $ m $ values.

$ m $	$ \frac{1}{10} \int_{0}^{10} \mathrm{e}_m(u(s))\, ds $	$ \frac{1}{10} \int_{0}^{10} \mathcal{E}_m(u(s)) \, ds $
3	3.3652e-03	8.6523e-05
6	1.0001e-03	2.1570e-04
9	2.1132e-04	1.8614e-04
15	5.4224e-06	5.5724e-06
20	4.3119e-07	2.5667e-07
25	5.6884e-08	1.0451e-08
30	1.1327e-09	3.2736e-10
35	2.0469e-11	2.6396e-12

| Show Table

DownLoad: CSV

Table 3. Case 2: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ d = 41 $, $ \Delta t = 10^{-3} $, $ 10000 $ equally spaced quadrature points, and different $ m $ values.

$ m $	$ \frac{1}{10} \int_{0}^{10} \mathrm{e}_m(u(s))ds $	$ \frac{1}{10} \int_{0}^{10} \mathcal{E}_m(u(s)) \, ds $
3	3.5296e-03	3.0515e-06
6	1.1402e-03	8.8746e-06
9	2.6359e-04	1.5090e-05
15	1.1126e-05	1.1156e-05
20	9.5913e-07	1.5183e-06
25	1.8140e-07	1.8704e-07
30	4.8765e-08	1.3491e-08
35	1.0291e-09	9.7629e-10
40	2.4680e-11	2.0416e-11

| Show Table

DownLoad: CSV

Table 4. Case 2: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ d = 43 $, $ \Delta t = 10^{-4} $, $ 10000 $ equally spaced quadrature points, and different $ m $ values.

$ m $	$ \frac{1}{10} \int_{0}^{10} \mathrm{e}_m(u(s))ds $	$ \frac{1}{10} \int_{0}^{10} \mathcal{E}_m(u(s)) \, ds $
3	3.5637e-03	2.7841e-06
6	1.1664e-03	8.1109e-06
9	2.7346e-04	1.3956e-05
15	1.2084e-05	1.1937e-05
20	1.1785e-06	1.6370e-06
25	2.2727e-07	1.2913e-07
30	6.2606e-08	9.5191e-09
35	2.1106e-09	6.1538e-10
40	1.0330e-10	2.0241e-11

| Show Table

DownLoad: CSV

Table 5. Case 2: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ d = 43 $, $ \Delta t = 2*10^{-5} $, $ 10000 $ equally spaced quadrature points, and different $ m $ values.

$ m $	$ \frac{1}{10} \int_{0}^{10} \mathrm{e}_m(u(s))ds $	$ \frac{1}{10} \int_{0}^{10} \mathcal{E}_m(u(s)) \, ds $
3	3.5668e-03	2.7594e-06
6	1.1688e-03	8.0390e-06
9	2.7436e-04	1.3843e-05
15	1.2172e-05	1.2002e-05
20	1.2030e-06	1.6405e-06
25	2.3339e-07	1.2463e-07
30	6.4120e-08	8.8971e-09
35	2.2654e-09	5.6187e-10
40	1.1211e-10	1.8216e-11

| Show Table

DownLoad: CSV

5.1.3. Case 3:

In this case, we use an even longer time interval, i.e., $ [0, T] = [0,100] $, and compare the time-averages for this time interval to those for the time intervals $ [0, T] = [0, 1] $ (Case 1) and $ [0, T] = [0, 10] $ (Case 2). For each time interval, we use the same time step values ($ \Delta t = 10^{-2} $) to generate the snapshots and all the subintervals in the composite trapezoidal rule utilized in the evaluation of the time average $ \frac{1}{T}\int_0^{T}\mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $. In Table 6, we list the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for all three time intervals and different $ m $ values. We note that the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ are positive for all time intervals and $ m $ values. Furthermore, the magnitude of the time-average of $ \mathcal{E}_m(u) $ is generally lower than the magnitude of the time-average of $ \mathrm{e}_m(u) $. Thus, we conclude that the time-average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ in (5.1) is positive for all time intervals and $ m $ values. Furthermore, we note that the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for the time intervals $ [0, T] = [0,100] $ and $ [0, T] = [0, 10] $ are close, whereas those for the time interval $ [0, T] = [0, 1] $ are slightly different. Thus, we conclude that the time interval $ [0, T] = [0, 10] $ is adequate for the approximation of the long time-average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $.

Table 6. Case 3: Time-averages of $ \mathrm{e}_m(u) $ for $ d = 36 $, $ \Delta t = 10^{-2} $, different $ m $ values, and all subintervals used in the composite trapezoidal rule.

$ m $	$ \frac{1}{100} \int_{0}^{100} \mathrm{e}_m(u(s)) \, ds $	$ \frac{1}{10} \int_{0}^{10} \mathrm{e}_m(u(s)) \, ds $	$ \int_{0}^{1} \mathrm{e}_m(u(s)) \, ds $
3	3.3687e-04	3.3683e-04	1.7634e-04
6	1.0001e-04	1.0001e-04	7.3142e-05
9	2.1146e-05	2.1147e-05	1.6701e-05
15	5.5621e-07	5.5742e-07	4.7129e-07
20	7.9605e-08	7.9605e-08	6.3799e-08
25	8.3483e-09	8.3489e-09	5.7426e-09
30	1.0839e-10	1.0840e-10	9.2430e-11
35	2.7710e-12	2.7718e-12	2.9575e-12

| Show Table

DownLoad: CSV

Table 7. Case 3: Time-averages of $ \mathcal{E}_{m}(u) $ for $ d = 36 $, $ \Delta t = 10^{-2} $, different $ m $ values, and all subintervals used in the composite trapezoidal rule.

$ m $	$ \frac{1}{100} \int_{0}^{100} \mathcal{E}_m(u(s)) \, ds $	$ \frac{1}{10} \int_{0}^{10} \mathcal{E}_m(u(s)) \, ds $	$ \int_{0}^{1} \mathcal{E}_m(u(s)) \, ds $
3	8.6270e-06	7.0320e-06	8.1057e-05
6	2.1568e-05	2.1520e-05	3.2405e-05
9	1.8614e-05	1.8599e-05	2.0255e-05
15	5.5718e-07	5.4232e-07	5.8207e-07
20	5.6044e-08	5.5874e-08	6.1423e-08
25	1.0936e-09	9.9763e-10	5.9386e-10
30	3.3046e-11	3.2478e-11	4.0949e-11
35	6.4170e-13	6.3739e-13	1.2125e-12

| Show Table

DownLoad: CSV

5.1.4. Case 4:

In this case, we use a much shorter time interval, i.e., $ [0, T] = [0, 0.1] $, and compare the time-averages for this time interval to those for the time intervals $ [0, T] = [0, 1], [0, T] = [0, 10], $ and $ [0, T] = [0,100] $ (Case 3). We use two different time step values to generate the snapshots, but the same (i.e., $ n = 5000 $) equally spaced subintervals in the composite trapezoidal rule utilized in the evaluation of the time average $ \frac{1}{T}\int_0^{T}\mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $. In Tables 8–9, we list the time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for two different time step values and different $ m $ values. We emphasize that, this time, the time-average of $ \mathrm{e}_m(u) $ is negative for some $ m $ values. Furthermore, the magnitude of the time-average of $ \mathcal{E}_m(u) $ is larger than the magnitude of the time-average of $ \mathrm{e}_m(u) $. This is in stark contrast with the previous cases.

Table 8. Case 4: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ \Delta t = 2*10^{-5} $, different $ m $ values, $ d = 18 $, and $ 5000 $ equally spaced subintervals used in the composite trapezoidal rule.

$ m $	$ \frac{1}{0.1} \int_{0}^{0.1} \mathrm{e}_m(u(s)) ds $	$ \frac{1}{0.1} \int_{0}^{0.1} \mathcal{E}_m(u(s)) ds $
3	-6.8687e-04	2.7378e-05
5	-2.6333e-05	1.7795e-05
7	-9.1458e-07	4.2432e-06
9	-9.8188e-09	4.1220e-07
13	2.3800e-10	1.5749e-09
15	8.5423e-12	4.8043e-11

| Show Table

DownLoad: CSV

Table 9. Case 4: Time-averages of $ \mathrm{e}_m(u) $ and $ \mathcal{E}_{m}(u) $ for $ \Delta t = 10^{-5} $, different $ m $ values, $ d = 18 $, and $ 5000 $ equally spaced subintervals used in the composite trapezoidal rule.

$ m $	$ \frac{1}{0.1} \int_{0}^{0.1} \mathrm{e}_m(u(s)) ds $	$ \frac{1}{0.1} \int_{0}^{0.1} \mathcal{E}_m(u(s)) ds $
3	-6.8807e-04	2.7338e-05
5	-2.6447e-05	1.7812e-05
7	-9.1691e-07	4.2755e-06
9	-9.5609e-09	4.1694e-07
13	2.5287e-10	1.5941e-09
15	1.0725e-11	4.7030e-11

| Show Table

DownLoad: CSV

6. Conclusions

In this preliminary study, we investigated theoretically and numerically the time-average of the exchange of energy among modes of reduced order models (ROMs) of fluid flows. In particular, we were interested in the statistical equilibrium problem, and especially in the long-time averaging of the ROM solutions. The main goal of the paper was to deduce the possible forward and backward average transfer of the energy among ROM basis functions (modes). We considered two types of ROM modes: Eigenfunctions of the Stokes operator and proper orthogonal decomposition (POD) modes. In Theorem 4.1 and Theorem 4.2, we proved analytical results for both types of ROM modes and we highlighted the differences between them, especially those stemming from the lack of orthogonality of the gradients of the POD basis functions.

In Section 5, we performed a preliminary numerical investigation aiming at determining whether the time-average energy exchange between POD modes (i.e., $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $) in Theorem 4.1 is positive. To this end, we used the one-dimensional Burgers equation as a mathematical model. We utilized a piecewise linear FE spatial discretization and a backward Euler temporal discretization. To compute the time-averages, we used the composite trapezoidal rule. We tested different time steps, different number of subintervals in the composite trapezoidal rule, and, most importantly, different time intervals, to ensure that the computed quantities are indeed approximations of the time-averages and not numerical artifacts. The main conclusion of our numerical study is that, for long enough time intervals (i.e., time intervals longer than $ [0, T] = [0, 10] $), the time-average $ \frac{1}{T}\int_0^{T} \mathrm{e}_m(u(s)) + \mathcal{E}_{m}(u(s)) \, ds $ in (5.1) is positive. Furthermore, the magnitude of the time-average of $ \mathcal{E}_m(u) $ is much lower than the magnitude of the time-average of $ \mathrm{e}_m(u) $.

There are several research directions that we plan to pursue. Probably the most important one is the numerical investigation of the theoretical results in three-dimensional, high Reynolds number flows, which could shed new light on the energy transfer among ROM modes. A related, but different numerical investigation was performed in ^[7].

Acknowledgments

Luigi C. Berselli and Roger Lewandowski for the research that led to the present paper were partially supported by a grant of the group GNAMPA of INdAM and by the project the University of Pisa within the grant PRA_2018_52 UNIPI Energy and regularity: New techniques for classical PDE problems. Traian Iliescu and Birgul Koc gratefully acknowledge the support provided by the NSF DMS-1821145 grant. We thank the two reviewers for the comments and suggestions that improved the paper.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	Makhalanyane TP, Valverde A, Gunnigle E, et al. (2015) Microbial ecology of hot desert edaphic systems. FEMS Microbiol Rev 39: 203–222. doi: 10.1093/femsre/fuu011
[2]	Sterflinger K, Tesei D, Zakharova K (2012) Fungi in hot and cold deserts with particular reference to microcolonial fungi. Fungal Ecol 5: 453–462. doi: 10.1016/j.funeco.2011.12.007
[3]	Atlas of Israel (1985) 3 Eds., Tel Aviv: Surveys of Israel.
[4]	Grishkan I, Nevo E (2010) Soil microfungi of the Negev Desert, Israel-adaptive strategies to climatic stress, In: Veress B, Szigethy J, Editors, Horizons in Earth Science Research, New York: Nova Science Publishers Inc., 313–333.
[5]	Abdullah SK, Al-Khesraji TO, Al-Edany TY (1986) Soil mycoflora of the Southern Desert of Iraq. Sydowia 39: 8–16.
[6]	Bates ST, Nash TH III, Garcia-Pichel F (2012) Patterns of diversity for fungal assemblages of biological soil crusts from the southwestern United States. Mycologia 104: 353–366. doi: 10.3852/11-232
[7]	Ciccarone C, Rambelli A (1998) A study on micro-fungi in arid areas. Notes on stress-tolerant fungi. Plant Biosyst 132: 17–20.
[8]	Conley CA, Ishkhanova G, McKay CP, et al. (2006) A preliminary survey of non-lichenized fungi cultured from the hyperarid Atacama Desert of Chile. Astrobiology 6: 521–526. doi: 10.1089/ast.2006.6.521
[9]	Halwagy R, Moustafa AF, Kamel S (1982) Ecology of the soil mycoflora in the desert of Kuwait. J Arid Environ 5: 109–125.
[10]	Mouchacca J (1993) Thermophilic fungi in desert soils: a neglected extreme environment, In: Allsopp D, Colwell RR, Hawksworth DL, Editors, Microbial Diversity and Ecosystem Function, Wallingford: CAB International, 265–288.
[11]	Mulder JL, El-Hendawy H (1999) Microfungi under stress in Kuwait's coastal saline depressions. Kuwait J Sci Eng 26: 157–172.
[12]	Ranzoni FV (1968) Fungi isolated in culture from soils of the Sonoran Desert. Mycologia 60: 356–371. doi: 10.2307/3757166
[13]	Romero-Olivares AL, Baptista-Rosas RC, Escalante AE, et al. (2013) Distribution patterns of Dikarya in arid and semiarid soils of Baja California, Mexico. Fungal Ecol 6: 92–101. doi: 10.1016/j.funeco.2012.09.004
[14]	Zhang T, Jia RL, Yu LY (2016) Diversity and distribution of soil fungal communities associated with biological soil crusts in the southeastern Tengger Desert (China) as revealed by 454 pyrosequencing. Fungal Ecol 23: 156–163. doi: 10.1016/j.funeco.2016.08.004
[15]	Zak J (2005) Fungal communities of desert ecosystems: links to climate change, In: Dighton J, White Jr. JF, Oudemans P, Editors, The Fungal Community, its Organization and Role in the Ecosystem, Boca Raton, FL: CRC Press, 659–681.
[16]	Wicklow DT (1981) Biogeography and conidial fungi, In: Cole GT, Kendrick B, Editors, Biology of Conidial Fungi, New York: Academic Press, 417–447.
[17]	Zak J, Sinsabaugh R, MacKay WP (1995) Windows of opportunity in desert ecosystems: their applications to fungal community development. Can J Bot 73: S1407–S1414.
[18]	Grishkan I, Beharav A, Kirzhner V, et al. (2007) Adaptive spatiotemporal distribution of soil microfungi in 'Evolution Canyon' III, Nahal Shaharut, extreme Southern Negev desert, Israel. Biol J Linn Soc 90: 263–277. doi: 10.1111/j.1095-8312.2007.00722.x
[19]	Grishkan I, Kidron GJ (2013) Biocrust-inhabiting cultured microfungi along a dune catena in the western Negev Desert, Israel. Eur J Soil Biol 56: 107–114. doi: 10.1016/j.ejsobi.2013.03.005
[20]	Grishkan I, Kidron GJ (2016) Vertical divergence of microfungal communities through the depth in different soil formations at Nahal Nizzana, western Negev desert, Israel. Geomicrobiol J 7: 564–577.
[21]	Grishkan I, Nevo E (2010) Spatiotemporal distribution of soil microfungi in the Makhtesh Ramon area, central Negev desert, Israel. Fungal Ecol 3: 326–337. doi: 10.1016/j.funeco.2010.01.003
[22]	Grishkan I, Zaady E, Nevo E (2006) Soil crust microfungi along a southward rainfall aridity gradient in the Negev desert, Israel. Eur J Soil Biol 42: 33–42. doi: 10.1016/j.ejsobi.2005.09.014
[23]	Goldreich Y, Karni O (2001) Climate and precipitation regime in the Arava Valley, Israel. Isr J Earth Sci 50: 53–59. doi: 10.1560/1V61-FPGF-Y5VK-ADAG
[24]	Domsch KH, Gams W, Anderson TH (2007) Compendium of Soil Fungi, 2 Eds., New York: Academic Press.
[25]	Noy-Meir I (1973) Desert ecosystems: environment and producers. Ann Rev Ecol Syst 4: 25–51. doi: 10.1146/annurev.es.04.110173.000325
[26]	Ginat H, Shlomi Y, Batarseh S, et al. (2011) Reduction in precipitation levels in the Arava Valley (southern Israel and Jordan), 1949–2009. J Dead-Sea Arava Res 1: 1–7.
[27]	Singer A (2007) The Soils of Israel, Berlin Heidelberg: Springer-Verlag.
[28]	Danin A (1983) Desert Vegetation of Israel and Sinai, Jerusalem: Cana Publishing House.
[29]	Rowell DL (1994) Soil Science: Methods and Applications, London: Longman Scientific and Technical.
[30]	Davet P, Rouxel F (2000) Detection and Isolation of Soil Microfungi, Enfield (NH), Plymouth: Science Publisher Inc.
[31]	Jeewon R, Hyde KD (2007) Detection and diversity of fungi from environmental samples: traditional versus molecular approaches, In: Varma A, Oelmuller R, Editors, Advanced Techniques in Soil Microbiology, Berlin Heidelberg: Springer-Verlag, 1–15.
[32]	Bills GF, Christensen M, Powell M, et al. (2004) Saprobic soil microfungi, In: Mueller GM, Bills GF, Foster MS, Editors, Biodiversity of Microfungi. Inventory and Monitoring Methods, San Francisco: Elsevier Academic Press, 271–302.
[33]	Kirk PM, The CABI Bioscience and CBS Database of Fungal Names. Available from www.indexfungorum.org.
[34]	Krebs CJ (1999) Ecological Methodology, Chicago: Addison Wesley Longman.
[35]	Smith JL, Doran JW (1996) Measurement and use of pH and electrical conductivity for soil quality analysis, In: Doran JW, Jones AJ, Editors, Methods for Assessing Soil Quality, Madison: Soil Science Society of America, 169–185.
[36]	Tripathi BM, Moroenyane I, Sherman C, et al. (2017) Trends in taxonomic and functional composition of soil microbiome along a precipitation gradient in Israel. Microb Ecol 74: 168–176. doi: 10.1007/s00248-017-0931-0
[37]	Valente DSM, Queiroz DM, Pinto C, et al. (2012) Definition of management zones in coffee production fields based on apparent soil electrical conductivity. Sci Agr 69: 173–179. doi: 10.1590/S0103-90162012000300001
[38]	Eigenberg RA, Doran JW, Nienaber JA, et al. (2002) Electrical conductivity monitoring of soil condition and available N with animal manure and a cover crop. Agr Ecosyst Environ 88: 183–193. doi: 10.1016/S0167-8809(01)00256-0
[39]	Grishkan I (2011) Ecological stress: Melanization as a response in fungi to radiation, In: Horikoshi K, Antranikian G, Bull A, et al., Editors, Extremophiles Handbook, Berlin Heidelberg: Springer-Verlag, 1135–1145.
[40]	Christensen M (1981) Species diversity and dominance in fungal community, In: Carroll GW, Wicklow DT, Editors, The Fungal Community, Its Organization and Role in the Ecosystem, New York: Marcell Dekker, 201–232.
[41]	Grishkan I, Jia RL, Kidron GJ, et al. (2015) Cultivable microfungal communities inhabiting biological soil crusts in the Tengger Desert, China. Pedosphere 25: 351–363. doi: 10.1016/S1002-0160(15)30003-5
[42]	Bates ST, Nash TH III, Sweat KG, et al. (2010) Fungal communities of lichen-dominated biological soil crusts: diversity, relative microbial biomass, and their relationship to disturbance and crust cover. J Arid Environ 74: 1192–1199. doi: 10.1016/j.jaridenv.2010.05.033
[43]	Porras-Alfaro A, Herrera J, Navig DO, et al. (2011) Diversity and distribution of soil fungal communities in a semiarid grassland. Mycologia 103: 10–21. doi: 10.3852/09-297
[44]	Ellis MB, Ellis JP (1997) Microfungi on Land Plant. An identification handbook, Slough, Berkshire: Richmond Publishing Co. Ltd.
[45]	Sutton BC (1980) The Coelomycetes. Kew, Surrey, England: Commonwealth Mycological Institute.
[46]	Grishkan I, Nevo E, Wasser SP (2003) Micromycete diversity in the hypersaline Dead Sea coastal area (Israel). Mycol Progr 2: 19–28. doi: 10.1007/s11557-006-0040-9
[47]	Bokhary HA (1998) Mycoflora of desert sand dunes of Riyadh Region, Saudi Arabia. J King Saud Univ 10: 15–29.
[48]	Oliveira LG, Cavalcanti MAQ, Fernandes MJS, et al. (2013) Diversity of filamentous fungi isolated from the soil in the semiarid area, Pernambuco, Brazil. J Arid Environ 95: 49–53. doi: 10.1016/j.jaridenv.2013.03.007
[49]	Powell AJ, Parchert KJ, Bustamante JM, et al. (2012) Thermophilic fungi in an arid land ecosystem. Mycologia 104: 813–825. doi: 10.3852/11-298

This article has been cited by:

1.	Ping Guo, Wei Tian, Huimin Li, Guangmin Zhang, Jianhui Li, Global characteristics and trends of research on construction dust: based on bibliometric and visualized analysis, 2020, 27, 0944-1344, 37773, 10.1007/s11356-020-09723-y
2.	Luigi C. Berselli, 2021, 9780128219546, 255, 10.1016/B978-0-12-821954-6.00011-9
3.	Shady E. Ahmed, Suraj Pawar, Omer San, Adil Rasheed, Traian Iliescu, Bernd R. Noack, On closures for reduced order models—A spectrum of first-principle to machine-learned avenues, 2021, 33, 1070-6631, 091301, 10.1063/5.0061577
4.	Huiping Wang, Yi Wang, Estimating CO2 emissions using a fractional grey Bernoulli model with time power term, 2022, 29, 0944-1344, 47050, 10.1007/s11356-022-18803-0
5.	Thierry Delisle, Peter A. Buhr, Advanced control‐flow and concurrency in C∀, 2021, 51, 0038-0644, 1005, 10.1002/spe.2925
6.	Junqi Zhu, Haotian Zheng, Li Yang, Shanshan Li, Liyan Sun, Jichao Geng, Evaluation of deep coal and gas outburst based on RS-GA-BP, 2023, 115, 0921-030X, 2531, 10.1007/s11069-022-05652-w
7.	Rahmat Ali Khan, Shaista Gul, Fahd Jarad, Hasib Khan, Existence results for a general class of sequential hybrid fractional differential equations, 2021, 2021, 1687-1847, 10.1186/s13662-021-03444-3
8.	Hongjun Han, Liyuan Wang, Diffraction of SV waves by a surface effective elliptic cavity, 2022, 17, 1559-3959, 65, 10.2140/jomms.2022.17.65
9.	Xingyue Deng, Jiuyuan Huo, Ling Wu, 2022, Chapter 20, 978-981-19-5208-1, 298, 10.1007/978-981-19-5209-8_20
10.	Pradeep Kumar Dabla, Kamal Upreti, Divakar Singh, Anju Singh, Jitender Sharma, Aashima Dabas, Damien Gruson, Bernard Gouget, Sergio Bernardini, Evgenija Homsak, Sanja Stankovic, Target association rule mining to explore novel paediatric illness patterns in emergency settings, 2022, 82, 0036-5513, 595, 10.1080/00365513.2022.2148121
11.	Mingsheng Yang, Jie Yang, Weihong Han, Jiawei Zhang, 2024, Chapter 24, 978-981-97-4518-0, 336, 10.1007/978-981-97-4519-7_24

Reader Comments

Your name:*

Email:*
© 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Microbiology

2.7 7.0

Metrics

Article views(5103) PDF downloads(757) Cited by(5)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(5) / Tables(4)

AIMS Microbiology

Spatiotemporal variations in soil cultivable mycobiota at the Arava desert (Israel) along latitudinal and elevational gradients

Related Papers:

Abstract

1. Introduction

2. Reduced order modeling

2.1. On the spectral decomposition

3. Preliminaries on long-time averages

4. Average transfer of energy at equilibrium

4.1. The POD case

4.2. The spectral case

5. Numerical results

5.1. Numerical results with step function initial condition

5.1.1. Case 1:

5.1.2. Case 2:

5.1.3. Case 3:

5.1.4. Case 4:

6. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

AIMS Microbiology

Spatiotemporal variations in soil cultivable mycobiota at the Arava desert (Israel) along latitudinal and elevational gradients

Related Papers:

Abstract

1. Introduction

2. Reduced order modeling

2.1. On the spectral decomposition

3. Preliminaries on long-time averages

4. Average transfer of energy at equilibrium

4.1. The POD case

4.2. The spectral case

5. Numerical results

5.1. Numerical results with step function initial condition

5.1.1. Case 1:

5.1.2. Case 2:

5.1.3. Case 3:

5.1.4. Case 4:

6. Conclusions

Acknowledgments

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog