5. Analysis of the Sectional Forces and the Flow Structure
- A.
Configuration with Fineness Ratio L/D = 7.5
Figure 7 plots with orange lines instantaneous sectional force coefficients for 20 instants chosen within the last 0.5 s of the calculation, for which the overall forces are represented in
Figure 6a. The thick blue line represents the discrete average curve, and the solid black circles are the experimental sectional force coefficients for the reference configuration of the fineness ratio of L/D = 15. From the visibility of the orange lines, we can conclude that Region 3 extends approximately up to x/D ~ 5. The rest of the body corresponds to Region 2, being responsible for the small oscillations that appear in
Figure 6a. The averaged side force reaches maxima or minima at x/D ≈ 0, 3.3 and 6.7.
Figure 8 shows iso-Q contours, colored with viscosity ratio
, computed at the last moment of the calculated period (t = 5 s).
Q is the second invariant of the velocity gradient tensor:
where
ω is the norm of the vorticity tensor, and
S is the norm of the strain rate tensor [
21]. These iso-Q contours help to identify coherent vortex structures. The first maximum of side force occurs at x/D = 3.3 (see
Figure 7), while the main vortex sheds between sections x/D = 4.5 and x/D = 4.75.
The solution is consistent with the analysis of Degani et al. [
16]: Region 3 is dominant, and Region 1 is virtually nonexistent for bodies with a lower fineness ratio (L/D < 16). There exists an effect of the base flow and an incipient Region 2, but the fluctuations of the flow are very small, leading to a practically stationary flow, as can be checked in the time history of the overall forces (see
Figure 6).
- B.
Configuration with Fineness Ratio L/D = 15
For the reference configuration of the fineness ratio of L/D = 15, an unsteady solution obtained over 2.1 s with the RSM model was used as the initial solution. The calculation was then resumed using the RSM-SAS model, computing four additional seconds. The solutions of the sectional side and normal force coefficients for several points in time over the last 0.1 s of the interval are shown in
Figure 9 (20 solutions taken at each time step
s). In this case, these two regions are clearly differentiated. Up to x/D ~ 7.0, all solutions coalesce (Region 3). In the rear part, the multiple lines indicate the unsteadiness of the flow (Region 2). Iso-surfaces of positive Q-function (up to 5000) are shown in
Figure 10 for the time t = 6.1 s. They are colored with the viscosity ratio magnitude (left) and vorticity (right).
The flow structure of the short configuration shown previously is similar to that of this reference configuration, as can be seen by comparing
Figure 8 and
Figure 10. The main difference lies in the length of Region 2, which is very small for the short configuration. The first maximum of side force occurs at the same location for both configurations, i.e., x/D = 3.3. The main vortex sheds in a similar location as that of the short configuration. In general, the flow structures of both configurations are similar up to x/D ≈ 7 (see
Figure 7 and
Figure 9), showing that the nose effect is dominant.
The red line marked in the body of
Figure 10 (left) shows approximately the transition between Regions 2 and 3. This corresponds to x/D ~ 7. There is a continuous vortex structure in the steady region, while the flow breaks into different scales in Region 2.
The blue marks with symbols ‘1, 2, and 3’ in the right figure (
Figure 10) indicate the approximate location of vortex shedding at the stationary region (Region 3), which does not correspond exactly with the maximum or minimum of the local side force coefficient. The minima or maxima (averaged side force coefficient) are located at x/D = 0.0, 3.33, 6.67, 9.24, 11.97 and 14. 1 (see
Figure 9). The vortex labeled 1 (
Figure 10) sheds at the nose approximately. However, the vortex labeled ‘2’ is shed at x/D = 4.50 approximately. This can be explained by looking at
Figure 11, which shows the pseudo-streamlines at four longitudinal sections and the Q-function contours taken at one of the 20 solutions analyzed. At this region, the flow is steady (see
Figure 9). At section x/D = 3.30 (local minimum of side force coefficient), the vortex labeled ‘2’ is attached, and its center is located at leeside and close to the y = 0 plane. There is an incipient vortex 3 attached at the starboard side but of lower strength. This balance leads to a maximum side force magnitude. This pair of vortices develop such that vortex 3 increases its strength while vortex 2 reduces it, thus contributing to reducing the side force; at section x/D = 4.30, vortex 3 is attached, and vortex 2 core is located at the plane of symmetry. In section x/D = 4.40, there is flow entrainment in vortex 2 from vortex 3. At section x/D = 4.50, vortex 2 is detached. The shedding occurs closer.
A similar pattern is reproduced with vortex 3 and a new one labeled ‘4’, leading to a vortex shedding at x/D = 6.67, which coincides approximately with another peak in the local side force coefficient. After this, from x/D = 7.0 on, there is an unsteady flow, and the flow breaks into several turbulent scales.
The Power Spectral Density (PSD) of the side and normal forces for the interval t = 2.1–6.1 s are shown in
Figure 12. There is a small content of energy above 20 Hz. The characteristic crossflow time and Strouhal number of the von Karman street are:
.
The theoretical value St = 0.2, characteristic of the infinite cylinder in crossflow, yields a frequency of 9.6 Hz, not concordant with the dominant frequency of 7.3 Hz in this calculation, which corresponds to a Strouhal of 0.15. It is worth noting that the experimental Strouhal for this configuration was 0.16, in accordance with the results by Prananta et al. [
14]. For the short configuration, there was no power spectral density analysis, as the flow was basically steady on the entire body.
Some other frequencies also display a significant content of energy. It is important to remark that for the correct prediction of the low-frequency part of the spectra of the forces it is needed a long physical time integration [
22]. For this reference configuration, the time integration was t = 4 s (8000 Δ
t). The dominant frequency is in accordance with the experimental Strouhal referred to in reference [
14], but this lower frequency of 4.8 Hz may be somewhat amplified for numerical reasons. Regarding the normal force coefficient, the amplitudes of the normal force oscillations are an order of magnitude lower. There are two dominant frequencies (11.5 and 16 Hz). It is important to remark that the normal force is independent of whether the vortex asymmetry is left or right-handed [
9]. The double frequency of the dominant frequency for the side force is 14.6 Hz, which lies between these two frequencies. The lower frequency with energy content is 3.1 Hz, similar to that of 2.8 Hz for the side force coefficient.
A deeper insight making a longer physical time integration to compare to the actual results is necessary to analyze the origin of this content of energy at those lower frequencies correctly.
In the test model, there were pressure data taken in several sections. Based on the crossflow dynamic pressure, we define the following pressure coefficient: .
Figure 13 plots the solutions of pressure coefficients in a circumferential direction for 20 instants over the last 0.1 s (magenta lines), together with the average values (blue lines) and the available experimental data (black dots). The figure includes four cross-sections. There is no available experimental information for the last section. It is important to remark that for this comparison, the origin (Φ = 0) of the roll angle corresponds to the leeside. Positive values denote locations at the starboard side, while negative Φ values represent the portside. The first and second sections (x/D = 1.67 and 3.67) are located in Region 3 of the flow. The third section (x/D = 6.29) is located just in the transition between Regions 2 and 3; it corresponds approximately to the maximum side force location: x/D = 6.7 (see
Figure 9). The flow in the leeside is slightly fluctuating. These fluctuations significantly increase at x/D = 9 (Region 2). Note the oscillation of the minimum pressure values. The unbalance of all the curves on the same side indicates that there is no alternate shedding.
The unsteadiness of the flow is located at the ogive and fore cylinder (up to x/D = 7.0), as shown by the force coefficients curves (
Figure 9), iso-Q contours (
Figure 10) and pressure coefficients (
Figure 13). Then, the oscillations occur at this part, although the different body parts have not been analyzed separately. In the very good detailed work shown in reference [
20] it is shown the time histories for the side force coefficients of the ogive and cylinder of a body of fineness ratio of L/D = 20. It is clear that the flow is steady in the nose region, while there is a fluctuating flow with a Strouhal number of St = 0.15 approximately. Although this calculation is done at laminar flow conditions (low Reynolds), this Strouhal is close to the value obtained for this configuration of lower fineness ratio (L/D = 15) at the same Mach number but at turbulent flow conditions.
The good comparison with experiments both in local sectional forces and pressures (see
Figure 9 and
Figure 13) indicates that the flow field structure shown in
Figure 10 with two differentiating regions—one steady in the forebody and one unsteady in the rear cylinder part—approximates to the real flow field close to the body at high incidence. The contours of the Q-function, together with the pseudo-streamlines of
Figure 11, help to identify the vortex structure and the alternating vortex shedding. A detailed explanation of this vortex structure in the forebody region is given in [
13].
This type of diagram, together with the averaged pressure coefficients at different sectional cross-sections, permits the study of the vortex structure.
Figure 14 plots the solutions of the averaged pressure coefficients in a circumferential direction at several cross-sections over the last 0.1 s.
The pressure coefficients at x/D = 2.8 to 4.0 are shown in the top left part of the Figure. They indicate a similar flow structure in this region, covering the ogive part of the body. The vortex core is located in the plane of symmetry at leeside, yielding a peak of suction at roll angle Φ = 0 deg. (measured from the top) whereas the port side separation line moves towards Φ = −90 deg. The presence of a secondary vortex is visible in the pressure coefficient at the port side. The minimum side force coefficient location is close to x/D = 3.2. It is worth noting that the main contribution to the side force comes from the attached flow region [
23]. In this section, the difference in minimum pressure coefficient between the port and starboard sides reaches its maximum. There is a dominant vortex 2 and an incipient counter-rotating vortex 3. The top right figure shows the pressure coefficients at cross sections x/D = 4.0 to 5.6. As shown in
Figure 11, at section x/D = 4.20 (between 4.00 and 4.4 of the figure), there is flow entrainment from vortex 3 to vortex 2. It can be seen that at region x/D = 4.4, the suction peak is reduced, leading to an increment of the minimum pressure at the leeside (Φ = 0 deg.), and the vortex 2 is detached at section x/D = 4.6. Between sections x/D = 5.2 and 5.6, vortex 2 is shed while vortex 3 evolves increasing its strength and moving its center towards the plane of symmetry at leeside (Φ = 0 deg., plane y = 0). The pressure coefficient at section x/D = 5.6 shows, as a result, a lateral force that decreases as the differences between minimum pressures at the port and starboard sides (close to Φ = −90 deg. and Φ = 90 deg.) are reduced significantly compared to that of section x/D = 3.2. The pressure coefficients shown in the bottom left figure (
Figure 14) are those corresponding to the cross sections x/D = 5.6 to 7.2. The change to positive side forces and to the transition to the unsteady flow region is shown in these plots. There is a dominant vortex 3 at the port side and an increasing strength vortex 4 at the starboard side, leading to a maximum side force between x/D = 6.6 and 7.2 (close to x/D = 7.0). In this section, the maximum difference in minimum pressure coefficients close to the separation lines (around Φ = −90 and Φ = 90 deg.) is produced.
The pressure coefficients shown in the bottom right of
Figure 14 are those corresponding to the cross sections x/D = 7.2 to 14. They show a more complex pattern, leading to a single plateau in the average pressure coefficient at leeside. This is similar to the pressure coefficient of the infinite circular cylinder at crossflow. The flow structure becomes very complex and unsteady. The Strouhal number for this configuration at incidence is St = 0.15. The contribution comes from this rear part of the cylinder.
Finally, we have analyzed the separation and attachment lines for this configuration and reference case [
23]. The information for 20 instants over the last 0.1 s is utilized. In a body at incidence, two main zones can be distinguished at a cross-section. On the windward, the boundary layer is attached. On the leeside, primary and secondary separations give rise to a complex flow field. This can be seen in
Figure 15 (left) at cross-section x/D = 3.67. The pressure contours and the pseudo streamlines are plotted. The attachment point is defined as the minimum velocity point, while the separation points are those where the boundary layer separates. In the figure, the angles Φ
1 and Φ
2 define the location in a circumferential direction of these points, one located at the port side and another at the starboard side. It has been checked in
Figure 14 that the attachment points are located close to Φ = 90 or Φ = −90 deg., according to the definition used for the pressure coefficient (Φ = 0 located at the leeside in the plane y = 0). It is interesting to calculate the attachment line location for the whole body and also the separations lines. The separation lines are calculated with the criterion of minimum shear stress magnitude, being this magnitude
. A plot of the separation lines (red color and left scale) together with the attachment line (blue color and right scale) is shown on the right side of
Figure 15. It is important to remark that for this figure, the origin of the angle is at the windward side (as seen in
Figure 15). An additional green line curve represents
(units in degrees). This curve represents the difference between the separation points on the port and starboard sides. This is an indicator of the level of asymmetry in each section. This way,
indicates symmetric separation points. The main conclusion obtained by looking at the figure is that the curves follow the same trend as the sectional side force coefficient curve (see
Figure 9). However, the attachment line evolves in an opposite phase. A clockwise turn of the separation points is accompanied by a counterclockwise movement of the attachment point. The flow passing through the starboard side feeds the strong central vortex: a low-pressure region compared to the port side values. As a consequence, a greater volume of the incoming flow is forced to pass through starboards, shifting the attachment point to the starboard side. In the unsteady region (X/D > 7.00), both the separation and the attachment lines tend to be symmetric scenarios in accordance with the decay of the cross-sectional side force amplitude. It is also worth noting that the attachment line is offset at x = 0 (about 1.5 deg. to the port side), which is an indicator that the flow is asymmetric from the nose (steady region) and then oscillates around this biased value. It is worth noting that at x/D = 7.0 approximately, both the attachment and separation lines are thicker lines, showing the fluctuation of the attachment and separation points at the cross sections of this rear zone. The unsteadiness of the flow is reflected by the fluctuation of these lines, which is greater as one moves toward the base region. This is another way of measuring and defining the unsteady flow Region 2.
- C.
Configuration with fineness ratio L/D = 30
The time history of the overall force coefficients for the long configuration can be seen in
Figure 6c. The period is 9.4 s. The averaged values and standard deviations are printed in
Table 1.
Figure 16 plots the sectional side and normal force coefficients for 20 instants, taken each
s, over a period of 0.1 s of calculation. As before, the blue line represents the discrete average values and the solid black circles, the experimental sectional force coefficients for the reference configuration. The numerical values in the ogive are similar to the experimental data of the reference configuration, indicating a small influence of the base flow in the nose. The steady zone stretches up to x/D ~ 11. Recall that this value was x/D ~ 6.5–7 for the reference configuration.
The averaged side force reaches maxima or minima at x/D ≈ 0, 3.3, 6.7, 9.4, 12.1, 14.6, 17 and 19.4 (see
Figure 16). After x/D = 20, it remains near zero, while the instantaneous force displays alternating positive and negative signs, the expected behavior of the vortex Karman street (Region 1). Thus, the three regions identified by Ramberg [
15] and Degani et al. [
16,
17] are visible in
Figure 16: Region 3 extends from the nose up to x/D = 11, approximately; Region 2, from x/D = 10 to 20, and Region 1, from x/D = 20 to the base.
Iso-surfaces of positive Q-function (up to 5000) are shown in
Figure 17 for the time t = 9 s. They are colored with the viscosity ratio magnitude.
The red lines marked in the body of
Figure 17 show approximately the transition between Regions 1 and 2 and Regions 2 and 3. This corresponds to x/D ~ 20 and x/D ~ 10 respectively. There is a continuous vortex structure in the steady region, while the flow breaks into different scales in Regions 1 and 2.
We can also appreciate the effect of the von Karman vortex street by looking at the pressure distributions along crossflow sections.
Figure 18 plots the pressure coefficient at a rear section at two points in time approximately half a cycle apart. The pressure at t = 2.55 s changes its pattern with respect to the pressure at t = 2.50 s. During this interval, the minimum pressure induced by the dominant vortex shifts from the port side to the starboard side.
Figure 19 shows the normalized Q-criterion and pseudo-streamlines at section x/D = 23.0 (rear body, flow Region 1), in which the flow has a similar behavior to that of section x/D = 25.0. These plots correspond to three different times (t = 2.505, 2555, and 2.605 s) within a period of 0.1 s again. There is a complex vortex structure with vortices of different scales. The nature of the fluctuating flow leading to a zero net side force within the period is shown here. However, the differences between the contours at t = 2.60 s (right side of the figure) compared to t = 2.505 s (left side of the figure) indicate the complexity of the flow at the rear part of the body.
Figure 20 shows 20 instantaneous pressure coefficients in circumferential direction during a period of 0.1 s and their average at several x/D cross sections. The first two sections (x/D = 3.00 and x/D = 9.00) belong to Region 3 (the steady flow region), the next two sections (x/D = 12.0 and 15.00) belong to the unsteady Region 2 (with oblique vortices), and the last two sections (x/D = 20.0 and 25.00) belong to the unsteady flow Region 1, in which the alternating vortex Karman street pattern is appreciated. The legend
is referred to the roll angle due to calculations at different roll angles being conducted. In this case, it is not relevant.
Pressures at x/D = 3 and x/D = 9 present the characteristic steady asymmetric behavior of Region 3. Note that these stations are close to the second and fourth peaks of the sectional side force depicted in
Figure 16a. Due to the spatial, though steady, alternation of vortices, the stronger depression occurs for negative values of Φ (port side). The pressures would be mirrored at stations near the third peak (x/D = 6.7). Pressures in Region 2 (x/D = 12 and x/D = 15) maintain the asymmetric behavior though introducing a level of unsteadiness. Note that the first station is close to a side-force peak (x/D = 12.1), whereas the last one is located between the peaks at x/D = 14.6 and x/D = 17, resulting in a lower level of asymmetry. Pressures at x/D = 20 and x/D = 25 exhibits a plateau in the leeside, whereas the instantaneous pressure coefficients display a change in the minimum value from one side to the other at a period of 0.1 s, indicating the von Karman street type of flow in this rear part of the body (Region 1).
PSD provides an alternative means to visualize the effect of Karman street.
Figure 21 shows the PSD of the side and normal force coefficients for the long configuration obtained over a period of 9 s.
There are two low frequencies for the side force with significant peaks, 9.7 Hz and 11.6 Hz (St = 0.194 and St = 0.23), the first of which is very close to the value of 10 Hz obtained with a Strouhal number St = 0.2. The normal force, as in the reference case, presents amplitudes an order of magnitude lower. The corresponding figure displays some content of energy at low frequencies (2.8 Hz), a peak at 20.5 Hz, which relates to the Karman shedding (the normal force is independent of whether the vortex asymmetry is left or right-handed), and a high-frequency peak at 87.5 Hz.
It is important to remark again that for the correct prediction of the low-frequency part of the spectra of the forces it is needed a long physical time integration [
22]. For this reference configuration, the time integration was t = 9.4 s (18800 Δ
t). This is a larger (more than double) period than that used for the calculation of the reference configuration (t = 4 s). However, the level of unsteadiness is larger. Therefore, the peaks at low frequencies observed in both the power spectral densities of the normal force coefficients of both configurations may be distorted by numerical noise. These low frequencies have little energy content for the side force coefficient of the long configuration (
Figure 21a), while there is some content for the reference configuration (
Figure 12a).
In order to obtain additional information, time histories for the pressure at four roll angles (Φ = 0, 180, 90, −90 deg.) and four cross sections (x/D = 6, 12, 18, 24) are shown in
Figure 22. The roll angle Φ = 0 corresponds to the windward side. It is interesting to notice that the high-frequency peak of the normal force at 87 Hz may be due to an oscillation of the stagnation point at the rear part of the body (Region 1). This oscillation appears in the last second of the evaluated period. An oscillation at 33 Hz appears in the stagnation point at all sections, affecting both the side and normal force coefficients. The oscillations of separation lines correspond to the pressure oscillations shown in
Figure 22 at the locations
and occur in the Karman street region. They exhibit the dominant frequency of 9.7 Hz in Region 1.
A plot of the separation lines (red color and left scale) together with the attachment line (blue color and right scale) is shown in
Figure 23. There are some spurious values for the separation line plots due to the type of calculation of the minimum shear stress. The main conclusion obtained by looking at the figure is that the curves follow the same trend as the sectional side force coefficient curve (see
Figure 16). However, there is apparently an important difference compared to the solution of the reference configuration (see
Figure 15). The attachment line for the reference configuration was a thin line for sections x/D > 7.00, reinforcing the steady nature of the flow there, as shown by the sectional force coefficients, pressure coefficients and Q-contours curves (
Figure 9,
Figure 10 and
Figure 13). In this case, the correspondent curves for the long configuration (
Figure 16,
Figure 17 and
Figure 20) indicate that a steady flow Region 3 exists from the nose up to x/D = 10.0 approximately. The separation lines curves indicate this up to x/D = 10.0. Nevertheless, the attachment line does not show the same behavior. There is a small fluctuation of the attachment line in this region which does not significantly affect the resultant pressures and then the resultant sectional and global forces. The explanation for this apparent fluctuation of the stagnation line follows.
First of all, the method of calculating this line is one seeking a minimum in a region that covers up to eight degrees (measured from
) where the values are very similar. The shear stress magnitude curve is flat. Small deviations can define a different position. It must be reminded that the cell size is 1.5 deg. and the attachment line can move from one cell to its neighbor cell. Then, a 1.5 deg. fluctuation observed in the nose region for the attachment line can be explained due to this method of calculation. Secondly, the pressure at the region close to
indicates a small oscillation with a frequency of 32–33 Hz on the entire body (
Figure 22), including the nose region. This is reflected in that the power spectral density of the side and normal force coefficients (
Figure 21) show some energy content at this frequency. Then, a small fluctuation of the stagnation line may occur also in the named steady flow Region 3 which covers the fore body (up to x/D = 10 approximately).
The pressure oscillations at
and
regions indicated in
Figure 22 occur at the rear body in Region 1 (From x/D = 20.0 to 30.0). The frequency corresponds to a St = 0.2, showing the Karman street-type region. The separation lines oscillate consequently, as shown in
Figure 23.
We have indicated that it was detected an oscillation of the stagnation point at the rear part of the body (Region 1). This oscillation appears in the last second of the evaluated period, and it is a high-frequency oscillation (f = 87.5 Hz, see
Figure 21 and
Figure 22). At the rear part, there are two main frequencies for the stagnation line: 32 Hz and 87.5 Hz, while in Regions 2 and 3, only 32 Hz with lower amplitudes.