Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape approximation”

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 361 records · Page 20

Shape Parameter for a Non-Axisymmetric Isothermal Dendrite

In previous work, we found approximate solutions for paraboloids having perturbations with four-fold axial symmetry in order to model dendritic growth in cubic materials. These solutions provide self-consistent corrections through second order in a shape parameter e to the Peclet number-supercooling relation of the Ivantsov solution. The parameter e is proportional to the amplitude of the four-fold correction to the dendrite shape, as measured from the Ivantsov paraboloid of revolution. We calculate e by comparing the dendrite tip shape to the portion of the equilibrium shape near the growth direction, (001), for anisotropic surface free energy, where the ni are components of the unit normal of the crystal surface. This comparison results in epsilon = -2(epsilon 4), independent of the Peclet number. From the experimental value of epsilon 4, we find epsilon approximately 0.011, in good agreement with the measured value epsilon approximately 0.008 of LaCombe et al.

McFadden, G. B.↗

Demonstration of Minimally Machined Honeycomb Silicon Carbide Mirrors

Honeycomb silicon carbide composite mirrors are made from a carbon fiber preform that is molded into a honeycomb shape using a rigid mold. The carbon fiber honeycomb is densified by using polymer infiltration pyrolysis, or through a reaction with liquid silicon. A chemical vapor deposit, or chemical vapor composite (CVC), process is used to deposit a polishable silicon or silicon carbide cladding on the honeycomb structure. Alternatively, the cladding may be replaced by a freestanding, replicated CVC SiC facesheet that is bonded to the honeycomb. The resulting carbon fiber-reinforced silicon carbide honeycomb structure is a ceramic matrix composite material with high stiffness and mechanical strength, high thermal conductivity, and low CTE (coefficient of thermal expansion). This innovation enables rapid, inexpensive manufacturing. The web thickness of the new material is less than 1 millimeter, and core geometries tailored. These parameters are based on precursor carbon-carbon honeycomb material made and patented by Ultracor. It is estimated at the time of this reporting that the HoneySiC(Trademark) will have a net production cost on the order of $38,000 per square meter. This includes an Ultracor raw material cost of about $97,000 per square meter, and a Trex silicon carbide deposition cost of $27,000 per square meter. Even at double this price, HoneySiC would beat NASA's goal of $100,000 per square meter. Cost savings are estimated to be 40 to 100 times that of current mirror technologies. The organic, rich prepreg material has a density of 56 kilograms per cubic meter. A charred carbon-carbon panel (volatile organics burnt off) has a density of 270 kilograms per cubic meter. Therefore, it is estimated that a HoneySiC panel would have a density of no more than 900 kilograms per cubic meter, which is about half that of beryllium and about onethird the density of bulk silicon carbide. It is also estimated that larger mirrors could be produced in a matter of weeks. Each cell is completely uniform, maintaining the shape of the inserted mandrel. Furthermore, the layup creates pressure that insures node bond strength. Each node is a composite laminate using only the inherent resin system to form the bond. This contrasts starkly with the other known method of producing composite honeycomb, in which individual corrugations are formed, cured, and then bonded together in a secondary process. By varying the size of the mandrels within the layup, varying degrees of density can be achieved. Typical sizes are 3/8 and 3/16 in. (approximately 10 and 5 millimeters). Cell sizes up to 1 in. (approximately 25 millimeters) have been manufactured. Similarly, the shape of the core can be altered for a flexible honeycomb structure.

Goodman, William↗

Polarimetric scattering from layered media with multiple species of scatterers

Geophysical media are usually heterogeneous and contain multiple species of scatterers. In this paper a model is presented to calculate effective permittivities and polarimetric backscattering coefficients of multispecies-layered media. The same physical description is consistently used in the derivation of both permittivities and scattering coefficients. The strong permittivity fluctuation theory is extended to account for the multiple species of scatterers with a general ellipsoidal shape whose orientations are randomly distributed. Under the distorted Born approximation, polarimetric scattering coefficients are obtained. These calculations are applicable to the special cases of spheroidal and spherical scatterers. The model is used to study effects of scatterer shapes and multispecies mixtures on polarimetric signatures of heterogeneous media. The multispecies model accounts for moisture content in scattering media such as snowpack in an ice sheet. The results indicate a high sensitivity of backscatter to moisture with a stronger dependence for drier snow and ice grain size is important to the backscatter. For frost-covered saline ice, model results for bare ice are compared with measured data at C band and then the frost flower formation is simulated with a layer of fanlike ice crystals including brine infiltration over a rough interface. The results with the frost cover suggest a significant increase in scattering coefficients and a polarimetric signature closer to isotropic characteristics compared to the thin saline ice case.

