Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic instability”

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 253 records · Page 14

Solar irradiance changes caused by g-modes and large scale convection

Solar irradiance measurements from the ACRIM experiment show a clear response to the rotation periods of g-mode oscillations (l = 1, 2, and 3) and their first harmonics. Peaks in the ACRIM spectrum at 16.6, 18.3, 20.7, 36.5, and about 71 days all lie within about 1 percent of periods arising from g-mode rotation. This means that the g-modes are a fundamental cause of irradiance fluctuations. On time scales of months and less they modulate the irradiance by means of transient flows of global scale which they stimulate in the sun's convective envelope. Dimensional arguments indicate that the flows carry up heat at an average rate of about 0.001 solar luminosities, which is not in conflict with observed changes in the irradiance. Five additional tests for g-modes and large-scale convection are given. An instability is described which undermines diffusion models of sunspot energy storage.

Wolff, C. L.↗

The dynamics and control of large flexible space structures. Part B: Development of continuum model and computer simulation

The equations of motion of an arbitrary flexible body in orbit were derived. The model includes the effects of gravity with all its higher harmonics. As a specific example, the motion of a long, slender, uniform beam in circular orbit was modelled. The example considers both the inplane and three dimensional motion of the beam in orbit. In the case of planar motion with only flexible vibrations, the pitch motion is not influenced by the elastic motion of the beam. For large values of the square of the ratio of the structural modal frequency to the orbital angular rate the elastic motion was decoupled from the pitch motion. However, for small values of the ratio and small amplitude pitch motion, the elastic motion was governed by a Hill's 3 term equation. Numerical simulation of the equation indicates the possibilities of instability for very low values of the square of the ratio of the modal frequency to the orbit angular rate. Also numerical simulations of the first order nonlinear equations of motion for a long flexible beam in orbit were performed. The effect of varying the initial conditions and the number of modes was demonstrated.

Bainum, P. M.↗

Phase Coupling Between Spectral Components of Collapsing Langmuir Solitons in Solar Type III Radio Bursts

We present the high time resolution observations of one of the Langmuir wave packets obtained in the source region of a solar type III radio burst. This wave packet satisfies the threshold condition of the supersonic modulational instability, as well as the criterion of a collapsing Langmuir soliton, i.e., the spatial scale derived from its peak intensity is less than that derived from its short time scale. The spectrum of t his wave packet contains an intense spectral peak at local electron plasma frequency, f(sub pe) and relatively weaker peaks at 2f(sub pe) and 3f(sub pe). We apply the wavelet based bispectral analysis technique on this wave packet and compute the bicoherence between its spectral components. It is found that the bicoherence exhibits two peaks at (approximately f(sub pe), approximately f(sub pe)) and (approximately f(sub pe) approximately 2f(sub pe)), which strongly suggest that the spectral peak at 2f(sub pe) probably corresponds to the second harmonic radio emission, generated as a result of the merging of antiparallel propagating Langmuir waves trapped in the collapsing Langmuir soliton, and, the spectral peak at 3f(sub pe) probably corresponds to the third harmonic radio emission, generated as a result of merging of a trapped Langmuir wave and a second harmonic electromagnetic wave.

Thejappa, G.↗

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.↗

In Situ Detection of Strong Langmuir Turbulence Processes in Solar Type III Radio Bursts

The high time resolution observations obtained by the WAVES experiment of the STEREO spacecraft in solar type III radio bursts show that Langmuir waves often occur as intense localized wave packets. These wave packets are characterized by short durations of only a few ms and peak intensities, which well exceed the supersonic modulational instability (MI) thresholds. These timescales and peak intensities satisfy the criterion of the solitons collapsed to spatial scales of a few hundred Debye lengths. The spectra of these wave packets consist of primary spectral peaks corresponding to beam-resonant Langmuir waves, two or more sidebands corresponding to down-shifted and up-shifted daughter Langmuir waves, and low frequency enhancements below a few hundred Hz corresponding to daughter ion sound waves. The frequencies and wave numbers of these spectral components satisfy the resonance conditions of the modulational instability (MI). Moreover, the tricoherences, computed using trispectral analysis techniques show that these spectral components are coupled to each other with a high degree of coherency as expected of the MI type of four wave interactions. The high intensities, short scale lengths, sideband spectral structures and low frequency spectral enhancements and, high levels of tricoherences amongst the spectral components of these wave packets provide unambiguous evidence for the supersonic MI and related strong turbulence processes in type III radio bursts. The implication of these observations include: (1) the MI and related strong turbulence processes often occur in type III source regions, (2) the strong turbulence processes probably play very important roles in beam stabilization as well as conversion of Langmuir waves into escaping radiation at the fundamental and second harmonic of the electron plasma frequency, fpe, and (3) the Langmuir collapse probably follows the route of MI in type III radio bursts.

