# gerris2D -DLEVEL=11 -DRE=100 circ-visc.gfs # gfsview end-11-100.gfs isolines.gfv Define A0 0.5 Define U0 1. 1 0 GfsAxi GfsBox GfsGEdge {} { Time { end = 100 } PhysicalParams { L = 50 } AdvectionParams { gc = 1 } Refine 4 Refine (LEVEL + 1./50.*(x*x + y*y)*(4. - LEVEL)) Solid (ellipse (0., 0., A0, A0)) SourceViscosity 1./RE Init {} { U = U0 } AdaptGradient { istep = 1 } { cmax = 5e-2 maxlevel = LEVEL } U AdaptGradient { istep = 1 } { cmax = 5e-2 maxlevel = LEVEL } V AdaptFunction { istep = 1 } { cmax = 1e-2 maxlevel = LEVEL } { return (fabs(dx("U"))+fabs(dy("U")))/fabs(U)*ftt_cell_size (cell); } EventStop { step = 0.1 } U 1e-3 DU OutputSimulation { start = end } end-LEVEL-RE.gfs } GfsBox { left = Boundary { BcDirichlet U U0 } right = BoundaryOutflow bottom = Boundary }