Magnetohydrodynamic Simulations Research Articles

Asymmetrical Looping Magnetic Fields and Marsward Flows on the Nightside of Mars

AbstractAs the interplanetary magnetic field (IMF) carried by the solar wind encounters the martian atmosphere, it tends to pile up and drape around the planet, forming looping magnetic fields and inducing marsward ion flows on the nightside. Previous statistical observations revealed asymmetrical distribution features within this morphology; however, the underlying physical mechanism remains unclear. In this study, utilizing a three‐dimensional multi‐fluid magnetohydrodynamic simulation model, we successfully reproduce the asymmetrical distributions of the looping magnetic fields and corresponding marsward flows on the martian nightside. Analyzing the magnetic forces resulting from the bending of the IMF over the polar area, we find that the asymmetry is guided by the orientation of the solar wind motional electric field (ESW). A higher solar wind velocity leads to enhanced magnetic forces, resulting in more tightly wrapped magnetic fields with an increased efficiency in accelerating flows as they approach closer to Mars.

Magnetohydrodynamic simulations of the space weather in Proxima b: Habitability conditions and radio emission

Magnetohydrodynamic simulations of the space weather in Proxima b: Habitability conditions and radio emission

Fast ion relaxation in ITER mediated by Alfvén instabilities

We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

Assessment of runaway electron beam termination and impact in ITER

The vertical motion and shrinking of the cold plasma column after a tokamak disruption leads to a natural decrease in the edge safety factor when most of the current is carried by runaway electrons (REs). Reaching a low edge safety factor can potentially cause a strong plasma instability. We present magnetohydrodynamic simulations of the termination of a post-disruption plateau-phase RE beam in ITER when the edge safety factor falls close to two. Growth of instabilities is observed to result in stochastization of the magnetic field and a prompt loss of REs. As RE impact must be mitigated in ITER, the effect of parameters that influence the final termination have been assessed. Higher background plasma resistivity is seen to cause larger mode magnitudes and stronger stochastization, leading to less remnant REs after the termination event. Lower ion-densities also project a qualitatively similar behavior although weaker in effect. Using computations from a wall collision model, the ensuing load distribution on the first-wall is also presented.

Co-evolution of dust grains and protoplanetary disks. II. Structure and evolution of protoplanetary disks: An analytical approach

Abstract In our previous study (Tsukamoto et al. 2023b, PASJ, 75, 835), we investigated the formation and early evolution of protoplanetary disks with 3D non-ideal magnetohydrodynamics simulations considering dust growth, and found that the modified equations of the conventional steady accretion disk model that consider magnetic braking, dust growth, and ambipolar diffusion reproduce the disk structure (such as density and vertical magnetic field) obtained from simulations very well. In this paper, as a sequel to our previous study, we analytically investigate the structure and evolution of protoplanetary disks corresponding to Class 0/I young stellar objects using the modified steady accretion disk model combining an analytical model of envelope accretion. We estimate that the disk radius is several astronomical units at the disk formation epoch and increases to several hundred astronomical units at the end of the accretion phase. The disk mass is estimated to be $0.01 \lesssim M_{\rm disk} \lesssim 0.1 \, M_\odot$ for a disk with a radius of several tens of astronomical units and a mass accretion rate of $\dot{M}_{\rm disk} \sim 10^{-6} \, M_\odot \,\, {\rm yr^{-1}}$. These estimates seems to be consistent with recent observations. We also found that, with typical disk ionization rates (ζ ≳ 10−19 s−1) and a moderate mass accretion rate ($\dot{M}_{\rm disk}\gtrsim 10^{-8} \, M_\odot \,\, {\rm yr^{-1}}$), magnetorotational instability is suppressed in the disk because of low plasma β and efficient ambipolar diffusion. We argue that the radial profile of specific angular momentum (or rotational velocity) at the disk outer edge should be continuously connected to that of the envelope if the disk evolves by magnetic braking, and should be discontinuous if the disk evolves by an internal angular momentum transport process such as gravitational instability or magnetorotational instability. Future detailed observations of the specific angular momentum profile around the disk outer edge are important for understanding the angular momentum transport mechanism of protoplanetary disks.

