Search

In gas evolving electrolysis, bubbles grow at electrodes due to a diffusive influx from oversaturation generated locally in the electrolyte by the electrode reaction. When considering electrodes of micrometer-size resembling catalytic islands, bubbles are found to approach dynamic equilibrium states at which they neither grow nor shrink. Such equilibrium states are found at low oversaturation for both, pinning and expanding wetting regimes of the bubbles and are based on the balance of local influx near the bubble foot and global outflux. Unlike the stability of pinned nano-bubbles studied earlier, the Laplace pressure plays a minor role only. We extend the analytical solution of Zhang & Lohse (2023) by taking into account the non-uniform distribution of dissolved gas around the bubble obtained from direct numerical simulation. This allows us to identify the parameter regions of bubble growth, dissolution and dynamic equilibrium as well as to analyze the stability of the equilibrium states. Finally, we draw conclusions on how to possibly enhance the efficiency of electrolysis., Comment: 16 pages, 13 figures

The distributed telescope array offers promise for conducting large-sky-area, high-frequency time domain surveys. Multiple telescopes can be deployed at each observation site, so intra-site observation task scheduling is crucial for enhancing observation efficiency and quality. Efficient use of observable time and rapid response to special situations are critical to maximize scientific discovery in time domain surveys. Besides, the competing scientific priorities, time-varying observation conditions, and capabilities of observation equipment, lead to a vast search space of the scheduling. So with the increasing number of telescopes and observation fields, balancing computational time with solution quality in observation scheduling poses a significant challenge. Informed by the seminal contributions of earlier studies on a multilevel scheduling model and global scheduler for time domain telescope array, this study is devoted to further exploring the site scheduler. Formulating the observation scheduling of multiple telescopes at the site as a cooperative decision-making problem, this paper proposes GRRIS, a real-time intra-site observation scheduling scheme for telescope array using graph and reinforcement learning. It employs a graph neural network to learn node features that can embed the spatial structure of the observation scheduling. An algorithm based on multi-agent reinforcement learning is designed to efficiently learn the optimum allocation policy of telescope agents to field nodes. Through numerical simulations with real-world scenarios, GRRIS can achieve up to a 22% solution improvement over the most competitive scheme. It offers better scalability and sub-second decision speed, meeting the needs of observation scheduling control for future distributed telescope arrays., Comment: Astronomical Journal

Large Language Models (LLMs) have demonstrated remarkable performance across diverse domains, prompting researchers to explore their potential for use in recommendation systems. Initial attempts have leveraged the exceptional capabilities of LLMs, such as rich knowledge and strong generalization through In-context Learning, which involves phrasing the recommendation task as prompts. Nevertheless, the performance of LLMs in recommendation tasks remains suboptimal due to a substantial disparity between the training tasks for LLMs and recommendation tasks and inadequate recommendation data during pre-training. This paper introduces RLRF4Rec, a novel framework integrating Reinforcement Learning from Recsys Feedback for Enhanced Recommendation Reranking(RLRF4Rec) with LLMs to address these challenges. Specifically, We first have the LLM generate inferred user preferences based on user interaction history, which is then used to augment traditional ID-based sequence recommendation models. Subsequently, we trained a reward model based on knowledge augmentation recommendation models to evaluate the quality of the reasoning knowledge from LLM. We then select the best and worst responses from the N samples to construct a dataset for LLM tuning. Finally, we design a structure alignment strategy with Direct Preference Optimization(DPO). We validate the effectiveness of RLRF4Rec through extensive experiments, demonstrating significant improvements in recommendation re-ranking metrics compared to baselines. This demonstrates that our approach significantly improves the capability of LLMs to respond to instructions within recommender systems.

We conducted a series of pore-scale numerical simulations on convective flow in porous media, with a fixed Schmidt number of 400 and a wide range of Rayleigh numbers. The porous media are modeled using regularly arranged square obstacles in a Rayleigh-B\'enard (RB) system. As the Rayleigh number increases, the flow transitions from a Darcy-type regime to an RB-type regime, with the corresponding $Sh-Ra_D$ relationship shifting from sublinear scaling to the classical 0.3 scaling of RB convection. Here, $Sh$ and $Ra_D$ represent the Sherwood number and Rayleigh-Darcy number, respectively. For different porosities, the transition begins at approximately $Ra_D = 4000$, at which point the characteristic horizontal scale of the flow field is comparable to the size of a single obstacle unit. When the thickness of the concentration boundary layer is less than about one-sixth of the pore spacing, the flow finally enters the RB regime. In the Darcy regime, the scaling exponent of $Sh$ and $Ra_D$ decreases as porosity increases. Based on the Grossman-Lohse theory (J. Fluid Mech., vol. 407, 2000; Phys. Rev. Lett., vol. 86, 2001), we provide an explanation for the scaling laws in each regime and highlight the significant impact of mechanical dispersion effects within the boundary layer. Our findings provide some new insights into the validity range of the Darcy model.

Objective: To explore and compare the performance of ChatGPT and other state-of-the-art LLMs on domain-specific NER tasks covering different entity types and domains in TCM against COVID-19 literature. Methods: We established a dataset of 389 articles on TCM against COVID-19, and manually annotated 48 of them with 6 types of entities belonging to 3 domains as the ground truth, against which the NER performance of LLMs can be assessed. We then performed NER tasks for the 6 entity types using ChatGPT (GPT-3.5 and GPT-4) and 4 state-of-the-art BERT-based question-answering (QA) models (RoBERTa, MiniLM, PubMedBERT and SciBERT) without prior training on the specific task. A domain fine-tuned model (GSAP-NER) was also applied for a comprehensive comparison. Results: The overall performance of LLMs varied significantly in exact match and fuzzy match. In the fuzzy match, ChatGPT surpassed BERT-based QA models in 5 out of 6 tasks, while in exact match, BERT-based QA models outperformed ChatGPT in 5 out of 6 tasks but with a smaller F-1 difference. GPT-4 showed a significant advantage over other models in fuzzy match, especially on the entity type of TCM formula and the Chinese patent drug (TFD) and ingredient (IG). Although GPT-4 outperformed BERT-based models on entity type of herb, target, and research method, none of the F-1 scores exceeded 0.5. GSAP-NER, outperformed GPT-4 in terms of F-1 by a slight margin on RM. ChatGPT achieved considerably higher recalls than precisions, particularly in the fuzzy match. Conclusions: The NER performance of LLMs is highly dependent on the entity type, and their performance varies across application scenarios. ChatGPT could be a good choice for scenarios where high recall is favored. However, for knowledge acquisition in rigorous scenarios, neither ChatGPT nor BERT-based QA models are off-the-shelf tools for professional practitioners., Comment: 22 pages with 2 figures

