Engineering topics

Ghia, U.

Publications and source records attributed to Ghia, U..

At least 19 records

NASA NTRS1996

Development of a Vorticity-Velocity Navier-Stokes Formulation for the Study of Compressibility Effects on Dynamic Stall

The first major area of this study was to develop a vorticity-velocity formulation and numerical solution algorithms suitable for the analyses of incompressible as well as low-to- moderate-speed compressible flows. Research performed towards contributing to the determination of the appropriate vorticity and dilation creation boundary conditions suggested to temporarily set aside this approach and use a primitive-variable approach other than the pseudo-compressibility approach used. The second major area of study was initiated to comprehensively examine the INS-2D and INS-3D programs from the point of view of boundary conditions. The research carried out was documented in the form of two technical papers which are included in Appendices A and B; the boundary-condition related issues for INS-3D are briefly mentioned.

Ghia, K. N.↗

NASA NTRS1992

Characterization of dynamic stall phenomenon using two-dimensional unsteady Navier-Stokes equations

Among the new significant aspects of the present work are: (1) the treatment of the far-field boundary; (2) the use of C-grid topology, with the branch-cut singularity treated analytically; (3) evaluation of the effect of the envelope of prevailing initial states, and finally; (4) the ability to employ streakline/pathline 'visualization' to probe the unsteady features prevailing in vortex-dominated flows. The far-field boundary is placed at infinity, using appropriate grid stretching. This contributes to the accuracy of the solutions, but raised a number of important issues which needed to be resolved; this includes determining the equivalent time-dependent circulation for the pitching airfoil. A secondary counter-clockwise vortex erupts from within the boundary layer and immediately pinches off the energetic leading-edge shear layer which then, through hydrodynamic instability, rolls up into the dynamic stall vortex. The streakline/pathline visualization serves to provide information for insight into the physics of the unsteady separated flow.

Ghia, K. N.↗

NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗

NASA NTRS1991

Examination of a vortex-ring interaction phenomenon in an axisymmetric flow

Unsteady separated flow in an axisymmetric configuration is investigated, with the objective of analyzing vortex-ring formation and break-up and the ensuing interactions believed to be leading to breakdown of laminar flow. The model problem considered is an idealized representation of a combustor with a centerbody. Direct solution of the time-dependent incompressible Navier-Stokes equations is employed, with central-differencing for all spatial derivatives, to determine the response of this flow for Re = 2000, 5000 and 10,000. The last case revealed a rather novel unsteady vortex-ring interaction phenomenon. This case has been further investigated using third-order accurate upwind differencing for the convective terms. The flow is interrogated carefully by examining the corresponding instantaneous vorticity, stream-function and velocity fields, as well as by tracking mass-less marker particles in some of the critical regions. Such detailed examination is deemed essential, prior to assigning fundamental physical relevance to the observed vortex-interaction phenomenon.

Ghia, U.↗

NASA NTRS1991

A semi-elliptic analysis for 2-D viscous flows through cascade configurations

A semielliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is negligible but which are significantly influenced by upstream interactions. A two-step alternating-direction-explicit numerical scheme is developed to solve the resulting governing equations efficiently. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Various simple channel and cascade flow configurations are examined for several values of Re, thickness ratio and Mach number, in order to establish the applicability of the IPNS model for these flows. With this model, solutions to compressible as well as nearly incompressible flows are obtained without any modification either in the analysis or in the solution procedure.

Ramamurti, R.↗

NASA NTRS1990

Low-dimensional modelling and prediction in fluid dynamics

Research into the development of a simplified, low-dimensional, state-variable model for a fluid dynamics system is discussed. The model would be capable of predicting the state of the system in as nearly real time as possible and would be suitable for implementation in a control loop. It is noted that state-variable modeling has an advantage for real-time health monitoring and control application over the Navier-Stokes equations since the latter require an enormous amount of real time. In this study, a dynamic procedure is developed for modeling and prediction using the CFD data generated for this model centerbody-combustor configuration. Some experiences with the NASA data available for the SSME are also reported. Application of this procedure to pressure and velocity measurements in an actual wind-tunnel experiment is described, and the applicability of the approach to health monitoring is considered.

Ghia, U.↗

NASA NTRS1990

Analysis and control of low-speed forced unsteady flow

A capability for numerically simulating 2-D flows in temporally deforming geometries is described, with emphasis on flow with forced unsteadiness, particularly on the simulation and analysis of these flows. The simulation of forced unsteady flows makes the examination of fundamental unsteady flow mechanisms, such as dynamic stall and unsteady separation, possible. A turbulence model is being incorporated into the analysis so as to obtain solutions for the higher Reynolds numbers used in the experiments. The analysis is also of utility in studying fluid-structure interactions, free surfaces, metal-forming, and bio-fluid mechanics involving flow through passages with flexible walls.

Ghia, U.↗

NASA NTRS1989

Analysis and control of unsteady separated flows

The active control of unsteady separated flows is considered, and the development of the separated flow following the formation and shedding of the starting vortex is characterized. In the method, the Navier-Stokes equations are formulated in time-dependent general orthogonal coordinates which remain boundary aligned for arbitrary boundary deformations. Results obtained for a backstep-channel flow with an oscillating flap demonstrate the effectiveness of the flap as a control mechanism.

