Английская Википедия:Boundary layer thickness

Материал из Онлайн справочника
Перейти к навигацииПерейти к поиску

Шаблон:Short description This page describes some of the parameters used to characterize the thickness and shape of boundary layers formed by fluid flowing along a solid surface. The defining characteristic of boundary layer flow is that at the solid walls, the fluid's velocity is reduced to zero. The boundary layer refers to the thin transition layer between the wall and the bulk fluid flow. The boundary layer concept was originally developed by Ludwig Prandtl[1] and is broadly classified into two types, bounded and unbounded.[2] The differentiating property between bounded and unbounded boundary layers is whether the boundary layer is being substantially influenced by more than one wall. Each of the main types has a laminar, transitional, and turbulent sub-type. The two types of boundary layers use similar methods to describe the thickness and shape of the transition region with a couple of exceptions detailed in the Unbounded Boundary Layer Section. The characterizations detailed below consider steady flow but is easily extended to unsteady flow.

The bounded boundary layer description

Bounded boundary layers is a name used to designate fluid flow along an interior wall such that the other interior walls induce a pressure effect on the fluid flow along the wall under consideration. The defining characteristic of this type of boundary layer is that the velocity profile normal to the wall often smoothly asymptotes to a constant velocity value denoted as ue(x). The bounded boundary layer concept is depicted for steady flow entering the lower half of a thin flat plate 2-D channel of height H in Figure 1 (the flow and the plate extends in the positive/negative direction perpendicular to the x-y-plane). Examples of this type of boundary layer flow occur for fluid flow through most pipes, channels, and wind tunnels. The 2-D channel depicted in Figure 1 is stationary with fluid flowing along the interior wall with time-averaged velocity u(x,y) where x is the flow direction and y is the normal to the wall. The H/2 dashed line is added to acknowledge that this is an interior pipe or channel flow situation and that there is a top wall located above the pictured lower wall. Figure 1 depicts flow behavior for H values that are larger than the maximum boundary layer thickness but less than thickness at which the flow starts to behave as an exterior flow. If the wall-to-wall distance, H, is less than the viscous boundary layer thickness then the velocity profile, defined as u(x,y) at x for all y, takes on a parabolic profile in the y-direction and the boundary layer thickness is just H/2.

At the solid walls of the plate the fluid has zero velocity (no-slip boundary condition), but as you move away from the wall, the velocity of the flow increases without peaking, and then approaches a constant mean velocity ue(x). This asymptotic velocity may or may not change along the wall depending on the wall geometry. The point where the velocity profile essentially reaches the asymptotic velocity is the boundary layer thickness. The boundary layer thickness is depicted as the curved dashed line originating at the channel entrance in Figure 1. It is impossible to define an exact location at which the velocity profile reaches the asymptotic velocity. As a result, a number of boundary layer thickness parameters, generally denoted as <math>\delta(x)</math>, are used to describe characteristic thickness scales in the boundary layer region. Also of interest is the velocity profile shape which is useful in differentiating laminar from turbulent boundary layer flows. The profile shape refers to the y-behavior of the velocity profile as it transitions to ue(x).

Файл:Bounded Boundary Layer Thickness Figure 1.jpg
Figure 1: Schematic drawing depicting fluid flow entering the bottom half of a 2-D channel with plate-to-plate spacing of H. The flow and channel extend perpendicular to the x-y-plane.

The 99% boundary layer thickness

The boundary layer thickness, <math>\delta</math>, is the distance normal to the wall to a point where the flow velocity has essentially reached the 'asymptotic' velocity, <math>u_e</math>. Prior to the development of the Moment Method, the lack of an obvious method of defining the boundary layer thickness led much of the flow community in the later half of the 1900s to adopt the location <math>y_{99}</math>, denoted as <math>\delta _{99}</math> and given by

<math>u(x,y_{99}) = 0.99u_e(x) \quad ,</math>

as the boundary layer thickness.

For laminar boundary layer flows along a flat plate channel that behave according to the Blasius solution conditions, the <math>\delta _{99}</math> value is closely approximated by[3]

<math> \delta _{99}(x)\approx 5.0 \sqrt{ {\nu x}\over u_0} = 5.0 {x \over \sqrt{\mathrm{Re}_x}} \quad ,</math>

where <math>u_e \approx u_0</math> is constant, and where

