Source author record

Giorgio Krstulovic

Giorgio Krstulovic 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

28works
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

28 published item(s)

preprint2022arXiv

Critical velocity for vortex nucleation and roton emission in a generalized model for superfluids

We study numerically the process of vortex nucleation at the wake of a moving object in superfluids using a generalized and non-local Gross-Pitaevskii model. The non-local potential is set to reproduce the roton minimum present in the excitation spectrum of superfluid helium. By applying numerically a Newton-Raphson method we determine the bifurcation diagram for different types of non-linearities and object sizes which allow for determining the corresponding critical velocities. In the case of a non-local potential, we observe that for small object sizes the critical velocity is simply determined by the Landau criterion for superfluidity whereas for large objects there is little difference between all models studied. Finally, we study dynamically in two and three dimensions how rotons and vortices are excited in the non-local model of superfluid.

preprint2022arXiv

Energy spectrum of two-dimensional acoustic turbulence

We report an exact unique constant-flux power-law analytical solution of the wave kinetic equation for the turbulent energy spectrum, $E(k)=C_1 \sqrt{\varepsilon\, a c_{\rm s} }/k$, of acoustic waves in 2D with almost linear dispersion law, $ω_k = c_{\rm s} k[1+(ak)^2]$, $ ak \ll 1$. Here $\varepsilon$ is the energy flux over scales, and $C_1$ is the universal constant which was found analytically. Our theory describes, for example, acoustic turbulence in 2D Bose-Einstein condensates (BECs). The corresponding 3D counterpart of turbulent acoustic spectrum was found over half a century ago, however, due to the singularity in 2D, no solution has been obtained until now. We show the spectrum $E(k)$ is realizable in direct numerical simulations of forced-dissipated Gross-Pitaevskii equation in the presence of strong condensate.

preprint2022arXiv

Testing wave turbulence theory for Gross-Pitaevskii system

We test the predictions of the theory of weak wave turbulence by performing numerical simulations of the Gross-Pitaevskii equation (GPE) and the associated wave-kinetic equation (WKE). We consider an initial state localized in Fourier space, and we confront the solutions of the WKE obtained numerically with GPE data for both the wave-action spectrum and the probability density functions (PDFs) of the Fourier mode intensities. We find that the temporal evolution of the GPE data is accurately predicted by the WKE, with no adjustable parameters, for about two nonlinear kinetic times. Qualitative agreement between the GPE and the WKE persists also for longer times with some quantitative deviations that may be attributed to the combination of breakdown of the theoretical assumptions underlying the WKE as well as numerical issues. Furthermore, we study how the wave statistics evolves toward Gaussianity in a time scale of the order of the kinetic time.The excellent agreement between direct numerical simulations of the GPE and the WKE provides a new and solid ground to the theory of weak wave turbulence.

preprint2022arXiv

Velocity circulation intermittency in finite-temperature turbulent superfluid helium

We study intermittency of circulation moments in turbulent superfluid helium by using experimental grid turbulence and numerical simulations of the Hall-Vinen-Bekarevich-Khalatnikov model. More precisely, we compute the velocity circulation $Γ_r$ in loops of size $r$ laying in the inertial range. For both, experimental and numerical data, the circulation variance shows a clear Kolmogorov scaling $\langle Γ_r^2 \rangle \sim r^{8/3}$ in the inertial range, independently of the temperature. Scaling exponents of high-order moments are comparable, within error bars, to previously reported anomalous circulation exponents in classical turbulence and low-temperature quantum turbulence numerical simulations.

preprint2021arXiv

Comment on "Theoretical analysis of quantum turbulence using the Onsager ideal turbulence theory''

In a recent paper [T. Tanogami Phys. Rev. E 103, 023106 ] proposes a scenario for quantum turbulence where the energy spectrum at scales smaller than the inter-vortex distance is dominated by a quantum stress cascade, in opposition to Kelvin wave cascade predictions. The purpose of the present comment is to highlight some physical issues in the derivation of the quantum stress cascade, in particular to stress that quantization of circulation has been ignored.

