Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “aspect ratio”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Radar Retrieval Evaluation and Investigation of Dendritic Growth Layer Polarimetric Signatures in a Winter Storm

Abstract This study evaluates ice particle size distribution and aspect ratio φ Multi-Radar Multi-Sensor (MRMS) dual-polarization radar retrievals through a direct comparison with two legs of observational aircraft data obtained during a winter storm case from the Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. In situ cloud probes, satellite, and MRMS observations illustrate that the often-observed K dp and Z DR enhancement regions in the dendritic growth layer can either indicate a local number concentration increase of dry ice particles or the presence of ice particles mixed with a significant number of supercooled liquid droplets. Relative to in situ measurements, MRMS retrievals on average underestimated mean volume diameters by 50% and overestimated number concentrations by over 100%. IWC retrievals using Z DR and K dp within the dendritic growth layer were minimally biased relative to in situ calculations where retrievals yielded −2% median relative error for the entire aircraft leg. Incorporating φ retrievals decreased both the magnitude and spread of polarimetric retrievals below the dendritic growth layer. While φ radar retrievals suggest that observed dendritic growth layer particles were nonspherical (0.1 ≤ φ ≤ 0.2), in situ projected aspect ratios, idealized numerical simulations, and habit classifications from cloud probe images suggest that the population mean φ was generally much higher. Coordinated aircraft radar reflectivity with in situ observations suggests that the MRMS systematically underestimated reflectivity and could not resolve local peaks in mean volume diameter sizes. These results highlight the need to consider particle assumptions and radar limitations when performing retrievals. significance statement Developing snow is often detectable using weather radars. Meteorologists combine these radar measurements with mathematical equations to study how snow forms in order to determine how much snow will fall. This study evaluates current methods for estimating the total number and mass, sizes, and shapes of snowflakes from radar using images of individual snowflakes taken during two aircraft legs. Radar estimates of snowflake properties were most consistent with aircraft data inside regions with prominent radar signatures. However, radar estimates of snowflake shapes were not consistent with observed shapes estimated from the snowflake images. Although additional research is needed, these results bolster understanding of snow-growth physics and uncertainties between radar measurements and snow production that can improve future snowfall forecasting.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Effects of film thickness on electrochemical properties of nanoscale polyethylenedioxythiophene (PEDOT) thin films grown by oxidative molecular layer deposition (oMLD)

Poly(3,4-ethylene dioxythiophene) (PEDOT) has a high theoretical charge storage capacity, making it of interest for electrochemical applications including energy storage and water desalination. Nanoscale thin films of PEDOT are particularly attractive for these applications to enable faster charging. Recent work has demonstrated that nanoscale thin films of PEDOT can be formed using sequential gas-phase exposures via oxidative molecular layer deposition, or oMLD, which provides advantages in conformality and uniformity on high aspect ratio substrates over other deposition techniques. But to date, the electrochemical properties of these oMLD PEDOT thin films have not been well-characterized. In this work, we examine the electrochemical properties of 5–100 nm thick PEDOT films formed using 20–175 oMLD deposition cycles. We find that film thickness of oMLD PEDOT films affects the orientation of ordered domains leading to a substantial change in charge storage capacity. Interestingly, we observe a minimum in charge storage capacity for an oMLD PEDOT film thickness of ~30 nm (60 oMLD cycles at 150 °C), coinciding with the highest degree of face-on oriented PEDOT domains as measured using grazing incidence wide angle X-ray scattering (GIWAXS). Thinner and thicker oMLD PEDOT films exhibit higher fractions of oblique (off-angle) orientations and corresponding increases in charge capacity of up to 120 mA h g -1 . Electrochemical measurements suggest that higher charge capacity in films with mixed domain orientation arise from the facile transport of ions from the liquid electrolyte into the PEDOT layer. Greater exposure of the electrolyte to PEDOT domain edges is posited to facilitate faster ion transport in these mixed domain films. In conclusion, these insights will inform future design of PEDOT coated high-aspect ratio structures for electrochemical energy storage and water treatment.

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

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗
DOE OSTI2024

Automated Manufacturing of Grid Stiffened Panels with Radically Reduced Tooling

