Engineering Papers⌕ Search

Engineering topics

Fung, K.-Y.

Publications and source records attributed to Fung, K.-Y..

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A compact solution to computational acoustics

This paper demonstrates that the linearized, dimensional Euler equations for acoustic computation can be accurately solved as a set of decoupled first-order wave equations, and that if ordered properly, this system of simple waves has unambiguous, easily implemented boundary conditions, allowing waves of same group speeds to pass through numerical boundaries or comply with wall conditions. Thus, the task of designing a complex multi-dimensional scheme with approximate far-field boundary conditions reduces to the design of higher order schemes for the one-dimensional simple wave equation. A compact finite-difference scheme and a characteristically exact but numerically n(th) order accurate boundary condition are introduced for solving the first order wave equation. Spanning a three-point two-level stencil, this low-dispersion implicit scheme has a third order spatial accuracy when used on nonuniform meshes, fourth order accurate on uniform meshes, and a temporal accuracy of second order due to the choice of trapezoidal integration for algorithmic simplicity. The robustness and accuracy of the scheme are demonstrated through a series of numerical experiments and comparisons with published results. When tested on the one-dimensional wave equation on a uniform grid, this scheme allows a Gaussian wave packet to pass through any finite domain with low numerical dispersion characteristic of a spatially fourth-order scheme and reflections at numerical boundaries maintained below truncation error. On highly stretched and irregular grids, only mild dispersions are found in the solution while solutions by other methods fail or are severely distorted. Yet, this scheme is no more sophisticated to solve or implement than the Crank-Nicolson scheme. This scheme has been tested on four categories of the ICASE/LaRC benchmark problems, which include propagation of acoustic and convective waves in Cartesian and cylindrical domains, reflection of acoustic wave at stationary/moving boundaries, and sound generation by gust-blade interaction.

Fung, K.-Y.↗
NASA NTRS1992

Effects Of Compressibility On Dynamic Stall

Report presents results of computations and measurements of compressible flow about an airfoil, angle of attack of which oscillates about static-stall angle. Study focuses on effects of compressibility on dynamic stall. Of particular interest are conditions determining onset of separation of flow, which leads to premature dynamic stall and consequent significant reduction of lift.

Carr, L. W.↗
NASA NTRS1992

Computed unsteady flows of airfoils at high incidence

The flow over an airfoil at an angle of attack above the static stall angle would ordinarily be massively separated. Under dynamic conditions, the onset of stall can be delayed to an angle, depending on the type of unsteadiness, much higher than that for static stall. The stall onset mechanisms under dynamic conditions are unclear. Due to extreme difficulties involved, experimental investigations, so far, have provided insufficient information about the flow field for the identification of the onset mechanisms. A course of classical boundary layer analysis augmented with numerical experiments and measured data is chosen here instead, with the hope that the identification of onset mechanisms can be more systematic and quantitative.

Fung, K.-Y.↗
NASA NTRS1989

Viscous-inviscid interaction and local grid refinement via truncation error injection

A methodology is presented which makes it possible to decouple a complex problem having multiple disparate length scales into problems of single length scale so that they can be solved more efficiently on a computer. The method is applied to a viscous transonic flow over an airfoil. It is found that accurate prediction of the flow over an airfoil can be obtained by solving the Euler equations on a relatively coarse global grid with viscous effects computed separately on a boundary-layer type grid and injected into the global grid solution as a combination of vorticity and trucation error.

Goble, Brian D.↗
NASA NTRS1988

The effects of compressibility on dynamic stall

In this paper, typical computational predictions and experimental measurements of compressible flow past an airfoil at dynamic stall conditions are studied and compared to develop an insight into the effect of compressibility on dynamic stall. The dependency of the critical Mach number on airfoil leading edge curvature, camber, and angle of attack is investigated. Evidence is presented to show that a local region of supersonic flow occurs on an oscillating airfoil, even for a freestream Mach number as low as 0.2, if the boundary layer remains attached and the angle of attack is sufficiently high; that a shock terminates this local supersonic bubble; and that the vorticity that this shock generates grows rapidly and becomes very unstable as the angle of attack increases beyond the value at which the maximum local flow speed first exceeds the speed of sound. It is suggested that these shock-induced effects compete with the dynamic viscous effects occurring in the boundary layer in determining the onset of separation, which can lead to premature dynamic stall and can significantly reduce the maximum dynamic lift that can otherwise be obtained.

Fung, K.-Y.↗
NASA NTRS1987