Nghiem, S. V.↗

A Boundary Condition Relaxation Algorithm for Strongly Coupled, Ablating Flows Including Shape Change

Implementations of a model for equilibrium, steady-state ablation boundary conditions are tested for the purpose of providing strong coupling with a hypersonic flow solver. The objective is to remove correction factors or film cooling approximations that are usually applied in coupled implementations of the flow solver and the ablation response. Three test cases are considered - the IRV-2, the Galileo probe, and a notional slender, blunted cone launched at 10 km/s from the Earth's surface. A successive substitution is employed and the order of succession is varied as a function of surface temperature to obtain converged solutions. The implementation is tested on a specified trajectory for the IRV-2 to compute shape change under the approximation of steady-state ablation. Issues associated with stability of the shape change algorithm caused by explicit time step limits are also discussed.

Gnoffo, Peter A.↗

Forebody and afterbody solutions of the Navier-Stokes equations for supersonic flow over blunt bodies in a generalized orthogonal coordinate system

A coordinate transformation, which can approximate many different two-dimensional and axisymmetric body shapes with an analytic function, is used as a basis for solving the Navier-Stokes equations for the purpose of predicting 0 deg angle of attack supersonic flow fields. The transformation defines a curvilinear, orthogonal coordinate system in which coordinate lines are perpendicular to the body and the body is defined by one coordinate line. This system is mapped in to a rectangular computational domain in which the governing flow field equations are solved numerically. Advantages of this technique are that the specification of boundary conditions are simplified and, most importantly, the entire flow field can be obtained, including flow in the wake. Good agreement has been obtained with experimental data for pressure distributions, density distributions, and heat transfer over spheres and cylinders in supersonic flow. Approximations to the Viking aeroshell and to a candidate Jupiter probe are presented and flow fields over these shapes are calculated.

Gnoffo, P. A.↗

X-ray observations of a flare in NGC 4151 from OSO 8

The 2-60-keV flux from NGC 4151 has been observed to change by a factor of 2 on a time scale of 1.5 days. No fluctuations in excess of a factor of 3 are detected on time scales less than 4 hours. During a total observation of approximately 11 days there were no statistically significant changes in spectral shape. The spectrum can be fitted by a power law with photon index of about 1.42 + or - 0.06 and a hydrogen column density of approximately 7.5 + or - 0.5 x 10 to the 22nd power atoms/sq cm. A 2-sigma residual to this fit implies fluorescent Fe line emission with an equivalent width of about 240 eV. Both synchrotron self-Compton and thermal Compton models are consistent with the X-ray data.

Mushotzky, R. F.↗

Estimating the Libration of Mercury by Remote Sensing of Gravity and Altimetry

One of the most intriguing scientific questions about Mercury is whether the planet possesses a liquid core. One way that this can be ascertained is through the measurement of the longitudinal libration of Mercury about its rotation axis. The MESSENGER (Mercury Surface, Space, Environment, Geochemistry, and Ranging) mission to Mercury will attempt this measurement by deriving the rotation of the planet from remotely sensed observations of the Mercury gravity field and the planet's shape. For a fully fluid core the libration is expected to be an approximately 350-m-amplitude oscillation (at the equator) about the mean rotation of the planet. Thus, the shape of the planet undergoes a librational oscillation in the same way that features on the surface participate in the motion. In this sense altimetric topography can be considered analogous to imaging, with the important exception that for altimetry, the long wavelengths corresponding to low degrees and orders in a spherical harmonic expansion, can be precisely determined from orbital measurements. The rotation of the planet thus introduces a time dependence of the topographic image in an inertial (i.e. fixed) reference frame and the libration is a variation in the planetary rotation rate, manifest as an oscillation of the topographic image. Similarly, gravity measures the distribution of mass within a planet and as such can be viewed as providing a vertically integrated 'image' of the internal density structure. For a librating planet the gravitational field will also reveal a variation in the rotation rate, manifest as an oscillation of the field with respect to an inertial reference frame. While the long-wavelength gravity field contains contributions from the radial distribution of mass with depth, the long wavelength terms are most sensitive to the mass distribution of the deep interior. Thus, both the topography and the gravity independently can be used to determine the (irregular) rotation of a planet, although the topographical method measures the rotation of the crust while the gravity measures the rotation of the distributed mass, and these are not necessarily identical. For example, for a planet that contains a fluid outer core that is effectively decoupled from the mantle the gravity field will not librate in the manner of a solid body. The libration obtained from the variable rotation of the gravity field will represent a combination of contributions due to a librating mantle and a differentially rotating fluid core. Differential rotation of the Earth's inner core had been theoretically predicted from three dimensional geodynamical models and subsequently reported from seismic observations. A difference between the libration of the lithosphere as determined from altimetry and the deep interior as determined from gravity, if it could be detected, might provide further insight into the nature of core-mantle coupling. Additional information is contained in the original extended abstract.

