Search NASA⌕ Search

SEARCH · Search NASA

Results for “low viscosity”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nonmonotonic Strain Rate Dependence on the Strain Hardening of Polymer Nanocomposites

This Letter investigates the external deformation on modifying the polymer–nanoparticle (NP) and NP–NP interactions as well as their influences on the macroscopic properties of polymer nanocomposites (PNCs). Specifically, the applied uniaxial extension brings together the NPs along the transverse stretching direction and catalyzes the percolation transition from the initially well-dispersed NPs in the PNCs. The percolated NP network thus increases the strength of PNCs and leads to a strong surge in the elongation viscosity at very low strain rates, exhibiting the unexpected nonmonotonic strain rate dependence on the strain hardening. These results reveal explicitly the unconventional roles of the external deformation on modulating the relative strength of the polymer–NP and the NP–NP interactions, enabling advanced structures and properties control of PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Improved single-particle phase-space distributions for viscous fluid dynamic models of relativistic heavy ion collisions

Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Present and Future of QCD: QCD Town Meeting White Paper – An Input to the 2023 NSAC Long Range Plan

It is currently understood that there are four fundamental forces in nature: gravitational, electromagnetic, weak and strong forces. The strong force governs the interactions between quarks and gluons, elementary particles whose interactions give rise to the vast majority of visible mass in the universe. The mathematical description of the strong force is provided by the non-Abelian gauge theory Quantum Chromodynamics (QCD). While QCD is an exquisite theory, constructing the nucleons and nuclei from quarks, and furthermore explaining the behavior of quarks and gluons at all energies, remain to be complex and challenging problems. Such challenges, along with the desire to understand all visible matter at the most fundamental level, position the study of QCD as a central thrust of research in nuclear science. Experimental insight into the strong force can be gained using large particle accelerator facilities, which are necessary to probe the very short distance scales over which quarks and gluons interact. The Long Range Plans (LRPs) exercise of 1989 and 1996 led directly to the construction of two world-class facilities: the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) that is focused on studying how the structure of hadrons emerges from QCD (cold QCD research), and the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL) that aims at the discovery and study of a new state of matter, the quark-gluon plasma (QGP), at extremely high temperatures (hot QCD research). These past investments have produced major advances. Nucleons and nuclei are being studied with increasing precision with a unified description of the partonic structure utilizing multi-dimensional imaging. Significant progress has been made, paving the way towards a complete picture of how quarks and gluons give rise to the mass, spin, and momentum of the nucleon. In hot QCD, the QGP is created in the collisions of nuclei at RHIC and the Large Hadron Collider (LHC) and is observed to behave like a fluid with very low specific shear viscosity; the current goals are to understand how the fluid behavior emerges from QCD and to characterize the temperature (and chemical potential) dependence of the properties of the QGP. As this White Paper is written, current experimental programs at CEBAF, RHIC and the LHC continue to provide exciting near term opportunities to capitalize on the investments in experimental equipment and accelerator operations. Most importantly, the QCD community looks forward to the construction of the Electron Ion Collider (EIC) as a major new facility to push forward QCD research in the next decades, with significant focus on exploring the properties of gluons, the mediators of the strong force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Modeling of High Solids Biomass Liquefaction, Techno-economic Assessment and Life Cycle Assessment for Industrial Scale UP

Biomass liquefaction improves corn-stover slurry transport within a biorefinery, and reduces challenges associated with solids handling. Slurries, formed with yield stresses (i.e., viscosities) that are low enough to facilitate pumping from one unit operation to the next, and that enable mixing at reasonable power inputs in large fermentation vessels, enable scaling of cellulose fermentations. Three liquefaction methods: enzyme, enzyme mimetic, and a hybrid approach were analyzed by a combination of techno-economic and life-cycle analyses, resulting in an assessment of cost-effectiveness and environmental impact. Of the three scenarios, enzyme-based processes emerged as the most efficient, both in cost and emissions reduction (gCO2-equivalent/kg slurry).

Ramirez, Jorge↗
DOE OSTI2023

Interpreting T-History Data on Supercooling and Phase Change Enthalpy of a Salt-Hydrogel Complex for Building Thermal Management