Infrastructure sensors installed at elevated positions offer a broader perception range and encounter fewer occlusions. Integrating both infrastructure and ego-vehicle data through V2X communication, known as vehicle-infrastructure cooperation, has shown considerable advantages in enhancing perception capabilities and addressing corner cases encountered in single-vehicle autonomous driving. However, cooperative perception still faces numerous challenges, including limited communication bandwidth and practical communication interruptions. In this paper, we propose CTCE, a novel framework for cooperative 3D object detection. This framework transmits queries with temporal contexts enhancement, effectively balancing transmission efficiency and performance to accommodate real-world communication conditions. Additionally, we propose a temporal-guided fusion module to further improve performance. The roadside temporal enhancement and vehicle-side spatial-temporal fusion together constitute a multi-level temporal contexts integration mechanism, fully leveraging temporal information to enhance performance. Furthermore, a motion-aware reconstruction module is introduced to recover lost roadside queries due to communication interruptions. Experimental results on V2X-Seq and V2X-Sim datasets demonstrate that CTCE outperforms the baseline QUEST, achieving improvements of 3.8% and 1.3% in mAP, respectively. Experiments under communication interruption conditions validate CTCE's robustness to communication interruptions., Comment: Accepted by IEEE ITSC 2024

Solidification of droplets is of great importance to various technological applications, drawing considerable attention from scientists aiming to unravel the fundamental physical mechanisms. In the case of multicomponent droplets undergoing solidification, the emergence of concentration gradients may trigger significant interfacial flows that dominate the freezing dynamics. Here, we experimentally investigate the fascinating snow-globe freezing dynamics of supercooled ethanol-water droplets. We reveal that these unique freezing dynamics are driven by solidification-induced solutal Marangoni flow within the droplets. We quantitatively characterize the concentration-dependent migration and growth dynamics of ice particles, tightly connecting them to the solutal Marangoni effect and the associated convective heat transfer. Moreover, we show that the final wrapping state of droplets can be modulated by the concentration of ethanol. Our findings may pave the way for novel insights into the physicochemical hydrodynamics of multicomponent liquids undergoing phase transitions., Comment: 5 pages, 4 figures

We investigate the coupling effects of the two-phase interface, viscosity ratio, and density ratio of the dispersed phase to the continuous phase on the flow statistics in two-phase Taylor-Couette turbulence at a system Reynolds number of 6000 and a system Weber number of 10 using interface-resolved three-dimensional direct numerical simulations with the volume-of-fluid method. Our study focuses on four different scenarios: neutral droplets, low-viscosity droplets, light droplets, and low-viscosity light droplets. We find that neutral droplets and low-viscosity droplets primarily contribute to drag enhancement through the two-phase interface, while light droplets reduce the system's drag by explicitly reducing Reynolds stress due to the density dependence of Reynolds stress. Additionally, low-viscosity light droplets contribute to greater drag reduction by further reducing momentum transport near the inner cylinder and implicitly reducing Reynolds stress. While interfacial tension enhances turbulent kinetic energy (TKE) transport, drag enhancement is not strongly correlated with TKE transport for both neutral droplets and low-viscosity droplets. Light droplets primarily reduce the production term by diminishing Reynolds stress, whereas the density contrast between the phases boosts TKE transport near the inner wall. Therefore, the reduction in the dissipation rate is predominantly attributed to decreased turbulence production, causing drag reduction. For low-viscosity light droplets, the production term diminishes further, primarily due to their greater reduction in Reynolds stress, while reduced viscosity weakens the density difference's contribution to TKE transport near the inner cylinder, resulting in a more pronounced reduction in the dissipation rate and consequently stronger drag reduction. Our findings provide new insights into the turbulence modulation in two-phase flow.

The radar signal processing algorithm is one of the core components in through-wall radar human detection technology. Traditional algorithms (e.g., DFT and matched filtering) struggle to adaptively handle low signal-to-noise ratio echo signals in challenging and dynamic real-world through-wall application environments, which becomes a major bottleneck in the system. In this paper, we introduce an end-to-end through-wall radar human detection network (TWP-CNN), which takes raw radar Analog-to-Digital Converter (ADC) signals without any preprocessing as input. We replace the conventional radar signal processing flow with the proposed DFT-based adaptive feature extraction (DAFE) module. This module employs learnable parameterized 3D complex convolution layers to extract superior feature representations from ADC signals, which is beyond the limitation of traditional preprocessing methods. Additionally, by embedding phase information from radar data within the network and employing multi-task learning, a more accurate detection is achieved. Finally, due to the absence of through-wall radar datasets containing raw ADC data, we gathered a realistic through-wall (RTW) dataset using our in-house developed through-wall radar system. We trained and validated our proposed method on this dataset to confirm its effectiveness and superiority in real through-wall detection scenarios., Comment: 11pages,13figures

Understanding the solidification dynamics of impacted water droplets is fundamental and crucial for applications, especially with the presence of frost and salt. Here, we experimentally investigate the spreading and freezing dynamics of saline droplets upon impact on a cold, frosty surface. Our findings demonstrate that the frost and salt can lead to an extended spreading of impacted droplets under specific conditions. In addition to the well-known 1/2 inertial spreading scaling law for droplets impact on a cold substrate, we observe a distinct transition to a 1/10 scaling law dominated by the capillary-viscous relation at low impacting velocities. We formulate the onset criterion for this extended spreading regime and derive a scaling dependence that captures the droplet's arrested diameter over various supercooling temperatures, by incorporating the effect of impact inertia, partial-wetting behavior and salinity. Based on the analysis, a unified model is proposed for predicting the droplet arrested diameter over a wide range of impact velocity and salinity., Comment: 10 pages, 10 figures

The fields of therapeutic application and drug research and development (R&D) both face substantial challenges, i.e., the therapeutic domain calls for more treatment alternatives, while numerous promising pre-clinical drugs have failed in clinical trials. One of the reasons is the inadequacy of Cross-drug Response Evaluation (CRE) during the late stages of drug R&D. Although in-silico CRE models bring a promising solution, existing methodologies are restricted to early stages of drug R&D, such as target and cell-line levels, offering limited improvement to clinical success rates. Herein, we introduce DeepCRE, a pioneering AI model designed to predict CRE effectively in the late stages of drug R&D. DeepCRE outperforms the existing best models by achieving an average performance improvement of 17.7% in patient-level CRE, and a 5-fold increase in indication-level CRE, facilitating more accurate personalized treatment predictions and better pharmaceutical value assessment for indications, respectively. Furthermore, DeepCRE has identified a set of six drug candidates that show significantly greater effectiveness than a comparator set of two approved drugs in 5/8 colorectal cancer organoids. This demonstrates the capability of DeepCRE to systematically uncover a spectrum of drug candidates with enhanced therapeutic effects, highlighting its potential to transform drug R&D.

Evaporation of multicomponent droplets is important in a wide range of applications, albeit complex, and requires a careful investigation. We experimentally and numerically investigate the evaporation characteristics of spherical, ethanol-water droplets with different initial concentration ratios in the acoustic levitation field. Imaging techniques and infrared thermometry are used for acquiring volume and surface temperature variations of droplets, reflecting their mass and heat transfer characteristics. Numerical simulations are conducted using modified parameters based on a theoretical model to consider the effect of the acoustic field. The calculation results show good agreement with the experimental data. The concentration and temperature distribution within the droplet is further investigated based on the numerical results., Comment: 14 pages, 6 figures

