Engineering topics

Davis, Sanford

Publications and source records attributed to Davis, Sanford.

17 records

NASA NTRS2002

Vorticity Transport and Wave Emission In A Protoplanetary Disk

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the disk to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived from the nonlinear simulation data and correlated with analytical dispersion relations from the linearized Euler and energy equations.

Davis, S. S.↗

NASA NTRS1999

Mid-Field Sonic Boom Extrapolation Methodology

In the design cycle of low-boom airplanes, sonic boom prediction must be accurate and efficient. The classical linear method, Whitham's F-function theory, has been widely applied to predict sonic boom signatures. However, linear theory fails to capture the nonlinear effects created by large civil transport. Computational fluid dynamics (CFD) has been used successfully to predict sonic boom signals at the near and mid fields. Nevertheless, it is computationally expansive in airplane design runs. In the present study, the method of characteristics is used to predict sonic boom signals in an efficient fashion. The governing equations are the axisymmetric Euler's equations with constant enthalpy. Since the method solves Euler's equations, it captures more nonlinear effects than the classical Whitham's F-function technique. Furthermore, the method of characteristics is an efficient marching scheme for initial value problems. In this study, we will first review the current CFD extrapolation technique and the work previously done in sonic boom extrapolation. Then, we will introduce the governing equations and the method of characteristics. Finally, we will show that the present method yields the same accurate results as previous CFD techniques, but with higher efficiency.

Cheung, Samson↗

NASA NTRS1998

Experimental and Computational Boundary-Layer Studies In a Supersonic Two-Dimensional Nozzle

Experiments and computations were carried out on the adiabatic laminar boundary layer developing along the surfaces of a two-dimensional supersonic nozzle, consisting of upper and lower contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes were employed. These codes were adapted to the characteristics of a specific wind tunnel nozzle, so that their numerical results could be directly compared with experimental data obtained in the same nozzle. Such comparisons were made for the boundary-layer growth on the nozzle contoured surfaces, and for the boundary-layer growth, surface streamlines and surface shear on the sidewalls. The three-dimensional Navier-Stokes code was found to be the only one to correctly predict the mean laminar boundary-layer flow on both the sidewalls and the contoured surfaces. Theory and experiment both indicated that the sidewall flow is highly three-dimensional, with non-uniform shear, corner vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Demetriades, Anthony↗

NASA NTRS1997

Quiet Flow Characteristics of the NASA-Ames Laminar Flow Supersonic Wind Tunnel for Mach 1.6 Transition Studies

Low-disturbance or 'quiet' wind tunnels are an essential part of any meaningful boundary layer transition research. In particular, the receptivity of supersonic boundary layers to wind tunnel disturbances can significantly alter the transition phenomena under investigation on a test model. Consequently, considerable study has gone into the design of a new wind tunnel to provide quiet (low-disturbance) flow, encompassing both theoretical and experimental efforts. Our pilot (eighth-scale) supersonic wind tunnel was reported in 1992. NASA-Ames Fluid Mechanics Laboratory (FML) commissioned a quiet wind tunnel in 1994 to support Supersonic Laminar Flow Control (SLFC) research. Known as the Laminar Flow Supersonic Wind Tunnel (LFSWT), this tunnel is designed to operate at potential cruise Mach numbers and unit Reynolds numbers (Re) of the High Speed Civil Transport (HSCT). The need to better understand the transition phenomena on the leading edge region of swept (HSCT) wings provided the impetus for building the LFSWT. Additional information is contained in the original extended abstract.

Wolf, Stephen W. D.↗

NASA NTRS1997

Boundary Layer Receptivity to Three-Dimensional Freestream Disturbances at Two-Dimensional Roughness

The receptivity of a laminar boundary layer to an isolated three-dimensional convected disturbance is investigated in a low-speed wind tunnel experiment. The disturbance is created by the short-duration pulsed displacement of a small low-aspect-ratio wing located upstream of a flat plate. The height of the wing is set so that the convected disturbance grazes the edge of the flat-plate boundary layer. A receptivity site is provided by a two-dimensional roughness strip on the surface of the plate. The different propagation speeds of acoustic, convected and instability waves cause the various wave packets from the pulsed displacement to arrive at a downstream measurement station at different times, separating the phenomena and allowing them to be studied independently. Ensemble- averaged measurements are made with and without roughness on the plate. Preliminary analysis of the measurements suggest the presence of a two-dimensional T-S wave packet arising from an interaction between an acoustic wave and the roughness, and a three-dimensional T-S wave packet arising from an interaction between the localized convected disturbance and the roughness strip. The growth rates and spatial characteristics of the disturbances and the instability wave packets are measured as they propagate downstream.

Dietz, Anthony↗

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 NTRS1994

NASA Ames Fluid Mechanics Laboratory research briefs

The Ames Fluid Mechanics Laboratory research program is presented in a series of research briefs. Nineteen projects covering aeronautical fluid mechanics and related areas are discussed and augmented with the publication and presentation output of the Branch for the period 1990-1993.

