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 5: | Line 5: | ||
# 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>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) | # 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} / \delta 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. | ||
# Perform a regression of <math>D_{ll}(n,\delta)</math> against <math>(\delta r_0)^{2/3}</math> for the appropriate range of bins and <math>\delta</math>r<sub>0</sub> separation distances. Be aware of [[Regressing structure function against bin separation | 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: <br /><br /> <math>D_{ll} = a_0 + a_1 (\delta r_0)^{2/3}</math>; or as <br /> <math>D_{ll} = a_0 + a_1 (\delta r_0)^{2/3}+a_3((\delta r_0)^{2/3})^3 </math> <br /><br /> the former being the [[canonical structure function method | canonical method]] that excludes non-turbulent velocity differences between bins, whereas the latter is a [[modified structure function method | 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). | # Perform a regression of <math>D_{ll}(n,\delta)</math> against <math>(\delta r_0)^{2/3}</math> for the appropriate range of bins and <math>\delta</math>r<sub>0</sub> separation distances. Be aware of [[Regressing structure function against bin separation | 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: <br /><br /> <math>D_{ll} = a_0 + a_1 (\delta r_0)^{2/3}</math>; or as <br /> <math>D_{ll} = a_0 + a_1 (\delta r_0)^{2/3}+a_3((\delta r_0)^{2/3})^3 </math> <br /><br /> the former being the [[canonical structure function method | canonical method]] that excludes non-turbulent velocity differences between bins, whereas the latter is a [[modified structure function method | 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). | ||
Revision as of 14:37, 27 May 2022
To calculate the dissipation rate at a specific range bin and a specific time ensemble:
- Extract or compute the along-beam bin center separation [[math]\displaystyle{ \delta r_0 }[/math]] based on the instrument geometry
- Calculate the along-beam velocity fluctuation time-series in each bin [math]\displaystyle{ n }[/math], where [[math]\displaystyle{ 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]\displaystyle{ 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]\displaystyle{ n_{\text{rmax}} = r_{max} / \delta r_0 }[/math]]
- Calculate the structure function [math]\displaystyle{ D_{ll} }[/math] for all possible bin separations [math]\displaystyle{ \delta }[/math] within [math]\displaystyle{ r_{max} }[/math] using either a bin-centred difference scheme or a forward-difference scheme.
- Perform a regression of [math]\displaystyle{ D_{ll}(n,\delta) }[/math] against [math]\displaystyle{ (\delta r_0)^{2/3} }[/math] for the appropriate range of bins and [math]\displaystyle{ \delta }[/math]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:
[math]\displaystyle{ D_{ll} = a_0 + a_1 (\delta r_0)^{2/3} }[/math]; or as
[math]\displaystyle{ D_{ll} = a_0 + a_1 (\delta r_0)^{2/3}+a_3((\delta r_0)^{2/3})^3 }[/math]
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). - Use the coefficient [math]\displaystyle{ a_1 }[/math] to calculate [math]\displaystyle{ \varepsilon }[/math] as
[math]\displaystyle{ \varepsilon = \left(\frac{a_1}{C_2}\right)^{2/3} }[/math]
where [math]\displaystyle{ C_2 }[/math] 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
