scholarly journals Towards a finite-time singularity of the Navier–Stokes equations. Part 2. Vortex reconnection and singularity evasion

2019 ◽  
Vol 870 ◽  
Author(s):  
H. K. Moffatt ◽  
Yoshifumi Kimura
Keyword(s):  
Dynamical System ◽  
Reynolds Number ◽  
Finite Time ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Time Singularity ◽  
Navier Stokes Equations ◽  
Reconnection Process ◽  
Vortex Reconnection ◽  
Finite Time Singularity

In Part 1 of this work, we have derived a dynamical system describing the approach to a finite-time singularity of the Navier–Stokes equations. We now supplement this system with an equation describing the process of vortex reconnection at the apex of a pyramid, neglecting core deformation during the reconnection process. On this basis, we compute the maximum vorticity $\unicode[STIX]{x1D714}_{max}$ as a function of vortex Reynolds number $R_{\unicode[STIX]{x1D6E4}}$ in the range $2000\leqslant R_{\unicode[STIX]{x1D6E4}}\leqslant 3400$, and deduce a compatible behaviour $\unicode[STIX]{x1D714}_{max}\sim \unicode[STIX]{x1D714}_{0}\exp [1+220(\log [R_{\unicode[STIX]{x1D6E4}}/2000])^{2}]$ as $R_{\unicode[STIX]{x1D6E4}}\rightarrow \infty$. This may be described as a physical (although not strictly mathematical) singularity, for all $R_{\unicode[STIX]{x1D6E4}}\gtrsim 4000$.

scholarly journals Towards a finite-time singularity of the Navier–Stokes equations Part 1. Derivation and analysis of dynamical system

2018 ◽  
Vol 861 ◽  
pp. 930-967 ◽  
Author(s):  
H. K. Moffatt ◽  
Yoshifumi Kimura
Keyword(s):  
Dynamical System ◽  
Finite Time ◽  
Stokes Equations ◽  
Rate Of Change ◽  
Navier Stokes ◽  
Tipping Points ◽  
Dimensionless Time ◽  
Time Singularity ◽  
Navier Stokes Equations ◽  
Finite Time Singularity

The evolution towards a finite-time singularity of the Navier–Stokes equations for flow of an incompressible fluid of kinematic viscosity$\unicode[STIX]{x1D708}$is studied, starting from a finite-energy configuration of two vortex rings of circulation$\pm \unicode[STIX]{x1D6E4}$and radius$R$, symmetrically placed on two planes at angles$\pm \unicode[STIX]{x1D6FC}$to a plane of symmetry$x=0$. The minimum separation of the vortices,$2s$, and the scale of the core cross-section,$\unicode[STIX]{x1D6FF}$, are supposed to satisfy the initial inequalities$\unicode[STIX]{x1D6FF}\ll s\ll R$, and the vortex Reynolds number$R_{\unicode[STIX]{x1D6E4}}=\unicode[STIX]{x1D6E4}/\unicode[STIX]{x1D708}$is supposed very large. It is argued that in the subsequent evolution, the behaviour near the points of closest approach of the vortices (the ‘tipping points’) is determined solely by the curvature$\unicode[STIX]{x1D705}(\unicode[STIX]{x1D70F})$at the tipping points and by$s(\unicode[STIX]{x1D70F})$and$\unicode[STIX]{x1D6FF}(\unicode[STIX]{x1D70F})$, where$\unicode[STIX]{x1D70F}=(\unicode[STIX]{x1D6E4}/R^{2})t$is a dimensionless time variable. The Biot–Savart law is used to obtain analytical expressions for the rate of change of these three variables, and a nonlinear dynamical system relating them is thereby obtained. The solution shows a finite-time singularity, but the Biot–Savart law breaks down just before this singularity is realised, when$\unicode[STIX]{x1D705}s$and$\unicode[STIX]{x1D6FF}/\!s$become of order unity. The dynamical system admits ‘partial Leray scaling’ of just$s$and$\unicode[STIX]{x1D705}$, and ultimately full Leray scaling of$s,\unicode[STIX]{x1D705}$and$\unicode[STIX]{x1D6FF}$, conditions for which are obtained. The tipping point trajectories are determined; these meet at the singularity point at a finite angle. An alternative model is briefly considered, in which the initial vortices are ovoidal in shape, approximately hyperbolic near the tipping points, for which there is no restriction on the initial value of the parameter$\unicode[STIX]{x1D705}$; however, it is still the circles of curvature at the tipping points that determine the local evolution, so the same dynamical system is obtained, with breakdown again of the Biot–Savart approach just before the incipient singularity is realised. The Euler flow situation ($\unicode[STIX]{x1D708}=0$) is considered, and it is conjectured on the basis of the above dynamical system that a finite-time singularity can indeed occur in this case.