Smith, D. E.↗

Quantitative Analysis of Spectral Interference of Spontaneous Raman Scattering in High-Pressure Fuel-Rich H2-Air Combustion

We present a theoretical study of the spectral interferences in the spontaneous Raman scattering spectra of major combustion products in 30-atm fuel-rich hydrogen-air flames. An effective methodology is introduced to choose an appropriate line-shape model for simulating Raman spectra in high-pressure combustion environments. The Voigt profile with the additive approximation assumption was found to provide a reasonable model of the spectral line shape for the present analysis. The rotational/vibrational Raman spectra of H2, N2, and H2O were calculated using an anharmonic-oscillator model using the latest collisional broadening coefficients. The calculated spectra were validated with data obtained in a 10-atm fuel-rich H2-air flame and showed excellent agreement. Our quantitative spectral analysis for equivalence ratios ranging from 1.5 to 5.0 revealed substantial amounts of spectral cross-talk between the rotational H2 lines and the N2 O-/Q-branch; and between the vibrational H2O(0,3) line and the vibrational H2O spectrum. We also address the temperature dependence of the spectral cross-talk and extend our analysis to include a cross-talk compensation technique that removes the nterference arising from the H2 Raman spectra onto the N2, or H2O spectra.

Kojima, Jun↗

Affordable Development and Demonstration of a Small NTR Engine and Stage: How Small is Big Enough?

The Nuclear Thermal Rocket (NTR) derives its energy from fission of uranium-235 atoms contained within fuel elements that comprise the engine's reactor core. It generates high thrust and has a specific impulse potential of approximately 900 seconds - a 100% increase over today's best chemical rockets. The Nuclear Thermal Propulsion (NTP) project, funded by NASA's AES program, includes five key task activities: (1) Recapture, demonstration, and validation of heritage graphite composite (GC) fuel (selected as the "Lead Fuel" option); (2) Engine Conceptual Design; (3) Operating Requirements Definition; (4) Identification of Affordable Options for Ground Testing; and (5) Formulation of an Affordable Development Strategy. During FY'14, a preliminary DDT&E plan and schedule for NTP development was outlined by GRC, DOE and industry that involved significant system-level demonstration projects that included GTD tests at the NNSS, followed by a FTD mission. To reduce cost for the GTD tests and FTD mission, small NTR engines, in either the 7.5 or 16.5 klbf thrust class, were considered. Both engine options used GC fuel and a "common" fuel element (FE) design. The small approximately 7.5 klbf "criticality-limited" engine produces approximately 157 megawatts of thermal power (MWt) and its core is configured with parallel rows of hexagonal-shaped FEs and tie tubes (TTs) with a FE to TT ratio of approximately 1:1. The larger approximately 16.5 klbf Small Nuclear Rocket Engine (SNRE), developed by LANL at the end of the Rover program, produces approximately 367 MWt and has a FE to TT ratio of approximately 2:1. Although both engines use a common 35 inch (approximately 89 cm) long FE, the SNRE's larger diameter core contains approximately 300 more FEs needed to produce an additional 210 MWt of power. To reduce the cost of the FTD mission, a simple "1-burn" lunar flyby mission was considered to reduce the LH2 propellant loading, the stage size and complexity. Use of existing and flight proven liquid rocket and stage hardware (e.g., from the RL10B-2 engine and Delta Cryogenic Second Stage) was also maximized to further aid affordability. This paper examines the pros and cons of using these two small engine options, including their potential to support future human exploration missions to the Moon, near Earth asteroids, and Mars, and recommends a preferred size. It also provides a preliminary assessment of the key activities, development options, and schedule required to affordably build, ground test and fly a small NTR engine and stage within a 10-year timeframe.