We perform the linear stability analysis and direct numerical simulations to study the effect of radius ratio on the instability and flow characteristics of the sheared annular centrifugal Rayleigh-B\'enard convection (ACRBC), where the cold inner cylinder and the hot outer cylinder rotate with a small angular velocity difference. With the shear enhancement, the thermal convection is suppressed and finally gets stable for different radius ratios $\eta\in[0.2, 0.95]$. Considering the inhomogeneous distribution of shear stresses in the base flow, a new global Richardson number $Ri_g$ is defined and the marginal-state curves for different radius ratios are successfully unified in the parameter domain of $Ri_g$ and the Rayleigh number $Ra$. The results are consistent with the marginal-state curve of the wall-sheared classical RBC in the streamwise direction, demonstrating that the basic stabilization mechanisms are identical. Moreover, systems with small radius ratios exhibit greater geometric asymmetry. On the one hand, this results in a smaller equivalent aspect ratio for the system, accommodating fewer convection roll pairs. Fewer roll pairs are more likely to cause a transition in the flow structure during shear enhancement. On the other hand, the shear distribution is more inhomogeneous, allowing for an outward shift of the convection region and the elevation of bulk temperature under strong shear.

Photonic integrated circuits based on erbium doped thin film lithium niobate on insulator has attracted broad interests with insofar various waveguide amplifiers and microlasers demonstrated. Wideband operation facilitated by the broadband absorption and emission of erbium ions necessitates the functional integration of wavelength filter and multiplexer on the same chip. Here a low-loss wavelength division multiplexer at the resonant pumping and emission wavelengths (~1480 nm and 1530~1560 nm) of erbium ions based on angled multimode interferometer, is realized in the erbium doped thin film lithium niobate on insulator fabricated by the photolithography assisted chemomechanical etching technique. The minimum on-chip insertion losses of the fabricated device are <0.7 dB for both wavelength ranges, and a 3-dB bandwidth of >20 nm is measured at the telecom C-band. Besides, direct visualization of the multimode interference pattern by the visible upconversion fluorescence of erbium ions compares well with the simulated light propagation in the multimode interferometer. Spectral tuning of the wavelength division multiplexer by structural design is also demonstrated and discussed., Comment: 11 pages, 5 figures

This work presents the first fabricated electrophysiology-optogenetic closed-loop bidirectional brain-machine interface (CL-BBMI) system-on-chip (SoC) with electrical neural signal recording, on-chip sleep staging and optogenetic stimulation. The first multiplexer with static assignment based table lookup solution (MUXnet) for multiplier-free NN processor was proposed. A state-of-the-art average accuracy of 82.4% was achieved with an energy consumption of only 0.2$\mu$J/class in sleep staging task., Comment: 2 pages, 6 figures. Accepted by IEEE Custom Integrated Circuits Conference (CICC) 2024. The codes for the MUXnet (constructing neural networks using multiplexers instead of multipliers) will be open-sourced after the Journal version of this work is accepted

This paper studies the distributed optimization problem under the influence of heavy-tailed gradient noises. Here, a heavy-tailed noise means that the noise does not necessarily satisfy the bounded variance assumption. Instead, it satisfies a more general assumption. The commonly-used bounded variance assumption is a special case of the considered noise assumption. A typical example of this kind of noise is a Pareto distribution noise with tail index within (1,2], which has infinite variance. Despite that there has been several distributed optimization algorithms proposed for the heavy-tailed noise scenario, these algorithms need a centralized server in the network which collects the information of all clients. Different from these algorithms, this paper considers that there is no centralized server and the agents can only exchange information with neighbors in a communication graph. A distributed method combining gradient clipping and distributed stochastic subgradient projection is proposed. It is proven that when the gradient descent step-size and the gradient clipping step-size meet certain conditions, the state of each agent converges to the optimal solution of the distributed optimization problem with probability 1. The simulation results validate the algorithm.

Superconducting materials exhibit unique physical properties and hold great scientific value and vast industrial application prospects. However, due to limitations such as critical temperature (TC) and critical current density (JC), the large-scale application of superconducting materials remains challenging. Chemical doping has been a commonly used method to enhance the superconductivity of B(P)SCCO. However, satisfactory enhancement results have been difficult to achieve. In this study, we introduced green light GaN p-n junction particles as inhomogeneous phases into B(P)SCCO polycrystalline particles to form a smart meta-superconductors (SMSCs) structure. Based on the electroluminescence properties of the p-n junction, the Cooper pairs were stimulated and strengthened to enhance the superconductivity of B(P)SCCO. Experimental results demonstrate that the introduction of inhomogeneous phases can indeed enhance the critical temperature TC, critical current density JC, and complete diamagnetism (Meissner effect) of B(P)SCCO superconductors. Moreover, When the particle size of raw material of B(P)SCCO is reduced from 30{\mu}m to 5{\mu}m, the grain size of the sintered samples also decreases, and the optimal doping concentration of the inhomogeneous phases increases from 0.15 wt.% to 0.2 wt.%, further improving the enhancement of superconductivity., Comment: 17 pages,7 figures

Understanding the dynamics of material objects advected by turbulent flows is a long standing question in fluid dynamics. In this perspective article we focus on the characterization of the statistical properties of non-interacting finite-sized massive spherical particles advected by a vigorous turbulent flow. We study the fluctuations and temporal correlations of particle accelerations and explore their behaviours with respect to the particle size and the particle mass density by means of fully-resolved numerical simulations. We observe that the measured trends can not be interpreted as the simple multiplicative combination of the two dominant effects: the spatial filtering of fluid accelerations and the added-mass-adjusted fluid-to-particle density ratio. We argue that other hydrodynamical forces or effects, e.g. preferential flow sampling, have still a significant role even at the largest particle sizes, which rich here the integral scale of turbulence., Comment: 7 pages, 4 figures

Direct numerical simulations are performed to explore the effects of rotating direction of the vertical asymmetric rough wall on the transport properties of Taylor-Couette (TC) flow up to a Taylor number of $\textit{Ta} = 2.39 \times 10^7$. It is shown that compared to the smooth wall, the rough wall with vertical asymmetric strips can enhance the dimensionless torque \textit{Nu}$_\omega$, and more importantly, at high \textit{Ta} clockwise rotation of the inner rough wall (the fluid is sheared by the steeper slope side of the strips) results in a significantly bigger torque enhancement as compared to the counter-clockwise rotation (the fluid is sheared by the smaller slope side of the strips) due to the larger convective contribution to the angular velocity flux, although the rotating direction has a negligible effect on the torque at low \textit{Ta}. The larger torque enhancement caused by the clockwise rotation of vertical asymmetric rough wall at high \textit{Ta} is then explained by the stronger coupling between the rough wall and the bulk due to the larger biased azimuthal velocity towards the rough wall at the mid-gap of TC system, the increased intensity of turbulence manifesting by larger Reynolds stress and thinner boundary layer, and the more significant contribution of the pressure force on the surface of rough wall to the torque., Comment: 17 pages,11 figures

