Search NASA⌕ Search

Engineering topics

Bradshaw, S. J.

Publications and source records attributed to Bradshaw, S. J..

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗

Laser-induced-fluorescence imaging of a spin-polarized ultracold neutral plasma in a magnetic field

In this work, we report Doppler-sensitive laser-induced-fluorescence (LIF) imaging of an ultracold neutral plasma in a magnetic field. Local values of ion density, hydrodynamic fluid velocity, temperature, and spin polarization are obtained using a fluorescence model based on velocity-resolved rate equations (REs) including the transfer of ions between states due to laser coupling and spontaneous emission. The RE approach captures optical pumping of ions into states that are not driven by the LIF excitation laser, and this is validated with experimental data. Combined molecular-dynamics and quantum-trajectories simulations verify that velocity-changing collisions have a negligible impact on the state population evolution for typical experimental conditions. Relative intensities of Zeeman components of the LIF spectra provide clear evidence that the ions are electron-spin-polarized when created by photoionization of magnetically trapped 88 Sr atoms. This probe opens many possibilities for studying thermal transport and the equilibration of neutral plasmas in overlapping regimes of strong coupling and magnetization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding Heating in Active Region Cores through Machine Learning. I. Numerical Modeling and Predicted Observables

To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow-up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real observations of NOAA 1158 in terms of the underlying heating frequency.

UV radiation↗

Patterns of Activity in A Global Model of A Solar Active Region

In this work we investigate the global activity patterns predicted from a model active region heated by distributions of nanoflares that have a range of frequencies. What differs is the average frequency of the distributions. The activity patterns are manifested in time lag maps of narrow-band instrument channel pairs. We combine hydrodynamic and forward modeling codes with a magnetic field extrapolation to create a model active region and apply the time lag method to synthetic observations. Our aim is not to reproduce a particular set of observations in detail, but to recover some typical properties and patterns observed in active regions. Our key findings are the following. (1) Cooling dominates the time lag signature and the time lags between the channel pairs are generally consistent with observed values. (2) Shorter coronal loops in the core cool more quickly than longer loops at the periphery. (3) All channel pairs show zero time lag when the line of sight passes through coronal loop footpoints. (4) There is strong evidence that plasma must be re-energized on a timescale comparable to the cooling timescale to reproduce the observed coronal activity, but it is likely that a relatively broad spectrum of heating frequencies are operating across active regions. (5) Due to their highly dynamic nature, we find nanoflare trains produce zero time lags along entire flux tubes in our model active region that are seen between the same channel pairs in observed active regions.

Bradshaw, S. J.↗

X-Ray Source Heights in a Solar Flare: Thick-Target Versus Thermal Conduction Front Heating

Observations of solar flares with RHESSI have shown X-ray sources traveling along flaring loops, from the corona down to the chromosphere and back up. The 2002 November 28 C1.1 flare, first observed with RHESSI by Sui et al. and quantitatively analyzed by O'Flannagain et al., very clearly shows this behavior. By employing numerical experiments, we use these observations of X-ray source height motions as a constraint to distinguish between heating due to a non-thermal electron beam and in situ energy deposition in the corona. We find that both heating scenarios can reproduce the observed light curves, but our results favor non-thermal heating. In situ heating is inconsistent with the observed X-ray source morphology and always gives a height dispersion with photon energy opposite to what is observed.

Reep, J. W.↗

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