Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental data”

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.

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

At least 1,207 records · Page 67

The NASA Langley laminar-flow-control experiment on a swept, supercritical airfoil - Drag equations

The Langley Research Center has designed a swept, supercritical airfoil incorporating Laminar Flow Control for testing at transonic speeds. Analytical expressions have been developed and an evaluation made of the experimental section drag, composed of suction drag and wake drag, using theoretical design information and experimental data. The analysis shows that, although the sweep-induced boundary-layer crossflow influence on the wake drag is too large to be ignored and there is not a practical method for evaluating these crossflow effects on the experimental wake data, the conventional unswept 2-D wake-drag computation used in the reduction of the experimental data is at worst 10 percent too high.

Brooks, Cuyler W., Jr.↗

Fundamental Study of Material Flow in Friction Stir Welds

The presented research project consists of two major parts. First, the material flow in solid-state, friction stir, butt-welds as been investigated using a marker insert technique. Changes in material flow due to welding parameter as well as tool geometry variations have been examined for different materials. The method provides a semi-quantitative, three-dimensional view of the material transport in the welded zone. Second, a FSW process model has been developed. The fully coupled model is based on fluid mechanics; the solid-state material transport during welding is treated as a laminar, viscous flow of a non-Newtonian fluid past a rotating circular cylinder. The heat necessary for the material softening is generated by deformation of the material. As a first step, a two-dimensional model, which contains only the pin of the FSW tool, has been created to test the suitability of the modeling approach and to perform parametric studies of the boundary conditions. The material flow visualization experiments agree very well with the predicted flow field. Accordingly, material within the pin diameter is transported only in the rotation direction around the pin. Due to the simplifying assumptions inherent in the 2-D model, other experimental data such as forces on the pin, torque, and weld energy cannot be directly used for validation. However, the 2-D model predicts the same trends as shown in the experiments. The model also predicts a deviation from the "normal" material flow at certain combinations of welding parameters, suggesting a possible mechanism for the occurrence of some typical FSW defects. The next step has been the development of a three-dimensional process model. The simplified FSW tool has been designed as a flat shoulder rotating on the top of the workpiece and a rotating, cylindrical pin, which extends throughout the total height of the flow domain. The thermal boundary conditions at the tool and at the contact area to the backing plate have been varied to fit experimental data such as temperature profiles, torque and tool forces. General aspects of the experimentally visualized material flow pattern are confirmed by the 3-D model.

Reynolds, Anthony P.↗

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electromagnetic backscattering by corner reflectors

The Geometrical Theory of Diffraction (GTD), which supplements Geometric Optics (GO), and the Physical Theory of Diffraction (PTD), which supplements Physical Optics (PO), are used to predict the backscatter cross sections of dihedral corner reflectors which have right, obtuse, or acute included angles. These theories allow individual backscattering mechanisms of the dihedral corner reflectors to be identified and provide good agreement with experimental results in the azimuthal plane. The advantages and disadvantages of the geometrical and physical theories are discussed in terms of their accuracy, usefulness, and complexity. Numerous comparisons of analytical results with experimental data are presented. While physical optics alone is more accurate and more useful than geometrical optics alone, the combination of geometrical optics and geometrical diffraction seems to out perform physical optics and physical diffraction when compared with experimental data, especially for acute angle dihedral corner reflectors.

Balanis, C. A.↗

Theoretical line shapes for rotational spectra of HCl in Ar

An interaction potential for HCl-Ar recently derived from spectra of van der Waals complexes (Hutson,1988) has been used to obtain converged close-coupling scattering S matrices from which pressure broadening (linewidth, line shift, and line coupling) cross sections have been computed within the Fano-Ben-Reuven formalism. Linewidths agree well with experimental data with the exception of the 0-1 line at low temperatures (an error of about 15 percent at 125 K). Line shifts, especially the largest, for the 0-1 line, are also in reasonable accord with experimental data. Line coupling is predicted to be quite small. The coupled-states molecular scattering approximation is found to be accurate except at very low collision energies. The infinite-order sudden approximations, on the other hand, is not reliable for this system. Inelastic scattering out of the spectroscopic levels accounts for only a fraction of the linewidth at thermal energies, especially in the lower rotational levels.

Green, Sheldon↗

Towards a methodology for robust parameter identification

Consideration is given to the problem of estimating, from experimental data, real parameters for a model with uncertainty in the form of both additive noise and norm-bounded perturbations. Such models frequently arise in robust control theory, and a framework is introduced for the consideration of experimental data in robust control analysis problems. If the analysis tools applied include robust stability tests for real parameter variations (real mu), the framework can be used to address the problem of robust parameter identification. While the techniques discussed can quickly become computationally overwhelming when applied to physical systems and real data, the approach introduces a novel way of looking at the identification problem and may be helpful in arriving at a more tractable methodology.

