Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “particle distributions”

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 793 records · Page 44

Phase space distribution functions and energy distributions of dark matter particles in haloes

ABSTRACT For a spherical dark matter halo with isotropic velocity distribution, the phase space distribution function (DF), the energy distribution, and the density profile form a set of self-consistent description of its equilibrium state, and knowing one is sufficient to determine the other two. The Navarro–Frenk–White density profile (NFW profile) is known to be a good approximation to the spherically averaged density distribution in simulated haloes. The DARKexp energy distribution is also known to compare well with the simulated energy distribution. We present a quantitative assessment of the NFW and DARKexp fits to the simulated DF and energy distribution for a wide range of haloes in a dark-matter-only simulation from the IllustrisTNG Project. As expected, we find that the NFW fits work well except at low energy when the density at small radii deviates from the NFW profile. Further, the NFW and DARKexp fits have comparable accuracy in the region where both fit well, but the DARKexp fits are better at low energy because they require matching of the central gravitational potential. We also find an approximate relation between the energy scale parametrizing the DARKexp energy distribution and that defined by the characteristic density and radius of the NFW profile. This relation may be linked to the relaxation process during halo formation.

Astronomy & Astrophysics↗

3002 Humidified Tandem Differential Mobility Analyzer (HTDMA) Instrument Handbook

The Brechtel Manufacturing Inc. (BMI) Humidified Tandem Differential Mobility Analyzer (HT-DMA Model 3002) (Brechtel and Kreidenweis 2000a,b, Henning et al. 2005, Xerxes et al. 2014) measures how aerosol particles of different initial dry sizes grow or shrink when exposed to changing relative humidity (RH) conditions. It uses two different mobility analyzers (DMA) and a humidification system to make the measurements. One DMA selects a narrow size range of dry aerosol particles, which are exposed to varying RH conditions in the humidification system. The second (humidified) DMA scans the particle size distribution output from the humidification system. Scanning a wide range of particle sizes enables the second DMA to measure changes in size or growth factor (growth factor = humidified size/dry size), due to water uptake by the particles. A Condensation Particle Counter (CPC) downstream of the second DMA counts particles as a function of selected size in order to obtain the number size distribution of particles exposed to different RH conditions.

54 ENVIRONMENTAL SCIENCES↗

Multi-Parameter Scattering Sensor and Methods

Methods, detectors and systems detect particles and/or measure particle properties. According to one embodiment, a detector for detecting particles comprises: a sensor for receiving radiation scattered by an ensemble of particles; and a processor for determining a physical parameter for the detector, or an optimal detection angle or a bound for an optimal detection angle, for measuring at least one moment or integrated moment of the ensemble of particles, the physical parameter, or detection angle, or detection angle bound being determined based on one or more of properties (a) and/or (b) and/or (c) and/or (d) or ranges for one or more of properties (a) and/or (b) and/or (c) and/or (d), wherein (a)-(d) are the following: (a) is a wavelength of light incident on the particles, (b) is a count median diameter or other characteristic size parameter of the particle size distribution, (c) is a standard deviation or other characteristic width parameter of the particle size distribution, and (d) is a refractive index of particles.

Greenberg, Paul S.↗

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗

Jupiter's radiation belts and upper atmosphere

Models of Jupiter's radiation belts have been constructed to determine the distribution of particles and their energies which will produce the observed decimetric radio emission. Data on the spectrum and the variation of emission with Jovian longitude have been used to show that the relativistic particles have a nearly isotropic distribution with high energies (of order 100 MeV) within 2 Jovian radii and a very flat distribution in the equatorial plane of low energy particles further out in the magnetosphere. Subtraction of the emission predicted by this model from the total radio emission shows that the thermal contribution in the frequency range between 3000 and 10,000 MHz is somewhat less than had been previously expected. (The brightness temperature of the planetary disk is 180 K at 3,000 MHz, for example.) This suggests that the ammonia mixing ratio in Jupiter's upper atmosphere may be as high as 0.002.

Degioanni, J. J.↗

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗

Polysiloxane coatings on microspheres based on Multiphase Flow with Interface Exchange-Discrete Element Modelling

