Iterative Brinkman penalization for remeshed vortex methods
Introduction
Direct Numerical simulations of bluff body flows using remeshed vortex methods have been traditionally associated with the computation of the vorticity flux, via boundary element methods (BEM), for the enforcement of the no-slip boundary condition [1]. In this context, the no-slip boundary condition is enforced in a fractional step algorithm. First, the no-through flow boundary condition is enforced by solving the potential flow problem using a vortex sheet, which is equivalent to a tangential slip velocity at the solid surface. The vortex sheet is then translated into a vorticity flux in order to enforce the no-slip boundary condition [2].
The particle velocity can be computed by solving the Poisson equation [1], [3], [4] using the Fast Multipole [5] or mesh based method (as in vortex-in-cell algorithms [6]). A combination of the boundary element method and the vortex-in-cell [6] method was presented in [7] using a particle–particle particle–mesh () algorithm [8] to correct the vortex-in-cell method with the vorticity obtained from the BEM.
The distortion of the particle locations can render vortex methods inaccurate. In order to remedy this problem, remeshing was introduced [1], [9] to regularize the particle locations. In remeshed vortex methods particles are advected with the velocity of the flow field. When particle locations get distorted their strength is interpolated onto mesh nodes, which become the particles to be advected in the next iteration. The Poisson equation as well as the viscous terms can be computed during remeshing [10] or after the remeshing step, using finite difference operators.
In remeshing, particle–mesh interpolations are computationally efficient by employing Cartesian grids [9] and interpolation functions that are tensorial products. The drawback of using a Cartesian mesh is that during remeshing spurious vorticity can be introduced in the interior of the body. A remedy is to use locally one-sided interpolation or body fitted grids [1], [11], [12]. However both of these approaches have limited applicability. In order to overcome this difficulty for complex, deforming geometries, in recent years a number of efforts have combined the Brinkman penalization techniques [13], [14] with remeshed vortex methods [15], [16], [17], [18], [19]. These simulations have demonstrated the capability of the penalization method to produce results in good agreement with benchmark simulations. A major drawback of these methods however, is that they impose a stringent time step on advancing the flow field to accurately capture vorticity generation at the boundary.
Here we propose a remedy to this situation by using an iterative scheme that ensures the accurate generation of vorticity at the boundary within each time step of the penalty method. We demonstrate the improvement in accuracy over the classical method, using a number of challenging benchmark problems, such as flow past an impulsively started cylinder and normal to a thin flat plate.
The paper is structured as follows: We introduce the governing equations and numerical discretization in Sections 2 Governing equations and numerical method, 3 Brinkman penalization of incompressible viscous flow, 4 An iterative Brinkman penalization method for incompressible viscous flow, followed by simulations (Sections 5 Impulsively started flow past a circular cylinder at, 6 Impulsively started flow normal to a flat plate of finite thickness at, 7 Flow past an impulsively started plate, inclined at 45° and, 8 Uniformly accelerated flow normal to a flat plate of finite thickness at an acceleration-formulated Reynolds number of) of the flow past an impulsively started cylinder at , the impulsively started flow normal to a flat plate at , the impulsively started flow at a 45 degree angle to a flat plate at , and the flow normal to a uniformly accelerated flat plate with an acceleration-based . We conclude by discussing the advantages of the present method over the classical penalization techniques.
Section snippets
Governing equations and numerical method
The velocity-vorticity form of the two-dimensional Navier–Stokes equations can be expressed in a Lagrangian frame as: where ν is the kinematic viscosity, and the vorticity ω is defined as the curl of the velocity field . The velocity field u can be represented by a Helmholtz decomposition where ψ is the stream function and ϕ is the velocity potential. For an incompressible flow () the stream function can be obtained by solving a Poisson equation:
We
Brinkman penalization of incompressible viscous flow
The penalization method enforces the no-slip boundary condition on the surface of a body in an incompressible flow [13] by introducing a source term localized around the surface of the body. The velocity of the flow u is modified by the penalization term as: where v denotes the velocity of the body and denotes the velocity field of the flow prior to penalization i.e. from the solution of Eq. (3). The penalization parameter λ has units of reciprocal time and is equivalent to a
An iterative Brinkman penalization method for incompressible viscous flow
We remark that while the velocity field is related to the vorticity through the elliptic Poisson equation, the penalization introduces only a local relation between velocity and vorticity. This inconsistency can be attributed to the non-consistent handling of the flow kinematics and the no-through flow boundary condition by the classical Brinkman penalization technique.
Here we propose to compute the vorticity field from the penalty term and use it in turn to compute an updated velocity field,
Impulsively started flow past a circular cylinder at
In this section we compare the new iterative penalization with the classical non-iterative explicit penalization method for the benchmark case of the impulsively started flow past a circular cylinder at . The Reynolds number, the non-dimensional time, and the drag coefficient for this case are defined by where D is the diameter of the cylinder. In order to compare with the BEM based vortex method of Koumoutsakos and Leonard [1] the simulation is performed
Impulsively started flow normal to a flat plate of finite thickness at
As mentioned earlier, the penalization method alone does not enforce a global correction of the flow field. A consequence of this is an inherent delay in producing proper surface vorticity to enforce the no-through boundary condition. Hence the flow normal to a flat plate presents a particularly difficult test case for the penalization method. For the impulsively started flow normal to a flat plate the Reynolds number, non-dimensional time, non-dimensional plate thickness, and the drag
Flow past an impulsively started plate, inclined at 45° and
In order to validate the method on an inherently asymmetrical flow, we study the flow past an impulsively started plate at 45° with a thickness and . The non-dimensionalized parameters are the same as those stated in Section 6 cf. Eq. (25). In the present test case, the spatial resolution , corresponding to 8 cell lengths across the plate thickness, and the time step is .
We observe in Fig. 11 that the penalization method does not conserve the zeroth moment of
Uniformly accelerated flow normal to a flat plate of finite thickness at an acceleration-formulated Reynolds number of
In this section we investigate a uniformly accelerated flow normal to a flat plate of finite thickness. The acceleration normalized Reynolds number β, the non-dimensional time, and the drag coefficient as given for this case by Koumoutsakos and Shiels [11] are: where a is the constant acceleration of the plate.
Following [11] we chose such that , and and to reach an acceleration-formulated Reynolds number of . The time step used
Conclusion
We demonstrate the need to extend the Brinkman penalization technique, by performing sub-iterations so as to enhance its accuracy in enforcing the no-through as well as the no-slip flow boundary conditions.
We have coupled this iterative Brinkman penalization technique with remeshed vortex methods. We demonstrate the improved accuracy of our method using benchmark problems computed by remeshed vortex methods with and without penalization. The calculated flow induced forces are in excellent
Acknowledgement
We would like to acknowledge the helpful discussions with Henrik Juul Spietz, Johannes Tophøj Rasmussen, Babak Hejazialhosseini and Wim van Rees. We are grateful for the lift and drag data provided by Jeff D. Eldredge. Finally, we thank the reviewers for their insightful and helpful comments.
References (30)
- et al.
Boundary conditions for viscous vortex methods
J. Comput. Phys.
(1994) - et al.
Vortex methods for direct numerical simulation of three-dimensional bluff body flows: applications to the sphere at , 500 and 1000
J. Comput. Phys.
(2002) - et al.
A fast algorithm for particle simulations
J. Comput. Phys.
(1987) - et al.
An immersed interface method for the vortex-in-cell algorithm
Comput. Struct.
(2007) An influence matrix particle–particle–mesh algorithm with exact particle–particle correction
J. Comput. Phys.
(2003)Inviscid axisymmetrization of an elliptical vortex ring
J. Comput. Phys.
(1997)- et al.
Modified interpolation kernels for treating diffusion and remeshing in vortex methods
J. Comput. Phys.
(2006) - et al.
Vortex methods with spatially varying cores
J. Comput. Phys.
(2000) - et al.
Computation of turbulent flow past an array of cylinders using a spectral method with Brinkman penalization
Eur. J. Mech. B, Fluids
(2001) - et al.
A vortex level set method for the two-way coupling of an incompressible fluid with colliding rigid bodies
J. Comput. Phys.
(2008)
GPU accelerated simulations of bluff body flows using vortex particle methods
J. Comput. Phys.
A multiresolution remeshed vortex-in-cell algorithm using patches
J. Comput. Phys.
Simulations of single and multiple swimmers with non-divergence free deforming geometries
J. Comput. Phys.
A high order solver for the unbounded Poisson equation
J. Comput. Phys.
A Fourier-based elliptic solver for vortical flows with periodic and unbounded directions
J. Comput. Phys.
Cited by (25)
An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies
2022, Journal of Computational PhysicsCitation Excerpt :Specifically, we focus on the vorticity evolution equation governing the 2D incompressible Navier-Stokes equation, written here in conservative form: Integrating this differential conservation law over a 2D region R leads to the integral form In this section, we discretize (6) with a conservative finite difference scheme, using numerical fluxes as described in Shu [41]. While Lagrangian advection schemes are common in vorticity-based discretizations of the 2D Navier-Stokes equations [5,7,42,43], including other immersed interface methods [39,35], we have chosen a conservative Eulerian transport scheme because it offers an explicit notion of the flux of a conserved quantity through a surface. This leads to a discrete form of the integral conservation law (7), which is essential to the discretization of Kelvin's theorem presented in section 5.1.
A remeshed vortex method for mixed rigid/soft body fluid–structure interaction
2021, Journal of Computational PhysicsA velocity-vorticity method for highly viscous 3D flows with application to digital rock physics
2021, Journal of Computational PhysicsImprovement of remeshed Lagrangian methods for the simulation of dissolution processes at pore-scale
2020, Advances in Water ResourcesA review on fluid dynamics of flapping foils
2020, Ocean EngineeringApplication of remeshed vortex method for the simulation of tip vortex at high Reynolds number
2019, Aerospace Science and TechnologyCitation Excerpt :While the Lagrangian advection is retained, the evaluation of variables on the mesh enhances the efficiency. The RVM has been used to simulate the vorticity around 2D and 3D bodies as well as the evolution of vortex structures [10–13]. Encouraging results were obtained and its suitability for vorticity simulation was verified.