Phase change materials (PCMs) are promising materials for storing and discharging latent heat in buildings to regulate the thermal environment. Measurement of salt hydrate PCM thermal performance is often challenging because of supercooling and phase segregation. Though differential scanning calorimetry (DSC) has been widely used to measure thermal properties of various phase change materials, DSC uses milligrams of sample and is unrealistic for predicting large-scale materials performance. The T-history method [1], in contrast, handles higher volumes but the analysis typically requires a lumped capacitance assumption. Here, we report and analyze T-History data on a novel salt-hydrogel complex with high viscosity, that renders low Biot number sample preparation extremely challenging. We re-examine the standard T-history method and compared against DSC data in this complex material. We report details of an experimental setup, built with necessary additions to a controlled environmental chamber. We carefully examine the validity of the lumped capacitance assumption and suggest a modified analysis using a computational approach for reliable data. Experimental parameters such as ramp rate and temperature range of the environmental chamber are considered via computational modeling to understand the effect of ramp rate on supercooling [2]. Finally, we report thermal cycling experiments on the salt-hydrogel complex. Our modified approach to measuring supercooling and enthalpy of fusion at large scales is important in better understanding the performance of phase-change materials at scale for building thermal storage. References: [1] Marin, Jose M, Belen Zalba, Luisa F Cabeza, and Harald Mehling. “Determination of Enthalpy Temperature Curves of Phase Change Materials with the Temperature-History Method: Improvement to Temperature Dependent Properties.” Measurement Science and Technology 14, no. 2 (February 1, 2003): 184–89. https://doi.org/10.1088/0957-0233/14/2/305. [2] Safari, A., R. Saidur, F.A. Sulaiman, Yan Xu, and Joe Dong. “A Review on Supercooling of Phase Change Materials in Thermal Energy Storage Systems.” Renewable and Sustainable Energy Reviews 70 (April 2017): 905–19. https://doi.org/10.1016/j.rser.2016.11.272

thermodynamics↗
DOE OSTI2024

Effect of glass forming additives on low-activity waste feed conversion to glass

A significant effort was invested in the past to develop and refine mathematical models that relate the composition of nuclear waste glasses with their properties, such as viscosity, electrical conductivity, or chemical durability. However, less attention has been paid to the formulation of the melter feed itself, such as the chemical form and the particle size of the glass forming and modifying additives (GFMA), which have a significant effect on the feed-to-glass conversion process during melting. To address this issue, we systematically changed the mineral composition of a simulated low-activity waste melter feed and inspected its melting behavior. When substituting minerals with corresponding oxides and hydroxides, we found that different alumina sources (kyanite, gibbsite, boehmite, or corundum) had the strongest effect on the feed melting process, whereas the sources of Ca, Mg, and Zr had little effect. Further, the x-ray diffraction analysis showed that the alumina sources differ in their dissolution kinetics: early dissolving alumina sources, such as gibbsite (Al(OH) 3 ) and boehmite (AlO(OH)), increase the transient glass-forming melt viscosity at early stages, when gases still evolve, causing extended foaming, whereas alumina sources that dissolve at high temperatures, such as kyanite (Al 2 SiO 5 ) and corundum (Al 2 O 3 ), keep the transient glass-forming melt viscosity low and lead to a faster foam collapse. Using viscosity-composition relationship to estimate the viscosity of transient glass-forming melts in the primary foaming range, we found that the primary foam began to collapse at 360 to 800 Pa s, and fully collapsed between 65 and 260 Pa s. This result agrees with our previous studies, according to which the glass-forming melt viscosity at the cold cap foam bottom ranged from 24 to 85 Pa s.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Predicted Melt Curve and Liquid Shear Viscosity of RDX up to 30 GPa

Recent grain-scale simulations of HMX and TATB have shown that predictions for hot spot formation in high explosives are particularly sensitive to accurate determinations of the pressure-dependent melt curve and the shear viscosity of the liquid phase. These physics terms are poorly constrained beyond ambient pressure for the explosive RDX. Here we adopt an all-atom modeling approach using molecular dynamics (MD) simulations to predict the melt curve of RDX near to detonation conditions (30 GPa) and determine the shear viscosity of the liquid as a function of temperature and pressure above the melt curve. Phase-coexistence simulations were used to determine the melt curve, which is predicted to vary by almost 1100 K as the pressure increases from 0 GPa to 30 GPa. Equilibrium MD simulations and the Green-Kubo formalism were used to obtain the pressure-temperature-dependent shear viscosity. The shear viscosity of RDX is predicted to be of similar magnitude to the viscosity of TATB at low GPa-range pressures, and to be roughly an order of magnitude lower than the viscosity of HMX. The temperature dependence of the shear viscosity is Arrhenius at a given pressure, and the exponential pre-factor and activation term exhibit a strong, yet complicated, pressure dependence. An empirical pressure-temperature dependent function for RDX shear viscosity is developed that simultaneously captures a wide range of MD predictions while taking an analytic form that extrapolates smoothly beyond the fitted regime. The relative strength of the pressure and temperature dependencies of these two physics terms is found to be of similar magnitude for RDX, HMX, and TATB, which motivates incorporating these results in future RDX grain-scale modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport Phenomena in Low Temperature Lithium-Ion Battery Electrolytes