Grid stiffened continuous fiber reinforced composite panels are an attractive option for creating lightweight structures due to the tailorability for various applications and the resulting high specific properties. However, the panel stiffeners and stiffener intersections result in high tooling complexity and correspondingly high cost of implementation. These factors have limited the impact of such structures in the composites industry. Previous research has demonstrated the ability to produce high quality, high aspect ratio beams, representative of individual grid stiffeners, using E-glass/PET comingled tow via direct digital manufacturing. Further, prior preliminary efforts have demonstrated the potential to use the same approach to manufacture grid intersections that have continuous fiber in both directions. To expand on the previous efforts in grid stiffeners produced by direct digital manufacture with radically reduced tooling requirements, this effort compares two methods of providing positioning and consolidation, nozzle vs. roller. Both processes are based on a commingled yarn feedstock. The extrusion through a nozzle has been shown to enable grid intersection control through local variations in applied consolidation and serves as the baseline process. However, this approach requires a continuous placement path to create the complete grid stiffened panel as no mechanism for cutting and restarting has been implemented. Alternatively, a newly developed placement head incorporating cut and refeed, mounted to a 6-axis robot, offers the potential of improved path placement efficiency. The two techniques are used to produce similar grid composite stiffeners to evaluate the effectiveness of producing the grid intersections. Rate of deposition of the two end effectors are compared and the quality of the associated grid stiffeners, and intersections, are determined through measurement of geometry, fiber volume fraction and void fraction.

Engineering↗
DOE OSTI2023

Metal hydride heat exchanger and method of use

A heat exchanger incorporates a metal hydride heat exchanger and mitigates the fluid mixing process, and thus greatly improves the heat transfer efficiency and heat recovery processes. The metal hydride heat exchanger has a container for the metal hydride that has a large aspect ratio. A plurality of high aspect container for the metal hydride may be coupled with a manifold.

Bahar, Bamdad↗
DOE OSTI2023

The Effects of Shallow Cumulus Cloud Shape on Interactions Among Clouds and Mixing With Near-Cloud Environments

Cloud lifecycles are influenced by interactions among clouds and mixing with cloud-free environments. Cloud shapes dictate the length of cloud boundaries, which affects these interactions. This study uses a large eddy simulation to examine the relationships between the shallow cumulus cloud shapes, processes at their boundaries, and lifecycles. To characterize the shape of clouds, two metrics are newly defined: one reflects the overall horizontal cloud shape, and the other defines the cloud edge irregularity. Simulated clouds experience increased irregularity, but little changes in aspect ratio. Irregularity-driven cloud perimeter growth strongly indicates cloud splitting, with a more pronounced effect than aspect ratio-driven perimeter growth. Smaller cloud-shell gradients of properties across cloud boundaries are associated with more irregular clouds, suggesting enhanced mixing of clouds with the cloud-free environment. These findings indicate that the shallow cumulus evolution is partly driven by the irregularity of the cloud edge through lateral mixing and cloud splitting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Research three-dimensional hierarchical graphite architectures for anodes for fast charging

With current lithium ion batteries optimized for performance under relatively low charge rate conditions, implementation of XFC has been hindered by drawbacks including Li plating, kinetic polarization, and heat dissipation. This project will utilize model-informed design of 3-D hierarchical electrodes to tune key XFCrelated variables like 1) bulk porosity/tortuosity 2) vertical pore diameter, spacing, and lattice 3) crystallographic orientation of graphite particles relative to exposed surfaces 4) interfacial chemistry of the graphite surfaces through "artificial sEr' formation using ALD 5) current collector surface roughness (aspect ratio, roughness factor, etc.). A key aspect of implementing novel electrodes is characterizing them in relevant settings. For this project, ultimately led out of University of Michigan by Neil Dasgupta, that includes both coin cell and 2+ Ah pouch cell testing, as well as comparison testing against baselines. Sandia National Labs will be conducting detailed cell characterization on iterative versions/improvements of the model-based hierarchical electrodes, as well as COTS cells for baseline comparisons. Key metrics include performance under fast charge conditions, as well as the absence or degree of lithium plating. Sandia will use their unique high precision cycling and rapid EIS capabilities to accurately characterize performance and any lithium plating during 6C charging and beyond, coupling electrochemical observations with cell teardown. Sandia will also design custom fixturing to cool cells during rapid charge, to decouple any kinetic effects brought about by cell heating and allow comparisons between different cells and charge rates. Using these techniques, Sandia will assess HOH electrodes from the University of Michigan, as well as aiding in iterative model and electrode design.

