SEARCH · Engineering Papers

Results for “Lifting map”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1971

Observational consequences of Martian wind regimes.

The connection with past and future observations of Mars of the high velocity relief winds deduced by Gierasch and Sagan (1971) is examined with the assistance of a large topographic map of Mars. A unimodal hypsometric curve is derived. Seasons and locales of wind velocities at the half surface pressure level below 80 m/sec, sufficient to lift dust at the surface, are identified. A number of suggestions are made that can be checked by observations with the Mariner Mars 1971 orbiters.

Sagan, C.↗

NASA NTRS1990

Aerothermodynamic measurements on a proposed assured crew return vehicle (ACRV) lifting-body configuration at Mach 6 and 10 in air

A 0.02-scale model of a lifting-body concept for possible application to the Assured Crew Return Vehicle from Space Station Freedom was tested at Mach 6 and 10 in air. Thermal mappings and surface streamline patterns were obtained at angles of attack ranging from 0 to 30 deg and unit Reynolds numbers Re from 2 to 8 x 10 to the 6th/ft. Areas that experienced the highest heating were near the model nose and tip-fin leading edges. The effect of Re on windward centerline heating coefficients was negligible, whereas increases in angles of attack produced increases in heating. At Mach 6 and the highest unit Re, turbulent heat at the windward centerline was three to four times the laminar level. Leeward crossflow separation and vortex reattachment along the centerline are evident across the Re and angle-of-attack ranges tested, indicative of a complex flowfield.

Horvath, Thomas J.↗

NASA NTRS1992

Finite wing aerodynamics with simulated glaze ice

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional wing is studied experimentally. The model used for these tests was a semi-span wing of effective aspect ratio five, mounted from the sidewall of the UIUC subsonic wind tunnel. The model has an NACA 0012 airfoil section on a rectangular, untwisted planform with interchangeable leading edges to allow for testing both the baseline and the iced wing geometry. A three-component sidewall balance was used to measure lift, drag and pitching moment on the clean and iced model. A four-beam two-color fiberoptic laser Doppler velocimeter (LDV) was used to map the flowfield along several spanwise cuts on the model. Preliminary results from LDV scans, which will be the bulk of this paper, are presented following the force balance measurement results. Initial comparison of LDV surveys compare favorably with inviscid theory results and 2D split hot-film measurements near the model surface.

Khodadoust, A.↗

NASA NTRS2015

Description, Usage, and Validation of the MVL-15 Modified Vortex Lattice Analysis Capability

