Lagrangian Statistics of Navier-Stokes- and MHD-Turbulence

We report on a comparison of high-resolution numerical simulations of Lagrangian particles advected by incompressible turbulent hydro- and magnetohydrodynamic (MHD) flows. Numerical simulations were performed with up to $1024^3$ collocation points and 10 million particles in the Navier-Stokes case and $512^3$ collocation points and 1 million particles in the MHD case. In the hydrodynamics case our findings compare with recent experiments from Mordant et al. [1] and Xu et al. [2]. They differ from the simulations of Biferale et al. [3] due to differences of the ranges choosen for evaluating the structure functions. In Navier-Stokes turbulence intermittency is stronger than predicted by a multifractal approach of [3] whereas in MHD turbulence the predictions from the multifractal approach are more intermittent than observed in our simulations. In addition, our simulations reveal that Lagrangian Navier-Stokes turbulence is more intermittent than MHD turbulence, whereas the situation is reversed in the Eulerian case. Those findings can not consistently be described by the multifractal modeling. The crucial point is that the geometry of the dissipative structures have different implications for Lagrangian and Eulerian intermittency. Application of the multifractal approach for the modeling of the acceleration PDFs works well for the Navier-Stokes case but in the MHD case just the tails are well described.Comment: to appear in J. Plasma Phy

Effect of turbulence on collisions of dust particles with planetesimals in protoplanetary disks

Planetesimals in gaseous protoplanetary disks may grow by collecting dust particles. Hydrodynamical studies show that small particles generally avoid collisions with the planetesimals because they are entrained by the flow around them. This occurs when $St$, the Stokes number, defined as the ratio of the dust stopping time to the planetesimal crossing time, becomes much smaller than unity. However, these studies have been limited to the laminar case, whereas these disks are believed to be turbulent. We want to estimate the influence of gas turbulence on the dust-planetesimal collision rate and on the impact speeds. We used three-dimensional direct numerical simulations of a fixed sphere (planetesimal) facing a laminar and turbulent flow seeded with small inertial particles (dust) subject to a Stokes drag. A no-slip boundary condition on the planetesimal surface is modeled via a penalty method. We find that turbulence can significantly increase the collision rate of dust particles with planetesimals. For a high turbulence case (when the amplitude of turbulent fluctuations is similar to the headwind velocity), we find that the collision probability remains equal to the geometrical rate or even higher for $St\geq 0.1$, i.e., for dust sizes an order of magnitude smaller than in the laminar case. We derive expressions to calculate impact probabilities as a function of dust and planetesimal size and turbulent intensity

Integration of new experiments into the reflood map

The reflood map, developed in the last decade, was designed as a tool to focus on preventive and mitigative measures to prevent severe core damage scenarios and to identify research areas. In the meanwhile, additional experiments were performed and their data are available, so that an update seems necessary. For homogeneous particulate debris configurations, several experimental programs shed light on coolability. For the transition from in-core early phase to in-core late phase i.e. the loss of the rod like geometry, experimental database is still rather scarce due to the difficulty to describe the state of the core components and the fluid adequately. According to the new insight, an update of the database and the reflood map is discussed, also with respect to assess available grace times

Automatic Neuron Detection in Calcium Imaging Data Using Convolutional Networks

