Jérémie Bec
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Jérémie Bec has a total of 69 co-authors. Here are the top co-authors by number of shared publications:
- Giorgio Krstulovic
Université Côte d'Azur
10 shared publications Holger Homann
Laboratoire J.-L. Lagrange, Université Côte d'Azur, OCA, CNRS, Bd. de l'Observatoire, 06300 Nice, France
9 shared publications- Massimo Cencini
Istituto dei Sistemi Complessi
9 shared publications - Federico Toschi
Department of Applied Physics and Science Education, Eindhoven University of Technology, Eindhoven, The Netherlands
6 shared publications - Luca Biferale
Department of Physics and INFN, University of Rome ‘Tor Vergata’, Via della Ricerca Scientifica 1, 00133 Rome, Italy
6 shared publications - Simon Thalabard
Institut de physique de Nice
6 shared publications
Latest Publications
- 1 March 2026
- Published by IOP Publishing in Europhysics Letters
- Vol. 153 (6) , 67001
Global manoeuvring of confined active particle swarms via emergent threshold-based control (a)
This study investigates the use of global control strategies to enhance the directed migration of swarms of interacting self-propelled particles confined in a channel. Uncontrolled dynamics naturally leads to wall accumulation, clogging, and band formation due to the interplay between self-organization, volume exclusion and confinement. This work explores whether a uniform global control, such as a magnetic field acting identically on all particles, can robustly enhance collective transport. Using a discrete Vicsek-like model, it is found that simple global alignment controls emerging from reinforcement learning, efficiently suppress unfavorable configurations and significantly increase the net particle flux along a prescribed channel direction. These results highlight that coarse, system-level observations are sufficient to achieve global manoeuvring by indirectly reshaping self-organized microscopic dynamics, even in regimes with strong fluctuations or partial ordering.
- 1 March 2026
- Published by IOP Publishing in Europhysics Letters
- Vol. 153 (6) , 67001
Global manoeuvring of confined active particle swarms via emergent threshold-based control (a)
This study investigates the use of global control strategies to enhance the directed migration of swarms of interacting self-propelled particles confined in a channel. Uncontrolled dynamics naturally leads to wall accumulation, clogging, and band formation due to the interplay between self-organization, volume exclusion and confinement. This work explores whether a uniform global control, such as a magnetic field acting identically on all particles, can robustly enhance collective transport. Using a discrete Vicsek-like model, it is found that simple global alignment controls emerging from reinforcement learning, efficiently suppress unfavorable configurations and significantly increase the net particle flux along a prescribed channel direction. These results highlight that coarse, system-level observations are sufficient to achieve global manoeuvring by indirectly reshaping self-organized microscopic dynamics, even in regimes with strong fluctuations or partial ordering.
- Open Access
- 15 October 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1021, A17
Surface quasi-geostrophic turbulence: the refined study of an active scalar
Surface quasi-geostrophic (SQG) theory describes the two-dimensional active transport of a scalar field, such as temperature, which – when properly rescaled – shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with fully developed three-dimensional turbulence. In particular, the Kraichnan – Leith – Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling for both scalar and velocity fields, and the presence of intermittency through multifractal scaling was pointed out by Sukhatme & Pierrehumbert (2002 Chaos12, 439–450), in unforced settings. In this work, we refine the discussion of these statistical analogies, using numerical simulations with up to $16\,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, coupling velocity increments with scalar differences, develop well-defined scaling ranges, highlighting the role of anomalous fluxes of all the scalar moments. However, the clean multiscaling properties of SQG transport are blurred when considering velocity and scalar fields separately. In particular, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires us to resort to generalised notions of scale invariance, such as extended self-similarity and a specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: they evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions that have been repeatedly observed in SQG numerics for the past 20 years.- 15 October 2025
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 1021, A17
Open AccessSurface quasi-geostrophic turbulence: the refined study of an active scalar
Surface quasi-geostrophic (SQG) theory describes the two-dimensional active transport of a scalar field, such as temperature, which – when properly rescaled – shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with fully developed three-dimensional turbulence. In particular, the Kraichnan – Leith – Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling for both scalar and velocity fields, and the presence of intermittency through multifractal scaling was pointed out by Sukhatme & Pierrehumbert (2002 Chaos12, 439–450), in unforced settings. In this work, we refine the discussion of these statistical analogies, using numerical simulations with up to $16\,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, coupling velocity increments with scalar differences, develop well-defined scaling ranges, highlighting the role of anomalous fluxes of all the scalar moments. However, the clean multiscaling properties of SQG transport are blurred when considering velocity and scalar fields separately. In particular, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires us to resort to generalised notions of scale invariance, such as extended self-similarity and a specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: they evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions that have been repeatedly observed in SQG numerics for the past 20 years. - Open Access
Surface quasigeostrophic turbulence: The refined study of an active scalar
SQG describes the 2D active transport of a scalar field, such as temperature, which -- when properly rescaled -- shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with 3D turbulence. In particular, the Kraichnan-Leith-Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling $\propto (\varepsilon \ell)^{1/3}$ for both scalar and velocity fields, and the presence of intermittency was pointed out by Sukhatme & Pierrehumbert ($Chaos$ $\boldsymbol{12}$, 439, 2002) in unforced settings. In this work, we refine these analogies using simulations up to $16,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, linking velocity increments with powers of scalar differences, exhibit clear scaling, revealing the role of anomalous fluxes of all the scalar moments. However, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires to resort to generalised notions of scale invariance, such as extended self-similarity and specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: They evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions, that have been repeatedly observed in SQG numerics for the past 20 years.Open AccessSurface quasigeostrophic turbulence: The refined study of an active scalar
SQG describes the 2D active transport of a scalar field, such as temperature, which -- when properly rescaled -- shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with 3D turbulence. In particular, the Kraichnan-Leith-Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling $\propto (\varepsilon \ell)^{1/3}$ for both scalar and velocity fields, and the presence of intermittency was pointed out by Sukhatme & Pierrehumbert ($Chaos$ $\boldsymbol{12}$, 439, 2002) in unforced settings. In this work, we refine these analogies using simulations up to $16,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, linking velocity increments with powers of scalar differences, exhibit clear scaling, revealing the role of anomalous fluxes of all the scalar moments. However, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires to resort to generalised notions of scale invariance, such as extended self-similarity and specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: They evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions, that have been repeatedly observed in SQG numerics for the past 20 years. - 2 December 2024
- Published by American Physical Society (APS) in Physical Review E
- Vol. 110 (6) , 064601
Effects of collective patterns, confinement, and fluid flow on active particle transport
The self-organization of active particles on a two-dimensional single-occupancy lattice is investigated, with an emphasis on the effects of boundary confinement and the influence of an external mean fluid flow. The study examines collective behaviors, particularly the transition from a disordered phase to the formation of orientationally ordered patterns, and their impact on particle transport and flux. In the absence of fluid flow, confinement causes particles to accumulate near the walls, leading to clogs or obstructions that hinder movement, or to the formation of bands aligned with the channel. Although these bands limit the particles' ability to freely self-propel, they still result in a net flux along the channel. The introduction of an external Poiseuille fluid flow induces vorticity, shifts the phase transition to higher alignment sensitivities, and promotes particle clustering at the channel center, significantly enhancing overall flux.
- 2 December 2024
- Published by American Physical Society (APS) in Physical Review E
- Vol. 110 (6) , 064601
Effects of collective patterns, confinement, and fluid flow on active particle transport
The self-organization of active particles on a two-dimensional single-occupancy lattice is investigated, with an emphasis on the effects of boundary confinement and the influence of an external mean fluid flow. The study examines collective behaviors, particularly the transition from a disordered phase to the formation of orientationally ordered patterns, and their impact on particle transport and flux. In the absence of fluid flow, confinement causes particles to accumulate near the walls, leading to clogs or obstructions that hinder movement, or to the formation of bands aligned with the channel. Although these bands limit the particles' ability to freely self-propel, they still result in a net flux along the channel. The introduction of an external Poiseuille fluid flow induces vorticity, shifts the phase transition to higher alignment sensitivities, and promotes particle clustering at the channel center, significantly enhancing overall flux.
- Open Access
- 25 November 2024
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 999, A83
Homogeneous turbophoresis of heavy inertial particles in turbulent flow
Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as ‘turbophoresis’, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space- and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number ${\textit {Pe}}_\ell$ , which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $\ell$ . Multifractal statistics of energy dissipation indicate that $ {\textit {Pe}}_\ell \sim \ell ^\delta /\tau _{p}$ with $\delta \approx 0.84$ . Numerical simulations support this behaviour and emphasise the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $\ell$ . This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes have a non-trivial distribution with a power-law tail $p(\mathcal {V}) \propto \mathcal {V}^{-\alpha }$ , with an exponent $\alpha$ that tends to 2 as ${\textit {Pe}}_\ell \to 0$ .- 25 November 2024
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 999, A83
Open AccessHomogeneous turbophoresis of heavy inertial particles in turbulent flow
Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as ‘turbophoresis’, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space- and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number ${\textit {Pe}}_\ell$ , which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $\ell$ . Multifractal statistics of energy dissipation indicate that $ {\textit {Pe}}_\ell \sim \ell ^\delta /\tau _{p}$ with $\delta \approx 0.84$ . Numerical simulations support this behaviour and emphasise the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $\ell$ . This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes have a non-trivial distribution with a power-law tail $p(\mathcal {V}) \propto \mathcal {V}^{-\alpha }$ , with an exponent $\alpha$ that tends to 2 as ${\textit {Pe}}_\ell \to 0$ . - 11 June 2024
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 9 (6) , L062501
Enhanced transport of flexible fibers by pole vaulting in turbulent wall-bounded flow
Long, flexible fibers transported by a turbulent channel flow sample nonlinear variations of the fluid velocity along their length. As the fibers tumble and collide with the boundaries, they bounce off with an impulse that propels them toward the center of the flow, similar to pole vaulting. As a result, the fibers migrate away from the walls, leading to depleted regions near the boundaries and more concentrated regions in the bulk. These higher concentrations in the center of the channel result in a greater net flux of fibers than what was initially imposed by the fluid. This effect becomes more pronounced as fiber length increases, especially when it approaches the channel height.
- 11 June 2024
- Published by American Physical Society (APS) in Physical Review Fluids
- Vol. 9 (6) , L062501
Enhanced transport of flexible fibers by pole vaulting in turbulent wall-bounded flow
Long, flexible fibers transported by a turbulent channel flow sample nonlinear variations of the fluid velocity along their length. As the fibers tumble and collide with the boundaries, they bounce off with an impulse that propels them toward the center of the flow, similar to pole vaulting. As a result, the fibers migrate away from the walls, leading to depleted regions near the boundaries and more concentrated regions in the bulk. These higher concentrations in the center of the channel result in a greater net flux of fibers than what was initially imposed by the fluid. This effect becomes more pronounced as fiber length increases, especially when it approaches the channel height.
- Open Access
- 6 May 2024
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 685, L4
Reduction of dust radial drift by turbulence in protoplanetary disks
Context. Dust particles in protoplanetary disks rotate at velocities exceeding those of the surrounding gas due to a lack of pressure support. Consequently, they experience a headwind from the gas that drives them toward the central star. Radial drift occurs on timescales much shorter than those inferred from disk observations or those required for dust to aggregate and form planets. Additionally, turbulence is often assumed to amplify the radial drift of dust in planet-forming disks when modeled through an effective viscous transport. However, the local interactions between turbulent eddies and particles are known to be significantly more intricate than in a viscous fluid. Aims. Our objective is to elucidate and characterize the dynamic effects of Keplerian turbulence on the mean radial and azimuthal velocities of dust particles. Methods. We employed 2D shearing-box incompressible simulations of the gas, which is maintained in a developed turbulent state while rotating at a sub-Keplerian speed. Dust is modeled as Lagrangian particles set at a Keplerian velocity, therefore experiencing a radial force toward the star through drag. Results. Turbulent eddies are found to reduce the radial drift, while simultaneously enhancing the azimuthal velocities of small particles. This dynamic behavior arises from the modification of dust trajectories due to turbulent eddies.
- 6 May 2024
- Published by EDP Sciences in Astronomy & Astrophysics
- Vol. 685, L4
Open AccessReduction of dust radial drift by turbulence in protoplanetary disks
Context. Dust particles in protoplanetary disks rotate at velocities exceeding those of the surrounding gas due to a lack of pressure support. Consequently, they experience a headwind from the gas that drives them toward the central star. Radial drift occurs on timescales much shorter than those inferred from disk observations or those required for dust to aggregate and form planets. Additionally, turbulence is often assumed to amplify the radial drift of dust in planet-forming disks when modeled through an effective viscous transport. However, the local interactions between turbulent eddies and particles are known to be significantly more intricate than in a viscous fluid. Aims. Our objective is to elucidate and characterize the dynamic effects of Keplerian turbulence on the mean radial and azimuthal velocities of dust particles. Methods. We employed 2D shearing-box incompressible simulations of the gas, which is maintained in a developed turbulent state while rotating at a sub-Keplerian speed. Dust is modeled as Lagrangian particles set at a Keplerian velocity, therefore experiencing a radial force toward the star through drag. Results. Turbulent eddies are found to reduce the radial drift, while simultaneously enhancing the azimuthal velocities of small particles. This dynamic behavior arises from the modification of dust trajectories due to turbulent eddies.
- Open Access
Homogeneous turbophoresis of heavy inertial particles in turbulent flow
Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as turbophoresis, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number $Pe_l$, which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $l$. Multifractal statistics of energy dissipation indicate that $Pe_l \sim l^δ/τ_{\rm p}$ with $δ\approx 0.84$. Numerical simulations support this behaviour and emphasises the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $l$. This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes has a non-trivial distribution with a power-law tail $p(V) \propto V^{-α}$, with an exponent $α$ that tends to 2 as $Pe_l\to0$.Open AccessHomogeneous turbophoresis of heavy inertial particles in turbulent flow
Heavy particles suspended in turbulent flow possess inertia and are ejected from violent vortical structures by centrifugal forces. Once piled up along particle paths, this small-scale mechanism leads to an effective large-scale drift. This phenomenon, known as turbophoresis, causes particles to leave highly turbulent regions and migrate towards calmer regions, explaining why particles transported by non-homogeneous flows tend to concentrate near the minima of turbulent kinetic energy. It is demonstrated here that turbophoretic effects are just as crucial in statistically homogeneous flows. Although the average turbulent activity is uniform, instantaneous spatial fluctuations are responsible for inertial-range inhomogeneities in the particle distribution. Direct numerical simulations are used to probe particle accelerations, specifically how they correlate to local turbulent activity, yielding an effective coarse-grained dynamics that accounts for particle detachment from the fluid and ejection from excited regions through a space and time-dependent non-Fickian diffusion. This leads to cast fluctuations in particle distributions in terms of a scale-dependent Péclet number $Pe_l$, which measures the importance of turbulent advection compared with inertial turbophoresis at a given scale $l$. Multifractal statistics of energy dissipation indicate that $Pe_l \sim l^δ/τ_{\rm p}$ with $δ\approx 0.84$. Numerical simulations support this behaviour and emphasises the relevance of the turbophoretic Péclet number in characterising how particle distributions, including their radial distribution function, depends on $l$. This approach also explains the presence of voids with inertial-range sizes, and the fact that their volumes has a non-trivial distribution with a power-law tail $p(V) \propto V^{-α}$, with an exponent $α$ that tends to 2 as $Pe_l\to0$. - 18 July 2023
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 967, R4
Velocity and acceleration statistics of heavy spheroidal particles in turbulence
Non-spherical particles transported by turbulent flow have a rich dynamics that combines their translational and rotational motions. Here, the focus is on small, heavy, inertial particles with a spheroidal shape fully prescribed by their aspect ratio. Such particles undergo an anisotropic, orientation-dependent viscous drag with the carrier fluid flow whose associated torque is given by the Jeffery equations. Direct numerical simulations of homogeneous, isotropic turbulence are performed to study systematically how the translational motion of such spheroidal particles depends on their shape and size. It is found that the Lagrangian statistics of both velocity and acceleration can be described in terms of an effective Stokes number obtained as an isotropic average over angles of the particle's orientation. Corrections to the translational motion of particles due to their non-sphericity and rotation can hence be recast as an effective radius obtained from such a mean.
- 18 July 2023
- Published by Cambridge University Press (CUP) in Journal of Fluid Mechanics
- Vol. 967, R4
Velocity and acceleration statistics of heavy spheroidal particles in turbulence
Non-spherical particles transported by turbulent flow have a rich dynamics that combines their translational and rotational motions. Here, the focus is on small, heavy, inertial particles with a spheroidal shape fully prescribed by their aspect ratio. Such particles undergo an anisotropic, orientation-dependent viscous drag with the carrier fluid flow whose associated torque is given by the Jeffery equations. Direct numerical simulations of homogeneous, isotropic turbulence are performed to study systematically how the translational motion of such spheroidal particles depends on their shape and size. It is found that the Lagrangian statistics of both velocity and acceleration can be described in terms of an effective Stokes number obtained as an isotropic average over angles of the particle's orientation. Corrections to the translational motion of particles due to their non-sphericity and rotation can hence be recast as an effective radius obtained from such a mean.
- 12 June 2023
- Published by Springer Nature in The European Physical Journal E
- Vol. 46 (6) , 1-22
Steering undulatory micro-swimmers in a fluid flow through reinforcement learning
This work aims at finding optimal navigation policies for thin, deformable microswimmers that progress in a viscous fluid by propagating a sinusoidal undul
- 12 June 2023
- Published by Springer Nature in The European Physical Journal E
- Vol. 46 (6) , 1-22
Steering undulatory micro-swimmers in a fluid flow through reinforcement learning
This work aims at finding optimal navigation policies for thin, deformable microswimmers that progress in a viscous fluid by propagating a sinusoidal undul