Asymptotic analysis of initial flow around an impulsively started circular cylinder using a Brinkman penalization method

The initial flow past an impulsively started translating circular cylinder is asymptotically analysed using a Brinkman penalization method on the Navier–Stokes equations. The asymptotic solution obtained shows that the tangential and normal slip velocities on the cylinder surface are of the order of $1/\sqrt {\lambda }$ and $1/\lambda$ , respectively, within the second approximation of the present asymptotic analysis, where $\lambda$ is the penalization parameter. This result agrees with the estimation of Carbou & Fabrie (Adv. Diff. Equ., vol. 8, 2003, pp. 1453–1480). Based on the asymptotic solution, the influence of the penalization parameter $\lambda$ is discussed on the drag coefficient that is calculated using the adopted three formulae. It can then be found that the drag coefficient calculated from the integration of the penalization term exhibits a half-value of the results of Bar-Lev & Yang (J. Fluid Mech., vol. 72, 1975, pp. 625–647) as $\lambda \to \infty$ .

The Coupled Volume of Fluid and Brinkman Penalization Methods for Simulation of Incompressible Multiphase Flows

In this work, we contribute to the development of numerical algorithms for the direct simulation of three-dimensional incompressible multiphase flows in the presence of multiple fluids and solids. The volume of fluid method is used for interface tracking, and the Brinkman penalization method is used to treat solids; the latter is assumed to be perfectly superhydrophobic or perfectly superhydrophilic, to have an arbitrary shape, and to move with a prescribed velocity. The proposed algorithm is implemented in the open-source software Basilisk and is validated on a number of test cases, such as the Stokes flow between a periodic array of cylinders, vortex decay problem, and multiphase flow around moving solids.

Utilization of the Brinkman Penalization to Represent Geometries in a High-Order Discontinuous Galerkin Scheme on Octree Meshes

We investigate the suitability of the Brinkman penalization method in the context of a high-order discontinuous Galerkin scheme to represent wall boundaries in compressible flow simulations. To evaluate the accuracy of the wall model in the numerical scheme, we use setups with symmetric reflections at the wall. High-order approximations are attractive as they require few degrees of freedom to represent smooth solutions. Low memory requirements are an essential property on modern computing systems with limited memory bandwidth and capability. The high-order discretization is especially useful to represent long traveling waves, due to their small dissipation and dispersion errors. An application where this is important is the direct simulation of aeroacoustic phenomena arising from the fluid motion around obstacles. A significant problem for high-order methods is the proper definition of wall boundary conditions. The description of surfaces needs to match the discretization scheme. One option to achieve a high-order boundary description is to deform elements at the boundary into curved elements. However, creating such curved elements is delicate and prone to numerical instabilities. Immersed boundaries offer an alternative that does not require a modification of the mesh. The Brinkman penalization is such a scheme that allows us to maintain cubical elements and thereby the utilization of efficient numerical algorithms exploiting symmetry properties of the multi-dimensional basis functions. We explain the Brinkman penalization method and its application in our open-source implementation of the discontinuous Galerkin scheme, Ateles. The core of this presentation is the investigation of various penalization parameters. While we investigate the fundamental properties with one-dimensional setups, a two-dimensional reflection of an acoustic pulse at a cylinder shows how the presented method can accurately represent curved walls and maintains the symmetry of the resulting wave patterns.