Search

The SARS-CoV-2 is transmitted not only through coughing, but also through breathing, speaking or singing. We perform direct numerical simulations of the turbulent transport of potentially infectious aerosols in short conversations, involving repetitive phrases separated by quiescent intervals. We estimate that buoyancy effects due to droplet evaporation are small, and neglect them. A two-way conversation is shown to significantly reduce the aerosol exposure compared with a relative monologue by one person and relative silence of the other. This is because of the 'cancelling' effect produced by the two interacting speech jets. Unequal conversation is shown to significantly increase the infection risk to the person who talks less. Interestingly, a small height difference is worse for infection spread, due to reduced interference between the speech jets, than two faces at the same level. For small axial separation, speech jets show large oscillations and reach the other person intermittently. We suggest a range of lateral separations between two people to minimize transmission risk. A realistic estimate of the infection probability is provided by including exposure through the eyes and mouth, in addition to the more common method of using inhaled virions alone. We expect that our results will provide useful inputs to epidemiological models and to disease management. [ABSTRACT FROM AUTHOR]

To make predictions for and design fusion reactors, a multitude of physical processes must be considered. In the edge and scrape-off layer (SOL), turbulent fluctuations intertwine with the plasma background, which is largely determined by neutral gas, and magnetic geometry plays an important role. A diffusive neutrals model has now been implemented in the global Braginskii edge turbulence code GRILLIX. The code is based on the flux-coordinate independent (FCI) approach, which allows efficient turbulence simulations in diverted equilibria. We validate simulations across the ASDEX Upgrade edge and SOL against measurements in discharge # 36190, and find much better agreement thanks to the neutrals. Disentangling the effects of the neutral gas, we find that it affects the plasma in several ways. Firstly, the ionization of neutrals modifies the radial profiles of plasma density and temperature, leading to a transition of the turbulence drive from the general ballooning type to the ion temperature gradient type. Secondly, strong poloidal asymmetries can be induced due to divertor recycling, depending on the ionization pattern. As ballooned perpendicular plasma transport is stronger at the low-field side, neutrals penetrate deeper into the plasma at the high-field side, leading to significant ionization and radiation there. With increasing divertor neutrals density, the targets cool down while plasma density increases, more strongly at the high-field side. Much of the dynamics take place directly around the X -point and along the separatrix, which can be resolved by the FCI approach. Potential remains in extending the model and the code, but our results build confidence that predictive capability is within reach for major design questions for fusion reactors, such as the near SOL fall-off length. [ABSTRACT FROM AUTHOR]

Abstract A systematic analysis of the discrete conservation properties of non-dissipative, central-difference approximations of the compressible Navier–Stokes equations is reported. A generalized splitting of the nonlinear convective terms is considered, and energy-preserving formulations are fully characterized by deriving a two-parameter family of split forms. Previously developed formulations reported in literature are shown to be particular members of this family; novel splittings are introduced and discussed as well. Furthermore, the conservation properties yielded by different choices for the energy equation (i.e. total and internal energy, entropy) are analyzed thoroughly. It is shown that additional preserved quantities can be obtained through a suitable adaptive selection of the split form within the derived family. Local conservation of primary invariants, which is a fundamental property to build high-fidelity shock-capturing methods, is also discussed in the paper. Numerical tests performed for the Taylor–Green Vortex at zero viscosity fully confirm the theoretical findings, and show that a careful choice of both the splitting and the energy formulation can provide remarkably robust and accurate results. Highlights • Energy-preserving split forms in compressible flow equations are studied. • A novel two-parameter family of energy-preserving splittings is derived. • A dynamic splitting procedure with optimal conservation properties is proposed. [ABSTRACT FROM AUTHOR]

The application of pseudo-symplectic Runge–Kutta methods to the incompressible Navier–Stokes equations is discussed in this work. In contrast to fully energy-conserving, implicit methods, these are explicit schemes of order p that preserve kinetic energy to order q , with q > p . Use of explicit methods with improved energy-conservation properties is appealing for convection-dominated problems, especially in case of direct and large-eddy simulation of turbulent flows. A number of pseudo-symplectic methods are constructed for application to the incompressible Navier–Stokes equations and compared in terms of accuracy and efficiency by means of numerical simulations. [ABSTRACT FROM AUTHOR]