The presence of a dispersed phase can significantly modulate the drag in turbulent systems. We derived a conserved quantity that characterizes the radial transport of azimuthal momentum in the fluid-fluid two-phase Taylor-Couette turbulence. This quantity consists of contributions from advection, diffusion, and two-phase interface, which are closely related to density, viscosity, and interfacial tension, respectively. We found that the presence of the two-phase interface consistently produces a positive contribution to the momentum transport and leads to drag enhancement, while decreasing the density and viscosity ratios of the dispersed phase to the continuous phase reduces the contribution of local advection and diffusion terms to the momentum transport, respectively, resulting in drag reduction. Therefore, we concluded that the decreased density ratio and the decreased viscosity ratio work together to compete with the presence of two-phase interface for achieving drag modulation in fluid-fluid two-phase turbulence.

Catastrophic phase inversion, the sudden breakdown of a dense emulsion, occurs when the dispersed majority phase irreversibly exchanges role with the continuous minority phase. This common process has been extensively studied over the past decades and yet its fundamental physical mechanism has remained largely unexplored. Here we experimentally and numerically study the dynamics of catastrophic phase inversion as it occurs when the volume fraction of the dispersed phase exceeds a critical volume fraction (typically around 92% in experiments). Our data accurately quantify the abrupt change of both the global torque and average droplet size at approaching and across the phase inversion point, exhibiting strong hysteresis. Most importantly, we reveal that the fluctuations in the global torque diverge as a power-law while approaching the critical volume-fraction and we connect their growth to the formation of highly heterogeneous spatial droplet structures. The present finding, unveiling the tight connection between fluctuations in dynamic heterogeneity and the critical divergence of torque fluctuation, paves the way to a quantitative description of catastrophic phase inversion as an out-of-equilibrium critical-like phenomena.

Wind-wave interactions impose wind forcing on wave surface and wave effects on turbulent wind structures, which essentially influences the wind-wave loading on structures. Existing research treats the wind and wave loading separately and ignores their interactions. The present study aims to characterize the turbulent airflow over wave surfaces and analyze the coupled wind-wave loading on offshore wind turbines. A high-fidelity two-phase model is developed to simulate highly turbulent wind-wave fields based on the open-source program OpenFOAM. A numerical case study is conducted to simulate extreme wind-wave conditions, where coupled wind-wave fields are applied. The intensity of wave induced turbulences and the height of wave influenced regions are affected by the wind velocity and wave heights. Higher wind velocities induce greater turbulence, which can be increased by over 100%. Then the combined wind-wave loading on offshore wind turbines is simulated under operational and extreme conditions. Under operational conditions, the wind-wave coupling effect on the combined loading is minimal. However, under extreme conditions, the coupled wind-wave fields lead to an increase in the average aerodynamic loading and a significant amplification of the fluctuation in the aerodynamic loading. Specifically, the maximum bending moment Fy at both the tower bottom and the monopile bottom experiences an increase of around 6%. The standard deviation of the shear force at the tower bottom increases by up to 45%. Also, the standard deviation of the bending moment at the tower bottom increases by approximately 27%. This study reveals the importance of considering the wind-wave coupling effect under extreme conditions, which provides valuable insights into the planning and design of offshore wind turbines., Comment: arXiv admin note: substantial text overlap with arXiv:2206.02085

Recent extreme tropical cyclones caused extensive damages to critical civil infrastructure globally. To better capture the unique hurricane wind characteristics, a large eddy simulation (LES) Hurricane Boundary Layer (HBL) model is developed by considering the variation of meso-scale kinematic and thermodynamic conditions. An asymmetric model is adopted to obtain the gradient wind velocity using the National Hurricane Center data. The meso-scale thermal conditions are obtained by extracting the hourly air temperature and relative humidity profiles from generated proxy soundings. Measurements recorded at the Aransas County airport during Hurricane Harvey and that at the City of Naples during Irma are used to validate the developed LES model. Research results show that the simulated 10-minute average wind speed and direction are consistent with the observations. The developed model can well predict the high wind turbulences, which are around 20% in Hurricane Harvey and 26% in Hurricane Irma. The 3-s gust wind speeds reach 62.4 m/s at 10-m elevation during Hurricane Harvey and 53.5 m/s at 15-m elevation during Hurricane Irma, close to the field observed data of 61.3 m/s and 54.2 m/s, respectively. The simulated 3-s gust factors are close to the observation except at some moments with significant variations because of the poorly understood physical phenomena. The simulated wind spectrum in longitudinal and lateral directions agrees well with the observed results. In summary, the developed LES-based HBL model can capture the main characteristics of hurricane structure and turbulence characteristics and is applicable for modeling civil infrastructure exposed to hurricanes.

Improving the critical temperature (TC), critical magnetic field (HC), and critical current (JC) of superconducting materials has always been one of the most significant challenges in the field of superconductivity, but progress has been slow over the years. Based on the concept of injecting energy to enhance electron pairing states, in this study, we have employed a solid-state sintering method to fabricate a series of smart meta-superconductors (SMSCs) consisting of p-n junction nanostructures with a wavelength of 550 nm, doped within an MgB2 matrix. Experimental results demonstrate that compared to pure MgB2 samples, the critical transition temperature (TC) has increased by 1.2 K, the critical current (JC) has increased by 52.8%, and the Meissner effect (HC) shows significant improvement in its diamagnetic properties. This phenomenon of enhanced superconducting performance can be explained by the coupling between superconducting electrons and evanescent waves.

Hypothesis: Surfactant-free microemulsion (SFME), an emerging phenomenology that occurs in the monophasic zone of a broad category of ternary mixtures 'hydrophobe/hydrotrope/water', has attracted extensive interests due to their unique physicochemical properties. The potential of this kind of ternary fluid for solubilization and drug delivery make them promising candidates in many industrial scenarios. Experiments: Here the thermodynamic behavior of these multiscale nanodomains formed in the ternary trans-anethol/ethanol/water system over a wide range of temperatures is explored. The macroscopic physical properties of the ternary solutions are characterized, with revealing the temperature dependence of refractive index and dynamic viscosity. Findings: With increasing temperature, the ternary system shows extended areas in the monophasic zone. We demonstrate that the phase behavior and the multiscale nanodomains formed in the monophasic zone can be precisely and reversibly tuned by altering the temperature. Increasing temperature can destroy the stability of the multiscale nanodomains in equilibrium, with an exponential decay in the scattering light intensity. Nevertheless, molecular-scale aggregates and mesoscopic droplets exhibit significantly different response behaviors to temperature stimuli. The temperature-sensitive nature of the ternary SFME system provides a crucial step forward exploring and industrializing its stability.

