Search

In this manuscript, it is analysed the onset and evolution of natural convection of an incompressible fluid between spherical shells. The shells are kept at a fixed temperature difference and aspect ratio, and the Rayleigh-Benard convection is driven by different radial gravity profiles. The analysis has been carried out by using a finite difference scheme to solve the three-dimensional Navier-Stokes equations in spherical coordinates. Numerical results are compared with theoretical predictions from linear and non-linear stability analysis, and differ from the expected critical Rayleigh number Ra_c = 1708 by less than 1 percent. In the range of Prandlt numbers Pr studied, and for all the different gravity profiles analysed, the system presents a dependence on its starting condition and flow history. Even in the region just above the onset of convection, two stable states are observed, with qualitative and quantitative differences, and exploring higher values of Ra introduces new modes and time dependency phenomena in the flow. These results are corroborated by spectral analysis.

Transport equations for heavy inertial particles in turbulent boundary layers may be derived from an underlying phase-space probability density function (PDF) equation. These equations, however, are unclosed, and the standard closure approach is to use a quasi-Normal approximation (QNA) in which the fourth moments are approximated as behaving as if the velocities were Normally distributed. Except for particles with weak inertia, the QNA leads to large quantitative errors, and is not consistent with the known asymptotic predictions of Sikovsky (Flow Turbulence Combust, vol. 92, 2014, pp. 41-64) for the moments of the PDF in the viscous sublayer. We derive a new closure approximation based on an asymptotic solution to the transport equations in regions where the effect of particle inertia is significant. The new closure is consistent with the asymptotic predictions of Sikovsky, but applies even outside the viscous sublayer. Comparisons with direct numerical simulations (DNS) show that the new closure gives similar results to the QNA (with the QNA results in slightly better agreement with the DNS) when the viscous Stokes number is $St<10$, but for $St>10$ the new model is in far better agreement with the DNS than the QNA. While the predictions from the new model leave room for improvement, the results suggest that this new closure strategy is a very effective alternative to the traditional QNA approach, and the closure could be refined in future work.

Inspired by the hemispherical asymmetry observed in the Earth's inner core, we perform direct numerical simulations to study the effect of the gravity center offset in spherical Rayleigh-B\'enard convection. We find that even a minimal shift of the gravity center has a pronounced influence on the flow structures. When the gravity center is shifted towards the South, the co-latitudinal buoyancy component creates an energetic jet on the Northern side of the inner sphere that is directed towards the outer sphere. As a result, a large-scale meridional circulation is formed. However, surprisingly, the global heat flux is not affected by the changes in the large-scale flow organization induced by the gravity center offset. Our results suggest that the hemispherical core asymmetry is key to model the flow phenomena in the Earth's outer core and mantle.

Direct numerical simulations are employed to reveal three distinctly different flow regions in rotating spherical Rayleigh-B\'enard convection. In the low-latitude region $\mathrm{I}$ vertical (parallel to the axis of rotation) convective columns are generated between the hot inner and the cold outer sphere. The mid-latitude region $\mathrm{II}$ is dominated by vertically aligned convective columns formed between the Northern and Southern hemispheres of the outer sphere. The diffusion-free scaling, which indicates bulk-dominated convection, originates from this mid-latitude region. In the equator region $\mathrm{III}$ the vortices are affected by the outer spherical boundary and are much shorter than in region $\mathrm{II}$. Thermally driven turbulence with background rotation in spherical Rayleigh-B\'enard convection is found to be characterized by three distinctly different flow regions. The diffusion-free scaling, which indicates the heat transfer is bulk-dominated, originates from the mid-latitude region in which vertically aligned vortices are stretched between the Northern and Southern hemispheres of the outer sphere.\ These results show that the flow physics in rotating convection are qualitatively different in planar and spherical geometries. This finding underlines that it is crucial to study convection in spherical geometries to better understand geophysical and astrophysical flow phenomena., Comment: 11pages, 6figures

