Numerical simulations of multiple scattering of the f−f-mode by flux tubes

We use numerial simulations to study the absorption and phase shift of surface-gravity waves caused by groups of magnetic flux tubes. The dependence of the scattering coefficients with the distance between the tubes and their positions is analyzed for several cases with two or three flux tubes embedded in a quiet Sun atmosphere. The results are compared with those obtained neglecting completely or partially multiple scattering effects. We show that multiple scattering has a significant impact on the absorption measurements and tends to reduce the phase shift. We also consider more general cases of ensembles of randomly distributed flux tubes, and we have evaluated the effects on the scattering measurements of changing the number of tubes included in the bundle and the average distance between flux tubes. We find that for the longest wavelength incoming waves multiple scattering enhances the absorption, and its efficiency increases with the number of flux tubes and the reduction of the distance between them.Comment: Accepted for publication in The Astrophysical Journa

Evaluation of the capability of local helioseismology to discern between monolithic and spaghetti sunspot models

The helioseismic properties of the wave scattering generated by monolithic and spaghetti sunspots are analyzed by means of numerical simulations. In these computations, an incident f or p1 mode travels through the sunspot model, which produces absorption and phase shift of the waves. The scattering is studied by inspecting the wavefield, computing travel-time shifts, and performing Fourier-Hankel analysis. The comparison between the results obtained for both sunspot models reveals that the differences in the absorption coefficient can be detected above noise level. The spaghetti model produces an steep increase of the phase shift with the degree of the mode at short wavelengths, while mode-mixing is more efficient for the monolithic model. These results provide a clue for what to look for in solar observations to discern the constitution of sunspots between the proposed monolithic and spaghetti models.Comment: Accepted for publication in The Astrophysical Journa

Helioseismic holography of simulated sunspots: magnetic and thermal contributions to travel times

Wave propagation through sunspots involves conversion between waves of acoustic and magnetic character. In addition, the thermal structure of sunspots is very different than that of the quiet Sun. As a consequence, the interpretation of local helioseismic measurements of sunspots has long been a challenge. With the aim of understanding these measurements, we carry out numerical simulations of wave propagation through sunspots. Helioseismic holography measurements made from the resulting simulated wavefields show qualitative agreement with observations of real sunspots. We use additional numerical experiments to determine, separately, the influence of the thermal structure of the sunspot and the direct effect of the sunspot magnetic field. We use the ray approximation to show that the travel-time shifts in the thermal (non-magnetic) sunspot model are primarily produced by changes in the wave path due to the Wilson depression rather than variations in the wave speed. This shows that inversions for the subsurface structure of sunspots must account for local changes in the density. In some ranges of horizontal phase speed and frequency there is agreement (within the noise level in the simulations) between the travel times measured in the full magnetic sunspot model and the thermal model. If this conclusion proves to be robust for a wide range of models, it would suggest a path towards inversions for sunspot structure.Comment: Accepted for publication in The Astrophysical Journa

Physical Baryon Resonance Spectroscopy from Lattice QCD

We complement recent advances in the calculation of the masses of excited baryons in quenched lattice QCD with finite-range regulated chiral effective field theory enabling contact with the physical quark mass region. We examine the P-wave contributions to the low-lying nucleon and delta resonances.Comment: Contributed paper at FB17, the 17th International Conference on Few-Body Problems in Physics, Durham, NC, June 5-10, 2003. 3 pages, 6 figure

Validating Forward Modeling and Inversions of Helioseismic Holography Measurements

Here we use synthetic data to explore the performance of forward models and inverse methods for helioseismic holography. Specifically, this work presents the first comprehensive test of inverse modeling for flows using lateral-vantage (deep-focus) holography. We derive sensitivity functions in the Born approximation. We then use these sensitivity functions in a series of forward models and inversions of flows from a publicly available magnetohydrodynamic quiet-Sun simulation. The forward travel times computed using the kernels generally compare favorably with measurements obtained by applying holography, in a lateral-vantage configuration, on a 15-hour time series of artificial Dopplergrams extracted from the simulation. Inversions for the horizontal flow components are able to reproduce the flows in the upper 3Mm of the domain, but are compromised by noise at greater depths.Comment: accepted for publication by the Astrophysical

An Optimal Control Formulation for Inviscid Incompressible Ideal Fluid Flow