AuriDESI: mock catalogues for the DESI Milky Way Survey

ABSTRACT The Dark Energy Spectroscopic Instrument Milky Way Survey (DESI MWS) will explore the assembly history of the Milky Way by characterizing remnants of ancient dwarf galaxy accretion events and improving constraints on the distribution of dark matter in the outer halo. We present mock catalogues that reproduce the selection criteria of MWS and the format of the final MWS data set. These catalogues can be used to test methods for quantifying the properties of stellar halo substructure and reconstructing the Milky Way’s accretion history with the MWS data, including the effects of halo-to-halo variance. The mock catalogues are based on a phase-space kernel expansion technique applied to star particles in the Auriga suite of six high-resolution lambda-cold dark matter magnetohydrodynamic zoom-in simulations. They include photometric properties (and associated errors) used in DESI target selection and the outputs of the MWS spectral analysis pipeline (radial velocity, metallicity, surface gravity, and temperature). They also include information from the underlying simulation, such as the total gravitational potential and information on the progenitors of accreted halo stars. We discuss how the subset of halo stars observable by MWS in these simulations corresponds to their true content and properties. These mock Milky Ways have rich accretion histories, resulting in a large number of substructures that span the whole stellar halo out to large distances and have substantial overlap in the space of orbital energy and angular momentum.

Generation and annihilation of three dimensional magnetic nulls in extrapolated solar coronal magnetic field: data-based Implicit Large Eddy simulation

Three-dimensional (3D) magnetic nulls are the points where magnetic field vanishes and are preferential sites for magnetic reconnection: a fundamental process which converts magnetic energy into kinetic energy, heat, and energy of non-thermal particles along with a rearrangement of magnetic field lines. Reconnection is ubiquitous in nature and plays a major role in various magnetically confined laboratory and space/astrophysical plasmas. In the solar corona, the reconnection manifests as coronal transients including solar flares, coronal mass ejections and coronal jets—often associated with 3D nulls. The nulls are generally found to be collocated with complex active regions on the solar photosphere and merits further attention, particularly in terms of their generation. A recent idealized magnetohydrodynamics simulation initiated with an analytically constructed preexisting proper radial null has identified magnetic reconnection to be responsible for spontaneous generation of these 3D nulls. It is then imperative to further explore the plausibility of spontaneous generation of nulls in naturally occurring plasmas, identify the mechanism and verify the outcome vis-à-vis observations. An apt test bed for such an initiative is the solar atmosphere, as abundant space and ground-based observations are available. In the above backdrop, the paper attempts to investigate 3D null generation by carrying out a data-based simulation of a C6.6 class flare associated with the photospheric active region NOAA 11 977. The simulation confirms spontaneous pairwise generation of 3D nulls with magnetic reconnections as the underlying cause. Importantly, magnetic field lines associated with the spontaneously generated nulls are found to trace observed chromospheric bright points—highlighting their observational relevance. Overall, such spontaneous generation and annihilation of nulls through magnetic reconnections opens up a new avenue for solar coronal and chromospheric heating.

Dynamics near the inner dead-zone edges in a proprotoplanetary disk

Abstract We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Elsässer numbers of the Ohmic resistivity and ambipolar diffusion is divided into two regions: “the transition zone” and “the coherent zone.” The coherent zone has the same properties as the conventional dead zone, and extends outside of the transition zone in the radial direction. Between the active and coherent zones, we discover the transition zone, the inner edge of which is identical to that of the conventional dead zone. The transition zone extends out over the regions where thermal ionization determines diffusion coefficients. The transition zone has completely different physical properties than the conventional dead zone, the so-called undead zone, and the zombie zone. The combination of amplification of the radial magnetic field owing to the ambipolar diffusion and a steep radial gradient of the Ohmic diffusivity causes the efficient evacuation of the net vertical magnetic flux from the transition zone within several rotations. Surface gas accretion occurs in the coherent zone but not in the transition zone. The presence of the transition zone prohibits mass and magnetic flux transport from the coherent zone to the active zone. Mass accumulation occurs at both edges of the transition zone as a result of mass supply from the active and coherent zones.

