SEARCH · Search NASA

Results for “constant flux”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Solutions of the equation of heat flow

The geometry of sunspots has been used to suggest a problem in heat flow. The equation of heat transport is solved for the case of a cylinder with a given thermal conductivity imbedded in an otherwise uniform medium with different conductivity. The surface of this region radiates heat with flux proportional to temperature. At a lower surface, either in heat flux or temperature is held constant. The cylinder can have an anisotropic thermal conductivity. The variations in temperature along the radiating surface have been determined. A simple approximation is noted which has been found to give a general solution with acceptable accuracy. This method may be of some use in other situations requiring the solution of Laplace's equation with a free surface. The analysis is used to set limits on the ratio of diameter to depth for cases which preserve the sharp surface temperature transition across the cylinder.

Margolis, S. H.↗

NASA NTRS1981

Bound oscillations on thin magnetic flux tubes - Convective instability and umbral oscillations

The possibility that 'tube waves' can be trapped on slender solar magnetic flux tubes is investigated. For rigid isothermal flux tubes, it is found that the flux tube geometry can by itself lead to waves which are trapped on the part of the tube that expands with height. Some geometries lead to trapped modes with eigenperiods near 180 s, if parameters appropriate to sunspot umbrae are chosen. It is possible that the umbral oscillations are a manifestation of such trapped waves, if sunspot umbrae consist of an assembly of slender flux tubes, as in the spaghetti model of Parker (1979). For flux tubes which have a constant ratio of Alfven speed to sound speed, it is found that it is primarily the variation of temperature with height which determines whether trapped waves can exist. Certain temperature profiles lead to disturbances for which omega squared is less than zero, corresponding to convective instability or Rayleigh-Taylor instability.

Hollweg, J. V.↗

NASA NTRS1981

Reference model atmosphere calculation - The Sunspot sunspot model

A composite temperature-density model of sunspots is developed from the results of five previous modeling attempts. The model is designed to describe the region extending from the deepest photospheric layers up through the chromosphere-corona transition region. Photospheric, chromospheric, and transition region temperature distributions proposed by other authors are employed to quantity the phenomena out to 2280 km. Additional calculations yield the various line profiles, the internal properties of the region above the sunspot including the chromospheric net radiative cooling rates, and the radiative flux as a function of depth in the photosphere. Observations of the hydrogen ion flux is used to account for the net cooling in the temperature minimum region. Furthermore, the integrated radiative flux is found to be constant with depth in the photosphere.

Avrett, E. H.↗

NASA NTRS1984

Very long baseline interferometry observations of the RS Canum Venaticorum system HR 5110 at 8.4 GHz

The RS CVn binary system HR 5110 was observed with a well-calibrated four-element VLBI array at a frequency of 8.4 GHz. The total flux density, which was nearly constant at 32 + or - 2 mJy, was monitored throughout the experiment with a 20 km baseline interferometer allowing accurate visibilities to be calculated. The observed visibilities are consistent with a source unresolved on all the baselines. The maximum size of the equivalent Gaussian source is 1.4 milli-arcsec (FWHM), corresponding to a linear size of 1.1 x 10 to the 12th cm at 52 pc. This maximum size is comparable with the overall size of the binary system. The inferred lower limit of the brightness temperature is 4 x 10 to the 8th K, which is higher than the coronal temperature of 10 to the 7th K measured in a soft X-ray survey. This is consistent with a nonthermal mechanism for the radiation process of the moderate radio outburst observed.

Mutel, R. L.↗

NASA NTRS1986

A model of the spatial and temporal variation of the Uranus thermal structure

Seasonal variability of the temperature structure of Uranus is modeled for all latitudes in the .0004 to 2 bar pressure range in anticipation of the Voyager encounter in January 1986. Atmospheric heating in the model results on the one hand from an internal heat source and, on the other hand, from absorption of solar energy by methane and by non-conservative aerosols located between the 0.5 and 2 bar levels. Various cases for the behavior of the internal heat flux are investigated, such as constant with latitude or constrained to yield a time-averaged thermal emission independent of latitude. Meridional transport of heat in the stably stratified atmosphere is not taken into account. The results indicate that the Voyager encounter time, very small north-south temperature asymmetry should be expected. Moreover, the northern hemisphere, although not illuminated, should emit as much energy (within one percent) as the southern hemisphere at this date. At a given latitude, extreme temperatures are reached at the equinoxes. At the poles, seasonal amplitudes of about 10 K in the upper stratosphere and 6 K at the 0.6 bar level are predicted, and the variation with time of the emission to space is found to be at most 20 percent. The atmosphere of Uranus appears to be characterized by very long radiative response times (mainly due to its cold temperature) which inhibit the large seasonal variations that one could otherwise expect in view of the high obliquity of the planet and its long orbital period.