A note on the solution of the Navier-Stokes equations for a spherically symmetric expansion into a very low pressure

1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar
Keyword(s):  
Shock Wave ◽  
Reynolds Number ◽  
Stokes Equation ◽  
Stokes Equations ◽  
Low Pressure ◽  
Navier Stokes ◽  
Spherically Symmetric ◽  
Navier Stokes Equation ◽  
Navier Stokes Equations ◽  
Type Flow

It is shown that, for a spherically symmetric expansion of a gas into a low pressure, the shock wave with area change region discussed earlier (Freeman & Kumar 1972) can be further divided into two parts. For the Navier–Stokes equation, these are a region in which the asymptotic zero-pressure behaviour predicted by Ladyzhenskii is achieved followed further downstream by a transition to subsonic-type flow. The distance of this final region downstream is of order (pressure)−2/3 × (Reynolds number)−1/3.

Free-stream coherent structures in parallel boundary-layer flows

2014 ◽  
Vol 752 ◽  
pp. 602-625 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Coherent Structures ◽  
Free Stream ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Homotopy Continuation ◽  
Boundary Layer Flows ◽  
Navier Stokes Equations ◽  
High Reynolds

AbstractOur concern in this paper is with high-Reynolds-number nonlinear equilibrium solutions of the Navier–Stokes equations for boundary-layer flows. Here we consider the asymptotic suction boundary layer (ASBL) which we take as a prototype parallel boundary layer. Solutions of the equations of motion are obtained using a homotopy continuation from two known types of solutions for plane Couette flow. At high Reynolds numbers, it is shown that the first type of solution takes the form of a vortex–wave interaction (VWI) state, see Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666), and is located in the main part of the boundary layer. On the other hand, here the second type is found to support an equilibrium solution of the unit-Reynolds-number Navier–Stokes equations in a layer located a distance of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}O(\ln \mathit{Re})$ from the wall. Here $\mathit{Re}$ is the Reynolds number based on the free-stream speed and the unperturbed boundary-layer thickness. The streaky field produced by the interaction grows exponentially below the layer and takes its maximum size within the unperturbed boundary layer. The results suggest the possibility of two distinct types of streaky coherent structures existing, possibly simultaneously, in disturbed boundary layers.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

scholarly journals Vortex reconnections in classical and quantum fluids

SeMA Journal ◽  
2021 ◽  
Author(s):  
Alberto Enciso ◽  
Daniel Peralta-Salas
Keyword(s):  
Stokes Equations ◽  
Quantum Fluids ◽  
Navier Stokes ◽  
Pitaevskii Equation ◽  
Global Approximation ◽  
Rigorous Results ◽  
Navier Stokes Equations ◽  
Localization Principle ◽  
Global Smooth Solutions ◽  
Vortex Reconnection

AbstractWe review recent rigorous results on the phenomenon of vortex reconnection in classical and quantum fluids. In the context of the Navier–Stokes equations in $$\mathbb {T}^3$$ T 3 we show the existence of global smooth solutions that exhibit creation and destruction of vortex lines of arbitrarily complicated topologies. Concerning quantum fluids, we prove that for any initial and final configurations of quantum vortices, and any way of transforming one into the other, there is an initial condition whose associated solution to the Gross–Pitaevskii equation realizes this specific vortex reconnection scenario. Key to prove these results is an inverse localization principle for Beltrami fields and a global approximation theorem for the linear Schrödinger equation.

scholarly journals Numerical Simulation of Slip effect on Lid-Driven Cavity Flow Problem for High Reynolds Number: Vorticity – Stream Function Approach

2021 ◽  
Vol 8 (3) ◽  
pp. 418-424
Author(s):  
Syed Fazuruddin ◽  
Seelam Sreekanth ◽  
G. Sankara Sekhar Raju
Keyword(s):  
Reynolds Number ◽  
Stream Function ◽  
Stokes Equations ◽  
Cavity Flow ◽  
Partial Slip ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Driven Cavity ◽  
Slip Conditions ◽  
Driven Cavity Flow

Incompressible 2-D Navier-stokes equations for various values of Reynolds number with and without partial slip conditions are studied numerically. The Lid-Driven cavity (LDC) with uniform driven lid problem is employed with vorticity - Stream function (VSF) approach. The uniform mesh grid is used in finite difference approximation for solving the governing Navier-stokes equations and developed MATLAB code. The numerical method is validated with benchmark results. The present work is focused on the analysis of lid driven cavity flow of incompressible fluid with partial slip conditions (imposed on side walls of the cavity). The fluid flow patterns are studied with wide range of Reynolds number and slip parameters.

Two-dimensional flow of a viscous fluid in a channel with porous walls

1991 ◽  
Vol 227 ◽  
pp. 1-33 ◽  
Author(s):  
Stephen M. Cox
Keyword(s):  
Finite Time ◽  
Stokes Equations ◽  
Pitchfork Bifurcation ◽  
Time Dependent ◽  
High Rate ◽  
Navier Stokes ◽  
Recent Theory ◽  
Navier Stokes Equations ◽  
Two Dimensional Flow ◽  
Steady Solutions

We consider the flow of a viscous incompressible fluid in a parallel-walled channel, driven by steady uniform suction through the porous channel walls. A similarity transformation reduces the Navier-Stokes equations to a single partial differential equation (PDE) for the stream function, with two-point boundary conditions. We discuss the bifurcations of the steady solutions first, and show how a pitchfork bifurcation is unfolded when a symmetry of the problem is broken.Then we describe time-dependent solutions of the governing PDE, which we calculate numerically. We analyse these unsteady solutions when there is a high rate of suction through one wall, and the other wall is impermeable: there is a limit cycle composed of an explosive phase of inviscid growth, and a slow viscous decay. The inviscid phase ‘almost’ has a finite-time singularity. We discuss whether solutions of the governing PDE, which are exact solutions of the Navier-Stokes equations, may develop mathematical singularities in a finite time.When the rates of suction at the two walls are equal so that the problem is symmetrical, there is an abrupt transition to chaos, a ‘homoclinic explosion’, in the time-dependent solutions as the Reynolds number is increased. We unfold this transition by perturbing the symmetry, and compare direct numerical integrations of the governing PDE with a recent theory for ‘Lorenz-like’ dynamical systems. The chaos is found to be very sensitive to symmetry breaking.

An Implicit Multigrid Scheme for the Compressible Navier-Stokes Equations With Low-Reynolds-Number Turbulence Closure

1998 ◽  
Vol 120 (2) ◽  
pp. 257-262 ◽  
Author(s):  
Peter Gerlinger ◽  
Dieter Bru¨ggemann
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Low Reynolds Number ◽  
Convergence Acceleration ◽  
Iterative Solvers ◽  
Turbulence Closure ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Low Reynolds ◽  
Turbulence Transport

A multigrid method for convergence acceleration is used for solving coupled fluid and turbulence transport equations. For turbulence closure a low-Reynolds-number q-ω turbulence model is employed, which requires very fine grids in the near wall regions. Due to the use of fine grids, convergence of most iterative solvers slows down, making the use of multigrid techniques especially attractive. However, special care has to be taken on the strong nonlinear turbulent source terms during restriction from fine to coarse grids. Due to the hyperbolic character of the governing equations in supersonic flows and the occurrence of shock waves, modifications to standard multigrid techniques are necessary. A simple and effective method is presented that enables the multigrid scheme to converge. A strong reduction in the required number of multigrid cycles and work units is achieved for different test cases, including a Mack 2 flow over a backward facing step.

Sign in / Sign up

Export Citation Format

Share Document