Source author record

Seth Davidovits

Seth Davidovits 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

4works
6topics
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

4 published item(s)

preprint2022arXiv

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal(rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multi-phase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(σ^{2 Γ}_{ρ/\langle ρ\rangle}-1)^{1/2} (1-σ_{ρ\langle ρ\rangle}^{-2})^{-1/2}\mathcal{M}^{-1}Γ^{-1/2}$ with the density dispersion $σ_{ρ/ \langle ρ\rangle}$, the turbulent Mach number $\mathcal{M}$, and the polytropic exponent $Γ$. Purely solenoidal and purely compressive driving correspond to $b \sim 1/3$ and $b \sim 1$, respectively. Using simulations in which a shock is driven into a multi-phase medium with structures of different sizes and $Γ< 1$, we find $b \sim 1$ for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

preprint2022arXiv

Turbulence Generation by Shock Interaction with a Highly Non-Uniform Medium

An initially planar shock wave propagating into a medium of non-uniform density will be perturbed, leading to the generation of post-shock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly-non-uniform density (order-unity normalized variance, $σ_ρ/\overlineρ \sim 1$) and strong shocks (shock Mach numbers $\overline{M}_s \gtrsim 10$). This leads to a highly disrupted shock and a turbulent post-shock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the post-shock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the pre-shock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The post-shock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

preprint2016arXiv

Compressing turbulence and sudden viscous dissipation with compression-dependent ionization state

Turbulent plasma flow, amplified by rapid 3D compression, can be suddenly dissipated under continuing compression. This effect relies on the sensitivity of the plasma viscosity to the temperature, $μ\sim T^{5/2}$. The plasma viscosity is also sensitive to the plasma ionization state. We show that the sudden dissipation phenomenon may be prevented when the plasma ionization state increases during compression, and demonstrate the regime of net viscosity dependence on compression where sudden dissipation is guaranteed. Additionally, it is shown that, compared to cases with no ionization, ionization during compression is associated with larger increases in turbulent energy, and can make the difference between growing and decreasing turbulent energy.

preprint2009arXiv

An Efficient Algorithm for Classical Density Functional Theory in Three Dimensions: Ionic Solutions

Classical density functional theory (DFT) of fluids is a valuable tool to analyze inhomogeneous fluids. However, few numerical solution algorithms for three-dimensional systems exist. Here we present an efficient numerical scheme for fluids of charged, hard spheres that uses $\mathcal{O}(N\log N)$ operations and $\mathcal{O}(N)$ memory, where $N$ is the number of grid points. This system-size scaling is significant because of the very large $N$ required for three-dimensional systems. The algorithm uses fast Fourier transforms (FFT) to evaluate the convolutions of the DFT Euler-Lagrange equations and Picard (iterative substitution) iteration with line search to solve the equations. The pros and cons of this FFT/Picard technique are compared to those of alternative solution methods that use real-space integration of the convolutions instead of FFTs and Newton iteration instead of Picard. For the hard-sphere DFT we use Fundamental Measure Theory. For the electrostatic DFT we present two algorithms. One is for the \textquotedblleft bulk-fluid\textquotedblright functional of Rosenfeld [Y. Rosenfeld. \textit{J. Chem. Phys.} 98, 8126 (1993)] that uses $\mathcal{O}(N\log N)$ operations. The other is for the \textquotedblleft reference fluid density\textquotedblright (RFD) functional [D. Gillespie et al., J. Phys.: Condens. Matter 14, 12129 (2002)]. This functional is significantly more accurate than the bulk-fluid functional, but the RFD algorithm requires $\mathcal{O}(N^{2})$ operations.