Search

Homoclinic and heteroclinic connections can form cycles and networks in phase space, which organize the global phenomena in dynamical systems. On the one hand, stability notions for (omni)cycles give insight into how many initial conditions approach the network along a single given (omni)cycle. On the other hand, the term switching is used to describe situations where there are trajectories that follow any possible sequence of heteroclinic connections along the network. Here we give a notion of asymptotic stability for general sequences along a network of homoclinic and heteroclinic connections. We show that there cannot be uncountably many aperiodic sequences that attract a set with nontrivial measure. Finally, we discuss examples where one may or may not expect aperiodic convergence towards a network and conclude with some open questions.

In so-called Lagrangian melting problems, a solid immersed in a fluid medium is free to rotate and translate in tandem with its phase-change from solid to liquid. Such configurations may be classified as a fluid-solid interaction (FSI) problem coupled to phase-change. Our present work proposes a numerical method capable of simulating these Lagrangian melting problems and adopts a front-tracking immersed-boundary (IB) method. We use the moving least squares IB framework, a well-established method for simulating a diverse range of FSI problems and extend this framework to accommodate melting by additionally imposing the Stefan condition at the interface. In the spirit of canonical front-tracking methods, the immersed solid is represented by a discrete triangulated mesh which is separate from the Eulerian mesh in which the governing flow equations are solved. A known requirement for these methods is the need for comparable Eulerian and Lagrangian grid spacings to stabilise interpolation and spreading operations between the two grids. For a melting object, this requirement is inevitably violated unless interventional remeshing is introduced. Our work therefore presents a novel dynamic remeshing procedure to overcome this. The remeshing is based on a gradual coarsening of the triangulated Lagrangian mesh and amounts to a negligible computational burden per timestep owing to the incremental and local nature of its operations, making it a scalable approach. Moreover, the coarsening is coupled to a volume-conserving smoothing procedure detailed by Kuprat et al. (2001), ensuring a zero net volume change in the remeshing step to machine precision. This added feature makes our present method highly specialised to the study of melting problems, where precise measurements of the melting solid's volume is often the primary predictive quantity of interest.

Sea ice is a key element of the global Earth system, with a major impact on global climate and regional weather. Unfortunately, accurate sea ice modeling is challenging due to the diversity and complexity of underlying physics happening there, and a relative lack of ground truth observations. This is especially true for the Marginal Ice Zone (MIZ), which is the area where sea ice is affected by incoming ocean waves. Waves contribute to making the area dynamic, and due to the low survival time of the buoys deployed there, the MIZ is challenging to monitor. In 2022-2023, we released 79 OpenMetBuoys (OMBs) around Svalbard, both in the MIZ and the ocean immediately outside of it. OMBs are affordable enough to be deployed in large number, and gather information about drift (GPS position) and waves (1-dimensional elevation spectrum). This provides data focusing on the area around Svalbard with unprecedented spatial and temporal resolution. We expect that this will allow to perform validation and calibration of ice models and remote sensing algorithms.

Electrolytically generated gas bubbles can significantly hamper the overall electrolysis efficiency. Therefore it is crucial to understand their dynamics in order to optimise water electrolyzer systems. Here we demonstrate a distinct transport mechanism where coalescence with microbubbles drives electrolyte droplets, resulting from the fragmentation of the Worthington jet, into the gas phase during hydrogen evolution reaction, both in normal and microgravity environments. This indicates that the H$_2$ bubble is not only composed of hydrogen gas and vapor but also includes electrolyte fractions. Reminiscent of bursting bubbles on a liquid-gas interface, this behavior results in a flow inside the bubble, which is further affected by Marangoni convection at the gas-electrolyte interface, highlighting interface mobility. In the case of electrode-attached bubbles, the sprayed droplets form electrolyte puddles at the bubble-electrode contact area, affecting the dynamics near the three-phase contact line and favoring bubble detachment from the electrode. The results of this work unravel important insights into the physicochemical aspects of electrolytic gas bubbles, integral for optimizing gas-evolving electrochemical systems. Besides, our findings are essential for studying the limits of jet formation and rupture relevant to acid mist formation in electrowinning, generation of sea spray aerosols, impact of droplets on liquid surfaces, etc., Comment: manuscript: 25 pages, 6 figures; SI: 12 pages, 5 figures, 1 table

The dynamics of drop impact on a rigid surface -- omnipresent in nature and technology -- strongly depends on the droplet's velocity, its size, and its material properties. The main characteristics are the droplet's force exerted on the surface and its maximal spreading radius. The crucial question is: How do they depend on the (dimensionless) control parameters, which are the Weber number $We$ (non-dimensionalized kinetic energy) and the Ohnesorge number $Oh$ (dimensionless viscosity)? Here we perform direct numerical simulations over the huge parameter range $1\le We \le 10^3$ and $10^{-3}\le Oh \le 10^2$ and in particular develop a unifying theoretical approach, which is inspired by the Grossmann-Lohse theory for wall-bounded turbulence [J. Fluid Mech. 407, 27 (2000); PRL 86, 3316 (2001)]. The key idea is to split the energy dissipation rate into the different phases of the impact process, in which different physical mechanisms dominate. The theory can consistently and quantitatively account for the $We$ and $Oh$ dependences of the maximal impact force and the maximal spreading diameter over the huge parameter space. It also clarifies why viscous dissipation plays a significant role during impact, even for low-viscosity droplets (low $Oh$), in contrast to what had been assumed in prior theories.

