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

• 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.) }

• 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 }

• 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 }

• 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.) }