We explore the role of gravitational settling on inertial particle concentrations in a wall-bounded turbulent flow. While it may be thought that settling can be ignored when the settling parameter $Sv\equiv v_s/u_\tau$ is small ($v_s$ - Stokes settling velocity, $u_\tau$ - fluid friction velocity), we show that even in this regime the settling may make a leading order contribution to the concentration profiles. This is because the importance of settling is determined, not by the size of $v_s$ compared with $u_\tau$ or any other fluid velocity scale, but by the size of $v_s$ relative to the other mechanisms that control the vertical particle velocity and concentration profile. We explain this in the context of the particle mean-momentum equation, and show that in general, there always exists a region in the boundary layer where settling cannot be neglected, no matter how small $Sv$ is (provided it is finite). Direct numerical simulations confirm the arguments, and show that the near-wall concentration is highly dependent on $Sv$ even when $Sv\ll 1$, and can reduce by an order of magnitude when $Sv$ is increased from $O(10^{-4})$ and $O(10^{-2})$. The results also show that the preferential sampling of ejection events in the boundary layer by inertial particles when $Sv=0$ is profoundly altered as $Sv$ is increased, and is replaced by a preferential sampling of sweep events due to the onset of the preferential sweeping mechanism.

We use theory and Direct Numerical Simulations (DNS) to explore the average vertical velocities and spatial distributions of inertial particles settling in a wall-bounded turbulent flow. The theory is based on the exact phase-space equation for the Probability Density Function describing particle positions and velocities. This allowed us to identify the distinct physical mechanisms governing the particle transport. We then examined the asymptotic behavior of the particle motion near the wall, revealing the fundamental differences to the near wall behavior that is produced when incorporating gravitational settling. When the average vertical particle mass flux is zero, the averaged vertical particle velocity is zero away from the wall due to the particles preferentially sampling regions where the fluid velocity is positive, which balances with the downward Stokes settling velocity. When the average mass flux is negative, the combined effects of turbulence and particle inertia lead to average vertical particle velocities that can significantly exceed the Stokes settling velocity, by as much as ten times. Sufficiently far from the wall, the enhanced vertical velocities are due to the preferential sweeping mechanism. However, as the particles approach the wall, the contribution from the preferential sweeping mechanism becomes small, and a downward contribution from the turbophoretic velocity dominates the behavior. Close to the wall, the particle concentration grows as a power-law, but the nature of this power law depends on the particle Stokes number. Finally, our results highlight how the Rouse model of particle concentration is to be modified for particles with finite inertia.

Direct numerical simulations two-way coupled with inertial particles are used to investigate the particle distribution and two-way coupling effect of low-inertia ($St_{LSM}=0.0625$, $St_{VLSM}=0.009$) and high-inertia ($St_{LSM}=0.475$, $St_{VLSM}=0.069$) particles associated with the large-scale motions (LSMs) and very-large-scale motions (VLSMs) in an open channel flow at a Reynolds number of $Re_\tau = 550$. One method of filtering the VLSMs from the flow is via artificial domain truncation, which alters the mean particle concentration profile and particle clustering due to the removal of VLSMs from a large domain simulation. In order to exclude possible correlation of the turbulence introduced by a small domain size with periodic boundary conditions, low- and high-pass filtering is performed during the simulation to isolate the particle interaction with different spatial scales. The results show that particle accumulation and turbophoresis are under-predicted without VLSMs, whereas the particle clustering and two-way coupling effects are mainly determined by particle coupling with LSMs. In the inner layer, the elongated streamwise anisotropic particle clustering can be reproduced by particles coupling solely with LSMs for low Stokes number ($St_{LSM}=0.0625$) particles. However, we do not observe similar particle clustering behavior in the outer layer as seen in the full simulation by coupling particles with either LSMs or VLSMs for high Stokes number ($St_{VLSM}=0.069$) particles. This indicates that the organized particle structures are formed by the joint action of LSMs and VLSMs, especially for high Stokes number particles in the outer layer.

This study is concerned with the statistics of vertical turbulent channel flow laden with inertial particles for two different volume concentrations ($\Phi_{V} = 3 \times 10^{-6}$ and $\Phi_{V} = 5 \times 10^{-5}$) at a Stokes number of $St^{+} = 58.6$ based on viscous units. Two independent direct numerical simulation models utilizing the point-particle approach are compared to recent experimental measurements, where all relevant nondimensional parameters are directly matched. While both numerical models are built on the same general approach, details of the implementations are different, particularly regarding how two-way coupling is represented. At low volume loading, both numerical models are in general agreement with the experimental measurements, with certain exceptions near the walls for the wall-normal particle velocity fluctuations. At high loading, these discrepancies are increased, and it is found that particle clustering is overpredicted in the simulations as compared to the experimental observations. Potential reasons for the discrepancies are discussed. As this study is among the first to perform one-to-one comparisons of particle-laden flow statistics between numerical models and experiments, it suggests that continued efforts are required to reconcile differences between the observed behavior and numerical predictions.

