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 1,117 records · Page 62

Recent advances in antenna synthesis.

A technique for the synthesis of planar aperture antennas has been developed. It applies equally well to continuous or array type apertures. Patterns with shaped main beams and complex sidelobe structures can be synthesized. The method involves applying corrections iteratively and converging toward the desired pattern. The iteration process terminates when the synthesized pattern approximates the desired pattern to within a specified tolerance. A computer program has been written to perform all calculations. The source distribution which is required to produce the synthesized pattern is provided as output. The method is such that patterns requiring complex source distributions will not be synthesized (such as superdirective patterns).

Stutzman, W. L.↗

Wide range multi-channel analog switch

A ten-channel analog switch has been developed as a sample and shape circuit for signals used as inputs to a logarithmic pulse height recording system. The circuit employs a field effect transistor as the switching element. Switching transients are reduced to approximately 2 mV by adding to the gated signal an inverted transient generated by an identical circuit with no signal applied. The circuit is linear to 2% from 3 mV to 10 V.

Levit, L. B.↗

Solutions to bi-Maxwellian transport equations for radial solar wind beyond 28 R(S)

This paper presents solar wind solutions for radial flow between 28 solar radii and 1 AU using the bi-Maxwellian-based 16-moment set of transport equations. In addition to the number density, drift velocity, and parallel and perpendicular temperatures, the 16-moment equations account for the transport of both longitudinal and transverse thermal energies as well as stress. Also, using the 16-moment approximation for the distribution function and assuming plasma parameter values characteristic of the solar wind, contour plots are generated for the proton velocity distribution function. It is shown how the shape of these plots depends on various macroscopic plasma parameters.

Demars, H. G.↗

A laboratory study of subjective response to sonic booms measured at White Sands Missile Range

The Sonic Boom Simulator of the Langley Research Center was used to quantify subjective loudness response to boom signatures consisting of: (1) simulator reproductions of booms recently recorded at White Sands Missile Range; (2) idealized N-waves; and (3) idealized booms having intermediate shocks. The booms with intermediate shocks represented signatures derived from CFD predictions. The recorded booms represented those generated by F15 and T38 aircraft flyovers and represented a variety of waveforms reflecting the effects of propagation through a turbulent atmosphere. These waveforms included the following shape categories: N-waves, peaked, rounded, and U-shaped. Results showed that Perceived Level and Zwicker Loudness Level were good estimators of the loudness of turbulence modified sonic booms. No significant differences were observed between loudness responses for the several shape categories when expressed in terms of Perceived Level. Thus, Perceived Level effectively accounted for waveform differences due to turbulence. Idealized booms with intermediate shocks, however, were rated as being approximately 2.7 dB(PL) less loud than the recorded signatures. This difference was not accounted for by PL.

Sullivan, Brenda M.↗

Gamma-ray pulsars: Emission zones and viewing geometries

There are now a half-dozen young pulsars detected in high-energy photons by the Compton Gamma-Ray Observatory (CGRO), showing a variety of emission efficiencies and pulse profiles. We present here a calculation of the pattern of high-energy emission on the sky in a model which posits gamma-ray production by charge-depleted gaps in the outer magnetosphere. This model accounts for the radio to gamma-ray pulse offsets of the known pulsars, as well as the shape of the high-energy pulse profiles. We also show that about one-third of emitting young radio pulsars will not be detected due to beaming effects, while approximately 2.5 times the number of radio-selected gamma-ray pulsars will be viewed only high energies. Finally we compute the polarization angle variation and find that the previously misunderstood optical polarization sweep of the Crab pulsar arises naturally in this picture. These results strongly support an outer magnetosphere location for the gamma-ray emission.

Romani, Roger W.↗

Computational Study of Surface Tension and Wall Adhesion Effects on an Oil Film Flow Underneath an Air Boundary Layer