MVL-15 is the most recent version of the Modified Vortex-Lattice (MVL) code developed within the Aerodynamics Systems Analysis Branch (ASAB) at NASA LaRC. The term "modified" refers to the primary modification of the core vortex-lattice methodology: inclusion of viscous aerodynamics tables that are linked to the linear solution via iterative processes. The inclusion of the viscous aerodynamics inherently converts the MVL-15 from a purely analytic linearized method to a semi-empirical blend which retains the rapid execution speed of the linearized method while empirically characterizing the section aerodynamics at all spanwise lattice points. The modification provides a means to assess non-linear effects on lift that occur at angles of attack near stall, and provides a means to determine the drag associated with the application of design strategies for lift augmentation such as the use of flaps or blowing. The MVL-15 code is applicable to the analyses of aircraft aerodynamics during cruise, but it is most advantageously applied to the analysis of aircraft operating in various high-lift configurations. The MVL methodology has been previously conceived and implemented; the initial concept version was delivered to the ASAB in 2001 (van Dam, C.), subsequently revised (Gelhausen, P. and Ozoroski, T. 2002 / AVID Inc., Gelhausen, P., and Roberts, M. 2004), and then overhauled (Ozoroski, T., Hahn, A. 2008). The latest version, MVL-15 has been refined to provide analysis transparency and enhanced to meet the analysis requirements of the Environmentally Responsible Aviation (ERA) Project. Each revision has been implemented with reasonable success. Separate applications of the methodology are in use, including a similar in-house capability, developed by Olson, E. that is tailored for structural and acoustics analyses. A central premise of the methodology is that viscous aerodynamic data can be associated with analytic inviscid aerodynamic results at each spanwise wing section, thereby providing a pathway to map viscous data to the inviscid results. However, a number of factors can sidetrack the analysis consistency during various stages of this process. For example, it should be expected that the final airplane lift curve and drag polar results depend strongly on the geometry and aerodynamics of the airfoil section; however, flap deflections and flap chord extensions change the local reference geometry of the input airfoil, the airplane wing, the tabulated non-dimensional viscous aerodynamics, and the spanwise links between the linear and the viscous aerodynamics. These changes also affect the bound circulation and therefore, calculation and integration of the induced angle of attack and induced drag. MVL-15 is configured to ensure these types of challenges are properly addressed. This report is a comprehensive manual describing the theory, use, and validation of the MVL-15 analysis tool. Section 3 summarizes theoretical, procedural, and characteristic features of MVL-15, and includes a list of the files required to setup, execute, and summarize an analysis. Section 4, Section 5, Section 6, and Section 7 combine to comprise the User's Guide portions of this report. The MVL-15 input and output files are described in Section 4 and Section 5, respectively; the descriptions are supplemented with example files and information about the file formats, parameter definitions, and typical parameter values. Section 6 describes the Wing Geometry Setup Utility and the 2d-Variants Utility files that simplify and assist setting up a consistent set of MVL-15 geometry and aerodynamics input parameters and input files. Section 7 describes the use of the 3d-Results Presentation Utility file that can be used to automatically create summary tables and charts from the MVL-15 output files. Section 8 documents the Validation Results of an extensive and varied validation test matrix, including results of an airplane analysis representative of the ERA Program. A start-to-finish example of the airplane analysis procedure is described in Section 7.

Ozoroski, Thomas A.↗

NASA NTRS1970

Calculation of incompressible fluid flow through cambered blades

Conformal mapping technique yields linear, approximate solutions for calculating flow of an incompressible fluid through staggered array of cambered blades for the cases of flow with partial cavitation and supercavitation. Lift and drag coefficients, cavitation number, cavity shape, and exit flow conditions can be determined.

Hsu, C. C.↗

NASA NTRS2010

Orion Entry Performance-Based Center-of-Gravity Box

The Orion capsule is designed both for Low Earth Orbit missions to the ISS and for missions to the moon. For ISS class missions, the capsule will use an Apollo-style direct entry. For lunar return missions, depending on the timing of the mission, the capsule could perform a direct entry or a skip entry of up to 4800 n.mi. in order to land in the coastal waters of California. The physics of atmospheric re-entry determine the capability of the Orion vehicle. For a given vehicle mass and shape, physics tells us that the driving parameters for an entry vehicle are the hypersonic lift-to-drag ratio (L/D) and the flight path angle at entry interface (gamma(sub EI)). The design of the Orion atmospheric re-entry must meet constraints during both nominal and dispersed flight conditions on landing accuracy, heating rate, total heat load, sensed acceleration, and proper disposal of the Service Module. These constraints define an entry corridor in the space of L/D-gamma(sub EI); if the vehicle falls within this corridor, then all constraints are met. The gamma(sub EI) dimension of the corridor can be further constrained by the gloads experienced during emergency entries. Thus, the entry performance for the Orion vehicle can be described completely by the L/D. Bounds on the hypersonic L/D necessary to achieve all the mission requirements can be defined for the given entry corridor. Landing accuracy performance drives the lower limit on L/D. In order to achieve the desired landing accuracy, a minimum L/D must be ensured. The design of the Thermal Protection System (TPS) drives the upper limit on L/D. A higher L/D can drive mass into the design of the TPS. Conversely, once the TPS is designed, the L/D must be ensured to stay below a certain limit in order for the TPS to stay within its design envelop. The L/D must stay within its upper and lower bounds during dispersed flight conditions. L/D is a function of both the aerodynamics and the center-of-gravity (CG) of the vehicle. The aerodynamics of the vehicle are determined by Computational Fluid Mechanics (CFD) and wind tunnel tests. However, the aerodynamics are not known precisely. Instead, an aerodynamic database has been developed where the aerodynamic coefficients are known to fall within a probabilistic band defined by upper and lower bounds. It is expected that the probabilistic band will shrink after the first missions are flown and real-world data is collected. Until that time, the Orion must be designed to the current aerodynamic database. Thus, for a given aerodynamic database with given uncertainties, the allowable range in L/D can be mapped to an allowable box for the CG location. The CG box is used to set requirements on the dispersions allowed for vehicle packaging and cargo storage. As the aerodynamic uncertainties decrease, the size of the CG box can increase. This paper discusses the technique used to map the minimum and maximum L/D bounds set by the entry performance requirements to the allowable dispersions in CG while accounting for aerodynamic uncertainties. The L/D is defined as the ratio of the lift force to the drag force. It is equivalent to the ratio of lift coefficient (C(sub L)) over drag coefficient (C(sub D)). C(sub L) and C(sub D) are functions of Mach number (M) and angle of attack (alpha). A Mach number of 25 is used as a measuring point of the hypersonic L/D. Variations in C(sub L), C(sub D) and alpha cause variations in L/D. Equation (1) shows the three contributions to the variation in L/D.