Bezard, B.↗

NASA NTRS1987

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗

NASA NTRS1987

An evaluation of candidate oxidation resistant materials

Ground based testing of materials considered for Kapton solar array blanket protection, graphite epoxy structural member protection, and high temperature radiators was performed in an RF plasma asher. Ashing rates for Kapton were correlated with rates measured on STS-8 to determine the exposure time equivalent to one year in low Earth orbit (LEO) at a constant density space station orbital flux. Protective coatings on Kapton from Tekmat, Andus Corporation, and LeRC were evaluated in the plasma asher and mass loss rates per unit area were measured for each sample. All samples evaluated provided some protection to the underlying surface but ion beam sputter deposited samples of SiO2 and SiO2 with 8% polytetrafluoroethylene (PTFE) showed no evidence of degradation after 47 hours of exposure. Mica paint was evaluated as a protective coating for graphite epoxy structural members. Mica appears to be resistant to attack by atomic oxygen but only offers some limited protection as a paint because the paint vehicles evaluated to date were not resistant to atomic oxygen. Four materials were selected for evaluation as candidate radiator materials: stainless steel, copper, niobium-1% zirconium, and titanium-6% aluminum-4% vanadium. These materials were surface textured by various means to improve their emittance. Emittances as high as 0.93 at 2.5 microns for stainless steel and 0.89 at 2.5 microns for Nb-1 Zr were obtained from surface texturing. There were no significant changes in emittance after asher exposure.

Rutledge, Sharon↗

NASA NTRS1990

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

NASA NTRS1989

The Li-7 and Be-7 deexcitation lines: Probes for accelerated particle transport models in solar flares

The photon energy spectrum of a spectral feature composed of the 429 and 478 keV gamma-ray lines from Li-7 and Be-7 (produced by interactions of flare-accelerated alpha particles with ambient He in the solar atmosphere) depends on the angular distribution of the interacting accelerated particles. This spectrum is calculated for limb and disc-centered flares using a loop model for the transport of the ions. In this model, the flux tube magnetic field is constant in the corona and converges in the chromosphere to the photosphere. Magnetic mirroring and MHD pitch-angle scattering are both taken into account. Comparison of these results with data from other experiments is presented.

Murphy, Ronald J.↗

NASA NTRS1992

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

NASA NTRS1993

On the stability of irradiated accretion disks

The effects of external heating on the stability of hot accretion disks are studied in some detail. It is known that geometrically thin, optically thick, nonirradiated accretion disks have two distinct branches in the surface density mass-accretion rate plane: the upper branch is radiation pressure-dominated and is unstable against thermal and secular perturbations, while the lower one is gas-pressure-dominated and is stable. We show quite generally that, even when disks are strongly irradiated, the upper branch remains unstable and the lower branch remains stable; the lower branch, however, can become radiation pressure dominated, if the irradiating flux, F(irr), is kept constant. A stable, radiation pressure dominated state thus appears. If F(irr) changes in proportion to the mass-accretion rate through the disk, the instabilities associated with radiation pressure dominated disks cannot be removed. Some observational implications are discussed in the context of long-term variations of low-mass X-ray binaries.

Mineshige, Shin↗

NASA NTRS1994

Spectral energy distributions of the brightest Palomar-Green quasars at intermediate redshifts