Davis, Sanford↗

NASA NTRS1993

Finite-Difference Algorithms For Computing Sound Waves

Governing equations considered as matrix system. Method variant of method described in "Scheme for Finite-Difference Computations of Waves" (ARC-12970). Present method begins with matrix-vector formulation of fundamental equations, involving first-order partial derivatives of primitive variables with respect to space and time. Particular matrix formulation places time and spatial coordinates on equal footing, so governing equations considered as matrix system and treated as unit. Spatial and temporal discretizations not treated separately as in other finite-difference methods, instead treated together by linking spatial-grid interval and time step via common scale factor related to speed of sound.

Davis, Sanford↗

NASA NTRS1993

An algorithm for solving the perturbed gas dynamic equations

The present application of a compact, higher-order central-difference approximation to the linearized Euler equations illustrates the multimodal character of these equations by means of computations for acoustic, vortical, and entropy waves. Such dissipationless central-difference methods are shown to propagate waves exhibiting excellent phase and amplitude resolution on the basis of relatively large time-steps; they can be applied to wave problems governed by systems of first-order partial differential equations.

Davis, Sanford↗

NASA NTRS1992

Scheme For Finite-Difference Computations Of Waves

Compact algorithms generating and solving finite-difference approximations of partial differential equations for propagation of waves obtained by new method. Based on concept of discrete dispersion relation. Used in wave propagation to relate frequency to wavelength and is key measure of wave fidelity.

Davis, Sanford↗

NASA NTRS1991

Low-dispersion finite difference methods for acoustic waves in a pipe

A new algorithm for computing one-dimensional acoustic waves in a pipe is demonstrated by solving the acoustic equations as an initial-boundary-value problem. Conventional dissipation-free second-order finite difference methods suffer severe phase distortion for grids with less that about ten mesh points per wavelength. Using the signal generation by a piston in a duct as an example, transient acoustic computations are presented using a new compact three-point algorithm which allows about 60 percent fewer mesh points per wavelength. Both pulse and harmonic excitation are considered. Coupling of the acoustic signal with the pipe resonant modes is shown to generate a complex transient wave with rich harmonic content.

Davis, Sanford↗

NASA NTRS1991

Computation Of Facilitated Transport of O2 In Hemoglobin

Report describes computations of unsteady facilitated transport of oxygen through liquid membrane of hemoglobin. Used here, "facilitated transport" means diffusion of permeant through membrane in which that diffusion enhanced by reversible chemical reaction between permeant and membrane. In this case, reversible reactions between hemoglobin and oxygen.

Davis, Sanford↗

NASA NTRS1991

Time-accurate finite difference solutions to the incompressible Navier-Stokes/energy equations

Two new algorithms for solving coupled Navier-Stokes and energy equations are presented. The algorithms are compared with available Navier-Stokes solutions for the model problem of a driven-cavity-flow over a range of Reynolds numbers. It is noted that the algorithms represent two different implementations of the fractional step approach to the solution of the coupled Navier-Stokes and energy equations using the Boussinesq approximation.

Mirfakhraee, Ali↗

NASA NTRS1990

Matrix-based, finite-difference algorithms for computational acoustics

A compact numerical algorithm is introduced for simulating multidimensional acoustic waves. The algorithm is expressed in terms of a set of matrix coefficients on a three-point spatial grid that approximates the acoustic wave equation with a discretization error of O(h exp 5). The method is based on tracking a local phase variable and its implementation suggests a convenient coordinate splitting along with natural intermediate boundary conditions. Results are presented for oblique plane waves and compared with other procedures. Preliminary computations of acoustic diffraction are also considered.

Davis, Sanford↗

NASA NTRS1990

Computation of the unsteady facilitated transport of oxygen in hemoglobin

The transport of a reacting permeant diffusing through a thin membrane is extended to more realistic dissociation models. A new nonlinear analysis of the reaction-diffusion equations, using implicit finite-difference methods and direct block solvers, is used to study the limits of linearized and equilibrium theories. Computed curves of molecular oxygen permeating through hemoglobin solution are used to illustrate higher-order reaction models, the effect of concentration boundary layers at the membrane interfaces, and the transient buildup of oxygen flux.

Davis, Sanford↗

NASA NTRS1989

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

NASA NTRS1987

Flow physics

A discussion is presented on the coupling of computational analysis and experiment. It is believed that this coupling is critical in developing new aerodynamic insights. Additionally, new methods for analyzing and interpreting data are discussed. These methods need to be developed in small-scale research studies and then applied to large-scale technology programs. The specific objectives of this program are threefold: (1) provide definitive data sets for the assessment of numerical simulations to the Navier-Stokes equations; (2) incorporate advanced instrumentation to measure the spatial and temporal structure of fluid flows; and (3) develop true parallelism between computational and experimental research using the 'scientific workstation' concept. The discussion is presented in viewgraph form.

Davis, Sanford↗