Bubble bursting and subsequent collapse of the open cavity at free surfaces of contaminated liquids can generate aerosol droplets, facilitating pathogen transport. After film rupture, capillary waves focus at the cavity base, potentially generating fast Worthington jets that are responsible for ejecting the droplets away from the source. While extensively studied for Newtonian fluids, the influence of non-Newtonian rheology on this process remains poorly understood. Here, we employ direct numerical simulations to investigate the bubble cavity collapse in viscoelastic media, such as polymeric liquids, examining how their elastic modulus $G$ and their relaxation time $\lambda$ affect jet and droplet formation. We show that the viscoelastic liquids yield Newtonian-like behavior as either parameter $G$ or $\lambda$ approaches zero, while increasing them suppresses jet formation due to elastic resistance to elongational flows. Intriguingly, for some cases with intermediate values of $G$ and $\lambda$, smaller droplets are produced compared to Newtonian fluids, potentially enhancing aerosol dispersal. By mapping the phase space spanned by the elastocapillary number (dimensionless $G$) and the Deborah number (dimensionless $\lambda$), we reveal three distinct flow regimes: (i) jets forming droplets, (ii) jets without droplet formation, and (iii) absence of jet formation. Our results elucidate the mechanisms underlying aerosol suppression versus fine spray formation in polymeric liquids, with implications for pathogen transmission and industrial processes involving viscoelastic fluids.

We offer a new model for the heat transfer and the turbulence intensity in strongly driven Rayleigh-Benard turbulence (the so-called ultimate regime), which in contrast to hitherto models is consistent with the new mathematically exact heat transfer upper bound of Choffrut et al. [J. Differential Equations 260, 3860 (2016)] and thus enables extrapolations of the heat transfer to geo- and astrophysical flows. The model distinguishes between four subregimes of the ultimate regime and well describes the measured heat transfer in various large-Ra experiments. In this new representation, which properly accounts for the Prandtl number dependence, the onset to the ultimate regime is seen in all available large-Ra data sets, though at different Rayleigh numbers, as to be expected for a non-normal-nonlinear instability., Comment: 14 pages, 2 figures

Evaluating the performance of causal discovery algorithms that aim to find causal relationships between time-dependent processes remains a challenging topic. In this paper, we show that certain characteristics of datasets, such as varsortability (Reisach et al. 2021) and $R^2$-sortability (Reisach et al. 2023), also occur in datasets for autocorrelated stationary time series. We illustrate this empirically using four types of data: simulated data based on SVAR models and Erd\H{o}s-R\'enyi graphs, the data used in the 2019 causality-for-climate challenge (Runge et al. 2019), real-world river stream datasets, and real-world data generated by the Causal Chamber of (Gamella et al. 2024). To do this, we adapt var- and $R^2$-sortability to time series data. We also investigate the extent to which the performance of score-based causal discovery methods goes hand in hand with high sortability. Arguably, our most surprising finding is that the investigated real-world datasets exhibit high varsortability and low $R^2$-sortability indicating that scales may carry a significant amount of causal information., Comment: Contribution for the Causal Inference for Time Series Data Workshop at the 40th Conference on Uncertainty in Artificial Intelligence (CI4TS 2024)

We argue for the epistemic and ethical advantages of pluralism in Reinforcement Learning from Human Feedback (RLHF) in the context of Large Language Models (LLM). Drawing on social epistemology and pluralist philosophy of science, we suggest ways in which RHLF can be made more responsive to human needs and how we can address challenges along the way. The paper concludes with an agenda for change, i.e. concrete, actionable steps to improve LLM development.

In piezo acoustic drop-on-demand inkjet printing a single droplet is produced for each piezo driving pulse. This droplet is typically multicomponent, including surfactants to control the spreading and drying of the droplet on the substrate. However, the role of these surfactants on the droplet formation process remains rather elusive. Surfactant concentration gradients may manifest across microsecond-to-second timescales, spanning both the rapid ejection of ink from the nozzle exit and the comparatively slower idling timescale governing the firing of successive droplets. In the present work, we study the influence of surfactants on droplet formation across 6 orders of magnitude in time. To this end, we visualize the microsecond droplet formation process using stroboscopic laser-induced fluorescence microscopy while we vary the nozzle idle time. Our results show that increasing the idle time up to O(1) s affects only the break-up dynamics of the inkjet but not its velocity. By contrast, for idle times $>$ O(1) s, both the break-up dynamics are altered and the velocity of the inkjet increases. We show that the increased velocity results from a decreased surface tension of the ejected droplet, which we extracted from the observed shape oscillations of the jetted droplets in flight. The measured decrease in surface tension is surprising as the $\mu$s timescale of droplet formation is much faster than the typical ms-to-s timescale of surfactant adsorption. By varying the bulk surfactant concentration, we show that the fast decrease in surface tension results from a local surfactant concentration increase to more than 200 times the CMC. Our results suggest that a local high concentration of surfactant allows for surfactant adsorption to the interface of an inkjet at the $\mu$s-to-ms timescale, which is much faster than the typical ms-to-s timescale of surfactant adsorption., Comment: 13 pages, 8 figures