The fringe-imaging skin friction (FISF) technique, which was originally developed by D. J. Monson and G. G. Mateer at Ames Research Center and recently extended to 3-D flows, is the most accurate skin friction measurement technique currently available. The principle of this technique is that the skin friction at a point on an aerodynamic surface can be determined by measuring the time-rate-of-change of the thickness of an oil drop placed on the surface under the influence of the external air boundary layer. Lubrication theory is used to relate the oil-patch thickness variation to shear stress. The uncertainty of FISF measurements is estimated to be as low as 4 percent, yet little is known about the effects of surface tension and wall adhesion forces on the measured results. A modified version of the free-surface Navier-Stokes solver RIPPLE, developed at Los Alamos National Laboratories, was used to compute the time development of an oil drop on a surface under a simulated air boundary layer. RIPPLE uses the volume of fluid method to track the surface and the continuum surface force approach to model surface tension and wall adhesion effects. The development of an oil drop, over a time period of approximately 4 seconds, was studied. Under the influence of shear imposed by an air boundary layer, the computed profile of the drop rapidly changes from its initial circular-arc shape to a wedge-like shape. Comparison of the time-varying oil-thickness distributions computed using RIPPLE and also computed using a greatly simplified numerical model of an oil drop equation which does not include surface tension and wall adhesion effects) was used to evaluate the effects of surface tension on FISF measurement results. The effects of surface tension were found to be small but not necessarily negligible in some cases.

Celic, Alan↗

Morphometric analysis of epidermal differentiation in primary roots of Zea mays

Epidermal differentiation in primary roots of Zea mays was divided into six cell types based on cellular shape and cytoplasmic appearance. These six cell types are: 1) apical protoderm, located at the tip of the root pole and characterized by periclinally flattened cells; 2) cuboidal protoderm, located approximately 230 microns from the root pole and characterized by cuboidal cells; 3) tabular epidermis, located approximately 450 microns from the root pole and characterized by anticlinally flattened cells; 4) cuboidal epidermis, located approximately 900 microns from the root pole and characterized by cuboidal cells having numerous small vacuoles; 5) vacuolate cuboidal epidermis, located approximately 1,500 microns from the root pole and characterized by cuboidal cells containing several large vacuoles; and 6) columnar epidermis, located approximately 2,200 microns from the root pole (i.e., at the beginning of the zone of elongation) and characterized by elongated cells. We also used stereology to quantify the cellular changes associated with epidermal differentiation. The quiescent center and the apical protoderm have significantly different ultrastructures. The relative volume of dictyosomes increases dramatically during the early stages of epidermal differentiation. This increase correlates inversely with the amount of coverage provided by the root cap and mucilage.

NASA Discipline Plant Biology↗

L Band Brightness Temperature from Forest: Comparison of Approximate Techniques

In this paper, three approximate physical microwave radiometry models have been used to calculate brightness temperatures from a forest canopy at L-band. These models are (1) tau-omega model (zero order scattering approximation to radiative transfer equations), (2) successive order of scattering model up to first order (first order scattering approximation to the radiative transfer equations), and (3) Peak technique utilizing the active solution obtained from the Distorted Born Approximation (DBA). These models are physically-based and treat vegetation as a layer of discrete scatterers over a rough surface. Vegetation components within the canopy are represented by canonical shapes such as dielectric discs and cylinders. The tau-omega model is based on a zero-order solution to the radiative transfer (RT) equations. The model ignores scattering except for the effect of the scatterers in the attenuation of the emission through the vegetation. Application of the tau-omega model to data acquired during airborne and ground-based campaigns over the years has solidified scientific understanding of microwave interactions with different landscapes. In particular, shrubland, grasslands, agricultural crops, and light to moderate vegetation have been investigated. Its applicability to areas with a significant tree fraction is unknown. The first order scattering model is based on an iterative solution of the RT equation up to the first order. The first order solution is obtained by substituting the zeroth-order solution into the scattering source term and then solving the resulting radiative transfer equations. This formulation adds a new scattering term to the tau-omega model. It represents emission by particles in the layer and emission by the ground that is scattered once by particles in the layer. The resulting model represents an improvement over the standard zero-order solution (the tau-omega model) since it accounts for the scattered vegetation and ground radiation that can have a pronounced effect on the observed brightness temperature. The third model is based on the Peake formulation in conjunction with the DBA. The procedure for calculation of forest emission is accomplished by first calculating the bistatic scattering cross section for each type of scatterer, then by using the DBA to calculate specular albedo of the ground and the diffused albedo of the layer. Once the albedos are determined, Peake s principle relating active and passive problems can be used to determine the effective emissivity of the forest layer.

Kurum, Mehmet↗