<math>\mathrm{Re}_x</math> is the Reynolds number,
<math>u_0</math> is the freestream velocity,
<math>u_e</math> is the asymptotic velocity,
<math>x</math> is the distance downstream from the start of the boundary layer, and
<math>\nu </math> is the kinematic viscosity.

For turbulent boundary layers along a flat plate channel, the boundary layer thickness, <math>\delta </math>, is given by[4]

<math> \delta (x) \approx 0.37 {x \over {\mathrm{Re}_x}^{1/5}} \quad .</math>

This turbulent boundary layer thickness formula assumes 1) the flow is turbulent right from the start of the boundary layer and 2) the turbulent boundary layer behaves in a geometrically similar manner[5] (i.e. the velocity profiles are geometrically similar along with the flow in the x-direction, differing only by scaling parameters in <math>y</math> and <math>u(x,y)</math>). Neither one of these assumptions is true for the general turbulent boundary layer case so care must be exercised in applying this formula.

Displacement thickness

The displacement thickness, <math>\delta_1</math> or <math>\delta^*</math>, is the normal distance to a reference plane representing the lower edge of a hypothetical inviscid fluid of uniform velocity <math>u_e</math> that has the same flow rate as occurs in the real fluid with the boundary layer.[6]

The displacement thickness essentially modifies the shape of a body immersed in a fluid to allow, in principle, an inviscid solution if the displacement thicknesses were known a priori.

The definition of the displacement thickness for compressible flow, based on mass flow rate, is

<math> {\delta_1(x)}= \int_0^{H/2} {\left(1-{\rho(x,y) u(x,y)\over \rho_e u_e(x)}\right) \,\mathrm{d}y} \quad ,</math>

where <math>\rho (x,y)</math> is the density. For incompressible flow, the density is constant so the definition based on volumetric flow rate becomes

<math> {\delta_1(x)}= \int_0^{H/2} {\left(1-{u(x,y)\over u_e(x)}\right) \,\mathrm{d}y} \quad .</math>

For turbulent boundary layer calculations, the time-averaged density and velocity are used.

For laminar boundary layer flows along a flat plate that behave according to the Blasius solution conditions, the displacement thickness is[7]

<math> \delta_1(x) \approx 1.72 \sqrt{{\nu x}\over u_0} \quad ,</math>

where <math>u_e \approx u_0</math> is constant.

The displacement thickness is not directly related to the boundary layer thickness but is given approximately as <math>\delta_1 \approx \delta /3</math>.[8] It has a prominent role in calculating the Shape Factor. It also shows up in various formulas in the Moment Method.

Momentum thickness

The momentum thickness, <math>\theta</math> or <math>\delta_2</math>, is the normal distance to a reference plane representing the lower edge of a hypothetical inviscid fluid of uniform velocity <math>u_e</math> that has the same momentum flow rate as occurs in the real fluid with the boundary layer.[9]

The momentum thickness definition for compressible flow based on the mass flow rate is[10][11][12]

<math> \delta_2(x) = \int_0^{H/2} {{\rho(x,y) u(x,y)\over \rho_e u_e(x)} {\left(1 - {\rho(x,y) u(x,y)\over \rho_e u_e(x)} \right)}} \,\mathrm{d}y \quad .</math>

For incompressible flow, the density is constant so that the definition based on volumetric flow rate becomes

<math> \delta_2(x) = \int_0^{H/2} {{u(x,y)\over u_e(x)} {\left(1 - {u(x,y)\over u_e(x)}\right)}} \,\mathrm{d}y \quad ,</math>

where <math>\rho</math> are the density and <math>u_e</math> is the 'asymptotic' velocity.

For turbulent boundary layer calculations, the time averaged density and velocity are used.

For laminar boundary layer flows along a flat plate that behave according to the Blasius solution conditions, the momentum thickness is[13]

<math> \delta_2(x) \approx 0.664 \sqrt{{\nu x}\over u_0} \quad ,</math>

where <math>u_e \approx u_0</math> is constant.

The momentum thickness is not directly related to the boundary layer thickness but is given approximately as <math>\delta_2 \approx \delta /6</math>.[14] It has a prominent role in calculating the Shape Factor.

A related parameter called the Energy Thickness[15] is sometimes mentioned in reference to turbulent energy distribution but is rarely used.

Shape factor