We present an experimental study on detachment characteristics of hydrogen bubbles during electrolysis. Using a transparent (Pt or Ni) electrode enables us to directly observe the bubble contact line and bubble size. Based on these quantities we determine other parameters such as the contact angle and volume through solutions of the Young-Laplace equation. We observe bubbles without ('pinned bubbles') and with ('spreading bubbles') contact line spreading, and find that the latter mode becomes more prevalent if the concentration of HClO4 is greater than or equal to 0.1 M. The departure radius for spreading bubbles is found to drastically exceed the value predicted by the well-known formula of W. Fritz (Physik. Zeitschr. 1935, 36, 379-384) for this case. We show that this is related to the contact line hysteresis, which leads to pinning of the contact line after an initial spreading phase at the receding contact angle. The departure mode is then similar to a pinned bubble and occurs once the contact angle reaches the advancing contact angle of the surface. A prediction for the departure radius based on these findings is found to be consistent with the experimental data., Comment: 12 figures

We present ExpIris, a separation logic framework for the (amortized) expected cost analysis of probabilistic programs. ExpIris is based on Iris, parametric in the language and the cost model, and supports both imperative and functional languages, concurrency, higher-order functions and higher-order state. ExpIris also offers strong support for correctness reasoning, which greatly eases the analysis of programs whose expected cost depends on their high-level behavior. To enable expected cost reasoning in Iris, we build on the expected potential method. The method provides a kind of "currency" that can be used for paying for later operations, and can be distributed over the probabilistic cases whenever there is a probabilistic choice, preserving the expected value due to the linearity of expectations. We demonstrate ExpIris by verifying the expected runtime of a quicksort implementation and the amortized expected runtime of a probabilistic binary counter.

Large Language Models (LLMs) have achieved impressive performance in text understanding and have become an essential tool for building smart assistants. Originally focusing on text, they have been enhanced with multimodal capabilities in recent works that successfully built visual instruction following assistants. As far as the graph modality goes, however, no such assistants have yet been developed. Graph structures are complex in that they represent relation between different features and are permutation invariant. Moreover, representing them in purely textual form does not always lead to good LLM performance even for finetuned models. As a result, there is a need to develop a new method to integrate graphs in LLMs for general graph understanding. This work explores the integration of the graph modality in LLM for general graph instruction following tasks. It aims at producing a deep learning model that enhances an underlying LLM with graph embeddings and trains it to understand them and to produce, given an instruction, an answer grounded in the graph representation. The approach performs significantly better than a graph to text approach and remains consistent even for larger graphs., Comment: This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

Understanding the stability mechanism of surface micro/nanobubbles adhered to gas-evolving electrodes is essential for improving the efficiency of water electrolysis, which is known to be hindered by the bubble coverage on electrodes. Using molecular simulations, the diffusion-controlled evolution of single electrolytic nanobubbles on wettability-patterned nanoelectrodes is investigated. These nanoelectrodes feature hydrophobic islands as preferential nucleation sites and allow the growth of nanobubbles in the pinning mode. In these simulations, a threshold current density distinguishing stable nanobubbles from unstable nanobubbles is found. When the current density remains below the threshold value, nucleated nanobubbles grow to their equilibrium states, maintaining their nanoscopic size. However, for the current density above the threshold value, nanobubbles undergo unlimited growth and can eventually detach due to buoyancy. Increasing the pinning length of nanobubbles increases the degree of nanobubble instability. By connecting the current density with the local gas oversaturation, an extension of the stability theory for surface nanobubbles [Lohse and Zhang, Phys. Rev. E, 2015, 91, 031003(R)] accurately predicts the nanobubble behavior found in molecular simulations, including equilibrium contact angles and the threshold current density. For larger systems that are not accessible to molecular simulations, continuum numerical simulations with the finite difference method combined with the immersed boundary method are performed, again demonstrating good agreement between numerics and theories.

A readout system based on the Timepix4 ASIC, is being developed for photon science. The TEMPUS detector can be operated in two distinct modes: a photon counting mode, which allows for conventional full-frame readout at rates up to 40 kfps; and an event-driven time-stamping mode, which allows excellent time resolution in the nanosecond regime in measurements with moderate X-ray flux. In this paper, we introduce the initial prototype, a single-chip system, and present the first results obtained at PETRA III and ESRF., Comment: 19 pages, 6 figures

This paper addresses the task of tracking Lissajous trajectories in the presence of quantized positioning sensors. To do so, theoretical results on tracking of continuous time periodic signals in the presence of output quantization are provided. With these results in hand, the application to Lissajous tracking is explored. The method proposed relies on the internal model principle and dispenses perfect knowledge of the system equations. Numerical results show that an arbitrary small scanning resolution is achievable despite large sensor quantization intervals.