In this paper we consider the Hamiltonian formulation of the equations of incompressible ideal fluid flow from the point of view of optimal control theory. The equations are compared to the finite symmetric rigid body equations analyzed earlier by the authors. We discuss various aspects of the Hamiltonian structure of the Euler equations and show in particular that the optimal control approach leads to a standard formulation of the Euler equations -- the so-called impulse equations in their Lagrangian form. We discuss various other aspects of the Euler equations from a pedagogical point of view. We show that the Hamiltonian in the maximum principle is given by the pairing of the Eulerian impulse density with the velocity. We provide a comparative discussion of the flow equations in their Eulerian and Lagrangian form and describe how these forms occur naturally in the context of optimal control. We demonstrate that the extremal equations corresponding to the optimal control problem for the flow have a natural canonical symplectic structure.Comment: 6 pages, no figures. To appear in Proceedings of the 39th IEEEE Conference on Decision and Contro

Scattering of the f-mode by small magnetic flux elements from observations and numerical simulations

The scattering of f-modes by magnetic tubes is analyzed using three-dimensional numerical simulations. An f-mode wave packet is propagated through a solar atmosphere embedded with three different flux tube models which differ in radius and total magnetic flux. A quiet Sun simulation without a tube present is also performed as a reference. Waves are excited inside the flux tube and propagate along the field lines, and jacket modes are generated in the surroundings of the flux tube, carrying 40% as much energy as the tube modes. The resulting scattered wave is mainly an f-mode composed of a mixture of m=0 and m=+/-1 modes. The amplitude of the scattered wave approximately scales with the magnetic flux. A small amount of power is scattered into the p_1-mode. We have evaluated the absorption and phase shift from a Fourier-Hankel decomposition of the photospheric vertical velocities. They are compared with the results obtained from the emsemble average of 3400 small magnetic elements observed in high-resolution MDI Doppler datacubes. The comparison shows that the observed dependence of the phase shift with wavenumber can be matched reasonably well with the simulated flux tube model. The observed variation of the phase-shifts with the azimuthal order $m$ appears to depend on details of the ensemble averaging, including possible motions of the magnetic elements and asymmetrically shaped elements.Comment: Accepted for publication in The Astrophysical Journa

Sublinear Estimation of Weighted Matchings in Dynamic Data Streams

This paper presents an algorithm for estimating the weight of a maximum weighted matching by augmenting any estimation routine for the size of an unweighted matching. The algorithm is implementable in any streaming model including dynamic graph streams. We also give the first constant estimation for the maximum matching size in a dynamic graph stream for planar graphs (or any graph with bounded arboricity) using $\tilde{O}(n^{4/5})$ space which also extends to weighted matching. Using previous results by Kapralov, Khanna, and Sudan (2014) we obtain a $\mathrm{polylog}(n)$ approximation for general graphs using $\mathrm{polylog}(n)$ space in random order streams, respectively. In addition, we give a space lower bound of $\Omega(n^{1-\varepsilon})$ for any randomized algorithm estimating the size of a maximum matching up to a $1+O(\varepsilon)$ factor for adversarial streams

Internal Flows in Free Drops (IFFD)

Within the framework of an Earth-based research task investigating the internal flows within freely levitated drops, a low-gravity technology development experiment has been designed and carried out within the NASA Glovebox facility during the STS-83 and STS-94 Shuttle flights (MSL-1 mission). The goal was narrowly defined as the assessment of the capabilities of a resonant single-axis ultrasonic levitator to stably position free drops in the Shuttle environment with a precision required for the detailed measurement of internal flows. The results of this entirely crew-operated investigation indicate that the approach is fundamentally sound, but also that the ultimate stability of the positioning is highly dependent on the residual acceleration characteristic of the Spacecraft, and to a certain extent, on the initial drop deployment of the drop. The principal results are: the measured dependence of the residual drop rotation and equilibrium drop shape on the ultrasonic power level, the experimental evaluation of the typical drop translational stability in a realistic low-gravity environment, and the semi-quantitative evaluation of background internal flows within quasi-isothermal drops. Based on these results, we conclude that the successful design of a full-scale Microgravity experiment is possible, and would allow accurate the measurement of thermocapillary flows within transparent drops. The need has been demonstrated, however, for the capability for accurately deploying the drop, for a quiescent environment, and for precise mechanical adjustments of the levitator