Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,405 records · Page 78

Compressing turbulence and sudden viscous dissipation with compression-dependent ionization state

Turbulent plasma flow, amplified by rapid three-dimensional compression, can be suddenly dissipated under continuing compression. Furthermore, this effect relies on the sensitivity of the plasma viscosity to the temperature, μ ~ T 5 / 2 . The plasma viscosity is also sensitive to the plasma ionization state. Here, we show that the sudden dissipation phenomenon may be prevented when the plasma ionization state increases during compression, and we demonstrate the regime of net viscosity dependence on compression where sudden dissipation is guaranteed. In addition, it is shown that, compared to cases with no ionization, ionization during compression is associated with larger increases in turbulent energy and can make the difference between growing and decreasing turbulent energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-dimensional hydrodynamic viscous electron flow in annular Corbino rings

The concept of fluidic viscosity is ubiquitous in condensed-matter systems hosting a continuum where macroscopic properties can emerge. While an important property of liquids and some solids, only recently was the viscosity of an electron shown to play a role in electronic transport experiments. In this Letter, we present nonlocal electronic transport measurements in concentric annular rings formed in high-mobility two-dimensional electron gases, and the resulting data show that viscous hydrodynamic flow can occur far away from the source-drain current region. Our conclusion of viscous electronic transport is further corroborated by simulations of the Navier-Stokes equations that are found to be in agreement with our measurements below T = 1 K . Finally, this work emphasizes the key role played by viscosity via electron-electron ( e − e ) interaction even when the electronic transport is restricted radially, and for which it should have played no major role. Published by the American Physical Society 2025

Vijayakrishnan, Sujatha (ORCID:0009000080933182)↗

Design of a High Toughness Epoxy for Superconducting Magnets and Its Key Properties

Nb 3 Sn accelerator magnets are poised to enable the luminosity upgrade of the Large Hadron Collider (LHC) at CERN, improving its potential for exploring physics beyond the standard model of particle physics. The prototype Nb 3 Sn magnets consistently need 10–25 quenches to achieve their best performance. A hypothesis is that the long training of these magnets may at least be partially induced by epoxy cracking and bonding failures. In 2018, we showed that several existing epoxies have a higher toughness and less tendency to crack at low temperatures than CTD-101K, the epoxy resin with which almost all Nb 3 Sn accelerator magnets have been impregnated. Here we explore a new high toughness formulation for reducing quench training of Nb 3 Sn accelerator magnets, through combining two amine curing agents to achieve a good glass transition temperature (T g ), adding a viscosity reducer to achieve low viscosity, and experimenting a coupling agent for improved bonding strengths. In this work, we report results of comprehensive materials tests, including thermal shock, tensile, compressive, shear tests, viscosity and T g tests.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dispersion and properties of zirconia suspensions for stereolithography

Abstract Stereolithography is an attractive technique for the fabrication of complex‐shaped ceramic components with high dimensional accuracy. One of the challenges in this technology is the development of high solid loading, low viscosity photosensitive ceramic suspension. In this study, the dispersion of zirconia in photocurable resin and the slurry properties were intensively investigated. Rheological measurements showed that DISPERBYK‐103 proved to be an effective dispersant. 42 vol% ZrO 2 suspension was successfully prepared using 3.5 wt% DISPERBYK‐103 as the dispersant, with a suitable viscosity (4.88 Pa·s) below the maximum allowable viscosity value (5 Pa·s) for stereolithography applications. The adsorption behavior of DISPERBYK‐103 on the surface of zirconia powders was characterized by TG and FT‐IR, confirming the dispersion effect of dispersant. Contact angle measurements were also conducted to show that the adsorption of DISPERBYK‐103 could help to improve the wettability between powder and photocurable resin. Results showed that DISPERBYK‐103 was effective for the preparation of suitable slurries for the development of ZrO 2 ceramics through stereolithography.

Li, Xingbang↗

Heat transfer from glass melt to cold cap: Melting rate correlation equation

Increasing the efficiency of all-electric melters is a hot topic not only in the waste glass melting industry, where the melting rate directly influences the life cycle of the cleanup process, but also in the commercial glass industry, where increasingly stringent environmental regulations and rapid development of renewable energy sources ask for transformation of traditional fossil fuel technology. Batch-to-glass conversion in an all-electric melter proceeds in the cold cap—a layer of reacting and melting material floating on molten glass. The heat transfer into the cold cap is controlled by convection and conduction in the thermal boundary layer on the melt side, and by the properties of the foam layer at the cold cap/melt boundary. An overview of factors affecting heat transfer is presented and assessed using data from laboratory and pilot-scale experiments. The heat flux can be improved by (i) decreasing glass melt viscosity, (ii) decreasing the temperature at which the foam starts to collapse at the cold cap bottom, and (iii) increasing melt convection by bubbling gas into the glass melt under the cold cap. Detailed understanding of the conversion process can guide the formulation of feeds that melt easily and quickly.