Rea, Jeremy R.↗

NASA NTRS1979

A summary of research on the NASA-Global Atmospheric Sampling Program performed by the Atmospheric Sciences Research Center

The annual variations of ozone near the tropopause are derived from aircraft exhibit year-to-year differences which are not explicitly accounted for by the simple, classical ozone transport theory. Phenomena such as tropopause lifting, interannual variations in the rates of stratospheric-tropospheric exchange and meridional mixing, contribute differently to the distribution of ozone in this altitude region. Ozone encounter climatologies have been represented by global maps which show the probabilities of exceeding ambient ozone levels of 200, 300, and 400 ppbV along flight routes during the year. Continuous ozone records obtained from the GASP system revealed the presence of gravity waves whose wavelength is of the order 20 km. The GASP data cannot, however, be utilized for the evaluation of horizontal fluxes of such quantities as ozone, sensible heat, and zonal momentum; the data are too sparsely and irregularly distributed for the computation of stable correlations. Multiple species data from the unique circumglobal flight of a Pan American airliner on 28-30 October 1977 are discussed with particular regard to the apparent interhemispheric differences in tropospheric species concentrations, variation between the Arctic and Antarctic stratospheres, to possible covariations between species, and to potential source regions for various constituents.

Falconer, P. D.↗

NASA NTRS1999

Coronal Dimmings and Energetic CMEs

We have analyzed the coronal dimmings for seven fast (> 600 km/s) coronal mass ejections (CMEs) occurring between 23 April and 9 May which were associated with flares from NOAA active region (AR) 8210. These dimming regions were identified by their strong depletion in coronal emission within a half hour of the estimated time of CME lift-off. They included areas which were as dark as quiescent coronal holes as well as other regions with weaker brightness depletions. We found that the extended dimming areas in these events generally mapped out the apparent "footprint" of the CME. In two of the seven cases, a pair of dimmings were more or less symmetrically positioned north and south of the flare site. In the five remaining cases, the dimmings were most prominent to the north of AR 8210 (approximately S15 latitude) and extended well north of the solar equator, consistent with the locations of the CMEs. We discuss the implications of these results for the sigmoid/double dimming/flux rope model of CMEs.

Thompson, B. J.↗

DOE OSTI2024

A Research Testbed for Experimental Connectivity and Automation in Cars

