Bin-centred difference scheme

From Atomix
Revision as of 01:22, 11 November 2021 by Jmmcmillan (talk | contribs) (Created page with "For a '''bin-centred difference''' scheme: # start at bin n = (n<sub>max</sub> / 2) + 1 ## start with <math>\delta</math> = 1 ## if <math>\delta</math> is '''''even''''' compu...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

For a bin-centred difference scheme:

  1. start at bin n = (nmax / 2) + 1
    1. start with δ = 1
    2. if δ is even compute the second order structure function D(n,δ) as the segment mean of the square of the velocity difference between the bins separated by distance δr0 centered around bin n:

      D(n, δ) = [v’(n+(δ/2), t) - v’(n-(δ/2), t)]2

      where the angled brackets indicate the mean across all t for the data segment yielding a velocity difference after the application of the Level 1 QC criteria
    3. if δ is odd compute the second order structure function D(n,δ) as the segment mean of the mean of the square of the velocity difference between the bins separated by distance δr0 centered on the upper and lower extent of bin n:

      dv'lo(n, δ, t) = v’(n+floor(δ/2), t) - v’(n-ceil(δ/2), t)
      dv'hi(n, δ, t) = v’(n+ceil(δ/2), t) - v’(n-floor(δ/2), t)

      where ceil and floor indicate the upper and lower integer value respectively, then

      D(n, δ) = [dv'lo(n, δ, t)2 + dv'hi(n, δ, t)2] / 2

      the angled brackets again indicating the mean across all t in the data segment yielding a velocity difference after the application of the Level 1 QC criteria
    4. increment δ and repeat steps until δ = nmax
  2. increment n and repeat steps until n + (nmax / 2) exceeds the bin number for which valid v’ are available