Ghia, U.↗

NASA NTRS1989

Direct method for solution of three-dimensional unsteady incompressible Navier-Stokes equations

A three-dimensional, unsteady analysis for the velocity-vorticity formulation of the Navier-Stokes (NS) equations is described. This formulation highlights the natural separation of the spin dynamics of a fluid particle from its translational kinematics. This separation is exploited to produce a single pass direct inversion of the full unsteady NS operator without the need to iterate either the future wall vorticity or the nonlinear vorticity convection terms. The method has uniform second-order accuracy in time and space. As an example, it is applied to the problem of the three-dimensional flow within a shear-driven cubical box. Results are obtained for Re = 100, 400, and 1000 using a preliminary (25, 25, 25) mesh.

Osswald, G. A.↗

NASA NTRS1988

Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids

A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils.

Ramamurti, R.↗

NASA NTRS1988

Solution of the Neumann pressure problem in general orthogonal coordinates using the multigrid technique

The multigrid (MG) technique has been advanced for use with Neumann boundary-value problems in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of incompressible flow using the velocity-pressure formulation of the Navier-Stokes equations. The finite-difference representation of the problem and the formulation of the restriction and coarse-grid correction operators are examined in detail. Maintaining consistency between these and the integral constraint associated with the Neumann problem is found to be critical for the success of the MG technique. The influence of the smoothing operator is examined by employing Gauss-Seidel, alternating-direction implicit, and strongly implicit techniques. The MG procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 14, depending on the type of smoothing operator employed and the values of the problem parameters.

Ghia, U.↗

NASA NTRS1987

A direct algorithm for solution of incompressible three-dimensional unsteady Navier-Stokes equations

A direct, implicit, numerical solution algorithm, with second-order accuracy in space and time, is constructed for the three-dimensional unsteady incompressible Navier-Stokes equations formulated in terms of velocity and vorticity, using generalized orthogonal coordinates to achieve the accurate solution of complex viscous flow configurations. A numerically stable, efficient, direct inversion procedure is developed for the computationally intensive divergence-curl elliptic velocity problem. This overdetermined partial differential operator is first formulated as a uniquely determined, nonsingular matrix-vector problem; this aspect of the procedure is a unique feature of the present analysis. The three-dimensional vorticity-transport equation is solved by a modified factorization technique which completely eliminates the need for any block-matrix inversions and only scalar tridiagonal matrices need to be inverted. The method is applied to the test problem of the three-dimensional flow within a shear-driven cubical box. Coherent streamwise vortex structures are observed within the steady-state flow at Re = 100.

Osswald, G. A.↗

NASA NTRS1987

Analysis of unsteady wake of a circular cylinder using Navier-Stokes equations

An unsteady Navier-Stokes analysis of the separated flow in the wake of a circular cylinder is presented in order to study the wake structure and the various interactions between its scales at low Re. The quasi-periodic and chaotic flows noted by Sreenivasan (1985) are searched for. The analysis is used to characterize the secondary vortex formation in the early stages of the flow development and its interaction with the primary vortex.

Ghia, K. N.↗

NASA NTRS1987

Simulation of separated flow past a bluff body using Navier-Stokes equations

Two-dimensional flow past a bluff body is presently simulated on the basis of an analysis that employs the incompressible, unsteady Navier-Stokes equations in terms of vorticity and stream function. The fully implicit, time-marching, alternating-direction, implicit-block Gaussian elimination used is a direct method with second-order spatial accuracy; this allows it to avoid the introduction of any artificial viscosity. Attention is given to the simulation of flow past a circular cylinder with and without symmetry, requiring the use of either the half or the full cylinder, respectively.

Ghia, K. N.↗

NASA NTRS1986

Simulation of self-induced unsteady motion in the near wake of a Joukowski airfoil

The unsteady Navier-Stokes analysis is shown to be capable of analyzing the massively separated, persistently unsteady flow in the post-stall regime of a Joukowski airfoil for an angle of attack as high as 53 degrees. The analysis has provided the detailed flow structure, showing the complex vortex interaction for this configuration. The aerodynamic coefficients for lift, drag, and moment were calculated. So far only the spatial structure of the vortex interaction was computed. It is now important to potentially use the large-scale vortex interactions, an additional energy source, to improve the aerodynamic performance.

Ghia, K. N.↗

NASA NTRS1986

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg., (2) Re =1000, alpha = 15 deg., (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗

NASA NTRS1985

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg, (2) Re = 1000, alpha = 15 deg, (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗

NASA NTRS1984

Multi-grid technique for the solution of incompressible Navier-Stokes equations

A computationally efficient method is proposed for obtaining fine-mesh solutions of the Navier-Stokes equations for high-Re flow, including separation. The method involves implementing a multi-grid solution procedure with suitably chosen elements, such that the resulting overall solution technique is efficient as well as robust. The robustness and efficiency of the solution technique are demonstrated by applying it to three model problems: flow in a driven cavity, downstream asymptotic flow in curved ducts of square and polar sections, and Newmann boundary-value problem in clustered curvilinear orthogonal coordinates.

Ghia, K. N.↗