Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSpontaneous Symmetry Breaking for Extreme Vorticity and Strain in the 3D Navier-Stokes Equations
96
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate the spatio-temporal structure of the most likely configurations realising extremely high vorticity or strain in the stochastically forced 3D incompressible Navier-Stokes equations. Most likely configurations are computed by numerically finding the highest probability velocity field realising an extreme constraint as solution of a large optimisation problem. High-vorticity configurations are identified as pinched vortex filaments with swirl, while high-strain configurations correspond to counter-rotating vortex rings. We additionally observe that the most likely configurations for vorticity and strain spontaneously break their rotational symmetry for extremely high observable values. Instanton calculus and large deviation theory allow us to show that these maximum likelihood realisations determine the tail probabilities of the observed quantities. In particular, we are able to demonstrate that artificially enforcing rotational symmetry for large strain configurations leads to a severe underestimate of their probability, as it is dominated in likelihood by an exponentially more likely symmetry broken vortex-sheet configuration.
rate research
Read More
Turbulent fluid flows are ubiquitous in nature and technology, and are mathematically described by the incompressible Navier-Stokes equations (INSE). A hallmark of turbulence is spontaneous generation of intense whirls, resulting from amplification of the fluid rotation-rate (vorticity) by its deformation-rate (strain). This interaction, encoded in the non-linearity of INSE, is non-local, i.e., depends on the entire state of the flow, constituting a serious hindrance in turbulence theory and in establishing regularity of INSE. Here, we unveil a novel aspect of this interaction, by separating strain into local and non-local contributions utilizing the Biot-Savart integral of vorticity in a sphere of radius R. Analyzing highly-resolved numerical turbulent solutions to INSE, we find that when vorticity becomes very large, the local strain over small R surprisingly counteracts further amplification. This uncovered self-attenuation mechanism is further shown to be connected to local Beltramization of the flow, and could provide a direction in establishing the regularity of INSE.
Whether the 3D incompressible Navier-Stokes equations can develop a finite time singularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the 3D incompressible axisymmetric Navier-Stokes equations with smooth initial data of finite energy develop nearly singular solutions at the origin. This nearly singular behavior is induced by a potential finite time singularity of the 3D Euler equations that we reported in cite{Hou-euler-2021}. One important feature of the potential Euler singularity is that the solution develops nearly self-similar scaling properties that are compatible with those of the 3D Navier-Stokes equations. We will present numerical evidence that the 3D Navier-Stokes equations develop nearly singular scaling properties with maximum vorticity increased by a factor of $10^7$. Moreover, the nearly self-similar profiles seem to be very stable to the small perturbation of the initial data. However, the 3D Navier-Stokes equations with our initial data do not develop a finite time singularity due to the development of a mild two-scale structure in the late stage, which eventually leads to viscous dominance over vortex stretching. To maintain the balance between the vortex stretching term and the diffusion term, we solve the 3D Navier-Stokes equations with a time-dependent viscosity roughly of order $O(|log(T-t)|^{-3})$ in the late stage. We present strong numerical evidence that the 3D Navier-Stokes equations with such time-dependent viscosity develop a finite time singularity.
This paper presents a low-communication-overhead parallel method for solving the 3D incompressible Navier-Stokes equations. A fully-explicit projection method with second-order space-time accuracy is adopted. Combined with fast Fourier transforms, the parallel diagonal dominant (PDD) algorithm for the tridiagonal system is employed to solve the pressure Poisson equation, differing from its recent applications to compact scheme derivatives computation (Abide et al. 2017) and alternating-direction-implicit method (Moon et al. 2020). The number of all-to-all communications is decreased to only two, in a 2D pencil-like domain decomposition. The resulting MPI/OpenMP hybrid parallel code shows excellent strong scalability up to $10^4$ cores and small wall-clock time per timestep. Numerical simulations of turbulent channel flow at different friction Reynolds numbers ($Re_{tau}$ = 550, 1000, 2000) have been conducted and the statistics are in good agreement with the reference data. The proposed method allows massively simulation of wall turbulence at high Reynolds numbers as well as many other incompressible flows.
In this paper, we present strong numerical evidences that the $3$D incompressible axisymmetric Navier-Stokes equations with degenerate diffusion coefficients and smooth initial data of finite energy develop a potential finite time locally self-similar singularity at the origin. The spatial part of the degenerate diffusion coefficient is a smooth function of $r$ and $z$ independent of the solution and vanishes like $O(r^2)+O(z^2)$ near the origin. This potential singularity is induced by a potential singularity of the $3$D Euler equations. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the property that the center of the traveling wave approaches the origin at a slower rate than the rate of the collapse of the singularity. The driving mechanism for this potential singularity is due to two antisymmetric vortex dipoles that generate a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop an additional small scale characterizing the thickness of the sharp front. On the other hand, the Navier-Stokes equations with a constant diffusion coefficient regularize the two-scale solution structure and do not develop a finite time singularity for the same initial data. The initial condition is designed in such a way that it generates a positive feedback loop that enforces a strong nonlinear alignment of vortex stretching, leading to a stable locally self-similar blowup at the origin. We perform careful resolution study and asymptotic scaling analysis to provide further support of the potential finite time locally self-similar blowup.
A conservative discretization of incompressible Navier-Stokes equations is developed based on discrete exterior calculus (DEC). A distinguishing feature of our method is the use of an algebraic discretization of the interior product operator and a combinatorial discretization of the wedge product. The governing equations are first rewritten using the exterior calculus notation, replacing vector calculus differential operators by the exterior derivative, Hodge star and wedge product operators. The discretization is then carried out by substituting with the corresponding discrete operators based on the DEC framework. Numerical experiments for flows over surfaces reveal a second order accuracy for the developed scheme when using structured-triangular meshes, and first order accuracy for otherwise unstructured meshes. By construction, the method is conservative in that both mass and vorticity are conserved up to machine precision. The relative error in kinetic energy for inviscid flow test cases converges in a second order fashion with both the mesh size and the time step.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
