Observational Signatures of Coronal Heating in Magnetohydrodynamic Simulations without Radiation or a Lower Atmosphere
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Poynting vector for electromagnetic energy arising from MHD disturbances, specifically for solar wind
We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .
To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.
We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.
We present a comprehensive analysis of super-Eddington black hole accretion simulations that solve the GRMHD equations coupled with angle-discretized radiation transport. The simulations span a range of accretion rates, two black hole spins, and two magnetic field topologies, and include resolution studies as well as comparisons with nonradiative models. Super-Eddington accretion flows consistently develop geometrically thick disks supported by radiation pressure, regardless of magnetic field configuration. Radiation generated in the inner disk drives substantial outflows, forming conical funnel regions that limit photon escape and result in very low radiation efficiency. The accretion flows are highly turbulent, with thermal energy transport dominated by radiation advection rather than diffusion. Angular momentum is primarily carried outward by Maxwell stress, with turbulent Reynolds stress playing a subdominant role. Both strong and weak jets are produced. Strong jets arise from sufficient net vertical magnetic flux and rapid black hole spin, and they can effectively evacuate the funnel, enabling radiation to escape through strong geometric beaming. In contrast, weak jets fail to clear the funnel, which becomes obscured by radiation-driven outflows and leads to distinct observational signatures. Spiral structures are observed in the plunging region, behaving like density waves. These super-Eddington models are applicable to a variety of astronomical systems, including ultraluminous X-ray sources, little red dots, and black hole transients.
Cosmic ray observations initiated by a solar flare are studied to measure the intensity and energy spectra of galactic cosmic ray protons in the 3-600 mev range by explorer xii
Dynamical behavior of cosmic ray gases, noting modes of hydromagnetic wave propagation and static equilibrium configurations for thermal interstellar gas and field
Stability of interstellar gas to perturbations in gas pressure, magnetic field and cosmic ray pressure from hydromagnetic viewpoint taking rotation into account
The properties of the synchrotron continuum are calculated on the basis of a steady spherically symmetric MHD model of the flow in the Crab nebula. The model accounts for the spectral distribution of the continuum from the optical to gamma-rays, and semiquantitatively for its observed spatial distribution. The model does not account for the radio portion of the Crab nebula's spectrum.
The Science Working Group for the High Resolution Solar Observatory (HRSO) laid plans beginning in 1984 for a series of workshops designed to stimulate a broadbased input from the scientific community to the HRSO mission. These workshops have the dual objectives of encouraging an early start on the difficult theoretical problems in radiative transfer, magnetohydrodynamics, and plasma physics that will be posed by the HRSO data, and maintaining current discussions of results in high resolution solar studies. This workshop was the second in the series. The workshop format presented invited review papers during the formal sessions and contributed poster papers for discussions during open periods. Both are presented.
In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.
Solar ionizing radiation critically depends on the level of the Sun’s magnetic activity. For robust physics-based forecasts, we employ the procedure of data assimilation, which combines theoretical modeling and observational data such that uncertainties in both the model and the observations are taken into account. Currently we are working in two major directions: 1) development of a new long-term forecast procedure on time-scales of the 11-year solar cycle, using a 2-dimensional mean-field dynamo model and synoptic magnetograms; 2) development of 3-dimensional radiative MHD (Magnetohydrodynamic) simulations to investigate the origin and precursors of local manifestations of magnetic activity, such as the formation of magnetic structures and eruptive dynamics.
In the present paper, we use Si IV 1393.755 Å spectral line observed by the Interface Region Imaging Spectrograph (IRIS) in the quiet-Sun to determine physical nature of the solar transition region (TR) oscillations. We analyze the properties of these oscillations using wavelet tools (e.g., power, cross-power, coherence, and phase difference) along with the stringent noise model (i.e., power-law + constant). We estimate the period of the intensity and Doppler velocity oscillations at each chosen location in the quiet-Sun (QS) and quantify the distribution of the statistically significant power and associated periods in one bright and two dark regions. In the bright TR region, the mean periods in intensity and velocity are 7 min, and 8 min respectively. In the dark region, the mean periods in intensity and velocity are 7 min, and 5.4 min respectively. We also estimate the phase difference between the intensity and Doppler velocity oscillations at each location. The statistical distribution of phase difference is estimated, which peaks at -119°± 13°, 33°± 10°, 102°± 10° in the bright region, while at -153°± 13°, 6°± 20°, 151°± 10° in the dark region. The statistical distribution reveals that the oscillations are caused by propagating slow magnetoacoustic waves encountered with the TR. Some of these locations may also be associated with the standing slow waves. Even, in the given time domain, several locations exhibit presence of both propagating and standing oscillations at different frequencies.
The interaction of the Crab pulsar with its surrounding nebula is modeled as a steady state, spherically symmetric relativistic wind which is confined by the outer layers of the supernova remnant. The pulsar wind is assumed to contain only positrons and electrons, and is terminated by a strong shock at a distance of 10 arcseconds from the pulsar, at the outer edge of the underluminous region. Down-stream of the shock, the flow decelerates and increases its pressure in order to match the boundary conditions imposed by outer nebula. Consistent flow solutions can be found if the magnetic luminosity of the pulsar wind is about 0.3 percent of its particle luminosity. The immediate postshock plasma is modeled as a power law distribution whose density and pressure moments satisfy the strong shock Rankine-Hugoniot relations. Using the downstream flow solution, the synchrotron continuum is calculated. The parameters of the pulsar wind and assumed power law spectral index are adjusted until the calculated total synchrotron power and the spectral distribution of the nebular continuum from optical to gamma-rays, are in reasonable accord with observations.
Previous studies of the magnetohydrodynamic (MHD) stability of solar coronal loops have not taken into account the effects of radiative energy loss in the energy equation. However, since coronal loops continuously lose energy by radiation and heat conduction, it is important to understand how these energy loss mechanisms affect MHD stability. We investigate the problem assuming that a magnetic loop has cylindrical geometry. As a first step, stability is studied for a localized mode, and the result is applied to a specific equilibrium. We find that the radiative energy loss effect not only changes the growth rate of ideally unstable modes, but also alters the stability boundary predicted by ideal MHD theory.
A solar oven heated by concentrated solar radiation as the heat source of a liquid metal magnetohydrodynamic (LMMHD) power generation system is proposed. The design allows the production of electric power in space, as well as on Earth, at high rates of efficiency. Two types of the solar oven suitable for the system are discussed.