Understanding Effective Diameter and Its Application to Terrestrial Radiation in Ice Clouds

The cloud property known as "effective diameter" or "effective radius", which in essence is the cloud particle size distribution (PSD) volume at bulk density divided by its projected area, is used extensively in atmospheric radiation transfer, climate modeling and remote sensing. This derives from the assumption that PSD optical properties can be uniquely described in terms of their effective diameter, D(sub e), and their cloud water content (CWC), henceforth referred to as the D(sub e)-CWC assumption. This study challenges this assumption, showing that while the D(sub e)-CWC assumption appears generally valid for liquid water clouds, it appears less valid for ice clouds in regions where (1) absorption is not primarily a function of either the PSD ice water content (IWC) or the PSD projected area, and (2) where wave resonance (i.e. photon tunneling) contributes significantly to absorption. These two regions often strongly coincide at terrestrial wavelengths when De less than ~60 m, which is where this D(sub e)-CWC assumption appears poorest. Treating optical properties solely in terms of D(sub e) and IWC may lead to errors up to 24%, 26% and 20% for terrestrial radiation in the window region regarding the absorption and extinction coefficients and the single scattering albedo, respectively. Outside the window region, errors may reach 33% and 42% regarding absorption and extinction. The magnitude and sign of these errors can change rapidly with wavelength, which may produce significant errors in climate modeling, remote sensing and other applications concerned with the wavelength dependence of radiation. Where the D(sub e)-CWC assumption breaks down, ice cloud optical properties appear to depend on D(sub e), IWC and the PSD shape. Optical property parameterizations in climate models and remote sensing algorithms based on historical PSD measurements may exhibit errors due to previously unknown PSD errors (i.e. the presence of ice artifacts due to the shattering of larger ice particles on the probe inlet tube during sampling). More recently developed cloud probes are designed to mitigate this shattering problem. Using realistic PSD shapes for a given temperature (and/or IWC) and cloud type may minimize errors associated with PSD shape in ice optics parameterizations and remote sensing algorithms. While this topic was investigated using two ice optics schemes (the Yang et al., 2005 database and the modified anomalous diffraction approximation, or MADA), a physical understanding of the limitations of the D(sub e)-IWC assumption was made possible by using MADA. MADA allows one to approximate the contribution of photon tunneling to absorption relative to other optical processes, which reveals that part of the error regarding the D(sub e)-IWC assumption can be associated with tunneling. By relating the remaining error to the radiation penetration depth in bulk ice (DELTA L) due to absorption, the domain where the D(sub e)-IWC assumption is weakest was described in terms of D(sub e) and DELTA L.

Mitchell, D. L.↗

A scaling law of the neutral penetration length and Balmer- α wing shape in high-temperature plasmas

Hydrogen atoms penetrating deep inside high-temperature magnetically confined plasmas by repetitive charge-exchange collisions result in a particle source of the plasma, which affects the plasma transport significantly. In this paper, we present an approximate solution of the fluid equations for neutral transport and an analytical representation of the neutral penetration length, in a simplified plasma geometry. This analysis predicts a power-law decay in the Balmer-α line wings which reflects the velocity distribution of the neutral atoms, with the power-law exponent analytically represented as well. These scaling laws are compared with a simple Monte–Carlo simulation and spectroscopic observations of Large Helical Device plasmas. Since the Balmer-α line wings are experimentally accessible, our formulation opens the possibility to quickly estimate the neutral penetration length from spectroscopic observations.

neutral opacity↗

Details of the computed flowfield over a circular cylinder at Reynolds number 1200

The application of an upwind-biased implicit approximate factorization Navier-Stokes algorithm to the unsteady impulsive start-up flow over a circular cylinder at Reynolds number 1200 is described. The complete form of the compressible Navier-Stokes equations is used, and the algorithm is second-order accurate in both space and time. The development with time of the shape and size of the separated vortical flow region is computed, as well as the time-variation of several boundary layer parameters and profile shapes. Computations, in general, show excellent agreement with experiment, although the present method predicts a more rapid onset of reversed flow on the cylinder than evidenced in experiment. The changes that the vortical region behind the cylinder undergoes as the symmetric flow transitions to periodic vortex shedding are discussed. The flow becomes periodic with a Strouhal frequency of 0.222, which compares well with the experimental value of approximately 0.21. The effect of grid density on the development of the unsteady flow is also shown.