Cyber-physical systems (CPS) are engineered systems in which the computational components and physical components of the system are interdependent and deeply linked. Discovery and development of CPS technologies are changing the way people interact with engineered systems. To demonstrate this point, consider the Global Positioning System (GPS) which revolutionized personal mobility. Its research origins stem from the US Department of Defense in 1973 for military use. By 2000 limitations on civilian consumption of GPS service were lifted. Immediately, personal GPS products were shipped for in-car navigation systems. Soon after, GPS receivers were embedded into mobile phones with the introduction of the Nokia N95 and Apple iPhone for maps and turn-by-turn navigation. With the benefit of hindsight we can see that GPS has become a cornerstone technology of our daily life. Satellite-based navigation has gone from a radical space technology, to a matter of uninteresting availability; it is the typical way that people navigate through roadways now. The combination of the satellite infrastructure, the integration of receiver modules, and powerful maps applications has significantly changed mobility in the 21st century.

42 ENGINEERING↗

NASA NTRS1983

Design analysis of Rayleigh-step floating-ring seals

The analysis and design of a 50-mm diameter floating-ring helium buffer seal are described. The seal rings incorporated Rayleigh-step lift pads to provide hydrodynamic forces to separate the rings from the shaft. Maximum surface speed is 183 m/s (600 fps) and maximum buffer gas pressure is 1389 kPa (200 psia). An operating range map was computed as a function of speed and pressure. Contradictory problems arise due to excessive friction preventing ring tracking at low-speed, high-pressure conditions and insufficent friction to retard inertia driven motions at high-speed, low-pressure conditions. Steady-state and dynamic analyses and performance are described, as well as the results of thermal studies.

Artiles, A.↗

DOE OSTI2023

Sulfur Lone Pairs Control Topology in Heterotrimetallic Complexes: An Experimental and Theoretical Study

