GfsEventStop

From Gerris

GfsEventStop is used to stop the simulation when steady-state is reached. This is done by monitoring the maximum change in a given variable at every locations between successive calls to GfsEventStop. If the absolute value of this change is smaller than or equal to a user-defined threshold the simulation is stopped.

The syntax in parameter files is:

[ GfsEvent ] V THRESHOLD ( DV )

where V is a variable, THRESHOLD a number and DV an optional variable name. If DV is given, the corresponding field contains the difference in V between successive calls.

Examples

• Lunar tides in Cook Strait, New Zealand

    EventStop { start = 100000 istep = 100 } A0 0.025

• Viscous flow past a sphere

    EventStop { step = 0.1 } U 1e-3 DU

• Boundary layer on a rotating disk

    EventStop { step = 0.1 } U 1e-3 DU

• Lid-driven cavity at Re=1000

  EventStop { istep = 10 } U 1e-4 DU

• Lid-driven cavity at Re=1000 (explicit scheme)

  EventStop { istep = 10 } U 1e-4 DU

• Lid-driven cavity with a non-uniform metric

  EventStop { istep = 10 } U 1e-4 DU

• Lid-driven cavity on an anisotropic mesh

  EventStop { istep = 10 } U 1e-4 DU

• Poiseuille flow

    EventStop { istep = 1 } U 1e-6 DU

• Bagnold flow of a granular material

    EventStop { istart = 10 istep = 10 } U 5e-6 DU

• Poiseuille flow with metric

    EventStop { istep = 1 } U 1e-6 DU

• Poiseuille flow with multilayer Saint-Venant

    EventStop { istep = 1 } U 1e-8 DU

• Creeping Couette flow of Generalised Newtonian fluids

  EventStop { istep = 1 } U 1e-4 DU

• Multi-layer Saint-Venant solver

    EventStop { step = 1 } U0 1e-9

• Circular droplet in equilibrium

  EventStop { istep = 10 } T DT

• Axisymmetric spherical droplet in equilibrium

  EventStop { istep = 10 } T DT

• Sessile drop

    EventStop { istep = 10 } K 1e-5 DK

• Transcritical flow over a bump

    EventStop { step = 1 } H 1e-5

• Transcritical flow with multiple layers

    EventStop { step = 1 } H 1e-3

• Creeping Couette flow between cylinders

    EventStop { step = 1e-2 } V 1e-5 DV

• Creeping Couette flow between eccentric cylinders

  EventStop { step = 1e-2 } U 1e-4 DU

• Flow between eccentric cylinders using bipolar coordinates

  EventStop { step = 1e-2 } U 1e-4 DU

• Flow between eccentric cylinders on a stretched grid

  EventStop { step = 1e-2 } U 1e-4 DU

• Charge relaxation in an axisymmetric insulated conducting column

    EventStop { istep = 10 } Ex 0.001

• Charge relaxation in a planar cross-section

    EventStop { istep = 10 } Ex 0.001

• Equilibrium of a droplet suspended in an electric field

    EventStop { istep = 10 } Rhoe 0.001 { relative = 1 }

• Gouy-Chapman Debye layer

    EventStop { istep = 10 } Phi 1e-4 DPhi

• Simple example of groundwater flow following Darcy's law

    EventStop { istep = 1 } U 1e-9

• Groundwater flow with piecewise constant permeability

    EventStop { istep = 1 } U 1e-9