Description of turbulent dynamics in the interstellar medium: Multifractal microcanonical analysis. II. Sparse filtering of Herschel observation maps and visualization of filamentarystructures at different length scales

We present significant improvements to our previous work on noise reduction in Herschel observation maps by defining sparse filtering tools capable of handling, in a unified formalism, a significantly improved noise reduction as well as a deconvolution in order to reduce effects introduced by the limited instrumental response (beam). We implement greater flexibility by allowing a wider choice of parsimonious priors in the noise-reduction process. More precisely, we introduce a sparse filtering and deconvolution approach approach of type $l^2$-$l^p$, with $p > 0$ variable and apply it to a larger set of molecular clouds using Herschel $m data in order to demonstrate their wide range of application. In the Herschel data, we are able to use this approach to highlight extremely fine filamentary structures and obtain singularity spectra that tend to show a significantly less log -normal behavior and a filamentary nature in the less dense regions. We also use high-resolution adaptive magneto-hydrodynamic simulation data to assess the quality of deconvolution in such a simulated beaming framework.

Sensitivity of magnetohydrodynamic simulations of Joule-heated conductors to the vaporization curve in equations of state.

Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equationsof state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. The inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor's expansion and diverted the metal into the warm dense matter regime.

Three-dimensional GRMHD simulations of rapidly rotating stellar core collapse

ABSTRACT We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M⊙ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Among our models, only those with cores having an initial magnetic field strength of 1012 G and rotation rates of 1 or 2 rad s−1 produce MHD jets. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same order-of-magnitudes. Due to the memory effect, the frequency of the neutrino GW from our full-fledged 3D-MHD models is in the range less than ∼10 Hz. Toward the future GW detection for a Galactic core-collapse supernova, if driven by the MR mechanism, the planned next-generation detector as DECIGO is urgently needed to catch the low-frequency signals.

Frozen-field Modeling of Coronal Condensations with MPI-AMRVAC. I. Demonstration in Two-dimensional Models

Large-scale coronal plasma evolutions can be adequately described by magnetohydrodynamics (MHD) equations. However, full multidimensional MHD simulations require substantial computational resources. Given the low plasma β in the solar corona, in many coronal studies, it suffices to approximate the magnetic field to remain topologically fixed and effectively conduct one-dimensional (1D) hydrodynamic (HD) simulations instead. This approach is often employed in studies of coronal loops and their liability to form condensations related to thermal instability. While 1D HD simulations along given and fixed field line shapes are convenient and fast, they are difficult to directly compare with multidimensional phenomena. Therefore, it is more convenient to solve volume-filling, multidimensional versions of the MHD equations where we freeze the magnetic field, transforming it into frozen-field HD (ffHD) equations for simulation. We have incorporated this ffHD module into our open-source MPI-AMRVAC code and tested it using a two-dimensional (2D) evaporation–condensation model to study prominence formation due to radiative losses. The 2D ffHD results are compared with those from actual 2D MHD and pseudo-2D HD simulations, analyzing the differences and their causes. Pseudo-2D studies account for the known flux tube expansion effects. Overall, the performance of 2D ffHD is close to that of 2D MHD and pseudo-2D HD. The 2D tests conducted in this paper will be extended in follow-up studies to 3D simulations based on analytical or observational approaches.

Numerical Calculations of Adiabatic Invariants From MHD‐Driven Magnetic Fields

AbstractThe adiabatic invariants (M, J, Φ) and the related invariants (M, K, L*) have been established as effective coordinate systems for describing radiation belt dynamics at a theoretical level, and through numerical techniques, can be paired with in situ observations to order phase‐space density. To date, methods for numerical techniques to calculate adiabatic invariants have focused on empirical models such the Tsyganenko models TS05, T96, and T89. In this work, we develop methods based on numerical integration and variable step size iteration for the calculation of adiabatic invariants, applying the method to the Lyon‐Fedder‐Mobarry (LFM) global magnetohydrodynamics (MHD) simulation code, with optional coupling to the Rice Convection Model (RCM). By opening the door to adiabatic invariant modeling with MHD magnetic fields, the opportunity for exploratory modeling work of radiation belt dynamics is enabled. Calculations performed using LFM are cross‐referenced with the same code applied to the T96 and TS05 Tsyganenko models evaluated on the LFM grid. Important aspects of the adiabatic invariant calculation are reviewed and discussed, including (a) sensitivity to magnetic field model used, (b) differences in the problem between quiet and disturbed geomagnetic states, and (c) the selection of key parameters, such as the magnetic local time step size for drift shell determination. The rigorous development and documentation of this algorithm additionally acts as preliminary step for future thorough reassessment of in situ phase‐space density results using alternative magnetic field models.

