Global stability analysis and direct numerical simulation of boundary layers with an isolated roughness element

Global stability analysis and direct numerical simulation (DNS) are performed to study boundary layer flows with an isolated roughness element. Wall-attached cuboids with aspect ratios $\eta=1$ and $\eta=0.5$ are investigated for fixed ratio of roughness height to displacement boundary layer thickness $h/\delta^*=2.86$. Global stability analysis is able to capture the frequency of the primary vortical structures. For $\eta=1$, only varicose instability is seen. For the thinner roughness element ($\eta=0.5$), the varicose instability dominates the sinuous instability, and the sinuous instability becomes more pronounced as $Re_h$ increases, due to increased spanwise shear in the near-wake region. The unstable modes mainly extract energy from the central streak, although the lateral streaks also contribute. The DNS results show that different instability features lead to different behavior and development of vortical structures in the nonlinear transition process. For $\eta=1$, the varicose mode is associated with the shedding of hairpin vortices. As $Re_h$ increases, the breakdown of hairpin vortices occurs closer to the roughness and sinuous breakdown behavior promoting transition to turbulence is seen in the farther wake. A fully-developed turbulent flow is established in both the inner and outer layers farther downstream when $Re_h$ is sufficiently high. For $\eta=0.5$, the sinuous wiggling of hairpin vortices is prominent at higher $Re_h$, leading to stronger interactions in the near wake, as a result of combined varicose and sinuous instabilities. A sinuous mode captured by dynamic mode decomposition (DMD) analysis, and associated with the `wiggling' of streaks persists far downstream

A variational volume-of-fluid approach for front propagation

A variational volume-of-fluid (VVOF) methodology is devised for evolving interfaces under curvature-dependent speed. The interface is reconstructed geometrically using the analytic relations of Scardovelli and Zaleski [1] and the advection of the volume fraction is performed using the algorithm of Weymouth and Yue (WY) [2] with a technique to incorporate a volume conservation constraint. The proposed approach has the advantage of simple implementation and straightforward extension to more complex systems. Canonical curves and surfaces traditionally investigated by the level set (LS) method are tested with the VVOF approach and results are compared with existing work in LS