Spacecraft design↗

The High-Energy Continuum Emission of the Gamma-Ray Blazar PKS 0528+134

We present Advanced Satellite for Cosmology and Astrophysics (ASCA) observations of the gamma-ray blazar PKS 0528 + 134, obtained at two separate epochs in 1994 August and 1995 March. These data represent the first measurement of the X-ray continuum emission of this source in the medium-hard X-ray band. Both ASCA spectra are consistent with a single power law with photon index GAMMA approximate 1.7-1.8 and column density N(sub H) approximately 5 x 10(exp 21) /sq cm, higher than Galactic. The X-ray flux increased by a factor of 4 in approximately 7 months without appreciable change of the spectral shape. During the lower state of 1994 August, PKS 0528 + 134 was observed simultaneously in the optical, X-rays, and at gamma-ray energies with EGRET. The gamma-ray intensity is the faintest detected thus far in the source, with a steep spectrum (GAMMA approximately 2.7). The extrapolation of the X-ray continuum to the gamma-ray range requires a sharp spectral break at approximately 10(exp 22) Hz. We discuss the radio through gamma-ray spectral energy distribution of PKS 0528 + 134, comparing the low state of 1994 August with the flare state of 1993 March. We show that in PKS 0528 + 134, a non-negligible contribution from the external radiation field is present and that, although synchrotron self-Compton scenarios cannot be ruled out, inverse Compton upscattering of thermal seed photons may be the dominant cooling process for the production of the high-energy continuum in this blazar.

Sambruna, Rita M.↗

Shapes for rotating airfoils

An airfoil which has particular application to the blade or blades of rotor aircraft and aircraft propellers is presented. The airfoil thickness distribution, camber and leading edge radius are shaped to locate the airfoil crest at a more aft position along the chord, and to increase the freestream Mach number at which sonic flow is attained at the airfoil crest. The reduced slope of the airfoil causes a reduction in velocity at the airfoil crest at lift coefficients from zero to the maximum lift coefficient. The leading edge radius is adjusted so that the maximum local Mach number at 1.25 percent chord and at the designed maximum lift coefficient is limited to about 0.48 when the Mach number normal to the leading edge is approximately 0.20. The lower surface leading edge radius is shaped so that the maximum local Mach number at the leading edge is limited to about 0.29 when the Mach number normal to the leading edge is approximately 0.20. The drag divergence Mach number associated with the airfoil is moved to a higher Mach number over a range of lift coefficients resulting in superior aircraft performance.

Bingham, G. J.↗

The peculiar type II supernova 1993J in M81: Transition to the nebular phase

We present optical spectra of the bright, peculiar Type II supernova 1993J in M81 spanning the first 14 months of its existence, revealing its transition to the nebular phase. Unlike the case in normal Type II supernovae, during the first 2-10 months the H-alpha emission line gradually becomes less prominent relative to other features such as (O I) lambda lambda 6300, 6364 and (Ca II) lambda lambda 7291, 7324, as we had predicted based on early-time (tau less than or approximately equal to 2 months) spectra. The nebular spectrum resembles those of the Type Ib/Ic supernovae 1985F and 1987M, although weak H-alpha emission is easily visible even at late times in SN 1993J. At tau = 8 months a close similarity is found with the spectrum of SN 1987K, the only other Type II supernova known to have undergone such a metamorphosis. The emission lines are considerably broader than those of normal Type II supernovae at comparable phases, consistent with the progenitor having lost a majority of its hydrogen envelope prior to exploding. Consequently, there is now little doubt that Type Ib, and probably Type Ic, supernovae result from core collapse in stripped, massive stars; models of the chemical evolution of galaxies in which these subtypes are ascribed to exploding white dwarfs must be appropriately modified. Although all of the emission lines in spectra of SN 1993J fade roughly exponentially for a considerable time, the fading of H-alpha begins to slow down at tau approximately = 8 months, and in the interval tau = 10-14 months its flux is constant, or even slightly rising in the wings of the line. This behavior, together with the box-like shape and great breadth (full width at half maximum (FWHM) approximately = 17 000 km/s) of the line profile, suggests that the H-alpha emission is being produced by the high-velocity outer layer of hydrogen ejecta interacting with circumstellar gas released by the progenitor prior to its explosion. A similar phenomenon has previously been seen at later phases in several Type II supernovae, most notably SN 1980K. Bumps (FWHM approximately = 1000 km/s, amplitude approximately = 20%) in the H-alpha profile are probably indicative of Rayleigh-Taylor instabilities in the cool gas behind the reverse shock. A very narrow component (unresolved, FWHM less than or approximately equal to 200 km/s) of H-alpha at the symmetric velocity of SN 1993J may instead be produced by a superposed H II region, or perhaps by recombination in a large circumstellar shell or ring that was ionized during the first few hours after outburst. In the near future the spectrum of SN 1993J should become increasingly dominated by broad H-alpha emission.