In this Letter, we systematically investigate the Leidenfrost temperature for hot solid substrates with various thermal diffusivities and surface roughnesses. Based on the experimental results, we build a phenomenological model that considers the thermal diffusivity of a solid substrate and derive a relationship between the surface roughness and the resulting vapor film thickness. The generality of this model is supported by experimental data for different liquids and solid substrates. Our model thus allows for a theoretical prediction of the Leidenfrost temperature and develops a comprehensive understanding of the Leidenfrost effect., Comment: 6 pages, 5 figures

Hypothesis: Multi-component droplet evaporation has received significant attention in recent years due to the broad range of applications such as material science, environmental monitoring, and pharmaceuticals. The selective evaporation induced by the different physicochemical properties of components is expected to influence the concentration distributions and the separation of mixtures, thereby leading to rich interfacial phenomena and phase interactions. Experiments: A ternary mixture system containing hexadecane, ethanol, and diethyl ether is investigated in this study. The diethyl ether exhibits both surfactant-like and co-solvent properties. Systematic experiments were performed using acoustic levitation technique to achieve a contact-less evaporation condition. The evaporation dynamics and temperature information are acquired in the experiments, using high-speed photography and infrared thermography technologies. Findings: Three distinct stages, namely, "Ouzo state", "Janus state", and "Encapsulating state", are identified for the evaporating ternary droplet in acoustic levitation. A self-sustaining periodic freezing & melting evaporation mode is reported. A theoretical model is developed to characterize the multi-stage evaporating behaviors. We demonstrate the capability to tune the evaporating behaviors by varying the initial droplet composition. This work provides a deeper understanding of the interfacial dynamics and phase transitions involved in multi-component droplets and proposes novel strategies for the design and control of droplet-based systems

Turbulent plane-Couette flow suspended with finite-size spheroidal particles is studied using fully particle-resolved direct numerical simulations. The effects of particle aspect ratio on turbulent arguments and particle statistics are explored, leading to the same conclusions as the previous experimental findings \citep{wang2022finite}. By performing stress analysis, we find that the presence of particles introduces extra stresses to the system and accounts for the global drag increases. The particle-laden flow cases exhibit spectra that are consistent with the scaling of $k^{-5/3}$ and $k^{-3}$ in the large and small scales, respectively. While the $k^{-3}$ scaling observed in the particle-laden flow is reminiscent of bubbly flow, an examination of the particle Reynolds number suggests that the mechanism responsible may not be attributable to the pseudo-turbulence induced by particles as in the case of bubbles. In the view of particle statistics, we observe that spherical and non-spherical particles preferentially cluster in the near-wall and the bulk region, respectively, and that the orientations of non-spherical particles are affected by their aspect ratios, especially in the near-wall region. The present numerical results, combined with previous experimental findings in \cite{wang2022finite}, provide in-depth information on both the fluid and the particle phase, contributing to a better understanding of particle suspension in shear flows., Comment: 22 pages, 15 figures

We perform a two-dimensional numerical study on the thermal effect of porous media on global heat transport and flow structure in Rayleigh-B\'enard (RB) convection, focusing on the role of thermal conductivity $\lambda$ of porous media, which ranges from $0.1$ to $50$ relative to the fluid. The simulation is carried out in a square RB cell with the Rayleigh number $Ra$ ranging from $10^7$ to $10^9$ and the Prandtl number $Pr$ fixed at $4.3$. The porosity of the system is fixed at $\phi=0.812$, with the porous media modeled by a set of randomly displayed circular obstacles. For a fixed $Ra$, the increase of conductivity shows a small effect on the total heat transfer, slightly depressing the Nusselt number. The limited influence comes from the small number of obstacles contacting with thermal plumes in the system as well as the counteraction of the increased plume area and the depressed plume strength. The study shows that the global heat transfer is insensitive to the conduction effect of separated porous media in the bulk region, which may have implications for industrial designs.

We investigate experimentally and theoretically the interactions between a cavitation bubble and a hemispherical pendant oil droplet immersed in water. In experiments, the cavitation bubble is generated by a focused laser pulse right below the pendant droplet with well-controlled bubble-wall distances and bubble-droplet size ratios. By high-speed imaging, four typical interactions are observed, namely, oil droplet rupture, water droplet entrapment, oil droplet large deformation, and oil droplet mild deformation. The bubble jetting at the end of collapse and the migration of the bubble centroid are particularly different in each bubble-droplet interaction. We propose theoretical models based on the method of images for calculating the Kelvin impulse and the anisotropy parameter which quantitatively reflects the migration of the bubble centroid at the end of the collapse. Finally, we explain that a combination of the Weber number and the anisotropy parameter determines the regimes of the bubble-droplet interactions., Comment: 26 pages, 14 figures

Sea ice is crucial in many natural processes and human activities. Understanding the dynamical couplings between the inception, growth and equilibrium of sea ice and the rich fluid mechanical processes occurring at its interface and interior is of relevance in many domains ranging from geophysics to marine engineering. Here we experimentally investigate the complete freezing process of water with dissolved salt in a standard natural convection system, i.e., the prototypical Rayleigh-B\'enard cell. Due to the presence of a mushy phase, the studied system is considerably more complex than the freezing of freshwater in the same conditions (Wang et al. 2021c). We measure the ice thickness and porosity at the dynamical equilibrium state for different initial salinities of the solution and temperature gaps across the cell. These observables are non-trivially related to the controlling parameters of the system as they depend on the heat transport mode across the cell. We identify in the experiments 5 out of the 6 possible modes of heat transport. We highlight the occurrence of brine convection through the mushy ice and of penetrative convection in stably stratified liquid underlying the ice. A one-dimensional multi-layer heat flux model built on the known scaling relations of global heat transport in natural convection systems in liquids and porous media is proposed. It allows, given the measured porosity of the ice, to predict the corresponding ice thickness, in a unified framework., Comment: 20 pages, 6 figures

Photolithography assisted chemo-mechanical etching (PLACE), a technique specifically developed for fabricating highquality large-scale photonic integrated circuits (PICs) on thin-film lithium niobate (TFLN), has enabled fabrication of a series of building blocks of PICs ranging from high-quality (high-Q) microresonators and low-loss waveguides to electrooptically (EO) tunable lasers and waveguide amplifiers. Aiming at high-throughput manufacturing of the PIC devices and systems, we have developed an ultra-high-speed high-resolution laser lithography fabrication system employing a high repetition rate femtosecond laser and a high-speed polygon laser scanner, by which a lithography fabrication efficiency of 4.8 cm2/h has been achieved at a spatial resolution of 200 nm. We demonstrate wafer-scale fabrication of TFLN-based photonic structures, optical phase masks as well as color printing