Smith, Roy S.↗

Progress in Validation of Wind-US for Ramjet/Scramjet Combustion

Validation of the Wind-US flow solver against two sets of experimental data involving high-speed combustion is attempted. First, the well-known Burrows- Kurkov supersonic hydrogen-air combustion test case is simulated, and the sensitively of ignition location and combustion performance to key parameters is explored. Second, a numerical model is developed for simulation of an X-43B candidate, full-scale, JP-7-fueled, internal flowpath operating in ramjet mode. Numerical results using an ethylene-air chemical kinetics model are directly compared against previously existing pressure-distribution data along the entire flowpath, obtained in direct-connect testing conducted at NASA Langley Research Center. Comparison to derived quantities such as burn efficiency and thermal throat location are also made. Reasonable to excellent agreement with experimental data is demonstrated for key parameters in both simulation efforts. Additional Wind-US feature needed to improve simulation efforts are described herein, including maintaining stagnation conditions at inflow boundaries for multi-species flow. An open issue regarding the sensitivity of isolator unstart to key model parameters is briefly discussed.

Engblom, William A.↗

Measurements and Theoretical Calculations of N2-broadening and N2-shift Coefficients in the v2 band of CH3D

In this paper, we report measured Lorentz N2-broadening and N2-induced pressure-shift coefficients of CH3D in the v2 fundamental band using a multispectrum fitting technique. These measurements were made by analyzing 11 laboratory absorption spectra recorded at 0.0056 cm(exp -1) resolution using the McMath-Pierce Fourier transform spectrometer located at the National Solar Observatory on Kitt Peak, Arizona. The spectra were obtained using two absorption cells with path lengths of 10.2 and 25 cm. The total sample pressures ranged from 0.98 to 402.25 Torr with CH3D volume mixing ratios of 0.01 in nitrogen. We have been able to determine the N2 pressure- broadening coefficients of 368 v2 transitions with quantum numbers as high as J"= 20 and K = 16, where K" = K' equivalent to K (for a parallel band). The measured N2-broadening coefficients range from 0.0248 to 0.0742 cm(exp -1) atm(exp -1) at 296 K. All the measured pressure-shifts are negative. The reported N2-induced pressure-shift coefficients vary from about 0.0003 to 0.0094 cm(exp -1) atm(exp -1). We have examined the dependence of the measured broadening and shift parameters on the J", and K quantum numbers and also developed empirical expressions to describe the broadening coefficients in terms of m (m = -J", J", and J" + 1 in the (sup Q)P-, (sup Q)Q-, and (sup Q)R-branch, respectively) and K. On average, the empirical expressions reproduce the measured broadening coefficients to within 4.7%. The N2-broadening and pressureshift coefficients were calculated on the basis of a semiclassical model of interacting linear molecules performed by considering in addition to the electrostatic contributions the atom atom Lennard-Jones potential. The theoretical results of the broadening coefficients are in good overall agreement with the experimental data (8.7%). The N2-pressure shifts whose vibrational contribution is derived from parameters fitted in the (sup Q)Q-branch of self-induced shifts of CH3D, are also in reasonable agreement with the scattered experimental data (20% in most cases).

Predoi-Cross, A.↗

Simulating Underexpanded Jets in Martian and Lunar Environments

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The GCD PSI Project, planned to be accomplished over a four year period, contains a Predictive Simulation Capability (PSC) Element focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the PSC Element, the PSI Project also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper describes the activities and accomplishments in the past year for the Prediction of Plume Flow in low pressure environments component of the PSI Project. While the Loci/Chem Computational fluid dynamics (CFD) application has been validated and used extensively for simulating launch environments in atmospheric conditions, use of this tool for simulating supersonic plumes at Mars-like ambient pressure requires further validation. CFD simulations of underexpanded jets have been performed using Loci/Chem for Mars-like conditions to predict several metrics for steady laminar, turbulent, and impinging plumes. The CFD results are compared with an experimental data set for low Reynolds number plumes at these conditions in order to evaluate the current capability of the Loci/Chem tool for these types of environments. The CFD validation results to date show reasonable agreement with the experimental data across all of the metrics of interest for the configurations considered. The Loci/Chem-Boltzmann CFD application is a hybrid continuum/rarefied flow solver which extends modeling capabilities to very low pressure environments such as those on the Moon. The Loci/Chem-Boltzmann solver uses a gradient-based continuum breakdown criterion to restrict solution of the computationally expensive Boltzmann equation to only a subset of the domain, while using the Navier-Stokes equations elsewhere. This CFD application is under active development, and current efforts toward establishing production ready capability for evaluating Lunar plume surface interactions are well under way. Initial simulations of an Apollo LEM indicate that several regions of the flow require solving the Boltzmann equations due to extreme rarefaction. Early simulations are promising, indicating reasonable overall computational time for a full 3D human scale lander simulation.