Polysiloxane coatings on yttria stabilized zirconia (YSZ) microspheres of 500 μm were simulated in a spout fluidized bed coater using Multiphase Flow with Interface Exchange-Discrete Element Modelling (MFiX-DEM). In this work, two different coater configurations were developed to study the influence of gas velocity and its distribution on particle dynamics. The presence of the Wurster tube not only enhances the distribution but also increases the overall residence time of the particles. Investigations were also carried out with different Wurster tube positions (normal, 10% and 20% lowered from its initial position). Among these, 20% lowered Wurster tube position demonstrated the most effective coating process. The effects of gas inlet pressure on the average gas velocity and the distribution of particles were analyzed. More than 97% of the particles can be retained. The derived results, including average gas velocity, particle retention percentage, and distribution of particles with gas velocity, are being used to guide the experimental work in obtaining defect-free coatings for YSZ microspheres.

36 MATERIALS SCIENCE↗

Decending motion of particle and its effect on ozone hole chemistry

Particle descending motion is one possible process which causes ozone loss near the tropopause in the Antarctic spring. However, this particle size distribution has not yet been measured. Particle settling is an important redistribution process of the chemical constituents contained in the particles. To understand particle settling effects on the Ozone Hole, information on the size distribution and the chemical composition of the particles is necessary.

Iwasaka, Y.↗

Experiments on excitation of Alfvén eigenmodes by alpha-particles with bump-on-tail distribution in JET DTE2 plasmas

Dedicated experiments were performed in JET DTE2 plasmas for obtaining an α-particle bump-on-tail (BOT) distribution aiming at exciting Alfvén eigenmodes (AEs). Neutral beam injection-only heating with modulated power was used so that fusion-born α-particles were the only ions present in the MeV energy range in these DT plasmas. The beam power modulation on a time scale shorter than the α-particle slowing down time was chosen for modulating the α-particle source and thus sustaining a BOT in the α-particle distribution. High-frequency modes in the toroidicity-induced Alfven eigenmode (TAE) frequency range and multiple short-lived modes in a wider frequency range have been detected in these DT discharges with interferometry, soft x-ray cameras, and reflectometry. The modes observed were localised close to the magnetic axis and were not seen in the Mirnov coils. Analysis with the TRANSP and Fokker-Planck FIDIT codes confirms that α-particle distributions with BOT in energy were achieved during some time intervals in these discharges though no clear correlation was found between the times of the high-frequency mode excitation and the BOT time intervals. The combined magneto-hydrodynamic (MHD) and kinetic modelling studies show that the high-frequency mode in the TAE frequency range is best fitted with a TAE of toroidal mode number n= 9. This mode is driven mostly by the on-axis beam ions while the smaller drive due to the pressure gradient of α-particles allows overcoming the marginal stability and exciting the mode. The observed multiple short-lived modes in a wider frequency range are identified as the on-axis kinetic AEs predicted in Rosenbluth and Rutherford.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Highly Durable and Active Pt/Sb-Doped SnO 2 Oxygen Reduction Reaction Electrocatalysts Produced by Atomic Layer Deposition

Platinum supported on mixed-metal oxides (MMOs) are a class of active and durable cathode catalysts for proton exchange membrane fuel cell (PEMFC) due to a combination of the high oxidative stability of the supports and strong-metal-support interactions (SMSI) that enable them to exceed the activity of Pt/C. Herein, we solve a significant remaining challenge with Pt/MMO systems, namely the relatively low surface area and porosity. This is achieved by dispersing nearly-uniform Pt clusters using atomic layer deposition (ALD) on highly conductive (6.2 S/cm) and stable antimony doped tin dioxide (ATO) support. ALD-Pt/ATO exhibited significantly higher electrochemically active surface area (ECSA) (74 m 2 /g) and oxygen reduction reaction (ORR) catalytic activity (102 mA/mg Pt at 0.9 V vs. RHE) compared to Pt/ATO synthesized using ethylene glycol (ECSA=31 m 2 /g Pt , mass activity=52 mA/mg Pt at 0.9 V vs. RHE) and formic acid reduction methods (ECSA=28 m 2 /g Pt , mass activity=46 mA/mg Pt at 0.9 V vs. RHE). Further characterization showed that wet chemical methods resulted in poorer Pt particle dispersion, poor control over Pt particle size distribution and chemical degradation of the support (during Pt deposition). Given the near-ideal Pt particle size distribution of the ALD-Pt/ATO, particle size growth and loss of ECSA was found to be minimal over the course of rigorous potential cycling. Thus, after 10,000 potential cycles between 1V and 1.5V vs. RHE, ALD-Pt/ATO and other Pt/ATOs were found to retain 100% of their initial ECSA compared to 57.6% retention for Pt/C. Upon testing in a H 2 /air PEMFC, following 1,000 potential cycles, the change in ALD-Pt/ATO performance was negligible while Pt/C exhibited a 68.2% loss of initial peak power density. Thus, ALD-Pt/ATO is an active and highly durable ORR electrocatalyst in PEMFCs under start-up-shut down conditions.