We investigate the coupling effect of buoyancy and shear based on an annular centrifugal Rayleigh-B\'enard convection (ACRBC) system in which two cylinders rotate with an angular velocity difference. Direct numerical simulations are performed in a Rayleigh number range 10^6 \le Ra \le 10^8, at fixed Prandtl number Pr=4.3, inversed Rossby number Ro^{-1}=20 and radius ratio \eta=0.5. The shear, represented by the non-dimensional rotational speed difference \Omega, varies from 0 to 10, corresponding to an ACRBC without shear and a radially heated Taylor-Couette flow with only the inner cylinder rotating, respectively. A stable regime is found in the middle part of the interval of \Omega, and divides the whole parameter space into three regimes: buoyancy-dominated regime, stable regime, and shear-dominated regime. Clear boundaries between the regimes are given by linear stability analysis. In the buoyancy-dominated regime, the flow is a quasi-two-dimensional flow on the r\varphi plane; as shear increases, both the growth rate of instability and the heat transfer is depressed. In the shear-dominated regime, the flow is mainly on the rz plane, and the heat transfer in this regime is greatly enhanced. The study shows shear can stabilize buoyancy-driven convection and reveals the complex coupling mechanism of shear and buoyancy, which may have implications for fundamental studies and industrial designs.

Deep neural networks (DNNs) have great potential to solve many real-world problems, but they usually require an extensive amount of computation and memory. It is of great difficulty to deploy a large DNN model to a single resource-limited device with small memory capacity. Distributed computing is a common approach to reduce single-node memory consumption and to accelerate the inference of DNN models. In this paper, we explore the "within-layer model parallelism", which distributes the inference of each layer into multiple nodes. In this way, the memory requirement can be distributed to many nodes, making it possible to use several edge devices to infer a large DNN model. Due to the dependency within each layer, data communications between nodes during this parallel inference can be a bottleneck when the communication bandwidth is limited. We propose a framework to train DNN models for Distributed Inference with Sparse Communications (DISCO). We convert the problem of selecting which subset of data to transmit between nodes into a model optimization problem, and derive models with both computation and communication reduction when each layer is inferred on multiple nodes. We show the benefit of the DISCO framework on a variety of CV tasks such as image classification, object detection, semantic segmentation, and image super resolution. The corresponding models include important DNN building blocks such as convolutions and transformers. For example, each layer of a ResNet-50 model can be distributively inferred across two nodes with five times less data communications, almost half overall computations and half memory requirement for a single node, and achieve comparable accuracy to the original ResNet-50 model. This also results in 4.7 times overall inference speedup.

Telescope arrays are receiving increasing attention due to their promise of higher resource utilization, greater sky survey area, and higher frequency of full space-time monitoring than single telescopes. Compared with the ordinary coordinated operation of several telescopes, the new astronomical observation mode has an order of magnitude difference in the number of telescopes. It requires efficient coordinated observation by large-domain telescopes distributed at different sites. Coherent modeling of various abstract environmental constraints is essential for responding to multiple complex science goals. Also, due to competing science priorities and field visibility, how the telescope arrays are scheduled for observations can significantly affect observation efficiency. This paper proposes a multilevel scheduling model oriented toward the problem of telescope-array scheduling for time-domain surveys. A flexible framework is developed with basic functionality encapsulated in software components implemented on hierarchical architectures. An optimization metric is proposed to self-consistently weight contributions from time-varying observation conditions to maintain uniform coverage and efficient time utilization from a global perspective. The performance of the scheduler is evaluated through simulated instances. The experimental results show that our scheduling framework performs correctly and provides acceptable solutions considering the percentage of time allocation efficiency and sky coverage uniformity in a feasible amount of time. Using a generic version of the telescope-array scheduling framework, we also demonstrate its scalability and its potential to be applied to other astronomical applications., Comment: Accepted for publication in the Astronomical Journal

Cross-matching operation, which is to find corresponding data for the same celestial object or region from multiple catalogues,is indispensable to astronomical data analysis and research. Due to the large amount of astronomical catalogues generated by the ongoing and next-generation large-scale sky surveys, the time complexity of the cross-matching is increasing dramatically. Heterogeneous computing environments provide a theoretical possibility to accelerate the cross-matching, but the performance advantages of heterogeneous computing resources have not been fully utilized. To meet the challenge of cross-matching for substantial increasing amount of astronomical observation data, this paper proposes Heterogeneous-computing-enabled Large Catalogue Cross-matcher (HLC2), a high-performance cross-matching framework based on spherical position deviation on CPU-GPU heterogeneous computing platforms. It supports scalable and flexible cross-matching and can be directly applied to the fusion of large astronomical cataloguesfrom survey missions and astronomical data centres. A performance estimation model is proposed to locate the performance bottlenecks and guide the optimizations. A two-level partitioning strategy is designed to generate an optimized data placement according to the positions of celestial objects to increase throughput. To make HLC2 a more adaptive solution, the architecture-aware task splitting, thread parallelization, and concurrent scheduling strategies are designed and integrated. Moreover, a novel quad-direction strategy is proposed for the boundary problem to effectively balance performance and completeness. We have experimentally evaluated HLC2 using public released catalogue data. Experiments demonstrate that HLC2 scales well on different sizes of catalogues and the cross-matching speed is significantly improved compared to the state-of-the-art cross-matchers., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

In modern solid-state drives (SSDs), the indexing of flash pages is a critical component in their storage controllers. It not only affects the data access performance, but also determines the efficiency of the precious in-device DRAM resource. A variety of address mapping schemes and optimization techniques have been proposed. However, most of them were developed with human-driven heuristics. They cannot automatically capture diverse data access patterns at runtime in SSD controllers, which leaves a large room for improvement. In this paper, we present a learning-based flash translation layer (FTL), named LeaFTL, which learns the address mapping to tolerate dynamic data access patterns via linear regression at runtime. By grouping a large set of mapping entries into a learned segment, it significantly reduces the memory footprint of the address mapping table, which further benefits the data caching in SSD controllers. LeaFTL also employs various optimization techniques, including out-of-band metadata verification to tolerate mispredictions, optimized flash allocation, and dynamic compaction of learned index segments. We implement LeaFTL with an SSD simulator and evaluate it with various storage workloads. LeaFTL saves the memory consumption of the mapping table by 2.9x on average and improves the storage performance by 1.4x on average, in comparison with state-of-the-art FTL schemes., Comment: This paper is accepted at the 28th Conference on Architectural Support for Programming Languages and Operating Systems (ASPLOS 2023)

The dispersed phase in turbulence can vary from almost inviscid fluid to highly viscous fluid. By changing the viscosity of the dispersed droplet phase, we experimentally investigate how the deformability of dispersed droplets affects the global transport quantity of the turbulent emulsion. Different kinds of silicone oil are employed to result in the viscosity ratio, $\zeta$, ranging from $0.53$ to $8.02$. The droplet volume fraction, $\phi$, is varied from 0\% to 10\% with a spacing of 2\%. The global transport quantity, quantified by the normalized friction coefficient $c_{f,\phi}/c_{f,\phi=0}$, shows a weak dependence on the turbulent intensity due to the vanishing finite-size effect of the droplets. The interesting fact is that, with increasing $\zeta$, the $c_{f,\phi}/c_{f,\phi=0}$ first increases and then saturates to a plateau value which is similar to that of the rigid particle suspension. By performing image analysis, this drag modification is interpreted from the aspect of droplet deformability, which is responsible for the breakup and coalescence effect of the droplets. The statistics of the droplet size distribution show that, with increasing $\zeta$, the stabilizing effect induced by interfacial tension comes to be substantial and the pure inertial breakup process becomes dominant. The measurement of the droplet distribution along the radial direction of the system shows a bulk-clustering effect, which can be attributed to the non-negligible coalescence effect of the droplet. It is found that the droplet coalescence effect could be suppressed as the $\zeta$ increases, thereby affecting the contribution of interfacial tension to the total stress, and accounting for the observed emulsion rheology., Comment: 17 pages, 8 figures