Rumsey, C. L.↗

Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping

A multi-beam, multi-color, polarized laser transmitter has been developed for mapping applications. It uses commercial off-the-shelf components for a lowcost approach for a ruggedized laser suitable for field deployment. The laser transmitter design is capable of delivering dual wavelengths, multiple beams on each wavelength with equal (or variable) intensities per beam, and a welldefined state of polarization. This laser transmitter has been flown on several airborne campaigns for the Slope Imaging Multi-Polarization Photon Counting Lidar (SIMPL) instrument, and at the time of this reporting is at a technology readiness level of between 5 and 6. The laser is a 1,064-nm microchip high-repetition-rate laser emitting energy of about 8 microjoules per pulse. The beam was frequency-doubled to 532 nm using a KTP (KTiOPO4) nonlinear crystal [other nonlinear crystals such as LBO (LiB3O5) or periodically poled lithium niobiate can be used as well, depending on the conversion efficiency requirements], and the conversion efficiency was approximately 30 percent. The KTP was under temperature control using a thermoelectric cooler and a feedback monitoring thermistor. The dual-wavelength beams were then spectrally separated and each color went through its own optical path, which consisted of a beam-shaping lens, quarterwave plate (QWP), and a birefringent crystal (in this case, a calcite crystal, but others such as vanadate can be used). The QWP and calcite crystal set was used to convert the laser beams from a linearly polarized state to circularly polarized light, which when injected into a calcite crystal, will spatially separate the circularly polarized light into the two linear polarized components. The spatial separation of the two linearly polarized components is determined by the length of the crystal. A second set of QWP and calcite then further separated the two beams into four. Additional sets of QWP and calcite can be used to further split the beams into multiple orders of two. The spatially separated beams had alternating linearly polarization states; a half-wave plate (HWP) array was then made to rotate the alternating states of A multi-beam, multi-color, polarized laser transmitter has been developed for mapping applications. It uses commercial off-the-shelf components for a lowcost approach for a ruggedized laser suitable for field deployment. The laser transmitter design is capable of delivering dual wavelengths, multiple beams on each wavelength with equal (or variable) intensities per beam, and a welldefined state of polarization. This laser transmitter has been flown on several airborne campaigns for the Slope Imaging Multi-Polarization Photon Counting Lidar (SIMPL) instrument, and at the time of this reporting is at a technology readiness level of between 5 and 6. The laser is a 1,064-nm microchip high-repetition-rate laser emitting energy of about 8 microjoules per pulse. The beam was frequency-doubled to 532 nm using a KTP (KTiOPO4) nonlinear crystal [other nonlinear crystals such as LBO (LiB3O5) or periodically poled lithium niobiate can be used as well, depending on the conversion efficiency requirements], and the conversion efficiency was approximately 30 percent. The KTP was under temperature control using a thermoelectric cooler and a feedback monitoring thermistor. The dual-wavelength beams were then spectrally separated and each color went through its own optical path, which consisted of a beam-shaping lens, quarterwave plate (QWP), and a birefringent crystal (in this case, a calcite crystal, but others such as vanadate can be used). The QWP and calcite crystal set was used to convert the laser beams from a linearly polarized state to circularly polarized light, which when injected into a calcite crystal, will spatially separate the circularly polarized light into the two linear polarized components. The spatial separation of the two linearly polarized components is determined by the length of the crystal. A cond set of QWP and calcite then further separated the two beams into four. Additional sets of QWP and calcite can be used to further split the beams into multiple orders of two. The spatially separated beams had alternating linearly polarization states; a half-wave plate (HWP) array was then made to rotate the alternating states of

Yu, Anthony W.↗

Acceleration of interstellar pickup ions in the disturbed solar wind observed on Ulysses