Uniturbulence statistics and analysis of factors influencing the energy spectrum

This paper investigates the dynamics of unidirectionally propagating surface Alfvén waves, employing magnetohydrodynamic numerical simulations and statistical methodologies. The primary goal of this work is to enhance our understanding of the nonlinear self-cascade of surface Alfvén waves, which we term as uniturbulence, by unraveling the complex relationships between various length scales and their interplay with turbulent energy transfer mechanisms. To achieve this, we extensively analyze the phenomenon of uniturbulence using methods such as power spectrum analysis, radially averaged Fourier transform, and kurtosis. We employ these techniques to investigate the spatiotemporal distributions of kinetic and magnetic energy in uniturbulent flows. We also reveal the crucial role of the density contrast's variations and the role of Yaglom's law in characterizing energy transfer mechanisms. Our findings reveal that the inertial range of the perpendicular kinetic energy and magnetic energy along the z-axis depicts a progressive change in slope values, ultimately approaching the often-observed values of −5/3 and −3/2, respectively. Furthermore, our kurtosis analysis highlights the non-Gaussian behavior of the flow field at different length scales and over time, offering a perspective on uniturbulence dynamics. The correlations observed among diverse statistical approaches emphasize the complex interplay between different length scales in the context of uniturbulence. Our findings contribute to understanding this phenomenon, establishing a basis for future investigations to clarify the connections regulating these turbulent dynamics.

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount Te2/Te2,ad−1≃0.0016Ms3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1/T i1, which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1/T i1 = 1, whereas heating by B-parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1/T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Modeling the secular evolution of embedded protoplanetary disks

Context. Protoplanetary disks are known to form around nascent stars from their parent molecular cloud as a result of angular momentum conservation. As they progressively evolve and dissipate, they also form planets. While a lot of modeling efforts have been dedicated to their formation, the question of their secular evolution, from the so-called class 0 embedded phase to the class II phase where disks are believed to be isolated, remains poorly understood. Aims. We aim to explore the evolution between the embedded stages and the class II stage. We focus on the magnetic field evolution and the long-term interaction between the disk and the envelope. Methods. We used the GPU accelerated code IDEFIX to perform a 3D, barotropic, non ideal magnetohydrodynamic (MHD) secular core collapse simulation that covers the system evolution from the collapse of the pre-stellar core until 100 kyr after the first hydrostatic core formation and the disk settling while ensuring sufficient vertical and azimuthal resolutions (down to 10−2 au) to properly resolve the disk internal dynamics and non axisymmetric perturbations. Results. The disk evolution leads to a power-law gas surface density in Keplerian rotation that extends up to a few 10 au. The magnetic flux trapped in the disk during the initial collapse decreases from 100 mG at disk formation down to 1 mG by the end of the simulation. After the formation of the first hydrostatic core, the system evolves in three phases. A first phase with a small (∼10 au), unstable, strongly accreting (∼ 10−5 M⊙ yr−1) disk that loses magnetic flux over the first 15 kyr, a second phase where the magnetic flux is advected with a smooth, expanding disk fed by the angular momentum of the infalling material, and a final phase with a gravitationally regulated ∼60 au disk accreting at at few 10−7 M⊙ yr−1. The initial isotropic envelope eventually feeds large-scale vertically extended accretion streamers, with accretion rates similar to that onto the protostar (∼ 10−6 M⊙ yr−1). Some of the streamer material collides with the disk’s outer edge and produces accretion shocks, but a significant fraction of the material lands on the disk surface without producing any noticeable discontinuity. Conclusions. While the initial disk size and magnetization are set by magnetic braking, self-gravity eventually drives accretion, so that the disk ends up in a gravitationally regulated state. This evolution from magnetic braking to self-gravity is due to the weak coupling between the gas and the magnetic field once the disk has settled. The weak magnetic field at the end of the class I phase (Bz ∼ 1 mG) is a result of the magnetic flux dilution in the disk as it expands from its initial relatively small size. This expansion should not be interpreted as a viscous expansion, as it is driven by newly accreted material from large-scale streamers with large specific angular momentum.