25 ENERGY STORAGE↗
DOE OSTI2022

Pathways towards break even for low convergence ratio direct-drive inertial confinement fusion

Following indirect-drive experiments which demonstrated promising performance for low convergence ratios (below 17), previous direct-drive simulations identified a fusion-relevant regime which is expected to be robust to hydrodynamic instability growth. This paper expands these results with simulated implosions at lower energies of 100 and 270 kJ, and ‘hydrodynamic equivalent’ capsules which demonstrate comparable convergence ratio, implosion velocity and in-flight aspect ratio without the need for cryogenic cooling, which would allow the assumptions of one-dimensional-like performance to be tested on current facilities. A range of techniques to improve performance within this regime are then investigated, including the use of two-colour and deep ultraviolet laser pulses. Finally, further simulations demonstrate that the deposition of electron energy into the hotspot of a low convergence ratio implosion through auxiliary heating also leads to significant increases in yield. Results include break even for 1.1 MJ of total energy input (including an estimated 370 kJ of short-pulse laser energy to produce electron beams for the auxiliary heating), but are found to be highly dependent upon the efficiency with which electron beams can be created and transported to the hotspot to drive the heating mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Optimization of Triply Periodic Minimal Surface Heat Exchanger to Achieve Compactness, High Efficiency, and Low-Pressure Drop

With advancements in additive manufacturing (AM) techniques, high-quality triply periodic minimal surface (TPMS) structures can now be produced. TPMS walled heat exchangers (HX) hold significant potential for industrial applications and are receiving increasing attention. This paper explores the impact of various TPMS design variables on flow and thermal performance to optimize TPMS heat exchangers for compactness, high efficiency, and low pressure drop. The design variables examined include the type of TPMS lattice, unit cell size, wall thickness, aspect ratio, TPMS orientation, and equivalent thickness. The study reveals that the flow and heat transfer performance of TPMS structures are significantly affected by these design variables. For the Gyroid, Diamond, and SplitP lattices, performance is nearly identical when the surface-to-volume ratio is kept constant. The average velocity of the fluid in the TPMS HX should be 0.3 m/s. The corresponding Re is between 300~800. Thin wall thickness, small equivalent thickness, and flat lattice configurations can significantly reduce pressure drop while maintaining the overall heat transfer coefficient. Additionally, the angle between the flow direction and TPMS orientation can increase pressure drop. Three aluminum heat exchangers were successfully printed using an AM machine, and testing results are comparable with theoretical prediction.

42 ENGINEERING↗
DOE OSTI2023

SiC Microchannel Heat Sinks for High Heat Flux Dissipation of 1 kW/cm 2

Thermal management has become a major barrier in the development of high-power micro-devices. Microchannel heat sinks (MCHS) are a promising method to dissipate high heat fluxes from high power density devices. Reducing the thermal resistance will reduce the temperature increase of the electronic device, allowing for higher powers to be achieved while maintaining a safe operating temperature. This paper reports on the use of silicon carbide (Sic) microchannels to dissipate high heat fluxes up to 1 kW/cm 2 over a footprint of 0.25 cm 2 . The use of high aspect ratio and decreased hydraulic diameter channels has allowed for very low thermal resistances to be achieved at higher heat duties compared to previous studies. The experiments performed resulted in heater temperatures of less than 60 °C and thermal resistances as low as 0.024 cm 2 . ° $C$/W between the heater surface and fluid inlet temperature. In conclusion, the demonstrated thermal performance of the SiC MCHS is an encouraging means for dissipating heat from high power density microelectronic devices.

42 ENGINEERING↗
DOE OSTI2021

An Asymmetric Dust Ring around a Very Low Mass Star ZZ Tau IRS