Langmuir↗

The Effects of Acoustic Treatment on Pressure Disturbances From a Supersonic Jet in a Circular Duct

The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined circular duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions at each axial location. These values are used to determine the growth rates and phase velocities of the instability waves and the near field pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, and wall impedance, and the effects of these changes on non-axisymmetric modes. The primary results indicate that the effects of the duct wall on stability characteristics diminish as the outer flow increases and as the jet azimuthal mode number increases. Also, wall reflections are reduced when using a finite impedance boundary condition at the wall; but in addition, reflections are reduced and growth rates diminished by keeping the imaginary part of the impedance negative when using the negative exponential for the harmonic dependence.

Dahl, Milo D.↗

Cyclotron maser and plasma wave growth in magnetic loops

Cyclotron maser and plasma wave growth which results from electrons accelerated in magnetic loops are studied. The evolution of the accelerated electron distribution is determined by solving the kinetic equation including Coulomb collisions and magnetic convergence. It is found that for modest values of the column depth of the loop the growth rates of instabilities are significantly reduced and that the reduction is much larger for the cyclotron modes than for the plasma wave modes. The large decrease in the growth rate with column depth suggests that solar coronal densities must be much lower than commonly accepted in order for the cyclotron maser to operate. The density depletion has to be similar to that which occurs during auroral kilometric radiation events in the magnetosphere. The resulting distributions are much more complicated than the idealized distributions used in many theoretical studies, but the fastest growing mode can still simply be determined by the ratio of electron plasma to gyrofrequency, U=omega(sub p)/Omega(sub e). However, the dominant modes are different than for the idealized situations with growth of the z-mode largest for U approximately less than 0.5, and second harmonic x-mode (s=2) or fundamental o-mode (s=1) the dominant modes for 0.5 approximately less than U approximately less than 1. The electron distributions typically contain more than one inverted feature which could give rise to wave growth. It is shown that this can result in simultaneous amplification of more than one mode with each mode driven by a different feature and can be observed, for example, by differences in the rise times of the right and left circularly polarized components of the associated spike bursts.

Hamilton, Russell J.↗

Amplification of Reynolds number dependent processes by wave distortion

A system using a hot-wire transducer as an analog of a liquid droplet of propellant was employed to investigate the ingredients of the acoustic instability of liquid-propellant rocket engines. It was assumed that the combustion process was vaporization-limited and that the combustion chamber was acoustically similar to a closed-closed right-circular cylinder. Before studying the hot-wire closed-loop system (the analog system), a microphone closed-loop system, which used the response of a microphone as the source of a linear feedback exciting signal, was investigated to establish the characteristics of self-sustenance of acoustic fields. Self-sustained acoustic fields were found to occur only at resonant frequencies of the chamber. In the hot-wire closed-loop system, the response of hot-wire anemometer was used as the source of the feedback exciting signal. The self-sustained acoustic fields which developed in the system were always found to be harmonically distorted and to have as their fundamental frquency a resonant frequency for which there also existed a second resonant frequency which was approximately twice the fundamental frequency.

Ventrice, M. B.↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Source of the Bursty Bulk Flow Diffuse Aurora: Electrostatic Cyclotron Harmonic and Whistler Waves in the Coupling of Bursty Bulk Flows to Auroral Precipitation

Electron cyclotron harmonic (ECH) and whistler chorus waves are recognized as the two mechanisms responsible for the resonant wave‐particle interactions necessary to precipitate plasma sheet electrons into the ionosphere, producing the diffuse Aurora. Previous work has demonstrated ECH waves dominate electron scattering at L shells >8, while whistler chorus dominates scattering at L shells L < 8. However, we find from Time History of Events and Macroscale (THEMIS) Interactions during Substorms observations of fast flows at L = 12 that oblique whistler chorus emissions play the dominant role in scattering electrons. Previous works have identified whistler‐mode waves within fast flows that are produced by an electron temperature anisotropy Te,⊥/Te,||> 1, consistent with electron betatron acceleration. Here, however, we find whistler chorus emissions throughout an interval of fast flows where Te,⊥/Te,||< 1. Parallel electron beams account for the enhanced parallel electron temperature and serve as the instability mechanism for the whistler chorus. The parallel electron beams and associated cigar‐shaped distributions are consistent with Fermi acceleration at dipolarizations in fast flows. We demonstrate that the scattering efficiency of the whistler chorus exceeds that of ECH waves, which THEMIS also detects during the fast flows. The obliquity of the whistler waves permits efficient scattering of lower‐energy electrons into the diffuse aurora. We conclude that Fermi acceleration of electrons provides one important free‐energy source for the wave‐particle interactions responsible for coupling plasma sheet electrons into the diffuse aurora during substorm conditions.