preprint2020arXiv

A new self-consistent approach of quantum turbulence in superfluid helium

We present the Fully cOUpled loCAl model of sUperfLuid Turbulence (FOUCAULT) that describes the dynamics of finite temperature superfluids. The superfluid component is described by the vortex filament method while the normal fluid is governed by a modified Navier-Stokes equation. The superfluid vortex lines and normal fluid components are fully coupled in a self-consistent manner by the friction force, which induces local disturbances in the normal fluid in the vicinity of vortex lines. The main focus of this work is the numerical scheme for distributing the friction force to the mesh points where the normal fluid is defined and for evaluating the velocity of the normalfluid on the Lagrangian discretization points along the vortex lines. In particular, we show that if this numerical scheme is not careful enough, spurious results may occur. The new scheme which we propose to overcome these difficulties is based on physical principles. Finally, we apply the new method to the problem of the motion of a superfluid vortex ring in a stationary normal fluid and in a turbulent normal fluid.

preprint2020arXiv

Active and finite-size particles in decaying quantum turbulence at low temperature

The evolution of a turbulent tangle of quantum vortices in presence of finite-size active particles is studied by means of numerical simulations of the Gross-Pitaevskii equation. Particles are modeled as potentials depleting the superfluid and described with classical degrees of freedom following a Newtonian dynamics. It is shown that particles do not modify the building-up and the decay of the superfluid Kolmogorov turbulent regime. It is observed that almost the totality of particles remains trapped inside quantum vortices, although they are occasionally detached and recaptured. The statistics of this process are presented and discussed. The particle Lagrangian dynamics is also studied. At large time scales, the velocity spectrum of particles is reminiscent of a classical Lagrangian turbulent behavior. At time-scales faster than the turnover time associated to the mean inter-vortex distance, the particle motion is dominated by oscillations due to Magnus effect. For light particles a non-classical scaling of the spectrum arises. The particle velocity and acceleration probability distribution functions are then studied. The decorrelation time of the particle acceleration is found to be shorter than in classical fluids, and related to the Magnus force experienced by the trapped particles.

preprint2020arXiv

How trapped particles interact with and sample superfluid vortex excitations

Particles have been used for more than a decade to visualize and study the dynamics of quantum vortices in superfluid helium. In this work we study how the dynamics of a collection of particles set inside a vortex reflects the motion of the vortex. We use a self-consistent model based on the Gross-Pitaevskii equation coupled with classical particle dynamics. We find that each particle oscillates with a natural frequency proportional to the number of vortices attached to it. We then study the dynamics of an array of particles trapped in a quantum vortex and use particle trajectories to measure the frequency spectrum of the vortex excitations. Surprisingly, due to the discreetness of the array, the vortex excitations measured by the particles exhibits bands, gaps and Brillouin zones, analogous to the ones of electrons moving in crystals. We then establish a mathematical analogy where vortex excitations play the role of electrons and particles that of the potential barriers constituting the crystal. We find that the height of the effective potential barriers is proportional to the particle mass and the frequency of the incoming waves. We conclude that large-scale vortex excitations could be in principle directly measured by particles and novel physics could emerge from particle-vortex interaction.

preprint2020arXiv

Inhomogeneous distribution of particles in co-flow and counterflow quantum turbulence

Particles are today the main tool to study superfluid turbulence and visualize quantum vortices. In this work, we study the dynamics and the spatial distribution of particles in co-flow and counterflow superfluid helium turbulence in the framework of the two-fluid Hall-Vinen-Bekarevich-Khalatnikov (HVBK) model. We perform three-dimensional numerical simulations of the HVBK equations along with the particle dynamics that depends on the motion of both fluid components. We find that, at low temperatures, where the superfluid mass fraction dominates, particles strongly cluster in vortex filaments regardless of their physical properties. At higher temperatures, as viscous drag becomes important and the two components become tightly coupled, the clustering dynamics in the coflowing case approach those found in classical turbulence, while under strong counterflow, the particle distribution is dominated by the quasi-two-dimensionalization of the flow.