We have combined low-dispersion International Ultraviolet Explorer (IUE) spectra with the optical/near-IR spectrophotometry of Neugebauer et al. (1987) in order to study the spectral energy distributions of seven of the brightest Palomar-Green (PG) quasars at intermediate redshifts (Z(sub em) greater than or equal to 0.9 and less than or equal to 1.5). Some of these PG quasars are barely detectable in long IUE exposures, so we have used the Gaussian Extraction (GEX) technique to maximize the signal-to-noise of the IUE data, and we have co-added all spectra available from the IUE archive for each QSO unless the ultraviolet spectra varied significantly from one exposure to the next. We have corrected the spectral energy distributions for Milky Way reddening using the observed neutral hydrogen column densities on each sight line and the gas-to-dust relation recently derived by Diplas & Savage. Six of the seven quasars are detected down to lambda much less than 700 A in the rest frame, and consequently continuum reddening due to dust in the immediate vicinity of the quasar can have a dramatic effect on the spectral energy distributions. In order to explore the possible importance of intrinsic continuum reddening, we have assembled a heuristic extinction curve which extends to lambda much less than 912 A. Using this heuristic extinction curve, we derive reasonable upper limits on the intrinsic E(B-V) for each quasar. We briefly discuss some of the implications of the derived intrinsic continuum reddening limits. We use geometrically thin accretion disk models to derive the black hole masses and accretion rates implied by the spectral energy distributions. Even if we neglect intrinsic reddening, we find that a large fraction of the quasars require super-Eddington accretion rates (which is not consistent with the thin disk assumption). Comparison of the data in this paper to a large body of data from the literature on the accretion disk M(sub BH) - M dot grid calculated by Wandel & Petrosian reveals that our quasars are among the brightest in the sky at 1450 A, and ostensibly suggests that the fraction of quasars which require super-Eddington accretion rates is much smaller than the fraction that we derive from our data alone. However, intrinsic continuum reddening has been ignored in this comparison, and a small amount of intrinsic reddening will push more of the quasars into the super-Eddington regime. We also plot the recent reverberation monitoring results on NGC 5548 and NGC 3783 on the Wandel & Petrosian grid, and we find that these Seyfert galaxies appear to vary along lines of constant M(sub BH). Continuum flux from two of the quasars in our main sample, PG 1338+416 and PG 1630+377, is detected at lambda(sub rest) less than 584 A. These quasars can in principle be used for the He I Gunn-Peterson test, but the S/N of IUE spectra of individual objects is usually too low to place interesting limits on the Gunn-Peterson optical depth. In order to improve the S/N, we have formed a composite spectrum from the spectra of five quasars detected with IUE at lambda(sub rest) less than 584 A, and we have used this composite spectrum to place a tighter limit on tau(sub GP, He I). We briefly discuss intermediate-redshift Lyman limit systems (Z(sub LL) greater than or equal to 0.5 and less than or equal to 1.6) detected in the IUE spectra of five quasars, including lower limits on N(H I) in each Lyman limit system.

Tripp, Todd M.↗

NASA NTRS1996

X-Ray Spectral Variability in NGC 7469

We present analyses of two Ginga observations and two observations from the ROSAT database of NGC 7469, focusing on the spectral variability observed on timescales of days and longer. During the 1988 Ginga observation, the hardness ratio (8-21 keV/3.4-5.7 keV) increased significantly as the total flux decreased by 30%. As the spectrum is well fit by the reflection model and since the spectra variability dominates the higher energy band, this could be explained by either a variation in the power law index or in the effective covering fraction of the reflecting material. This ambiguity is inherent in reflection modeling of Ginga spectra from moderate flux Seyfert 1 galaxies. Assuming that the power law index did not change, we find that the reflected flux is consistent with being constant, suggesting that much of the reflecting material may be located more than 3 light-days from the continuum source with the molecular torus being a plausible site. This scenario is also supported by the report of a narrow rather than broad iron K-alpha line in the ASCA data by Guainazzi et al. NGC 7469 was faint during the 1989 Ginga observation, but variability was observed with doubling timescale of 5 hr, and the spectrum was harder. A reflection component could not be constrained, and the change in the spectrum could be explained by an increase in neutral absorption. The brighter of two ROSAT spectra was significantly softer, and in both spectra there was evidence of spectral complexity, as has been previously reported by Turner, George, & Mushotzky and Brandt et al. The spectrum could be fit by a variety of two-component models, including a warm absorber model, an ionized disk model, and a thermal model with single-component blackbody spectrum, but joint fitting of the 1988 average Ginga spectrum and the nonsimultaneous ROSAT spectra favored thermal models, and other models required an anomalously high reflection ratio. This model is supported by the observation of a soft excess component and the lack of ionized absorption edges in the ASCA spectrum by Guainazzi et al. The long-term spectral variability could be explained by relative variability between the power-law and soft excess component normalizations, perhaps implying that hard X-ray reprocessing in thermal material does not dominate on long timescales.