We use three-dimensional direct numerical simulations of homogeneous isotropic turbulence in a cubic domain to investigate the dynamics of heavy, chiral, finite-size inertial particles and their effects on the flow. Using an immersed-boundary method and a complex collision model, four-way coupled simulations have been performed, and the effects of particle-to-fluid density ratio, turbulence strength and particle volume fraction have been analysed. We find that freely falling particles on the one hand add energy to the turbulent flow but, on the other hand, they also enhance the flow dissipation: depending on the combination of flow parameters, the former or the latter mechanism prevails, thus yielding enhanced or weakened turbulence. Furthermore, particle chirality entails a preferential angular velocity which induces a net vorticity in the fluid phase. As turbulence strengthens, the energy introduced by the falling particles becomes less relevant and stronger velocity fluctuations alter the solid phase dynamics, making the effect of chirality irrelevant for the large-scale features of the flow. Moreover, comparing the time history of collision events for chiral particles and spheres (at the same volume fraction) suggests that the former tend to entangle, in contrast to the latter which rebound impulsively.

Three-dimensional direct numerical simulations (DNS) are used to revisit the experiments on multibubble cavitation performed by Bremond et al. (https://doi.org/10.1063/1.2396922, Phys. Fluids 18, 121505 (2006), https://doi.org/10.1103/PhysRevLett.96.224501, Phys. Rev. Lett. 96, 224501 (2006)). In particular, we aim at understanding the asymmetry observed therein during the expansion and collapse of bubble clusters subjected to a pressure pulse. Our numerical simulations suggest that the asymmetry is due to the force applied by the imposed pressure pulse and it is a consequence of the finite effective speed of sound in the liquid. By comparing our numerical results to the experiments, we found that the effective speed of sound under the experimental conditions was smaller than that of degassed water due to microbubbles in the system which resulted from prior cavitation experiments in the same setup. The estimated values of the effective speed of sound are consistent with those derived from the classical theory of wave propagation in liquids with small amounts of gas. To support this theory, we also present evidence of tiny bubbles remaining in the liquid bulk as a result of the fragmentation of large bubbles during the prior cavitation experiments. Furthermore, we find that this asymmetry also alters the direction of the liquid jet generated during the last stages of bubble collapse.

Turbulent shear flows driven by a combination of a pressure gradient and buoyancy forcing are investigated using direct numerical simulations. Specifically, we consider the setup of a differentially heated vertical channel subject to a Poiseuille-like horizontal pressure gradient. We explore the response of the system to its three control parameters: the Grashof number $Gr$, the Prandtl number $Pr$, and the Reynolds number $Re$ of the pressure-driven flow. From these input parameters, the relative strength of buoyancy driving to the pressure gradient can be quantified by the Richardson number $Ri=Gr/Re^2$. We compare the response of the mixed vertical convection configuration to that of mixed Rayleigh-B\'enard convection and find a nearly identical behaviour, including an increase in wall friction at higher $Gr$ and a drop in the heat flux relative to natural convection for $Ri=O(1)$. This closely matched response is despite vastly different flow structures in the systems. No large-scale organisation is visible in visualisations of mixed vertical convection - an observation that is quantitatively confirmed by spectral analysis. This analysis, combined with a statistical description of the wall heat flux, highlights how moderate shear suppresses the growth of small-scale plumes and reduces the likelihood of extreme events in the local wall heat flux. Vice versa, starting from a pure shear flow, the addition of thermal driving enhances the drag due to the emission of thermal plumes., Comment: 31 pages, 19 figures, submitted to J. Fluid Mech

Ensembles of identical atoms exhibit peculiar collective properties in their interaction with radiation depending on geometry and environment where they are embedded in. A remarkably clean and versatile platform to study collective effects in resonant light scattering are M\"ossbauer nuclei placed in planar x-ray waveguides. Here we conceive and demonstrate experimentally distinct temporal emission characteristics in these systems, ranging from a tunable accelerated exponential decay all the way to a pronounced oscillatory emission pattern, depending on the waveguide geometry and mode of excitation. The observed temporal and spatial emission characteristics of the collectively excited nuclear state in the waveguide -- the nuclear exciton -- are well reproduced by a unified theoretical model. Our findings pave the way for applications ranging from fundamental studies of cooperative emission at hard x-ray frequencies up to new methods of narrowband x-ray control via the engineering of collective radiation patterns.

We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh-B\'enard convection, i.e., in Rayleigh-B\'enard flow with Couette or Poiseuille type shear forcing, by extending the Grossmann & Lohse (2000,2001,2002,2004) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$, the Prandtl number $Pr$, and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that in turbulent Rayleigh-B\'enard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (1932) logarithmic friction law, suggesting some kind of universality between purely shear driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$, $0.5 \leq Pr \leq 5.0$, and $0 \leq Re_S \leq 10^4$., Comment: 22 pages, 9 figures, 2 tables. Submitted to J. Fluid Mech

Over the past decade, there has been a growing interest in the study of multicomponent drops. These drops exhibit unique phenomena, as the interplay between hydrodynamics and the evolving physicochemical properties of the mixture gives rise to distinct and often unregulated behaviors. Of particular interest is the complex dynamic behavior of the drop contact line, which can display self-lubrication effect. The presence of a slipping contact line in self-lubricating multicomponent drops can suppress the coffee-stain effect, conferring valuable technological applications. This review will explain the current understanding of the self-lubrication effect of drops, and cover an analysis of fundamental concepts and recent advances in colloidal assembly. The potential applications of self-lubricating drops across different fields will also be highlighted.

X-ray near field holography has proven to be a powerful 2D and 3D imaging technique with applications ranging from biomedical research to material sciences. To reconstruct meaningful and quantitative images from the measurement intensities, however, it relies on computational phase retrieval which in many cases assumes the phase-shift and attenuation coefficient of the sample to be proportional. Here, we demonstrate an efficient phase retrieval algorithm that does not rely on this homogeneous-object assumption and is a generalization of the well-established contrast-transfer-function (CTF) approach. We then investigate its stability and present an experimental study comparing the proposed algorithm with established methods. The algorithm shows superior reconstruction quality compared to the established CTF-based method at similar computational cost. Our analysis provides a deeper fundamental understanding of the homogeneous object assumption and the proposed algorithm will help improve the image quality for near-field holography in biomedical applications

Effective communication between humans and collaborative robots is essential for seamless Human-Robot Collaboration (HRC). In noisy industrial settings, nonverbal communication, such as gestures, plays a key role in conveying commands and information to robots efficiently. While existing literature has thoroughly examined gesture recognition and robots' responses to these gestures, there is a notable gap in exploring the design of these gestures. The criteria for creating efficient HRC gestures are scattered across numerous studies. This paper surveys the design principles of HRC gestures, as contained in the literature, aiming to consolidate a set of criteria for HRC gesture design. It also examines the methods used for designing and evaluating HRC gestures to highlight research gaps and present directions for future research in this area.

A uniform solidification front undergoes non-trivial deformations when encountering an insoluble dispersed particle in a melt. For solid particles, the overall deformation characteristics are primarily dictated by heat transfer between the particle and the surroundings, remaining unaffected by the rate of approach of the solidification front. In this Letter, we show that, conversely, when interacting with a droplet or a bubble, the deformation behaviour exhibits entirely different and unexpected behaviour. It arises from an interfacial dynamics which is specific to particles with free interfaces, namely thermal Marangoni forces. Our study employs a combination of experiments, theory, and numerical simulations to investigate the interaction between the droplet and the freezing front and unveils its surprising behaviour. In particular, we quantitatively understand the dependence of the front deformation $\Delta$ on the front propagation velocity, which, for large front velocities, can even revert from attraction ($\Delta < 0$) to repulsion ($\Delta > 0$).

The evaporation of water/1,2-hexanediol binary drops shows remarkable segregation dynamics, with hexanediol-rich spots forming at the rim, thus breaking axisymmetry. While the segregation of hexanediol near the rim can be attributed to the preferential evaporation of water, the symmetry-breaking and spot formation could not yet be successfully explained. With three-dimensional simulations and azimuthal stability analysis of a minimal model, we investigate the flow and composition in the drop. We show that a slightly non-monotonic surface tension causes the emergence of a counter-rotating Marangoni vortex in the hexanediol-rich rim region, which subsequently becomes azimuthally unstable and forms the observed spots. Accurate measurements reveal that the surface tension is indeed non-monotonic (~0.3 mN/m). This work provides valuable insight for applications like inkjet printing or spray cooling.

Water usually contains dissolved gases, and because freezing is a purifying process these gases must be expelled for ice to form. Bubbles appear at the freezing front and are trapped in ice, making pores. These pores come in a range of sizes from microns to millimeters and their shapes are peculiar; never spherical but elongated, and usually fore-aft asymmetric. We show that these remarkable shapes result of a delicate balance between freezing, capillarity, and mass diffusion. A non-linear ordinary differential equation suffices to describe the bubbles, with only two non-dimensional parameters representing the supersaturation and the freezing rate. Our experiments provide us with a large variety of pictures of bubble shapes. We show that all of them have their rounded tip well described by an asymptotic regime of the differential equation, and that most of them can have their full shape quantitatively matched by a full solution. This enables the measurement of the freezing conditions of ice samples, and the design of freeze-cast porous materials. Furthermore, the equation exhibits a bifurcation that explains why some bubbles grow indefinitely and make long cylindrical ``ice worms'', well known to glaciologists., Comment: 14 pages, 9 figures

The evolution of electrogenerated gas bubbles during water electrolysis can significantly hamper the overall process efficiency. Promoting the departure of electrochemically generated bubbles during (water) electrolysis is therefore beneficial. For a single bubble, a departure from the electrode surface occurs when buoyancy wins over the downward-acting forces (e.g. contact, Marangoni, and electric forces). In this work, the dynamics of a pair of H$_2$ bubbles produced during hydrogen evolution reaction in 0.5 M H$_2$SO$_4$ using dual platinum micro-electrode system is systematically studied by varying the electrode distance and the cathodic potential. By combining high-speed imaging and electrochemical analysis, we demonstrate the importance of bubble-bubble interactions for the departure process. We show that bubble coalescence may lead to substantially earlier bubble departure as compared to buoyancy effects alone, resulting in considerably higher reaction rates at constant potential. However, due to continued mass input and conservation of momentum repeated coalescence events with bubbles close to the electrode may drive departed bubbles back to the surface beyond a critical current, which increases with the electrode spacing. The latter leads to the resumption of bubble growth near the electrode surface, followed by buoyancy-driven departure. While less favourable at small electrode spacing, this configuration proves to be very beneficial at larger separations increasing the mean current up to 2.4 times compared to a single electrode under the conditions explored in this study., Comment: manuscript: 27 pages, 7 figures; SI: 13 pages, 9 figures

Frozen water might appear opaque since gas bubbles can get trapped in the ice during the freezing process. They nucleate and then grow near the advancing solidification front, due to the formation of a gas supersaturation region in its vicinity. A delicate interplay between the rate of mass transfer and the rate of freezing dictates the final shapes and sizes of the entrapped gas bubbles. In this work, we experimentally and numerically investigate the initial growth of such gas bubbles that nucleate and grow near the advancing ice front. We show that the initial growth of these bubbles is governed by diffusion and is enhanced due to a combination of the presence of the background gas concentration gradient and the motion of the approaching front. Additionally, we recast the problem into that of mass transfer to a moving spherical object in a homogeneous concentration field, finding good agreement between our experimental data and the existing scaling relations for that latter problem. Lastly, we address how fluid flow around the bubble might further affect this growth and qualitatively explore this through numerical simulations., Comment: 19 pages, 10 figures

SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons., Comment: Submitted 20th Apr. 2023, published 25th January 2024 (accepted version)

We report on the melting dynamics of ice suspended in fresh water and subject to natural convective flows. Using direct numerical simulations we investigate the melt rate of ellipsoidal objects for $2.32\times 10^4 \leq \text{Ra} \leq 7.61\times 10^8$, where \text{Ra} is the Rayleigh number defined with the temperature difference between the ice and the surrounding water. We reveal that the system exhibits non-monotonic behavior in three control parameters. As a function of the aspect ratio of the ellipsoidal, the melting time shows a distinct minimum that is different from a disk which has the minimum perimeter. Furthermore, also with \text{Ra} the system shows a non-monotonic trend, since for large \text{Ra} and large aspect ratio the flow separates, leading to distinctly different dynamics. Lastly, since the density of water is non-monotonic with temperature, the melt rate depends non-monotonically also on the ambient temperature, as for intermediate temperatures ($\unit{4}{\celsius}$--$\unit{7}{\celsius}$) the flow is (partially) reversed. In general, the shape which melts the slowest is quite distinct from that of a disk.

Solidification of complex liquids is pertinent to numerous natural and industrial processes. Here, we examine the freezing of a W/O/W double-emulsion, i.e., water-in-oil compound droplets dispersed in water. We show that the solidification of such hierarchical emulsions can trigger a topological transition; for example, in our case, we observe the transition from the stable W/O/W state to a (frozen) O/W single-emulsion configuration. Strikingly, this transition is characterised by sudden expulsion of the inner water drop from the encapsulating oil droplet. We propose that this topological transition is triggered by the freezing of the encapsulating oil droplet from the outside in, putting tension on the inner water drop thus, destabilizing the W/O/W configuration. Using high-speed imaging we characterize the destabilization process. Interestingly, we find that below a critical size of the inner drop, $R_{\mathrm{in,crit}} \approx 19 \, \mu\mathrm{m}$, the topological transition does not occur any more and the double-emulsion remains stable, in line with our interpretation., Comment: 10 pages, 7 figures

Particles at liquid interfaces have the tendency to cluster due to capillary forces competing to gravitational buoyancy (i.e., normal to the distorted free surface). This is known as the Cheerios effect. Here we experimentally and theoretically study the interaction between two submerged particles near an advancing water-ice interface during the freezing process. Particles that are thermally more conductive than water are observed to attract each other and form clusters once frozen. We call this feature the frozen Cheerios effect, where interactions are driven by alterations to the direction of the experienced repelling force (i.e., normal to the distorted isotherm). On the other hand, particles less conductive than water separate, highlighting the importance of thermal conduction during freezing. Based on existing models for single particle trapping in ice, we develop an understanding of multiple particle interaction. We find that the overall efficacy of the particle-particle interaction critically depends on the solidification front velocity. Our theory explains why the thermal conductivity mismatch between the particles and water dictates the attractive/repulsive nature of the particle-particle interaction., Comment: 15 pages, 10 figures

A liquid drop impacting a rigid substrate undergoes deformation and spreading due to normal reaction forces, which are counteracted by surface tension. On a non-wetting substrate, the drop subsequently retracts and takes off. Our recent work (Zhang et al., \textit{Phys. Rev. Lett.}, vol. 129, 2022, 104501) revealed two peaks in the temporal evolution of the normal force $F(t)$ -- one at impact and another at jump-off. The second peak coincides with a Worthington jet formation, which vanishes at high viscosities due to increased viscous dissipation affecting flow focusing. In this article, using experiments, direct numerical simulations, and scaling arguments, we characterize both the peak amplitude $F_1$ at impact and the one at take off ($F_2$) and elucidate their dependency on the control parameters: the Weber number $We$ (dimensionless impact kinetic energy) and the Ohnesorge number $Oh$ (dimensionless viscosity). The first peak amplitude $F_1$ and the time $t_1$ to reach it depend on inertial timescales for low viscosity liquids, remaining nearly constant for viscosities up to 100 times that of water. For high viscosity liquids, we balance the rate of change in kinetic energy with viscous dissipation to obtain new scaling laws: $F_1/F_\rho \sim \sqrt{Oh}$ and $t_1/\tau_\rho \sim 1/\sqrt{Oh}$, where $F_\rho$ and $\tau_\rho$ are the inertial force and time scales, respectively, which are consistent with our data. The time $t_2$ at which the amplitude $F_2$ appears is set by the inertio-capillary timescale $\tau_\gamma$, independent of both the viscosity and the impact velocity of the drop. However, these properties dictate the magnitude of this amplitude., Comment: accepted version, minor changes from the last version, Supplementary Videos: https://youtube.com/playlist?list=PLf5C5HCrvhLH8BwyeMSusdvyFmQYGP2qt&si=h9cQv-2luiqKC9Ar

Direct numerical simulations are utilised to investigate mass transfer processes at gas-evolving electrodes that experience successive formation and detachment of bubbles. The gas-liquid interface is modeled employing an Immersed Boundary Method. We simulate the growth phase of the bubbles followed by their departure from the electrode surface in order to study the mixing induced by these processes. We find that the growth of the bubbles switches from a diffusion-limited mode at low to moderate fractional bubble-coverages of the electrode to reaction-limited growth dynamics at high coverages. Furthermore, our results indicate that the net transport within the system is governed by the effective buoyancy driving induced by the rising bubbles and that mechanisms commonly subsumed under the term `microconvection' do not significantly affect the mass transport. Consequently, the resulting gas transport for different bubble sizes, current densities, and electrode coverages can be collapsed onto one single curve and only depends on an effective Grashof number. The same holds for the mixing of the electrolyte when additionally taking the effect of surface blockage by attached bubbles into account. For the gas transport to the bubble, we find that the relevant Sherwood numbers also collapse onto a single curve when accounting for the driving force of bubble growth, incorporated in an effective Jakob number. Finally, linking the hydrogen transfer rates at the electrode and the bubble interface, an approximate correlation for the gas-evolution efficiency has been established. Taken together, these findings enable us to deduce parametrizations for all response parameters of the systems.

Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a few gigaelectronvolt (GeV). Using the H.E.S.S. array of Cherenkov telescopes, we discovered a novel radiation component emerging beyond this generic GeV cutoff in the Vela pulsar's broadband spectrum. The extension of gamma-ray pulsation energies up to at least 20 teraelectronvolts (TeV) shows that Vela pulsar can accelerate particles to Lorentz factors higher than $4\times10^7$. This is an order of magnitude larger than in the case of the Crab pulsar, the only other pulsar detected in the TeV energy range. Our results challenge the state-of-the-art models for high-energy emission of pulsars while providing a new probe, i.e. the energetic multi-TeV component, for constraining the acceleration and emission processes in their extreme energy limit., Comment: 38 pages, 6 figures. This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, Nat Astron (2023), and is available online at https://doi.org/10.1038/s41550-023-02052-3

We consider the process of convective dissolution in homogeneous and isotropic porous media. The flow is unstable due to the presence of a solute that induces a density difference responsible for driving the flow. The mixing dynamics is thus driven by a Rayleigh-Taylor instability at the pore scale. We investigate the flow at the scale of the pores using Hele-Shaw type experiment with bead packs, 2D DNS and physical models. Experiments and simulations have been specifically designed to mimic the same flow conditions, namely matching porosities, high Schmidt numbers, and linear dependency of fluid density with solute concentration. In addition, the solid obstacles of the medium are impermeable to fluid and solute. We characterise the evolution of the flow via the mixing length, which quantifies the extension of the mixing region and grows linearly in time. The flow structure, analysed via the centre-line mean wavelength, is observed to grow in agreement with theoretical predictions. Finally, we analyse the dissolution dynamics of the system, quantified through the mean scalar dissipation, and three mixing regimes are observed: (i) The evolution is controlled by diffusion, which produces solute mixing across the initial horizontal interface; (ii) when the interfacial diffusive layer is sufficiently thick, it becomes unstable, forming finger-like structures and driving the system into a convection-dominated phase; (iii) when the fingers have grown sufficiently to touch the horizontal boundaries of the domain, the mixing reduces dramatically due to the absence of fresh unmixed fluid. With simple physical models, we explain the physics of the results obtained numerically and experimentally. The solute evolution presents a self-similar behaviour, and it is controlled by different length scales in each stage of the mixing process, namely length scale of diffusion, pore size, and domain height.

Iceberg melting is a critical factor for climate change, contributing to rising sea levels and climate change. However, the shape of an iceberg is an often neglected aspect of its melting process. Our study investigates the influence of different ice shapes and ambient flow velocities on melt rates by conducting direct numerical simulations. Our study focuses on the ellipsoidal shape, with the aspect ratio as the control parameter. It plays a crucial role in the melting process, resulting in significant variations in the melt rate between different shapes. Without flow, the optimal shape for a minimal melt rate is the disk (2D) or sphere (3D), due to the minimal surface area. However, as the ambient flow velocity increases, the optimal shape changes with the aspect ratio. We find that ice with an elliptical shape (when the long axis is aligned with the flow direction) can melt up to 10\% slower than a circular shape when exposed to flowing water. We provide a quantitative theoretical explanation for this optimal shape, based on the competition between surface area effects and convective heat transfer effects. Our findings provide insight into the interplay between phase transitions and ambient flows, contributing to our understanding of the iceberg melting process and highlighting the need to consider the aspect ratio effect in estimates of iceberg melt rates.

The presence of bubbles in a turbulent flow changes the flow drastically and enhances the mixing. Adding salt to the bubbly aqueous flow changes the bubble coalescence properties as compared to pure water. Here we provide direct experimental evidence that also the turbulent thermal energy spectra are strongly changed. Experiments were performed in the Twente Mass and Heat Transfer water tunnel,in which we can measure the thermal spectra in bubbly turbulence in salty water. We find that the mean bubble diameter decreases with increasing concentration of salt (NaCl), due to the inhibition of bubble coalescence. With increasing salinity, the transition frequency from the classical $-5/3$ scaling of the thermal energy spectrum to the bubble induced $-3$ scaling shifts to higher frequencies, thus enhancing the overall thermal energy. We relate this frequency shift to the smaller size of the bubbles for the salty bubbly flow. Finally we measure the heat transport in the bubbly flow, and show how it varies with changing void fraction and salinity: Increases in both result into increases in the number of extreme events., Comment: 18 pages, 10 figures, submitted to International Journal of Multiphase Flow

While artificial intelligence (AI) has made many successful applications in various domains, its adoption in healthcare lags a little bit behind other high-stakes settings. Several factors contribute to this slower uptake, including regulatory frameworks, patient privacy concerns, and data heterogeneity. However, one significant challenge that impedes the implementation of AI in healthcare, particularly in orthopedics, is the lack of explainability and interpretability around AI models. Addressing the challenge of explainable AI (XAI) in orthopedics requires developing AI models and algorithms that prioritize transparency and interpretability, allowing clinicians, surgeons, and patients to understand the contributing factors behind any AI-powered predictive or descriptive models. The current contribution outlines several key challenges and opportunities that manifest in XAI in orthopedic practice. This work emphasizes the need for interdisciplinary collaborations between AI practitioners, orthopedic specialists, and regulatory entities to establish standards and guidelines for the adoption of XAI in orthopedics., Comment: This paper was accepted at The 2023 World Congress in Computer Science, Computer Engineering, and Applied Computing (CSCE'23)

Two-dimensional horizontally periodic Rayleigh-B\'enard convection between stress-free boundaries displays two distinct types of states, depending on the initial conditions. Roll states are composed of pairs of counter-rotating convection rolls. Windy states are dominated by strong horizontal wind (also called zonal flow) that is vertically sheared, precludes convection rolls, and suppresses heat transport. Windy states occur only when the Rayleigh number $Ra$ is sufficiently above the onset of convection. At intermediate $Ra$ values, windy states can be induced by suitable initial conditions, but they undergo a transition to roll states after finite lifetimes. At larger $Ra$ values, where windy states have been observed for the full duration of simulations, it is unknown whether they represent chaotic attractors or only metastable states that would eventually undergo a transition to roll states. We study this question using direct numerical simulations of a fluid with a Prandtl number of 10 in a layer whose horizontal period is 8 times its height. At each of seven $Ra$ values between $9\times10^6$ and $2.25\times10^7$ we have carried out 200 or more simulations, all from initial conditions leading to windy convection with finite lifetimes. The lifetime statistics at each $Ra$ indicate a memoryless process with survival probability decreasing exponentially in time. The mean lifetimes grow with $Ra$ approximately as $Ra^4$. This analysis provides no $Ra$ value at which windy convection becomes stable; it might remain metastable at larger $Ra$ with extremely long lifetimes.

Automatic segmentation of knee bony anatomy is essential in orthopedics, and it has been around for several years in both pre-operative and post-operative settings. While deep learning algorithms have demonstrated exceptional performance in medical image analysis, the assessment of fairness and potential biases within these models remains limited. This study aims to revisit deep learning-powered knee-bony anatomy segmentation using plain radiographs to uncover visible gender and racial biases. The current contribution offers the potential to advance our understanding of biases, and it provides practical insights for researchers and practitioners in medical imaging. The proposed mitigation strategies mitigate gender and racial biases, ensuring fair and unbiased segmentation results. Furthermore, this work promotes equal access to accurate diagnoses and treatment outcomes for diverse patient populations, fostering equitable and inclusive healthcare provision., Comment: This paper has been accepted by IEEE BHI 2023

Rotational effects are commonly neglected when considering the dynamics of freely rising or settling isotropic particles. Here, we demonstrate that particle rotations play an important role for rising as well as for settling cylinders in situations when mass eccentricity, and thereby a new pendulum timescale, is introduced to the system. We employ two-dimensional simulations to study the motion of a single cylinder in a quiescent unbounded incompressible Newtonian fluid. This allows us to vary the Galileo number, density ratio, relative moment of inertia, and Centre-Of-Mass offset (COM) systematically and beyond what is feasible experimentally. For certain buoyant density ratios, the particle dynamics exhibit a resonance mode, during which the coupling via the Magnus lift force causes a positive feedback between translational and rotational motions. This mode results in vastly different trajectories with significantly larger rotational and translational amplitudes and an increase of the drag coefficient easily exceeding a factor two. We propose a simple model that captures how the occurrence of the COM offset induced resonance regime varies, depending on the other input parameters, specifically the density ratio, the Galileo number, and the relative moment of inertia. Remarkably, depending on the input parameters, resonance can be observed for centre-of-mass offsets as small as a few percent of the particle diameter, showing that the particle dynamics can be highly sensitive to this parameter., Comment: 37 pages, 17 figures