Results from numerical simulations of the mixing of two puffs of scalars released in a turbulent flow channel are used to introduce a measure of mixing quality, and to investigate the effectiveness of turbulent mixing as a function of the location of the puff release and the molecular diffusivity of the puffs. The puffs are released from instantaneous line sources in the flow field with Schmidt numbers that range from 0.7 to 2400. The line sources are located at different distances from the channel wall, starting from the wall itself, the viscous wall layer, the logarithmic layer, and the channel center. The mixing effectiveness is quantified by following the trajectories of individual particles with a Lagrangian approach and carefully counting the number of particles from both puffs that arrive at different locations in the flow field as a function of time. A new measure, the mixing quality index Ø, is defined as the product of the normalized fraction of particles from the two puffs at a flow location. The mixing quality index can take values from 0, corresponding to no mixing, to 0.25, corresponding to full mixing. The mixing quality in the flow is found to depend on the Schmidt number of the puffs when the two puffs are released in the viscous wall region, while the Schmidt number is not important for the mixing of puffs released outside the logarithmic region.

We investigate statistically steady states of turbulent flows when molecular degrees of freedom, in particular vibration, are taken into account. Unlike laminar flows initially in thermal non-equilibrium which asymptotically relax towards thermal equilibrium, turbulent flows present persistent departures from thermal equilibrium. This is due to fluctuations in temperature and other thermodynamic variables, which are known to increase with turbulent Mach number. Analytical results are compared to direct numerical simulations at a range of Reynolds and Mach numbers as well as molecular parameters such as relaxation times. Turbulent fluctuations are also shown to disrupt the distribution of energy between translational-rotational-vibrational modes even if thermal equilibrium is attained instantaneously relative to turbulence time scales, an effect that increases with characteristic relaxation times. Because of the nonlinear relation between temperature and vibrational energy in equilibrium, the fluctuation of the latter could be strongly positively skewed with long tails in its probability density function. This effect is stronger in flows with strong temperature fluctuations and when vibrational modes are partially excited. Because of the finite-time relaxation of vibration, departures from equilibrium result in very strong transfers of energy from the translational-rotational mode to the vibrational mode. A simple spectral model can explain the stronger departures from thermal equilibrium observed at the small scales. The spectral behaviour of the instantaneous vibrational energy can be described by classical phenomenology for passive scalars. [ABSTRACT FROM AUTHOR]

The effects of rough surfaces on turbulent channel flow are modelled by an extra force term in the Navier–Stokes equations. This force term contains two parameters, related to the density and the height of the roughness elements, and a shape function, which regulates the influence of the force term with respect to the distance from the channel wall. This permits a more flexible specification of a rough surface than a single parameter such as the equivalent sand grain roughness. The effects of the roughness force term on turbulent channel flow have been investigated for a large number of parameter combinations and several shape functions by direct numerical simulations. It is possible to cover the full spectrum of rough flows ranging from hydraulically smooth through transitionally rough to fully rough cases. By using different parameter combinations and shape functions, it is possible to match the effects of different types of rough surfaces. Mean flow and standard turbulence statistics have been used to compare the results to recent experimental and numerical studies and a good qualitative agreement has been found. Outer scaling is preserved for the streamwise velocity for both the mean profile as well as its mean square fluctuations in all but extremely rough cases. The structure of the turbulent flow shows a trend towards more isotropic turbulent states within the roughness layer. In extremely rough cases, spanwise structures emerge near the wall and the turbulent state resembles a mixing layer. A direct comparison with the study of Ashrafian, Andersson & Manhart (Intl J. Heat Fluid Flow, vol. 25, 2004, pp. 373–383) shows a good quantitative agreement of the mean flow and Reynolds stresses everywhere except in the immediate vicinity of the rough wall. The proposed roughness force term may be of benefit as a wall model for direct and large-eddy numerical simulations in cases where the exact details of the flow over a rough wall can be neglected. [ABSTRACT FROM PUBLISHER]

The Lagrangian block hydrodynamics is formulated based on the block advection of fluid. By enforcing the mass and momentum conservations on the Lagrangian mesh, the numerical oscillation problem encountered in the classical Eulerian computational methods is circumvented. A large number of the previously computationally difficult problems in environmental fluid mechanics are successfully simulated using the method. Examples of these simulations are described in this paper. [ABSTRACT FROM AUTHOR]

We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer zones away from the midplane in order to assess the effects of boundaries. The dynamical α quenching phenomenology requires that the magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the small-scale magnetic helicity is found to reach a steady state. In agreement with the earlier work, the magnetic helicity fluxes of small-scale fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that 'softer' (halo) boundary conditions give a lower energy of the saturated mean magnetic field. [ABSTRACT FROM AUTHOR]