25 ENERGY STORAGE↗

Intercomparison and Evaluation of Global Aerosol Microphysical Properties Among Aerocom Models of a Range of Complexity

Many of the next generation of global climate models will include aerosol schemes which explicitly simulate the microphysical processes that determine the particle size distribution. These models enable aerosol optical properties and cloud condensation nuclei (CCN) concentrations to be determined by fundamental aerosol processes, which should lead to a more physically based simulation of aerosol direct and indirect radiative forcings. This study examines the global variation in particle size distribution simulated by 12 global aerosol microphysics models to quantify model diversity and to identify any common biases against observations. Evaluation against size distribution measurements from a new European network of aerosol supersites shows that the mean model agrees quite well with the observations at many sites on the annual mean, but there are some seasonal biases common to many sites. In particular, at many of these European sites, the accumulation mode number concentration is biased low during winter and Aitken mode concentrations tend to be overestimated in winter and underestimated in summer. At high northern latitudes, the models strongly underpredict Aitken and accumulation particle concentrations compared to the measurements, consistent with previous studies that have highlighted the poor performance of global aerosol models in the Arctic. In the marine boundary layer, the models capture the observed meridional variation in the size distribution, which is dominated by the Aitken mode at high latitudes, with an increasing concentration of accumulation particles with decreasing latitude. Considering vertical profiles, the models reproduce the observed peak in total particle concentrations in the upper troposphere due to new particle formation, although modelled peak concentrations tend to be biased high over Europe. Overall, the multimodel- mean data set simulates the global variation of the particle size distribution with a good degree of skill, suggesting that most of the individual global aerosol microphysics models are performing well, although the large model diversity indicates that some models are in poor agreement with the observations. Further work is required to better constrain size-resolved primary and secondary particle number sources, and an improved understanding of nucleation an growth (e.g. the role of nitrate and secondary organics) will improve the fidelity of simulated particle size distributions.

Cloud physics↗

Application of Aerosol Hygroscopicity Measured at the Atmospheric Radiation Measurement Program's Southern Great Plains Site to Examine Composition and Evolution

A Differential Mobility Analyzer/Tandem Differential Mobility Analyzer (DMA/TDMA) was used to measure submicron aerosol size distributions, hygroscopicity, and occasionally volatility during the May 2003 Aerosol Intensive Operational Period (IOP) at the Central Facility of the Atmospheric Radiation Measurement Program's Southern Great Plains (ARM SGP) site. Hygroscopic growth factor distributions for particles at eight dry diameters ranging from 0.012 micrometers to 0.600 micrometers were measured throughout the study. For a subset of particle sizes, more detailed measurements were occasionally made in which the relative humidity or temperature to which the aerosol was exposed was varied over a wide range. These measurements, in conjunction with backtrajectory clustering, were used to infer aerosol composition and to gain insight into the processes responsible for evolution. The hygroscopic growth of both the smallest and largest particles analyzed was typically less than that of particles with dry diameters of about 0.100 micrometers. It is speculated that condensation of secondary organic aerosol on nucleation mode particles is largely responsible for the minimal hygroscopic growth observed at the smallest sizes considered. Growth factor distributions of the largest particles characterized typically contained a nonhygroscopic mode believed to be composed primarily of dust. A model was developed to characterize the hygroscopic properties of particles within a size distribution mode through analysis of the fixed size hygroscopic growth measurements. The performance of this model was quantified through comparison of the measured fixed size hygroscopic growth factor distributions with those simulated through convolution of the size-resolved concentration contributed by each of the size modes and the mode-resolved hygroscopicity. This transformation from sizeresolved hygroscopicity to mode-resolved hygroscopicity facilitated examination of changes in the hygroscopic properties of particles within a size distribution mode that accompanied changes in the sizes of those particles. This model was used to examine three specific cases in which the sampled aerosol evolved slowly over a period of hours or days.

