ترغب بنشر مسار تعليمي؟ اضغط هنا

Maximum Amplification of Enstrophy in 3D Navier-Stokes Flows

59   0   0.0 ( 0 )
 نشر من قبل Bartosz Protas
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

This investigation concerns a systematic search for potentially singular behavior in 3D Navier-Stokes flows. Enstrophy serves as a convenient indicator of the regularity of solutions to the Navier Stokes system --- as long as this quantity remains finite, the solutions are guaranteed to be smooth and satisfy the equations in the classical (pointwise) sense. However, there are no estimates available with finite a priori bounds on the growth of enstrophy and hence the regularity problem for the 3D Navier-Stokes system remains open. In order to quantify the maximum possible growth of enstrophy, we consider a family of PDE optimization problems in which initial conditions with prescribed enstrophy $mathcal{E}_0$ are sought such that the enstrophy in the resulting Navier-Stokes flow is maximized at some time $T$. Such problems are solved computationally using a large-scale adjoint-based gradient approach derived in the continuous setting. By solving these problems for a broad range of values of $mathcal{E}_0$ and $T$, we demonstrate that the maximum growth of enstrophy is in fact finite and scales in proportion to $mathcal{E}_0^{3/2}$ as $mathcal{E}_0$ becomes large. Thus, in such worst-case scenario the enstrophy still remains bounded for all times and there is no evidence for formation of singularity in finite time. We also analyze properties of the Navier-Stokes flows leading to the extreme enstrophy values and show that this behavior is realized by a series of vortex reconnection events.

قيم البحث

اقرأ أيضاً

The goal of this study is to analyze the fine structure of nonlinear modal interactions in different 1D Burgers and 3D Navier-Stokes flows. This analysis is focused on preferential alignments characterizing the phases of Fourier modes participating i n triadic interactions, which are key to determining the nature of energy fluxes between different scales. We develop novel diagnostic tools designed to probe the level of coherence among triadic interactions realizing different flow scenarios. We consider extreme 1D viscous Burgers flows and 3D Navier-Stokes flows which are complemented by singularity-forming inviscid Burgers flows as well as viscous Burgers flows and Navier-Stokes flows corresponding to generic turbulent and simple unimodal initial data, such as the Taylor-Green vortex. The main finding is that while the extreme viscous Burgers and Navier-Stokes flows reveal the same relative level of enstrophy amplification by nonlinear effects, this behaviour is realized via modal interactions with vastly different levels of coherence. In the viscous Burgers flows the flux-carrying triads have phase values which saturate the nonlinearity thereby maximizing the energy flux towards small scales. On the other hand, in 3D Navier-Stokes flows with the extreme initial data the energy flux to small scales is realized by a very small subset of helical triads. The second main finding concerns the role of initial coherence. Comparison of the flows resulting from the extreme and generic initial conditions shows striking similarities between these two types of flows, for the 1D viscous Burgers equation as well as the 3D Navier-Stokes equation.
172 - Thomas Y. Hou 2021
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 incompre ssible 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.
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.
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, th e 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.
We introduce a model of interacting singularities of Navier-Stokes, named pin,cons. They follow a Hamiltonian dynamics, obtained by the condition that the velocity field around these singularities obeys locally Navier-Stokes equations. This model can be seen of a generalization of the vorton model of Novikov, that was derived for the Euler equations. When immersed in a regular field, the pin,cons are further transported and sheared by the regular field, while applying a stress onto the regular field, that becomes dominant at a scale that is smaller than the Kolmogorov length. We apply this model to compute the motion of a dipole of pin,cons. When the initial relative orientation of the dipole is inside the interval (0, pi/2), a dipole made of pin,con of same intensity exhibits a transient collapse stage, following a scaling with dipole radius tending to 0 like (tc - t) power 0.63. For long time, the dynamics of the dipole is however repulsive, with both components running away from each other to infinity.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا