Abstract:SUMMARYThis paper presents a numerical study of a two-dimensional time-dependent ow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coe cient at the cylinder and for the pressure di erence between the front and the back of the cylinder at the ÿnal time. In addition, the accuracy of these values obtained with di erent time stepping schemes and di erent ÿnite element methods is studied.
“…In [18], this boundary condition is derived in section 3 and appears two lines below (15). Equation (4) of [18] is a discretized version of (32) with q n+1 replaced by ∆t φ n+1 .…”
“…The computational mesh is shown in Figure 7 and the contour plot of the stream function at t= [2,4,5,6,7,8] is shown in Figure 11. For comparison with [15] we also calculate the drag and lift coefficients, denoted by c d (t) and c l (t), which are the x and y components of the quantity…”
Section: Benchmark Tests With Finite Elementsmentioning
How to properly specify boundary conditions for pressure is a longstanding problem for the incompressible Navier-Stokes equations with no-slip boundary conditions. An analytical resolution of this issue stems from a recently developed formula for the pressure in terms of the commutator of the Laplacian and Leray projection operators. Here we make use of this formula to (a) improve the accuracy of computing pressure in two kinds of existing time-discrete projection methods implicit in viscosity only, and (b) devise new higher-order accurate time-discrete projection methods that extend a slip-correction idea behind the well-known finite-difference scheme of Kim and Moin. We test these schemes for stability and accuracy using various combinations of C 0 finite elements. For all three kinds of time discretization, one can obtain 3rd-order accuracy for both pressure and velocity without a time-step stability restriction of diffusive type. Furthermore, two kinds of projection methods are found stable using piecewise linear elements for both velocity and pressure.
“…In [18], this boundary condition is derived in section 3 and appears two lines below (15). Equation (4) of [18] is a discretized version of (32) with q n+1 replaced by ∆t φ n+1 .…”
“…The computational mesh is shown in Figure 7 and the contour plot of the stream function at t= [2,4,5,6,7,8] is shown in Figure 11. For comparison with [15] we also calculate the drag and lift coefficients, denoted by c d (t) and c l (t), which are the x and y components of the quantity…”
Section: Benchmark Tests With Finite Elementsmentioning
How to properly specify boundary conditions for pressure is a longstanding problem for the incompressible Navier-Stokes equations with no-slip boundary conditions. An analytical resolution of this issue stems from a recently developed formula for the pressure in terms of the commutator of the Laplacian and Leray projection operators. Here we make use of this formula to (a) improve the accuracy of computing pressure in two kinds of existing time-discrete projection methods implicit in viscosity only, and (b) devise new higher-order accurate time-discrete projection methods that extend a slip-correction idea behind the well-known finite-difference scheme of Kim and Moin. We test these schemes for stability and accuracy using various combinations of C 0 finite elements. For all three kinds of time discretization, one can obtain 3rd-order accuracy for both pressure and velocity without a time-step stability restriction of diffusive type. Furthermore, two kinds of projection methods are found stable using piecewise linear elements for both velocity and pressure.
“…We will use the Q 2 finite element for the velocity and the P disc 1 (discontinuous linears) finite element for the pressure. This pair of finite element spaces has been proven to be among the best performing ones for discretizing the incompressible Navier-Stokes equations [6,9,10,13]. Since laminar flows will be considered in the numerical studies, a stabilization of the spatial discretization of the Navier-Stokes equations, [1], or the application of a turbulence model, [26], is not necessary.…”
Section: The Navier-stokes Equationsmentioning
confidence: 99%
“…In numerical studies for incompressible flows, the Crank-Nicolson scheme has shown a good relation of accuracy to efficiency [15]. In particular, it was considerably more accurate than the backward Euler scheme (θ 1 = θ 4 = 1, θ 2 = θ 3 = 0), see also [9]. Next, the system (12), (13) is linearized by a fixed point iteration: Given (u…”
Population balance systems are models for processes in nature and industry which lead to a coupled system of equations (Navier-Stokes equations, transport equations, . . .) where the equations are defined in domains with different dimensions. This paper will study the impact of using different schemes for solving the three-dimensional equation of a precipitation process in a two-dimensional flow domain. The numerical schemes for the three-dimensional equation are assessed with respect to the median of the volume fraction of the particle size distribution and the computational costs. It turns out that in the case of a structured flow field with small variations in time all schemes give qualitatively the same results. For a highly time-dependent flow field, the evolution of the median of the volume fraction differs considerably between first order and higher order schemes.
“…The unsteady flow around a cylinder problem [33] is governed by the unsteady Navier-Stokes equations with Re = 100 and f = 0. The problem is given with homogeneous initial value u(x, y, 0) = 0 and boundary conditions at inflow (x = 0) and outflow (x = 2.2) boundaries as u(0, y, t) = u(2.2, y, t) = 0.41 −2 sin( t/8)(6y(0.41− y), 0), 0 y 0.41 (24) No-slip boundary conditions are described at the other boundaries (y = 0, 0.41 and around the cylinder).…”
Section: Unsteady Flow Around a Cylindermentioning
SUMMARYWe consider the Galerkin finite element method for the incompressible Navier-Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite-dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual-free bubble functions. To find the bubble part of the solution, a two-level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier-Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.