Filippenko, Alexei V.↗

A mobile work station concept for mechanically aided astronaut assembly of large space trusses

This report presents results of a series of truss assembly tests conducted to evaluate a mobile work station concept intended to mechanically assist astronaut manual assembly of erectable space trusses. The tests involved assembly of a tetrahedral truss beam by a pair of test subjects with and without pressure (space) suits, both in Earth gravity and in simulated zero gravity (neutral buoyancy in water). The beam was assembled from 38 identical graphite-epoxy nestable struts, 5.4 m in length with aluminum quick-attachment structural joints. Struts and joints were designed to closely simulate flight hardware. The assembled beam was approximately 16.5 m long and 4.5 m on each of the four sides of its diamond-shaped cross section. The results show that average in-space assembly rates of approximately 38 seconds per strut can be expected for struts of comparable size. This result is virtually independent of the overall size of the structure being assembled. The mobile work station concept would improve astronaut efficiency for on-orbit manual assembly of truss structures, and also this assembly-line method is highly competitive with other construction methods being considered for large space structures.

Heard, W. L., Jr.↗

Exact Integrations of Polynomials and Symmetric Quadrature Formulas over Arbitrary Polyhedral Grids

This paper is concerned with two important elements in the high-order accurate spatial discretization of finite volume equations over arbitrary grids. One element is the integration of basis functions over arbitrary domains, which is used in expressing various spatial integrals in terms of discrete unknowns. The other consists of quadrature approximations to those integrals. Only polynomial basis functions applied to polyhedral and polygonal grids are treated here. Non-triangular polygonal faces are subdivided into a union of planar triangular facets, and the resulting triangulated polyhedron is subdivided into a union of tetrahedra. The straight line segment, triangle, and tetrahedron are thus the fundamental shapes that are the building blocks for all integrations and quadrature approximations. Integrals of products up to the fifth order are derived in a unified manner for the three fundamental shapes in terms of the position vectors of vertices. Results are given both in terms of tensor products and products of Cartesian coordinates. The exact polynomial integrals are used to obtain symmetric quadrature approximations of any degree of precision up to five for arbitrary integrals over the three fundamental domains. Using a coordinate-free formulation, simple and rational procedures are developed to derive virtually all quadrature formulas, including some previously unpublished. Four symmetry groups of quadrature points are introduced to derive Gauss formulas, while their limiting forms are used to derive Lobatto formulas. Representative Gauss and Lobatto formulas are tabulated. The relative efficiency of their application to polyhedral and polygonal grids is detailed. The extension to higher degrees of precision is discussed.

Liu, Yen↗

Unusual locations of Earth's bow shock on September 24 - 25, 1987: Mach number effects

