Source author record

S. D. Baalrud

S. D. Baalrud appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

6 published item(s)

preprint2022arXiv

Strong Coulomb Coupling Influences Ion and Neutral Temperatures in Atmospheric Pressure Plasmas

Molecular dynamics simulations are used to model ion and neutral temperature evolution in partially-ionized atmospheric pressure plasma at different ionization fractions. Results show that ion-ion interactions are strongly coupled at ionization fractions as low as 10^-5 and that the temperature evolution is influenced by effects associated with the strong coupling. Specifically, disorder-induced heating is found to rapidly heat ions on a timescale of the ion plasma period (~10s ps) after an ionization pulse. This is followed by the collisional relaxation of ions and neutrals, which cools ions and heats neutrals on a longer (~ns) timescale. Slight heating then occurs over a much longer (~ 100s ns) timescale due to ion-neutral three-body recombination. An analytic model of the temperature evolution is developed that agrees with the simulation results. A conclusion is that strong coupling effects are important in atmospheric pressure plasmas.

preprint2020arXiv

Diffusion coefficients in the envelopes of white dwarfs

The diffusion of elements is a key process in understanding the unusual surface composition of white dwarfs stars and their spectral evolution. The diffusion coefficients of Paquette et al. (1986) have been widely used to model diffusion in white dwarfs. We perform new calculations of the coefficients of inter-diffusion and ionic thermal diffusion with 1) a more advanced model that uses a recent modification of the calculation of the collision integrals that is more suitable for the partially ionized, partially degenerate and moderately coupled plasma, and 2) classical molecular dynamics. The coefficients are evaluated for silicon and calcium in white dwarf envelopes of hydrogen and helium. A comparison of our results with Paquette et al. shows that the latter systematically underestimates the coefficient of inter-diffusion yet provides reliable estimates for the relatively weakly coupled plasmas found in nearly all types of stars as well as in white dwarfs with hydrogen envelopes. In white dwarfs with cool helium envelopes (Teff < 15000K), the difference grows to more than a factor of two. We also explored the effect of the ionization model used to determine the charges of the ions and found that it can be a substantial source of discrepancy between different calculations. Finally, we consider the relative diffusion time scales of Si and Ca in the context of the pollution of white dwarf photospheres by accreted planetesimals and find factor of > 3 differences between calculations based on Paquette et al. and our model.

preprint2013arXiv

Kinetic Theory of the Presheath and the Bohm Criterion

A kinetic theory of the Bohm criterion is developed that is based on positive-exponent velocity moments of the plasma kinetic equation. This result is contrasted with the conventional kinetic Bohm criterion that is based on a v^{-1} moment of the Vlasov equation. The salient difference between the two results is that low velocity particles dominate in the conventional theory, but are essentially unimportant in the new theory. It is shown that the derivation of the conventional kinetic Bohm criterion is flawed. Low velocity particles can cause unphysical divergences in the conventional theory. These divergent contributions are avoided with this new approach. The two theories are compared using example distribution functions from previous presheath models. The importance of ion-ion and electron-electron collisions to determining the particle distribution functions throughout the presheath is also discussed. A kinetic equation that accounts for wave-particle scattering by convective instabilities is used to show that ion-acoustic instabilities in the presheath of low temperature plasmas (where T_e \gg T_i) can cause both ions and electrons to obtain Maxwellian distribution functions near the sheath.

preprint2013arXiv

Particle-in-cell simulations of a current-free double layer

Current-free double layers of the type reported in plasmas in the presence of an expanding magnetic field [C. Charles and R. W. Boswell, Appl. Phys. Lett. 82, 1356 (2003)] are modeled theoretically and with particle-in-cell/Monte Carlo simulations. Emphasis is placed on determining what mechanisms affect the electron velocity distribution function (EVDF) and how the EVDF influences the double layer. A theoretical model is developed based on depletion of electrons in certain velocity intervals due to wall losses and repletion of these intervals due to ionization and elastic electron scattering. This model is used to predict the range of neutral pressures over which a double layer can form and the electrostatic potential drop of the double layer. These predictions are shown to compare well with simulation results.

preprint2011arXiv

Hall magnetohydrodynamic reconnection in the plasmoid unstable regime

A set of reduced Hall magnetohydrodynamic (MHD) equations are used to evaluate the stability of large aspect ratio current sheets to the formation of plasmoids (secondary islands). Reconnection is driven by resistivity in this analysis, which occurs at the resistive skin depth $d_η\equiv S_L^{-1/2} \sqrt{L v_A/γ}$, where $S_L$ is the Lundquist number, $L$ the length of the current sheet, $v_A$ the Alfvén speed, and $γ$ the growth rate. Modifications to a recent resistive MHD analysis [N.\ F.\ Loureiro, A.\ A.\ Schekochihin, and S.\ C. Cowley, Phys.\ Plasmas {\bf 14}, 100703 (2007)] arise when collisions are sufficiently weak that $d_η$ is shorter than the ion skin depth $d_i \equiv c/ω_{pi}$. Secondary islands grow faster in this Hall MHD regime: the maximum growth rate scales as $(d_i/L)^{6/13} S_L^{7/13} v_A/L$ and the number of plasmoids as $(d_i/L)^{1/13} S_L^{11/26}$, compared to $S_L^{1/4} v_A/L$ and $S^{3/8}$, respectively, in resistive MHD.

preprint2011arXiv

Reduced magnetohydrodynamic theory of oblique plasmoid instabilities

The three-dimensional nature of plasmoid instabilities is studied using the reduced magnetohydrodynamic equations. For a Harris equilibrium with guide field, represented by $\vc{B}_o = B_{po} \tanh (x/λ) \hat{y} + B_{zo} \hat{z}$, a spectrum of modes are unstable at multiple resonant surfaces in the current sheet, rather than just the null surface of the polodial field $B_{yo} (x) = B_{po} \tanh (x/λ)$, which is the only resonant surface in 2D or in the absence of a guide field. Here $B_{po}$ is the asymptotic value of the equilibrium poloidal field, $B_{zo}$ is the constant equilibrium guide field, and $λ$ is the current sheet width. Plasmoids on each resonant surface have a unique angle of obliquity $θ\equiv \arctan(k_z/k_y)$. The resonant surface location for angle $θ$ is $x_s = - λ\arctanh (\tan θB_{zo}/B_{po})$, and the existence of a resonant surface requires $|θ| < \arctan (B_{po} / B_{zo})$. The most unstable angle is oblique, i.e. $θ\neq 0$ and $x_s \neq 0$, in the constant-$ψ$ regime, but parallel, i.e. $θ= 0$ and $x_s = 0$, in the nonconstant-$ψ$ regime. For a fixed angle of obliquity, the most unstable wavenumber lies at the intersection of the constant-$ψ$ and nonconstant-$ψ$ regimes. The growth rate of this mode is $γ_{\textrm{max}}/Γ_o \simeq S_L^{1/4} (1-μ^4)^{1/2}$, in which $Γ_o = V_A/L$, $V_A$ is the Alfvén speed, $L$ is the current sheet length, and $S_L$ is the Lundquist number. The number of plasmoids scales as $N \sim S_L^{3/8} (1-μ^2)^{-1/4} (1 + μ^2)^{3/4}$.