A shape factor is used in boundary layer flow to help to differentiate laminar and turbulent flow. It also shows up in various approximate treatments of the boundary layer including the Thwaites method for laminar flows. The formal definition is given by

<math>H_{12}(x) = \frac {\delta_1(x)}{\delta_2(x)} \quad ,</math>

where <math>H_{12}</math> is the shape factor, <math>\delta_1</math> is the displacement thickness and <math>\delta_2</math> is the momentum thickness.

Conventionally, <math>H_{12}</math> = 2.59 (Blasius boundary layer) is typical of laminar flows, while <math>H_{12}</math> = 1.3 - 1.4 is typical of turbulent flows near the laminar-turbulent transition.[16] For turbulent flows near separation, <math>H_{12} \approx </math>2.7.[17] The dividing line defining laminar-transitional and transitional-turbulent <math>H_{12}</math> values is dependent on a number of factors so it is not always a definitive parameter for differentiating laminar, transitional, or turbulent boundary layers.

Moment method

A relatively new method[18][19] for describing the thickness and shape of the boundary layer uses the mathematical moment methodology which is commonly used to characterize statistical probability functions. The boundary layer moment method was developed from the observation that the plot of the second derivative of the Blasius boundary layer for laminar flow over a plate looks very much like a Gaussian distribution curve. The implication of the second derivative Gaussian-like shape is that the velocity profile shape for laminar flow is closely approximated as a twice integrated Gaussian function.[20]

The moment method is based on simple integrals of the velocity profile that use the entire profile, not just a few tail region data points as does <math>\delta_{99}</math>. The moment method introduces four new parameters that help describe the thickness and shape of the boundary layer. These four parameters are the mean location, the boundary layer width, the velocity profile skewness, and the velocity profile excess. The skewness and excess are true shape parameters as opposed to the simple ratio parameters like the H12. Applying the moment method to the first and second derivatives of the velocity profile generates additional parameters that, for example, determine the location, shape, and thickness of the viscous forces in a turbulent boundary layer. A unique property of the moment method parameters is that it is possible to prove that many of these velocity thickness parameters are also similarity scaling parameters. That is, if similarity is present in a set of velocity profiles, then these thickness parameters must also be similarity length scaling parameters.[21]

It is straightforward to cast the properly scaled velocity profile and its first two derivatives into suitable integral kernels.

The central moments based on the scaled velocity profiles are defined as

<math> {\zeta_n(x)}= \int_0^{H/2} {(y-m(x))^n {1\over\delta_1(x)} \left(1-{u(x,y)\over u_e(x)}\right) \mathrm{d}y} \quad ,</math>

where <math>\delta_1(x)</math> is the displacement thickness and the mean location, <math>m(x)</math> is given by

<math> m(x) = \int_0^{H/2} { y {1\over\delta_1(x)} \left(1-{u(x,y)\over u_e(x)}\right) \mathrm{d}y} \quad .</math>

There are some advantages to also include descriptions of moments of the boundary layer profile derivatives with respect to the height above the wall. Consider the first derivative velocity profile central moments given by

<math> {\kappa_n(x)} = \int_0^{H/2} { (y-{\delta_1(x)})^n {d \{u(x,y)/u_e(x) \} \over dy} \mathrm{d}y} \quad ,</math>

where the first derivative mean location is the displacement thickness <math>\delta_1(x)</math>.

Finally the second derivative velocity profile central moments are given by

<math> {\lambda_n(x)} = \int_0^{H/2} { (y-{\mu_1(x)})^n {d^2 \{-\mu_1(x) u(x,y)/u_e(x) \} \over dy^2} \mathrm{d}y} \quad ,</math>

where the second derivative mean location, <math>\mu_1(x)</math>, is given by

<math> {\mu_1(x)} = {u_{e}(x) \over \left.{\fracШаблон:Du(x,y)Шаблон:Dy} \right|_{y = 0}} = {\upsilon u_e(x) \over \tau_w(x)} \quad ,</math>

where <math>\upsilon</math> is the viscosity and where <math>\tau_w(x)</math> is the wall shear stress. The mean location, <math>\mu_1</math>, for this case is formally defined as ue(x) divided by the area under the second derivative curve.

The above equations work for both laminar and turbulent boundary layers as long as the time-averaged velocity is used for the turbulent case.

