Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity”

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 379 records · Page 21

Contaminant Leaching From Intact Saltstone Monoliths - 17517

At the Department of Energy's (DOE's) Savannah River Site (SRS) chemically reducing materials, such as blast furnace slag (BFS), are added to grout formulations mixed with low-level radioactive salt solution in order to enhance the attenuation of redox sensitive contaminants (e.g., technetium-99 (Tc-99)). The resulting cementitious material, known as saltstone, is deposited in a series of concrete vaults for permanent disposal at the Saltstone Disposal Facility (SDF). Under oxidizing conditions, Tc persists as an anion in the +7 oxidation state, i.e., pertechnetate (TcO{sub 4}{sup -}), with very limited retention generally associated with amphoteric soil oxides. However, Tc(VII) is subject to chemical reduction to the +4 oxidation state, which is less soluble and mobile in the environment. Thus, chemically reducing grouts provide both a physical (i.e., low saturated hydraulic conductivities (K{sub sat}) that limit H{sub 2}O turnover and O{sub 2} exposure) and chemical barrier (i.e., residual reductive capacity) to contaminant release. However, many of the previous experiments evaluating the ability of saltstone to chemically reduce and immobilize Tc have been conducted using ground saltstone materials as sorbents, a practice that is likely to reduce the moisture level in the material, enhance exposure to O{sub 2} and alter sorbent properties in an unpredictable manner. Therefore, the objective of the current study was to evaluate contaminant leaching from intact monoliths that better represent the initial physical and chemical state of saltstone within the SDF disposal units. For the current study, saltstone simulants were produced utilizing SRR prescribed formulations: initial batches were spiked with rhenium (Re) to serve as nonradioactive analog for Tc-99, and subsequent batches were spiked with Tc-99 for comparison. The relative concentrations of I, Re and Tc-99 in the salt-waste simulants were consistent with the average concentrations of I-129 and Tc-99 in the low-level feed waste. In addition, SRR retrieved intact monoliths from Saltstone Disposal Unit (SDU) Cell 2A in 2015 for testing. Contaminant mass transfer rates for Tc-99 and other contaminants from spiked saltstone simulants and actual SDF saltstone samples were assessed using EPA Method 1315, Mass Transfer Rates of Constituents in Monolithic or Compacted Granular Materials Using a Semi-Dynamic Tank Leaching Procedure. This method was recently adopted for evaluating contaminant leaching from intact monolithic materials. Given the importance of redox speciation in controlling the mobility of Tc-99, the EPA mass transfer tests for the spiked simulants were initially conducted using leachates in equilibrium with three different test atmospheres: (1) oxic, (2) anoxic, and (3) anoxic reducing. Nitrate (NO{sub 3}{sup -}) leaching was evaluated as a poorly retained waste constituent. The leaching behavior of cesium-137 (Cs-137), representing about 98% of the initial radioactivity of the saltstone, was also evaluated for the SDU-extracted samples. Results from Method 1315 were also compared to a novel Dynamic Leaching Method (DLM) in which a flexible-wall permeameter was used to achieve saturated leaching under an elevated hydraulic gradient in an effort to evaluate the persistence of reductive capacity and subsequent changes in contaminant partitioning within intact saltstone monoliths. The composition of the chemical leachates from both tests was thoroughly analyzed in an effort to identify potential solid phases controlling contaminant partitioning through geochemical modeling. For the EPA 1315 tests using spiked simulants, Re leaching rates (and other poorly sorbing contaminants like NO{sub 3}{sup -} and iodine), as indicated by high effective diffusivities (D{sub e}) and low leachability indices (i.e. LI = - log[D{sub e}]), were much higher than Tc-99, which was attributed to poor initial reduction and/or rapid oxidation of Re, indicating that Re is a poor chemical analog for evaluating Tc partitioning under reducing conditions. Technetium-99 leaching rates for the spiked saltstone samples also appeared to be sensitive to the estimated reduction capacity of the dry feed materials used in making grout, a characteristic that can vary between different BFS sources. Technetium-99 leaching rates for the intact saltstone samples in the EPA 1315 test were quite similar, with less variability between the three test samples. This may be indicative of the longer curing time before the SDF samples were tested (∼3 years). Somewhat surprisingly, the leaching rates for Cs-137 from the SDF samples were generally lower than the poorly sorbing contaminants, NO{sub 3}{sup -} and iodine. For the DLM tests, Re leaching from spiked saltstone was much higher than observed for Tc-99, with rates similar to that of non-reactive grout constituents, consistent with EPA 1315 leaching rates. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Solid Secondary Waste Form Testing Activities for the Integrated Disposal Facility Performance Assessment Support - 17567