Lithium-ion batteries face low temperature performance issues, limiting the adoption of technologies ranging from electric vehicles to stationary grid storage. This problem is thought to be exacerbated by slow transport within the electrolyte, which in turn may be influenced by ion association, solvent viscosity, and cation transference number. How these factors collectively impact low temperature transport phenomena, however, remains poorly understood. Here we show using all-atom classical molecular dynamics (MD) simulations that the dominant factor influencing low temperature transport in LP57 (1 M LiPF 6 in 3:7 ethylene carbonate (EC)/ethyl methyl carbonate (EMC)) is solvent viscosity, rather than ion aggregation or cation transference number. We find that ion association decreases with decreasing temperature, while the cation transference number is positive and roughly independent of temperature. In an effort to improve low temperature performance, we introduce γ -butyrolactone (GBL) as a low viscosity co-solvent to explore two alternative formulations: 1 M LiPF 6 in 15:15:70 EC/GBL/EMC and 3:7 GBL/EMC. While GBL reduces solution viscosity, its low dielectric constant results in increased ion pairing, yielding neither improved bulk ionic conductivity nor appreciably altered ion transport mechanisms. We expect that these results will enhance understanding of low temperature transport and inform the development of superior electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2022

Optimization of matrix viscosity improves polypropylene fiber dispersion and properties of engineered cementitious composites

Engineered cementitious composites (ECC) is a durable cementitious material with high tensile ductility and strain-hardening characteristics. Although considered as a cost-effective fiber for ECC, polypropylene (PP) fiber is reportedly difficult to disperse in mortar matrix due to its high aspect ratio and hydrophobicity. In this study, the matrix viscosity of a low-carbon ECC based on limestone calcined clay cement was tailored as a variable to improve PP fiber dispersion under a pre-determined mixing protocol. The effect of matrix viscosity on the composite fresh and hardened properties was investigated experimentally. Results suggested an optimal range of matrix viscosity (10.3–11.5 Pa . s) favors the composite tensile strength and strain capacity at 28 days. At the optimal state with a 0.1 % viscosity modifying admixture (VMA)-to-binder mass ratio, PP-ECC achieved 7.0 % tensile strain capacity and 3.5 MPa ultimate tensile strength. When matrix viscosity falls outside the desired range, both ultimate tensile strength and strain capacity were diminished. By tailoring the VMA dosage, the matrix viscosity can be adjusted for desired fiber dispersion, workability, and mechanical properties. Finally, the findings of this study provide a technical reference for the practical design and application of PP-ECC.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting

In this paper, we develop high-order nodal discontinuous Galerkin (DG) methods for hyperbolic conservation laws that satisfy invariant domain preserving properties using subcell flux corrections and convex limiting. These methods are based on a subcell flux corrected transport (FCT) methodology that involves blending a high-order target scheme with a robust, low-order invariant domain preserving method that is obtained using a graph viscosity technique. Furthermore, the new low-order discretizations are based on sparse stencils which do not increase with the polynomial degree of the high-order DG method. As a result, the accuracy of the low-order method does not degrade when used with high-order target methods. The method is applied to both scalar conservation laws, for which the discrete maximum principle is naturally enforced, and to systems of conservation laws such as the Euler equations, for which positivity of density and a minimum principle for specific entropy are enforced. Numerical results are presented on a number of benchmark test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Enhancing membrane performance for CO 2 capture from flue gas with ultrahigh MW polyvinylamine

