Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).
more »
« less
Organizing bacterial vortex lattices by periodic obstacle arrays
Abstract Recent experiments have shown that the complex spatio-temporal vortex structures emerging in active fluids are susceptible to weak geometrical constraints. This observation poses the fundamental question of how boundary effects stabilize a highly ordered pattern from seemingly turbulent motion. Here we show, by a combination of continuum theory and experiments on a bacterial suspension, how artificial obstacles guide the flow profile and reorganize topological defects, which enables the design of bacterial vortex lattices with tunable properties. To this end, the continuum model is extended by appropriate boundary conditions. Beyond the stabilization of square and hexagonal lattices, we also provide a striking example of a chiral, antiferromagnetic lattice exhibiting a net rotational flow, which is induced by arranging the obstacles in a Kagome-like array.
more »
« less
- Award ID(s):
- 1707900
- PAR ID:
- 10214376
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Communications Physics
- Volume:
- 3
- Issue:
- 1
- ISSN:
- 2399-3650
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The present study examines the northern stratosphere during April 2020, when the polar vortex split into two cyclonic vortices during a winter‐early spring period with the strongest ozone depletion on record. We investigate the dynamical evolution leading to the split at middle stratospheric levels, including the fate of fluid parcels on the vortex boundary during its rupture and the distribution of ozone between the vortices resulting from the split. We also illustrate the vertical structure of the vortices after the split. The findings obtained with Lagrangian methods confirm the key role for the split played by a flow with a special configuration of barriers to the motion of parcels. A trajectory analysis clarifies how the ozone distribution between vortices was such that ozone poorest air remained in the main vortex. The offspring vortex had a deep structure from the troposphere and later decayed to vanish by the end of April.more » « less
-
Active turbulence, or chaotic self-organized collective motion, is often observed in concentrated suspensions of motile bacteria and other systems of self-propelled interacting agents. To date, there is no fundamental understanding of how geometrical confinement orchestrates active turbulence and alters its physical properties. Here, by combining large-scale experiments, computer modeling, and analytical theory, we have identified a generic sequence of transitions occurring in bacterial suspensions confined in cylindrical wells of varying radii. With increasing the well’s radius, we observed that persistent vortex motion gives way to periodic vortex reversals, four-vortex pulsations, and then well-developed active turbulence. Using computational modeling and analytical theory, we have shown that vortex reversal results from the nonlinear interaction of the first three azimuthal modes that become unstable with the radius increase. The analytical results account for our key experimental findings. To further validate our approach, we reconstructed equations of motion from experimental data. Our findings shed light on the universal properties of confined bacterial active matter and can be applied to various biological and synthetic active systems.more » « less
-
Vortex breakdown is an instability which occurs in swirling flows and is seen in both incompressible and compressible flow regimes. The choice of numerical algorithm to com- pute this flow is critical to resolving the relevant physical phenomena. In previous work, we described the procedure of combining a Barely Implicit Correction (BIC) algorithm with a fourth-order Flux-Corrected Transport (FCT) algorithm to eliminate the sound speed limit for explicit calculations of low-speed flows, and showed that BIC-FCT could predict the three major vortex breakdown modes (spiral, bubble and double-helix). Here, we show a direct comparison of the implicit BIC and explicit FCT algorithms by using the same time step within the explicit stability limit. The comparison is conducted on simulations of a three-dimensional, swirling jet flow with vortex breakdown, by examining the flow structures and examining the property fields. Differences are observed in the phase of the downstream spiral and double-helix modes. Based on a pressure wave analysis, we attribute this phase shift to the difference of how the explicit FCT and implicit BIC treat an open outflow boundary.more » « less
-
van_Kessel, Julia C (Ed.)ABSTRACT Bacterial motility plays a crucial role in biofilm development, yet the underlying mechanism remains not fully understood. Here, we demonstrate that the flagellum-driven motility ofPseudomonas aeruginosaenhances biofilm formation by altering the orientation of bacterial cells, an effect controlled by shear stress rather than shear rate. By tracking wild-typeP. aeruginosaand its non-motile mutants in a microfluidic channel, we demonstrate that while non-motile cells align with the flow, many motile cells can orient toward the channel sidewalls, enhancing cell surface attachment and increasing biofilm cell density by up to 10-fold. Experiments with varying fluid viscosities further demonstrate that bacterial swimming speed decreases with increasing fluid viscosity, and the cell orientation scales with the shear stress rather than shear rate. Our results provide a quantitative framework to predict the role of motility in the orientation and biofilm development under different flow conditions and viscosities.IMPORTANCEBiofilms are ubiquitous in rivers, water pipes, and medical devices, impacting the environment and human health. While bacterial motility plays a crucial role in biofilm development, a mechanistic understanding remains limited, hindering our ability to predict and control biofilms. Here, we reveal how the motility ofPseudomonas aeruginosa, a common pathogen, influences biofilm formation through systematically controlled microfluidic experiments with confocal and high-speed microscopy. We demonstrate that the orientation of bacterial cells is controlled by shear stress. While non-motile cells primarily align with the flow, many motile cells overcome the fluid shear forces and reorient toward the channel sidewalls, increasing biofilm cell density by up to 10-fold. Our findings provide insights into how bacterial transition from free-swimming to surface-attached states under varying flow conditions, emphasizing the role of cell orientation in biofilm establishment. These results enhance our understanding of bacterial behavior in flow environments, informing strategies for biofilm management and control.more » « less
An official website of the United States government