Solid secondary waste (SSW) that will be generated at the Waste Treatment and Immobilization Plant (WTP) and from other activities at the Hanford site is destined for disposal at the Integrated Disposal Facility (IDF). Construction of the IDF was undertaken based on a Performance Assessment (PA) performed in 2001. In 2016, a SSW data package was prepared to support the development of an updated PA (2017 IDF PA). There are many inputs required to support the IDF PA. Data needs were identified based on experience from the previous PAs and from the development of conceptual models and modeling approaches to be applied for the 2017 PA. Interactions with the IDF PA team and the inventory data package developer helped to focus detailed efforts on a few key waste streams, contaminants of concern, and disposal configurations that are expected to be the primary contributors to the dose resulting from SSW disposed at the IDF. The key waste streams include: High Efficiency Particulate Air (HEPA) filters, ion exchange resins, granular activated carbon, and Ag-mordenite. The key parameters addressed in the 2016 data package include: saturated hydraulic conductivity, moisture characteristic curves, effective diffusion coefficient, distribution coefficient and porosity. The IDF PA team also identified a set of key contaminants that are expected to be the primary concern for the PA calculations (Tc-99, I-129, Cs-137, Sr-90, uranium isotopes (and total uranium), chromium, mercury and nitrate). These contaminants were the focus of the SSW data package. Specific formulations had not been identified for cementitious materials that will be used to encapsulate or stabilize SSW, and no IDF-specific experiments had been conducted to obtain the data for the PA. Thus, for the 2016 SSW data package, recommended property values were provided for a range of representative candidate materials based on a review of existing literature. However, starting in FY17, IDF-specific experiments will be performed for Hanford SSW expected to contribute most significantly to the dose for the IDF. This paper will summarize the information obtained for the 2016 data package and the plans for further testing to refine key inputs for the IDF PA. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling the Long-Term Behavior of Cementitious Barriers - 18188

Cementitious materials are commonly included as backfill/barrier materials in geological repositories for radioactive waste. An understanding of the long-term behavior of cement-based barriers is therefore required to provide confidence in their safety functions within an environmental safety case. Geochemical modelling has been undertaken on several different radioactive waste repository concepts (and their natural/industrial analogues), to explore the nature of cement degradation over timescales relevant to postclosure safety (many tens of thousands of years to c. 1 Ma). In radionuclide transport calculations undertaken as part of performance assessments, it is typically considered that the properties of cementitious engineered barrier materials change as they degrade, potentially resulting in an increase in porosity hydraulic conductivity, and effective diffusion coefficients. However, reactive transport modelling (supported by some analogue data) suggests that, in some repository systems, a degree of clogging of the cement pore space could occur and that some advective flow paths (cracks in the barriers, fractures in the host rock) could become (to some extent) infilled by cement and rock alteration products. This could result in a reduction in water flow and contaminant transport in the vicinity of the cement-rock interface, but some uncertainties remain (especially process couplings and the nature of water flow in complex fracture networks that are becoming clogged). In addition, the armoring of fracture surfaces could result in decreased potential for diffusion of solutes into the rock matrix. It is envisaged that ongoing work will continue to improve reactive-transport model capabilities, such that the implications cement degradation has for contaminant transport can be predicted with greater confidence. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Oxidation of Advanced Steel Cladding Alloys in Extreme Environments

