An interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of-magnitude analysis of the subsonic three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first-order continuity effects is accomplished using a penalty differential constraint concept within a finite-element solution algorithm. A Reynolds stress constitutive equation, with low-turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at MM =0.08 and #e/C=0.6xl0/m. Additional results document the turbulence closure model.