preprint2020arXiv

Quantum vortex reconnections mediated by trapped particles

Reconnections between quantum vortex filaments in presence of trapped particles are investigated using numerical simulations of the Gross--Pitaevskii equation. Particles are described with classical degrees of freedom and modeled as highly repulsive potentials which deplete the superfluid. First, the case of a vortex dipole with a single particle trapped inside one of the vortices is studied. It is shown that the reconnection takes place at the position of the particle as a consequence of the symmetry breaking induced by it. The separation rate between the reconnecting points is compatible with the known dynamics of quantum vortex reconnections and it is independent of the particle mass and size. After the reconnection, the particle is pushed away with a constant velocity and its trajectory is deflected because of the transverse momentum exchange with the vortex filaments. The momentum exchanges between the particle, the vortex, and a density pulse are characterized. Finally, the reconnection of two linked rings, each of them with several initially randomly distributed particles is studied. It is observed that generically, reconnections take place at the location of trapped particles. It is shown that reconnection dynamics is unaffected for light particles.

preprint2020arXiv

Stochastic motion of finite-size immiscible impurities in a dilute quantum fluid at finite temperature

The dynamics of an active, finite-size and immiscible impurity in a dilute quantum fluid at finite temperature is characterized by means of numerical simulations of the projected Gross--Pitaevskii equation. The impurity is modeled as a localized repulsive potential and described with classical degrees of freedom. It is shown that impurities of different sizes thermalize with the fluid and undergo a stochastic dynamics compatible with an Ornstein--Uhlenbeck process at sufficiently large time-lags. The velocity correlation function and the displacement of the impurity are measured and an increment of the friction with temperature is observed. Such behavior is phenomenologically explained in a scenario where the impurity exchanges momentum with a dilute gas of thermal excitations, experiencing an Epstein drag.

preprint2016arXiv

A Vortex Filament Tracking Method for the Gross-Pitaevskii Model of a Superfluid

We present an accurate and robust numerical method to track quantized vortex lines in a superfluid described by the Gross-Pitaevskii equation. By utilizing the pseudo-vorticity field of the associated complex scalar order parameter of the superfluid, we are able to track the topological defects of the superfluid and reconstruct the vortex lines which correspond to zeros of the field. Throughout, we assume our field is periodic to allow us to make extensive use of the Fourier representation of the field and its derivatives in order to retain spectral accuracy. We present several case studies to test the precision of the method which include the evaluation of the curvature and torsion of a torus vortex knot, and the measurement of the Kelvin wave spectrum of a vortex line and a vortex ring. The method we present makes no a-priori assumptions on the geometry of the vortices and is therefore applicable to a wide range of systems such as a superfluid in a turbulent state that is characterised by many vortex rings coexisting with sound waves. This allows us to track the positions of the vortex filaments in a dense turbulent vortex tangle and extract statistical information about the distribution of the size of the vortex rings and the inter-vortex separations. In principle, the method can be extended to track similar topological defects arising in other physical systems.

preprint2016arXiv

Evolution of a superfluid vortex filament tangle driven by the Gross-Pitaevskii equation

The development and decay of a turbulent vortex tangle driven by the Gross-Pitaevskii equation is studied. Using a recently-developed accurate and robust tracking algorithm, all quantised vortices are extracted from the fields. The Vinen's decay law for the total vortex length with a coefficient that is in quantitative agreement with the values measured in Helium II is observed. The topology of the tangle is then studied showing that linked rings may appear during the decay. The tracking also allows for determining the statistics of small-scales quantities of vortex lines, exhibiting large fluctuations of curvature and torsion. Finally, the temporal evolution of the Kelvin wave spectrum is obtained providing evidence of the development of a weak-wave turbulence cascade.

preprint2016arXiv

Grid superfluid turbulence and intermittency at very low temperature

Low-temperature grid generated turbulence is investigated by using numerical simulations of the Gross-Pitaevskii equation. The statistics of regularized velocity increments are studied. Increments of the incompressible velocity are found to be skewed for turbulent states. Results are later confronted with the (quasi) homogeneous and isotropic Taylor-Green flow, revealing the universality of the statistics. For this flow, the statistics are found to be intermittent and a Kolmogorov constant close to the one of classical fluid is found for the second order structure function.