Wendel, D. E.↗

On the unique structure of the magnetic fields of Uranus and Neptune

The magnetic fields of Uranus and Neptune, which have comparable dipole, quadrupole, and octupole harmonics, are unique in the present-day solar system, but they resemble the geomagnetic field at the epochs of excursions and reversals known from paleomagnetic data. The precession dynamo model, in which the dominant role in the generation of the planetary magnetic fields is played by external gravitational forces, allows us to propose two scenarios for the formation of the unique topology of the magnetic fields of Uranus and Neptune. In the first case, tidal flows in the 'oceans' of these two planets extend down to the depths where the matter has a noticeable electric conductivity and velocity. A hydromagnetic interaction of the moving conducting fluid with the planetary magnetic field outside the generation region results in the deformation of the field and the deceleration of the motion under the action of the radial magnetic field. In the second case, the deformation of the field facilitates drastic changes in cyclonic cells within the generation region causing instabilities that result in a multi-polar field structure, excursions, and inversions. This paper considers this problem in greater detail by using the Neptune-Triton system as an example.

Dolginov, Sh. SH.↗

A Regional Nuclear Conflict Would Compromise Global Food Security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 ◦C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30◦N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

Jaegermeyr, Jonas↗

Variability of accretion disks surrounding black holes: The role of inertial-acoustic mode instabilities

The global nonlinear time-dependent evolution of the inertial-acoustic mode instability in accretion disks surrounding black holes has been investigated. The viscous stress is assumed to be proportional to the gas pressure only, i.e., tau = alphap(sub g). It is found that an oscillatory nonsteady behavior exists in the inner regions of disks (r is less than 10r(sub g) where r(sub g) is the Schwarzschild radius) for sufficiently large alpha(greater than or approximately equal to 0.2) and for mass accretion rates less than about 0.3 times the Eddington value. The variations of the integrated bolometric luminosity from the disk, Delta L/L, are less than 3%. A power spectrum analysis of these variations reveals a power spectrum which can be fitted to a power-law function of the frequency Pis proportional to f(exp -gamma), with index gamma = 1.4-2.3 and a low-frequency feature at about 4 Hz in one case. In addition, a narrow peak centered at a frequency corresponding to the maximum epicyclic frequency of the disk at approximately 100-130 Hz and its first harmonic is also seen. The low-frequency modulations are remarkably similar to those observed in black hole candidate systems. The possible existence of a scattering corona in the inner region of the disk and/or other processes contributing to the power at high frequencies in the inner region of the accretion disk may make the detection of the high-frequency component difficult.

Chen, Xingming↗

Harmonic Balance Computations of Fan Aeroelastic Stability

A harmonic balance (HB) aeroelastic analysis, which has been recently developed, was used to determine the aeroelastic stability (flutter) characteristics of an experimental fan. To assess the numerical accuracy of this HB aeroelastic analysis, a time-domain aeroelastic analysis was also used to determine the aeroelastic stability characteristics of the same fan. Both of these three-dimensional analysis codes model the unsteady flowfield due to blade vibrations using the Reynolds-averaged Navier-Stokes (RANS) equations. In the HB analysis, the unsteady flow equations are converted to a HB form and solved using a pseudo-time marching method. In the time-domain analysis, the unsteady flow equations are solved using an implicit time-marching approach. Steady and unsteady computations for two vibration modes were carried out at two rotational speeds: 100 percent (design) and 70 percent (part-speed). The steady and unsteady results obtained from the two analysis methods compare well, thus verifying the recently developed HB aeroelastic analysis. Based on the results, the experimental fan was found to have no aeroelastic instability (flutter) at the conditions examined in this study.

Bakhle, Milind A.↗

Key Observations and Discoveries from the Vector Electric Field Investigation on the C/NOFS Satellite