With the moments and the mean locations defined, the boundary layer thickness and shape can be described in terms of the boundary layer widths (variance), skewnesses, and excesses (excess kurtosis). Experimentally, it is found that the thickness defined as <math>\delta_m = m + 3\sigma_m</math> where <math>\sigma_m=\zeta_2^{1/2}</math>, tracks the <math>\delta_{99}</math> very well for turbulent boundary layer flows.[22]

Taking a cue from the boundary layer momentum balance equations, the second derivative boundary layer moments, <math>{\lambda_n}</math> track the thickness and shape of that portion of the boundary layer where the viscous forces are significant. Hence the moment method makes it possible to track and quantify the laminar boundary layer and the inner viscous region of turbulent boundary layers using <math>{\lambda_n}</math> moments whereas the boundary layer thickness and shape of the total turbulent boundary layer is tracked using <math>{\zeta_n}</math> and <math>{\kappa_n}</math> moments.

Calculation of the 2nd derivative moments can be problematic since under certain conditions the second derivatives can become positive in the very near-wall region (in general, it is negative). This appears to be the case for interior flow with an adverse pressure gradient (APG). Integrand values do not change sign in standard probability framework so the application of the moment methodology to the second derivative case will result in biased moment measures. A simple fix[23] is to exclude the problematic values and define a new set of moments for a truncated second derivative profile starting at the second derivative minimum. If the width, <math>{\sigma_v}</math>, is calculated using the minimum as the mean location, then the viscous boundary layer thickness, defined as the point where the second derivative profile becomes negligible above the wall, can be properly identified with this modified approach.

For derivative moments whose integrands do not change sign, the moments can be calculated without the need to take derivatives by using integration by parts to reduce the moments to simply integrals based on the displacement thickness kernel given by

<math> {\alpha_n}(x)= \int_0^{H/2} {y^n \left(1-{u(x,y)\over u_e(x)}\right) \mathrm{d}y} \quad .</math>

For example, the second derivative <math>\sigma_v</math> value is <math>\sigma_v = \sqrt{{-\mu_1}^2+2\mu_1 \alpha_0}</math> and the first derivative skewness, <math>\gamma_1</math>, can be calculated as

<math>\gamma_{1}(x) = \kappa_3/\kappa_{2}^{3/2} = (2\delta_1^{3} - 6\delta_1 \alpha_1 + 3\alpha_2)/(2\alpha_1 - \delta_1^{2})^{3/2} \quad .</math>

This parameter was shown to track the boundary layer shape changes that accompany the laminar to turbulent boundary layer transition.[24]

Numerical errors encountered in calculating the moments, especially the higher-order moments, are a serious concern. Small experimental or numerical errors can cause the nominally free stream portion of the integrands to blow up. There are certain numerical calculation recommendations[25] that can be followed to mitigate these errors.

The unbounded boundary layer description

Unbounded boundary layers, as the name implies, are typically exterior boundary layer flows along walls (and some very large gap interior flows in channels and pipes). Although not widely appreciated, the defining characteristic of this type of flow is that the velocity profile goes through a peak near the viscous boundary layer edge and then slowly asymptotes to the free stream velocity u0. An example of this type of boundary layer flow is near-wall air flow over a wing in flight. The unbounded boundary layer concept is depicted for steady laminar flow along a flat plate in Figure 2. The lower dashed curve represents the location of the maximum velocity umax(x) and the upper dashed curve represents the location where u(x,y) essentially becomes u0, i.e. the boundary layer thickness location. For the very thin flat plate case, the peak is small resulting in the flat plate exterior boundary layer closely resembling the interior flow flat channel case. This has led much of the fluid flow literature to incorrectly treat the bounded and unbounded cases as equivalent. The problem with this equivalence thinking is that the maximum peak value can easily exceed 10-15% of u0 for flow along a wing in flight.[26] The differences between the bounded and unbounded boundary layer was explored in a series of Air Force Reports.[27][28][29]

The unbounded boundary layer peak means that some of the velocity profile thickness and shape parameters that are used for interior bounded boundary layer flows need to be revised for this case. Among other differences, the laminar unbounded boundary layer case includes viscous and inertial dominated regions similar to turbulent boundary layer flows.

Файл:Unbounded Boundary Layer Thickness.jpg
Figure 2: The depiction of the laminar “unbounded” boundary layer along a 2-D flat plate with the flow and plate extending perpendicular to the x-y-plane.

Moment method

