Work with the complex representation and take its real part to obtain the physical magnetic field. For constant magnetic diffusivity and the prescribed circular velocity , the resistive induction equation and give
To see the last equation directly, . Taking its curl reproduces the induction equation if , up to a time-dependent additive constant removable by the gauge transformation. This is the scalar form of advection and magnetic diffusion in cylindrical coordinates.
Put . The material derivative of is , because . A radial derivative instead acts on both amplitude and phase:
Also . Substitution gives the shearing-coordinate magnetic flux equation, here for azimuthal mode number one:
Writing it with on the left is important when considering zero magnetic diffusivity: one must not literally divide by . In the ideal limit, , so . This is magnetic flux freezing: the scalar flux pattern is carried with each rotating fluid ring. It does not make the magnetic field stationary. Indeed,
so differential rotation winds the field and can make its azimuthal component grow linearly in time. Rigid rotation, for which , only rotates the pattern.
For the specified interior angular velocity, and away from the circulation edge. Let , and . The amplitude equation becomes
For a smooth initial amplitude on the cell scale, at bounded away from zero and the edge, the term dominates when . Its accumulated damping is order one at , whereas the integrated terms proportional to are only and the terms independent of are . The resulting cubic-time resistive damping in a differentially rotating cell is
The magnetic field has this leading exponential decay envelope. Its components also contain the algebraic winding factor visible in , so the claim concerns the cubic exponent, not an exact amplitude with no prefactor. Equivalently the radial wavenumber grows as and the integrated resistive damping is .
This is phase mixing in magnetic flux expulsion. The circulation creates ever finer gradients, allowing weak magnetic diffusion to act much sooner than its unsheared time :
Rapid winding followed by resistive smoothing expels the imposed field from the cell interior. The maintained exterior magnetic field is redistributed into a circulation-edge boundary layer. The result is an interior high-magnetic Reynolds number transient, not a uniform assertion at the axis and the nonsmooth edge, nor the exact solution with a continuously imposed far field. The diffusion length at is , identifying where a local bulk approximation can fail. The linear profile also needs smoothing at a regular rotation axis and at if used as a fully smooth global flow.
In photospheric magnetoconvection, circulating and diverging convection transports weak magnetic flux towards converging downflow regions. Solar granulation concentrates field in intergranular lanes, while larger-scale supergranulation helps form the solar magnetic network. Flux expulsion explains why a highly conducting photosphere can have relatively weak field in cell interiors and intermittent concentrations at their edges: high electrical conductivity delays diffusion on the original scale, but does not prevent diffusion on the fine scales created by advection. A long-lived circulation must last long enough for the expulsion time; rapidly changing cells need not reach this idealized state. Once the field becomes strong, the Lorentz force modifies the convection, so this prescribed-flow calculation alone neither determines the final field strength nor describes strong-field sunspots.
Solar granule 2026-10-14
A solar granule is a photospheric convection cell, with relatively hot rising material in its interior and cooler descending material near its edges. The pattern formed by their bright interiors and darker edges is solar granulation. Their evolving horizontal outflows can sweep weak magnetic flux toward surrounding lanes.
A motivating measurement used the McMath solar optical telescope at Kitt Peak. The Sun's image drifted across a -high, -wide optical slit; a double-pass spectrograph selected wavelength intervals from wide at to at . An EMI 9750 QB photodetector or silicon photodiode supplied a digitized signal: samples were summed into 4096 positions per crossing, including dark and off-limb sky readings. [ref]
Controls included daily realignment of optical telescope mirrors, wavelength settings between spectral lines insensitive to the rotational Doppler shift, selection of detectors without detectable hysteresis, and measured electronic nonlinearity corrections below of maximum signal. Limb-position calibration combined centres of gravity of the edge derivative profiles with almanac semidiameters; the preliminary edge shifts were consistent with astronomical seeing of approximately arcseconds full width at half maximum. [ref]
The decreasing centre-to-limb readout gives a specific intensity profile, constraining stratification of the photosphere under a radiative transfer model. Astronomical seeing, instrumental scattering of light, and solar granulation limit this inference; brightness alone does not uniquely determine temperature versus geometric height. These observations need not follow the stipulated linear limb-darkening law.