Turbulence affects the structure and motions of nearly all temperature and density regimes in the interstellar gas. This two-part review summarizes the observations, theory, and simulations of interstellar turbulence and their implications for many fields of astrophysics. The first part begins with diagnostics for turbulence that have been applied to the cool interstellar medium and highlights their main results. The energy sources for interstellar turbulence are then summarized along with numerical estimates for their power input. Supernovae and superbubbles dominate the total power, but many other sources spanning a large range of scales, from swing-amplified gravitational instabilities to cosmic ray streaming, all contribute in some way. Turbulence theory is considered in detail, including the basic fluid equations, solenoidal and compressible modes, global inviscid quadratic invariants, scaling arguments for the power spectrum, phenomenological models for the scaling of higher-order structure functions, the direction and locality of energy transfer and cascade, velocity probability distributions, and turbulent pressure. We emphasize expected differences between incompressible and compressible turbulence. Theories of magnetic turbulence on scales smaller than the collision mean free path are included, as are theories of magneto-hydrodynamic turbulence and their various proposals for power spectra. Numerical simulations of interstellar turbulence are reviewed. Models have reproduced the basic features of the observed scaling relations, predicted fast decay rates for supersonic MHD turbulence, and derived probability distribution functions for density. Thermal instabilities and thermal phases have a new interpretation in a supersonically turbulent medium. Large-scale models with various combinations of self-gravity, magnetic fields, supernovae, and star formation are beginning to resemble the observed interstellar medium in morphology and... [ABSTRACT FROM AUTHOR]

We study a one-dimensional discrete analog of the von Kármán flow widely investigated in turbulence, made of a lattice of anharmonic oscillators excited by both ends in the presence of a dissipative term proportional to the second-order finite difference of the velocities, similar to the viscous term in a fluid. The dynamics of the model shows striking similarities with an actual turbulent flow, both at local and global scales. Calculations of the probability distribution function of velocity increments, extensively studied in turbulence, with a very large number of points in order to determine accurately the statistics of rare events, allow us to provide a meaningful comparison of different theoretical expressions of the PDFs. [Copyright &y& Elsevier]

• Methods to impose virtual origins on the velocities in turbulent direct simulations are studied. • Robin conditions for all three components are identified as a suitable method to impose such origins. • The origin for near-wall turbulence is controlled by the spanwise and wall-normal virtual origins. • The origin for the mean velocity profile, in turn, is set by the streamwise origin. The relative wall-normal displacement of the origin perceived by different components of near-wall turbulence is known to produce a change in drag. This effect is produced for instance by drag-reducing surfaces of small texture-size like riblets and superhydrophobic surfaces. To facilitate the research on how these displacements alter near-wall turbulence, this paper studies different strategies to model such displacement effect through manipulated boundary conditions. Previous research has considered the effect of offsetting the virtual origins perceived by the tangential components of the velocity from the reference, boundary plane, where the wall-normal velocity was set to zero. These virtual origins are typically characterised by slip-length coefficients in Robin, slip-like boundary conditions. In this paper, we extend this idea and explore several techniques to define and implement virtual origins for all three velocity components on direct numerical simulations (DNSs) of channel flows, with special emphasis on the wall-normal velocity. The aim of this work is to provide a suitable foundation to extend the existing understanding on how these virtual origins affect the near-wall turbulence, and ultimately aid in the formulation of simplified models that capture the effect of complex surfaces on the overlying flow and on drag, without the need to resolve fully the turbulence and the surface texture. From the techniques tested, Robin boundary conditions for all three velocities are found to be the most satisfactory method to impose virtual origins, relating the velocity components to their respective wall-normal gradients linearly. Our results suggest that the effect of virtual origins on the flow, and hence the change in drag that they produce, can be reduced to an offset between the virtual origin perceived by the mean flow and that perceived by the overlying turbulence, and that turbulence remains otherwise smooth-wall-like, as proposed by Luchini (1996). The origin for turbulence, however, would not be set by the spanwise virtual origin alone, but by a combination of the spanwise and wall-normal origins. These observations suggest the need for an extension of Luchini's virtual-origin theory to predict the change in drag, accounting for the wall-normal transpiration when its effect is not negligible. [ABSTRACT FROM AUTHOR]

• B-spline-based filters are provably dissipative for periodic and bounded domains. • A dynamic procedure outperforms a priori selection of filtering parameters. • An efficient implementation algorithm makes cost comparable to compact (implicit) finite difference filters. A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. These simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy. [ABSTRACT FROM AUTHOR]