Gasparini, Roberto↗

Scanning Mobility Particle Sizer (SMPS)-Aerodynamic Particle Sizer (APS) Merged Size Distribution (mergedsmpsaps) Value-Added Product Report

Aerosol particles influence the Earth’s radiation balance directly by absorbing and scattering light and indirectly by influencing cloud formation, properties, and lifetimes. Measurements of aerosol particle optical properties, mass loading, size distributions, microphysical properties, cloud formation properties, and chemical composition are important for understanding the aerosol life cycle and for validating earth system models that predict these quantities. The Atmospheric Radiation Measurement (ARM) user facility’s Aerosol Observing System (AOS) is a highly instrumented platform designed to house instruments that measure many of these aerosol properties in situ. Currently, ARM operates at least four different instruments that measure a portion of the ambient aerosol size distribution. Most users are interested in the entire size distribution or a portion of the size distribution that extends across the measurement range of multiple instruments. However, merging these distributions is not trivial as the instruments employ different measurement principals and, in most cases, report data as a function of different representations of the aerosol diameter.

54 ENVIRONMENTAL SCIENCES↗

Characterization of Iron and Organic Carbon Colloids in Boreal Rivers and Their Fate at High Salinity

Riverine colloids are important carriers of macronutrients, trace metals, and pollutants into marine waters. The aim of the current study was to extend the understanding of iron (Fe) and organic carbon (OC) colloids in boreal rivers and their fate at higher salinities. X-ray absorbance spectroscopy (XAS) and dynamic light scattering (DLS) were combined to explore Fe speciation and colloidal characteristics such as size and surface charge and how these are affected at increasing salinity. XAS confirmed the presence of two Fe phases in the river waters—Fe-organic matter (OM) complexes and Fe(oxy)hydroxides. From DLS measurements on filtered and unfiltered samples, three particle size distributions were identified. The smallest particles (10–40 nm) were positively charged and suggested to consist of essentially bare Fe (oxy)hydroxide nanoparticles. The largest particles (300–900 nm) were dominated by Fe(oxy)hydroxides associated with chromophoric molecular matter. An intermediate size distribution (100–200 nm) with a negative surface charge was presumably dominated by OM and containing Fe-OM complexes. Increasing the salinity resulted in a removal of the smallest distribution. Unexpectedly, both the intermediate and largest size distributions were still detected at high salinity. The collective results suggest that Fe(oxy) hydroxides and Fe-OM complexes are both found across the wide size range studied and that colloidal size does not necessarily reflect either Fe speciation or stability toward salinity-induced aggregation. The findings further demonstrate that also particles beyond the typically studied <0.45-μm size range should be considered to fully understand the riverine transport and fate of macronutrients, trace metals, and pollutants.

54 ENVIRONMENTAL SCIENCES↗

Optical measurement of interplanetary particulates from Pioneer 10

The spatial concentration and size distribution for particulates measured in situ by the asteroid/meteoroid detector on Pioneer 10 between 1.0 and 3.5 AU are presented. The size distribution includes particles of radii from about 35 microns to 10 cm. Extrapolation from the smallest particle sizes measured shows good agreement with the results of the particle penetration detector carried on the same spacecraft. Within the uncertainties of this experiment, a single size distribution seems appropriate for all but the smallest particles measured. There is evidence for the existence of a planetary sweeping effect in the vicinity of the Mars orbit. As calculated from the particle spatial distribution, the zodiacal light brightness is found to vary approximately as the inverse square of solar distance out to about 2.25 AU and then decrease more rapidly. The absolute value of the zodiacal light brightness as calculated from the particle spatial concentration is found to be too high by a factor of 10.

Soberman, R. K.↗

Analysis of hematite attrition in a grid jet apparatus

