Computational aspects of the prediction of multidimensional transonic flows in turbomachinery

The analytical prediction and description of transonic flow in turbomachinery is complicated by three fundamental effects: (1) the fluid equations describing the transonic regime are inherently nonlinear, (2) shock waves may be present in the flow, and (3) turbomachine blading is geometrically complex, possessing large amounts of curvature, stagger, and twist. A three-dimensional computation procedure for the study of transonic turbomachine fluid mechanics is described. The fluid differential equations and corresponding difference operators are presented, the boundary conditions for complex blade shapes are described, and the computational implementation and mapping procedures are developed. Illustrative results of a typical unthrottled transonic rotor are also presented

Leading-edge vortex research: Some nonplanar concepts and current challenges

Some background information is provided for the Vortex Flow Aerodynamics Conference and that current slender wing airplanes do not use variable leading edge geometry to improve transonic drag polar is shown. Highlights of some of the initial studies combining wing camber, or flaps, with vortex flow are presented. Current vortex flap studies were reviewed to show that there is a large subsonic data base and that transonic and supersonic generic studies have begun. There is a need for validated flow field solvers to calculate vortex/shock interactions at transonic and supersonic speeds. Many important research opportunities exist for fundamental vortex flow investigations and for designing advanced fighter concepts

Finite element analysis of transonic flows in cascades: Importance of computational grids in improving accuracy and convergence

The finite element method is applied for the solution of transonic potential flows through a cascade of airfoils. Convergence characteristics of the solution scheme are discussed. Accuracy of the numerical solutions is investigated for various flow regions in the transonic flow configuration. The design of an efficient finite element computational grid is discussed for improving accuracy and convergence

Extension of transonic flow computational concepts in the analysis of cavitated bearings

An analogy between the mathematical modeling of transonic potential flow and the flow in a cavitating bearing is described. Based on the similarities, characteristics of the cavitated region and jump conditions across the film reformation and rupture fronts are developed using the method of weak solutions. The mathematical analogy is extended by utilizing a few computational concepts of transonic flow to numerically model the cavitating bearing. Methods of shock fitting and shock capturing are discussed. Various procedures used in transonic flow computations are adapted to bearing cavitation applications, for example, type differencing, grid transformation, an approximate factorization technique, and Newton's iteration method. These concepts have proved to be successful and have vastly improved the efficiency of numerical modeling of cavitated bearings

The application of Riegels' rule and time-like damping to transonic flow calculations

Finite difference relaxation solutions of the nonlinear small perturbation equations have proven reliable and successful in determining the transonic flowfields about thin airfoils. However, application of the small perturbation approach to thick airfoils usually results in an accuracy less than desirable. The incorporation of Riegels' Rule and time-like damping into the small perturbation approach and their application to thick and thin airfoils in transonic flow are discussed. Studies for thick and thin airfoils are presented. It is concluded that Riegels' Rule and damping should both be included in small perturbation transonic flow calculations

Solution of steady and unsteady transonic-vortex flows using Euler and full-potential equations

Two methods are presented for inviscid transonic flows: unsteady Euler equations in a rotating frame of reference for transonic-vortex flows and integral solution of full-potential equation with and without embedded Euler domains for transonic airfoil flows. The computational results covered: steady and unsteady conical vortex flows; 3-D steady transonic vortex flow; and transonic airfoil flows. The results are in good agreement with other computational results and experimental data. The rotating frame of reference solution is potentially efficient as compared with the space fixed reference formulation with dynamic gridding. The integral equation solution with embedded Euler domain is computationally efficient and as accurate as the Euler equations

The role of flow geometry in influencing the stability criteria for low angular momentum axisymmetric black hole accretion

Using mathematical formalism borrowed from dynamical systems theory, a complete analytical investigation of the critical behaviour of the stationary flow configuration for the low angular momentum axisymmetric black hole accretion provides valuable insights about the nature of the phase trajectories corresponding to the transonic accretion in the steady state, without taking recourse to the explicit numerical solution commonly performed in the literature to study the multi-transonic black hole accretion disc and related astrophysical phenomena. Investigation of the accretion flow around a non rotating black hole under the influence of various pseudo-Schwarzschild potentials and forming different geometric configurations of the flow structure manifests that the general profile of the parameter space divisions describing the multi-critical accretion is roughly equivalent for various flow geometries. However, a mere variation of the polytropic index of the flow cannot map a critical solution from one flow geometry to the another, since the numerical domain of the parameter space responsible to produce multi-critical accretion does not undergo a continuous transformation in multi-dimensional parameter space. The stationary configuration used to demonstrate the aforementioned findings is shown to be stable under linear perturbation for all kind of flow geometries, black hole potentials, and the corresponding equations of state used to obtain the critical transonic solutions. Finally, the structure of the acoustic metric corresponding to the propagation of the linear perturbation studied are discussed for various flow geometries used.Comment: 13 pages. 5 figure

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented

A Green's function formulation for a nonlinear potential flow solution applicable to transonic flow

Routine determination of inviscid subsonic flow fields about wing-body-tail configurations employing a Green's function approach for numerical solution of the perturbation velocity potential equation is successfully extended into the high subsonic subcritical flow regime and into the shock-free supersonic flow regime. A modified Green's function formulation, valid throughout a range of Mach numbers including transonic, that takes an explicit accounting of the intrinsic nonlinearity in the parent governing partial differential equations is developed. Some considerations pertinent to flow field predictions in the transonic flow regime are discussed

Numerical calculation of transonic boattail flow

A viscid-inviscid interaction procedure for the calculation of subsonic and transonic flow over a boattail was developed. This method couples a finite-difference inviscid analysis with an integral boundary-layer technique. Results indicate that the effect of the boundary layer is as important as an accurate inviscid method for this type of flow. Theoretical results from the solution of the full transonic-potential equation, including boundary layer effects, agree well with the experimental pressure distribution for a boattail. Use of the small disturbance transonic potential equation yielded results that did not agree well with the experimental results even when boundary-layer effects were included in the calculations