# Title: Air-Water capillary wave
#
# Description:
#
# Same test as before but with density and viscosity ratio
# corresponding to an air/water interface.
#
# \begin{table}[htbp]
# \caption{\label{convergence}Convergence of the relative error between the analytical
# solution and simulation results.}
# \begin{center}
# \begin{tabular}{|l|ccccc|} \hline
# Method & $8^2$ & $16^2$ & $32^2$ & $64^2$ & $128^2$ \\ \hline
# \input{convergence.tex} & 0.00313 \\
# \end{tabular}
# \end{center}
# \end{table}
#
# \begin{figure}[htbp]
# \caption{\label{amplitude}Evolution of the amplitude of the capillary wave as a
# function of non-dimensional time $\tau=\omega_0 t$.}
# \begin{center}
# \includegraphics[width=\hsize]{amplitude.eps}
# \end{center}
# \end{figure}
#
# Author: St\'ephane Popinet
# Command: sh ../capwave.sh air-water.gfs
# Version: 1.2.0
# Required files: convergence.ref prosperetti
# Generated files: convergence.tex amplitude.eps
#
# Theoretical solution generated using:
#  ~/local/src/laplace/oscillations planar 1 1 0.0182571749236 0.273038508 1 1 0.0012 0.01 0 0 | awk '{print $1*15.7402, ($2 > 0. ? $2 : -$2)}' > prosperetti
#
3 5 GfsSimulation GfsBox GfsGEdge {} {
  Time { end = 1.58928694288774963184 }
  ApproxProjectionParams { tolerance = 1e-6 }
  ProjectionParams { tolerance = 1e-6 }
  Refine floor(LEVEL + 1 - (LEVEL - 2)*fabs(y)/1.5)
  VariableTracerVOF T
  VariableFiltered T1 T 1
  VariableCurvature K T
  SourceTension T 1 K
  VariablePosition Y T y
  Global {
      #define VAR(T,min,max)   (min + CLAMP(T,0,1)*(max - min))
      #define RHO(T)            VAR(T, 1.2/1000., 1.)
      #define MU(T)             VAR(T, 1.8e-5/1.003e-3, 1.)
  }
  PhysicalParams { alpha = 1./RHO(T1) }
  SourceViscosity 0.0182571749236*MU(T1)
  InitFraction T (y - 0.01*cos (2.*M_PI*x))
  OutputScalarNorm { step = 0.00198785108553814829 } {
      awk '{printf ("%g %g\n", $3*15.7402, $9); fflush(stdout); }' > wave-LEVEL
  } { v = (T > 0. && T < 1. ? Y : 0.) }
}
GfsBox {}
GfsBox {}
GfsBox {}
1 1 right
2 2 right
3 3 right
1 2 top
1 3 bottom