Monday, December 29, 2014

Radome Induced Vortex Shedding Part 2


Introduction

In a previous ITI-RCS post, we illustrated the nature of vortices that are shed from behind a sphere-cone radome in the presence of a 100 mph wind.  We did this numerically using the Detached Eddy Simulation (DES) scheme present in ANSYS CFX for simulating highly turbulent transient flows.  The simulation was a quick attempt to recreate one of the wind tunnel tests conducted by D’Amato and Fanning in Reference 1.  We used a small fluid domain and coarse mesh (relative to the anticipated turbulence scales) in order to predict the major characteristics of the flow.  This was sufficient to capture the basic response of large eddies shedding and rolling away from behind the radome, however it did not capture the full nature of the flow.

In the current post, we present the results of three different simulation methods for the same wind tunnel test, two of which are of significantly higher fidelity than the previous post.  In addition, we assess the accuracy of the methods for predicting the flow.  The first simulation is the traditional Reynolds Averaged Navier-Stokes (RANS) method.  RANS is a steady state method that averages turbulence in time and models the smallest turbulence scales so that they do not have to be resolved through mesh discretization.  The second method is Scale Adaptive Simulation (SAS).  SAS is a variant of unsteady RANS where the Von Karman length-scale has been included allowing the solver to resolve eddies similar to Large Eddy Simulation (LES).  The last method is DES, which is a combination of RANS and LES.  DES uses RANS in regions where turbulence scales are small (such as near walls) and LES in the rest of the fluid domain.    Both DES and SAS are transient methods.  All of the methods utilize the Shear Stress Transport (SST) turbulence model.


Simulation Meshes

The meshes for the three simulations are shown in Figure 1.  The upper mesh is a half domain model that assumes symmetry exists down the length of the wind tunnel.  A region approximately 15 radome diameters ahead of the radome is included in the fluid domain.  This is done in order to give the boundary layer in the entering fluid time to develop sufficiently.  Approximately 13 diameters are included behind the radome for the development of a wake.  The unstructured mesh contains 1 million cells ranging in size from 0.5 in near the dome to 20 in elsewhere.  An inflation layer is present along the radome wall with a first cell thickness of 5 mils.  A slighter courser inflation layer is present along the wind tunnel walls.
 
We use the mesh shown in Figure 1 (b) & (c) for both the DES and SAS simulations.  The dimensions of the fluid domain are the same as for the RANS model except that no symmetry is assumed.  The mesh is unstructured with 27 million cells ranging in size between 0.25 in near the radome to 20 in far from the radome.  Cells in a region anticipated to be the wake zone are no larger than 0.35 in.  An inflation layer with identical settings as the RANS model is present along the radome wall.  A similar inflation layer is present along the other walls.


(a)
(b)

(c)

Figure 1:  (a) Mesh for RANS model.  (b) Mesh for  DES and SAS models.  (c) Cross sectional view of DES and SAS mesh scoped to the wake zone around the radome.


Vortex Structures

Results of the RANS simulation are shown in Figure 2 and Figure 3.  In both figures the results are mirrored about the symmetry plane to provide a complete picture.   The pressure distribution on the radome is shown along with an iso-surface of the swirl strength around the radome.  This surface is colored based on the flow velocity.  One half is also semi-transparent to reveal the radome surface pressure contours and flow streamlines.  The iso-surface is an indication of the presence of vortices.  We see two major types of vortices in the flow.  The first is horseshoe shaped and arcs around the front and sides of the radome.  Behind the radome are a set of long wake vortices.  These are formed as the flow separates on the backside of the radome.  They are then turned and lengthened as the flow rolls them against the ground.


Figure 2:  RANS simulation results.  Close-up view.
Figure 3:  RANS simulation results.  Broad view.
The steady-state RANS model is only capable of producing a time-averaged vortex response.  One would certainly expect the true transient response to be much more dynamic.

Video 1 shows the transient response that is predicted using the SAS model.  The motion has been slowed by roughly a factor of 30.  As before, vortices around the radome are represented using swirl strength iso-surfaces, colored by the velocity magnitude.  We see that the flow is much richer than predicted using RANS.  The location of flow separation on the radome wanders with time and tightly wound vortices shed rapidly from the radome.  They then roll over one another and coalesce into larger vortices that spread outward and continue to roll downstream. As they move they change shape, spinning off and consuming smaller vortices in a chaotic fashion.  We note that the path of the vortices roughly matches that of the tubular vortices from the RANS model.



Video 2 shows the predicted response of the DES model.  As with the SAS model, time is slowed by a factor of roughly 30.  While the flow characteristics are still rich, the structure of the vortices is less chaotic.  Larger, more stable vortices are shed from the radome surface.  They combine quickly to form vortices that tend to hold their structure better as they grow and roll downstream.




Error Comparison

D’Amato and Fanning measured the time-averaged pressure at over 200 points on the radome surface during their wind tunnel test.  From this they determined the overall lift and drag forces.  Table 1 lists the measured data as well as the predicted values from three simulations.  The table also shows the errors in the maximum and minimum pressures from the simulations, the root mean squares (RMS) of their errors, as well as the lift and drag errors.  We see that the method with the lowest error is DES.  It achieves RMS and lift errors of 5.4% and 1.7%, respectively, which are roughly half that of the RANS method.  Unfortunately, none of the methods predict the drag force well.

The quality of the three methods is easier to assess in Figure 4.  The absolute error in the pressure prediction on the surface of the radome is plotted with identical color scales.  It is clear to see that the DES solution is the most accurate.


Table 1: Radome and simulation data and error.


Figure 4:  Plots of pressure error for RANS (left), SAS (middle) and DES (right) solutions
Reference 1:  D’Amato, R., Fanning, W. R., “Pressure Distributions on Sphere-Cone Radomes in Uniform and Gradient Flows,” Techincal Note, ESD-TR-68-84, Lincoln Library, MIT, 1968.