Hrma, Pavel↗

X-ray tomography of feed-to-glass transition of simulated borosilicate waste glasses

The feed composition of a high level nuclear waste (HLW) glass melter affects the overall melting rate by influencing the chemical, thermophysical, and morphological properties of a relatively insulating cold cap layer over the molten phase where the primary feed vitrification reactions occur. Data from X ray computed tomography imaging of melting pellets comprised of a simulated high-aluminum HLW feed heated at a rate of 10°C/min reveal the distribution and morphology of bubbles, collectively known as primary foam, within this layer for various SiO 2 /(Li 2 CO 3 +H 3 BO 3 +Na 2 CO 3 ) mass fractions at temperatures between 600°C and 1040°C. To track melting dynamics, cross-sections obtained through the central profile of the pellet were digitally segmented into primary foam and a condensed phase. Pellet dimensions were extracted using Photoshop CS6 tools while the DREAM.3D software package was used to calculate pellet profile area, average and maximum bubble areas, and two-dimensional void fraction. The measured linear increase in the pellet area expansion rates – and therefore the increase in batch gas evolution rates – with SiO 2 /(Li 2 CO 3 +H 3 BO 3 +Na 2 CO 3 ) mass fraction despite an exponential increase in viscosity of the final waste glass at 1050°C and a lower total amount of gas-evolving species suggest that the retention of primary foam with large average bubble size at higher temperatures results from faster reaction kinetics rather than increased viscosity. However, viscosity does affect the initial foam collapse temperature by supporting the growth of larger bubbles. Because the maximum bubble size is limited by the pellet dimensions, larger scale studies are needed to understand primary foam morphology at high temperatures. This temperature-dependent morphological data can be used in future investigations to synthetically generate cold cap structures for use in models of heat transfer within a HLW glass melter.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Evaluating Transfer and Pumping of Slurries from Pulsed Jet Mixed Vessels

Due to gravity, solids in slurries will settle if density differences between the solids and liquid are positive (i.e., particle has a negative buoyant force) unless rheological properties and flow conditions are adequate to overcome the gravitational effects. The rate of settling depends on the force balance of the particle, which includes the surface forces associated with fluid rheology. Given the same fluid and solid properties, particles that are larger and denser tend to settle faster. When pumping slurry into a vessel at concentrations precluding hindered settling with insufficient mixing, particle and density distributions can result in preferential settling, creating stratification in the solids concentration within the vessel. For vessels with transfer line inlets located in the lower portion of the tank, the stratified solids concentration may be detrimental to the transfer system performance. Elevated concentrations of solids in the slurry entrained at the inlet to the transfer line can result in the effective viscosity or slurry bulk density exceeding the design limits of the pump. These conditions could result in plugging of the transfer line or onset of cavitation of the pumps because of excessive pressure drop. These conditions can be exacerbated with periodic inlet conditions existing at the transfer line inlet. Periodic conditions can result when vessel mixing is intermittent such as with pulsed jet mixers (PJM). The transfer line inlet conditions are impacted by the periodic nature of the PJM operations with respect to suspension of solids and their transport to the inlet of the transfer line. A scaling approach is presented, and corresponding test requirements are developed for assessing the prevention of plugging the pipeline. Line plugging mechanisms are addressed that do not include plugging due to steady-state high-density slurry entering the transfer line and reducing the net positive suction head available (NPSHA) at the pump inlet to below that required for pump operation. Items considered include the transition to reduced relative flow velocities, such that the critical pipe velocity for solids deposition, Ucd, is not maintained, and segregation of heavy solids during the transport. The recommended requirements to prevent plugging include: • Limits for viscosity and density for entrained slurry to prevent the pressure drop in the pipeline from exceeding pump capacity. • Limits for viscosity and density for entrained slurry to prevent the net positive suction head available (NPSHA) from falling below the net positive suction head required (NPSHR) for operating the pump. • Transfer line velocity and flow rate requirements to maintain solids in suspension, while avoiding line plugging that results from deposition of solids within the transfer line. This paper describes the development of the scaling and testing requirements to verify that proposed approaches for transfer and pump out are appropriately developed

