Search

The French Riviera is very often threatened by flash floods. These hydro-meteorological events, which are fast and violent, have catastrophic consequences on life and property. The development of forecasting tools may help to limit the impacts of these extreme events. Our purpose here is to demonstrate the possibility of using b-flood (a subset of the Basilisk library http://basilisk.fr/, last access: 8 November 2021), which is a 2D tool based on the shallow-water equations and adaptive mesh refinement. The code is first validated using analytical test cases describing different flow regimes. It is then applied to the Toce river valley physical model produced by ENEL-HYDRO in the framework of the CADAM project and on a flash-flood case over the urbanized Toce area produced during the IMPACT project. Finally, b-flood is applied to the flash flood of October 2015 in Cannes in south-eastern France, which demonstrates the feasibility of using software based on the shallow-water equations and mesh refinement for flash-flood simulation in small watersheds (less than 100 km2) and on a predictive computational timescale.

We study the coalescence of two drops of an ideal fluid driven by surface tension. The velocity of approach is taken to be zero and the dynamical effect of the outer fluid (usually air) is neglected. Our approximation is expected to be valid on scales larger than $\ell_{\nu} = \rho\nu^2/\sigma$, which is 10 nm for water. Using a high-precision boundary integral method, we show that the walls of the thin retracting sheet of air between the drops reconnect in finite time to form a toroidal enclosure. After the initial reconnection, retraction starts again, leading to a rapid sequence of enclosures. Averaging over the discrete events, we find the minimum radius of the liquid bridge connecting the two drops to scale like $r_b \propto t^{1/2}$. [ABSTRACT FROM AUTHOR]

Superhydrophobic surfaces generate very high contact angles as a result of their microstructure. The impact of a water drop on such a surface shows unusual features, such as total rebound at low impact speed. We report experimental and numerical investigations of the impact of approximately spherical water drops. The axisymmetric free surface problem, governed by the Navier–Stokes equations, is solved numerically with a front-tracking marker-chain method on a square grid. Experimental observations at moderate velocities and capillary wavelength much less than the initial drop radius show that the drop evolves to a staircase pyramid and eventually to a torus. Our numerical simulations reproduce this effect. The maximal radius obtained in numerical simulations precisely matches the experimental value. However, the large velocity limit has not been reached experimentally or numerically. We discuss several complications that arise at large velocity: swirling motions observed in the cross-section of the toroidal drop and the appearance of a thin film in the centre of the toroidal drop. The numerical results predict the dry-out of this film for sufficiently high Reynolds and Weber numbers. When the drop rebounds, it has a top-heavy shape. In this final stage, the kinetic energy is a small fraction of its initial value. [ABSTRACT FROM AUTHOR]

We propose the first continuous model with long range screening (shadowing) that describes columnar growth in one space dimension, as observed in plasma sputter deposition. It is based on a new continuous partial derivative equation with non-linear diffusion and where the shadowing effects apply on all the different processes. [ABSTRACT FROM AUTHOR]

The effect of freezing on contact line motion is a scientific challenge in the understanding of the solidification of capillary flows. In this Letter, we experimentally investigate the spreading and freezing of a water droplet on a cold substrate. We demonstrate that solidification stops the spreading because the ice crystals catch up with the advancing contact line. Indeed, we observe the formation and growth of ice crystals along the substrate during the drop spreading, and show that their velocity equals the contact line velocity when the drop stops. Modeling the growth of the crystals, we predict the shape of the crystal front and show that the substrate thermal properties play a major role on the frozen drop radius.

The COVID-19 pandemic has generated many concerns about cross-contamination risks, particularly in hospital settings and Intensive Care Units (ICU). Virus-laden aerosols produced by infected patients can propagate throughout ventilated rooms and put medical personnel entering them at risk. Experimental results found with a schlieren optical method have shown that the air flows generated by a cough and normal breathing were modified by the oxygenation technique used, especially when using High Flow Nasal Canulae, increasing the shedding of potentially infectious airborne particles. This study also uses a 3D Computational Fluid Dynamics model based on a Lattice Boltzmann Method to simulate the air flows as well as the movement of numerous airborne particles produced by a patient's cough within an ICU room under negative pressure. The effects of different mitigation scenarii on the amount of aerosols potentially containing SARS-CoV-2 that are extracted through the ventilation system are investigated. Numerical results indicate that adequate bed orientation and additional air treatment unit positioning can increase by 40% the number of particles extracted and decrease by 25% the amount of particles deposited on surfaces 45s after shedding. This approach could help lay the grounds for a more comprehensive way to tackle contamination risks in hospitals, as the model can be seen as a proof of concept and be adapted to any room configuration.

Thin liquid or gas films are everywhere in nature, from foams to submillimetric bubbles at a free surface, and their rupture leaves a collection of small drops and bubbles. However, the mechanisms at play responsible for the bursting of these films is still in debate. The present study thus aims at understanding the drainage dynamics of the thin air film squeezed by gravity between a millimetric droplet and a smooth solid or a liquid thin film. Solving coupled lubrication equations and analyzing the dominant terms in the solid- and liquid-film cases, we explain why the drainage is much faster in the liquid-film case, leading often to a shorter coalescence time, as observed in recent experiments., Competing Interests: The authors declare no competing interest.

Following the idea that dissipation in turbulence at high Reynolds number is dominated by singular events in space-time and described by solutions of the inviscid Euler equations, we draw the conclusion that in such flows, scaling laws should depend only on quantities appearing in the Euler equations. This excludes viscosity or a turbulent length as scaling parameters and constrains drastically possible analytical pictures of this limit. We focus on the drag law deduced by Newton for a projectile moving quickly in a fluid at rest. Inspired by this Newton's drag force law (proportional to the square of the speed of the moving object in the limit of large Reynolds numbers), which is well verified in experiments when the location of the detachment of the boundary layer is defined, we propose an explicit relationship between the Reynolds stress in the turbulent wake and quantities depending on the velocity field (averaged in time but depending on space). This model takes the form of an integrodifferential equation for the velocity which is eventually solved for a Poiseuille flow in a circular pipe.

In head-on collisions between two droplets, reflexive separation is frequently formed, showing tentative coalescence followed by disintegration into two primary drops. With higher impact inertia relative to surface tension, characterized by a Weber number (We), more satellite droplets are created between the primary drops. In the symmetric configuration, the existing phenomenological models indicate the absence of satellite droplets at the onset We when the coalesced drops start to break. Supported by experimental and simulation evidence, here we demonstrate the exclusive formation of at least one droplet after pinch of the thread connecting the colliding drops. In accordance with the universal features of a thinning liquid filament approaching singularity as predicted by scaling theories of pinch-off, the mechanism of satellite droplet formation in the symmetrical impact of droplets is clarified. Via slight breaking of the symmetry, no satellite droplet can be observed, thus providing a possible interpretation for the discrepancy in the literature and implications for controlling undesirable drop formation.

When an object is placed at the surface of a liquid, its weight deforms the interface. For two identical spherical objects, such a deformation creates an attractive force, leading to the aggregation of the two-body system. Here, we experimentally investigate the interaction between two granular rafts, formed by the aggregation of dense millimeter-sized beads placed at an oil-water interface. The interfacial deformation created by such a two-dimensional object exceeds by at least an order of magnitude the deformation of a single bead. This leads to unusually high capillary forces which strongly depend on the number of particles. Likewise, because the raft grows in size as more particles are added, the viscous drag experienced increases along with the capillary attraction, leading to a non-trivial dependence of the balance of forces on the number of beads. By studying the relative motion of two granular rafts in relation with the interfacial deformation they generate, we derive a model for the observed speed profiles. With this work, we generalize how the capillary interaction between two non-identical complex structures evolves with their respective geometry.

We investigate the impact of liquid drops on deep pools of aqueous glycerol solutions with variable pool viscosity and air pressure both experimentally and numerically. With this approach, we are able to address drop impacts on substrates that continuously transition from low-viscosity liquids to almost solids. We show that the generic corolla spreading out from the impact point consists of two distinct sheets, namely an ejecta sheet fed by the drop liquid and a second sheet fed by the substrate liquid, which evolve on separated timescales. These two sheets contribute to a varying extent to the corolla overall dynamics and splashing, depending on the viscosity ratio between the two liquids.

Many differential equations involved in natural sciences show singular behaviors; i.e., quantities in the model diverge as the solution goes to zero. Nonetheless, the evolution of the singularity can be captured with self-similar solutions, several of which may exist for a given system. How to characterize the transition from one self-similar regime to another remains an open question. By studying the classic example of the pinch-off of a viscous liquid thread, we show experimentally that the geometry of the system and external perturbations play an essential role in the transition from a symmetric to an asymmetric solution. Moreover, this transient regime undergoes unexpected log-scale oscillations that delay dramatically the onset of the final self-similar solution. This result sheds light on the strong impact external constraints can have on predictions established to explain the formation of satellite droplets or on the rheological tests applied on a fluid, for example., Competing Interests: The authors declare no conflict of interest.

Solid particles can modify the properties of liquid interfaces and are therefore widely used to coat drops, bubbles, and stabilize emulsions and foams. Here, we propose a new, easy, and affordable method to produce millimetric to centimetric water-in-water capsules using solid particles. We prevent the coalescence of a water drop at an oil-water interface using a monolayer of large, dense, and hydrophobic particles: a "granular raft". The capsule is then formed by a mechanical instability occurring when the interface collapses under the combined load of the floating drop and particle weight. During the destabilization, the water drop sinks into the water subphase through an oil-particle film which covers it to produce the armored capsule. By modeling the raft as a heavy membrane, we predict the floating drop shape, the raft deformation, its destabilization and highlight the complex dual nature (solid- and liquid-like) of the capsule shell. Because armored capsules' content is isolated, transportable, and easily releasable, they are great candidates for applications requiring transport of water-soluble compounds in aqueous systems such as green chemistry or cell biology.

The Sagdeev-Zaslavski (SZ) equation for wave turbulence is analytically derived, both in terms of a generating function and of a multipoint probability density function (PDF), for weakly interacting waves with initial random phases. When the initial amplitudes are also random, a one-point PDF equation is derived. Such analytical calculations remarkably agree with results obtained in totally different fashions. Numerical investigations of the two-dimensional nonlinear Schrödinger equation (NLSE) and of a vibrating plate prove the following: (i) Generic Hamiltonian four-wave systems rapidly attain a random distribution of phases independently of the slower dynamics of the amplitudes, vindicating the hypothesis of initially random phases. (ii) Relaxation of the Fourier amplitudes to the predicted stationary distribution (exponential) happens on a faster time scale than relaxation of the spectrum (Rayleigh-Jeans distribution). (iii) The PDF equation correctly describes dynamics under different forcings: The NLSE has an exponential PDF corresponding to a quasi-Gaussian solution, as the vibrating plates, that also shows some intermittency at very strong forcings.

Instruments that belong to the gong family exhibit nonlinear dynamics at large amplitudes of vibration. In the specific case of the xiaoluo gong, this nonlinear behavior results in a pitch glide of several modes of the instrument in addition to harmonic distortion and internal resonances. This study applies a linear modal active control to a xiaoluo gong in an attempt to change its sound properties. First, a modal damping control of the fundamental mode based on a linear identification and a state space controller is applied in the small amplitude regime (no pitch glide). Results indicate that modal control influences not only the controlled mode but also the frequency components involved in distortion or internal resonance phenomena. Second, a modal damping control is performed in the large amplitude regime (in the presence of pitch glide). Results show that modal control does not affect the pitch glide. However, the controller becomes effective at a time trigger which is related to the instantaneous frequency.

We report experiments that yield new insights on the behavior of granular rafts at an oil-water interface. We show that these particle aggregates can float or sink depending on dimensionless parameters taking into account the particle densities and size and the densities of the two fluids. We characterize the raft shape and stability and propose a model to predict its shape and maximum length to remain afloat. Finally we find that wrinkles and folds appear along the raft due to compression by its own weight, which can trigger destabilization. These features are characteristics of an elastic instability, which we discuss, including the limitations of our model.

We investigate experimentally the quenching of a liquid pancake, obtained through the impact of a water drop on a cold solid substrate (0 °C to -60 °C). We show that, below a certain substrate temperature, fractures appear on the frozen pancake and the crack patterns change from a 2D fragmentation regime to a hierarchical fracture regime as the thermal shock increases. The different regimes are discussed and the transition temperatures are estimated through classical fracture scaling arguments. Finally, a phase diagram presents how these regimes can be controlled by the drop impact parameters.

We investigate the onset of intermittency for vibrating elastic plate turbulence in the framework of the weak wave turbulence theory using a numerical approach. The spectrum of the displacement field and the structure functions of the fluctuations are computed for different forcing amplitudes. At low forcing, the spectrum predicted by the theory is observed, while the fluctuations are consistent with Gaussian statistics. When the forcing is increased, the spectrum varies at large scales, corresponding to the oscillations of nonlinear structures made of ridges delimited by d cones. In this regime, the fluctuations exhibit small-scale intermittency that can be fitted via a multifractal model. The analysis of the nonlinear frequency shows that the intermittency is linked to the breakdown of the weak turbulence at large scales only.

The dynamics of random weakly nonlinear waves is studied in the framework of vibrating thin elastic plates. Although it has been previously predicted that no stationary inverse cascade of constant wave action flux could exist in the framework of wave turbulence for elastic plates, we present substantial evidence of the existence of a time-dependent inverse cascade, opening up the possibility of self-organization for a larger class of systems. This inverse cascade transports the spectral density of the amplitude of the waves from short up to large scales, increasing the distribution of long waves despite the short-wave fluctuations. This dynamics appears to be self-similar and possesses a power-law behavior in the short-wavelength limit which significantly differs from the exponent obtained via a Kolmogorov dimensional analysis argument. Finally, we show explicitly a tendency to build a long-wave coherent structure in finite time.

We study the dynamical regime of wave turbulence of a vibrated thin elastic plate based on experimental and numerical observations. We focus our study on the strongly nonlinear regime described in a previous Letter by Yokoyama and Takaoka. At small forcing, a weakly nonlinear regime is compatible with the weak turbulence theory when the dissipation is localized at high wave number. When the forcing intensity is increased, a strongly nonlinear regime emerges: singular structures dominate the dynamics at large scales whereas at small scales the weak turbulence is still present. A turbulence of singular structures with folds and D cones develops that alters significantly the energy spectra and causes the emergence of intermittency.

A study of the complex drying dynamics of polymeric mixtures with optical microscopy and gravimetric measurement is presented. Droplet formation is observed, followed by a collapse that leads to the residual craters in the dried film. The process is followed in situ under well-defined temperature and hygrometric conditions to determine the origin and nature of these droplets and craters. The drying process is usually completed within 1 h. The observations are explained using a simple diffusion model based on experimental results collected from mass and optical measurements as well as Raman confocal microspectrometry. Although the specific polymeric mixtures used here are of interest to the cosmetic industry, the general conclusions reached can apply to other polymeric aqueous solutions with applications to commercial and artistic painting.

We investigate theoretically and numerically a model of a supersolid in a dipole-blockaded Bose-Einstein condensate. The dependence of the superfluid fraction with an imposed thermal bath and a uniform boost velocity on the condensate is considered. Specifically, we observe a critical velocity for the nucleation of vortices in our system that is strongly linked to a steplike decrease in the superfluid fraction. We are able to use a scaling argument based on the energy required to activate a vortex, relating the critical temperature to the critical velocity, and find that this relationship is in good agreement with the numerical simulations carried out on the nonlocal Gross-Pitaevskii equation.

The splashing of a drop impacting onto a liquid pool produces a range of different sized microdroplets. At high impact velocities, the most significant source of these droplets is a thin liquid jet emerging at the start of the impact from the neck that connects the drop to the pool. We use ultrahigh-speed video imaging in combination with high-resolution numerical simulations to show how this ejecta gives way to irregular splashing. At higher Reynolds numbers, its base becomes unstable, shedding vortex rings into the liquid from the free surface in an axisymmetric von Kármán vortex street, thus breaking the ejecta sheet as it forms.

A drop impacting a target cutout in a thin polymer film is wrapped by the film in a dynamic sequence involving both capillary forces and inertia. Different 3D structures can be produced from a given target by slightly varying the impact parameters. A simplified model for a nonlinear dynamic Elastica coupled with a drop successfully explains this shape selection and yields detailed quantitative agreement with experiments. This first venture into the largely unexplored dynamics of elastocapillary assemblies opens up the perspective of mass production of 3D packages with individual shape selection.

We describe experimentally the formation of a pattern for drop impacts on thin liquid films for a large range of impact parameters. Using the shallow-water approximation, we are able to explain the main mechanisms leading to these patterns: it consists in the linear instability of the self-similar axisymmetric radial solution of the equations. Agreement between the experiments and the theory is remarkably good, leading, in particular, to the prediction that the most unstable fold number scales like (We/h∞)2/7.

We simulate the impact of a viscous liquid drop onto a smooth dry solid surface. As in experiments, when ambient air effects are negligible, impact flattens the falling drop without producing a splash. The no-slip boundary condition at the wall produces a boundary layer inside the liquid. Later, the flattening surface of the drop traces out the boundary layer. As a result, the eventual shape of the drop is a "pancake" of uniform thickness except at the rim, where surface tension effects are significant. The thickness of the pancake is simply the height where the drop surface first collides with the boundary layer.

Constant flux atom deposition into a porous medium is shown to generate a dense overlayer and a diffusion profile. Scaling analysis shows that the overlayer acts as a dynamic control for atomic diffusion in the porous substrate. This is modeled by generalizing the porous diffusion equation with a time-dependent diffusion coefficient equivalent to a nonlinear rescaling of time.

A two-dimensional model for the erosion generated by viscous free-surface flows, based on the shallow-water equations and the lubrication approximation, is presented. It has a family of self-similar solutions for straight erodible channels, with an aspect ratio that increases in time. It is also shown, through a simplified stability analysis, that a laminar river can generate various bar instabilities very similar to those observed in natural rivers. This theoretical similarity reflects the meandering and braiding tendencies of laminar rivers indicated by F. Métivier and P. Meunier [J. Hydrol. 27, 22 (2003)]. Finally, we propose a simple scenario for the transition between patterns observed in experimental erodible channels.

The mechanical behavior of a supersolid is studied in the framework of a fully explicit model derived from the Gross-Pitaevskii equation without assuming any defect or vacancy. A set of coupled nonlinear partial differential equations plus boundary conditions is derived. The conditions of mechanical equilibrium are studied under external constraints such as steady rotation and external stress. Our model explains the experimentally observed paradoxical behavior: a nonclassical rotational inertia fraction in the limit of small rotation speed but a solidlike elastic response to small stress or an external force field.

We study the long-time evolution of waves of a thin elastic plate in the limit of small deformation so that modes of oscillations interact weakly. According to the theory of weak turbulence (successfully applied in the past to plasma, optics, and hydrodynamic waves), this nonlinear wave system evolves at long times with a slow transfer of energy from one mode to another. We derive a kinetic equation for the spectral transfer in terms of the second order moment. We show that such a theory describes the approach to an equilibrium wave spectrum and represents also an energy cascade, often called the Kolmogorov-Zakharov spectrum. We perform numerical simulations that confirm this scenario.

We show that when water droplets gently impact on a hydrophobic surface, the droplet shoots out a violent jet, the velocity of which can be up to 40 times the drop impact speed. As a function of the impact velocity, two different hydrodynamic singularities are found that correspond to the collapse of the air cavity formed by the deformation of the drop at impact. It is the collapse that subsequently leads to the jet formation. We show that the divergence of the jet velocity can be understood using simple scaling arguments. In addition, we find that very large air bubbles can remain trapped in the drops. The surprising occurrence of the bubbles for low-speed impact is connected with the nature of the singularities, and can have important consequences for drop deposition, e.g., in ink-jet printing.

We study the formation of a large-scale coherent structure (a condensate) in classical wave equations by considering the defocusing nonlinear Schrödinger equation as a representative model. We formulate a thermodynamic description of the classical condensation process by using a wave turbulence theory with ultraviolet cutoff. In three dimensions the equilibrium state undergoes a phase transition for sufficiently low energy density, while no transition occurs in two dimensions, in complete analogy with standard Bose-Einstein condensation in quantum systems. On the basis of a modified wave turbulence theory, we show that the nonlinear interaction makes the transition to condensation subcritical. The theory is in quantitative agreement with the numerical integration of the nonlinear Schrödinger equation.

We investigate long time numerical simulations of the inviscid Rayleigh-Taylor instability at Atwood number one using a boundary integral method. We are able to attain the asymptotic behavior for the spikes predicted by Clavin and Williams for which we give a simplified demonstration. In particular, we observe that the spike's curvature evolves as t(3), while the overshoot in acceleration shows good agreement with the suggested 1/t(5) law. Moreover, we obtain consistent results for the prefactor coefficients of the asymptotic laws. Eventually we exhibit the self-similar behavior of the interface profile near the spike.

We show how giant vortices can be stabilized for strong external potentials in Bose-Einstein condensates. We illustrate the formation of these vortices thanks to the Ginzburg-Landau dissipative dynamics for two typical potentials in two spatial dimensions. The giant vortex stability is studied for the particular case of a rotating cylindrical hard wall. Due to axial symmetry the minimization of the perturbed energy is simplified into a one dimensional relaxation dynamics. Solving this 1D minimization problem, we observe that giant vortices are either never stable, or only stable in a finite frequency range. Finally we obtain the marginal curve for the minimum frequency needed to observe a giant vortex., (Copyright 2004 American Institute of Physics)

We propose a simple continuum model to interpret the shearing motion of dense, dry and cohesion-less granular media. Compressibility, dilatancy and Coulomb-like friction are the three basic ingredients. The granular stress is split into a rate-dependent part representing the rebound-less impacts between grains and a rate-independent part associated with long-lived contacts. Because we consider stationary flows only, the grain compaction and the grain velocity are the two main variables. The predicted velocity and compaction profiles are in apparent qualitative agreement with most of the experimental or numerical results concerning free-surface shear flows as well as confined shear flows.

We study the growth of a periodic pattern in one dimension for a model of spinodal decomposition, the Cahn-Hilliard equation. We particularly focus on the intermediate region, where the nonlinearity cannot be neglected anymore, and before the coalescence dominates. The dynamics is captured through the standard technique of a solubility condition performed over a particular family of quasistatic solutions. The main result is that the dynamics along this particular class of solutions can be expressed in terms of a simple ordinary differential equation. The density profile of the stationary regime found at the end of the nonlinear growth is also well characterized. Numerical simulations correspond satisfactorily to the analytical results through three different methods and asymptotic dynamics are well recovered, even far from the region where the approximations hold.

We present a combined experimental and theoretical study of memory effects in vibration-induced compaction of granular materials. In particular, the response of the system to an abrupt change in shaking intensity is measured. At short times after the perturbation a granular analog of aging in glasses is observed. Using a simple two-state model, we are able to explain this short-time response. We also discuss the possibility for the system to obey an approximate pseudo-fluctuation-dissipation theorem relationship and relate our work to earlier experimental and theoretical studies of the problem.

A theoretical and experimental study of patterns in vibrated granular layers is presented. An order parameter model based on the parametric Ginzburg-Landau equation is used to describe strongly nonlinear excitations including hexagons, interfaces between flat antiphase domains, and new localized objects, superoscillons. The experiments confirm the existence of superoscillons and bound states of superoscillons and interfaces. On the basis of the order parameter model we predict analytically and confirm experimentally that additional subharmonic driving results in the controlled motion of interfaces.

We present a statistical equilibrium model of self-organization in a class of focusing, nonintegrable nonlinear Schrodinger (NLS) equations. The theory predicts that the asymptotic-time behavior of the NLS system is characterized by the formation and persistence of a large-scale coherent solitary wave, which minimizes the Hamiltonian given the conserved particle number (L2-norm squared), coupled with small-scale random fluctuations, or radiation. The fluctuations account for the difference between the conserved value of the Hamiltonian and the Hamiltonian of the coherent state. The predictions of the statistical theory are tested against the results of direct numerical simulations of NLS, and excellent qualitative and quantitative agreement is demonstrated. In addition, a careful inspection of the numerical simulations reveals interesting features of the transitory dynamics leading up to the long-time statistical equilibrium state starting from a given initial condition. As time increases, the system investigates smaller and smaller scales, and it appears that at a given intermediate time after the coalescense of the soliton structures has ended, the system is nearly in statistical equilibrium over the modes that it has investigated up to that time.

We discuss the problem of an explosion in the cubic-quintic superfluid model, in relation to some experimental observations. We show numerically that an explosion in such a model might induce a cavitation bubble for large enough energy. This gives a consistent view for rebound bubbles in superfluid and we identify the loss of energy between the successive rebounds as radiated waves. We compute self-similar solution of the explosion for the early stage, when no bubbles have been nucleated. The solution also gives the wave number of the excitations emitted through the shock wave.