preprint2016arXiv

Statistical steady state in turbulent droplet condensation

Motivated by systems in which droplets grow and shrink in a turbulence-driven supersaturation field, we investigate the problem of turbulent condensation in a general manner. Using direct numerical simulations we show that the turbulent fluctuations of the supersaturation field offer different conditions for the growth of droplets which evolve in time due to turbulent transport and mixing. Based on that, we propose a Lagrangian stochastic model for condensation and evaporation of small droplets in turbulent flows. It consists of a set of stochastic integro-differential equations for the joint evolution of the squared radius and the supersaturation along the droplet trajectories. The model has two parameters fixed by the total amount of water and the thermodynamic properties, as well as the Lagrangian integral timescale of the turbulent supersaturation. The model reproduces very well the droplet size distributions obtained from direct numerical simulations and their time evolution. A noticeable result is that, after a stage where the squared radius simply diffuses, the system converges exponentially fast to a statistical steady state independent of the initial conditions. The main mechanism involved in this convergence is a loss of memory induced by a significant number of droplets undergoing a complete evaporation before growing again. The statistical steady state is characterised by an exponential tail in the droplet mass distribution. These results reconcile those of earlier numerical studies, once these various regimes are considered.

preprint2015arXiv

A lattice model for the Eulerian description of heavy particle suspensions

Modeling dispersed solid phases in fluids still represents a computational challenge when considering a small-scale coupling in wide systems, such as the atmosphere or industrial processes at high Reynolds numbers. A numerical method is here introduced for simulating the dynamics of diffusive heavy inertial particles in turbulent flows. The approach is based on the position/velocity phase-space particle distribution. The discretization of velocities is inspired from lattice Boltzmann methods and is chosen to match discrete displacements between two time steps. For each spatial position, the time evolution of particles momentum is approximated by a finite-volume approach. The proposed method is tested for particles experiencing a Stokes viscous drag with a prescribed fluid velocity field in one dimension using a random flow, and in two dimensions with the solution to the forced incompressible Navier-Stokes equations. Results show good agreements between Lagrangian and Eulerian dynamics for both spatial clustering and the dispersion in particle velocities. This demonstrates the suitability of the proposed approach at large Stokes numbers or for situations where details of collision processes are important.

preprint2014arXiv

Clustering, fronts, and heat transfer in turbulent suspensions of heavy particles

Heavy inertial particles transported by a turbulent flow are shown to concentrate in the regions where an advected passive scalar, such as temperature, displays very strong front-like discontinuities. This novel effect is responsible for extremely high levels of fluctuations for the passive field sampled by the particles that impacts the heat fluxes exchanged between the particles and the surrounding fluid. Instantaneous and averaged heat fluxes are shown to follow strongly intermittent statistics and anomalous scaling laws.

preprint2014arXiv

Structures and Lagrangian statistics of the Taylor-Green Dynamo

The evolution of a Taylor-Green forced magnetohydrodynamic (MHD) system showing dynamo activity is analyzed via direct numerical simulations. The statistical properties of the velocity and magnetic field in Eulerian coordinates and along trajectories of fluid elements (Lagrangian coordinates) are found to change between the kinematic, non-linear and saturated regime. Fluid element (tracer) trajectories change from chaotic quasi-isotropic (kinematic phase) to mean magnetic field aligned (saturated phase). The probability density functions (PDFs) of the magnetic field change from strongly non-Gaussian in the kinematic to quasi-Gaussian PDFs in the saturated regime so that their flatness give a precise handle on the definition of the limiting points of the three regimes. Also the statistics of the fluctuations of the kinetic and magnetic energy along fluid trajectories change. All this goes along with a dramatic increase of the correlation time of velocity and magnetic field fluctuations experienced by tracers significantly exceeding one turbulent large-eddy turn-over time. A remarkable consequence is an intermittent scaling regime of the Lagrangian magnetic field structure functions at unusually long time scales.

