We study an individual based model describing competition in space between two different alleles. Although the model is similar in spirit to classic models of spatial population genetics such as the ...stepping stone model, here however space is continuous and the total density of competing individuals fluctuates due to demographic stochasticity. By means of analytics and numerical simulations, we study the behavior of fixation probabilities, fixation times, and heterozygosity, in a neutral setting and in cases where the two species can compete or cooperate. By concluding with examples in which individuals are transported by fluid flows, we argue that this model is a natural choice to describe competition in marine environments.
We present a collection of eight data sets from state-of-the-art experiments and numerical simulations on turbulent velocity statistics along particle trajectories obtained in different flows with ...Reynolds numbers in the range R{lambda}in120:740. Lagrangian structure functions from all data sets are found to collapse onto each other on a wide range of time lags, pointing towards the existence of a universal behavior, within present statistical convergence, and calling for a unified theoretical description. Parisi-Frisch multifractal theory, suitably extended to the dissipative scales and to the Lagrangian domain, is found to capture the intermittency of velocity statistics over the whole three decades of temporal scales investigated here.
We compare experimental data and numerical simulations for the dynamics of inertial particles with finite density in turbulence. In the experiment, bubbles and solid particles are optically tracked ...in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The probability density functions (PDF) of particle accelerations and their auto-correlation in time are computed. Numerical results are obtained from a direct numerical simulation in which a suspension of passive pointwise particles is tracked, with the same finite density and the same response time as in the experiment. We observe a good agreement for both the variance of acceleration and the autocorrelation time scale of the dynamics; small discrepancies on the shape of the acceleration PDF are observed. We discuss the effects induced by the finite size of the particles, not taken into account in the present numerical simulations.
The achievement of Bose–Einstein condensation in ultra-cold vapours of alkali atoms has given enormous impulse to the study of dilute atomic gases in condensed quantum states inside magnetic traps ...and optical lattices. High-purity and easy optical access make them ideal candidates to investigate fundamental issues on interacting quantum systems. This review presents some theoretical issues which have been addressed in this area and the numerical techniques which have been developed and used to describe them, from mean-field models to classical and quantum simulations for equilibrium and dynamical properties. After an introductory overview on dilute quantum gases, both in the homogeneus state and under harmonic or periodic confinement, the article is organized in three main sections. The first concerns Bose-condensed gases at zero temperature, with main regard to the properties of the ground state in different confinements and to collective excitations and transport in the condensate. Bose–Einstein-condensed gases at finite temperature are addressed in the next section, the main emphasis being on equilibrium properties and phase transitions and on dynamical and transport properties associated with the presence of the thermal cloud. Finally, the last section is focused on theoretical and computational issues that have emerged from the efforts to drive gases of fermionic atoms and boson–fermion mixtures deep into the quantum degeneracy regime, with the aim of realizing novel superfluids from fermion pairing. The attention given in this article to methods beyond standard mean-field approaches should make it a useful reference point for future advances in these areas.
We investigate and compare the accuracy and efficiency of different numerical approaches to model the dynamics of finite-size particles using the lattice Boltzmann method (LBM). This includes the ...standard bounce-back (BB) and the equilibrium interpolation (EI) schemes. To accurately compare the different implementations, we first introduce a boundary condition to approximate the flow properties of an unbounded fluid in a finite simulation domain, taking into account the perturbation induced by a moving particle. We show that this boundary treatment is efficient in suppressing detrimental effects on the dynamics of spherical and ellipsoidal particles arising from the finite size of the simulation domain. We then investigate the performances of the BB and EI schemes in modeling the dynamics of a spherical particle settling under Stokes conditions, which can now be reproduced with great accuracy thanks to the treatment of the exterior boundary. We find that the EI scheme outperforms the BB scheme in providing a better accuracy scaling with respect to the resolution of the settling particle, while suppressing finite-size effects due to the particle discretization on the lattice grid. Additionally, in order to further increase the capability of the algorithm in modeling particles of sizes comparable to the lattice spacing, we propose an improvement to the EI scheme, the complete equilibrium interpolation (CEI). This approach allows us to accurately capture the boundaries of the particle also when located between two fluid nodes. We evaluate the CEI performance in solving the dynamics of an under-resolved particle under analogous Stokes conditions and also for the case of a rotating ellipsoid in a shear flow. Finally, we show that EI and CEI are able to recover the correct flow solutions also at small, but finite, Reynolds number. Adopting the CEI scheme it is not only possible to detect particles with zero lattice occupation, but also to increase up to one order of magnitude the accuracy of the dynamics of particles with a size comparable to the lattice spacing with respect to the BB and the EI schemes.
An important aspect in numerical simulations of particle-laden turbulent flows is the interpolation of the flow field needed for the computation of the Lagrangian trajectories. The accuracy of the ...interpolation method has direct consequences for the acceleration spectrum of the fluid particles and is therefore also important for the correct evaluation of the hydrodynamic forces for almost neutrally buoyant particles, common in many environmental applications. In order to systematically choose the optimal tradeoff between interpolation accuracy and computational cost we focus on comparing errors: the interpolation error is compared with the discretization error of the flow field. In this way one can prevent unnecessary computations and still retain the accuracy of the turbulent flow simulation. From the analysis a practical method is proposed that enables direct estimation of the interpolation and discretization error from the energy spectrum. The theory is validated by means of direct numerical simulations (DNS) of homogeneous, isotropic turbulence using a spectral code, where the trajectories of fluid tracers are computed using several interpolation methods. We show that B-spline interpolation has the best accuracy given the computational cost. Finally, the optimal interpolation order for the different methods is shown as a function of the resolution of the DNS simulation.