BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Bradley W. Hindman

Bradley W. Hindman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-phastro-ph.EPphysics.flu-dyn

Catalog footprint

What is connected

18works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2024arXiv

Confinement of the Solar Tachocline by a Non-Axisymmetric Dynamo

We recently presented the first 3D numerical simulation of the solar interior for which tachocline confinement was achieved by a dynamo-generated magnetic field. In this followup study, we analyze the degree of confinement as the magnetic field strength changes (controlled by varying the magnetic Prandtl number) in a coupled radiative zone (RZ) and convection zone (CZ) system. We broadly find three solution regimes, corresponding to weak, medium, and strong dynamo magnetic field strengths. In the weak-field regime, the large-scale magnetic field is mostly axisymmetric with regular, periodic polarity reversals (reminiscent of the observed solar cycle), but fails to create a confined tachocline. In the strong-field regime, the large-scale field is mostly non-axisymmetric with irregular, quasi-periodic polarity reversals, and creates a confined tachocline. In the medium-field regime, the large-scale field resembles a strong-field dynamo for extended intervals, but intermittently weakens to allow temporary epochs of strong differential rotation. In all regimes, the amplitude of poloidal field strength in the RZ is very well explained by skin-depth arguments, wherein the oscillating field that gives rise to the skin depth (in the medium- and strong-field cases) is a non-axisymmetric field structure rotating with respect to the RZ. These new simulations reaffirm that tachocline confinement by the solar dynamo (the so-called fast magnetic confinement scenario) is possible, but suggest a new picture in which non-axisymmetric field components rotating with respect to the RZ play the primary role, instead of the regularly reversing axisymmetic field associated with the 22-year cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraining Global Solar Models through Helioseismic Analysis

Global hydrodynamic simulations of internal solar dynamics have focused on replicating the conditions for solar-like differential rotation and meridional circulation using the results of helioseismic inversions as a constraint. Inferences of meridional circulation, however, have provided controversial results showing the possibility of one, two, or multiple cells along the radius. To resolve this controversy and develop a more robust understanding of global flow regimes in the solar interior, we apply a "forward-modeling" approach to the analysis of helioseismic signatures of meridional circulation profiles obtained from numerical simulations. We employ the global acoustic modeling code GALE to simulate the propagation of acoustic waves through regimes of mean mass flows generated by global hydrodynamic and magnetohydrodynamic models: EULAG, the Pencil Code, and the Rayleigh code. These models are used to create synthetic dopplergram data products, used as inputs for local time-distance helioseismology techniques. Helioseismic travel-time signals from solutions obtained through global numerical simulations are compared directly with inferences from solar observations, in order to set additional constraints on global model parameters in a direct way. We show that even though these models are able to replicate solar-like differential rotation, the resulting rotationally-constrained convection develops a multi-cell global meridional circulation profile that is measurably inconsistent with local time-distance inferences of solar observations. However, we find that the development of rotationally-unconstrained convection close to the model surface is able to maintain solar-like differential rotation, while having a significant impact on the helioseismic travel-time signal, replicating solar observations within one standard deviation of the error due to noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revisiting the Sun's Strong Differential Rotation along Radial Lines