Launched in 2008 and operating for 7.5 years, the Air Force Communication /Navigation Outage Forecasting System (C/NOFS) satellite included the Vector Electric Field Investigation (VEFI) designed and built at NASA’s Goddard Space Flight Center. VEFI successfully met its objectives: 1) investigate the role of ambient electric fields in initiating nighttime density depletions and turbulence; 2) determine the quasi-DC electric fields associated with abrupt density depletions, and 3) quantify the spectrum of the irregularities associated with density depletions, providing many key observations and discoveries including: -- Global (low latitude) vector DC electric fields at 16 s/sec revealing variations with local time, longitude, and season between extremely low and moderate solar activity-- Reversed zonal E x B drifts below the F-region ledge at sunset and simultaneous observations of large scale Kelvin-Helmholtz instabilities as seeds of spread-F -- Large expanses of quasi-coherent kilometer-scale vector wave observations (electric field and density) below 450 km and their discovery as a source of scintillations using the C/NOFS GPS -- Reversed zonal DC electric fields and simultaneous observations of afternoon counter electrojet-- Enhanced zonal DC electric fields at sunrise and full vector plasma flow continuity at the terminator-- Vector electric field and density irregularities extending to meter-scales within equatorial plasma depletions -- Intense electric field structure within the equatorial ionosphere at night in the absence of density depletions -- First observations of the westward equatorial electrojet in post-midnight ionosphere and possible association with downward meridional drifts-- Spaceborne Dst observations and implications for ring current asymmetries-- Vector electric and magnetic field Poynting flux within depletions and TIDs-- Measurements of ionospheric reflectance and Alfvenic waves associated with TIDs-- Observations of Alfven resonators -- ULF magnetic field structure within the nightime equatorial ionosphere-- Discovery of Schumman resonances in space and implications of ELF radiowave propagation-- Vector observations of 50-60 Hz powerline radiation without harmonics-- Ion cyclotron resonance absorption lines associated with ELF hiss and the identification of ambient ions-- Parallel electric fields associated with lower hybrid waves driven by thunderstorm lightning-related sferics-- Discovery of Z-mode radiation in the equatorial ionosphere including its associated with density depletions-- Optical lightning detector waveforms and electric field sferics observed up to altitudes of 800 km

Robert F Pfaff↗

Acoustic Receptivity of a Blasius Boundary Layer with 2-D and Oblique Surface Waviness

An experimental investigation was conducted to examine acoustic receptivity and subsequent boundary-layer instability evolution for a Blasius boundary layer formed on a flat plate in the presence of two-dimensional (2-D) and oblique (3-D) surface waviness. The effect of the non-localized surface roughness geometry and acoustic wave amplitude on the receptivity process was explored. The surface roughness had a well defined wavenumber spectrum with fundamental wavenumber k (sub w). A planar downstream traveling acoustic wave was created to temporally excite the flow near the resonance frequency of an unstable eigenmode corresponding to k (sub ts) = k (sub w). The range of acoustic forcing levels, epsilon, and roughness heights, DELTA h, examined resulted in a linear dependence of receptivity coefficients; however, the larger values of the forcing combination epsilon dot DELTA h resulted in subsequent nonlinear development of the Tollmien-Schlichting (T-S) wave. This study provided the first experimental evidence of a marked increase in the receptivity coefficient with increasing obliqueness of the surface waviness in excellent agreement with theory. Detuning of the 2-D and oblique disturbances was investigated by varying the streamwise wall-roughness wavenumber a,, and measuring the T-S response. For the configuration where laminar-to-turbulent breakdown occurred, the breakdown process was found to be dominated by energy at the fundamental and harmonic frequencies, indicative of K-type breakdown.

King, Rudolph A.↗

Test Cases for Flutter of the Benchmark Models Rectangular Wings on the Pitch and Plunge Apparatus

The supercritical airfoil was chosen as a relatively modem airfoil for comparison. The BOO12 model was tested first. Three different types of flutter instability boundaries were encountered, a classical flutter boundary, a transonic stall flutter boundary at angle of attack, and a plunge instability near M = 0.9 and for zero angle of attack. This test was made in air and was Transonic Dynamics Tunnel (TDT) Test 468. The BSCW model (for Benchmark SuperCritical Wing) was tested next as TDT Test 470. It was tested using both with air and a heavy gas, R-12, as a test medium. The effect of a transition strip on flutter was evaluated in air. The B64AOlO model was subsequently tested as TDT Test 493. Some further analysis of the experimental data for the BOO12 wing is presented. Transonic calculations using the parameters for the BOO12 wing in a two-dimensional typical section flutter analysis are given. These data are supplemented with data from the Benchmark Active Controls Technology model (BACT) given and in the next chapter of this document. The BACT model was of the same planform and airfoil as the BOO12 model, but with spoilers and a trailing edge control. It was tested in the heavy gas R-12, and was instrumented mostly at the 60 per cent span. The flutter data obtained on PAPA and the static aerodynamic test cases from BACT serve as additional data for the BOO12 model. All three types of flutter are included in the BACT Test Cases. In this report several test cases are selected to illustrate trends for a variety of different conditions with emphasis on transonic flutter. Cases are selected for classical and stall flutter for the BSCW model, for classical and plunge for the B64AOlO model, and for classical flutter for the BOO12 model. Test Cases are also presented for BSCW for static angles of attack. Only the mean pressures and the real and imaginary parts of the first harmonic of the pressures are included in the data for the test cases, but digitized time histories have been archived. The data for the test cases are available as separate electronic files. An overview of the model and tests is given, the standard formulary for these data is listed, and some sample results are presented.

Robert M Bennett↗