Evaluation of Flow Solver Accuracy using Five Simple Unsteady Validation Cases
|
|
- Lawrence Casey
- 6 years ago
- Views:
Transcription
1 Evaluation of Flow Solver Accuracy using Five Simple Unsteady Validation Cases Bradford E. Green, Ryan Czerwiec Naval Air Systems Command, Patuxent River, MD, Chris Cureton, Chad Lillian, Sergey Kernazhitskiy, Tim Eymann, Jason Torres Air Force SEEK Eagle Office, Eglin Air Force Base, FL, Keith Bergeron, Robert Decker United States Air Force Academy, USAF Academy, CO, The goal of this study was to evaluate the time accuracy of five different Navier-Stokes flow solvers using five simple two-dimensional validation cases. The five flow solvers included Cobalt, USM3D, FUN3D, Beggar and NASCART-GT. The five validation cases included the laminar flow over a circular cylinder, the laminar shock tube, the inviscid convecting vortex, the 18%-thick circular-arc airfoil and the NACA 0015 pitching airfoil. This study was conducted in an effort to better understand the capabilities of each flow solver for time-accurate calculations. The eventual goal of the project was to use one or more of these flow solvers to analyze the unsteady transonic flow over a fighter aircraft with stores in captive carriage. Evaluating the capabilities of each flow solver with simple validation cases rather than a complicated fighter aircraft was desirable. a c c a c d c l c m c n Ch e f k L M p Re St T t u x y z = speed of sound = chord of airfoil = 2D axial force coefficient = 2D drag coefficient = 2D lift coefficient = 2D pitching-moment coefficient = 2D normal force coefficient = wall heat coefficient = internal energy = frequency = reduced frequency = reference length = Mach number = static pressure = Reynolds number = Strouhal number = static temperature = time = velocity = x dimension = y dimension = z dimension = angle of attack (degrees) = pitch angle (degrees) = vortex strength = density = thickness Nomenclature 1
2 = vorticity = freestream I. Introduction With recent advances in computers and turbulence models, computational fluid dynamics (CFD) is being used more than ever before to compute the complicated flow over military aircraft. The advances in computers allow for larger grids to be generated and the final CFD solutions to be obtained faster. Steady-state solutions at a single condition on unstructured grids with approximately 50 million cells can be obtained in 12 hours or less. Complete polars with 15 angles of attack can be completed in a few days. Depending on the availability of computers, sometimes these polars can be obtained in only one day. The computer times for structured flow solvers and grids are even less. The above estimates are for steady-state solutions. Sometimes, though, unsteady calculations must be performed to obtain accurate results. These time-accurate calculations add to the complexity of the problem, since the time step and number of Newton subiterations must be chosen properly to generate an accurate solution. Therefore, when doing time-accurate calculations, it is important to understand the impact that time step and Newton subiterations have on the flow solution. Despite the ability to generate larger grids and use improved turbulence models, caution must still be taken when generating flow solutions. Flow solver, grid, turbulence model, time step and Newton subiteration studies should still be conducted. In addition, the results should be compared to known wind-tunnel or flight-test data where possible to ensure that accurate solutions are being generated. The eventual goal of this work is to predict the unsteady flow over a fighter aircraft with stores in captive carriage. Before embarking on this challenging problem, it was desired to validate the CFD tools with several simple two-dimensional validation cases. CFD could only be expected to generate accurate data on the unsteady, three-dimensional fighter-aircraft problem if CFD were first able to accurately predict some basic validation cases. Five unsteady validation cases were selected for this study. These cases include the laminar flow over a circular cylinder, the laminar shock tube, the inviscid convecting vortex, the 18%-thick circular-arc airfoil and the NACA 0015 pitching airfoil. The goal of this work was to analyze these five validation cases with five different flow solvers. These flow solvers include Cobalt, USM3D, FUN3D, Beggar and NASCART-GT. While the developers of all flow solvers use a basic set of validation cases to prove the accuracy of the particular flow solver, these validation cases often vary from one flow solver to another. While these validation cases may be similar between different flow solvers, they are not always identical. In this study, each of the flow solvers will be run on the same set of validation cases allowing for the results to be compared. The flow solvers will be introduced in the following section. The results of the validation cases will then be presented, followed by a brief summary. II. Flow Solvers A. Cobalt The commercial flow solver Cobalt [1] is a Navier-Stokes flow solver that solves the Navier-Stokes equations using a cell-centered, finite-volume approach applicable to arbitrary cell topologies. Since Cobalt accepts arbitrary cell topologies and can be used on both two- and three-dimensional grids, it will be possible to use Cobalt to analyze any of the grids that are generated during this project. Cobalt is second-order accurate in space and time. A Newton sub-iteration scheme is employed for time-accurate calculations. Several different turbulence models have been coded in Cobalt, including Spalart-Allmaras (SA) [2], SA with rotation correction, Wilcox k- and Menter s shear stress transport (SST) [3] model. The SA and SST turbulence models can be combined with the Delayed Detached- Eddy Simulation (DDES) for time-accurate calculations. Cobalt is also capable of analyzing overset grids and cases employing rigid body motion. B. USM3D USM3D [4], which was developed at NASA Langley Research Center and is part of the Tetrahedral Unstructured Software System (TetrUSS), is a three-dimensional, cell-centered, finite-volume Navier-Stokes flow solver for unstructured, tetrahedral meshes. USM3D uses Roe s flux-difference splitting technique to compute the inviscid flux quantities across each cell face. The MinMod flux limiter is also implemented within USM3D. The implicit Gauss-Seidel time-stepping scheme is used to advance the solution to a steady-state condition. USM3D is 2
3 second-order in space and time. The one-equation SA turbulence model and the two-equation Mentor s SST model are available within USM3D. Spalart-Allmaras is used to form a Detached-Eddy Simulation (DES) model for timeaccurate calculations. It is important to mention that the DDES model available in Cobalt is a variation of the DES model available in USM3D. C. FUN3D The FUN3D [5] flow solver was developed at NASA Langley Research Center. FUN3D solves the Navier- Stokes equations on unstructured two- and three-dimensional grids. The flow solver is node-based and uses a finitevolume discretization. FUN3D employs multigrid with implicit time stepping where the linear system is solved using either point-implicit, line-implicit or Newton-Krylov schemes. The FUN3D options for upwind flux functions include Roe flux-difference splitting, Van Leer flux-vector splitting, AUFS and HLLC. The SA and SST turbulence models are also available within FUN3D. The time-accurate calculations within FUN3D are third-order accurate with error controllers. FUN3D is also capable of analyzing problems with time-varying translation and rotation. D. Beggar The Beggar code [6], which originated at the Wright Laboratory, is currently being maintained and developed at the Air Force Seek Eagle Office (AFSEO) at Eglin Air Force Base in Fort Walton Beach, FL. Beggar is a computational tool that can predict the flow over stores as they separate from aircraft and can calculate the trajectories and orientations of the stores as they separate. As a result, grid assembly, flow solver and six degree-offreedom (DOF) components exist within Beggar. For this project, only the capabilities of Beggar as a flow solver will be evaluated. Beggar solves the Euler and Navier-Stokes equations on blocked and overset structured grids. Beggar solves the governing equations using a finite-volume, cell-centered discretization. The Steger-Warming scheme or the Roe scheme can be used to calculate the inviscid fluxes. The Baldwin-Lomax [7], SA and k-epsilon turbulence models are used by Beggar for turbulent flows. Beggar can be used in both steady-state and time-accurate modes. DES is also available for use with SA for time-accurate calculations. E. NASCART-GT NASCART-GT, which was developed at Georgia Tech, is capable of solving the Euler equations, the Euler equations with an integral boundary layer and the Navier-Stokes equations on two- and three-dimensional Cartesian grids. NASCART-GT uses adaptive mesh refinement in high flow gradient regions to improve the accuracy and efficiency of the solution. Before solving the governing equations, NASCART-GT generates a Cartesian grid on the specified geometry in an automated fashion, reducing the grid generation time usually required before a flow solver can be executed. NASCART-GT is up to fifth-order accurate in space and uses Roe s Approximate Riemann solver, AUSM and AUSMPW+ schemes. NASCART-GT is up to second-order accurate in time. The k-epsilon turbulence model is available within NASCART-GT. III. Summary of Validation Cases The goal of this work was to use five validation cases to evaluate the accuracy of each of the five flow solvers. These cases include the laminar flow over a circular cylinder, the laminar shock tube, the inviscid convecting vortex, the NACA 0015 pitching airfoil, and the 18%-thick circular-arc airfoil. All of these cases are viscous, with the exception of the inviscid convecting vortex case. The NACA 0015 pitching airfoil and the 18%-thick circular-arc airfoil are turbulent cases while the circular cylinder and shock tube cases are laminar cases. Table 1 summarizes the flow conditions, reference length and time-steps that were used with Cobalt for each validation case. In this table, the Reynolds number is based on the reference length shown in the table. For some of these cases, the free-stream temperature and pressure are arbitrary, but the Reynolds number shown in the table must be maintained with each flow solver. The flow conditions and reference length for the shock tube and convecting vortex cases are shown in metric units in Table 1. These cases can also be evaluated in English units. In doing so, the grid remains the same and the reference length becomes 1 ft. The pressure and temperature must be set to give the correct value for Reynolds number. The time steps given in Table 1 are non-dimensional. The dimensional time steps can be calculated using the equation: 3
4 t t L a Table 1. Summary of flow conditions for validation cases. Validation Case Mach Re * T + Circular cylinder R psf 1 ft Shock tube , K Pa 1 m P + L t Convecting vortex 0.5 N/A K 101,325 Pa 1 m 0.04 Pitching airfoil x10 6 Prescribed 500 R psf 1 ft Circular-arc airfoil x R psf 1 ft * Reynolds number is based on the reference length L + Pressure and temperature are arbitrary, but the Reynolds number must be maintained In the following sections, each of the validation cases will be discussed. Each of these sections contains a description of the validation case, a description of the computational grids that were used and a summary of results. IV. Laminar Flow over Circular Cylinder A. Description The laminar flow over a circular cylinder was used to evaluate the ability of each flow solver to predict a periodic unsteady viscous flow. For this problem, the laminar Navier-Stokes equations were solved. The flow conditions for this validation case were taken from Nichols et al. [9] and are shown in Table 1. Two different time steps were used for the calculations. The number of Newton sub-iterations was varied to ensure temporal accuracy. The Strouhal number and amplitude of the cylinder lift and drag coefficients were evaluated as well as the average value of the drag coefficient. These same quantities were evaluated by Nichols. The Strouhal number is defined as f St u B. Computational Grids For the circular cylinder validation case, each flow solver analyzed one grid. For Cobalt and BEGGAR, the structured grid of Nichols was used. This grid is shown in Figure 1. This grid had 401 points on the cylinder and 201 points normal to the cylinder. The normal spacing at the surface of the cylinder was diameters. The farfield boundary was placed 98.7 diameters from the surface of the cylinder. For USM3D and FUN3D, an unstructured tetrahedral grid of similar grid spacing was generated. This grid is shown in Figure 2. C. Results The results for the circular cylinder validation case are shown in Figures 3 through 9. In Figure 3, the lift coefficient from the Cobalt flow solver is plotted versus time over one cycle of the flow. On this plot, arrows illustrate the value of the lift coefficient at nine different times along this cycle. The times are labeled with letters A through I. In Figure 4, the Mach contours in the vicinity of the cylinder are shown at the times labeled in Figure 3. 4
5 Using Figure 4, one is able to visualize how the flow is changing with time. The Mach contours illustrate the shedding of vorticity from the cylinder. In Figure 5, the lift amplitude is plotted as a function of the number of sub-iterations for the small and large time steps for Cobalt, USM3D, FUN3D and Beggar. This plot confirms that the lift amplitude converges for both time steps as the number of sub-iterations is increased. The converged values of lift amplitude vary between 0.73 and The lift amplitude is a measure of how damped the solutions are for each of the flow solvers. In Figure 6, the lift Strouhal number is plotted for each of the flow solvers. This figure shows that the lift Strouhal number converges for each of the flow solvers as the number of sub-iterations is increased. Based on the available experimental data, the lift Strouhal number should be approximately to [9, 10 and 11]. The lift Strouhal number for each of the flow solvers was in the vicinity of the experimental data. In Figure 7, the average drag value is plotted for each of the flow solvers. Once again, the average drag value converges as the number of sub-iterations is increased. The values for the converged average drag coefficient predicted by the flow solvers range from 1.29 to Based on experimental data, the average drag coefficient should be approximately 1.34 [9]. In Figure 8, the drag amplitude is plotted for each of the flow solvers. In these plots, the converged drag amplitude varies between and As in the case of the lift amplitude, the drag amplitude is a measure of how damped the solutions are. In Figure 9, the drag Strouhal number is plotted for each of the flow solvers. Once again, the drag Strouhal number converges as the number of sub-iterations is increased. The values of drag Strouhal number in the plot range from 0.34 to approximately Attempts to run the NASCART-GT flow solver on this validation case were not successful. V. Inviscid Convecting Vortex A. Description The inviscid convecting vortex validation case was also used to evaluate each of the flow solvers. For this case, the Euler (inviscid) equations were solved. This validation case presented the opportunity to evaluate the ability of the flow solvers to maintain the shape and strength of a vortex while it is convected downstream. The ability of a flow solver to convect a vortex downstream when analyzing a complex military aircraft is very important. In flight, vortices from the wing on an aircraft can impact the flow over horizontal tails. If the vortex is not properly convected downstream in the CFD simulation, then the resulting forces and moments on the tail will be inaccurate. The boundary conditions for this problem were taken from Nichols [8]. The initial vortex strength was prescribed as: where (2 R)exp(0.5(1 R)) 2 R ( x x z, ) ( z 0) x 0 and z 0 represent the location of the vortex center. A vortex strength of 5 was chosen for this study. Rather than initializing the flow solution with the traditional free-stream values at the beginning of the CFD calculation, each flow solver had to be modified to prescribe the vortex described above. This task alone was very challenging. In several cases, the vortex was not prescribed properly and erroneous results were generated. Caution must be taken when prescribing the initial vortex. In the case of the commercial flow solver Cobalt, the developers of the code made the modifications to the code since the source code is not distributed to the users of the software. The flow conditions used within Cobalt for this problem are shown in Table 1. The vortex was allowed to complete 5 cycles (500 time steps/cycle). The minimum pressure of the vortex was tracked to examine the ability of the flow solvers to maintain the vortex strength. B. Computational Grids For the inviscid convecting vortex validation case, three structured grids were used as the baseline grids. These grids were those of Nichols. These grids were 10 units in height by 10 units in width. The number of points in each of these grids was varied in an effort to determine the impact of grid spacing on the convecting vortex. The coarse 5
6 grid was 41 x 41, which corresponds to a grid spacing of 0.25 units. The medium grid was 81 x 81, which corresponds to a grid spacing of units. The fine grid was 161 x 161, which corresponds to a grid spacing of units. These structured grids are shown in Figure 10. To generate the solutions with Cobalt, the baseline set of structured grids were used. The vortex was placed in the center of the grid. Through the use of a periodic boundary condition, the vortex was progressed five cycles on each grid. Unfortunately, none of the other flow solvers used were capable of using a periodic boundary condition. For the Beggar flow solver, the baseline set of grids (those shown in Figure 10) were modified to eliminate the need for a periodic boundary condition. Each of the grids were modified by extrapolating the grids to the right (in the positive x-direction) while maintaining the grid spacing. The extrapolated grids had a height of 10 units and a width of 50 units. Unstructured tetrahedral grids with the same dimensions and similar grid spacings were generated for use with USM3D and FUN3D. The coarse, medium and fine unstructured grids are shown in Figure 11. C. Results In Figures 12 through 14, the results of the Cobalt, Beggar and FUN3D flow solvers are shown. In these figures, the contours of the vortex after five cycles for the coarse, medium and fine grids are compared to the contours of the initial vortex. The impact of grid spacing on the vortex is obvious. For the fine grid, the contours of the vortex are nearly maintained after five cycles with each flow solver. As the grid becomes more coarse, the contours of the vortex are not maintained. For the ideal flow solver and grid, the contours of the initial vortex would be maintained. In Figure 15, the results of the USM3D flow solver are shown for the fine grid after five cycles of the vortex. The contours of the vortex are maintained on the fine grid. Plots of the medium and coarse grids were not available. In Figure 16, the minimum pressure as a function of the number of vortex cycles and grid size is shown for Cobalt, USM3D, FUN3D and Beggar. The results are consistent across all four flow solvers. Fine grids are better able to convect the vortex than coarse grids. The minimum pressure is plotted in Figure 17 for the NASCART-GT flow solver. In this case, the initial minimum pressure is not maintained, even for the fine grid. A lot of effort was expended to find the source of this problem, which also plagued several of the other flow solvers during the course of this validation case. The problem for all of the flow solvers was traced to an incorrect initialization of the vortex within the flow solver. The other flow solvers were corrected and meaningful results were generated. However, due to time and budget limitations, the NASCART-GT was dropped from the project at this point. VI. Laminar Shock Tube A. Description A laminar shock tube problem was used to evaluate the ability of each flow solver to properly calculate the shock motion as a function of time. For this problem, the laminar Navier-Stokes equations were solved. Including the viscous effects presented the opportunity to calculate the heat transfer characteristics for this problem. The initial conditions for this problem were taken from Nichols [8] and are shown in Table 2. Table 2. Initial conditions for shock tube validation case. / ref p/p ref T wall /T ref Left state Right state The flow conditions for this problem are shown in Table 1. The free-stream values from Table 1 were used for the reference conditions shown in Table 2. Isothermal, no-slip boundary conditions were used on the upper and lower walls and a temperature of degrees Kelvin was specified on the two walls. The density, pressure, velocity and internal energy on the centerline of the shock tube and the wall heat transfer coefficient were calculated at a non-dimensional time of 0.2. The density, pressure, velocity and internal energy were compared to the inviscid theoretical results and the wall heat transfer coefficient was compared with the results of Nichols. B. Computational Grids The structured grid used for the shock tube validation case was that of Nichols. This structured grid, which is shown in Figure 18, was used by Cobalt and Beggar. For FUN3D and USM3D, unstructured tetrahedral grids were 6
7 generated. These grids are shown in Figure 19. Each grid was 1 unit in width and 0.1 units in height. The upper and lower walls were viscous with a normal spacing of approximately units. C. Results The results for the laminar shock tube problem are presented in Figures 20 through 24. All of the results in these figures show the state of the solution at a non-dimensional time of 0.2. Where applicable, the results are compared to the exact solution. The results from all three time steps are shown in the plots. In Figure 20, the density is plotted as a function of x for the Cobalt, USM3D, FUN3D and Beggar flow solvers. The FUN3D results tend to overshoot the exact solution. This is probably a result of the fact that FUN3D is a nodebased solver. A finer grid might yield a better correlation with the exact solution. In Figures 21, 22 and 23, the pressure, velocity and internal energy are plotted as a function of x. The same overshoot with FUN3D is present in each of the figures. In Figure 24, the wall heat transfer coefficient is plotted as a function of x for Cobalt. The other flow solvers were unable to output the wall heat transfer coefficient. The Cobalt results for wall heat transfer coefficient are compared to those of Nichols, which are also included in the plot. The Beggar flow solver was not used for the remaining validation cases in this paper. Time and budget issues dictated that it be dropped from the suite of flow solvers being used during this project. VII. 18%-Thick Circular-Arc Airfoil A. Description The 18%-thick circular-arc airfoil (also known as the bi-convex airfoil) was used to evaluate the flow solvers. For this problem, the turbulent Navier-Stokes equations were solved. The SA and SST turbulence models were evaluated for each flow solver. The CFD results were compared to wind-tunnel data obtained by McDevitt [12]. Shock-induced boundary-layer separation and a hysteresis effect are present in the transonic flow over the bi-convex airfoil. The goal of this validation exercise was to determine if each of the flow solvers could predict these phenomena. The axial force, normal force and pitching-moment coefficients, as well as the reduced frequency, were evaluated during these calculations. As shown in Table 1, the majority of these calculations were performed at Mach The approach of Rumsey et al. [13] was used during the current study. B. Computational Grids The two-dimensional structured grid of Rumsey was used by Cobalt for the 18%-thick circular-arc airfoil validation case. This grid is shown in Figure 25. The C-type grid has 257 points around the airfoil and 81 points normal to the airfoil. The grid had 177 points on the airfoil itself. The grid extends 20 chord lengths in each direction to the farfield. The viscous spacing normal to the airfoil is approximately chord lengths. Unstructured tetrahedral grids of similar size and grid spacing were generated for FUN3D and USM3D. This unstructured grid is shown in Figure 26. C. Results The results for the 18%-thick circular-arc airfoil validation case are presented in Figures 27 through 35. In Figure 27, the normal force coefficient as a function of time is plotted for one cycle of the flow at Mach The Cobalt flow solver with the SST turbulence model was used to generate these results. In this plot, nine instances in time are labeled with the letters A through I. In Figure 28, contours of Mach number are shown in the vicinity of the 18%-thick circular-arc airfoil for the times labeled in Figure 27. This figure shows the unsteadiness of the flow, which results from the shock-induced flow separation which occurs on the airfoil. It is important to point out that the airfoil was stationary during these calculations. The unsteadiness that is seen in the figure is not due to any motion of the airfoil. The snapshots in Figure 28 show that the flow separation on the airfoil oscillates between the upper and lower surfaces. The results in Figures 29 and 30 are for static cases ran at Mach These plots show the impact of time step, sub-iterations and turbulence model on the reduced frequency, which is defined as 2 f (0.5c) k u When calculating the reduced frequency, the pressure component of the lift coefficient is used. 7
8 In Figure 29, the reduced frequency is plotted as a function of time step for the Cobalt, FUN3D and USM3D flow solvers and the SA and SST turbulence models. In general, the reduced frequency converges as the time step is reduced. Converged values of the reduced frequency range from 0.47 to Based on experimental data, the reduced frequency is In Figure 30, the reduced frequency is plotted as a function of the number of sub-iterations. The reduced frequency for Cobalt is converged after 10 sub-iterations. The reduced frequency for USM3D and FUN3D do not appear to be converged after 20 sub-iterations. Once again, based on experimental data the reduced frequency is In the second phase of this study, the ability of the flow solvers to predict the hysteresis effect was evaluated. To do this, an initial time-accurate solution was obtained at Mach 0.7. Once converged, the airfoil was analyzed at Mach 0.71 starting from the final solution of the Mach 0.7 calculations. Once converged, the Mach 0.72 calculations were started using the final solution of the Mach 0.71 case. This process was continued for Mach 0.73, 0.74, 0.75 and Once the Mach 0.76 calculations were converged, the Mach number was reduced to Mach 0.7 at 0.01 increments. Each case at a new Mach number restarted from the solution at the previous Mach number. Time-accurate solutions were generated for each Mach number. The results for the hysteresis calculations are shown in Figures 31 through 33. In these figures, the axial-force, normal-force and pitching-moment coefficients are plotted as a function of time and Mach number for Cobalt, USM3D and FUN3D. The results for both the SA and SST turbulence models are shown in the plots. For increasing Mach number, Cobalt, FUN3D and USM3D predict the onset of unsteadiness between Mach 0.73 and Mach For decreasing Mach number, Cobalt and FUN3D predict the onset of steadiness between Mach 0.73 and Mach 0.72 while USM3D predicts the onset of steadiness between Mach 0.74 and Mach For USM3D, the change between steady flow and unsteady flow is between Mach 0.73 and Mach 0.74 regardless of whether the Mach number is increasing or decreasing. In Figure 34, the results for USM3D at Mach are shown. At this Mach number, the flow is steady if the Mach number is increasing (as shown in the left plot of Figure 34) and unsteady if the Mach number is decreasing (as shown in the right plot of Figure 34). Therefore, USM3D predicts the hysteresis effect, but at a slightly different Mach number than Cobalt and FUN3D. For Cobalt and FUN3D, the flow is steady at Mach 0.73 if the Mach number is increasing and unsteady at Mach 0.73 if the Mach number is decreasing. During a wind-tunnel experiment of this airfoil, unsteady flow began at Mach 0.76 for increasing Mach number and unsteady flow ceased at Mach 0.73 for decreasing Mach number. [14] In Figure 35, the free-stream Mach numbers where regions of unsteady flow were predicted by CFD are plotted for increasing and decreasing Mach number. These results are compared with the unsteady regions seen during the WT test. When analyzing this figure, it is important to note that CFD results were generated between Mach 0.7 and Mach No CFD results were generated above Mach However, WT data was produced above Mach For increasing Mach number, CFD predicts the onset of unsteady flow at a lower Mach number than was seen in the WT test. For CFD, the onset of unsteady flow begins at Mach For the WT test, the onset begins at Mach The reason for this difference between CFD and WT is unknown. For decreasing Mach number, unsteady flow is predicted by CFD to Mach 0.73 or Mach 0.735, depending on the flow solver. Unsteady flow continues during the WT test to Mach There is general agreement between CFD and the WT test in this regard. VIII. NACA 0015 Pitching Airfoil A. Description The low angle of attack pitching NACA 0015 airfoil case of Piziali [15] was also included as a validation case for this project. For this case, the turbulent Navier-Stokes equations were solved. The SA and SST turbulence models were evaluated. This case provided the opportunity to validate each of the CFD codes for a simple rigid body motion. For the angle of attack range being studied here, the flow is attached throughout the majority of the cycle. This allows some of the flow solvers being used during this project to be run in 2-D mode. The flow conditions for this case are shown in Table 1. The angle of attack was prescribed as sin(2 ft) where f is the frequency of oscillation. A frequency of oscillation of 10 cycles per second was used for the calculations. The airfoil is rotated about the quarter-chord point. 8
9 B. Computational Grids The structured grid of Nichols was used by the Cobalt flow solver. This grid is shown in Figure 36. This is a C- type grid with 441 points around the surface and 71 points normal to the surface. The grid had 241 points on the airfoil itself. The grid spacing at the surface in the normal direction is chord lengths. The farfield boundaries of the grid are approximately 25 chord lengths from the airfoil surface. An unstructured tetrahedral grid was generated for use with the USM3D and FUN3D flow solvers. This grid is shown in Figure 37. C. Results In Figure 38, contours of Mach number are shown in the vicinity of the airfoil as the airfoil is pitched from 4 to 8.2 to -0.2 to 4. The results of the Cobalt flow solver and the SST turbulence model were used for this figure. The acceleration of the flow around the leading edge at the higher pitch angles is clearly evident in the Mach contours. The Mach contours also indicate that the flow is attached throughout the majority of the pitch cycle. The lift, drag and pitching-moment coefficients are plotted as a function of pitch angle in Figures 39, 40 and 41. The results for the wind tunnel data and the Cobalt, FUN3D and USM3D flow solvers are included. Each flow solver was run with both the SA and SST turbulence models. In each case, only the pressure component of the forces and moments are being compared. In Figure 39, the lift coefficient is plotted. The Cobalt results compare well to the WT data. The FUN3D results compare well to the WT data at the lower pitch angles. At the higher pitch angles, FUN3D does not predict the hysteresis effect as the results for the up-sweep and down-sweep are nearly identical. The hysteresis effect is predicted with USM3D, but the results under-predict the lift at the higher pitch angles. The differences between the SA and SST results are not significant. The drag coefficient is plotted in Figure 40. Each flow solver accurately predicts the shape of the data. Cobalt and USM3D appear to under-predict the drag, while FUN3D over-predicts the drag. Once again, the differences between the SA and SST models are subtle. In Figure 41, the pitching-moment coefficient is plotted as a function of pitch angle. Each flow solver is able to predict the shape of the plot. While Cobalt appears to be closest to the WT data, the results of each flow solver are rotated from the WT data. It s interesting to note that the FUN3D/SST results predict the roughness of the data. The curves from the other plots are smooth. IX. Summary The goal of this study was to evaluate the time accuracy of five different Navier-Stokes flow solvers using five simple two-dimensional validation cases. The five flow solvers included Cobalt, USM3D, FUN3D, Beggar and NASCART-GT. The five validation cases included the laminar flow over a circular cylinder, the laminar shock tube, the inviscid convecting vortex, the 18%-thick circular-arc airfoil and the NACA 0015 pitching airfoil. This study was conducted in an effort to better understand the capabilities of each flow solver for time-accurate calculations. The eventual goal of the project was to use one or more of these flow solvers to analyze the unsteady transonic flow over a fighter aircraft with stores in captive carriage. Evaluating the capabilities of each flow solver with simple validation cases rather than a complicated fighter aircraft was desirable. While the goal of the project was to evaluate all five flow solvers on all five validation cases, the NASCART-GT flow solver was not successfully validated on any of the validation cases. The remaining four flow solvers were successfully used on the circular cylinder, laminar shock tube and inviscid convecting vortex validation cases. The Cobalt, USM3D and FUN3D flow solvers were successfully used on the 18%-thick circular-arc airfoil and the NACA 0015 pitching airfoil validation cases. Funding, time and flow solver limitations were the main reasons that all five flow solvers were not successfully used on all five validation cases. For the calculations that were completed, the CFD results compare well to the validation data that was provided. There are cases, however, where the CFD results do not compare well to the data. An example of this can be seen in the results of the pitching-moment coefficient for the NACA 0015 pitching airfoil. Since the Cobalt, FUN3D and USM3D flow solvers were successfully used on all five validation cases, these three flow solvers would be the best choices for accurately predicting the transonic, unsteady flow over a fighter aircraft with stores in captive carriage. This conclusion is based on the results of this study alone. Evaluating the abilities of these flow solvers to accurately predict the flow over the fighter aircraft alone and the stores/weapons alone would provide more information on making a decision regarding the flow solver or solvers that could accurately predict the complex problem at hand. 9
10 Acknowledgments The authors would like to thank the Institute for HPC (High Performance Computing) Applications to Air Armament (IHAAA) program for funding this project. The computer time used for the CFD calculations during this project was obtained from the HPC Modernization Program (HPCMP). The authors would like to acknowledge Mr. B.L. Chesser, AFSEO, and Dr. Robert Nichols, UAB, for their help and guidance during this project. References [1] Strang, W. Z., Tomaro, R. F., and Grismer, M. J., The Defining Methods of Cobalt60: a Parallel, Implicit, Unstructured Euler/Navier-Stokes Flow Solver, AIAA , January [2] Spalart P., and Allmaras S. A., One-Equation Turbulence Model for Aerodynamic Flows, AIAA , January [3] Menter, F.R., Improved Two-Equation k-omega Turbulence Models for Aerodynamic Flows, NASA TM , October [4] Frink, N. T., Tetrahedral Unstructured Navier-Stokes Method for Turbulent Flows, AIAA Journal, Vol. 36, No. 11, Nov. 1998, pp [5] Anderson, W. K., Rausch, R. D., and Bonhaus, D. L., Implicit/Multigrid Algorithms for Incompressible Turbulent Flows on Unstructured Grids, Journal of Computational Physics, Vol. 128, Issue 2, 1996, pp [6] Rizk, M., Ellison, S., Prewitt, N., Beggar A Store Separation Predictive Tool, AIAA , June [7] Baldwin, B. S., and Lomax, H., Thin Layer Approximation and Algebraic Model Separated Turbulent Flows, AIAA , January [8] Nichols, R. H., and Heikkinen, B. D., Validation of Implicit Algorithms for Unsteady Flows Including Moving and Deforming Grids, AIAA , January [9] Tritton, D. J., Experiments on the Flow Past a Circular Cylinder at Low Reynolds Numbers, Journal of Fluid Mechanics, Vol. 6, 1959, pp [10] Roshko, A., On the Development of Turbulent Wakes from Vortex Streets, NACA Report 1191, [11] Friehe, C. A., Vortex Shedding from Cylinders at Low Reynolds Numbers, Journal of Fluid Mechanics, Vol. 100, Part 2, 1980, pp [12] McDevitt, J. B., Supercritical Flow about a Thick Circular-Arc Airfoil, NASA TM 78549, January [13] Rumsey, C. L., Sanetrik, M. D., Biedron, R. T., Melson, N. D., and Parlette, E. B., Efficiency and Accuracy of Time- Accurate Turbulent Navier-Stokes Computations, Computers and Fluids, Vol. 25, No. 2, pp , [14] McDevitt, J. B., Levy, L. L., and Deiwert, G. S., Transonic Flow about a Thick Circular-Arc Airfoil, AIAA Journal, Vol. 14, No. 5, May 1976, pp [15] Piziali, R. A., 2-D and 3-D Oscillating Wing Aerodynamics for a Range of Angles of Attack Including Stall, NASA TM 4632, Sep
11 Figures Figure 1. Structured grid of the circular cylinder used for the Cobalt and Beggar calculations Figure 2. Unstructured grid of the circular cylinder used for the USM3D and FUN3D calculations. A B H I C G c l 0 D F E t Figure 3. Lift coefficient as a function of time for the laminar circular cylinder at Mach 0.2, Re c = 150, = 0 for the Cobalt flow solver. 11
12 Time A Time F Time B Time G Time C Time H Time D Time I Time E Figure 4. Mach number contours as a function of time for the laminar circular cylinder at Mach 0.2, Re c = 150, = 0 for the Cobalt flow solver. 12
13 Cobalt USM3D c l Amp FUN3D Beggar c l Amp Sub-iterations Small Time Step Sub-iterations Large Time Step Figure 5. Amplitude of lift coefficient versus number of subiterations for the circular cylinder validation case. Cobalt USM3D St (lift) FUN3D Beggar St (lift) Sub-iterations Small Time Step Sub-iterations Large Time Step Figure 6. Lift Strouhal number versus number of subiterations for the circular cylinder validation case. 13
14 Cobalt USM3D c d Avg FUN3D Beggar c d Avg Sub-iterations Small Time Step Sub-iterations Large Time Step Figure 7. Average drag coefficient versus number of subiterations for the circular cylinder validation case. Cobalt USM3D c d Amp FUN3D Beggar c d Amp Sub-iterations Small Time Step Sub-iterations Large Time Step Figure 8. Amplitude of drag coefficient versus number of subiterations for the circular cylinder validation case. 14
15 Cobalt USM3D St (drag) FUN3D Beggar St (drag) Sub-iterations Small Time Step Sub-iterations Large Time Step Figure 9. Drag Strouhal number versus number of subiterations for the circular cylinder validation case. 15
16 Figure 10. Coarse Grid 41x41 Medium Grid 81x81 Fine Grid 161x161 Coarse, medium and fine structured grids used for the Cobalt convecting vortex calculations. Coarse Grid Medium Grid Fine Grid Figure 11. Coarse, medium and fine unstructured grids used for the USM3D and FUN3D convecting vortex calculations. 16
17 Initial Vortex After 5 Cycles: Fine Grid Medium Grid Coarse Grid Figure 12. Initial Vortex Results for Cobalt for the inviscid convecting vortex. After 5 Cycles: Fine Grid Medium Grid Coarse Grid Figure 13. Results for Beggar for the inviscid convecting vortex. Initial Vortex After 5 Cycles: Fine Grid Medium Grid Coarse Grid Figure 14. Results for FUN3D for the inviscid convecting vortex. Initial Vortex After 5 Cycles: Fine Grid Figure 15. Results for USM3D for the inviscid convecting vortex. 17
18 Cobalt USM3D Min p/p FUN3D Beggar Min p/p Cycles Cycles Coarse Grid Medium Grid Fine Grid Figure 16. Minimum pressure versus cycles for the inviscid convecting vortex. NASCART Min p/p Cycles Coarse Grid Medium Grid Fine Grid Figure 17. Minimum pressure versus cycles for the inviscid convecting vortex for NASCART. 18
19 Figure 18. Structured grid of the shock tube used for the Cobalt and Beggar calculations. Figure 19. Unstructured grid of the shock tube used for the USM3D and FUN3D calculations. 19
20 Cobalt USM3D FUN3D Beggar x x Exact Large t Medium t Small t Figure 20. Density in the shock tube after a non-dimensional time of 0.2. Cobalt USM3D p FUN3D Beggar p x x Exact Large t Medium t Small t Figure 21. Pressure in the shock tube after a non-dimensional time of
21 Cobalt USM3D u FUN3D Beggar u x x Exact Large t Medium t Small t Figure 22. Velocity in the shock tube after a non-dimensional time of 0.2. Cobalt USM3D e FUN3D Beggar e x x Exact Large t Medium t Small t Figure 23. Energy in the shock tube after a non-dimensional time of
22 Cobalt Nichols in AIAA Ch x x Figure 24. Wall heat coefficient in the shock tube after a non-dimensional time of 0.2. Figure 25. Structured grid of the 18%-thick circular-arc airfoil used for the Cobalt calculations. 22
23 Figure 26. Untructured grid of the 18%-thick circular-arc airfoil used for the USM3D and FUN3D calculations. D E C F c n A B G H I t Figure 27. Normal force coefficient as a function of time for the 18%-thick circular-arc airfoil at Mach 0.76, Re c = 12 million, = 0 for the Cobalt flow solver and the SST turbulence model. 23
24 Time A Time F Time B Time G Time C Time H Time D Time I Time E Figure 28. Mach number contours as a function of time for the 18%-thick circular-arc airfoil at Mach 0.76, Re c = 12 million, = 0 with the Cobalt flow solver and the SST turbulence model. 24
25 Cobalt FUN3D USM3D k t t t SA SST Figure 29. Reduced frequency versus time step for the 18%-thick circular-arc airfoil at Mach 0.76, Re c = 12 million, = 0. Cobalt FUN3D USM3D k Sub-iterations Sub-iterations Sub-iterations SA SST Figure 30. Reduced frequency versus subiterations for the 18%-thick circular-arc airfoil at Mach 0.76, Re c = 12 million, = 0. 25
26 Cobalt with SA FUN3D with SA USM3D with SA c a M Cobalt with SST FUN3D with SST USM3D with SST c a M Figure 31. t Axial Force t Mach Number Axial-force coefficient versus time and Mach number for the 18%-thick circulararc airfoil at = 0. t Cobalt with SA FUN3D with SA USM3D with SA c n M Cobalt with SST FUN3D with SST USM3D with SST c n M Figure 32. t Normal Force t Mach Number Normal-force coefficient versus time and Mach number for the 18%-thick circulararc airfoil at = 0. t 26
27 Cobalt with SA FUN3D with SA USM3D with SA c m M Cobalt with SST FUN3D with SST USM3D with SST c m M Figure 33. t Pitching Moment t Mach Number Pitching-moment coefficient versus time and Mach number for the 18%-thick circular-arc airfoil at = 0. t Mach Increasing Mach Mach Sweep Mach Decreasing Mach c a M t t t Axial Force Mach Number Figure 34. Axial-force coefficient versus time and Mach number for the 18%-thick circulararc airfoil for USM3D with the SA turbulence model at = 0. Increasing Mach Decreasing Mach USM3D Cobalt FUN3D WT Data Mach Number Mach Number Figure 35. Regions of unsteady flow predicted by wind-tunnel experiment and the Cobalt, FUN3D, USM3D flow solvers for the 18%-thick circular-arc airfoil at = 0. 27
28 Figure 36. Structured grid of the NACA 0015 pitching airfoil that was used for the Cobalt calculations. Figure 37. Unstructured grid of the NACA 0015 pitching airfoil that was used for the USM3D and FUN3D calculations. 28
29 = 4, increasing = 1.03, decreasing = 6.97, increasing = -0.2 = 8.2 = 1.03, increasing = 6.97, decreasing = 4, increasing = 4, decreasing Figure 38. Mach number contours as a function of pitch angle for the NACA 0015 pitching airfoil at Mach 0.29, Re c = 1.95 million with the Cobalt flow solver and the SST turbulence model. 29
30 Cobalt with SA FUN3D with SA USM3D with SA c l Cobalt with SST FUN3D with SST USM3D with SST c l Figure 39. Pitch Angle Pitch Angle Pitch Angle WT Data Small Time Step Large Time Step Lift coefficient as a function of pitch angle for the NACA 0015 pitching airfoil. Cobalt with SA FUN3D with SA USM3D with SA c d Cobalt with SST FUN3D with SST USM3D with SST c d Figure 40. Pitch Angle Pitch Angle Pitch Angle WT Data Small Time Step Large Time Step Drag coefficient as a function of pitch angle for the NACA 0015 pitching airfoil. Cobalt with SA FUN3D with SA USM3D with SA c m Cobalt with SST FUN3D with SST USM3D with SST c m Pitch Angle WT Data Pitch Angle Small Time Step Pitch Angle Large Time Step Figure 41. Pitching-moment coefficient as a function of pitch angle for the NACA 0015 pitching airfoil. 30
Introduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
Application of Wray-Agarwal Turbulence Model for Accurate Numerical Simulation of Flow Past a Three-Dimensional Wing-body
Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 4-28-2016 Application
Introduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
Debojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons Stores Separation Introduction Flight Test Expensive, high-risk, sometimes catastrophic
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
High-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder
High-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder Aerospace Application Areas Aerodynamics Subsonic through Hypersonic Aeroacoustics Store release & weapons bay analysis High lift devices
LES Applications in Aerodynamics
LES Applications in Aerodynamics Kyle D. Squires Arizona State University Tempe, Arizona, USA 2010 Tutorial School on Fluid Dynamics: Topics in Turbulence Center for Scientific Computation and Mathematical
A Computational Study of the Abrupt Wing Stall (AWS) Characteristics for Various Fighter Jets: Part I, F/A-18E and F-16C
AIAA 2003-0746 A Computational Study of the Abrupt Wing Stall (AWS) Characteristics for Various Fighter Jets: Part I, F/A-18E and F-16C Paresh Parikh NASA Langley Research Center, Hampton, VA James Chung
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 -
SPC 307 Aerodynamics. Lecture 1. February 10, 2018
SPC 307 Aerodynamics Lecture 1 February 10, 2018 Sep. 18, 2016 1 Course Materials drahmednagib.com 2 COURSE OUTLINE Introduction to Aerodynamics Review on the Fundamentals of Fluid Mechanics Euler and
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN ABSTRACT
AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN Adeel Khalid *, Daniel P. Schrage + School of Aerospace Engineering, Georgia Institute of Technology
ANSYS FLUENT. Airfoil Analysis and Tutorial
ANSYS FLUENT Airfoil Analysis and Tutorial ENGR083: Fluid Mechanics II Terry Yu 5/11/2017 Abstract The NACA 0012 airfoil was one of the earliest airfoils created. Its mathematically simple shape and age
COMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY. Abstract. I. Introduction
COMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY John R. Cipolla 709 West Homeway Loop, Citrus Springs FL 34434 Abstract A series of computational fluid dynamic (CFD)
Modeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
The Spalart Allmaras turbulence model
The Spalart Allmaras turbulence model The main equation The Spallart Allmaras turbulence model is a one equation model designed especially for aerospace applications; it solves a modelled transport equation
Application of STAR-CCM+ to Helicopter Rotors in Hover
Application of STAR-CCM+ to Helicopter Rotors in Hover Lakshmi N. Sankar and Chong Zhou School of Aerospace Engineering, Georgia Institute of Technology, Atlanta, GA Ritu Marpu Eschol CD-Adapco, Inc.,
Investigation of the Influence of the Turbulent Transition on the Transonic Periodic Flow
Investigation of the Influence of the Turbulent Transition on the Transonic Periodic Flow C. Tulita *, E. Turkbeyler and S. Raghunathan Department of Aeronautical Engineering, The Quenn s University of
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.
Estimating Vertical Drag on Helicopter Fuselage during Hovering
Estimating Vertical Drag on Helicopter Fuselage during Hovering A. A. Wahab * and M.Hafiz Ismail ** Aeronautical & Automotive Dept., Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310
Computational Fluid Dynamics for Engineers
Tuncer Cebeci Jian P. Shao Fassi Kafyeke Eric Laurendeau Computational Fluid Dynamics for Engineers From Panel to Navier-Stokes Methods with Computer Programs With 152 Figures, 19 Tables, 84 Problems and
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson
AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS
AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS Emre Alpman Graduate Research Assistant Aerospace Engineering Department Pennstate University University Park, PA, 6802 Abstract A new methodology
Ail implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations
Ail implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations Vimala Dutta Computational and Theoretical Fluid Dynamics Division National Aerospace
Profile Catalogue for Airfoil Sections Based on 3D Computations
Risø-R-58(EN) Profile Catalogue for Airfoil Sections Based on 3D Computations Franck Bertagnolio, Niels N. Sørensen and Jeppe Johansen Risø National Laboratory Roskilde Denmark December 26 Author: Franck
Usage of CFX for Aeronautical Simulations
Usage of CFX for Aeronautical Simulations Florian Menter Development Manager Scientific Coordination ANSYS Germany GmbH Overview Elements of CFD Technology for aeronautical simulations: Grid generation
Abstract. Introduction
EULER SOLUTIONS AS LIMIT OF INFINITE REYNOLDS NUMBER FOR SEPARATION FLOWS AND FLOWS WITH VORTICES Wolfgang Schmidt and Antony Jameson Dornier GmbH, D-7990 Friedrichshafen, FRG and Princeton University,
Detached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads
Detached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads Matt Knapp Chief Aerodynamicist TLG Aerospace, LLC Presentation Overview Introduction to TLG Aerospace
Rotorcraft Noise Prediction with Multi-disciplinary Coupling Methods. Yi Liu NIA CFD Seminar, April 10, 2012
Rotorcraft Noise Prediction with Multi-disciplinary Coupling Methods Yi Liu NIA CFD Seminar, April 10, 2012 Outline Introduction and Background Multi-disciplinary Analysis Approaches Computational Fluid
NSU3D Results for the Fourth AIAA Drag Prediction Workshop
NSU3D Results for the Fourth AIAA Drag Prediction Workshop Dimitri Mavriplis Mike Long Department of Mechanical Engineering University of Wyoming, Laramie, Wyoming, 82072-3295 Simulation results for the
FAR-Wake Workshop, Marseille, May 2008
Wake Vortices generated by an Aircraft Fuselage : Comparison of Wind Tunnel Measurements on the TAK Model with RANS and RANS-LES Simulations T. Louagie, L. Georges & P. Geuzaine Cenaero CFD-Multiphysics
Missile External Aerodynamics Using Star-CCM+ Star European Conference 03/22-23/2011
Missile External Aerodynamics Using Star-CCM+ Star European Conference 03/22-23/2011 StarCCM_StarEurope_2011 4/6/11 1 Overview 2 Role of CFD in Aerodynamic Analyses Classical aerodynamics / Semi-Empirical
Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit
Washington University in St. Louis Washington University Open Scholarship Mechanical Engineering and Materials Science Independent Study Mechanical Engineering & Materials Science 12-21-2016 Computational
Introduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
A DRAG PREDICTION VALIDATION STUDY FOR AIRCRAFT AERODYNAMIC ANALYSIS
A DRAG PREDICTION VALIDATION STUDY FOR AIRCRAFT AERODYNAMIC ANALYSIS Akio OCHI, Eiji SHIMA Kawasaki Heavy Industries, ltd Keywords: CFD, Drag prediction, Validation Abstract A CFD drag prediction validation
Computation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables
Computation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables By Veer N. Vatsa NASA Langley Research Center, Hampton, VA 23681 Mail Stop 128, email: v.n.vatsa@larc.nasa.gov
Driven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
The viscous forces on the cylinder are proportional to the gradient of the velocity field at the
Fluid Dynamics Models : Flow Past a Cylinder Flow Past a Cylinder Introduction The flow of fluid behind a blunt body such as an automobile is difficult to compute due to the unsteady flows. The wake behind
APPROACH FOR NUMERICAL MODELING OF AIRFOIL DYNAMIC STALL
APPROACH FOR NUMERICAL MODELING OF AIRFOIL DYNAMIC STALL Cvetelina Velkova, Ivan Dobrev, Michael Todorov, Fawaz Massouh To cite this version: Cvetelina Velkova, Ivan Dobrev, Michael Todorov, Fawaz Massouh.
Numerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
Keywords: flows past a cylinder; detached-eddy-simulations; Spalart-Allmaras model; flow visualizations
A TURBOLENT FLOW PAST A CYLINDER *Vít HONZEJK, **Karel FRAŇA *Technical University of Liberec Studentská 2, 461 17, Liberec, Czech Republic Phone:+ 420 485 353434 Email: vit.honzejk@seznam.cz **Technical
Second Symposium on Hybrid RANS-LES Methods, 17/18 June 2007
1 Zonal-Detached Eddy Simulation of Transonic Buffet on a Civil Aircraft Type Configuration V.BRUNET and S.DECK Applied Aerodynamics Department The Buffet Phenomenon Aircraft in transonic conditions Self-sustained
Aerodynamic Analysis of Forward Swept Wing Using Prandtl-D Wing Concept
Aerodynamic Analysis of Forward Swept Wing Using Prandtl-D Wing Concept Srinath R 1, Sahana D S 2 1 Assistant Professor, Mangalore Institute of Technology and Engineering, Moodabidri-574225, India 2 Assistant
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary
Aerodynamic Design Optimization of UAV Rotor Blades using a Genetic Algorithm
Aerodynamic Design Optimization of UAV Rotor Blades using a Genetic Algorithm Hak-Min Lee 1), Nahm-Keon Hur 2) and *Oh-Joon Kwon 3) 1), 3) Department of Aerospace Engineering, KAIST, Daejeon 305-600, Korea
Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon Kim, Hogeon Kim
Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon Kim, Hogeon Kim Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon
Grid. Apr 09, 1998 FLUENT 5.0 (2d, segregated, lam) Grid. Jul 31, 1998 FLUENT 5.0 (2d, segregated, lam)
Tutorial 2. Around an Airfoil Transonic Turbulent Flow Introduction: The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a non-zero angle of attack. You will use the
Analysis of an airfoil
UNDERGRADUATE RESEARCH FALL 2010 Analysis of an airfoil using Computational Fluid Dynamics Tanveer Chandok 12/17/2010 Independent research thesis at the Georgia Institute of Technology under the supervision
An Upwind Multigrid Method for Solving Viscous Flows on Unstructured Triangular Meshes
An Upwind Multigrid Method for Solving Viscous Flows on Unstructured Triangular Meshes by Daryl Lawrence Bonhaus B.S. June 1990, University of Cincinnati A Thesis submitted to The Faculty of The School
Team 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud
Team 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud CFD Support s OpenFOAM and UberCloud Containers enable efficient, effective, and easy access and use of MEET THE TEAM End-User/CFD
Efficient Unstructured Mesh Flow Solver Using Hamiltonian Paths and Strand Grids
Efficient Unstructured Mesh Flow Solver Using Hamiltonian Paths and Strand Grids *Yong Su Jung 1), Bharath Govindarajan 2) and James Baeder 3) 1), 2), 3) Department of Aerospace Engineering, University
Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
Progress and Future Prospect of CFD in Aerospace
Progress and Future Prospect of CFD in Aerospace - Observation from 30 years research - Kozo Fujii Institute of Space and Astronautical Science (ISAS) Japan Aerospace Exploration Agency (JAXA) Japan JAXA:
Challenges in Boundary- Layer Stability Analysis Based On Unstructured Grid Solutions
Challenges in Boundary- Layer Stability Analysis Based On Unstructured Grid Solutions Wei Liao National Institute of Aerospace, Hampton, Virginia Collaborators: Mujeeb R. Malik, Elizabeth M. Lee- Rausch,
CFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence
CFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence Kavya H.P, Banjara Kotresha 2, Kishan Naik 3 Dept. of Studies in Mechanical Engineering, University BDT College of Engineering,
A MESH ADAPTATION METHOD FOR SIMULATION OF UNSTEADY FLOWS
A MESH ADAPTATION METHOD FOR SIMULATION OF UNSTEAD FLOWS C. H. Zhou* * Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing, 6, China Keywords: mesh adaptation, unsteady
39th AIAA Aerospace Sciences Meeting and Exhibit January 8 11, 2001/Reno, NV
AIAA 1 717 Static Aero-elastic Computation with a Coupled CFD and CSD Method J. Cai, F. Liu Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697-3975 H.M. Tsai,
AIAA A Solution on the F-18C for Store Separation Simulation Using Cobalt 60 R. F. Tomaro, F. C. Witzeman and W. Z. Strang Computational
AIAA 99-0122 A Solution on the F-18C for Store Separation Simulation Using Cobalt 60 R. F. Tomaro, F. C. Witzeman and W. Z. Strang Computational Sciences Branch Air Vehicles Directorate Air Force Research
Strömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
Axisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
Simulation of Turbulent Flow in an Asymmetric Diffuser
Simulation of Turbulent Flow in an Asymmetric Diffuser 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 3 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University of Iowa C.
The Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
NUMERICAL SIMULATION OF 3D FLAPPING WING BASED ON CHIMERA METHOD
26 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES NUMERICAL SIMULATION OF 3D FLAPPING WING Wenqing Yang, Bifeng Song, Wenping Song School of Aeronautics, Northwestern Polytechnical University,
Solution Recording and Playback: Vortex Shedding
STAR-CCM+ User Guide 6663 Solution Recording and Playback: Vortex Shedding This tutorial demonstrates how to use the solution recording and playback module for capturing the results of transient phenomena.
An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment
The Aeronautical Journal November 2015 Volume 119 No 1221 1451 An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment J. A. Camberos
AIAA A Computational Study of Subsonic Flows over A Medium Range Cargo Aircraft
S AIAA 23-366 A Computational Study of Subsonic Flows over A Medium Range Cargo Aircraft 2st Applied Aerodynamic Conference 23-26 June 23 Orlando, Florida For permission to copy or to republish, contact
Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software
Reports of Research Institute for Applied Mechanics, Kyushu University No.150 (71 83) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 3: For the Case
Keywords: CFD, aerofoil, URANS modeling, flapping, reciprocating movement
L.I. Garipova *, A.N. Kusyumov *, G. Barakos ** * Kazan National Research Technical University n.a. A.N.Tupolev, ** School of Engineering - The University of Liverpool Keywords: CFD, aerofoil, URANS modeling,
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI
Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI François Thirifay and Philippe Geuzaine CENAERO, Avenue Jean Mermoz 30, B-6041 Gosselies, Belgium Abstract. This paper reports
APPLICATION OF A NAVIER-STOKES SOLVER TO THE ANALYSIS OF MULTIELEMENT AIRFOILS AND WINGS USING MULTIZONAL GRID TECHNIQUES
APPLICATION OF A NAVIER-STOKES SOLVER TO THE ANALYSIS OF MULTIELEMENT AIRFOILS AND WINGS USING MULTIZONAL GRID TECHNIQUES Kenneth M. Jones* NASA Langley Research Center Hampton, VA Robert T. Biedron Analytical
THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS
March 18-20, 2013 THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS Authors: M.R. Chiarelli, M. Ciabattari, M. Cagnoni, G. Lombardi Speaker:
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
HPCMP CREATE TM -AV Kestrel Dual Mesh Computations on the ROBIN Fuselage
HPCMP CREATE TM -AV Kestrel Dual Mesh Computations on the ROBIN Fuselage Jennifer Abras NAVAIR Applied Aerodynamics and Store Separation Branch Nathan Hariharan HPCMP CREATE-AV 12 th Overset Grid Symposium
Research Directions in. Unsteady Aerodynamics
APS DFD 2004 Seattle Research Directions in Unsteady Aerodynamics P. Spalart* Boeing Commercial Airplanes *with Strelets and Squires Outline Theory: )-: Computational Fluid Dynamics Resolved and modeled
Numerical Study of Blending Hybrid RANS/LES Method and γ-re θ Transition Model for Unsteady Turbulent Flow Analysis
Numerical Study of Blending Hybrid RANS/LES Method and γ-re θ Transition Model for Unsteady Turbulent Flow Analysis *Jeong Hwan Sa 1), Kum Won Cho 2) and Soo Hyung Park 3) 1) Department of Computational
Multidisciplinary Applications of Detached-Eddy Simulation to Separated Flows at High Reynolds Numbers (Challenge 92)
Multidisciplinary Applications of Detached-Eddy Simulation to Separated Flows at High Reynolds Numbers (Challenge 92) Scott A. Morton USAF Academy Scott.Morton@usafa.af.mil Denis B. Kholodar USAF Academy
Numerical Investigation of Transonic Shock Oscillations on Stationary Aerofoils
Numerical Investigation of Transonic Shock Oscillations on Stationary Aerofoils A. Soda, T. Knopp, K. Weinman German Aerospace Center DLR, Göttingen/Germany Symposium on Hybrid RANS-LES Methods Stockholm/Sweden,
ITU/FAA Faculty of Aeronautics and Astronautics
S. Banu YILMAZ, Mehmet SAHIN, M. Fevzi UNAL, Istanbul Technical University, 34469, Maslak/Istanbul, TURKEY 65th Annual Meeting of the APS Division of Fluid Dynamics November 18-20, 2012, San Diego, CA
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS
COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS Dr W. Malalasekera Version 3.0 August 2013 1 COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE
Direct numerical simulations of flow and heat transfer over a circular cylinder at Re = 2000
Journal of Physics: Conference Series PAPER OPEN ACCESS Direct numerical simulations of flow and heat transfer over a circular cylinder at Re = 2000 To cite this article: M C Vidya et al 2016 J. Phys.:
Modelling Flow Phenomena in Time Dependent Store Release from Transonic Aircraft
Modelling Flow Phenomena in Time Dependent Store Release from Transonic Aircraft D A MacLucas and I M A Gledhill Defence, Peace, Safety and Security Unit, CSIR, PO Box 395, Pretoria 0001 South Africa E-mail:
EXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS
EXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS Edisson Sávio de Góes Maciel, edissonsavio@yahoo.com.br Mechanical
Simulation of Turbulent Flow around an Airfoil
1. Purpose Simulation of Turbulent Flow around an Airfoil ENGR:2510 Mechanics of Fluids and Transfer Processes CFD Lab 2 (ANSYS 17.1; Last Updated: Nov. 7, 2016) By Timur Dogan, Michael Conger, Andrew
EVALUATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE URBAN FLOW
EVALATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE RBAN FLOW Jarkko Saloranta and Antti Hellsten Helsinki niversity of Technology, Laboratory of Aerodynamics, Finland INTRODCTION In this work we study
Introduction to Aerodynamic Shape Optimization
Introduction to Aerodynamic Shape Optimization 1. Aircraft Process 2. Aircraft Methods a. Inverse Surface Methods b. Inverse Field Methods c. Numerical Optimization Methods Aircraft Process Conceptual
Numerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Numerical Simulation Study on Aerodynamic Characteristics of the High Speed Train under Crosswind Fan Zhao,
OpenFOAM GUIDE FOR BEGINNERS
OpenFOAM GUIDE FOR BEGINNERS Authors This guide has been developed by: In association with: Pedro Javier Gamez and Gustavo Raush The Foam House Barcelona ETSEIAT-UPC June 2014 2 OPENFOAM GUIDE FOR BEGINNERS
Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software
Reports of Research Institute for Applied Mechanics, Kyushu University, No.150 (60-70) March 2016 Reproducibility of Complex Turbulent Flow Using Commercially-Available CFD Software Report 2: For the Case
Comparison of CFD Simulation of a Hyundai I20 Model with Four Different Turbulence Models
RESEARCH ARTICLE OPEN ACCESS Comparison of CFD Simulation of a Hyundai I20 with Four Different Turbulence s M. Vivekanandan*, R. Sivakumar**, Aashis. S. Roy*** *(Uttam Industrial Engg. Pvt. Ltd., Tiruchirapalli,
Finite Element Flow Simulations of the EUROLIFT DLR-F11 High Lift Configuration
Finite Element Flow Simulations of the EUROLIFT DLR-F11 High Lift Configuration Kedar C. Chitale MANE Dept., Rensselaer Polytechnic Institute, NY 12180 Michel Rasquin Leadership Computing Facility, Argonne
Estimation of Flow Field & Drag for Aerofoil Wing
Estimation of Flow Field & Drag for Aerofoil Wing Mahantesh. HM 1, Prof. Anand. SN 2 P.G. Student, Dept. of Mechanical Engineering, East Point College of Engineering, Bangalore, Karnataka, India 1 Associate
Effect of Position of Wall Mounted Surface Protrusion in Drag Characteristics At Low Reynolds Number
ISSN (e): 2250 3005 Volume, 07 Issue, 11 November 2017 International Journal of Computational Engineering Research (IJCER) Effect of Position of Wall Mounted Surface Protrusion in Drag Characteristics
Research Article A Computational Investigation of Unsteady Aerodynamics of Insect-Inspired Fixed Wing Micro Aerial Vehicle s 2D Airfoil
Advances in Aerospace Engineering, Article ID 5449, 7 pages http://dx.doi.org/1.1155/214/5449 Research Article A Computational Investigation of Unsteady Aerodynamics of Insect-Inspired Fixed Wing Micro
Case C3.1: Turbulent Flow over a Multi-Element MDA Airfoil
Case C3.1: Turbulent Flow over a Multi-Element MDA Airfoil Masayuki Yano and David L. Darmofal Aerospace Computational Design Laboratory, Massachusetts Institute of Technology I. Code Description ProjectX
ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent
MEGR 7090-003, Computational Fluid Dynamics :1 7 Spring 2015 Computational Study of Laminar Flowfield around a Square Cylinder using Ansys Fluent Rahul R Upadhyay Master of Science, Dept of Mechanical
Algorithmic Developments in TAU
Algorithmic Developments in TAU Ralf Heinrich, Richard Dwight, Markus Widhalm, and Axel Raichle DLR Institute of Aerodynamics and Flow Technology, Lilienthalplatz 7, 38108, Germany ralf.heinrich@dlr.de,