More Like this

1. Finite Time Blowup of 2D Boussinesq and 3D Euler Equations with $$C^{1,\alpha}$$ Velocity and Boundary

   https://doi.org/https://doi.org/10.1007/s00220-021-04067-1  

   Chen, Jiajie; Hou, Thomas Y (April 2021, Communications in mathematical physics)

   null (Ed.)

   Inspired by the numerical evidence of a potential 3D Euler singularity by Luo- Hou [30,31] and the recent breakthrough by Elgindi [11] on the singularity formation of the 3D Euler equation without swirl with $$C^{1,\alpha}$$ initial data for the velocity, we prove the finite time singularity for the 2D Boussinesq and the 3D axisymmetric Euler equations in the presence of boundary with $$C^{1,\alpha}$$ initial data for the velocity (and density in the case of Boussinesq equations). Our finite time blowup solution for the 3D Euler equations and the singular solution considered in [30,31] share many essential features, including the symmetry properties of the solution, the flow structure, and the sign of the solution in each quadrant, except that we use $$C^{1,\alpha}$$ initial data for the velocity field. We use a dynamic rescaling formulation and follow the general framework of analysis developed by Elgindi in [11]. We also use some strategy proposed in our recent joint work with Huang in [7] and adopt several methods of analysis in [11] to establish the linear and nonlinear stability of an approximate self-similar profile. The nonlinear stability enables us to prove that the solution of the 3D Euler equations or the 2D Boussinesq equations with $$C^{1,\alpha}$$ initial data will develop a finite time singularity. Moreover, the velocity field has finite energy before the singularity time. 

   more » « less
2. Singularity formation in 3D Euler equations with smooth initial data and boundary

   https://doi.org/10.1073/pnas.2500940122  

   Chen, Jiajie; Hou, Thomas_Y (June 2025, Proceedings of the National Academy of Sciences)

   A long-standing fundamental open problem in mathematical fluid dynamics and nonlinear partial differential equations is to determine whether solutions of the 3D incompressible Euler equations can develop a finite-time singularity from smooth, finite-energy initial data. Leonhard Euler introduced these equations in 1757 [L. Euler,Mémoires de l’Académie des Sci. de Berlin11, 274–315 (1757).], and they are closely linked to the Navier–Stokes equations and turbulence. While the general singularity formation problem remains unresolved, we review a recent computer-assisted proof of finite-time, nearly self-similar blowup for the 2D Boussinesq and 3D axisymmetric Euler equations in a smooth bounded domain with smooth initial data. The proof introduces a framework for (nearly) self-similar blowup, demonstrating the nonlinear stability of an approximate self-similar profile constructed numerically via the dynamical rescaling formulation. 

   more » « less
3. Stability of a Subcritical Fluid Model for Fair Bandwidth Sharing with General File Size Distributions

   https://doi.org/10.1287/stsy.2019.0058  

   Fu, Yingjia; Williams, Ruth J. (September 2020, Stochastic Systems)

   null (Ed.)

   This work concerns the asymptotic behavior of solutions to a (strictly) subcritical fluid model for a data communication network, where file sizes are generally distributed and the network operates under a fair bandwidth-sharing policy. Here we consider fair bandwidth-sharing policies that are a slight generalization of the [Formula: see text]-fair policies introduced by Mo and Walrand [Mo J, Walrand J (2000) Fair end-to-end window-based congestion control. IEEE/ACM Trans. Networks 8(5):556–567.]. Since the year 2000, it has been a standing problem to prove stability of the data communications network model of Massoulié and Roberts [Massoulié L, Roberts J (2000) Bandwidth sharing and admission control for elastic traffic. Telecommunication Systems 15(1):185–201.], with general file sizes and operating under fair bandwidth sharing policies, when the offered load is less than capacity (subcritical conditions). A crucial step in an approach to this problem is to prove stability of subcritical fluid model solutions. In 2012, Paganini et al. [Paganini F, Tang A, Ferragut A, Andrew LLH (2012) Network stability under alpha fair bandwidth allocation with general file size distribution. IEEE Trans. Automatic Control 57(3):579–591.] introduced a Lyapunov function for this purpose and gave an argument, assuming that fluid model solutions are sufficiently smooth in time and space that they are strong solutions of a partial differential equation and assuming that no fluid level on any route touches zero before all route levels reach zero. The aim of the current paper is to prove stability of the subcritical fluid model without these strong assumptions. Starting with a slight generalization of the Lyapunov function proposed by Paganini et al., assuming that each component of the initial state of a measure-valued fluid model solution, as well as the file size distributions, have no atoms and have finite first moments, we prove absolute continuity in time of the composition of the Lyapunov function with any subcritical fluid model solution and describe the associated density. We use this to prove that the Lyapunov function composed with such a subcritical fluid model solution converges to zero as time goes to infinity. This implies that each component of the measure-valued fluid model solution converges vaguely on [Formula: see text] to the zero measure as time goes to infinity. Under the further assumption that the file size distributions have finite pth moments for some p > 1 and that each component of the initial state of the fluid model solution has finite pth moment, it is proved that the fluid model solution reaches the measure with all components equal to the zero measure in finite time and that the time to reach this zero state has a uniform bound for all fluid model solutions having a uniform bound on the initial total mass and the pth moment of each component of the initial state. In contrast to the analysis of Paganini et al., we do not need their strong smoothness assumptions on fluid model solutions and we rigorously treat the realistic, but singular situation, where the fluid level on some routes becomes zero, whereas other route levels remain positive. 

   more » « less