For exterior unbounded boundary layer flows, it is necessary to modify the moment equations to achieve the desired goal of estimating the various boundary layer thickness locations. The peaking behavior of the velocity profile means the area normalization of the <math>\zeta_n(x)</math> moments becomes problematic. To avoid this problem, it has been suggested[30] that the unbounded boundary layer be divided into viscous and inertial regions and that the boundary layer thickness can then be calculated using separate moment integrals specific to that region. That is, the inner viscous region of laminar and turbulent unbounded boundary layer regions can be tracked using modified <math>{\lambda_n}</math> moments whereas the inertial boundary layer thickness can be tracked using modified <math>{\zeta_n}</math> and <math>{\kappa_n}</math> moments. The slow rate at which the peak asymptotes to the free stream velocity means that the calculated boundary layer thickness values are typically much larger than the bounded boundary layer case.

The modified <math>{\zeta_n}</math> and <math>{\kappa_n}</math> moments for the inertial boundary layer region are created by: 1) replacing the lower integral limit by the location of the velocity peak designated by <math>{\delta_{max}}</math>, 2) changing the upper integral limit to h where h is located deep in the free stream, and 3) changing the velocity scale from <math>u_e</math> to <math>u_0</math>. The displacement thickness in the modified moments must be calculated using the same integral limits as the modified moment integrals. By taking <math>\delta_{max}</math> as the mean location, the modified 3-sigma boundary layer thickness becomes <math>\delta_m = \delta_{max} + 3\sigma_i</math> where <math>\sigma_i</math> is the modified <math>{\zeta_2^{1/2}}</math> width.

The modified <math>{\lambda_n}</math> second derivative moments can be calculated using the same integrals as defined above but with <math>\delta_{max}</math> replacing H/2 for the upper integral limit. To avoid numerical errors, certain calculation recommendations[31] should be followed. The same concerns for the second derivative moments in regards to APG bounded boundary layers for the bounded case above also apply to the modified moments for the unbounded case.

An example of the modified moments are shown for unbounded boundary layer flow along a wing section in Figure 3.[32] This figure was generated from a 2-D simulation[33] for laminar airflow over a NACA_0012 wing section. Included in this figure are the modified 3-sigma <math>\delta _{m}</math>, the modified 3-sigma <math>\delta _{v}</math>, and the <math>\delta _{99}</math> locations. The modified <math>\delta _{m}/\delta _{99}</math> ratio value is 311, the modified <math>\delta _{v}/\delta _{99}</math> ratio value is ~2, and the <math>u_{max}</math> value is 9% higher than the <math>u_{0}</math> value. The large difference between the <math>\delta _{m}</math> and <math>\delta _{v}</math> compared to the <math>\delta _{99}</math> value demonstrates the inadequacy of the <math>\delta _{99}</math> boundary layer thickness. Furthermore, the large velocity peak demonstrates the problem with treating interior bounded boundary layers as equivalent to exterior unbounded boundary layers.

Файл:A Boundary Layer Profile for Air Flow over a NACA0012 Wing.jpg
Figure 3: The velocity profile from a NACA0012 airfoil simulation at x/c = 0.3.[34]

δmax thickness

The location of the velocity peak, denoted as <math>\delta_{max}</math> is an obvious demarcation location for the unbounded boundary layer. The main appeal of this choice is that this location is approximately the dividing location between the viscous and inertial regions. For the laminar flow simulation along a wing,[35] umax located at δmax is found to approximate the viscous boundary layer thickness given as <math>\delta_{max} \approx \delta_v^{4.3} = \mu_1 </math>+<math> 4.3\sigma_v</math> indicating the velocity peaks just above the viscous boundary layer thickness δv. For the inertial regions of both laminar and turbulent flows, <math>\delta_{max}</math> is a convenient lower boundary for the moment integrals. If the width, <math>{\sigma_i}</math>, is calculated using <math>\delta_{max}</math> as the mean location then the boundary layer thickness, defined as the point where the velocity essentially becomes u0 above the wall, can then be properly identified.

The 99% boundary layer thickness

A significant implication of the peaking behavior is that the 99% thickness, <math>\delta _{99}</math>, is NOT recommended[36] as a thickness parameter for the exterior flow, unbounded boundary layer since it no longer corresponds to a boundary layer location of consequence. It is only useful for unbounded laminar flow along a very thin flat plate at zero incidence angle to the flow direction since the peak for this case will be very small and the velocity profile will be closely approximated as the bounded boundary layer case. For thick plates-walls, non-zero incidence angles, or flow around most solid surfaces, the excess flow due to form drag results in a near-wall peak in the velocity profile making <math>\delta _{99}</math> not useful.