International Sun Earth Explorer 1 (ISEE 1) and Interplanetary Monitoring Platform 8 (IMP 8) data are used to identify 19 crossings of Earth's bow shock during a 30-hour period following 0000 UT on September 24, 1987. Apparent standoff distances for the shock are calculated for each crossing using two methods and the spacecraft location; one method assumes the average shock shape, while the other assumes a ram pressure-dependent shock shape. The shock's apparent standoff distance, normally approximately 14 R(sub E), is shown to increase from near 10 R(sub E) initially to near 19 R(sub E) during an 8-hour period, followed by an excursion to near 35 R(sub E) (where two IMP 8 shock crossings occur) and an eventual return to values smaller than 19 R(sub E). The Alfven M(sub A) and fast magnetosonic M(sub ms). Mach numbers remain above 2 and the number density above 4/cu cm for almost the entire period. Ram pressure effects produce the initial near-Earth shock location, whereas expansions and contractions of the bow shock due to low Mach number effects account, qualitatively and semiquantitatively, for the timing and existence of almost all the remaining ISEE crossings and both IMP 8 crossings. Significant quantitative differences exist between the apparent standoff distances for the shock crossings and those predicted using the observed plasma parameters and the standard model based on Spreiter et al.'s (1966) gasdynamic equation. These differences can be explained in terms of either a different dependence of the standoff distance on Mach number at low M(sub A) and M(sub ms), or variations in shock shape with M(sub A) and M(sub ms) (becoming increasingly "puffed up" with decreasing M(sub A) and M(sub ms), as expected theoretically), or by a combination of both effects.

Cairns, Iver H.↗

Shape optimization for aerodynamic efficiency and low observability

Field methods based on the finite-difference approximations of the time-domain Maxwell's equations and the potential-flow equation have been developed to solve the multidisciplinary problem of airfoil shaping for aerodynamic efficiency and low radar cross section (RCS). A parametric study and an optimization study employing the two analysis methods are presented to illustrate their combined capabilities. The parametric study shows that for frontal radar illumination, the RCS of an airfoil is independent of the chordwise location of maximum thickness but depends strongly on the maximum thickness, leading-edge radius, and leadingedge shape. In addition, this study shows that the RCS of an airfoil can be reduced without significant effects on its transonic aerodynamic efficiency by reducing the leading-edge radius and/or modifying the shape of the leading edge. The optimization study involves the minimization of wave drag for a non-lifting, symmetrical airfoil with constraints on the airfoil maximum thickness and monostatic RCS. This optimization study shows that the two analysis methods can be used effectively to design aerodynamically efficient airfoils with certain desired RCS characteristics.

Vinh, Hoang↗

Application of Newton's method to the postbuckling of rings under pressure loadings

The postbuckling response of circular rings (or long cylinders) is examined. The rings are subjected to four types of external pressure loadings; each type of pressure is defined by its magnitude and direction at points on the buckled ring. Newton's method is applied to the nonlinear differential equations of the exact inextensional theory for the ring problem. A zeroth approximation for the solution of the nonlinear equations, based on the mode shape corresponding to the first buckling pressure, is derived in closed form for each of the four types of pressure. The zeroth approximation is used to start the iteration cycle in Newton's method to compute numerical solutions of the nonlinear equations. The zeroth approximations for the postbuckling pressure-deflection curves are compared with the converged solutions from Newton's method and with similar results reported in the literature.

Thurston, Gaylen A.↗

Theta(1) and gamma Tauri: Understanding the coronal structure of Hyades giants. EUV spectroscopy of BY Dra systems

These programs involve: (1) analysis and interpretation of EUVE spectrometer observations of the active giant beta Cet in comparison to the Hyades giant theta(1) Tau, and (2) analysis and interpretation of EUVE spectrometer observations of the BY Dra systems FK Aqr, DH Leo, and BH Lyn. EUVE carried out observations of beta Cet, but has yet to perform an observation of theta(1) Tau. In beta Cet, a number of Fe lines from high ionization species were observed, up to Fe XXIV. The spectrum overall resembles that seen in the active binary Capella (alpha Aur). All three BY Dra systems were observed by EUVE, and show evidence of high temperature (approximately 10(exp 7) K) plasma; FK Aqr and DH Leo show significant variability in their Deep Survey lightcurves. In FK Aqr, spectral differences between its 'quiescent' and 'active' states suggest possible differences in the plasma density. In DH Leo, the Deep Survey lightcurve, taken over nearly 8 days, shows a distinct period of approximately 1.05 days, similar to the photometric period. The emission measure distributions of all three systems are rather similar in shape, and can be well-represented by a power law with slope approximately 1.5 from 6.2-7.0 in log T.

Stern, Robert A.↗