We present Atacama Large Millimeter/submillimeter Array (ALMA) gas and dust observations at band 7 (339 GHz: 0.89 mm) of the protoplanetary disk around a very low mass star ZZ Tau IRS with a spatial resolution of 025. The {sup 12}CO J = 3 → 2 position–velocity diagram suggests a dynamical mass of ZZ Tau IRS of ∼0.1–0.3 M {sub ☉}. The disk has a total flux density of 273.9 mJy, corresponding to an estimated mass of 24–50 M {sub ⊕} in dust. The dust emission map shows a ring at r = 58 au and an azimuthal asymmetry at r = 45 au with a position angle of 135°. The properties of the asymmetry, including radial width, aspect ratio, contrast, and contribution to the total flux, were found to be similar to the asymmetries around intermediate mass stars (∼2 M {sub ☉}) such as MWC 758 and IRS 48. This implies that the asymmetry in the ZZ Tau IRS disk shares a similar origin with others, despite the star being ∼10 times less massive. Our observations also suggest that the inner and outer parts of the disk may be misaligned. Overall, the ZZ Tau IRS disk shows evidence of giant planet formation on a ∼10 au scale at a few megayears. If confirmed, it will challenge existing core accretion models in which such planets have been predicted to be extremely hard to form around very low mass stars.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Hydro-mechanical coupled modeling of hydraulic fracturing using the mesh fragmentation technique

The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.

58 GEOSCIENCES↗
DOE OSTI2024

Quasilinear modelling of collisional trapped electron modes

The simplification of collision operators is necessary for quasilinear turbulence modelling used with integrated modelling frameworks, such as the gyrokinetic code QuaLiKiz. The treatment of collisions greatly impacts the accuracy of trapped electron mode (TEM) modelling, which is necessary to predict the electron heat flux and the balance between inward and outward particle fluxes. In particular, accurate particle flux predictions are necessary to successfully model density peaking in the tokamak core. We explored two ways of improving collisional TEM model reduction for tokamak core plasmas. First, we carried out linear GENE simulations to study the complex interplay between collisions and trapped electrons. We then used these simulations to define an effective trapped fraction to characterize the collisional TEM based on two key parameters, the local inverse-aspect ratio $\epsilon$ and the collisionality $\nu ^\ast$. One aspect missing from analytical TEM research is that the collisional frequency and the bounce-transit frequency are both velocity dependent; this effective trapped fraction takes both into account. In doing so, we determined that two parameters are not enough to model the collisional TEM, as an additional third free parameter was necessary. We determined that this model, as currently formulated, is not suitable for integrated modelling purposes. Second, we directly improved QuaLiKiz's Krook operator, which relies on two free parameters. We determined that these parameters required adjustments against higher-fidelity collisional models. In order to improve density profile predictions when paired with integrated models, we refined the Krook operator by using GENE simulations as a higher-fidelity point of comparison. We then demonstrate strong improvement of density peaking predictions of QuaLiKiz within the integrated modelling framework JETTO.

fusion plasma↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An asymptotic Grad–Shafranov equation for quasisymmetric stellarators

A first-order model is derived for quasisymmetric stellarators where the vacuum field due to coils is dominant, but plasma-current-induced terms are not negligible and can contribute to magnetic differential equations, with $\beta$ of the order of the ratio induced to vacuum fields. Under these assumptions, it is proven that the aspect ratio must be large and a simple expression can be obtained for the lowest-order vacuum field. The first-order correction, which involves both vacuum and current-driven fields, is governed by a Grad–Shafranov equation and the requirement that flux surfaces exist. These two equations are not always consistent, and so this model is generally overconstrained, but special solutions exist that satisfy both equations simultaneously. One family of such solutions is the set of first-order near-axis solutions. Thus, the first-order near-axis model is a subset of the model presented here. Several other solutions outside the scope of the near-axis model are also found. A case study comparing one such solution to a VMEC-generated solution shows good agreement.

Nikulsin, Nikita (ORCID:0000000318611777)↗
DOE OSTI2022

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2 /M 1 ) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r $\in$ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ~0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ~10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ~ 0.01 and sustain deep, asymmetric cavities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Evolution of Accreting Binaries: From Brown Dwarfs to Supermassive Black Holes

Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios (q ≡ M 2 /M 1 ) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios H/r ϵ {0.1, 0.05, 0.03} and vary both the magnitude and spatial dependence of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with 0.01 ≤ q ≤ 0.04 outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outward. However, in a number of cases, the disks accreting onto binaries with mass ratios ~0.07 fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the ~10% level even for mass ratios of 0.01. We demonstrate that the spatial dependence of the viscosity (e.g., α vs. constant ν) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at q ~ 0.01 and sustain deep, asymmetric cavities.

79 ASTRONOMY AND ASTROPHYSICS↗