Engineering topics
Publications and source records attributed to Ekaterinaris, John A..
The prediction of separation bubbles on NACA 65-213 and NACA 0012 using a modified Chen-Thyson transition model is presented. The contents include: 1) Background; 2) Analysis of NACA 65-213 separation bubble using cebeci's viscous-inviscid interaction method; 3) Analysis of NACA 0012 separation bubble using navier-stokes method; and 4) Comparison with experiment.
The ability of one- and two-equation turbulence models to predict unsteady separated flows over airfoils is evaluated. An implicit, factorized, upwind-biased numerical scheme is used for the integration of the compressible, Reynolds averaged Navier-Stokes equations. The turbulent eddy viscosity is obtained from the computed mean flowfield by integration of the turbulent field equations. The two-equation turbulence models are discretized in space with an upwind-biased, second order accurate total variation diminishing scheme. One and two-equation turbulence models are first tested for a separated airfoil flow at fixed angle of incidence. The same models are then applied to compute the unsteady flowfields about airfoils undergoing oscillatory motion at low subsonic Mach numbers. Experimental cases where the flow has been tripped at the leading edge and where natural transition was allowed to occur naturally are considered. The more recently developed field-equation turbulence models capture the physics of unsteady separated flow significantly better than the standard kappa-epsilon and kappa-omega models. However, certain differences in the hysteresis effects are obtained. For an untripped high-Reynolds-number flow, it was found necessary to take into account the leading edge transitional flow region in order to capture the correct physical mechanism that leads to dynamic stall.
A Navier-Stokes/potentia/flow interactive solution method suitable for the solution of steady-state and unsteady flowfields around airfoils has been developed. The Navier-Stokes equations are solved in the close proximity of an airfoil and in its wake. The inviscid flowfield surrounding the viscous flow regions is assumed to be irrotational and isentropic. In the inviscid flow region, the potential flow equations, which are based on distributed sources and vortices, are solved. The two solutions are strongly coupled through the application of their boundary conditions. The strongly coupled Navier-Stokes/potential flow interactive solution method confines computations to a small domain in the proximity of an airfoil. For steady flow solutions, the computational domain may be confined to a region which extends less than one-fifth of a chord length distance away from the airfoil surface. Computed solutions have the same accuracy as the full domain NavierStokes solutions. Yet, as a result of the significantly reduced computational domain and increased convergence rates, the Navier-Stokes/potential flow interactive solution method is about 40 percent more efficient computationally.
Euler and Navier-Stokes solutions for airfoil flows on zonal grids are presented. The governing equations are solved with an implicit, iterative, factorized numerical scheme. The inviscid fluxes are examined with a third-order accurate upwind method. Zonal grid solutions are compared with experimental measurements for flows over airfoils at fixed angles of incidence. The computed unsteady solutions for rapidly pitching and oscillating airfoils are in good agreement with experiments.
An upwind-biased implicit scheme is used to investigate steady-state and unsteady Navier-Stokes solutions of the vortical flow over a double-delta wing configuration. The governing equations are solved numerically with a fully upwind, implicit, iterative, and factorized numerical scheme. Steady-state solutions for fixed angles of attack and unsteady solutions for a sinusoidal oscillatory motion are obtained. The steady-state solutions on the baseline grid are in agreement with the experiment, and grid refinements show some improvements of the predictions. The higher-order accuracy of the present scheme yields equivalent solutions on smaller grid densities compared to solutions obtained with a second-order accurate method on larger grids. As the angle of attack increases, the grid resolution requirements for adequate resolution of the leeward-side vortical flowfield become very severe. The unsteady solutions are in general agreement with the measurements and show a qualitative correlation with the experiment.
The vortical flowfield over a double-delta wing configuration, consisting of a sharp leading edge 76 deg sweep strake and a 40 deg sweep wing section is investigated. The governing equations are solved numerically with a partially upwind, finite-difference, factorized algorithm. The leeward side vortex system resulting from the strake and wing vortices is investigated for subsonic, high-Reynolds-number flow at various angles of incidence. At low angles of attack the strake and wing vortices remain separate over the wing section, while for flows at higher angles of attack the two vortices merge and vortex breakdown develops. Vortex breakdown appears initially at the trailing edge region of the wing section. As the angle of attack increases bursting occurs further upstream closer to the strake section. The effect of numerical grid density is investigated, and the solutions are compared with available experimental data.