GfsAdaptGradient

From Gerris

GfsAdaptGradient uses a "cell cost" defined as the norm of the local gradient of a given variable or function multiplied by the cell size.

The syntax in parameter files is:

[ GfsAdapt ] [ GfsFunction ]

Examples

• B\'enard--von K\'arm\'an Vortex Street for flow around a cylinder at Re=160

  AdaptGradient { istep = 1 } { maxlevel = 6 cmax = 1e-2 } T

• Vortex street around a "heated" cylinder

  AdaptGradient { istep = 1 } { maxlevel = 6 cmax = 1e-2 } T

• Parallel simulation on four processors

  AdaptGradient { istep = 1 } { maxlevel = 6 cmax = 1e-2 } T

• Rayleigh-Taylor instability

  AdaptGradient { istep = 1 } { maxlevel = 7 cmax = 1e-2 } T

• Boussinesq flow generated by a heated cylinder

  AdaptGradient { istep = 1 } { maxlevel = 8 cmax = 5e-2 } T

• Collapse of a column of grains

    AdaptGradient { istep = 1 } {
	cmax = 0
	maxlevel = LEVEL
    } T

• Small amplitude solitary wave interacting with a parabolic hump

    AdaptGradient { istep = 1 } { 
	cmax = 1e-4
	cfactor = 2
	maxlevel = 8
	minlevel = 6
    } (P + Zb)

• Shock reflection by a circular cylinder

    AdaptGradient { istep = 1 } {
        cmax = 0.1
        cfactor = 2
        maxlevel = 9
    } P

• Tsunami runup onto a complex three-dimensional beach

    AdaptGradient { istep = 1 } {
        cmax = 1e-3
        cfactor = 2
        maxlevel = 9
        minlevel = 6
    } (P < DRY ? 0. : P + Zb)

• The 2004 Indian Ocean tsunami

    AdaptGradient { start = 212 istep = 1 end = 272 } {
	cmax = 0.05 cfactor = 2
	minlevel = 5 maxlevel = LEVEL 
    } D

    AdaptGradient { start = 528 istep = 1 end = 588 } {
	cmax = 0.05 cfactor = 2
	minlevel = 5 maxlevel = LEVEL 
    } D

    AdaptGradient { start = 853 istep = 1 end = 913 } {
	cmax = 0.05 cfactor = 2
	minlevel = 5 maxlevel = LEVEL 
    } D	

    AdaptGradient { start = 1213 istep = 1 end = 1273 } {
	cmax = 0.05 cfactor = 2
	minlevel = 5 maxlevel = LEVEL 
    } D	

    AdaptGradient { istep = 1 } {
	cmax = 0.05 cfactor = 2
	minlevel = 5 maxlevel = LEVEL 
    } (P < DRY ? 0. : P + Zb)

• "Garden sprinkler effect" in wave model

    AdaptGradient { istep = 1 } { cmax = 0.04 minlevel = MINLEVEL maxlevel = 6 } Hs

• Time-reversed VOF advection in a shear flow

    AdaptGradient { istart = 1 istep = 1 } { cmax = 0 maxlevel = 8 } T

• Time-reversed advection of a VOF concentration

    AdaptGradient { istart = 1 istep = 1 } { cmax = 0 maxlevel = LEVEL } T

• Conservation of diffusive tracer

    AdaptGradient { istep = 1 } { minlevel = 3 maxlevel = 5 cmax = 1e-2 } T

• Viscous flow past a sphere

    AdaptGradient { istep = 1 } { cmax = 5e-2 maxlevel = LEVEL } U

    AdaptGradient { istep = 1 } { cmax = 5e-2 maxlevel = LEVEL } V

• Mass conservation

    AdaptGradient { istep = 1 } { cmax = 1e-3 minlevel = 4 maxlevel = (x < 0.25 ? 6 : 7) } T1

• Mass conservation with solid boundary

    AdaptGradient { istep = 1 } { cmax = 1e-3 minlevel = 4 maxlevel = (x < 0. ? 7 : 8) } T1

• Momentum conservation for large density ratios

    AdaptGradient { istep = 1 } { cmax = 1e-3 maxlevel = level } T

• Circular droplet in equilibrium

  AdaptGradient { istep = 1 } { cmax = 1e-6 maxlevel = LEVEL } T

• Axisymmetric spherical droplet in equilibrium

  AdaptGradient { istep = 1 } { cmax = 1e-6 maxlevel = LEVEL } T

• Scalings for Plateau--Rayleigh pinchoff

   AdaptGradient { istep = 1 } { maxlevel = 6 cmax = 0 } T

• Circular dam break on a ``cubed sphere''

    AdaptGradient { istep = 1 } { cmax = 1e-2 maxlevel = 8 } P

• Advection of a cosine bell around the sphere

  AdaptGradient { istep = 1 } { cmax = 1e-4 maxlevel = LEVEL } T

• Charge relaxation in an axisymmetric insulated conducting column

    AdaptGradient { istep = 1 } { cmax = 1e-4 minlevel = 5 maxlevel = LEVEL } T

• Charge relaxation in a planar cross-section

    AdaptGradient { istep = 1 } { cmax = 1e-4 minlevel = 5 maxlevel = LEVEL + 1 } T

• Equilibrium of a droplet suspended in an electric field

    AdaptGradient { istep = 1 } { cmax = 1e-4 minlevel = 4 maxlevel = 7 } T

• Gouy-Chapman Debye layer

    AdaptGradient { istep = 5 } { cmax = 0.02 minlevel = 3 maxlevel = 5 } Cneg