Accident tolerant claddings must withstand the loss of coolant for longer times than conventional zirconium-based reactor claddings. Advanced steel alloys (such as Alloy 33, APMT) offer enhanced resistance to high-temperature oxidation in steam environments due to the formation of protective oxide layers of chromia and/or alumina oxides on the surface of the alloy, which enhance their stability at high temperatures by acting as an effective diffusion barrier, so resulting in low rate of oxidation. These passivating oxide layers, resistant to accelerated oxidation, offer an improvement over the current designs by reducing the production of hydrogen gas, and by delaying the release of fission products into the environment. Hence, these alloys present a great opportunity to explore alternative cladding materials with better accident tolerance for commercial light water reactors. Alloy 33 and APMT showed enhanced resistance to corrosion under the normal operation and off-nominal conditions of light water reactors.

36 MATERIALS SCIENCE↗

End Effects and Load Diffusion in Composite Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Specific problems recently considered were focussed on end effects in sandwich structures and for functionally graded materials. Both linear and nonlinear (geometric and material) problems have been addressed. Our goal is the development of readily applicable design formulas for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses. The decay behavior of stresses and other field quantities provides a significant aid towards this process. The analysis is also amenable to parameter study with a large parameter space and should be useful in structural tailoring studies.

Horgan, Cornelius O.↗

The effect of helium diffusion on the ages of globular clusters

Evolutionary tracks for halo stars were calculated, and isochrones which include the effects of microscopic diffusion of helium were constructed. The isochrones were fitted to a metal poor (M92) and a moderately metal rich (NGC 288) globular cluster using an updated version of the Revised Yale Isochrone color calibration. Ages of the two clusters were also determined using the difference between the turnoff magnitude and horizontal branch magnitude, and the difference in color between the main-sequence turnoff and lower giant branch. Considering all methods and constraints, diffusion is argued to reduce the derived ages of M92 and NGC 288 by 0.5-1 Gyr. The maximum age reduction that diffusion could cause is 3 Gyr. Age estimates including diffusion indicate that M92 is 16 +/- 2 Gyr old, and that M92 is about 3 Gyr older than NGC 288, assuming that the clusters have the same O/Fe of +0.4.

Chaboyer, Brian↗

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser-induced pressure-wave and barocaloric effect during flash diffusivity measurements

We report laser-induced pressure-wave and barocaloric effect captured by an infrared detector during thermal diffusivity measurements. Very fast (< 1 ms) and negative transients during laser flash measurements were captured by the infrared detector on thin, high thermal conductivity samples. Standard thermal diffusivity analysis only focuses the longer time scale thermal transient measured from the back surface due to thermal conduction. These negative spikes are filtered out and ignored as noise or anomaly from instrument. This study confirmed that the initial negative signal was indeed a temperature drop induced by the laser pulse. The laser pulse induced instantaneous volume expansion and the associated cooling in the specimen can be explained by the barocaloric effect. The initial cooling (< 100 microsecond) is also known as thermoelastic effect in which a negative temperature change is generated when the material is elastically deformed by volume expansion. A subsequent temperature oscillation in the sample was observed and only lasted about one millisecond. The pressure-wave induced thermal signal was systematically studied and analyzed. In conclusion, the underlying physics of photon-mechanical-thermal energy conversions and the potential of using this signal to study barocaloric effects in solids are discussed.

36 MATERIALS SCIENCE↗

Effects of Hydrodynamics of Soot Formation in Laminar Opposed-Jet Diffusion Flames