Membranes for post-combustion CO 2 capture are required to have a high CO 2 permeance due to the limited driving force. For facilitated transport membranes synthesized with polyvinylamine (PVAm), a defect-free selective layer of <200 nm is usually required to render sufficient permeance with high CO 2 /N 2 selectivity. In order to meet such a demand through the knife-coating process, the coating solution needs to have a high viscosity at a relatively low concentration to minimize its penetration into the substrate. The demand was met by synthesizing PVAm with an ultrahigh molecular weight (MW) via inverse emulsion polymerization (IEP) in this study. Compared to solution polymerization, IEP isolates the reaction in inverse micelles suspended in a continuous organic phase, which allows excellent dissipation of the heat generated by the reaction and reduces gel formation drastically. The polymerization parameters, including monomer concentration, initiator concentration, and reaction temperature, were investigated to obtain the optimal MW of PVAm for membrane performance. As compared to solution polymerization, IEP enhanced the MW of PVAm from 1.2 to 12.7 MDa, which minimized the penetration of the coating solution into the substrate and hence the extra mass transfer resistance. The effect of MW and degree of hydrolysis of PVAm on the membrane transport properties were studied. Furthermore, by employing PVAm with a MW of 12.7 MDa to strengthen the polymer matrix, the loading of piperazine glycinate in the membrane was increased up to 85 wt.%. In conclusion, the resultant membrane achieved a CO 2 permeance of 839 GPU and a CO 2 /N 2 selectivity of 161 at the typical flue gas temperature of 57°C.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Machine learning-enabled discovery of ionic liquid–solvent electrolytes exhibiting high ionic conductivity

Ionic liquids (ILs), which are a class of materials with versatile nature and growing popularity, are facing impediments toward widespread usage as electrolytes due to various factors such as low ionic conductivity, high viscosity, high market price etc. One of the ways these limitations can be addressed is by mixing ILs with a molecular solvent. In a combinatorial sense, there exists an immense number of specific IL–solvent combinations. An exhaustive experimental or even simulation-based investigation of the chemical space spanned by such combinations can be extremely time-consuming, expensive, and nearly impossible. An alternative approach is to employ machine learning-based models developed from available databases. Although there exists prior literature that integrates machine learning to investigate mixtures of specific solvents with ILs, these models lack generalization necessitating development of a large number of ML models to handle various solvents. To remedy this shortcoming, as a part of designing green electrolytes with high ionic conductivity that can have potential applications in next-generation batteries and solar cells, this work aims to develop a unified machine learning model to predict ionic conductivity of any IL–solvent mixture system. In this regard, three models, namely, Random Forest, extreme gradient boosting (XGBoost), and artificial neural network (ANN) were formulated using the NIST ILThermo database. The dataset contained 549 unique ionic liquids from 16 cation families and 81 unique solvents, representing a total of 23 712 datapoints. SHAPLEY additive explanation (SHAP) method was used to assess the impact of various features on model prediction and their significance was compared with literature to gain physical insight about the model behavior. Finally, using the developed models, approximately 2.5 million IL–solvent mixtures at five different compositions were screened at room temperature. The high-throughput screening yielded nearly 19 000 IL–solvent mixtures for which ionic conductivity was found to exceed the ionic conductivity of conventional Li-ion battery electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2024

Laser‐Sintered Silver Metallization for Silicon Heterojunction Photovoltaic Cells

Herein, a novel metallization technique is reported for crystalline silicon heterojunction (SHJ) solar cells in which silver (Ag) fingers are printed on the SHJ substrates by dispensing Ag nanoparticle‐based inks through a needle and then sintered with a continuous‐wave carbon dioxide (CO 2 ) laser. The impact of the Ag ink viscosity on the line quality and the line resistance is investigated on three Ag inks with different viscosities. Increasing ink viscosity yields higher Ag contact heights, larger aspect ratios, and lower line resistance values. The Ag line height increases from less than a micrometer to ≈18.62 ± 3.48 μm with the increasing viscosity. Photoluminescence imaging shows that the low‐resistance Ag metal contacts obtained do not result in any passivation damage of the SHJ substrate. This is because the wavelength of light emitted from the CO 2 laser (i.e., 10.6 μm) leads to optical absorption in the Ag, but this light is effectively transparent to the transparent conductive oxide film, amorphous silicon films, and crystalline silicon substrate. Bulk resistivity values as low as 6.5 μΩ cm are obtained for the laser‐sintered Ag contact and printed using the Ag ink with the highest viscosity in this work.

Mousumi, Jannatul Ferdous↗
DOE OSTI2020

Viscous Properties of a Degenerate One-Dimensional Fermi Gas

