Skip to content
Snippets Groups Projects
Commit 1d2dcbc6 authored by Pierre Augier's avatar Pierre Augier
Browse files

2022strat_turb_toro: Add def Ivelo and Idiss

parent 80f36124
No related branches found
No related tags found
1 merge request!82022strat_turb_toro...
......@@ -398,8 +398,9 @@
\includegraphics[width=0.48\textwidth]{%
../tmp/fig_isotropy_velo_vs_R}
}
\caption{Large scale isotropy coefficient $\Ivelo$. \label{fig:large-scale-isotropy}}
\caption{Large scale isotropy coefficient $\Ivelo$ versus $F_h$ in (a) and $\R$ in (b).
\label{fig:large-scale-isotropy}}
\end{figure}
Figure~\ref{fig:large-scale-isotropy} ...
......@@ -402,7 +403,11 @@
\end{figure}
Figure~\ref{fig:large-scale-isotropy} ...
\begin{equation}
\Ivelo \equiv \frac{3 \mean{{v_z}^2}}{\mean{{v_x}^2 + {v_y}^2 + {v_z}^2}}
\end{equation}
%% Small scale isotropy
\begin{figure}
......@@ -412,8 +417,9 @@
\includegraphics[width=0.48\textwidth]{%
../tmp/fig_isotropy_diss_vs_R}
}
\caption{Small scale isotropy coefficient $\Idiss$. \label{fig:small-scale-isotropy}}
\caption{Small scale isotropy coefficient $\Idiss$ versus $F_h$ in (a) and $\R$ in (b).
\label{fig:small-scale-isotropy}}
\end{figure}
Figure~\ref{fig:small-scale-isotropy} ...
......@@ -416,7 +422,15 @@
\end{figure}
Figure~\ref{fig:small-scale-isotropy} ...
\begin{equation}
\Idiss \equiv \frac{1 - r_{\text{diss}}}{1 - r_{\text{diss,\,iso}}},
\end{equation}
where $r_{\text{diss}}$ is the ratio of the kinetic energy dissipation associated with
vertical gradients over the total kinetic energy dissipation $\epsK$, and
$r_{\text{diss,\,iso}} = 1/3$ is the value for this ratio for isotropic turbulence.
%% Isotropy coefficient: summary
\begin{figure}
......
......@@ -20,6 +20,7 @@
xs = np.linspace(1e-2, 4e-2, 2)
ax.plot(xs, 1.8e1 * xs**1, "b")
ax.text(1.5e-2, 0.21, "${F_h}^1$")
xs = np.linspace(2e-2, 9.0e-2, 2)
ax.plot(xs, 3.5e0 * xs**0.5, "r")
......@@ -23,6 +24,7 @@
xs = np.linspace(2e-2, 9.0e-2, 2)
ax.plot(xs, 3.5e0 * xs**0.5, "r")
ax.text(2e-2, 0.75, "${F_h}^{0.5}$")
ax.set_xlabel("$F_h$")
ax.set_ylabel(r"$L_v / L_h$")
......
......@@ -21,6 +21,7 @@
xs = np.linspace(1e-2, 5.0e-2, 2)
ax.plot(xs, 1e1 * xs**1, "b")
ax.text(1.8e-2, 0.15, "${F_h}^1$")
xs = np.linspace(1e-2, 1.0e-1, 2)
ax.plot(xs, 2e0 * xs**0.5, "r")
......@@ -24,6 +25,7 @@
xs = np.linspace(1e-2, 1.0e-1, 2)
ax.plot(xs, 2e0 * xs**0.5, "r")
ax.text(1.2e-2, 0.35, "${F_h}^{0.5}$")
ax.set_xlabel("$F_h$")
ax.set_ylabel(r"$U_v / U_h$")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment