2025-04-03 | | Total: 9
Understanding and solving fluid dynamics equations efficiently remains a fundamental challenge in computational physics. Traditional numerical solvers and physics-informed neural networks struggle to capture the full range of frequency components in partial differential equation solutions, limiting their accuracy and efficiency. Here, we propose the Multi-stream Physics Hybrid Network, a novel neural architecture that integrates quantum and classical layers in parallel to improve the accuracy of solving fluid dynamics equations, namely Kovasznay flow problem. This approach decomposes the solution into separate frequency components, each predicted by independent Parallel Hybrid Networks, simplifying the training process and enhancing performance. We evaluated the proposed model against a comparable classical neural network, the Multi-stream Physics Classical Network, in both data-driven and physics-driven scenarios. Our results show that the Multi-stream Physics Hybrid Network achieves a reduction in root mean square error by 36% for velocity components and 41% for pressure prediction compared to the classical model, while using 24% fewer trainable parameters. These findings highlight the potential of hybrid quantum-classical architectures for advancing computational fluid dynamics.
Microswimming cells and robots exhibit diverse behaviours due to both their swimming and their environment. One of the core environmental features impacting inertialess swimming is background flows. While the influence of select flows, particularly shear flows, have been extensively investigated, these are special cases. Here, we examine inertialess swimmers in more general flows, specifically general linear planar flows that may also possess rapid oscillations. Relatively weak symmetry constraints are imposed on the swimmer to ensure planarity and to reduce complexity. A further constraint reflecting common observation is imposed, namely that the swimmer is inefficient, which we suitably define. This introduces two separate timescales: a fast timescale associated with swimmer actuation, and a second timescale associated with net swimmer movement, with inefficiency dictating that this latter timescale is much slower, allowing for a multiple timescale simplification of the governing equations. With the exception of mathematically precise edge cases, we find that the behaviour of the swimmer is dictated by two parameter groupings, both of which measure balances between the angular velocity and rate of strain of the background flow. While the measures of flow angular velocity and strain rates that primarily govern the rotational dynamics are modulated by swimmer properties, the primary features of the translational motion are determined solely by a ratio of flow angular velocity to strain rate. Hence, a simple classification of the swimmer dynamics emerges. For example, this illustrates the limited extent to which, and how, microswimmers may control their orientations and trajectories in flows.
When do mucus films plug lung airways? Using reduced-order simulations of a large ensemble of randomly perturbed films, we show that the answer is not determined by just the film's volume. While very thin films always stay open and very thick films always plug, we find a range of intermediate films for which plugging is uncertain. The fastest-growing linear mode of the Rayleigh-Plateau instability ensures that the film's volume is divided among multiple humps. However, the nonlinear growth of these humps can occur unevenly, due to spontaneous axial sliding -- a lucky hump can sweep up a disproportionate share of the film's volume and so form a plug. This sliding-induced plugging is robust and prevails with or without gravitational and ciliary transport.
We conduct linear analyses of convection in domains larger than the temperature scale height. We employ both analytical and numerical methods in these analyses. In the case excluding all dissipation, the typical time scale of convection is determined by the free fall time over the temperature scale height. We quantitatively show the condition for the Boussinesq and Wentzel-Kramers-Brillouin (WKB) approximations to be applicable. We provide a reassessment of the critical Rayleigh number, a key indicator of convection, and show that WKB approximation tends to underestimate the critical Rayleigh number, particularly when the temperature scale height is comparable to or smaller than the domain height. We show clear explanation why both thermal conduction and viscosity are required for stabilizing negative entropy gradient medium.
This experimental study aims to investigate the effect of different base blowing configurations on the aerodynamics of a squareback Ahmed body. Four different slot configurations through which air is injected were studied. Each configuration had the same blowing area, equivalent to 10% of the base area of the body, and was designated according to its geometric shape: square (S), vertical (V), cross (C) and horizontal (H). The same range of injected flow rates, C_{q}, was tested for each slot geometry at a Reynolds number Re= 65000, with corresponding wind tunnel measurements. Our experiments revealed the effect of blowing on the near wake and consequently on the base pressure and drag of the Ahmed body. In particular, the geometry of the slots was shown to be a crucial factor in influencing aerodynamics, especially at blowing flow rates close to C_{q,opt}, which is the blowing flow rate that provides the minimum drag coefficient. The centered square slot (S) is the configuration that achieves a greater drag reduction. This behavior is attributed to an elongation in the recirculation region behind the Ahmed body, a reduction in the backflow inside the recirculation bubble, and a decrease in the wake asymmetry associated with the Reflexional Symmetry Breaking (RSB) mode. Conversely, the vertically oriented slot geometries, such as the cross (C) and the vertical (V) configurations, showed limited drag reduction capability while maintaining or even intensifying the wake asymmetry. The horizontal (H) slot represented an intermediate case, mitigating the wake asymmetry to a large extent, but proving to be less effective in reducing the drag than the square case. The hierarchy of the blowing configurations was dictated by the modifications induced in the near wake behind the Ahmed body, which influenced the asymmetry of the wake and the filling/emptying of the recirculation region.
Microscale turbulent flow in porous media is conducive to the development of flow instabilities due to strong vortical and shearing flow occurring within the pore space. When the flow instabilities around individual solid obstacles interact with numerous others within the porous medium, unique symmetry-breaking phenomena emerge as a result. This paper focuses on investigations of the vortex dynamics and flow instabilities behind solid obstacles in porous media, emphasizing how solid obstacle geometry and porosity influence both microscale and macroscale flow behavior. Two distinct symmetry-breaking mechanisms were identified in different porosity ranges. In low porosity media (< 0.8), a "deviatory flow" phenomenon occurs, where the macroscale flow deviates from the direction of applied pressure gradient at Reynolds numbers above 500. Deviatory flow is a source of macroscale Reynolds stress anisotropy, which is counterbalanced by a diminished vortex core size. In the intermediate porosity regime (0.8-0.95), a "jetting flow" mechanism creates asymmetric microscale velocity channels in the pore space through temporally biased vortex shedding, occurring during the transition to turbulence. Both symmetry-breaking phenomena are critically influenced by solid obstacle shape, porosity, and Reynolds number. Circularity of solid obstacle geometry and an adequately high Reynolds number provide critical conditions for symmetry-breaking, whereas porosity can be used to parametrize the degree of symmetry-breaking. This paper provides fundamental insights into the intricate flow dynamics in porous media, offering a comprehensive understanding of how microscale vortex interactions generate macroscale flow asymmetries across different geometric configurations.
We present a novel approach for simulating acoustic (pressure) wave propagation across different media separated by a diffuse interface through the use of a weak compressibility formulation. Our method builds on our previous work on an entropy-stable discontinuous Galerkin spectral element method for the incompressible Navier-Stokes/Cahn-Hilliard system \cite{manzanero2020entropyNSCH}, and incorporates a modified weak compressibility formulation that allows different sound speeds in each phase. We validate our method through numerical experiments, demonstrating spectral convergence for acoustic transmission and reflection coefficients in one dimension and for the angle defined by Snell's law in two dimensions. Special attention is given to quantifying the modeling errors introduced by the width of the diffuse interface. Our results show that the method successfully captures the behavior of acoustic waves across interfaces, allowing exponential convergence in transmitted waves. The transmitted angles in two dimensions are accurately captured for air-water conditions, up to the critical angle of 13^\circ. This work represents a step forward in modeling acoustic propagation in incompressible multiphase systems, with potential applications to marine aeroacoustics.
Small-scale gravity waves play a crucial role in atmospheric tracer transport, yet their effects remain unresolved in climate models and must be parameterized. This study investigates how gravity waves influence large-scale tracer distributions, utilizing a multiple-scale analysis to systematically identify the governing terms of gravity wave-induced tracer fluxes. The analysis reveals both leading-order and next-order impacts: the former being the inertia-gravity wave-induced tracer Stokes drift, which acts perpendicular to both the large-scale tracer gradient and the wave number, while the latter becomes significant at lower latitudes where Coriolis effects diminish. A numerical framework is developed to incorporate these fluxes into a gravity wave parameterization model, potentially enhancing climate model accuracy without requiring explicit resolution of small-scale wave dynamics. Model validation against high-resolution wave-resolving simulations confirms the effectiveness of this approach. By improving the representation of gravity wave-induced tracer transport, this research advances the accuracy of climate simulations, particularly in their depiction of microphysics and radiative processes.
Direct numerical simulation of three-dimensional acoustic turbulence has been performed for both weak and strong regimes. Within the weak turbulence, we demonstrate the existence of the Zakharov-Sagdeev spectrum \propto k^{-3/2} not only for weak dispersion but in the non-dispersion (ND) case as well. Such spectra in the k-space are accompanied by jets in the form of narrow cones. These distributions are realized due to small nonlinearity compared with both dispersion/diffraction. Increasing pumping in the ND case due to dominant nonlinear effects leads to the formation of shocks. As a result, the acoustic turbulence turns into an ensemble of random shocks with the Kadomtsev-Petviashvili spectrum.