In this paper, we study the viscous properties of a system of weakly interacting spin-1/2 fermions in one dimension. Accounting for the effect of interactions on the quasiparticle energy spectrum, we obtain the bulk viscosity of this system at low temperatures. Our result is valid for frequencies that are small compared with the rate of fermion backscattering. For frequencies larger than this exponentially small rate, the excitations of the system become decoupled from the center of mass motion, and the fluid is described by two-fluid hydrodynamics. We calculate the three transport coefficients required to describe viscous dissipation in this regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Global Simulations of Phase State and Equilibration Time Scales of Secondary Organic Aerosols with GEOS-Chem

The phase state of secondary organic aerosols (SOA) can range from liquid through amorphous semisolid to glassy solid, which is important to consider as it influences various multiphase processes including SOA formation and partitioning, multiphase chemistry, and cloud activation. In this study, we simulate the glass transition temperature and viscosity of SOA over the globe using the global chemical transport model, GEOS-Chem. The simulated spatial distributions show that SOA at the surface exist as liquid over equatorial regions and oceans, semisolid in the midlatitude continental regions, and glassy solid over lands with low relative humidity. The predicted SOA viscosities are mostly consistent with the available measurements. In the free troposphere, SOA particles are mostly predicted to be semisolid at 850 hPa and glassy solid at 500 hPa, except over tropical regions including Amazonia, where SOA are predicted to be low viscous. Phase state also exhibits seasonal variation with a higher frequency of semisolid and solid particles in winter compared to warmer seasons. We calculate equilibration time scales of SOA partitioning (τ eq ) and effective mass accommodation coefficient (α eff ), indicating that τ eq is shorter than the chemical time step of GEOS-Chem of 20 min and αeff is close to unity for most locations at the surface level, supporting the application of equilibrium SOA partitioning. However, τ eq is prolonged and α eff is lowered over drylands and most regions in the upper troposphere, suggesting that kinetically limited growth would need to be considered for these regions in future large-scale model studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Stability, Combustion, and Compatibility of High-Viscosity Heavy Fuel Oil Blends with a Fast Pyrolysis Bio-Oil

Properties related to the combustion, stability, and compatibility of blends composed of high-viscosity heavy fuel oil (HFO) and highly acidic pyrolysis bio-oil were determined to assess the utility of bio-oil as a marine fuel. The addition of bio-oil was shown to be fully stable with HFO at blend levels up to 50 mass % for up to 2 weeks. Bio-oil concentrations as low as 5 mass % significantly reduced the viscosity of HFO at 25 and 50 °C. Aging studies at 50 and 90 °C showed that the HFO inhibited the polymerization of bio-oil. The heating value and lubricity showed a linear dependency with bio-oil content, and combustion quality was acceptable for blends containing up to 15% bio-oil. The highly acidic bio-oil was found to be corrosive to carbon steel, 2.25Cr-1Mo steel, and 409 stainless steels, but not 304L and 316L. When blended into HFO at levels less than 19 mass %, no measurable corrosion was observed on any of the steel materials, but a 50 mass % concentration produced low-to-moderate corrosion in the carbon steel, 2.25Cr-1Mo steel, and 409 stainless steel grades. The combination of good blend stability, polymerization inhibition, reduced viscosity, and acceptable compatibility for low blend levels suggests that bio-oils may be suitable for use as a marine fuel.

09 BIOMASS FUELS↗
DOE OSTI2025

Bulk viscosity of nuclear matter with pions in the neutrino-trapped regime

Recent work has shown that a population of thermal pions could modify the equation of state and transport properties of hot and dense neutron-rich matter and introduce new reaction pathways to change the proton fraction. In this article we study their impact on the bulk viscosity of dense matter, focusing on the neutrino-trapped regime that would be realized in neutron star mergers and supernovae. We find that the presence of a thermal population of pions alters the bulk viscosity by modifying the EoS (via the susceptibilities) and by providing new reaction pathways to achieve beta equilibrium. In neutron star merger conditions, the bulk viscosity in neutrino-trapped npeμ matter (without pions) has its peak at temperatures of at most a couple MeV and is quite small at temperatures of tens of MeV. We find that thermal pions enhance the low-temperature peak of the bulk viscosity by a factor of a few and shift it to slightly lower temperatures. At higher temperatures, where the pion abundance is large but the bulk viscosity is traditionally small, pions can increase the bulk viscosity by an order of magnitude or more, although it is still orders of magnitude smaller than its peak value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