Leighly, Karen↗

NASA NTRS2000

Protein Crystal Growth Dynamics and Impurity Incorporation

The general concepts and theories of crystal growth are proven to work for biomolecular crystallization. This allowed us to extract basic parameters controlling growth kinetics - free surface energy, alpha, and kinetic coefficient, beta, for steps. Surface energy per molecular site in thermal units, alpha(omega)(sup 2/3)/kT approx. = 1, is close to the one for inorganic crystals in solution (omega is the specific molecular volume, T is the temperature). Entropic restrictions on incorporation of biomolecules into the lattice reduce the incorporation rate, beta, by a factor of 10(exp 2) - 10(exp 3) relative to inorganic crystals. A dehydration barrier of approx. 18kcal/mol may explain approx. 10(exp -6) times difference between frequencies of adding a molecule to the lattice and Brownian attempts to do so. The latter was obtained from AFM measurements of step and kink growth rates on orthorhombic lysozyme. Protein and many inorganic crystals typically do not belong to the Kossel type, thus requiring a theory to account for inequivalent molecular positions within its unit cell. Orthorhombic lysozyme will serve as an example of how to develop such a theory. Factors deteriorating crystal quality - stress and strain, mosaicity, molecular disorder - will be reviewed with emphasis on impurities. Dimers in ferritin and lysozyme and acetylated lysozyme, are microheterogeneous i.e. nearly isomorphic impurities that are shown to be preferentially trapped by tetragonal lysozyme and ferritin crystals, respectively. The distribution coefficient, K defined as a ratio of the (impurity/protein) ratios in crystal and in solution is a measure of trapping. For acetylated lysoyzme, K = 2.15 or, 3.42 for differently acetylated forms, is independent of both the impurity and the crystallizing protein concentration. The reason is that impurity flux to the surface is constant while the growth rate rises with supersaturation. About 3 times lower dimer concentration in space grown ferritin and lysozyme crystals might be examples explaining higher quality of the space grown protein crystal. Depletion of solution with respect to isomorphic impurities around a growing crystal may be K times deeper than with respect to the crystallizing protein.

Chernov, Alex A.↗

NASA NTRS2006

The Boiling eXperiment Facility (BXF) for the Microgravity Science Glovebox (MSG)

Boiling is an effective means of cooling by removing heat from surfaces through vaporization of a working fluid. It is also affected by both the magnitude and direction of gravity. By conducting pool boiling tests in microgravity, the effect of buoyancy n the overall boiling process and the relative magnitude of other phenomena can be assessed. The Boiling eXperiment Facility (BXF) is being built for the Microgravity Science Glovebox. This facility will conduct two pool boiling studies. The first study the Microheater Array Boiling Experiment (MABE) uses two 96 element microheater arrays, 2.7 mm and 7.0 mm in size, to measure localized hear fluxes while operating at a constant temperature. The other experiment, the Nucleate Pool Boiling eXperiment (NPBX) uses a 85 mm diameter heater wafer that has been "seeded" with five individually-controlled nucleation sites to study bubble nucleation, growth, coalescence and departure. The BXF uses normal-perfluorohexane as the test fluid and will operate between pressures of 60 to 244 Pa. and temperatures of 35 to 60 C. Both sets of experimental heaters are highly instrumented. Pressure and bulk fluid temperature measurements will be made with standard rate video. A high speed video system will be used to visualize the boiling process through the bottom of the MABE heater arrays. The BXF is currently scheduled to fly on Utilization Flight-13A.1 to the ISS with facility integration into the MSG and operation during Increment 15

McQuillen, John↗

NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to replace the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model produces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed.

Secondary Ejecta↗

DOE OSTI2024

Accretion onto a Supermassive Black Hole Binary before Merger

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20 M (G = c = 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate $\dot{M}$ onto the black holes first decreases and then plateaus, dropping by only a factor of ~3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, "sloshing" streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless level $\phi$ ~ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PATCHWORK MHD, decreasing the computational expense to ~3% of conventional single-grid methods' cost.

79 ASTRONOMY AND ASTROPHYSICS↗