Impact of edge harmonic oscillations on the divertor heat flux in NSTX
All the data was uploaded with .cvs file, we have not uploaded the figure 1 data since it is just photo show field of view of IR and GPI diagnostic.
edge harmonic oscillation↗
SEARCH · Engineering Papers
Results for “heat flux width”
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.
Performance assessment of long-legged tightly-baffled divertor geometries in the ARC reactor concept
Extremely intense power exhaust channels are projected for tokamak-based fusion power reactors; a means to handle them remains to be demonstrated. Advanced divertor configurations have been proposed as potential solutions. Recent modelling of tightly baffled, long-legged divertor geometries for the divertor test tokamak concept, ADX, has shown that these concepts may access passively stable, fully detached regimes over a broad range of parameters. The question remains as to how such divertors may perform in a reactor setting. To explore this, numerical simulations are performed with UEDGE for the long-legged divertor geometry proposed for the ARC pilot plant conceptual design - a device with projected heat flux power width (λq||) of 0.4 mm and power exhaust of 93 MW - first for a simplified Super-X divertor configuration (SXD) and then for the actual X-point target divertor (XPTD) being proposed. It is found that the SXD, combined with 0.5% fixed-fraction neon impurity concentration, can produce passively stable, detached divertor regimes for power exhausts in the range of 80-108 MW - fully accommodating ARC's power exhaust. The XPTD configuration is found to reduce the strike-point temperature by a factor of ~10 compared to the SXD for small separations (~1.4λq||) between main and divertor X-point magnetic flux surfaces. Even greater potential reductions are identified for reducing separations to ~1λq|| or less. The power handling response is found to be insensitive to the level of cross-field convective or diffusive transport assumed in the divertor leg. By raising the separatrix density by a factor of 1.5, stable fully detached divertor solutions are obtained that fully accommodate the ARC exhaust power without impurity seeding. To our knowledge, this is the first time an impurity-free divertor power handling scenario has been obtained in edge modelling for a tokamak fusion power reactor with λq|| of 0.4 mm.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
New heat flux model for non-axisymmetric divertor infrared structures
A convective heat flux model for perturbed plasmas, based on guiding center ion drift in vacuum elds (A. Wingen, et al., Phys. Plasmas 21 (2014) 012509), has been updated. The old model only considered ion heat flux, while here also electron heat flux is included. The updated model predicts divertor heat flux distributions in non-axisymmetric (3D) plasmas with applied Resonant Magnetic Perturbation (RMP) fields, and includes electric scalar potentials. It is found that a radial electric field in the near Scrape-o Layer (SOL) can considerably shift the footprints toroidally, leading to a smearing out eff ect of the incident heat flux, while a simple model for sheath potential has little impact on footprints. Various approaches to model electron heat flux are studied. A convective electron model, based on collisionless free streaming, is found to yield the best agreement with measurements, while a conductive model requires a at temperature gradient inside lobes to yield acceptable peak heat ux values. A heuristic heat flux layer approach, based on a fixed layer width also requires a limited heat flux inside the last closed flux surface (LCFS); by selecting various locations of the LCFS, the results of the conductive or convective model can be recovered respectively. The sum of ion and electron heat fluxes, both obtained by the convective model, is compared to experimental data for multiple time slices in DIII-D. Strike point splitting is observed with peak heat fluxes and layer widths that compare well to infrared camera (IR) measurements.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effect of Lateral Conduction in a Rectangular Fuel Plate and Optimum Stripe Widths for the OPAL Reactor
The impetus for this work is the Proliferation Resistance Optimization (PRO-X) program, which was created by the National Nuclear Security Administration (NNSA) of the US Department of Energy (DOE) to provide a framework for developing reactor designs that minimize the production of special nuclear materials while maximizing performance for peaceful uses. The Argentine company INVAP and Argonne National Laboratory collaborated to share knowledge and expertise regarding the thermal analysis of research reactors. Working together, they designed a series of test problems based on the Open Pool Australian Light Water (OPAL) reactor, whose fuel is in the form of flat plates cooled by water flowing through narrow rectangular channels. ANL used their PLTEMP/ANL code. INVAP used their TERMIC code (and its new multi-plate version, TERMIC-MP). Then, INVAP and ANL compared the two sets of results. 1-D models were used during the collaboration since for research reactor analysis, largely 1-D thermal-hydraulics models are commonly used because they are simple to apply and technically defendable. The collaboration between the two organizations aims for improvements in thermal-hydraulic models. Improvements to methods and models can reduce modeling uncertainties and excessive conservatism, allowing greater reactor performance without reducing predicted safety margins. Therefore, these improvements potentially expand the design space, making additional designs feasible, which, in turn, may achieve better fuel utilization and proliferation resistance in the designs. This report provides justification for the lateral node size (stripe width) chosen for the 1-D computer models during the collaboration. The lateral heat conduction in OPAL fuel plates is evaluated using the two-dimensional conduction capability of PLTEMP/ANL to determine a representative hot-stripe width that can be conservatively used for one-dimensional steady-state thermal-hydraulic analysis. Basically, a representative hot-stripe width needs to be selected so that the actual heat flux averaged over this stripe width is greater than or equal to the maximum heat flux obtained by the two-dimensional analysis. Detailed results in this report show that a hot-stripe width of 8.125 mm can be conservatively used for the one-dimensional thermal-hydraulics analysis of all plates. Therefore, eight equal stripes (8.125 mm each) are recommended to be modeled in the MCNP and PLTEMP/ANL full core models of the OPAL reactor.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
The role of plasma–atom and molecule interactions on power & particle balance during detachment on the MAST Upgrade Super-X divertor
This paper shows first quantitative analysis of the detachment processes in the MAST Upgrade Super-X divertor (SXD). We identify an unprecedented impact of plasma-molecular interactions involving molecular ions (likely $D^+_2$), resulting in strong ion sinks (Molecular Activated Recombination—MAR), leading to a reduction of ion target flux. The MAR ion sinks exceed the divertor ion sources before electron-ion recombination (EIR) starts to occur, suggesting that significant ionisation occurs outside of the divertor chamber. In the EIR region, $T_e \ll 0.2$ eV is observed and MAR remains significant in these deep detached phases. The total ion sink strength demonstrates the capability for particle (ion) exhaust in the Super-X Configuration. Molecular Activated Dissociation is the dominant volumetric neutral atom creation process can lead to an electron cooling of 20% of $P_\mathrm{SOL}$. The measured total radiative power losses in the divertor chamber are consistent with inferred hydrogenic radiative power losses. This suggests that intrinsic divertor impurity radiation, despite the carbon walls, is minor in the divertor chamber. This contrasts previous TCV results, which may be associated with enhanced plasma-neutral interactions and reduced chemical erosion in the detached, tightly baffled SXD. The above observations have also been observed in higher heat flux (narrower SOL width) type I ELMy H-mode discharges. This provides evidence that the characterisation in this paper may be general.
plasma-molecular interactions↗
Reassessing energy deposition for the ITER 5 MA vertical displacement event with an improved DINA model
The beryllium (Be) main chamber wall interaction during a 5 MA/1.8 T upward, unmitigated VDE scenario, first analysed in [J. Coburn et al., Phys. Scr. T171 (2020) 014076] for ITER, has been re-evaluated using the latest energy deposition analysis software. Updates to the DINA disruption model are summarized, including an improved numerical convergence for the 0D power balance, limitations on the safety factor within the plasma core, and the choice to maintain a constant plasma +halo poloidal cross-section. Such updates result in a broad halo region and higher radiated power fractions compared to previous models. The new scenario lasts for ~75 ms and deposits ~29 MJ of energy, with the radial distribution of parallel heat flux q‖(r)resembling an exponential falloff with an effective λ_q=75-198 mm. A maximum halo width w_h of 0.52 m at the outboard midplane is observed. SMITER field line tracing and energy deposition simulations calculate a q_{⊥,max} of ~83 MW/m^2 on the upper first wall panels (FWP). Heat transfer calculations with the MEMOS-U code show that the FWP surface temperature reaches ~1000 K, well below the Be melt threshold. Variations of this 5 MA scenario with Be im-purity densities from 0 to 3∙10^19 m^-3 also remain below the melt threshold despite differences in energy deposition and duration. These results are in contrast to the early study which predicted melt damage to the first wall [J. Coburn et al., Phys. Scr. T171 (2020) 014076], and emphasize the importance of accurate models for the halo width w_h and the heat flux distribution q‖(r)within that halo width. The 2020 halo model in DINA has been compared with halo current experiments on COMPASS, JET, and Alcator C-Mod, and the preliminary results build confidence in the broad halo width predictions. Results for the 5 MA VDE are compared with those for a 15 MA equivalent, generated using the new DINA model. At the higher current, significant melting of the upper FWP is to be expected.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The Effects of Sensor Field-of-View on the Geometrical Characteristics of Sea Ice Leads and Implications for Large-Area Heat Flux Estimates
The release of heat from sea ice fractures ('leads') is an important component of the heat budget in the Arctic, but their impact on regional scale climate is difficult to assess without more information on their distribution in both space and time. Remote sensing of leads using satellite data, specifically AVHRR thermal and Landsat visible-band imagery, is examined empirically with respect to lead width, orientation, and area fraction. The geometrical aspects of the sensor are simulated so that the effect of sensor field-of-view on retrieved lead width statistics can be assessed. This is done using Landsat data and simulated lead networks degraded to AVHRR pixel sizes. The analyses illustrate how leads of sufficiently high contrast tend to 'grow' with increasing pixel size and how small or low contrast leads disappear. The relationship between lead contrast and the width/field-of-view ratio is also examined in order to determine the limits of lead detectability, and illustrates the multivalued nature of the problem of lead width retrieval. To help quantify the importance of changes in lead statistics, turbulent heat flux is calculated as a function of lead width and lead fraction. It is shown that pixel size has a substantial effect on estimates of turbulent heat transfer from leads to the atmosphere.
Key, J.↗
Turbulence driven widening of the near-SOL power width in ASDEX Upgrade H-Mode discharges
Divertor heat flux data from infra red (IR) from various tokamaks in H-Mode regime show that the power width, $\lambda_{q}$, scales approximately like 1.6 $(R/a) ρ_{s,pol}$. However, the IR based scaling comes with the restriction that only discharges were considered with attached divertor conditions. Here we set up a new data base covering low edge density data but also including high density plasmas reaching the H-mode density limit. We use Thomson-scattering to measure the electron temperature decay length which will set the near-SOL power width through parallel heat conduction. We present scaling results for the density, temperature and pressure separatrix fall-off lengths. As the principal result we derive a generalized power width scaling which reads as $\lambda_{q}\propto\,\rho_{s,pol}\cdot(1+2.1\alpha_{t}^{1.7})$ where $\alpha_{t}$ describes a normalized collisionality ($\alpha_{t} = 3\cdot\,10^{-18}\,R\,q^2\,n\,Z_{eff}T^{-2}$). The parameter $\alpha_{t}$ describes the relative importance of the interchange effect on drift-wave turbulence as proposed by Scott and is found to be closely linked to the diamagnetic parameter αd introduced by Rogers, Drake and Zeiler (RDZ). The edge operational space for ASDEX Upgrade displays a boundary qualitatively in agreement with the region labelled as density limit in the edge plasma phase space proposed by RDZ. The data base shows in the limit of low edge densities ($\alpha_{t}\simeq\,0.2$) accurate agreement to the IR based power width scaling. At elevated separatrix densities ($\alpha_{t}\simeq\,1$) the power width is broadened by a factor of up to about three accompanied by a reduction of the $\text{H}_{98,y2}$ confinement factor. Finally, plasmas with higher triangularity show higher values for $\text{H}_{98,y2}$ in the same range of normalized collisionality.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
On treatments of fetch and stability sensitivity in large-area estimates of sensible heat flux over sea ice
Bulk transfer coefficients estimated as a function of atmospheric stability and sea-ice lead width are combined with Arctic meteorological observations and ice thickness data to calculate the sensitivity of turbulent flux estimates to changes in lead width, wind speed, air temperature, and ice thickness for a high-concentration ice pack. These results are considered in terms of bulk transfer parameterizations that use a fixed transfer coefficient r that address atmospheric stability only. On the basis of the fetch-sensitive parameterizations considered here, differences in lead width for widths up to about 200 m can exert a substantial influence on sensible heat transfer coefficients and heat flux from leads under typical Arctic conditions. Fluxes from an open water lead decrease by 34% if fetch increases from 10 m to 100 m. This effect is greatest for open water leads, decreases considerably as leads refreeze, and is negligible for ice thicker than about 0.3 m. If open or newly refrozen leads make up 2% of the ice cover, than an increase in mean fetch from 10 m to 100 m yields a decrease of about 2 W/sq m in areally averaged flux from the ice pack. Calculations using observed and theoretical lead width distributions suggest that parameterizing lead widths in a sea ice model can be done effectively using a single, representative lead width rather than requiring a full distribution of widths. When coupled to the lowest atmospheric boundary layer using a bulk similarity theory model, this sensitivity of heat transfer to fetch results in substantially higher near-surface air temperatures over narrow leads, with equilibrium air temperatures decreasing by about 50% as fetch increases from 10 to 100 m.
Maslanik, J. A.↗
How fluctuation intensity flux drives SOL expansion
Abstract Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) magnetohydrodynamic (MHD) turbulence broadens the stable scrape-off layer (SOL) by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes. Pressure fluctuation skewness correlates well with the spreading flux–with the zero crossing of skewness and Γ ε spatially coincident–suggesting the role of coherent fluctuation structures and the presence of intermittency in λ q broadening. λ q ∼ B p − 1 scaling persists for the broadened SOL. We show that the spreading flux increases for increasing pedestal pressure gradient ∇ P 0 and for decreasing pedestal collisionality υ ped ∗ . This trend is due to the dominance of peeling modes for large ∇ P 0 and low υ ped ∗ . Ultimately, we see that a state of weak MHD turbulence, as for small ELMs, is very attractive for heat load management. Our findings have transformative implications for future fusion reactor designs and call for experimental investigations to validate the observed trends.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Dependence of Deep Convective Cell Properties on Meteorological and Aerosol Conditions during TRACER
Deep convective cells significantly influence Earth’s energy balance and water cycle. However, their accurate representation in numerical models remains challenging due to their small spatiotemporal scales and limited observational constraints. This study examines over ∼400 deep convective cells near Houston, observed by a dual-polarization C-band radar during the Tracking Aerosol Convection Interactions Experiment (TRACER) intensive observation period (June–September 2022). Cells are categorized by lifetime into short-lived (<40 min), intermediate-lived (40–80 min), and long-lived (80+ min) groups. Long-lived cells were broader (∼13.2 km at 2–4-km height) and deeper (∼11.4 km) than short-lived cells (∼6.4-km width, ∼7.31-km height). Using random forest (RF) modeling and correlation analyses, precipitable water vapor (PWV), 2–6-km lapse rate, 0–8-km bulk shear, and fine aerosol mass concentration (Mass_f) are identified as key predictors of cell lifetime. Higher PWV is associated with significantly longer convective cell lifetimes compared to the low-PWV group, particularly within low 2–6-km temperature lapse rate (LR_26km), moderate-to-higher 0–8-km bulk shear (BS_08km), and low-to-moderate Mass_f environments. RF analysis also identifies low-level (0–2 km) equivalent potential temperature, PWV, Mass_f, and surface latent heat flux as key predictors for cell width and height. Short-lived cells have higher aerosol number concentrations (500–1000-nm size range), linked to onshore wind conditions and marine aerosols; however, their low concentration suggests the sensitivity may reflect associated meteorological regimes rather than a direct aerosol effect. Long-lived cells have higher concentrations of organic and sulfate aerosols, while short-lived cells exhibit higher black carbon concentrations. These results highlight the intricate dependence of convective cell lifetimes and structure on environmental moisture, thermodynamics, wind shear, and aerosol characteristics.
54 ENVIRONMENTAL SCIENCES↗
Gyrokinetic prediction of core tungsten peaking in a WEST plasma with nitrogen impurities
Tungsten peaking is predicted in the core of a WEST plasma with total-f gyrokinetic simulations, including both collisional and turbulent transport. This prediction is validated with a synthetic diagnostic of the bolometry. Although nitrogen impurities are shown to reduce the neoclassical peaking of tungsten on-axis, the overall tungsten peaking increases when nitrogen impurities are present, as they reduce the turbulence screening off-axis. This finding is important for the plasma current ramp-up phase of ITER, where light impurities seeding will be desirable to achieve low temperatures at the plasma-facing components and reduce tungsten sputtering. It provides further argument for applying early ECRH heating to maintain margins on the core power balance. The neoclassical peaking factor is cross-verified between XGC and FACIT. The heat flux at separatrix and the heat load width are modeled by XGC and compared to WEST data.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Effects of compression and expansion ramp fuel injector configuration on scramjet combustion and heat transfer
A scramjet combustor with four wall-ramp injectors containing Mach-1.7 fuel jets in the base of the ramps was investigated experimentally. During the test program, two swept ramp injector designs were evaluated. One swept-ramp model had 10-deg compression-ramps and the other had 10-deg expansion cavities between flush wall ramps. The scramjet combustor model was instrumented with pressure taps and heat-flux gages. The pressure measurements indicated that both injector configurations were effective in promoting mixing and combustion. Autoignition occurred for the compression-ramp injectors, and the fuel began to burn immediately downstream of the injectors. In tests of the expansion ramps, a pilot was required to ignite the fuel, and the fuel did not burn for a distance of at least two gaps downstream of the injectors. Once initiated, combustion was rapid in this configuration. Heat transfer measurements showed that the heat flux differed greatly both across the width of the combustor and along the length of the combustor.
Stouffer, Scott D.↗
Parameterized Study of Heat Load Trends in a Subscale Rotating Detonation Rocket Engine
NASA has developed a subscale rotating detonation rocket engine platform, enabling rapid parameterization of components and their associated performances. A major performance metric of interest for RDREs is the total heat absorbed at a given operating condition and wave mode operation. To investigate this, several sets of hardware were produced to have variation of the contraction ratio, injector geometry, and length. A wide range of total mass flow rates and mixture ratios were also explored using gaseous methane/oxygen. All hardware was produced using laser powder bed fusion GRCop-42, GRX-810, or C-103 alloys depending on the component cooling requirements. Outer and inner body hardware were made from GRCop-42 where heat fluxes were expected to be high and were water cooled through integrated coolant channels. Nozzles were made from C-103 which are radiatively cooled. Finally, the injectors were produced using GRCop-42 or GRX-810 as passive cooling via propellant injection has previously found to be sufficient given the small gap width of exposed injection surface. Total heat load and bulk heat flux measurements to the hot walls are reported. Trends in heat load and heat flux are shown with relation to parameters such as chamber pressure, area ratio, and wave mode. Heat loads between the inner and outer bodies and the chamber and shroud section are compared.
Heat Flux↗
Parameterized Study of Heat Load Trends in a Subscale Rotating Detonation Rocket Engine
NASA has developed a subscale rotating detonation rocket engine platform, enabling rapid parameterization of components and their associated performances. A major performance metric of interest for RDREs is the total heat absorbed at a given operating condition and wave mode operation. To investigate this, several sets of hardware were produced to have variation of the contraction ratio, injector geometry, and length. A wide range of total mass flow rates and mixture ratios were also explored using gaseous methane/oxygen. All hardware was produced using laser powder bed fusion GRCop-42, GRX-810, or C-103 alloys depending on the component cooling requirements. Outer and inner body hardware were made from GRCop-42 where heat fluxes were expected to be high and were water cooled through integrated coolant channels. Nozzles were made from C-103 which are radiatively cooled. Finally, the injectors were produced using GRCop-42 or GRX-810 as passive cooling via propellant injection has previously found to be sufficient given the small gap width of exposed injection surface. Total heat load and bulk heat flux measurements to the hot walls are reported. Trends in heat load and heat flux are shown with relation to parameters such as chamber pressure, area ratio, and wave mode. Heat loads between the inner and outer bodies and the chamber and shroud section are compared.
Heat Loads↗
Nonlinear axially symmetric circulations in a nearly inviscid atmosphere
The structure of certain axially symmetric circulations in a stably stratified, differentially heated, rotating Boussinesq fluid on a sphere is analyzed. A simple approximate theory (similar to that introduced by Schneider (1977)) is developed for the case in which the fluid is sufficiently inviscid that the poleward flow in the Hadley cell is nearly angular momentum conserving. The theory predicts the width of the Hadley cell, the total poleward heat flux, the latitude of the upper level jet in the zonal wind, and the distribution of surface easterlies and westerlies. Fundamental differences between such nearly inviscid circulations and the more commonly studied viscous axisymmetric flows are emphasized. The theory is checked against numerical solutions to the model equations.
Held, I. M.↗
Solar Flare Heating with Turbulent Suppression of Thermal Conduction
During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppmission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function.
Joel C. Allred↗
Solar Flare Heating with Turbulent Suppression of Thermal Conduction
During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function
Solar flares↗