Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls

An Extension of Essentially Non-Oscillatory Shock-Capturing Schemes to Multi-Dimensional Systems of Conservation Laws

In recent years, a class of numerical schemes for solving hyperbolic partial differential equations has been developed which generalizes the first-order method of Godunov to arbitrary order of accuracy. High-order accuracy is obtained, wherever the solution is smooth, by an essentially non-oscillatory (ENO) piecewise polynomial reconstruction procedure, which yields high-order pointwise information from the cell averages of the solution at a given point in time. When applied to piecewise smooth initial data, this reconstruction enables a flux computation that provides a time update of the solution which is of high-order accuracy, wherever the function is smooth, and avoids a Gibbs phenomenon at discontinuities. The promising results of Harten et al., in the use of ENO schemes in solving the one-dimensional Euler equations of gas dynamics, have aroused considerable interest in the aerodynamic community. However, the application of these schemes to areas of scientific and industrial interest, such as aircraft configuration, obviously requires compressible flow solutions in more than one spatial dimension. It is this extension of ENO schemes to multi-dimensional application to which this study is dedicated. In particular, a two-dimensional extension is proposed for the Euler equations of gas dynamics. Among the issues to be considered in this extension are achieving formal high-order two-dimensional spatial accuracy, the implementation of boundary conditions and applications to general curvilinear coordinates

An Investigation of High-Order Shock-Capturing Methods for Computational Aeroacoustics

Topics covered include: Low-dispersion scheme for nonlinear acoustic waves in nonuniform flow; Computation of acoustic scattering by a low-dispersion scheme; Algorithmic extension of low-dispersion scheme and modeling effects for acoustic wave simulation; The accuracy of shock capturing in two spatial dimensions; Using high-order methods on lower-order geometries; and Computational considerations for the simulation of discontinuous flows

A New Time Domain Formulation for Broadband Noise Predictions

A new analytic result in acoustics called "Formulation 1B," proposed by Farassat, is used to compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the loading source term. The formulation contains a far field surface integral that depends on the time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The unsteady surface pressure for this test case is analytically specied from a result based on linear airfoil theory. This test case is used to examine the velocity scaling properties of Formulation 1B and to demonstrate its equivalence to Formulation 1A of Farassat. The new acoustic formulation, again with an analytic surface pressure, is then used to predict broadband noise radiated from an airfoil immersed in homogeneous, isotropic turbulence. The results are compared with experimental data previously reported by Paterson and Amiet. Good agreement between predictions and measurements is obtained. Finally, an alternative form of Formulation 1B is described for statistical analysis of broadband noise

Broadband Noise Prediction When Turbulence Simulation Is Available - Derivation of Formulation 2B and Its Statistical Analysis

We show that a simple modification of Formulation 1 of Farassat results in a new analytic expression that is highly suitable for broadband noise prediction when extensive turbulence simulation is available. This result satisfies all the stringent requirements, such as permitting the use of the exact geometry and kinematics of the moving body, that we have set as our goal in the derivation of useful acoustic formulas for the prediction of rotating blade and airframe noise. We also derive a simple analytic expression for the autocorrelation of the acoustic pressure that is valid in the near and far fields. Our analysis is based on the time integral of the acoustic pressure that can easily be obtained at any resolution for any observer time interval and digitally analyzed for broadband noise prediction. We have named this result as Formulation 2B of Farassat. One significant consequence of Formulation 2B is the derivation of the acoustic velocity potential for the thickness and loading terms of the Ffowcs Williams-Hawkings (FW-H) equation. This will greatly enhance the usefulness of the Fast Scattering Code (FSC) by providing a high fidelity boundary condition input for scattering predictions

Computational considerations for the simulation of shock-induced sound

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order-accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order-accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks in a quasi-one-dimensional nozzle flow. The design order of accuracy is achieved in the smooth regions of a steady-state test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order-accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustics simulations of flows with shocks that are raised by these results are discussed

Permeable Surface Corrections for Ffowcs Williams and Hawkings Integrals

The acoustic prediction methodology discussed herein applies an acoustic analogy to calculate the sound generated by sources in an aerodynamic simulation. Sound is propagated from the computed flow field by integrating the Ffowcs Williams and Hawkings equation on a suitable control surface. Previous research suggests that, for some applications, the integration surface must be placed away from the solid surface to incorporate source contributions from within the flow volume. As such, the fluid mechanisms in the input flow field that contribute to the far-field noise are accounted for by their mathematical projection as a distribution of source terms on a permeable surface. The passage of nonacoustic disturbances through such an integration surface can result in significant error in an acoustic calculation. A correction for the error is derived in the frequency domain using a frozen gust assumption. The correction is found to work reasonably well in several test cases where the error is a small fraction of the actual radiated noise. However, satisfactory agreement has not been obtained between noise predictions using the solution from a three-dimensional, detached-eddy simulation of flow over a cylinder

A comparison of two formulations for high-order accurate essentially non-oscillatory schemes

The finite-volume and finite-difference implementations of high-order accurate essentially non-oscillatory shock-capturing schemes are discussed and compared. Results obtained with fourth-order accurate algorithms based on both formulations are examined for accuracy, sensitivity to grid irregularities, resolution of waves that are oblique to the mesh, and computational efficiency. Some algorithm modifications that may be required for a given application are suggested. Conclusions that pertain to the relative merits of both formulations are drawn, and some circumstances for which each might be useful are noted

Airframe Noise Prediction by Acoustic Analogy: Revisited

The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill s acoustic analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry-independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic literature would suggest some interest in the use of tailored or exact Green s function (EGF) for aerodynamic noise problems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill s equation with either the Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of development. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and use of these numerical Green functions are then discussed

HTLV-III Serology in Hemophilia: Relationship with Immunologic Abnormalities

We investigated the relationship of the presence of antibodies to HTLV-III and immunologic abnormalities in patients with hemophilia. Serum antibodies to HTLV-III were analyzed by ELISA assay, immunoprecipitation of labeled cell extracts, and immunoprecipitation of purified HTLV-III p24. Thirty-four (61%) of the total group (n = 56) had antibody to HTLV-III; 34 (76%) of 45 patients given commercial factor VIII preparations were seropositive, compared with none of 11 patients treated exclusively with cryoprecipitate obtained from volunteer blood donors. Of patients who were seropositive for HTLV-III antibody, 94% had abnormal T4/T8 ratios, and 33% of those whose serum was antibody negative had abnormal T4/T8 ratios; five patients, each antibody positive, have lymphadenopathy syndrome. Sequential studies in a subset of patients indicate that there is a changing pattern of antibody production to HTLV-III antigens after seroconversion