periodic flow, mixing, pipeline plugging, scaling,↗

Optical CT imaging of solid radiochromic dosimeters in mismatched refractive index solutions using a scanning laser and large area detector

Purpose: The practical use of the PRESAGE® solid plastic dosimeter is limited by the inconvenience of immersing it in high-viscosity oils to achieve refractive index matching for optical computed tomography (CT) scanning. The oils are slow to mix and difficult to clean from surfaces, and the dosimeter rotation can generate dynamic Schlieren inhomogeneity patterns in the reference liquid, limiting the rotational and overall scan speed. Therefore, it would be beneficial if lower-viscosity, water-based solutions with slightly unmatched refractive index could be used instead. The purpose of this work is to demonstrate the feasibility of allowing mismatched conditions when using a scanning laser system with a large acceptance angle detector. A fiducial-based ray path measurement technique is combined with an iterative CT reconstruction algorithm to reconstruct images. Methods: A water based surrounding liquid with a low viscosity was selected for imaging PRESAGE® solid dosimeters. Liquid selection was optimized to achieve as high a refractive index as possible while avoiding rotation-induced Schlieren effects. This led to a refractive index mismatch of 6% between liquid and dosimeters. Optical CT scans were performed with a fan-beam scanning-laser optical CT system with a large area detector to capture most of the refracted rays. A fiducial marker placed on the wall of a cylindrical sample occludes a given light ray twice. With knowledge of the rotation angle and the radius of the cylindrical object, the actual internal path of each ray through the dosimeter can be calculated. Scans were performed with 1024 projections of 512 data samples each, and rays were rebinned to form 512 parallel-beam projections. Reconstructions were performed on a 512 × 512 grid using 100 iterations of the SIRT iterative CT algorithm. Proof of concept was demonstrated with a uniformly attenuating solution phantom. PRESAGE® dosimeters (11 cm diameter) were irradiated with Cobalt-60 irradiator to achieve either a uniform dose or a 2-level “step-dose” pattern. Results: With 6% refractive index mismatching, a circular field of view of 85% of the diameter of a cylindrical sample can be reconstructed accurately. Reconstructed images of the test solution phantom were uniform (within 3%) inside this radius. However, the dose responses of the PRESAGE® samples were not spatially uniform, with variations of at least 5% in sensitivity. The variation appears as a “cupping” artifact with less sensitivity in the middle than at the periphery of the PRESAGE® cylinder. Polarization effects were also detected for these samples. Conclusions: The fiducial-based ray path measurement scheme, coupled with an iterative reconstruction algorithm, enabled optical CT scanning of PRESAGE® dosimeters immersed in mismatched refractive index solutions. However, improvements to PRESAGE® dose response uniformity are required.

60 APPLIED LIFE SCIENCES↗

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics

In this article, strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Concentration Dependence of VO2+ Crossover of Nafion for Vanadium Redox Flow Batteries

The VO2+ crossover, or permeability, through Nafion in a vanadium redox flow battery (VRFB) was monitored as a function of sulfuric acid concentration and VO2+ concentration. A vanadium rich solution was flowed on one side of the membrane through a flow field while symmetrically on the other side a blank or vanadium deficit solution was flowed. The blank solution was flowed through an electron paramagnetic resonance (EPR) cavity and the VO2+ concentration was determined from the intensity of the EPR signal. Concentration values were fit using a solution of Fick s law that allows for the effect of concentration change on the vanadium rich side. The fits resulted in permeability values of VO2+ ions across the membrane. Viscosity measurements of many VO2+ and H2SO4 solutions were made at 30 60 C. These viscosity values were then used to determine the effect of the viscosity of the flowing solution on the permeability of the ion. 2013 The Electrochemical Society. [DOI: 10.1149/2.004306jes] All rights reserved.

Lawton, Jamie↗

Index tests for geosynthetic clay liners containing bentonite–polymer composites