Heterotrimetallic complexes with (N 2 S 2 )M metallodithiolates, M = Ni 2+ , [Fe(NO)] 2+ , and [Co(NO)] 2+ , as bidentate chelating ligands to a central trans-Cr(NO)(MeCN) unit were characterized as the first members of a new class, NiCrNi, FeCrFe, CoCrCo. The complexes exhibit a cisoid structural topology, ascribed to the stereoactivity of the available lone pair(s) on the sulfur donors, resulting in a dispersed, electropositive pocket from the N/N and N/S hydrocarbon linkers wherein the Cr-NO site is housed. Computational studies explored alternative isomers (transoid and inverted cisoid) that suggest a combination of electronic and steric effects govern the geometrical selectivity. Electrostatic potential maps readily display the dominant electronegative potential from the sulfurs which force the NO to the electropositive pocket. The available S lone pairs work in synergy with the π-withdrawing ability of NO to lift Cr out of the S 4 plane toward the NO and stabilize the geometry. The metallodithiolate ligands bound to Cr(NO) thus find structural consistency across the three congeners. Although the dinitrosyl [(bme-dach)Co(NO)-Mo(NO)(MeCN)-(bme-dach)Co(MeCN)][PF 6 ] 2 (CoMoCo') analogue displays chemical noninnocence and a partial Mo–Co bond toward (N 2 S 2 )Co'(NCCH 3 ) in an “asymmetric butterfly” topology, the stability of the {Cr(NO)} 5 unit prohibits such bond rearrangement. Magnetism and EPR studies illustrate spin coupling across the sulfur thiolate sulfur bridges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

Estimate of loads during wing-vortex interactions by Munk's transverse-flow method

The inviscid, incompressible interaction of a wing with a vortex is studied by use of Munk's transverse-flow method. The method assumes that the loading on the wing is such that the local circulatory flow of the vortex is turned so that the wing and its vortex wake act as a barrier to the flow. This permits the analysis to be carried out by mapping the transverse flowfield into the flow about a circle to find the vorticity distribution in the wake. Closed-form expressions are then derived for the bound circulation in the wind and for the lift and rolling moment induced by the vortex on the encountering wing. Comparisons of the loads predicted by these relationships with those of vortex-lattice theory for a flat wing of a rectangular planform indicate that they accurately represent the various parameters when the aspect ratio of the encountering wing is less than about two. When flat rectangular wings of higher aspect ratios are considered, some sort of correction is needed. Examples are then presented to illustrate some applications of the results.

Rossow, Vernon J.↗

DOE OSTI2024

Liquid Interfacial Electron Microscopy Identifies Nanogalvanic Corrosion in Pearlitic Steel

The nanoscale mechanisms of localized corrosion in low carbon steels have remained elusive due to the complexity of studying the degradative material behavior at nanoscale solid-liquid interfaces. We identified various steps in the nanogalvanic corrosion processes using in-situ liquid-cell scanning transmission electron microscopy (STEM) using a microfluidic holder by Hummingbird Scientific. Initial work, performed at low magnification, identified the initiation point on a 1018 low-carbon steel surface. This initiation point was determined to be a triple junction of two ferrite grains bridging a cementite grain in contact with a baseline electrolyte of 6 uM CO2 dissolved in a buffered (2.78 uM Na2SO4) aqueous solution, pH 6.1. The pre-etched low-carbon steel surface was prepared using focused ion beam lift-out procedures to extract a cross-section of the low-carbon steel surface, which then was thinned to about 150 nm and transferred to a SiN membrane microfluidic window. The transfer was made using a lift-out needle to attach the low-carbon steel lamella to the corner of the SiN window, and then Pt/C deposition held the lamella in contact with the window while it was released from the lift out needle. To identify the triple point on the low carbon steel lamella, prior to attachment on the SiN window, the sample was characterized for compositional variations with energy dispersive x-ray spectroscopy mapping, grain orientation and phase mapping with precession electron diffraction, and thickness mapping with energy filtered transmission electron microscopy. This pre-characterization prior to the in-situ experiment provided a map of the multiphase and multigrain structure, where the in-situ liquid cell imaging provided a clear understanding of the initiation point on the sample. These data were cross-correlated to paint a holistic picture of the triple junction site, enabling low electron-fluence in-situ snapshot imaging to avoid dominating the native corrosion reactions with effects from the incident electron beam. This initial result identified that localized, nanogalvanic corrosion at the phase interface was the dominant corrosion process in the low-carbon steel, so we next targeted the observation of an array of these nanogalvanic features phase boundaries in a pearlite grain. Near-surface ferrite/cementite phase interfaces that typify pearlitic low-carbon steel were extracted, pre-characterized, and imaged for the in-situ corrosion processes. The sample was a cross-section from a pearlite grain, with alternating ferrite and cementite grains that extended microns down from the pre-etched low-carbon steel pipe surface. After contact with a buffered aqueous solution, the phase boundaries between the ferrite and cementite began to dissolve, with observable material loss and thickness changes in the dark-field and bright-field STEM images. Within minutes, the corrosion front proceeded deeper into the material, claiming a thin layer of ferrite around all exposed phase boundaries before progressing laterally into the ferrite matrix, converting the ferrite to corrosion product normal to each buried cementite grain. Formation of the corrosion product causes a volumetric expansion, creating a lateral wedging force that mechanically ejects the cementite grains from their grooves and leaves behind percolation channels into the steel substructure. Rapid and deleterious, this nanogalvanic corrosion pathway represents an important target for understanding and preventing run-away degradation in this common building material. Observation of this corrosion mechanism was enabled by the combination of pre-characterization using standard structural, grain, and compositional analysis in the TEM, which provides maps for understanding the reaction propagation captured in low-dose, in-situ, liquid-cell STEM.

corrosion↗

DOE OSTI2026

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation) [SWR-26-095]

