Processing your ADCP data using structure function techniques: Difference between revisions

From Atomix
Yuengdjern (talk | contribs)
No edit summary
Yuengdjern (talk | contribs)
No edit summary
Line 4: Line 4:


# Extract or compute the [[along-beam bin center separation]] [<math>\delta r_0</math>] based on the instrument geometry
# Extract or compute the [[along-beam bin center separation]] [<math>\delta r_0</math>] based on the instrument geometry
# Calculate the [[along-beam velocity fluctuation]] time-series in each bin <math>n</math>, where [<math>v’(n, t)</math>] from the along-beam velocity data that has met the QC criteria (i.e. the data in Level 2 of the netcdf file)
# Calculate the [[along-beam velocity fluctuation]] time-series in each bin <math>n</math>, where [<math>b’(n, t)</math>] from the along-beam velocity data that has met the QC criteria (i.e. the data in Level 2 of the netcdf file)
# Select the maximum distance (<math>r_{max}</math>) over which to compute the structure function based on conditions of the flow (e.g., expected max overturn). The corresponding number of bins is [<math>n_{\text{rmax}} = r_{max} / r_0</math>]
# Select the maximum distance (<math>r_{max}</math>) over which to compute the structure function based on conditions of the flow (e.g., expected max overturn). The corresponding number of bins is [<math>n_{\text{rmax}} = r_{max} / r_0</math>]
# Calculate the structure function <math>D_{ll}</math> for all possible bin separations <math>\delta</math> within <math>r_{max}</math>  using either a [[bin-centred difference scheme]] or a [[forward-difference]] scheme.  
# Calculate the structure function <math>D_{ll}</math> for all possible bin separations <math>\delta</math> within <math>r_{max}</math>  using either a [[bin-centred difference scheme]] or a [[forward-difference]] scheme.  

Revision as of 09:49, 23 May 2022

To calculate the dissipation rate at a specific range bin and a specific time ensemble:

Schematic showing along-beam distance r and radial velocities.
  1. Extract or compute the along-beam bin center separation [δr0] based on the instrument geometry
  2. Calculate the along-beam velocity fluctuation time-series in each bin n, where [Failed to parse (syntax error): {\displaystyle b’(n, t)} ] from the along-beam velocity data that has met the QC criteria (i.e. the data in Level 2 of the netcdf file)
  3. Select the maximum distance (rmax) over which to compute the structure function based on conditions of the flow (e.g., expected max overturn). The corresponding number of bins is [nrmax=rmax/r0]
  4. Calculate the structure function Dll for all possible bin separations δ within rmax using either a bin-centred difference scheme or a forward-difference scheme.
  5. Perform a regression of Dll(n,δ) against (δr0)2/3 for the appropriate range of bins and δr0 separation distances. Be aware of special considerations for forward-difference, center-difference schemes in setting up the regression calculation. The regression is typically done as a least-squares fit, either as:

    Dll=a0+a1(δr0)2/3; or as
    Dll=a0+a1(δr0)2/3+a3((δr0)2/3)3

    the former being the canonical method that excludes non-turbulent velocity differences between bins, whereas the latter is a modified method that includes non-turbulent velocity differences between bins due to any oscillatory signal (e.g. surface waves, motion of the ADCP on a mooring).
  6. Use the coefficient a1 to calculate ε as

    ε=(a1C2)2/3

    where C2 is an empirical constant, typically taken as 2.0 or 2.1.




Next step: Apply quality-control on dissipation rates (QA2)

Previous step: Apply quality-control on velocity time series data (QA1)

Return to ADCP Flow Chart front page