The hydraulic conductivity of geosynthetic clay liners (GCLs) containing bentonite–polymer composites (BPCs) is controlled by a pore clogging mechanism related to the swelling of sodium bentonite (NaB) granules and formation of polymer gel in the pore space. Three index testing procedures that are influenced by bentonite swelling and rheology of the polymer gel were evaluated in this study for screening BPC GCLs for hydraulic conductivity: viscosity, liquid limit by fall-cone penetration and active water content. Experiments were conducted using polyacrylamide (PAM) polymer, dry-blended mixtures of PAM and NaB and a BPC from a commercially available GCL. These materials were contacted with inorganic solutions representing leachates from a range of waste containment scenarios. For the BPC that was evaluated, hydraulic conductivities less than 10 -10 m/s were obtained when the viscosity exceeded 220 mPa s and the active water content exceeded 300%. Hydraulic conductivities less than 10 -11 m/s were achieved for viscosity exceeding 350 mPa s, fall-cone liquid limit exceeding 190 and active water content exceeding 300%. Additional testing is needed to assess the generality of these index tests and the identified thresholds when applied to other BPC GCLs and other permeant liquids.

Engineering↗

Polymer Retention Evaluation in a Heavy Oil Sand for a Polymer Flooding Application on Alaska’s North Slope

For a polymer flooding field trial in a heavy oil reservoir on Alaska’s North Slope, polymer retention is a key parameter. Because of the economic impact of retention, this parameter was extensively studied using field core material and conditions. Several important findings were noted. Polymer retention based on effluent viscosity measurements can be overestimated unless the correct (non-linear) relation between polymer concentration and viscosity is used. Polymer degradation (either mechanical or oxidative) can also lead viscosity-based measurements to overestimate retention. Inaccessible pore volume (IAPV) can be overestimated if insufficient brine is flushed through the sand between polymer banks. Around 100 pore volumes (PV) of brine may be needed to displace mobile polymer to approach a true residual resistance factor and properly measure IAPV. Even for a sand pack with kwsor = 20 mD, IAPV was zero for HPAM with molecular weight (Mw) of 18 million g/mol. Fine-grain particles (<20 μm) strongly impacted polymer retention values. Native NB#1 sand with a significant component of particles <20 μm exhibited 290 μg/g, while the same sand exhibited 28 μg/g after these small particles were removed. Polymer retention did not necessarily correlate with mineral composition. The NB#1, NB#3, and OA sands had similar elemental and clay compositions, but the NB#1 sand exhibited ~10 times higher retention than the NB#3 sand. Polymer retention did not necessarily correlate with permeability. NB#1 sand exhibited much higher retention than OA sand, even though NB#1 sand was twice as permeable as OA sand. No evidence of chromatographic separation of HPAM molecular weights was found in our experiments. Although retention tended to be greater without a residual oil saturation (than at Sor), the effect was not strong. Aging a core (with high oil saturation) at 60°C reduced HPAM retention by a factor of two. As expected, under similar conditions, polymer retention was greater for a higher Mw HPAM (18 million g/mol) than for a lower Mw HPAM (10-12 million g/mol). In many cases with high polymer retention values (e.g., 240 μg/g), polymer arrival at the end of the core was relatively quick, but achieving the injected concentration occurred gradually over many pore volumes. This effect was not due chromatographic separation of polymer molecular weights. Results from modeling of this behavior was consistent with concentrationdependent polymer retention. The form assumed for the retention function in a simulator can have an important impact on the timing and magnitude of the oil response from a polymer flood. Field-based observations can underestimate polymer retention, depending on when the tracer and polymer concentrations were measured and the assumptions made about reservoir heterogeneity.

Wang, Dongmei↗

Ionic Liquids: Breakthrough Absorption Technology for Post-Combustion CO 2 Capture

This is the final report for DE-FC26-07NT43091 Ionic Liquids: Breakthrough Absorption Technology for Post-Combustion CO 2 Capture. A detailed summary is provided of the ionic liquid (IL) discovery process, synthesis and testing results, process / systems modeling, lab-scale operational testing, corrosion testing and commercialization possibilities. The work resulted in the discovery of a new class of ionic liquids (ILs) that efficiently react with CO 2 in a 1:1 stoichiometry with no water present and no increase in viscosity. The enthalpy of reaction was tuned to optimize process economics. The IL was found to have excellent corrosion behavior with and without CO 2 present. In lab-scale tests, the IL was able to effectively remove CO 2 from a simulated flue gas stream, although mass transfer was slower than with aqueous monoethanolamine (MEA) due to higher viscosities. The non-volatile nature of the solvent and its high thermal stability, however, make it an intriguing option. An independent systems analysis indicates that the economics of using the best IL discovered to date (NDIL0157), are at least comparable to and potentially slightly better than - the Fluor Econamine FG PlusTM process (DOE Case 12). Further work should be directed at improving mass transfer / lowering viscosity and developing commercial synthesis routes to make these ILs at scale in an inexpensive manner. Demonstration of the process at larger scales is also warranted, as is the exploration of other process configurations that leverage the anhydrous nature of the solvent and its extremely low volatility.

