sandbox/lopez/remove.c

    Test for withdrawal of droplets

    //#include "grid/cartesian.h"
#include "run.h"
#include "tag.h"
#include "src/droplet_stat.h"

scalar c[];

int main()
{
  origin (-0.5, -0.5);
  init_grid (64);
  run ();
}

    A set of circles of decreasing radius are distributed in the computational domain.

    A proper initialization would need the Vofi library to compute the correct volume fractions but we do a rough initizalization setting to 1 these circles.

    event init (i = 0) {  
  foreach()
    c[] = 0.;
  
  int nx = 3, ny = 2;
  double radius = 0.1;
  for (int j = 1; j <= ny; j++) 
    for (int i = 1; i <= nx; i++) {
      double xc = -0.5 + (double) i/(nx + 1); 
      double yc = -0.5 + (double) j/(ny + 1);
      radius *= 0.8;
      foreach()
	if(sq(x-xc)+sq(y-yc) < sq(radius))
	  c[] = 1.;
    }
}

event drop_remove (t = 0) {
  droplet_remove (c, 3, true); 
  output_field ({c}, stderr);    
}
    set pm3d map 
set size ratio -1
splot 'log' u 1:2:3
    The three largest droplets remain (three removed) (script)

    The three largest droplets remain (three removed) (script)