Acceleration of interstellar pickup H(+) and He(+) as well as of solar wind protons and alpha particles has been observed on Ulysses during the passage of a corotating interaction region (CIR) at approximately 4.5 AU. Injection efficiencies for both the high thermal speed interstellar pickup ions (H(+) and He(+)) and the low thermal speed solar wind ions (H(+) and He(++) are derived using velocity distribution functions of protons, pickup He(+) and alpha particles from less than 1 to 60 keV/e and of ions (principally protons) above approximately 60 keV. The observed spatial variations of the few keV and the few hundred keV accelerated pickup protons across the forward shock of CIR indicate a two stage acceleration mechanism. Thermal ions are first accelerated to speeds of 3 to 4 times the solar wind speed inside the CIR, presumably by some statistical mechanism, before reaching higher energies by a shock acceleration process. Our results also indicate that (1) the injection efficiencies for pickup ions are almost 100 times higher than they are for solar wind ions, (2) pickup H(+) and He(+) are the two most abundant suprathermal ion species and they carry a large fraction of the particle thermal pressure, (3) the injection efficiency is highest for protons, lowest for He(+), and intermediate for alpha particles, (4) both H(+) and He(+) have identical spectral shapes above the cutoff speed for pickup ions, and (5) the solar wind frame velocity distribution function of protons has the form F(w) = F(sub o)w(sup -4) for 1 is less than w is less than approximately 5, where w is the ion speed divided by the solar wind speed. Above w approximately 5-10 the proton spectrum becomes steeper. These results have important implications concerning acceleration of ions by shocks and CIRs, acceleration of anomalous cosmic rays, and particle dynamics in the outer heliosphere.

Gloeckler, G.↗

Near-Net-Shape Processing of Sintered Fibrous Ceramics Achieved

A variety of sintered fibrous ceramic (SFC) materials have been developed over the last 50 years as thermal barrier materials for reentry applications. SFC materials typically exhibit very low thermal conductivities combined with low densities and good thermal stability up to 2500 F. These materials have flown successfully on the space shuttle orbiters since the 1960's. More recently, the McDonnell Douglas Corporation successfully used SFC tiles as a heat shield on the underside of its DC X test vehicle. For both of these applications, tiles are machined from blocks of a specific type of SFC called an alumina-enhanced thermal barrier (AETB). The sizes of these blocks have been limited by the manufacturing process. In addition, as much as 80 to 90 percent of the material can be lost during the machining of tiles with significant amounts of curvature. To address these problems, the NASA Glenn Research Center at Lewis Field entered a cooperative contract with the Boeing Company to develop a vacuum-assisted forming process that can produce large (approximately 4 square feet), severely contoured panels of AETB while saving costs in comparison to the conventional cast-and-machine billet process. For shuttle use, AETB is slurry cast, drained, and fired to form square billets conforming to the shape of the filtration box. The billets are then cut into tiles of the appropriate size for thermally protecting the space shuttle. Processing techniques have limited the maximum size of AETB billets to 21.5 square inches by 6.5-in. thick, but the space shuttles use discrete heat shield tiles no more than 8 to 12 square inches. However, in other applications, large, complex shapes are needed, and the tiling approach is undesirable. For such applications, vacuum-assisted forming can produce large parts with complex shapes while reducing machining waste and eliminating cemented joints between bonded billets. Because it allows contoured shapes to be formed, material utilization is inherently high. Initial estimates show that the amount of material lost during machining can be reduced by 50 percent or more. In addition, a fiber alignment favorable for minimum heat transfer is maintained for all panel shapes since the fibers are aligned parallel to the contoured surface of the forming tool or mold. The vacuum-assisted forming process can complete the entire forming operation in a matter of minutes and can produce multiple parts whose size is limited only by the size of the forming tool. To date, panels as large as 2 square feet have been demonstrated The vacuum-assisted forming process starts with the fabrication of a permeable forming tool, or mold, with the proper part contour. This reusable tool is mounted over an internal rib support structure, as depicted in the diagram, such that a vacuum can be pulled on the bottom portion of the tool. AETB slurry is then poured over and around the tool, liquid is drawn from the slurry, and the part forms over the tool surface. The part is then dried, fired, and finished machined. Future plans include an evaluation of the need for additional coatings and surface-toughness treatments to extend the durability and performance of this material.

Angel, Paul W.↗

A survey of nonlinear robust optimization

Robust optimization (RO) has attracted great attention from the optimization community over the past decade. RO is dedicated to solving optimization problems subject to uncertainty: design constraints must be satisfied for all the values of the uncertain parameters within a given uncertainty set. Uncertainty sets may be modeled as deterministic sets (boxes, polyhedra, ellipsoids), in which case the RO problem may be reformulated via worst-case analysis, or as families of distributions. The challenge of RO is to reformulate or approximate robust constraints so that the uncertain optimization problem is transformed into a tractable deterministic optimization problem. Most reformulation methods assume linearity of the robust constraints or uncertainty sets of favorable shape, which represents only a fraction of real-world applications. This survey addresses nonlinear RO and includes problem formulations and applications, solution approaches, and available software with code samples.

97 MATHEMATICS AND COMPUTING↗

Moments-based interface reconstruction, remap and advection

We present a new moment-of-fluid (MOF 2 ) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to reconstruct a convex material polygon or a union of convex polygons that approximate the respective material fragment. The new method requires information about the material moments only for the cell under consideration. Furthermore, the MOF 2 method allows to exactly reproduce several convex shapes: corners, filaments, and some concave shapes: cell-complements to corners and filaments.

97 MATHEMATICS AND COMPUTING↗

Financial Analysis of the High Flow Experiment conducted at the Glen Canyon Dam during Water Year 2023

The Glen Canyon Dam (GCD) is a Colorado River Storage Project (CRSP) power resource that is a component of the Salt Lake City Area Integrated Projects (SLCA/IP). The 2016 record of decision (ROD) for the GCD long-term experimental and management plan (LTEMP) final Environmental Impact Statement (EIS) specified criteria for GCD monthly water releases, daily and hourly operating limits, and experimental releases. This report examines the financial implications of the high flow experiment (HFE) conducted at GCD during the spring of Water Year (WY) 2023 as required by the LTEMP HFE Protocol. This report is part of a series of reports that describe the financial costs of LTEMP experimental releases since the 2016 ROD was adopted in January 2017. Previous reports analyzed the impact of several past HFEs and Bug Flow Experiments. This report focuses on the HFE conducted in April 2023. For this experimental release, financial costs of approximately $1.33 million were incurred because the HFE required sustained water releases exceeding the power plant’s maximum turbine flow rate. In addition, during the experiment, operators were not allowed to shape GCD power production, either to follow Firm Electric Service (FES) customer day-ahead energy deliveries or to respond to market prices. This study identifies the main factors contributing to the HFE costs and examines the interdependencies among these factors. It applies an integrated set of tools to estimate Western Area Power Administration (WAPA) financial impacts by simulating GCD under two types of cases; namely, (1) a “With Experiment” case that mimics the operations that actually occurred and (2) a “Without Experiment” case that simulates operations under the assumption that the HFE did not occur. The “With Experiment” case mimics operations during the HFE and the entire month the HFE occurred. It complies with LTEMP hourly and daily operating criteria. The “Without Experiment” case assumes that the HFE did not occur. The monthly water release volume is assumed to be identical under both cases. The Colorado River Storage Project Python-based model (CRiSPPy) model was the main modeling tool used to simulate the dispatch of the GCD hydropower plant and associated water releases from Lake Powell. In the modeling process, the research team used extensive data sets and historical information on SLCA/IP power plant characteristics, hydrologic conditions, and WAPA’s power purchases and sales prices. In addition to estimating the financial impact of the HFE, the team used the CRiSPPy model to gain insights into the interplay among ROD operating criteria, exceptions made to criteria to accommodate the HFE, and WAPA operating practices.

13 HYDRO ENERGY↗

Analysis of the ultrafine fraction of the Apollo 14 regolith

Analyses were obtained on more than 2400 randomly selected particles from the sub-37 micron (ultrafine) fraction of ten Apollo 14 regolith samples. The analyses were conducted with an energy dispersive electron microprobe system. The semiquantitative data were used to group the particles into ten categories. The pyroxene/plagioclase and olivine/plagioclase ratios are inconsistent with those ratios in the Apollo 14 breccias and rocks. The data suggest that fragmented basalts similar to Apollo 12 olivine basalts may have made significant contributions to the ultrafine fraction of the Fra Mauro regolith. Among a number of unusual particles encountered are brown, birefringent lath-shaped grains with 60 wt % SiO2 and 34 wt % FeO(FeSi2O5) and a glass with 20 to 25 wt % CaO, 0 to 8 wt % MgO, 40 to 45 wt % Al2O3 and approximately 30 wt % SiO2.

Finkelman, R. B.↗