54 ENVIRONMENTAL SCIENCES↗

Carbon Dioxide-Water Emulsions for Enhanced Oil Recovery and Permanent Sequestration of Carbon Dioxide

This project involves the use of an innovative new invention – Particle Stabilized Emulsions (PSEs) of Carbon Dioxide-in-Water and Water-in-Carbon Dioxide for Enhanced Oil Recovery (EOR) and Permanent Sequestration of Carbon Dioxide. The EOR emulsion would be injected into a semi-depleted oil reservoir such as Dover 33 in Otsego County, Michigan. It is expected that the emulsion would dislocate the stranded heavy crude oil from the rock granule surfaces, reduce its viscosity, and increase its mobility. The advancing emulsion front should provide viscosity control which drives the reduced-viscosity oil toward the production wells. The make-up of the emulsion would be subsequently changed so it interacts with the surrounding rock minerals in order to enhance mineralization, thereby providing permanent sequestration of the injected CO 2 .

02 PETROLEUM↗

CO2-Binding Organic Liquids Gas Capture with Polarity-Swing-Assisted Regeneration Full Technology Feasibility Study B1 - Solvent-based Systems

PNNL, Fluor Corporation and Queens University (Kingston, ON) successfully completed a three year comprehensive study of the CO2BOL water-lean solvent platform with Polarity Swing Assisted Regeneration (PSAR). This study encompassed solvent synthesis, characterization, environmental toxicology, physical, thermodynamic and kinetic property measurements, Aspen Plus™ modeling and bench-scale testing of a candidate CO2BOL solvent molecule. Key Program Findings The key program findings are summarized as follows: • PSAR favorably reduced stripper duties and reboiler temperatures with little/no impact to absorption column • >90% CO2 capture was achievable at reasonable liquid-gas ratios in the absorber • High rich solvent viscosities (up to 600 cP) were successfully demonstrated in the bench-scale system. However, the projected impacts of high viscosity to capital cost and operational limits compromised the other levelized cost of electricity benefits. • Low thermal conductivity of organics significantly increased the required cross exchanger surface area, and potentially other heat exchange surfaces. • CO2BOL had low evaporative losses during bench-scale testing • There was no evidence of foaming during bench scale testing • Current CO2BOL formulation costs project to be $35/kg • Ecotoxicity (Water Daphnia) was comparable between CO2BOL and MEA (169.47 versus 103.63 mg/L) • Full dehydration of the flue gas was determined to not be economically feasible. However, modest refrigeration (13 MW for the 550 MW reference system) was determined to be potentially economically feasible, and still produce a water-lean condition for the CO2BOLs (5 wt% steady-state water loading). • CO2BOLs testing with 5 wt% water loading did not compromise anhydrous performance behavior, and showed actual enhancement of CO2 capture performance. • Mass transfer of CO2BOLs was not greatly impeded by viscosity • Facile separation of antisolvent from lean CO2BOL was demonstrated on the bench cart • No measurable solvent degradation was observed over 4 months of testing – even with 5 wt% water present

Heldebrant, David J↗

Low-Severity Hydroprocessing to Stabilize Bio-oil: TechnoEconomic Assessment

The impetus for this study was the suggestion that recent developments in fast pyrolysis (FP) bio-oil production had indicated instability of the bio-oil in storage which might lead to unacceptable viscosity increases. Commercial operation of FP in Finland began in 2014 and the distribution of the bio-oil to isolated users has been proposed as the long-term plan. Stability of the shipped bio-oil therefore became a concern. Experimental results at PNNL with low-severity hydroprocessing of bio-oil for stabilization has validated a process in which the stability of the bio-oil could be improved, as measured by viscosity increase following storage of the product at 80 °C for 24h. In the work reported here the assessed process configuration consists of fast pyrolysis followed by low temperature and pressure hydroprocessing to produce a stable fuel oil product. The product could then be stored for an extended period of time without significant viscosity increase. This work was carried out as part of a collaborative project between Technical Research Centre of Finland (VTT) and Pacific Northwest National Laboratory (PNNL). The public funding agents for the work were Tekes in Finland and the Bioenergy Technologies Office of the U.S. Department of Energy. The effort was proposed as an evaluation of the process developed in earlier collaboration and jointly invented by VTT and PNNL researchers.