Feedback of Efficient Shock Acceleration on Magnetic-field Structure Inside Young Type Ia Supernova Remnants

Using an effective adiabatic index γ eff to mimic the feedback of efficient shock acceleration, we simulate the temporal evolution of a young type Ia supernova remnant (SNR) with two different background magnetic field (BMF) topologies: a uniform and a turbulent BMF. The density distribution and magnetic-field characteristics of our benchmark SNR are studied with two-dimensional cylindrical magnetohydrodynamic simulations. When γ eff is considered, we find that: (1) the two-shock structure shrinks and the downstream magnetic-field orientation is dominated by the Rayleigh–Taylor instability structures; (2) there exists more quasi-radial magnetic fields inside the shocked region; and (3) inside the intershock region, both the quasi-radial magnetic energy density and the total magnetic energy density are enhanced: in the radial direction, with γ eff = 1.1, they are amplified about 10–26 times more than those with γ eff = 5/3. While in the angular direction, the total magnetic energy densities could be amplified about 350 times more than those with γ eff = 5/3, and there are more grid cells within the intershock region where the magnetic energy density is amplified by a factor greater than 100.

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

Accretion Funnel Reconfiguration during an Outburst in a Young Stellar Object: EX Lupi

EX Lupi, a low-mass young stellar object, went into an accretion-driven outburst in 2022 March. The outburst caused a sudden phase change of ∼112° ± 5° in periodically oscillating multiband lightcurves. Our high-resolution spectra obtained with the High Resolution Spectrograph (HRS) on board the Southern African Large Telescope also revealed a consistent phase change in the periodically varying radial velocities (RVs), along with an increase in the RV amplitude of various emission lines. The phase change and increase in RV amplitude morphologically translates to a change in the azimuthal and latitudinal location of the accretion hotspot over the stellar surface, which indicates a reconfiguration of the accretion funnel geometry. Our three-dimensional magnetohydrodynamic simulations reproduce the phase change for EX Lupi. To explain the observations, we explored the possibility of forward shifting of the dipolar accretion funnel as well as the possibility of the emergence of a new accretion funnel. During the outburst, we also found evidence of the hotspot’s morphology extending azimuthally asymmetrically with a leading hot edge and cold tail along the stellar rotation. Further, our high-cadence photometry showed that the accretion flow has clumps. We also detected possible clumpy accretion events in the HRS spectra that showed episodically highly blueshifted wings in the Ca ii IR triplet and Balmer H lines.

Radiation MHD Simulations of Soft X-Ray Emitting Regions in Changing Look AGN

Strong soft X-ray emission called soft X-ray excess is often observed in luminous active galactic nuclei (AGN). It has been suggested that the soft X-rays are emitted from a warm (T = 106 ∼ 107 K) region that is optically thick for the Thomson scattering (warm Comptonization region). Motivated by the recent observations that soft X-ray excess appears in changing look AGN (CLAGN) during the state transition from a dim state without broad emission lines to a bright state with broad emission lines, we performed global three-dimensional radiation magnetohydrodynamic simulations, assuming that the mass accretion rate increases and becomes around 10% of the Eddington accretion rate. The simulation successfully reproduces a warm, Thomson-thick region outside the hot radiatively inefficient accretion flow near the black hole. The warm region is formed by efficient radiative cooling due to inverse Compton scattering. The calculated luminosity 0.01−0.08 L Edd is consistent with the luminosity of CLAGN. We also found that the warm Comptonization region is well described by the steady model of magnetized disks supported by azimuthal magnetic fields. When the antiparallel azimuthal magnetic fields supporting the radiatively cooled region reconnect around the equatorial plane of the disk, the temperature of the region becomes higher by releasing the magnetic energy transported to the region.