The decay constants of the low lying S-wave $B_c$ mesons, i.e. $B_c(nS)$ and $B^*_c(nS)$ with $n\leq 3$, are calculated in the nonrelativistic quark model. The running coupling of the strong interaction is taken into account, and the uncertainties due to varying parameters and losing Lorentz covariance are considered carefully. As a byproduct, the decay constants of the low lying S-wave charmonium and bottomium states are given in the appendixes., Comment: 9 pages, 2 figures

Emulsions are common in many natural and industrial settings. Recently, much attention has been put on understanding the dynamics of turbulent emulsions. This paper reviews some recent studies of emulsions in turbulent Taylor-Couette flow, mainly focusing on the statistics of the dispersed phase and the global momentum transport of the system. We first study the size distribution and the breakup mechanism of the dispersed droplets for turbulent emulsions with a low volume-fraction (dilute) of the dispersed phase. For systems with a high volume-fraction of the dispersed phase (dense), we address the detailed response of the global transport (effective viscosity) of the turbulent emulsion and its connection to the droplet statistics. Finally, we will discuss catastrophic phase inversions, which can happen when the volume fraction of the dispersed phase exceeds a critical value during dynamic emulsification. We end the manuscript with a summary and an outlook including some open questions for future research. This article is part of the theme issue `Taylor-Couette and Related Flows on the Centennial of Taylor's Seminal Philosophical Transactions Paper'., Comment: 15 pages, 4 figures

In this letter, the merging-vs-bouncing response of obliquely-oriented colliding jets under elevated and reduced gaseous environment pressures was experimentally examined. Experiments with water and n-tetradecane confirmed that the collision outcome transitions from merging to bouncing, and then to merging again, when the impact velocity was increased. This behavior which was previously reported for atmospheric pressure, has now also been observed at elevated and reduced pressures. New results also show that there exists a critical pressure (0.9 bar for tetradecane and 5 bar for water) below which increasing pressure promotes bouncing (expands the bouncing regime), while beyond this, merging is promoted (reduces the bouncing regime) instead. This leads to a non-monotonic influence of pressure on the non-coalescence outcomes of collisional jets, which was not previously reported. The study provides evidence of new behaviors in colliding jets at reduced and elevated pressures, which differs from well-studied droplet-droplet collisions., Comment: Corrected the abstract, results unchanged, a few minor edits based on reviews. This is the published version in Phy. Fluids

Water droplet freezing is a common phenomenon in our daily life. In both natural scenarios and industrial production, different surface inclinations bring distinctive deformation and freezing dynamics to frozen droplets. We explore the freezing of pendent and sessile droplets at different Bond number regimes. The effect of gravity on the droplet freezing process is analyzed by considering droplet morphology, freezing front dynamics, and freezing time. It is found that gravity can significantly influence droplet freezing processes via shaping the initial droplet, resulting in the flattening or elongation of pendent and sessile droplets, respectively. We show that the droplet initial geometry is the most important parameter and it completely controls the droplet freezing. Despite the significant difference in the initial droplet shape several remarkable similarities have been found for pendent and sessile droplets at small and large Bond numbers. The final height of a frozen droplet is found to be linearly proportional to its initial height. The time evolution of the ice-liquid-air contact line is found to reproduce the power-law $t^{0.5}$, but noticeably faster than the Stefan 1-D icing front propagation. As a consequence, the time to freeze a droplet is faster than predicted by the Stefan model and it is found to be dependent on the initial droplet height and base radius through a simple power-law.

We introduce MABe22, a large-scale, multi-agent video and trajectory benchmark to assess the quality of learned behavior representations. This dataset is collected from a variety of biology experiments, and includes triplets of interacting mice (4.7 million frames video+pose tracking data, 10 million frames pose only), symbiotic beetle-ant interactions (10 million frames video data), and groups of interacting flies (4.4 million frames of pose tracking data). Accompanying these data, we introduce a panel of real-life downstream analysis tasks to assess the quality of learned representations by evaluating how well they preserve information about the experimental conditions (e.g. strain, time of day, optogenetic stimulation) and animal behavior. We test multiple state-of-the-art self-supervised video and trajectory representation learning methods to demonstrate the use of our benchmark, revealing that methods developed using human action datasets do not fully translate to animal datasets. We hope that our benchmark and dataset encourage a broader exploration of behavior representation learning methods across species and settings., Comment: To appear in ICML 2023, Project website: https://sites.google.com/view/computational-behavior/our-datasets/mabe2022-dataset

We report on the modification of the spectrum of a passive scalar inside a turbulent flow by the injection of large bubbles. While the spectral modification through bubbles is well known and well analyzed for the velocity fluctuations, little is known on how bubbles change the fluctuations of an approximately passive scalar, in our case temperature. Here we uncover the thermal spectral scaling behavior of a turbulent multiphase thermal mixing layer. The development of a $-3$ spectral scaling is triggered. By injecting large bubbles ($\text{Re}_{\text{bub}} = \mathcal{O}(10^2)$) with gas volume fractions $\alpha$ up to 5\%. For these bubbly flows, the $-5/3$ scaling is still observed at intermediate frequencies for low $\alpha$ but becomes less pronounced when $\alpha$ further increases and it is followed by a steeper $-3$ slope for larger frequencies. This $-3$ scaling range extends with increasing gas volume fraction. The $-3$ scaling exponent coincides with the typical energy spectral scaling for the velocity fluctuations in high Reynolds number bubbly flows. We identify the frequency scale of the transition from the $-5/3$ scaling to the $-3$ scaling and show how it depends on the gas volume fraction.

Wind-wave interaction involves wind forcing on wave surface and wave effects on the turbulent wind structures, which essentially influences the wind and wave loading on structures. Existing research on wind-wave interaction modeling ignores the inherent strong turbulences of wind. The present study aims to characterize the turbulent airflow over wave surfaces and wave dynamics under wind driving force. A high-fidelity two-phase model is developed to simulate highly turbulent wind-wave fields. Instead of using uniform wind, inherent wind turbulences are prescribed at the inlet boundary using the turbulent spot method. The developed model is validated by comparing the simulated wind-wave flow characteristics with experimental data. With the validated model, a numerical case study is conducted on a 10^2 m scale under extreme wind and wave conditions. The result shows that when inherent wind turbulences are considered, the resultant turbulence is strengthened and is the summation of the inherent turbulence and the wave-induced turbulence. In addition, the wave coherent velocities and shelter effect are enhanced because of the presence of wind inherent turbulence. The regions of intense turbulence depend on the relative speed between wind velocity and wave phase speed. Higher wind velocities induce greater turbulence intensities, which can be increased by up to 17%. The different relative speed between wind and wave can induce opposite positive-negative patterns of wave coherent velocities. The wave-coherent velocity is approximately proportional to the wind velocity, and the influenced region mainly depends on the wave heights.

Energy carrier evolution is crucial for material performance. Ultrafast microscopy has been widely applied to visualize the spatiotemporal evolution of energy carriers. However, direct imaging of small amounts of energy carriers on nanoscale remains difficult due to extremely weak transient signals. Here we present a method for ultrasensitive and high-throughput imaging of energy carrier evolution in space and time. This method combines femtosecond pump-probe techniques with interferometric scattering microscopy (iSCAT), named Femto-iSCAT. The interferometric principle and unique spatially-modulated contrast enhancement increase the transient image contrast by >2 orders of magnitude and enable the exploration of new science. We address three important and challenging problems: transport of different energy carriers at various interfaces, heterogeneous hot electron distribution and relaxation in single plasmonic resonators, and distinct structure-dependent edge state dynamics of carriers and excitons in optoelectronic semiconductors. Femto-iSCAT holds great potential as a universal tool for ultrasensitive imaging of energy carrier evolution in space and time.

We study the dynamics of gyrotactic swimmers in turbulence, whose orientation is governed by gravitational torque and local fluid velocity gradient. The gyrotaxis strength is measured by the ratio of the Kolmogorov time scale to the reorientation time scale due to gravity, and a large value of this ratio means the gyrotaxis is strong. By means of direct numerical simulations, we investigate the effects of swimming velocity and gyrotactic stability on spatial accumulation and alignment. Three-dimensional Vorono{\"\i} analysis is used to study the spatial distribution and time evolution of the particle concentration. We study spatial distribution by examing the overall preferential sampling and where clusters and voids (subsets of particles that have small and large Vorono{\"\i} volumes respectively) form. Compared with the ensemble particles, the preferential sampling of clusters and voids is found to be more pronounced. The clustering of fast swimmers lasts much longer than slower swimmers when the gyrotaxis is strong and intermediate, but an opposite trend emerges when the gyrotaxis is weak. In addition, we study the preferential alignment with the Lagrangian stretching direction, with which passive slender rods have been known to align. We show that the Lagrangian alignment is reduced by the swimming velocity when the gyrotaxis is weak, while the Lagrangian alignment is enhanced for the regime in which gyrotaxis is strong.

We study the translational and rotational dynamics of neutrally-buoyant finite-size spheroids in hydrodynamic turbulence by means of fully resolved numerical simulations. We examine axisymmetric shapes, from oblate to prolate, and the particle volume dependences. We show that the accelerations and rotations experienced by non-spherical inertial-scale particles result from volume filtered fluid forces and torques, similar to spherical particles. However, the particle orientations carry signatures of preferential alignments with the surrounding flow structures, which is reflected in distinct axial and lateral fluctuations for accelerations and rotation rates. The randomization of orientations does not occur even for particles with volume equivalent diameter size in the inertial range, here up to 60 $\eta$ at $Re_{\lambda}=120$. Additionally, we demonstrate that the role of fluid boundary layers around the particles cannot be neglected to reach a quantitative understanding of particle statistical dynamics, as they affect the intensities of angular velocities, and the relative importance of tumbling with respect to spinning rotations. This study brings to the fore the importance of inertial-scale flow structures in homogeneous and isotropic turbulence and their impacts on the transport of neutrally-buoyant bodies with size in the inertial range., Comment: 21 pages, 11 figures

By varying the oil volume fraction, the microscopic droplet size and the macroscopic rheology of emulsions are investigated in a Taylor-Couette (TC) turbulent shear flow. Although here oil and water in the emulsions have almost the same physical properties (density and viscosity), unexpectedly, we find that oil-in-water (O/W) and water-in-oil (W/O) emulsions have very distinct hydrodynamic behaviors, i.e., the system is clearly asymmetric. By looking at the micro-scales, the average droplet diameter hardly changes with the oil volume fraction neither for O/W nor for W/O. However, for W/O it is about 50% larger than that of O/W. At the macro-scales, the effective viscosity of O/W is higher when compared to that of W/O. These asymmetric behaviors can be traced back to the presence of surface-active contaminants in the system. By introducing an oil-soluble surfactant at high concentration, remarkably, we recover the symmetry (droplet size and effective viscosity) between O/W and W/O emulsions. Based on this, we suggest a possible mechanism responsible for the initial asymmetry. Next, we discuss what sets the droplet size in turbulent emulsions. We uncover a -6/5 scaling dependence of the droplet size on the Reynolds number of the flow. Combining the scaling dependence and the droplet Weber number, we conclude that the droplet fragmentation, which determines the droplet size, occurs within the boundary layer and is controlled by the dynamic pressure caused by the gradient of the mean flow, as proposed by Levich (1962), instead of the dynamic pressure due to turbulent fluctuations, as proposed by Kolmogorov (1949). The present findings provide an understanding of both the microscopic droplet formation and the macroscopic rheological behaviors in dynamic emulsification, and connects them.

We experimentally investigate the drag modification by neutrally buoyant finite-size particles with various aspect ratios in a Taylor-Couette (TC) turbulent flow. The current Reynolds number, $Re$, ranges from $6.5\times10^3$ to $2.6\times10^4$, and the particle volume fraction, $\Phi$, is up to $10\%$. Particles with three kinds of aspect ratio, $\lambda$, are used: $\lambda=1/3$ (oblate), $\lambda=1$ (spherical) and $\lambda=3$ (prolate). Unlike the case of bubbly TC flow, we find that the suspended finite-size particles increase the drag of the TC system regardless of their aspect ratios. In addition, the normalized friction coefficient, $c_{f,\Phi}/c_{f,\Phi=0}$, decreases with increasing $Re$, the reason could be that in the current low volume fractions the turbulent stress becomes dominant at higher $Re$. As $Re$ increases, the particles distribute more evenly in the entire system, which results from both the greater turbulence intensity and the more pronounced finite-size effects of the particles at higher $Re$. Moreover, it is found that the variation of the particle aspect ratios leads to different particle collective effects. The suspended spherical particles, which tend to cluster near the walls and form a particle layer, significantly affect the boundary layer and result in maximum drag modification. The minimal drag modification is found in the oblate case, where the particles preferentially cluster in the bulk region, and thus the particle layer is absent. Based on the optical measurement results, it can be concluded that, in the low volume fraction ranges ($\Phi=0.5\%$ and $\Phi = 2\%$ here), the larger drag modification is connected to the near-wall particle clustering. The present findings suggest that the particle shape plays a significant role in drag modification.