FOILPOLARS (Grassmannian Foil Shape Sweeps for Polar Generation): Multifidelity aerodynamic polar data generation for hydrofoil/tidal-turbine airfoil sections. Foilpolars ties together three pieces: *AeroSandbox supplies the baseline airfoil coordinates (UIUC database). *G2Aero parameterizes those shapes on a Grassmannian manifold (Karcher mean + PGA basis) and samples new perturbed shapes around that basis. *XFoil (panel method) and NeuralFoil (neural-network surrogate, shipped with AeroSandbox) each solve the resulting shapes for lift, drag, moment, and pressure at the swept angles of attack, Reynolds numbers, and n_crit values. Design optimization of foil shapes in a computationally efficient way requires polars data across many candidate shapes, not just a handful of baseline foils. However, high-fidelity CFD at that scale is too costly, and naive shape perturbation strays from realistic geometries. FOILPOLARS addresses this by loading baseline airfoils (via AeroSandbox) and mapping them onto a Grassmannian manifold (via G2Aero), computing a Karcher mean and principal geodesic analysis (PGA) basis. New shapes are sampled by perturbing PGA coefficients, keeping them close to the manifold of realistic foils. Each sampled shape is evaluated across a configurable sweep of angle of attack, Reynolds number, and critical amplification factor using two solvers: XFoil (panel method) and NeuralFoil (neural-network surrogate), producing a paired dataset of lift, drag, moment, pressure, convergence, and confidence, indexed alongside each shape's PGA coefficients and shared Grassmannian basis in a single xarray dataset. From this, FOILPOLARS produces convergence summaries and comparison plots per shape, Reynolds number, and n_crit. A command-line interface exposes each pipeline stage independently, supporting data-driven design, optimization, and machine-learning workflows for foils.

Sandhu, Rimple [National Laboratory of the Rockies↗

NASA NTRS1973

A compilation of studies from Atmospheric Variability Experiment (AVE)

Five methods for obtaining the sign of vertical motion at various levels in the atmosphere were investigated to determine which gave the best explanation for areas of rain and no-rain in the southeastern United States during the period February 19-22, 1964. The methods used were the terrain-induced vertical motion, the kinematic method including the terrain effect, the adiabatic method, the omega equation, and the vorticity equation combined with Dines' Compensation Principle. Stability and moisture availability were considered but not as extensively as vertical motion. Values of vertical motion obtained by the kinetic method, including orographic lifting produced the best agreement with the observed areas. When terms in the omega equation were added through use of constant multipliers, results comparable to the adiabatic method were obtained. Without this addition large uncertainties occurred. Maps showing areas where terms of the omega equation would indicate positive vertical motion did not correlate well with the observed rainfall patterns.

Scoggins, J. R.↗

NASA NTRS1986

Analytical observations on the aerodynamics of a delta wing with leading edge flaps

The effect of a leading edge flap on the aerodynamics of a low aspect ratio delta wing is studied analytically. The separated flow field about the wing is represented by a simple vortex model composed of a conical straight vortex sheet and a concentrated vortex. The analysis is carried out in the cross flow plane by mapping the wing trace, by means of the Schwarz-Christoffel transformation into the real axis of the transformed plane. Particular attention is given to the influence of the angle of attack and flap deflection angle on lift and drag forces. Both lift and drag decrease with flap deflection, while the lift-to-drag ratioe increases. A simple coordinate transformation is used to obtain a closed form expression for the lift-to-drag ratio as a function of flap deflection. The main effect of leading edge flap deflection is a partial suppression of the separated flow on the leeside of the wing. Qualitative comparison with experiments is presented, showing agreement in the general trends.

Oh, S.↗

DOE OSTI2022

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1972

Numerical calculation of transonic flow about swept wings.

Description of a mixed elliptic-hyperbolic relaxation algorithm which calculates solutions to the three-dimensional, nonlinear transonic small disturbance potential equation for flows about thin swept lifting wings with free-stream Mach number less than 1. The algorithm is designed to treat supercritical flows, including cases with embedded shock waves. Nonrectangular planform shapes, including sweep and taper, are treated by a coordinate transformation which maps the wing planform into a rectangle. Computed results at angles of attack of 0 and 2 deg for a AR = 4, constant chord, 23.75 deg sweptback planform model with a Lockheed C141 airfoil section are compared with data obtained experimentally for both subcritical and supercritical flows. Subcritical results are also compared with those obtained by a subsonic 'panel' method.

Ballhaus, W. F.↗