GfsOutputScalarNorm

From Gerris

GfsOutputScalarNorm is used to write the volume-weighted norms over the whole domain of a given scalar.

The statistics are written using the following formatting:

DESCRIPTION time: T first:  FIRST second:  SECOND infty:  INFTY

with:

DESCRIPTION

a description of the scalar field (without any spaces),

T

the physical time,

FIRST

the 1-norm i.e. the average of the absolute values,

SECOND

the 2-norm i.e. the root-mean-square norm,

INFTY

the infinity-norm i.e. the maximum absolute value.

The syntax in parameter files is:

[ GfsOutputScalar ]

Examples

• Collapse of a column of grains

    OutputScalarNorm { istep = 10 } X-H0-LEVEL { v = (T > 0.1 ? X : G_MAXDOUBLE) }

    OutputScalarNorm { istep = 10 } Y-H0-LEVEL { v = (T > 0.1 ? Y : G_MAXDOUBLE) }

    OutputScalarNorm { istep = 10 } Yc-H0-LEVEL { v = (x < 0.1 ? Y : G_MAXDOUBLE) }

• Savart--Plateau--Rayleigh instability of a water column

    OutputScalarNorm { istep = 1 } v { v = Velocity }

• Lunar tides in Cook Strait, New Zealand

    OutputScalarNorm { istep = 1 } p { v = P }

    OutputScalarNorm { istep = 1 } u { v = Velocity }

    OutputScalarNorm { istep = 10 } a0 { v = sqrt(A0*A0 + B0*B0) }

• Dam break on complex topography

    OutputScalarNorm { istep = 10 } u { v = (P > 0. ? U/P : 0.) }

• Tsunami runup onto a complex three-dimensional beach

    OutputScalarNorm { istep = 1 } u { v = Velocity }

• The 2004 Indian Ocean tsunami

    OutputScalarNorm { istep = 1 } u { v = Velocity }

    OutputScalarNorm { istep = 1 } U { v = U }

    OutputScalarNorm { istep = 1 } V { v = V }

    OutputScalarNorm { istep = 1 } hwet { v = Hwet }

• Conservation of diffusive tracer

    OutputScalarNorm { istep = 1 } diff { v = (T - Te) }

• Estimation of the numerical viscosity

  OutputScalarNorm { istep = 1 } divLEVEL { v = Divergence }

• Estimation of the numerical viscosity with refined box

  OutputScalarNorm { istep = 1 } divLEVEL { v = Divergence }

• Numerical viscosity for the skew-symmetric scheme

  OutputScalarNorm { istep = 1 } divLEVEL { v = Divergence }

• Numerical viscosity for the skew-symmetric scheme with refined box

  OutputScalarNorm { istep = 1 } div-LEVEL { v = Divergence }

• Numerical viscosity for vorticity/streamfunction formulation

  OutputScalarNorm { istep = 1 } divLEVEL { v = (dx("U") + dy("V")) }

• Lid-driven cavity at Re=1000

  OutputScalarNorm { istep = 10 } du { v = DU }

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

  OutputScalarNorm { istep = 10 } du { v = DU }

• Lid-driven cavity with a non-uniform metric

  OutputScalarNorm { istep = 10 } du { v = DU }

• Lid-driven cavity on an anisotropic mesh

  OutputScalarNorm { istep = 10 } du { v = DU }

• Bagnold flow of a granular material

    OutputScalarNorm { istep = 10 } du-LEVEL { v = DU }

• Creeping Couette flow of Generalised Newtonian fluids

  OutputScalarNorm { istep = 1 } du-MODEL { v = DU }

• Hydrostatic balance with solid boundaries and viscosity

    OutputScalarNorm { istep = 1 } v { v = V }

• Hydrostatic balance with quadratic pressure profile

    OutputScalarNorm { istep = 1 } v { v = V }

• Potential flow around a thin plate

  OutputScalarNorm { start = end } stdout { v = Velocity } 

• Circular droplet in equilibrium

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $9*sqrt(0.8) }' > La-LAPLACE-LEVEL
  } { v = Velocity }

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9 }' > E-LAPLACE-LEVEL
  } { v = (Tref - T) }

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9 }' > EK-LAPLACE-LEVEL
  } { v = (T > 0 && T < 1 ? K - 2.5 : 0) }

• Axisymmetric spherical droplet in equilibrium

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $9*sqrt(0.8); fflush (stdout); }' > La-LAPLACE-LEVEL
  } { v = Velocity }

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9; fflush (stdout); }' > E-LAPLACE-LEVEL
  } { v = (Tref - T) }

  OutputScalarNorm { istep = 1 } {
    awk '{ print MU*$3/(0.8*0.8), $5, $7, $9; fflush (stdout); }' > EK-LAPLACE-LEVEL
  } { v = (T > 0 && T < 1 ? (K - 5.)/2. : 0) }

• Planar capillary waves

  OutputScalarNorm { step = 3.04290519077e-3 } {
      awk '{printf ("%g %g\n", $3*11.1366559937, $9); fflush(stdout); }' > wave-LEVEL
  } { v = (T > 0. && T < 1. ? Y : 0.) }

• Air-Water capillary wave

  OutputScalarNorm { step = 0.00198785108553814829 } {
      awk '{printf ("%g %g\n", $3*15.7402, $9); fflush(stdout); }' > wave-LEVEL
  } { v = (T > 0. && T < 1. ? Y : 0.) }

• Fluids of different densities

  OutputScalarNorm { step = .00225584983639310905 } {
      awk '{printf ("%g %g\n", $3*15.016663878457, $9); fflush(stdout); }' > wave-LEVEL
  } { v = (T > 0. && T < 1. ? Y : 0.) }

• Pure gravity wave

  OutputScalarNorm { step = .00225584983639310905 } {
      awk '{printf ("%g %g\n", $3*16.032448313657, $9); fflush(stdout); }' > wave-LEVEL
  } { v = (T > 1e-6 && T < 1. - 1e-6 ? Y : 0.) }

• Scalings for Plateau--Rayleigh pinchoff

   OutputScalarNorm { istep = 1 } u { v = U }

• Sessile drop

    OutputScalarNorm { istep = 10 } v-ANGLE { v = Velocity }

• Lake-at-rest balance in an inclined domain with cut cells

    OutputScalarNorm { start = end } u { v = U }

• Lake-at-rest balance in an inclined domain with bipolar metric

    OutputScalarNorm { start = end } u { v = U }

• Creeping Couette flow between cylinders

    OutputScalarNorm { istep = 1 } dv { v = DV }

• Creeping Couette flow between eccentric cylinders

  OutputScalarNorm { istep = 1 } du { v = DU }

• Flow between eccentric cylinders using bipolar coordinates

  OutputScalarNorm { istep = 1 } du { v = DU }

• Flow between eccentric cylinders on a stretched grid

  OutputScalarNorm { istep = 1 } du { v = DU }

• Rossby--Haurwitz wave

    OutputScalarNorm { istep = 10 } v-LEVEL { v = V }

• Rossby--Haurwitz wave with a free surface

    OutputScalarNorm { istep = 10 } v-LEVEL { v = V }

• Rossby--Haurwitz wave with Saint-Venant

    OutputScalarNorm { istep = 10 } v-LEVEL { v = V }

• Dielectric-dieletric planar balance

    OutputScalarNorm { start = end } v { v = V }

• Balance with solid boundaries

    OutputScalarNorm { start = end } v { v = V }