preprint2014arXiv

Turbulent pair dispersion as a continuous-time random walk

The phenomenology of turbulent relative dispersion is revisited. A heuristic scenario is proposed, in which pairs of tracers undergo a succession of independent ballistic separations during time intervals whose lengths fluctuate. This approach suggests that the logarithm of the distance between tracers self-averages and performs a continuous-time random walk. This leads to specific predictions for the probability distribution of separations, that differ from those obtained using scale-dependent eddy-diffusivity models (e.g. in the framework of Richardson's approach). Such predictions are tested against high-resolution simulations and shed new lights on the explosive separation between tracers.

preprint2012arXiv

Kelvin-wave cascade and dissipation in low-temperature superfluids vortices

We study the statistical properties of the Kelvin waves propagating along quantized superfluid vortices driven by the Gross-Pitaevskii equation. No artificial forcing or dissipation is added. Vortex positions are accurately tracked. This procedure directly allows us to obtain the Kevin-waves occupation-number spectrum. Numerical data obtained from long time integration and ensemble-average over initial conditions supports the spectrum proposed in [L'vov and Nazarenko, JETP Lett 91, 428 (2010)]. Kelvin wave modes in the inertial range are found to be Gaussian as expected by weak-turbulence predictions. Finally the dissipative range of the Kelvin-wave spectrum is studied. Strong non-Gaussian fluctuations are observed in this range.

preprint2012arXiv

Mass fluctuations and diffusion in time-dependent random environments

A mass ejection model in a time-dependent random environment with both temporal and spatial correlations is introduced. When the environment has a finite correlation length, individual particle trajectories are found to diffuse at large times with a displacement distribution that approaches a Gaussian. The collective dynamics of diffusing particles reaches a statistically stationary state, which is characterized in terms of a fluctuating mass density field. The probability distribution of density is studied numerically for both smooth and non-smooth scale-invariant random environments. A competition between trapping in the regions where the ejection rate of the environment vanishes and mixing due to its temporal dependence leads to large fluctuations of mass. These mechanisms are found to result in the presence of intermediate power-law tails in the probability distribution of the mass density. For spatially differentiable environments, the exponent of the right tail is shown to be universal and equal to -3/2. However, at small values, it is found to depend on the environment. Finally, spatial scaling properties of the mass distribution are investigated. The distribution of the coarse-grained density is shown to posses some rescaling properties that depend on the scale, the amplitude of the ejection rate, and the Hölder exponent of the environment.

preprint2011arXiv

Alfvén waves and ideal two-dimensional Galerkin truncated magnetohydrodynamics

We investigate numerically the dynamics of two-dimensional Euler and ideal magnetohydrodynamics (MHD) flows in systems with a finite number of modes, up to $4096^2$, for which several quadratic invariants are preserved by the truncation and the statistical equilibria are known. Initial conditions are the Orszag-Tang vortex with a neutral X-point centered on a stagnation point of the velocity field in the large scales. In MHD, we observe that the total energy spectra at intermediate times and intermediate scales correspond to the interactions of eddies and waves, $E_T(k)\sim k^{-3/2}$. Moreover, no dissipative range is visible neither for Euler nor for MHD in two dimensions; in the former case, this may be linked to the existence of a vanishing turbulent viscosity whereas in MHD, the numerical resolution employed may be insufficient. When imposing a uniform magnetic field to the flow, we observe a lack of saturation of the formation of small scales together with a significant slowing-down of their equilibration, with however a cut-off independent partial thermalization being reached at intermediate scales.

preprint2011arXiv

Anomalous vortex ring velocities induced by thermally-excited Kelvin waves and counterflow effects in superfluids

Dynamical counterflow effects on vortex evolution under the truncated Gross-Pitaevskii equation are investigated. Standard longitudinal mutual friction effects are produced and a dilatation of vortex rings is obtained at large counterflow. A strong temperature-dependent anomalous slowdown of vortex rings is observed and attributed to the presence of thermally exited Kelvin waves. This generic effect of finite-temperature superfluids is estimated using energy equipartition and orders of magnitude are given for weakly interacting Bose-Einstein condensates and superfluid $^4{\rm He}$.

preprint2011arXiv

Axial dipolar dynamo action in the Taylor-Green vortex

We present a numerical study of the magnetic field generated by the Taylor-Green vortex. We show that periodic boundary conditions can be used to mimic realistic boundary conditions by prescribing the symmetries of the velocity and magnetic fields. This gives insight in some problems of central interest for dynamos: the possible effect of velocity fluctuations on the dynamo threshold, the role of boundary conditions on the threshold and on the geometry of the magnetic field generated by dynamo action. In particular, we show that an axial dipolar dynamo similar to the one observed in a recent experiment can be obtained with an appropriate choice of the symmetries of the magnetic field. The nonlinear saturation is studied and a simple model explaining the magnetic Prandtl number dependence of the super/sub critical nature of the dynamo transition is given.

preprint2011arXiv

Dispersive bottleneck delaying thermalization of turbulent Bose-Einstein Condensates

A new mechanism of thermalization involving a direct energy cascade is obtained in the truncated Gross-Pitaevskii dynamics. A long transient with partial thermalization at small-scales is observed before the system reaches equilibrium. Vortices are found to disappear as a prelude to final thermalization. A bottleneck that produces spontaneous effective self-truncation and delays thermalization is characterized when large dispersive effects are present at the truncation wavenumber. Order of magnitude estimates indicate that self-truncation takes place in turbulent Bose-Einstein condensates. This effect should also be present in classical hydrodynamics and models of turbulence.

preprint2011arXiv

Energy cascade with small-scales thermalization, counterflow metastability and anomalous velocity of vortex rings in Fourier-truncated Gross-Pitaevskii equation

The statistical equilibria of the dynamics of the Gross-Pitaevskii Equation (GPE) with a finite range of spatial Fourier modes are characterized using a new algorithm, based on a stochastically forced Ginzburg-Landau equation (SGLE), that directly generates grand canonical distributions.. The SGLE-generated distributions are validated against finite-temperature GPE-thermalized states and exact low-temperature results. A standard second-order $λ$-transition is exhibited. A new mechanism of GPE thermalization through a direct cascade of energy is found using initial conditions with mass and energy distributed at large scales. A long transient with partial thermalization at small-scales is observed. Vortices are shown to disappear as a prelude to final thermalization and their annihilation is related to the contraction of vortex rings due to mutual friction. Increasing the amount of dispersion at truncation wavenumber is shown to slowdown thermalization and vortex annihilation. A bottleneck that produces spontaneous effective self truncation with partial thermalization is characterized in the limit of large dispersive effects. Metastable counter-flow states are generated using the SGLE algorithm. Spontaneous nucleation of vortex ring is observed and the corresponding Arrhenius law is characterized. Dynamical counter-flow effects on vortex evolution are investigated. Longitudinal effects are produced and measured. A dilatation of vortex rings is obtained for larger counterflows. The vortex ring longitudinal velocity has a strong dependence on temperature, an effect that is related to the presence of finite-amplitude Kelvin waves. This anomalous vortex ring velocity is quantitatively reproduced by assuming equipartition of energy of the Kelvin waves. Orders of magnitude are given for the predicted effects in weakly interacting Bose-Einstein condensates and superfluid $^4{\rm He}$.

preprint2009arXiv

Comment on "Superfluid turbulence from quantum Kelvin wave to classical Kolmogorov cascade". [arXiv:0905.0159v1]

In this comment we point out that the high wavenumber $k^{-3}$ power-law observed by the PRL, [v. 103, 084501 (2009) by J. Yepez, G. Vahala, L.Vahala and M. Soe, arXiv:0905.0159] is an artifact stemming from the definition of the kinetic energy spectra and is thus not directly related to a Kelvin wave cascade. We also clarify a confusion about the wavenumber intervals on which Kolmogorov and Kelvin wave cascades are expected to take place.