Displacement thickness, momentum thickness, and shape factor

The displacement thickness, momentum thickness, and shape factor can, in principle, all be calculated using the same approach described above for the bounded boundary layer case. However, the peaked nature of the unbounded boundary layer means the inertial section of the displacement thickness and momentum thickness will tend to cancel the near wall portion. Hence, the displacement thickness and momentum thickness will behave differently for the bounded and unbounded cases. One option to make the unbounded displacement thickness and momentum thickness approximately behave as the bounded case is to use umax as the scaling parameter and δmax as the upper integral limit.

Notes

Шаблон:Reflist

References

  • Ludwig Prandtl (1904), “Über Flüssigkeitsbewegung bei sehr kleiner Reibung,” Verhandlungen des Dritten Internationalen Mathematiker-Kongresses in Heidelberg 1904, A. Krazer, ed., Teubner, Leipzig, 484–491(1905).
  • Hermann Schlichting (1979), Boundary-Layer Theory, 7th ed., McGraw Hill, New York, U.S.A.
  • Swanson, R. Charles and Langer, Stefan (2016), “Comparison of NACA 0012 Laminar Flow Solutions: Structured and Unstructured Grid Methods,” NASA/TM-2016-219003.
  • Wang, Xia, William K. George and Luciano Castillo (2004), "Separation Criterion for Turbulent Boundary Layers Via Similarity Analysis," J. of Fluids Eng., vol. 126, pp. 297-304.
  • Weyburne, David (2006). "A mathematical description of the fluid boundary layer," Applied Mathematics and Computation, vol. 175, pp.  1675–1684
  • Weyburne, David (2014). "New thickness and shape parameters for the boundary layer velocity profile," Experimental Thermal and Fluid Science, vol. 54, pp. 22–28
  • Weyburne, David (2017), "Inner/Outer Ratio Similarity Scaling for 2-D Wall-bounded Turbulent Flows," arXiv:1705.02875 [physics.flu-dyn].
  • Weyburne, David (2020a). "A Boundary Layer Model for Unbounded Flow Along a Wall," Air Force Tech Report: AFRL-RY-WP-TR-2020-0004,DTIC Accession # AD1091170.
  • Weyburne, David (2020b). "The Unbounded and Bounded Boundary Layer Models for Flow Along a Wall," Air Force Tech Report: AFRL-RY-WP-TR-2020-0005, DTIC Accession # AD1094086.
  • Weyburne, David (2020c). "A New Conceptual Model for Laminar Boundary Layer Flow," Air Force Tech Report: AFRL-RY-WP-TR-2020-0006, DTIC Accession # AD1091187.
  • Whitfield, David (1978). "Integral Solution of Compressible Turbulent Boundary Layers Using Improved Velocity Profiles," AEDO-TR-78-42.

Further reading

  1. L. Prandtl, 1904
  2. Weyburne, 2017
  3. Schlichting, p.140
  4. Schlichting, p. 638
  5. Schlichting, p.152
  6. Schlichting, p. 140
  7. Schlichting, p. 141
  8. Schlichting, p. 28
  9. Schlichting, p. 141
  10. Schlichting, p. 354
  11. Whitfield, p. 13
  12. Schlichting, p. 258
  13. Schlichting, p. 141
  14. Schlichting, p. 161
  15. Schlichting, p. 354
  16. Schlichting, p. 454.
  17. X. Wang, W. George, L. Castillo, 2004
  18. Weyburne, 2006
  19. Weyburne, 2014
  20. Weyburne, 2006, p. 1678
  21. Weyburne, 2017
  22. Weyburne, 2014, p. 26
  23. Weyburne, 2020a
  24. Weyburne, 2014, p. 25
  25. Weyburne, 2014
  26. Weyburne, 2020a
  27. Weyburne, 2020a
  28. Weyburne, 2020b
  29. Weyburne, 2020c
  30. Weyburne, 2020a
  31. Weyburne, 2014
  32. Weyburne, 2020a
  33. R. Swanson and S. Langer, 2016
  34. R. Swanson and S. Langer, 2016
  35. Weyburne, 2020a
  36. Weyburne, 2020a