src/test/large-ns.c

    Large-amplitude standing wave

    See large.c for a more detailed description.

    #include "grid/multigrid.h"
    #include "navier-stokes/centered.h"
    #include "two-phase.h"
    #include "navier-stokes/conserving.h"
    #include "reduced.h"
    
    int main()
    {
      rho1 = 1./1000.;
      rho2 = 1.;
      G.y = -9.81;
      DT = 1e-3 [0,1];
      Y0 = - L0/2.;
      N = 256;
      run();
    }
    
    event init (i = 0)
    {
      double a = 0.07, k = 2.*pi;
      fraction (f, - (a*cos(k*x) - y));
    }
    
    event profiles (t = 0.1; t += 0.1; t <= 0.5)
    {
      output_facets (f, stderr);
    }