Current state-of-the-art models of the solar convection zone consist of solutions to the Navier-Stokes equations in rotating, 3D spherical shells. Such models are highly sensitive to the choice of boundary conditions. Here, we present two suites of simulations differing only in their outer thermal boundary condition, which is either one of fixed-entropy or fixed-entropy-gradient. We find that the resulting differential rotation is markedly different between the two sets. The fixed-entropy-gradient simulations have strong differential rotation contrast and isocontours tilted along radial lines (in good agreement with the Sun's interior rotation revealed by helioseismology), whereas the fixed-entropy simulations have weaker contrast and contours tilted in the opposite sense. We examine in detail the force balances in our models and find that the poleward transport of heat by Busse columns drives a thermal wind responsible for the different rotation profiles. We conclude that the Sun's strong differential rotation along radial lines may result from the solar emissivity being invariant with latitude (which is similar to the fixed-entropy-gradient condition in our models) and the poleward transport of heat by Busse columns. In future work on convection in the solar context, we strongly advise modelers to use a fixed-gradient outer boundary condition.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An improved multi-ridge fitting method for ring-diagram helioseismic analysis

Context: There is a wide discrepancy in current estimates of the strength of convection flows in the solar interior obtained using different helioseismic methods applied to observations from SDO/HMI. The cause for these disparities is not known. Aims: As one step in the effort to resolve this discrepancy, we aim to characterize the multi-ridge fitting code for ring-diagram helioseismic analysis that is used to obtain flow estimates from local power spectra of solar oscillations. Methods: We updated the multi-ridge fitting code developed by Greer et al.(2014) to solve several problems we identified through our inspection of the code. In particular, we changed the merit function to account for the smoothing of the power spectra, model for the power spectrum, and noise estimates. We used Monte Carlo simulations to generate synthetic data and to characterize the noise and bias of the updated code by fitting these synthetic data. Results: The bias in the output fit parameters, apart from the parameter describing the amplitude of the p-mode resonances in the power spectrum, is below what can be measured from the Monte-Carlo simulations. The amplitude parameters are underestimated; this is a consequence of choosing to fit the logarithm of the averaged power. We defer fixing this problem as it is well understood and not significant for measuring flows in the solar interior. The scatter in the fit parameters from the Monte-Carlo simulations is well-modeled by the formal error estimates from the code. Conclusions: We document and demonstrate a reliable multi-ridge fitting method for ring-diagram analysis. The differences between the updated fitting results and the original results are less than one order of magnitude and therefore we suspect that the changes will not eliminate the aforementioned orders-of-magnitude discrepancy in the amplitude of convective flows in the solar interior.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Emergence of Solar Supergranulation as a Natural Consequence of Rotationally-Constrained Interior Convection

We investigate how rotationally-constrained, deep convection might give rise to supergranulation, the largest distinct spatial scale of convection observed in the solar photosphere. While supergranulation is only weakly influenced by rotation, larger spatial scales of convection sample the deep convection zone and are presumably rotationally influenced. We present numerical results from a series of nonlinear, 3-D simulations of rotating convection and examine the velocity power distribution realized under a range of Rossby numbers. When rotation is present, the convective power distribution possesses a pronounced peak, at characteristic wavenumber $\ell_\mathrm{peak}$, whose value increases as the Rossby number is decreased. This distribution of power contrasts with that realized in non-rotating convection, where power increases monotonically from high to low wavenumbers. We find that spatial scales smaller than $\ell_\mathrm{peak}$ behave in analogy to non-rotating convection. Spatial scales larger than $\ell_\mathrm{peak}$ are rotationally-constrained and possess substantially reduced power relative to the non-rotating system. We argue that the supergranular scale emerges due to a suppression of power on spatial scales larger than $\ell\approx100$ owing to the presence of deep, rotationally-constrained convection. Supergranulation thus represents the largest non-rotationally-constrained mode of solar convection. We conclude that the characteristic spatial scale of supergranulation bounds that of the deep convective motions from above, making supergranulation an indirect measure of the deep-seated dynamics at work in the solar dynamo. Using the spatial scale of supergranulation in conjunction with our numerical results, we estimate an upper bound of 10 m s$^{-1}$ for the Sun's bulk rms convective velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Helioseismic Imaging of Fast Convective Flows Throughout the Near-Surface Shear Layer

Using a new implementation of ring-diagram helioseismology, we ascertain the strength and spatial scale of convective flows throughout the near-surface shear layer. Our ring-diagram technique employs highly overlapped analysis regions and an efficient method of 3D inversion to measure convective motions with a resolution that ranges from $3 \ \mathrm{Mm}$ at the surface to $80 \ \mathrm{Mm}$ at the base of the layer. We find the rms horizontal flow speed to peak at $427 \ \mathrm{m \ s^{-1}}$ at the photosphere and fall to a minimum of $124 \ \mathrm{m \ s^{-1}}$ between $20 \ \mathrm{Mm}$ and $30 \ \mathrm{Mm}$. From the velocity amplitude and the dominant horizontal scales seen at each depth, we infer the level of rotational influence on convection to be low near the surface, but transition to a significant level at the base of the near-surface shear layer with a Rossby number varying between 2.2 to as low as 0.1.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Spectral Amplitude of Stellar Convection and its Scaling in the High-Rayleigh-Number Regime

Convection plays a central role in the dynamics of any stellar interior, and yet its operation remains largely-hidden from direct observation. As a result, much of our understanding concerning stellar convection necessarily derives from theoretical and computational models. The Sun is, however, exceptional in that regard. The wealth of observational data afforded by its proximity provides a unique testbed for comparing convection models against observations. When such comparisons are carried out, surprising inconsistencies between those models and observations become apparent. Both photospheric and helioseismic measurements suggest that convection simulations may overestimate convective flow speeds on large spatial scales. Moreover, many solar convection simulations have difficulty reproducing the observed solar differential rotation due to this apparent overestimation. We present a series of 3-dimensional (3-D) stellar convection simulations designed to examine how the amplitude and spectral distribution of convective flows are established within a star's interior. While these simulations are non-magnetic and non-rotating in nature, they demonstrate two robust phenomena. When run with sufficiently high Rayleigh number, the integrated kinetic energy of the convection becomes effectively independent of thermal diffusion, but the spectral distribution of that kinetic energy remains sensitive to both of these quantities. A simulation that has converged to a diffusion-independent value of kinetic energy will divide that energy between spatial scales such that low-wavenumber power is overestimated, and high-wavenumber power is underestimated relative to a comparable system possessing higher Rayleigh number. We discuss the implications of these results in light of the current inconsistencies between models and observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Five-minute oscillation power within magnetic elements in the solar atmosphere

It has long been known that magnetic plage and sunspots are regions in which the power of acoustic waves is reduced within the photospheric layers. Recent observations now suggest that this suppression of power extends into the low chromosphere and is also present in small magnetic elements far from active regions. In this paper we investigate the observed power supression in plage and magnetic elements, by modelling each as a collection of vertically aligned magnetic fibrils and presuming that the velocity within each fibril is the response to buffeting by incident $p$ modes in the surrounding field-free atmosphere. We restrict our attention to modeling observations made near solar disk center, where the line-of-sight velocity is nearly vertical and hence, only the longitudinal component of the motion within the fibril contributes. Therefore, we only consider the excitation of axisymmetric sausage waves and ignore kink oscillations as their motions are primarily horizontal. We compare the vertical motion within the fibril with the vertical motion of the incident $p$ mode by constructing the ratio of their powers. In agreement with observational measurements we find that the total power is suppressed within strong magnetic elements for frequencies below the acoustic cut-off frequency. We also find that the magnitude of the power deficit increases with the height above the photosphere at which the measurement is made. Further, we argue that the area of the solar disk over which the power suppression extends increases as a function of height.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multi-Ridge Fitting for Ring-Diagram Helioseismology

Inferences of sub-surface flow velocities using local domain ring-diagram helioseismology depend on measuring the frequency splittings of oscillation modes seen in acoustic power spectra. Current methods for making these measurements utilize maximum-likelihood fitting techniques to match a model of modal power to the spectra. The model typically describes a single oscillation mode, and each mode in a given power spectrum is fit independently. We present a new method that produces measurements with greater reliability and accuracy by fitting multiple modes simultaneously. We demonstrate how this method permits measurements of sub-surface flows deeper into the Sun while providing higher uniformity in data coverage and velocity response closer to the limb of the solar disk. While the previous fitting method performs better for some measurements of low-phase-speed modes, we find this new method to be particularly useful for high phase-speed modes and small spatial areas.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An Interpretation of Flare-Induced and Decayless Coronal-Loop Oscillations as Interference Patterns

We present an alternative model of coronal-loop oscillations which considers that the waves are trapped in a 2D waveguide formed by the entire arcade of field lines. This differs from the standard 1D model which treats the waves as the resonant oscillations of just the visible bundle of field lines. Within the framework of our 2D model, the two types of oscillations that have been observationally identified, flare-induced waves and "decayless" oscillations, can both be attributed to MHD fast waves. The two components of the signal differ only because of the duration and spatial extent of the source that creates them. The flare-induced waves are generated by strong localized sources of short duration, while the decayless background can be excited by a continuous, stochastic source. Further, the oscillatory signal arising from a localized, short-duration source can be interpreted as a pattern of interference fringes produced by waves that have traveled diverse routes of various pathlengths through the waveguide. The resulting amplitude of the fringes slowly decays in time with an inverse square root dependence. The details of the interference pattern depend on the shape of the arcade and the spatial variation of the Alfvén speed. The rapid decay of this wave component, which has previously been attributed to physical damping mechanisms that remove energy from resonant oscillations, occurs as a natural consequence of the interference process without the need for local dissipation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Equilibrium models of coronal loops that involve curvature and buoyancy

We construct magnetostatic models of coronal loops in which the thermodynamics of the loop is fully consistent with the shape and geometry of the loop. This is achieved by treating the loop as a thin, compact, magnetic fibril that is a small departure from a force-free state. The density along the loop is related to the loop's curvature by requiring that the Lorentz force arising from this deviation is balanced by buoyancy. This equilibrium, coupled with hydrostatic balance and the ideal gas law, then connects the temperature of the loop with the curvature of the loop without resorting to a detailed treatment of heating and cooling. We present two example solutions: one with a spatially invariant magnetic Bond number (the dimensionless ratio of buoyancy to Lorentz forces) and the other with a constant radius of curvature of the loop's axis. We find that the density and temperature profiles are quite sensitive to curvature variations along the loop, even for loops with similar aspect ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kink Oscillations of a Curved, Gravitationally Stratified Coronal Loop

Loops of magnetic field in the corona are observed to oscillate and these oscillations have been posited to be the superposition of resonant kink waves. To date, most analyses of these oscillations have concentrated on calculating the frequency shifts that result from spatial variation in the kink wave speed. Further, most have ignored gravity and treated the loop as a straight tube. Here we ignore spatial variation in the wave speed, but self-consistently include the effects of gravity and loop curvature in both the equilibrium loop model and in the wave equation. We model a coronal loop as an isolated, thin, magnetic fibril that is anchored at two points in the photosphere. The equilibrium shape of the loop is determined by a balance between magnetic buoyancy and magnetic tension, which is characterized by a Magnetic Bond Number ε, that is typically small |ε| << 1. This balance produces a loop that has a variable radius of curvature. The resonant kink waves of such a loop come in two polarizations that are decoupled from each other: waves with motion completely within the plane of the loop (normal oscillations) and waves with motions that are completely horizontal, perpendicular to the plane of the loop (binormal oscillations). We solve for the eigensolutions of both polarizations using perturbation theory for small Magnetic Bond Number. For modes of the same order, normal oscillations have smaller eigenfrequencies than binormal oscillations. The additional forces of buoyancy and magnetic tension from the curvature of the loop increase and decrease the mode frequencies, respectively. The ratio of the frequencies of the first overtone to the fundamental mode is modified by the inclusion of buoyancy and curvature. We find that the normal polarization possesses a frequency ratio that exceeds the canonical value of 2, whereas the binormal polarization has a ratio less than 2.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Improved Method for Fitting p-Mode Profile Asymmetries

In a power spectrum of the Sun's acoustic waves, the p modes have distinctly skewed frequency profiles. Furthermore, the asymmetry is observed to have the opposite sign in power spectra made from line-of-sight velocity and continuum intensity. The asymmetry and its reversal in sign has previously been explained using a combination of mechanisms that involve the acoustic source. A localized acoustic source within an acoustic cavity naturally generates asymmetric profiles through wave interference; however, the sign of the asymmetry due to this mechanism is identical for all observables. The reversal of the asymmetry between velocity and intensity observations has been attributed to the visibility of the source itself (i.e., "correlated noise"). In this paper, I will show that asymmetry generated by a localized source can be interpreted as either a wave interference effect in physical space, or a mode interference effect in spectral space. I advocate a new mode-fitting procedure based on this new interpretation, whereby the complex phases of all the modes determine the mode asymmetries. Further, I suggest that information about the acoustic source function is encapsulated in the amplitude of each mode, and present a scheme by which the source function can be obtained from measured mode amplitudes by standard helioseismic inversion techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Axisymmetric Scattering of p Modes by Thin Magnetic Tubes

We examine the scattering of acoustic p-mode waves from a thin magnetic fibril embedded in a gravitationally stratified atmosphere. The scattering is mediated through the excitation of slow sausage waves on the magnetic tube, and only the scattering of the monopole component of the wavefield is considered. Since such tube waves are not confined by the acoustic cavity and may freely propagate along the field lines removing energy from the acoustic wavefield, the excitation of fibril oscillations is a source of acoustic wave absorption as well as scattering. We compute the mode mixing that is achieved and the absorption coefficients and phase shifts. We find that for thin tubes the mode mixing is weak and the absorption coefficient is small and is a smooth function of frequency over the physically relevant band of observed frequencies. The prominent absorption resonances seen in previous studies of unstratified tubes are absent. Despite the relatively small absorption, the phase shift induced can be surprisingly large, reaching values as high as 15 degree for f modes. Further, the phase shift can be positive or negative depending on the incident mode order and the frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gyroscopic Pumping in the Solar Near-Surface Shear Layer

We use global and local helioseismic inversions to explore the prevailing dynamical balances in the solar Near-Surface Shear Layer (NSSL). The differential rotation and meridional circulation are intimately linked, with a common origin in the turbulent stresses of the upper solar convection zone. The existence and structure of the NSSL cannot be attributed to the conservation of angular momentum by solar surface convection, as is often supposed. Rather, the turbulent angular momentum transport accounts for the poleward meridional flow while the often overlooked meridional force balance maintains the mid-latitude rotational shear. We suggest that the base of the NSSL is marked by a transition from baroclinic to turbulent stresses in the meridional plane which suppress Coriolis-induced circulations that would otherwise establish a cylindrical rotation profile. The turbulent angular momentum transport must be non-diffusive and directed radially inward. Inferred mean flows are consistent with the idea that turbulent convection tends to mix angular momentum but only if the mixing efficiency is inhomogeneous and/or anisotropic. The latitudinal and longitudinal components of the estimated turbulent transport are comparable in amplitude and about an order of magnitude larger than the vertical component. We estimate that it requires 2--4% of the solar luminosity to maintain the solar NSSL against the inertia of the mean flow. Most of this energy is associated with the turbulent transport of angular momentum out of the layer, with a spin-down time scale of $\sim$ 600 days. We also address implications of these results for numerical modeling of the NSSL.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Helioseismic detection of deep meridional flow

Steady meridional flow makes no first-order perturbation to the frequencies of helioseismic normal modes. It does, however, Doppler shift the local wavenumber, thereby distorting the eigenfunctions. For high-degree modes, whose peaks in a power spectrum are blended into continuous ridges, the effect of the distortion is to shift the locations of those ridges. From this blended superposition of modes, one can isolate oppositely directed wave components with the same local horizontal wavenumber and measure a frequency difference which can be safely used to infer the subsurface background flow. But such a procedure fails for the components of the more-deeply-penetrating low-degree modes that are not blended into ridges. Instead, one must analyze the spatial distortions explicitly. With a simple toy model, we illustrate one method by which that might be accomplished by measuring the spatial variation of the oscillation phase. We estimate that by this procedure it might be possible to infer meridional flow deep in the solar convection zone.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Subsurface Circulations within Active Regions

Using high-resolution ring analysis (HRRA) we deduce subsurface flows within magnetic active regions and within quiet sun. With this procedure we are capable of measuring flows with a horizontal spatial resolution of 2 degrees in heliographic angle (or roughly 20 Mm). From the resulting flow fields we deduce mean inflow rates into active regions, mean circulation speeds around active regions, and probability density functions (PDFs) of properties of the flow field. These analyses indicate that active regions have a zonal velocity that exceeds that of quiet sun at the same latitude by 20 m/s, yet active regions advect poleward at the same rate as quiet sun. We also find that almost all active regions possess a mean inflow (20-30 m/s) and a cyclonic circulation (~5 m/s) at their peripheries, whereas their cores, where the sunspots are located, are zones of strong anticyclonic outflow (~50 m/s). From the PDFs, we find that active regions modify the structure of convection with a scale greater than that of supergranulation. Instead of possessing an asymmetry between inflows and outflows (with a larger percentage of the surface occupied by outflows), as is seen in quiet sun, active regions possess symmetric distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Generation of Coronal Loop Waves below the Photosphere by p-Mode Forcing

Recent observations of coronal-loop waves by TRACE and within the corona as a whole by CoMP clearly indicate that the dominant oscillation period is 5 minutes, thus implicating the solar p modes as a possible source. We investigate the generation of tube waves within the solar convection zone by the buffeting of p modes. The tube waves--in the form of longitudinal sausage waves and transverse kink waves--are generated on the many magnetic fibrils that lace the convection zone and pierce the solar photosphere. Once generated by p-mode forcing, the tube waves freely propagate up and down the tubes, since the tubes act like light fibers and form a waveguide for these magnetosonic waves. Those waves that propagate upward pass through the photosphere and enter the upper atmosphere where they can be measured as loop oscillations and other forms of propagating coronal waves. We treat the magnetic fibrils as vertically aligned, thin flux tubes and compute the energy flux of tube waves that can generated and driven into the upper atmosphere. We find that a flux in excess of 10^5 ergs/cm^2/s can be produced, easily supplying enough wave energy to explain the observations. Furthermore, we compute the associated damping rate of the driving p modes and find that the damping is significant compared to observed line widths only for the lowest order p modes.

0 citations0 reviewsNo code linkedOpen access