4. Fractional Voigt-Regularization of the 3D Navier–Stokes and Euler Equations: Global Well-Posedness and Limiting Behavior

   https://doi.org/10.1007/s00021-025-00948-w  

   Brzeźniak, Zdzisław; Larios, Adam; Safarik, Isabel (August 2025, Journal of Mathematical Fluid Mechanics)

   The Voigt regularization is a technique used to model turbulent flows, offering advantages such as sharing steady states with the Navier-Stokes equations and requiring no modification of boundary conditions; however, the parabolic dissipative character of the equation is lost. In this work we propose and study a generalization of the Voigt regularization technique by introducing a fractional power $$r$$ in the Helmholtz operator, which allows for dissipation in the system, at least in the viscous case. We examine the resulting fractional Navier-Stokes-Voigt (fNSV) and fractional Euler-Voigt (fEV) and show that global well-posedness holds in the 3D periodic case for fNSV when the fractional power $$r \geq \frac{1}{2}$$ and for fEV when $$r>\frac{5}{6}$$. Moreover, we show that the solutions of these fractional Voigt-regularized systems converge to solutions of the original equations, on the corresponding time interval of existence and uniqueness of the latter, as the regularization parameter $$\alpha \to 0$$. Additionally, we prove convergence of solutions of fNSV to solutions of fEV as the viscosity $$\nu \to 0$$ as well as the convergence of solutions of fNSV to solutions of the 3D Euler equations as both $$\alpha, \nu \to 0$$. Furthermore, we derive a criterion for finite-time blow-up for each system based on this regularization. These results may be of use to researchers in both pure and applied fluid dynamics, particularly in terms of approximate models for turbulence and as tools to investigate potential blow-up of solutions. 

   more » « less
5. Squirmers with swirl at low Weissenberg number

   https://doi.org/10.1017/jfm.2020.987  

   Housiadas, Kostas D.; Binagia, Jeremy P.; Shaqfeh, Eric S.G. (March 2021, Journal of Fluid Mechanics)

   We investigate aspects of the spherical squirmer model employing both large-scale numerical simulations and asymptotic methods when the squirmer is placed in weakly elastic fluids. The fluids are modelled by differential equations, including the upper-convected Maxwell (UCM)/Oldroyd-B, finite-extensibility nonlinear elastic model with Peterlin approximation (FENE-P) and Giesekus models. The squirmer model we examine is characterized by two dimensionless parameters related to the fluid velocity at the surface of the micro-swimmer: the slip parameter $$\xi $$ and the swirl parameter $$\zeta $$ . We show that, for swimming in UCM/Oldroyd-B fluids, the elastic stress becomes singular at a critical Weissenberg number, Wi , that depends only on $$\xi$$ . This singularity for the UCM/Oldroyd-B models is independent of the domain exterior to the swimmer, or any other forces considered in the momentum balance for the fluid – we believe that this is the first time such a singularity has been explicitly demonstrated. Moreover, we show that the behaviour of the solution at the poles is purely extensional in character and is the primary reason for the singularity in the Oldroyd-B model. When the Giesekus or the FENE-P models are utilized, the singularity is removed. We also investigate the mechanism behind the speed and rotation rate enhancement associated with the addition of swirl in the swimmer's gait. We demonstrate that, for all models, the speed is enhanced by swirl, but the mechanism of enhancement depends intrinsically on the rheological model employed. Special attention is paid to the propulsive role of the pressure and elucidated upon. We also study the region of convergence of the perturbation solutions in terms of Wi . When techniques that accelerate the convergence of series are applied, transformed solutions are derived that are in very good agreement with the results obtained by simulations. Finally, both the analytical and numerical results clearly indicate that the low- Wi region is more important than one would expect and demonstrate all the major phenomena observed when swimming with swirl in a viscoelastic fluid. 

   more » « less