Very large-scale motions (VLSMs) and large-scale motions (LSMs) coexist at moderate Reynolds numbers in a very long open channel flow. Direct numerical simulations two-way coupled with inertial particles are analysed using spectral information to investigate the modulation of VLSMs. In the wall-normal direction, particle distributions (mean/preferential concentration) exhibit two distinct behaviors in the inner flow and outer flow, corresponding to two highly anisotropic turbulent structures, LSMs and VLSMs. This results in particle inertia's non-monotonic effects on the VLSMs: low inertia (based on the inner scale) and high inertia (based on the outer scale) both strengthen the VLSMs whereas moderate and very high inertia have little influence. Through conditional tests, low and high inertia particles enhance VLSMs following two distinct routes. Low inertia particles promote VLSMs indirectly through the enhancement of the regeneration cycle (the self-sustaining mechanism of LSMs) in the inner region whereas high inertia particles enhance the VLSM directly through contribution to the Reynolds shear stress at similar temporal scales in the outer region. This understanding also provides more general insight into inner-outer interaction in high Reynolds number, wall-bounded flows.

Two-way coupled direct numerical simulations are used to investigate the effects of inertial particles on self-sustained, turbulent coherent structures (i.e. the so-called the regeneration cycle) in plane Couette flow at low Reynolds number just above the onset of transition. Tests show two limiting behaviors with increasing particle inertia, similar to the results from the linear stability analysis of saffman, 1962: low-inertia particles trigger the laminar-to-turbulent instability whereas high-inertia particles tend to stabilize turbulence due to the extra dissipation induced by particle-fluid coupling. Furthermore, it is found that the streamwise coupling between phases is the dominant factor in damping the turbulence and is highly related to the spatial distribution of the particles. The presence of particles in different turbulent coherent structures (large scale vortices or large scale streaks) determines the turbulent kinetic energy of particulate phase, which is related to the particle response time scaled by the turnover time of large scale vortices. By quantitatively investigating the periodic character of the whole regeneration cycle and the phase difference between linked sub-steps, we show that the presence of inertial particles does not alter the periodic nature of the cycle or the relative length of each of the sub-steps. Instead, high-inertia particles greatly weaken the large-scale vortices as well as the streamwise vorticity stretching and lift-up effects, thereby suppressing the fluctuating amplitude of the large scale streaks. The primary influence of low-inertia particles, however, is to strengthen the large scale vortices, which fosters the cycle and ultimately reduces the critical Reynolds number.

Direct numerical simulations of turbulent suspension flows are carried out with the Force-Coupling Method in plane Couette and pressure-driven channel configurations. Dilute to moderately concentrated suspensions of neutrally buoyant finite-size (L_y /d = 20) spherical particles are considered when the Reynolds number is slightly above the laminar-turbulent transition. Tests performed with synthetic streaks, in both turbulent channel and Couette flows, show clearly that particles trigger the instability in channel flow whereas the plane Couette flow becomes laminar. Moreover, we have shown that particles have a pronounced impact on pressure-driven flow through a detailed temporal and spatial analysis whereas they have no significant impact on plane Couette flow configuration. The substantial difference between both flows is related to the spatial preferential distribution of particles in the large-scale rolls (inactive motion) in Couette flow, whereas they are accumulated in the ejection (active motion) regions in the pressure-driven flow. Through investigation of particle modification on two distinct flow configurations, we were able to show the specific response of turbulent structures and the modulation of the fundamental mechanisms composing the regeneration cycle in the buffer layer of near-wall turbulence. Especially for pressure-driven flow, the particles enhance the lift- up and let it act continuously whereas the particles do not significantly alter the streak breakdown process. The reinforcement of the streamwise vortices is attributed to the vorticity stretching term by the wavy streaks. The smaller and more numerous wavy streaks enhance the vorticity stretching and consequently strengthen the circulation of the large-scale streamwise vortex in suspension flow., Comment: Submitted to JFM and under consideration