A truncation error injection approach to viscous-inviscid interaction

An approach to viscous-inviscid interaction which is based on truncation error injection is presented in the context of solving flow over an airfoil. A two-dimensional interpolation scheme is used to restrict the fine grid solutions to the global coarse grid. Details on the current implementation of the approach are given, and the boundary conditions being used are discussed. Inviscid results from a NACA0012 airfoil test case and the viscous results are presented.

Goble, B. D.↗
NASA NTRS1985

Refined numerical solution of the transonic flow past a wedge

A numerical procedure combining the ideas of solving a modified difference equation and of adaptive mesh refinement is introduced. The numerical solution on a fixed grid is improved by using better approximations of the truncation error computed from local subdomain grid refinements. This technique is used to obtain refined solutions of steady, inviscid, transonic flow past a wedge. The effects of truncation error on the pressure distribution, wave drag, sonic line, and shock position are investigated. By comparing the pressure drag on the wedge and wave drag due to the shocks, a supersonic-to-supersonic shock originating from the wedge shoulder is confirmed.

Liang, S.-M.↗
NASA NTRS1985

Computation of unsteady transonic aerodynamics with steady state fixed by truncation error injection

A novel technique is introduced for efficient computations of unsteady transonic aerodynamics. The steady flow corresponding to body shape is maintained by truncation error injection while the perturbed unsteady flows corresponding to unsteady body motions are being computed. This allows the use of different grids comparable to the characteristic length scales of the steady and unsteady flows and, hence, allows efficient computation of the unsteady perturbations. An example of typical unsteady computation of flow over a supercritical airfoil shows that substantial savings in computation time and storage without loss of solution accuracy can easily be achieved. This technique is easy to apply and requires very few changes to existing codes.

Fung, K.-Y.↗
NASA NTRS1985

A simple, accurate and efficient algorithm for unsteady transonic flow

The derivation of an algorithm for computing an airfoil's response to small unsteady perturbations is described. The analysis of the far-field boundary conditions for the flow field, and the application of the time-linearization technique to develop an accurate description of the flow field are discussed. The treatment of moving shock waves is examined. The alternating direction implicit scheme of Ballhaus and Steger (1975) is utilized for the discretization of the governing equations and boundary conditions. Numerical examples demonstrating the applicability of the time-linearization method to determine the correlation between airfoil and pressure distribution are provided. The effect of wind tunnel walls on unsteady transonic testing is studied.

Fung, K.-Y.↗
NASA NTRS1982

A model for unsteady transonic indicial responses

A method of characterizing the indicial curve for the response frequencies of flutter in the transonic regime by means of a time scale parameter is illustrated. Time linearization of the unsteady potential equation is assumed adequate and the perturbed potential about a steady flow is defined with appropriate boundary conditions. An indicial motion can then be obtained for a given mode, with the lift and drag coefficients which result from a step change in incidence or control surface deflection available to describe any kind of motion. The lift and drag coefficients are shown to result analytically by addition of the multiplicand involving the time scale parameter. If the coefficients are maximized, then at reduced frequencies only two parameters are necessary to form the indicial response curve.

Fung, K.-Y.↗
NASA NTRS1978

A new method for designing shock-free transonic configurations

A new method for the design of shock-free supercritical airfoils, wings, and three-dimensional configurations is described. Results illustrating this procedure in two and three dimensions are given. They include modifications to part of the upper surface of an NACA 64A410 airfoil that will maintain shock-free flow over a range of Mach numbers for a fixed lift coefficient, and the modifications required on part of the upper surface of a swept wing with an NACA 64A410 root section to achieve shock-free flow. While the results are given for inviscid flow, the same procedures can be employed iteratively with a boundary layer calculation in order to achieve shock-free viscous designs. With a shock-free pressure field the boundary layer calculation will be reliable and not complicated by the difficulties of shock-wave boundary-layer interaction.

Sobieczky, H.↗
NASA NTRS1977

Unsteady transonic flow computations

The numerical procedures previously developed for computing nonlinear and time-linearized small-perturbation unsteady transonic flows are briefly reviewed, and the effects of unsteady modes of motion on two-dimensional transonic flows are evaluated. The numerical procedure used comprises an alternating-direction implicit scheme and treats shock waves as discontinuities in the flow. Comparison of the time-linearized results with fully nonlinear calculations delineates their range of applicability. The unsteady behavior due to harmonic pitching and flap oscillations of an NACA airfoil is also examined.

Seebass, A. R.↗