Effects of flow (hydrodynamic properties on the presence of soot in hydrocarbon-fueled laminar opposed-jet diffusion flames were studied experimentally at atmospheric pressure, emphasizing effects of velocities normal to the Rom sheet. These velocities were varied for conditions corresponding to combustion in air by transferring nitrogen from the oxidizer stream to the fuel stream which increases the stoichiometric mixture fraction of the flame and causes the stagnation plane of the flow to shift toward the fuel-rich side of the flame sheet. Fuels considered included acetylene, ethylene, ethane, propylene, propane, and 1-3 butadiene. Present measurements consisted of the critical strain rates for the flames to contain soot (the soot extinction limit) and for the flames to extinguish (the flame extinction limit) It was found that increasing the stoichiometric mixture fraction causes a progressive reduction of the critical strain rates for both flame and soot extinction; however, their ratio increases, and even becomes unbounded in most instances to yield a permanently-blue-flame regime. The results suggest that soot formation in nonpremixed flames can be controlled by varying velocities normal to the flame sheet. Nevertheless, definitive conclusions along these lines require evaluation of effects of corresponding variations of fuel and oxygen concentrations on soot formation when velocities normal to the flame sheet are changed by varying stoichiometric mixture fractions for laminar opposed-jet diffusion flames.

Lin, K.-C.↗

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗

NASA HECC Geometry and Performance Review Part 3: A Numerical and Experimental Investigation of Tip Clearance Effects on the Vaneless Diffuser Configuration

Tip clearance effects in centrifugal compressors have been extensively investigated to understand the losses associated with the flow in the impeller tip clearance gap between the rotating blades and the stationary shroud. In Part 3 of this multipart investigation, experimental data and validated numerical simulations from Parts 1 and 2 were used to analyze the effects of the size of the tip clearance gap on the High Efficiency Centrifugal Compressor performance and aerodynamics. Four tip clearance gaps ranging from 0.012-in to 0.030-in (2.0% to 4.9% of the impeller exit blade height) were considered at both design and off-design operating conditions. The total pressure ratio and efficiency of the stage are found to decay linearly with increasing tip gaps. The sensitivity of the impeller performance to the tip gap was found to vary with the rotational speed of the compressor. Spanwise surveys of flow angle, total pressure, and total temperature collected at the impeller exit at design speed are used to validate numerical simulations at each experimental tip gap condition. Numerical simulations show that increased turbulence kinetic energy near the shroud at larger tip gaps leads to spanwise mixing of high entropy fluid near the impeller trailing edge which decreases the useful work input by the impeller. The data presented have been made available to the public in the HECC Data Archive located at https://storage.googleapis.com/hecc-data/NASA-HECC-Data-Archive.zip.

centrifugal compressor↗

A numerical study of the effects of ambipolar diffusion on the collapse of magnetic gas clouds

The gravitational collapse of isothermal, nonrotating magnetic gas clouds have been calculated numerically, including the effects of ambipolar diffusion. The fractional ionization in the clouds is approximated by a power-law function of the gas density, f = K/n to the q-power, where K and q are adjustable parameters. Eleven numerical experiments were run, and the results indicate that the asymptotic character of collapse is determined mainly by the value of q and is largely independent of the other parameters characterizing a cloud (e.g., K, cloud mass). In particular, there is nearly a one-to-one correspondence between q and the slope, x, of the central magnetic field strength-gas density relationship. If q is no more than 0.8, a cloud collapses asymptotically, as though the magnetic field were 'frozen' to the neutral matter. The magnetic field strength at the center of a collapsing cloud is strongly amplified during collapse even for values of q of about 1, despite extremely low values of fractional ionization. A discussion of the theoretical basis for this unexpected behavior is given. Possible implications of our results for the problems of magnetic braking of rotating protostars and star formation in general are also presented.

Black, D. C.↗

The effect of velocity overshoot on the performance of magnetohydrodynamic subsonic diffusers

The effect of velocity overshoot, a velocity distribution encountered in magnetohydrodynamic channels with peak velocity near the wall instead of at the center, on the performance of two-dimensional subsonic diffusers was investigated. For the simplified conditions analyzed in this paper it was found that the ratio of peak-to-centerline velocity increases along the diffuser for a diffuser half angle greater than some critical value. It was also found that irrespective of the accompanying inlet temperature distribution the wall shear stress and the wall heat flux is substantially larger when the inlet velocity profile has an overshoot than that for a fully developed turbulent inlet velocity profile.

Greywall, M. S.↗

Evolutionary effects of helium diffusion in population II stars

Evolutionary calculations continuing until well past turnoff are presented for models of low-mass Population II stars which take into account the effects of the diffusion of helium relative to hydrogen. Evolutionary tracks, cluster isochrones and hydrogen distributions were obtained for stellar masses in the range 0.75 to 1.01 solar masses, both in the presence and absence of diffusion. It is found that for a star of a given mass, diffusion speeds up the evolutionary process on the main sequence, although after turnoff evolution is slowed with respect to the case without diffusion. As the stars ascend the red giant branch, their outer regions are remixed so that evidence of helium diffusion is erased, and the evolutionary tracks of the models with and without diffusion converge. Thus, if the age of a globular cluster is determined from the absolute magnitude at turnoff or from fitting isochrones, diffusion results in a 25% reduction in the derived age at a turnoff magnitude of 4.23, and a 14% reduction at a turnoff magnitude of 3.45.

Stringfellow, G. S.↗

Proton irradiation effects on minority carrier diffusion length and defect introduction in homoepitaxial and heteroepitaxial n-GaN [Proton irradiation effects on minority carrier diffusion length and defect introduction in homoepitaxial n-GaN]

Here, inherent advantages of wide bandgap materials make GaN-based devices attractive for power electronics and applications in radiation environments. Recent advances in the availability of wafer-scale, bulk GaN substrates have enabled the production of high quality, low defect density GaN devices, but fundamental studies of carrier transport and radiation hardness in such devices are lacking. Here, we report measurements of the hole diffusion length in low threading dislocation density (TDD), homoepitaxial n-GaN, and high TDD heteroepitaxial n-GaN Schottky diodes before and after irradiation with 2.5 MeV protons at fluences of 4–6 × 10 13 protons/cm 2 . We also characterize the specimens before and after irradiation using electron beam-induced-current (EBIC) imaging, cathodoluminescence, deep level optical spectroscopy (DLOS), steady-state photocapacitance, and lighted capacitance-voltage (LCV) techniques. We observe a substantial reduction in the hole diffusion length following irradiation (50%–55%) and the introduction of electrically active defects which could be attributed to gallium vacancies and associated complexes (V Ga -related), carbon impurities (C-related), and gallium interstitials (Ga i ). EBIC imaging suggests long-range migration and clustering of radiation-induced point defects over distances of ~500 nm, which suggests mobile Ga i . Following irradiation, DLOS and LCV reveal the introduction of a prominent optical energy level at 1.9 eV below the conduction band edge, consistent with the introduction of Ga i .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗

Effects of turbulent diffusion and back-reaction on the dust distribution around two resonant planets

In evolved and dusty circumstellar discs, two planets with masses comparable to Jupiter and Saturn that migrate outwards while maintaining an orbital resonance can produce distinctive features in the dust distribution. Dust accumulates at the outer edge of the common gas gap, which behaves as a dust trap, where the local dust concentration is significantly enhanced by the planets’ outward motion. Concurrently, an expanding cavity forms in the dust distribution inside the planets’ orbits, because dust does not filter through the common gaseous gap and grain depletion in the region continues via inward drifting. There is no cavity in the gas distribution because gas can filter through the gap, although ongoing gas accretion on the planets can reduce the gas density in the inner disc. Such behaviour was demonstrated by means of simulations neglecting the effects of dust diffusion due to turbulence and of dust backreaction on the gas. Both effects may alter the formation of the dust peak at the gap outer edge and of the inner dust cavity, by letting grains filter through the dust trap. We performed high-resolution hydrodynamical simulations of the coupled evolution of gas and dust species, the latter treated as pressureless fluids, in the presence of two giant planets. We show that diffusion and backreaction can change some morphological aspects of the dust distribution but do not alter some main features, such as the outer peak and the expanding inner cavity. Furthermore, these findings are confirmed for different parametrizations of gas viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