Calcium imaging is an important technique for monitoring the activity of thousands of neurons simultaneously. As calcium imaging datasets grow in size, automated detection of individual neurons is becoming important. Here we apply a supervised learning approach to this problem and show that convolutional networks can achieve near-human accuracy and superhuman speed. Accuracy is superior to the popular PCA/ICA method based on precision and recall relative to ground truth annotation by a human expert. These results suggest that convolutional networks are an efficient and flexible tool for the analysis of large-scale calcium imaging data.Comment: 9 pages, 5 figures, 2 ancillary files; minor changes for camera-ready version. appears in Advances in Neural Information Processing Systems 29 (NIPS 2016

Dispersion Relations for Thermally Excited Waves in Plasma Crystals

Thermally excited waves in a Plasma crystal were numerically simulated using a Box_Tree code. The code is a Barnes_Hut tree code proven effective in modeling systems composed of large numbers of particles. Interaction between individual particles was assumed to conform to a Yukawa potential. Particle charge, mass, density, Debye length and output data intervals are all adjustable parameters in the code. Employing a Fourier transform on the output data, dispersion relations for both longitudinal and transverse wave modes were determined. These were compared with the dispersion relations obtained from experiment as well as a theory based on a harmonic approximation to the potential. They were found to agree over a range of 0.9<k<5, where k is the shielding parameter, defined by the ratio between interparticle distance a and dust Debye length lD. This is an improvement over experimental data as current experiments can only verify the theory up to k = 1.5.Comment: 8 pages, Presented at COSPAR '0

Soil Nitrogen Dynamics Following Harvesting and Conversion of Red Alder and Douglas-Fir Stands

Drastic reductions in NO3- leaching have been observed after harvesting of mature red alder (Alnus rubra Bong.) stands. Our objective was to examine whether these reduction were linked to changes in soil N dynamics. Adjacent alder and Douglas fir [Pseudotsuga menziessii (Mirbel.) Franco] stands on young glacial soils (Alderwood; a loamy-skeletal, mixed, mesic, ortstein Aquic Haplorthod) in western Washington were harvested and replanted with either alder or Douglas fir seedlings; reference plots were established in nearby undisturbed stands. Three years after site conversion, when NO3- leaching declined most drastically in the harvested alder plots, net N mineralization and net nitrification in the upper soil were determined seasonally by in situ incubation using the buried-bag technique. There were no significant increases in soil NH4-N during incubation, indicating that nitrification rates closely followed net N mineralization in all plots. In the N-poor soils of the harvested and uncut Douglas fir plots, nitrification was consistently \u3c2 \u3eµg N g–1 per 30 d. In the N-rich alder reference plot, nitrification was highest in summer (52 µg N g–1 per 30 d), lowest in winter (4 µg N g–1), and intermediate in spring (20 µg N g–1). Seasonal fluctuations in nitrification diminished in the harvested alder plots. Nitrification rates were significantly lower than in the uncut alder plot in summer only (5–9 µg N g–1 per 30 d), but this difference was insufficient to account for the lower NO3- leaching rates in the conversion plots. Soils in the harvested plots were drier and experienced greater temperature extremes, but these changes did not correlate well with the variability in N-mineralization rates between uncut and harvested alder plots. Nitrification was positively correlated to soil temperature in the uncut alder plot only. The establishment of red alder seedlings did not appear to have influenced the N dynamics in either N-poor or N-rich soils

Analytical pre-test support of boil-down test QUENCH-11

Analytische Unterstützung zur Vorbereitung des Ausdampf-Versuchs QUENCH-11 Im QUENCH-Vorhaben des Forschungszentrums Karlsruhe soll das Fluten eines teilweise zerstörten Kerns untersucht werden. Der zweite LACOMERA Versuch Q-L2 (QUENCH-11) beginnt mit einer Ausdampfphase des Bündels, bis der Wasserspiegel das untere Bündel¬ende erreicht hat. Ein derartiger Versuch wurde bislang noch nicht in der QUENCH-Anlage durchgeführt, so dass mit SCDAP/RELAP5 mod3.2.irs eine Machbarkeitsstudie erforderlich war. Die Ergebnisse zeigen, dass eine Zusatzheizung im unteren Plenum notwendig ist, um den Wasserstand und die Verdampfungsrate (Dampfmassenstrom in der Ausdampfphase) unabhängig von der angestrebten Maximaltemperatur im Bündel zu regeln. Für eine verläss¬liche Versuchsplanung sowie zur Erstellung der Energiebilanz muss die Zusatzheizung in¬nerhalb des unteren Plenums unterhalb der Wasseroberfläche installiert werden, damit die Heizleistung vollständig in das Wasser eingekoppelt wird. Um die Verdampfungsrate über längere Zeit aufrecht zu erhalten, muss zusätzlich Wasser in das untere Plenum eingespeist werden. Anhand dieser Rechnungen wird der Testablauf im Detail diskutiert. Eine entsprechende Studie zeigte die Durchführbarkeit eines solchen Ausdampftests und war die Grundlage für die oben erwähnten Änderungen in der Anlage und der Versuchs-Durchführung gegenüber früheren Tests. Eine Reihe von Vorversuche wurde durchgeführt, um die Brauchbarkeit der Änderungen an der Anlage und der geplanten Versuchsführung zu prüfen und um Daten für das thermohydraulische Verhalten der Anlage zu bekommen, an denen die Code-Modelle für die Voraus- und Nachrechnungen von QUENCH-11 getestet werden können. Im Anschluss an die Vorversuche wurden wie bei früheren QUENCH-Tests detaillierte Vorausrechnungen mit verschiedenen Codes zu Versuchsablauf und -steuerung durchgeführt. Drei Forschungs¬einrichtungen in der EU waren beteiligt. Die berechneten Ergebnisse reagieren empfindlich auf Änderungen der Versuchsparameter wie das anfängliche axiale Temperaturprofil und die eingespeiste elektrische Leistung, wie es auch für die untersuchten physikalischen Bedin¬gungen im Versuch erwartet werden kann