Particulate attrition is of interest for novel carbon-capture processes such as chemical looping combustion because the makeup cost of oxygen carrier is a significant portion of operating cost. As such, models to study and predict attrition of various oxygen carriers in fluidized bed systems are being developed. One of the regions of concern in fluidized bed systems is the high-velocity jet region near gas distributors in a fluid bed. This work studies the attrition of hematite particles using a modified ASTM apparatus to measure the particle size distribution throughout the experiment. Bed weight and gas velocity were varied. Hematite particles above the corresponding threshold value had a severe variation of particle size distributions which decreased with time. Weight fractions of the sieves over time were fit to a linear, time-variant population balance model to offer insight into particle attrition. The first-order rate constant was modified as a decaying exponential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particle Size Measurements from the 2018 Honeywell Uncertified Research Engine Test in the NASA Propulsion Systems Laboratory

The particle size measurements from the Honeywell Uncertified Research Engine Icing Test in the NASA Propulsion Systems Laboratory conducted in 2018 are presented. This work focuses on describing the experimental arrangement, the processing and analysis methods, and final results for select cases acquired during ice crystal cloud conditions. The measurements presented will include data acquired by two High Speed Imaging instruments located both upstream and downstream of the fan, in the engine bypass section. Experimental results demonstrate the expected downward shift in the particle size distribution as the cloud passes through the engine fan and the ice particles break-up. This work also demonstrates the feasibility of acquiring particle size distributions in the environment of a running engine.

Engine icing↗

Development of New Lunar Highland Regolith Simulant, NUW-LHT-5M

Introduction: NASA has a need for large quantities of lunar simulants that closely match future manned lunar missions at landing sites near the Lunar South Pole. The current simulant was designed to approximate NU-LHT-2M, and -4M, except using natural minerals and a fully synthetic, non-basaltic, high calcium glass. Washington Mills in Niagara Falls, NY was contracted to produce this new simulant due to their electric arc furnace (EAF) fusion technology and capabilities for crushing and sizing ceramics. The primary objectives of this work were: (1) to create a commercially produced simulant that contains a high fidelity, high calcium (An90+), low Mg/Fe glass without basaltic constituents; (2) to create a simulant that uses existing, terrestrially mined, and readily available minerals blended with a synthetic glass, and (3) accomplish objectives (1) and (2) using conventional processing methods capable of production levels that can meet substantial current and future demands of NASA projects. Another consideration was that highland simulants closest to meeting these requirements, NU-LHT-2M and -4M types, are out of production and largely out of stock. Design: The design chosen for development and production for this work is NU-LHT-2M, de-signed by Doug Stoeser and Douglas Rickman, which uses Stillwater anorthosite and norite rocks and synthetic glass. The composition targeted the average composition, glass content and particle size distribution of Apollo 16 samples [1]. The design utilizes two Stillwater minerals, anorthosite (37.7 wt.%), norite (17.6 wt.%), a commercially sourced olivine (4.7 wt.%), and a synthetic glass (40 wt.%). Processing: Stillwater rocks were hand collected and initially crushed by the USGS. Further crushing and milling was performed at Washing-ton Mills. For the glass, oxides were blended in the ap-propriate ratios, blended in a v-blender, and placed in a graphite lined, water cooled pot. The mixture was fused in Washington Mills’ pilot scale, 500 kW EAF. The molten glass was then poured into water and quenched to inhibit crystallization. Several pour/quench cycles were required to produce the total quantity of glass needed. Anorthosite, norite and olivine were milled together to the target particle size in a ball mill. The glass was ball milled separately. The complete simulant was produced by blending the milled materials. Particle size distribution, phase and chemical analyses of the materials was were per-formed at Washington Mills. Results: The glass composition was found to match the Apollo 16 oxide content average [1] quite well, with the exception of SiO2 and Fe2O3. The deviation of SiO2 and Fe2O3 from batched content was believed to be due to limitations inherent in the use of the arc furnace technique for these types of materials. A particle size distribution of the full simulant formulation closely matching the average for the Apollo 16 samples was also achieved. Fig. 1 shows the morphology of the final simulant particles. The size distribution is compared to the target in Fig. 2. Conclusions: A new lunar highlands type simulant containing a high-quality synthetic glass, closely matching an average composition of the Apollo 16 regolith was produced. Processing utilized easily scalable mineral and glass making methods.

simulant↗