Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userParallel-in-Space-and-Time Simulation of the Three-Dimensional, Unsteady Navier-Stokes Equations for Incompressible Flow
65
0
0.0
(
0
)
Added by
Daniel Ruprecht
Publication date
2017
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper we combine the Parareal parallel-in-time method together with spatial parallelization and investigate this space-time parallel scheme by means of solving the three-dimensional incompressible Navier-Stokes equations. Parallelization of time stepping provides a new direction of parallelization and allows to employ additional cores to further speed up simulations after spatial parallelization has saturated. We report on numerical experiments performed on a Cray XE6, simulating a driven cavity flow with and without obstacles. Distributed memory parallelization is used in both space and time, featuring up to 2,048 cores in total. It is confirmed that the space-time-parallel method can provide speedup beyond the saturation of the spatial parallelization.
rate research
Read More
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 introduce and analyze a space-time least-squares method associated to the unsteady Navier-Stokes system. Weak solution in the two dimensional case and regular solution in the three dimensional case are considered. From any initial guess, we construct a minimizing sequence for the least-squares functional which converges strongly to a solution of the Navier-Stokes system. After a finite number of iterates related to the value of the viscosity constant, the convergence is quadratic. Numerical experiments within the two dimensional case support our analysis. This globally convergent least-squares approach is related to the damped Newton method when used to solve the Navier-Stokes system through a variational formulation.
We show that in bounded domains with no-slip boundary conditions, the Navier-Stokes pressure can be determined in a such way that it is strictly dominated by viscosity. As a consequence, in a general domain we can treat the Navier-Stokes equations as a perturbed vector diffusion equation, instead of as a perturbed Stokes system. To illustrate the advantages of this view, we provide a simple proof of the unconditional stability of a difference scheme that is implicit only in viscosity and explicit in both pressure and convection terms, requiring no solution of stationary Stokes systems or inf-sup conditions.
We accomplish two major tasks. First, we show that the turbulent motion at large scales obeys Gaussian statistics in the interval 0 < Rlambda < 8.8, where Rlambda is the microscale Reynolds number, and that the Gaussian flow breaks down to yield place to anomalous scaling at the universal Reynolds number bounding the inequality above. In the inertial range of turbulence that emerges following the breakdown, the effective Reynolds number based on the turbulent viscosity, Rlambda* assumes this same constant value of about 9. This scenario works also for the emergence of turbulence from an initially non-turbulent state. Second, we derive expressions for the anomalous scaling exponents of structure functions and moments of spatial derivatives, by analyzing the Navier-Stokes equations in the form developed by Hopf. We present a novel procedure to close the Hopf equation, resulting in expressions for zetan in the entire range of allowable moment-order, n, and demonstrate that accounting for the temporal dynamics changes the scaling from normal to anomalous. For large n, the theory predicts the saturation of zetan with n, leading to two inferences: (a) the smallest length scale etan = LRe-1 << LRe-3/4, where Re is the large-scale Reynolds number, and (b) velocity excursions across even the smallest length scales can sometimes be as large as the large scale velocity itself. Theoretical predictions for each of these aspects are shown to be in quantitative agreement with available experimental and numerical data.
The approximation of the value function associated to a stabilization problem formulated as optimal control problem for the Navier-Stokes equations in dimension three by means of solutions to generalized Lyapunov equations is proposed and analyzed. The specificity, that the value function is not differentiable on the state space must be overcome. For this purpose a new class of generalized Lyapunov equations is introduced. Existence of unique solutions to these equations is demonstrated. They provide the basis for feedback operators, which approximate the value function, the optimal states and controls, up to arbitrary order.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