Thomas Shurtz↗

Identification of Best Practices for Predicting Inlet Performance Using FUN3D Part 2: Installed Inlets.

A series of studies were performed to assess the impacts of boundary condition type and placement, grid refinement, and modeling parameters such as turbulence model and flux limiter on the predicted inlet performance for installed inlet configurations using the FUN3D flow solver. Two configurations were considered for the study; a wall-mounted Boundary Layer Ingestion (BLI) inlet and the C607 propulsion model tested in the 8x6 Supersonic Wind Tunnel at the NASA Glenn Research Center. The results of the studies were to be used to recommend best practices, as well as to assess the accuracy of FUN3D for inlet predictions. The results of BLI inlet stud-ies showed a minimal impact of grid refinement on the predicted inlet performance for a constant mass flow rate through the inlet. For the C607 propulsion model, the results showed that while the FUN3D predictions at the Aerodynamic Inter-face Plane (AIP) qualitatively agree with the experimental data, FUN3D showed a tendency to overpredict the circumferential distortion metric (IDCmax) and un-derpredict both the radial distortion metric (IDRmax) and the pressure recovery at the AIP (PRAIP ), with the differences between FUN3D and the experimental data increasing with increasing grid refinement and Mach number. Additionally, the outflow boundary location studies performed for both geometries showed that the solution at the AIP was not significantly impacted by the outflow boundary location as long as it was not placed at the location of the AIP. The modeling parameter studies did not indicate a path forward for improved predictions for either inlet con-figuration. Finally, comparisons between the mass flow plug and outflow geometry versions of the C607 propulsion model illustrated favorable agreement, which indi-cates that the differences observed are not caused by the outflow boundary model for this problem. This problem poses significant challenges to Reynolds-averaged Navier-Stokes (RANS) solvers due to the presence of shocks and flow separation in the inlet.

FUN3D↗

Simulating Underexpanded Plumes in Martian and Lunar Environments

While the Loci/Chem Computational fluid dynamics (CFD) application has been validated and used extensively for simulating launch environments in atmospheric conditions, use of this tool for simulating supersonic plumes at much lower ambient pressure requires further validation. Simulations of underexpanded plumes have been performed using Loci/Chem at Martian pressure to predict several metrics for laminar, turbulent, and impinging plumes. The CFD results are compared with an experimental data set for low Reynolds number plumes at these conditions in order to evaluate the current capability of the Loci/Chem tool for these types of environments. The CFD validation results to date show reasonable agreement with the experimental data across all the metrics of interest for the configurations considered. The Loci/Chem-Boltzmann tool is a hybrid continuum/rarefied flow solver which extends modeling capabilities to very low-pressure environments such as those on the Moon. The Loci/Chem-Boltzmann solver uses a gradient-based continuum breakdown criterion to restrict solution of the computationally expensive Boltzmann equation to only a subset of the domain, while using the Navier-Stokes equations elsewhere. This CFD application is under active development, and current efforts toward establishing production ready capability for evaluating Lunar plume surface interactions are well under way. Initial simulations of an Apollo Lunar excursion module indicate that several regions of the flow require solving the Boltzmann equations due to extreme rarefaction. Early simulations are promising, indicating reasonable overall computational time for a full 3D human scale lander simulation.

Plume Surface Interaction↗

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Baseline Validation of Unstructured Grid Reynolds-Averaged Navier-Stokes Toward Flow Control

The value of the use of the Reynolds-averaged Navier-Stokes methodology for active flow control applications is assessed. An experimental flow control database exists for a NACA0015 airfoil modified at the leading edge to implement a fluidic actuator; hence, this configuration is used. Computational results are documented for the baseline wing configuration (no control) with the experimental results and assumes two-dimensional flow. The baseline wing configuration has discontinuities at the leading edge, trailing edge, and aft of midchord on the upper surface. A limited number of active flow control applications have been tested in the laboratory and in flight. These applications include dynamic stall control using a deformable leading edge, separation control for takeoff and landing flight conditions using piezoelectric devices, pulsed vortex generators, zero-net-mass oscillations, and thrust vectoring with zero-net-mass piezoelectric-driven oscillatory actuation. As yet, there is no definitive comparison with experimental data that indicates current computational capabilities can quantitatively predict the large aerodynamic performance gains achieved with active flow control in the laboratory. However, one study using the Reynolds-averaged Navier-Stokes (RANS) methodology has shown good quantitative agreement with experimental results for an isolated zero-net-mass actuator. In addition, some recent studies have used RANS to demonstrate qualitative performance gains compared with the experimental data for separation control on an airfoil. Those quantitative comparisons for both baseline and flow control cases indicated that computational results were in poor quantitative agreement with the experiments. The current research thrust will investigate the potential use of an unstructured grid RANS approach to predict aerodynamic performance for active flow control applications building on the early studies. First the computational results must quantitatively match experiments for the no-control case before proceeding to the time-dependent flow control case. This paper documents the baseline (no-control) case using an unswept airfoil configuration. The next section describes the configurations used for the computations and the experimentals. The computational approach is then described followed by results and concluding remarks.

Joslin, Ronald D.↗

The Paradox of a Wet (High H2O) and Dry (Low H2O/Ce) Mantle: High Water Concentrations in Mantle Garnet Pyroxenites from Hawaii

Water dissolved as trace amounts in anhydrous minerals has a large influence on the melting behavior and physical properties of the mantle. The water concentration of the oceanic mantle is inferred from the analyses of Mid-Ocean Ridge Basalt (MORB) and Oceanic Island Basalt (OIB). but there is little data from actual mantle samples. Moreover, enriched mineralogies (pyroxenites, eclogites) are thought as important sources of heterogeneity in the mantle, but their water concentrations and their effect on the water budget and cycling in the mantle are virtually unknown. Here, we analyzed by FTIR water in garnet clinopyroxenite xenoliths from Salt Lake Crater, Oahu, Hawaii. These pyroxenites are high-pressure (>20kb) crystal fractionates from alkalic melts. The clinopyroxenes (cpx) have 260 to 576 ppm wt H2O, with the least differentiated samples (Mg#>0.8) in the 400-500 ppm range. Orthopyroxene (opx) contain 117-265 ppm H2O, about half of that of cpx, consistent with other natural sample studies, but lower than cpx/opx equilibrium from experimental data. The pyroxenite cpx and opx H2O concentrations are at the high-end of on-and off-craton peridotite xenolith concentrations and those of Hawaiian spinel peridotites. In contrast, garnet has extremely low water contents (<5ppm H2O). There is no correlation between H2O in cpx and lithophile element concentrations. Phlogopite is present in some samples, and its modal abundance shows a positive correlation in Mg# with cpx, implying equilibrium. However, there is no correlation between H2O concentrations and or the presence of phlogopite. These data imply that cpx and opx may be at water saturation, far lower than experimental data suggest. Reconstructed bulk rock pyroxenite H2O ranges from 200-460 ppm (average 331 +/- 75 ppm), 2 to 8 times higher than H2O estimates for the MORB source (50-200 ppm), but in the range of E-MORB, OIB and the source of rejuvenated Hawaiian magmas. The average bulk rock pyroxenite H2O/Ce is 69 +/-35, lower than estimates of the MORB source (approx 150) or FOZO, C (200-250) mantle component, but consistent with "dry" EM sources (<100). These data suggest that a metasomatized, refertilized oceanic lithosphere that contains pyroxenitic veins (e.g. the lower part of an oceanic plate, where ascending melts can become trapped and crystallize), will have both higher water concentrations and low H2O/Ce, and may contribute to EM-type OIB sources, like that of Samoa basalts. Therefore, a low H2O/Ce mantle source may not necessarily be "dry".

Peslier, Anne H.↗

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Neural Networks for Supercapacitor Modeling

This software implements a one-dimensional model of a supercapacitor and calibrates it to experimental data. A system of coupled partial differential equation models the salt concentration and the ionic and electronic potentials. The term in the equation modeling the capacitance contains a neural network, whose parameters are optimized to fit the provided experimental data. The software is written in python using the open-source library JAX to calculate the gradients needed for the optimization.

Salazar de Troya, Miquel↗

Three-dimensional nearfield characterization of a VSTOL jet in turbulent crossflow

This paper documents the measurement and prediction of the near-field evolution of the steady, subsonic, turbulent three-dimensional external aerodynamic flow-field associated with a high speed jet issued from an orifice perpendicular to an imposed crossflow velocity field. A parabolic Navier-Stokes CFD solver is employed, utilizing a Kappa-epsilon closure system with algebraic Reynolds stress equation, in concert with suitably defined porous far-field boundary conditions and a virtual source initial condition procedure. The CFD results correlate well with NASA generated near-field LDA experimental data on mean velocity and rms fluctuation velocity (interpolated to turbulent kinetic energy) distributions. The numerical results in the farther field predict the vortex roll-up and jet deflection/entrainment characterization verified by pitot-static mean velocity experimental data. Detailed results are presented for a circular jet for velocity ratios ranging 4 to 8.

Baker, A. J.↗