No Arabic abstract
In fluid mechanics, a lot of authors have been reporting analytical solutions of Euler and Navier-Stokes equations. But there is an essential deficiency of non-stationary solutions indeed. In our presentation, we explore the case of non-stationary flows of the Euler equations for incompressible fluids, which should conserve the Bernoulli-function to be invariant for the aforementioned system. We use previously suggested ansatz for solving of the system of Navier-Stokes equations (which is proved to have the analytical way to present its solution in case of conserving the Bernoulli-function to be invariant for such the type of the flows). Conditions for the existence of exact solution of the aforementioned type for the Euler equations are obtained. The restrictions at choosing of the form of the 3D nonstationary solution for the given constant Bernoulli-function B are considered. We should especially note that pressure field should be calculated from the given constant Bernoulli-function B, if all the components of velocity field are obtained.
Incompressible 3D Euler equations develop high vorticity in very thin pancake-like regions from generic large-scale initial conditions. In this work we propose an exact solution of the Euler equations for the asymptotic pancake evolution. This solution combines a shear flow aligned with an asymmetric straining flow, and is characterized by a single asymmetry parameter and an arbitrary transversal vorticity profile. The analysis is based on detailed comparison with numerical simulations performed using a pseudo-spectral method in anisotropic grids of up to 972 x 2048 x 4096.
The incompressible three-dimensional ideal flows develop very thin pancake-like regions of increasing vorticity. These regions evolve with the scaling $omega_{max}(t)proptoell(t)^{-2/3}$ between the vorticity maximum and pancake thickness, and provide the leading contribution to the energy spectrum, where the gradual formation of the Kolmogorov interval $E_{k}propto k^{-5/3}$ is observed for some initial flows [Agafontsev et. al, Phys. Fluids 27, 085102 (2015)]. With the massive numerical simulations, in the present paper we study the influence of initial conditions on the processes of pancake formation and the Kolmogorov energy spectrum development.
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.
We review the continuous symmetry approach and apply it to find the solution, via the construction of constants of motion and infinitesimal symmetries, of the 3D Euler fluid equations in several instances of interest, without recourse to Noethers theorem. We show that the vorticity field is a symmetry of the flow and therefore one can construct a Lie algebra of symmetries if the flow admits another symmetry. For steady Euler flows this leads directly to the distinction of (non-)Beltrami flows: an example is given where the topology of the spatial manifold determines whether the flow admits extra symmetries. Next, we study the stagnation-point-type exact solution of the 3D Euler fluid equations introduced by Gibbon et al. (Physica D, vol.132, 1999, pp.497-510) along with a one-parameter generalisation of it introduced by Mulungye et al. (J. Fluid Mech., vol.771, 2015, pp.468-502). Applying the symmetry approach to these models allows for the explicit integration of the fields along pathlines, revealing a fine structure of blowup for the vorticity, its stretching rate, and the back-to-labels map, depending on the value of the free parameter and on the initial conditions. Finally, we produce explicit blowup exponents and prefactors for a generic type of initial conditions.
The interplay between incompressibility and stratification can lead to non-conservation of horizontal momentum in the dynamics of a stably stratified incompressible Euler fluid filling an infinite horizontal channel between rigid upper and lower plates. Lack of conservation occurs even though in this configuration only vertical external forces act on the system. This apparent paradox was seemingly first noticed by Benjamin (J. Fluid Mech., vol. 165, 1986, pp. 445-474) in his classification of the invariants by symmetry groups with the Hamiltonian structure of the Euler equations in two dimensional settings, but it appears to have been largely ignored since. By working directly with the motion equations, the paradox is shown here to be a consequence of the rigid lid constraint coupling through incompressibility with the infinite inertia of the far ends of the channel, assumed to be at rest in hydrostatic equilibrium. Accordingly, when inertia is removed by eliminating the stratification, or, remarkably, by using the Boussinesq approximation of uniform density for the inertia terms, horizontal momentum conservation is recovered. This interplay between constraints,action at a distance by incompressibility, and inertia is illustrated by layer-averaged exact results, two-layer long-wave models, and direct numerical simulations of the incompressible Euler equations with smooth stratification.