biomass↗

Rheological and Mechanical Property Measurements of PMDI Foam at Elevated Temperatures

A study was undertaken to determine the viscosity of liquefied 20 lb/ft 3 poly methylene diisocyanate (PMDI) foam and the stress required to puncture solid PMDI foam at elevated temperatures. For the rheological measurements the foam was a priori liquefied in a pressure vessel such that the volatiles were not lost in the liquefaction process. The viscosity of the liquefied PMDI foam was found to be Newtonian with a power law dependence on temperature log 10 (μ/Pa s) = 20.6 – 9.5 log 10 (T/°C) for temperatures below 170 °C. Above 170 °C, the viscosity was in the range of 0.3 Pa s which is close to the lower measurement limit (≈ 0.1 Pa s) of the pressurized rheometer. The mechanical pressure required to break through 20lb/ft 3 foam was 500-800 psi at temperatures from room temperature up to 180 °C. The mechanical pressure required to break through 10 lb/ft 3 was 170-300 psi at temperatures from room temperature up to 180 °C. We have not been able to cause gas to break through the 20 lb/ft 3 PMDI foam at gas pressures up to 100 psi.

36 MATERIALS SCIENCE↗

Plasma Science and Innovation Center at Washington, Wisconsin, and Utah State: Final Scientific Report for the University of Wisconsin-Madison

The University of Wisconsin-Madison component of the Plasma Science and Innovation Center (PSI Center) contributed to modeling capabilities and algorithmic efficiency of the Non-Ideal Magnetohydrodynamics with Rotation (NIMROD) Code, which is widely used to model macroscopic dynamics of magnetically confined plasma. It also contributed to the understanding of direct-current (DC) injection of electrical current for initiating and sustaining plasma in three spherical torus experiments: the Helicity Injected Torus-II (HIT-II), the Pegasus Toroidal Experiment, and the National Spherical Torus Experiment (NSTX). The effort was funded through the PSI Center's cooperative agreement with the University of Washington and Utah State University over the period of March 1, 2005 - August 31, 2016. In addition to the computational and physics accomplishments, the Wisconsin effort contributed to the professional education of four graduate students and two postdoctoral research associates. The modeling for HIT-II and Pegasus was directly supported by the cooperative agreement, and contributions to the NSTX modeling were in support of work by Dr. Bickford Hooper, who was funded through a separate grant. Our primary contribution to model development is the implementation of detailed closure relations for collisional plasma. Postdoctoral associate Adam Bayliss implemented the temperature-dependent effects of Braginskii's parallel collisional ion viscosity. As a graduate student, John O'Bryan added runtime options for Braginskii's models and Ji's K2 models of thermal conduction with magnetization effects and thermal equilibration. As a postdoctoral associate, O'Bryan added the magnetization effects for ion viscosity. Another area of model development completed through the PSI-Center is the implementation of Chodura's phenomenological resistivity model. Finally, we investigated and tested linear electron parallel viscosity, leveraged by support from the Center for Extended Magnetohydrodynamic Modeling (CEMM). Work on algorithmic efficiency improved NIMROD's element-based computations. We reordered arrays and eliminated a level of looping for computations over the data points that are used for numerical integration over elements. Moreover, the reordering allows fewer and larger communication calls when using distributed-memory parallel computation, thereby avoiding a data starvation problem that limited parallel scaling over NIMROD's Fourier components for the periodic coordinate. Together with improved parallel preconditioning, work that was supported by CEMM, these developments allowed NIMROD's first scaling to over 10,000 processor cores. Another algorithm improvement supported by the PSI Center is nonlinear numerical diffusivities for implicit advection. We also developed the Stitch code to enhance the flexibility of NIMROD's preprocessing. Our simulations of HIT-II considered conditions with and without fluctuation-induced amplification of poloidal flux, but our validation efforts focused on conditions without amplification. A significant finding is that NIMROD reproduces the dependence of net plasma current as the imposed poloidal flux is varied. The modeling of Pegasus startup from localized DC injectors predicted that development of a tokamak-like configuration occurs through a sequence of current-filament merger events. Comparison of experimentally measured and numerically computed cross-power spectra enhance confidence in NIMROD's simulation of magnetic fluctuations; however, energy confinement remains an open area for further research. Our contributions to the NSTX study include adaptation of the helicity-injection boundary conditions from the HIT-II simulations and support for linear analysis and computation of 3D current-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