manus_continuum_granular1

manuscript files for first continuum-till paper
git clone git://src.adamsgaard.dk/manus_continuum_granular1
Log | Files | Refs

commit 255caafff86199c794ab3806cd2239eff32a331c
parent 4b45b2c20ee90c421e248df474d70892c3490039
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Wed,  4 Sep 2019 13:19:39 +0200

Remove fig7 of normalized shear velocities now that fig6 is rate controlled

Diffstat:
Mcontinuum-granular-manuscript1.tex | 16++++------------
Mexperiments/fig1.pdf | 0
Mexperiments/fig2.pdf | 0
Mexperiments/fig3.pdf | 0
Mexperiments/fig6.pdf | 0
5 files changed, 4 insertions(+), 12 deletions(-)

diff --git a/continuum-granular-manuscript1.tex b/continuum-granular-manuscript1.tex @@ -119,7 +119,7 @@ Future research will investigate how wide grain-size distributions affect strain \label{sub:numerical_solution_procedure} The presented formulation is applicable to any spatial dimensionality. For the purposes of this study we apply it in a 1D spatial reference system. -The axis $z$ is pointed upwards with a domain length of $L_z$. +Axis $x$ is along flow and axis $z$ is pointed upwards with a domain thickness of $L_z$. Shear deformation is restricted to occur in horizontal (x) shear zones. We assign depth coordinates $z_i$ and fluidity $g_i$ to a regular grid with ghost nodes and cell spacing $\Delta z$. The normal stress is assumed to increase with depth due to lithostatic pressure from the overburden ($\sigma_\text{n}(z) = \int^{z'=L_z}_{z'=z} \rho_\text{s} \phi G dz' + \sigma_\text{n,t}$), where G is the magnitude of gravitational acceleration and $\sigma_\text{n,t}$ is the normal stress applied on the top of the domain. @@ -166,6 +166,9 @@ We use a representative grain size of $d = 0.04$ m, a grain material density of Dimensionless material parameters $A$ and $b$ from Eq.~\ref{eq:g_local} and~\ref{eq:cooperativity} are 0.4 and 0.9377, respectively. These values are constrained from experiments on glass beads \citep{Damsgaard2013, Henann2016}. +Importantly, the resultant shear velocities are in this setup not limited by anything but sediment kinematics. +The simulated velocities are for the most part far greater than any glacial setting, where horizontal stresses keep ice masses in place over weak beds. + For the first experiment with variable water pressure, we apply a water-pressure forcing amplitude of 50 kPa that modulates effective stress at the top around 100 kPa (Fig.~\ref{fig:stick_slip}). @@ -281,17 +284,6 @@ The shear velocities during the first cycle ($t<1$ d) is slightly different from \end{center} \end{figure*} -\begin{figure*}[htbp] - \begin{center} - \includegraphics[width=15.0cm]{experiments/fig7.pdf} - \caption{\label{fig:stick_slip_depth_normalized}% - Pore-pressure diffusion and normalized strain distribution with depth with a sinusoidal water-pressure forcing from the top (Fig.~\ref{fig:stick_slip}). - The forcing has a daily periodocity, and plot lines are one hour in simulation time apart. - The horizontal magenta line marks skin depth from Eq.~\ref{eq:skin_depth}. - } - \end{center} -\end{figure*} - \begin{figure}[htbp] \begin{center} \includegraphics[width=7.5cm]{experiments/fig8.pdf} diff --git a/experiments/fig1.pdf b/experiments/fig1.pdf Binary files differ. diff --git a/experiments/fig2.pdf b/experiments/fig2.pdf Binary files differ. diff --git a/experiments/fig3.pdf b/experiments/fig3.pdf Binary files differ. diff --git a/experiments/fig6.pdf b/experiments/fig6.pdf Binary files differ.