JC

J. Casacuberta Puig

info

Please Note

15 records found

Journal article (2026) - Koen J. Groot, Jordi Casacuberta, Stefan Hickel
A detailed derivation, analysis, and verification is given for the non-orthogonal, plane-marching Parabolized Stability Equations (PSE) approach. In applying the approach to a flow distorted by a medium-amplitude crossflow vortex, we determine its linear secondary instability mechanisms. We show that converged solutions can be achieved for a broad frequency range with an existing stabilization method for the line-marching PSE approach. We verify that 1) solutions converge versus grid size in all dimensions, 2) primary disturbance solutions agree with line-marching PSE results, and 3) secondary disturbance solutions match amplitude and growth-rate evolution of reference Direct Numerical Simulation (DNS) results. We show how and why the type-II instability displays a delayed neutral point when modeled with the plane-marching approach versus the considered local stability approaches, whether the streamwise evolution of the distorted base flow is accounted for or not. This may explain why the type-II disturbance is scarcely captured by DNS in the literature. ...
This work explores the use of a shallow surface hump for passive control and stabilisation of stationary crossflow (CF) instabilities. Wind tunnel experiments are conducted on a spanwise-invariant swept-wing model. The influence of the hump on the boundary layer stability and laminar-turbulent transition is assessed through infrared thermography and particle image velocimetry measurements. The results reveal a strong dependence of the stabilisation effect on the amplitude of the incoming CF disturbances, which is conditioned via discrete roughness elements at the wing leading edge. At a high forcing amplitude, weakly nonlinear stationary CF vortices interact with the hump and result in an abrupt anticipation of transition, essentially tripping the flow. In contrast, at a lower forcing amplitude, CF vortices interact with the hump during linear growth. Notable stabilisation of the primary CF disturbance and considerable transition delay with respect to the reference case (i.e. without hump) is then observed. The spatial region just downstream of the hump apex is shown to be key to the stabilisation mechanism. In this region, the primary CF disturbances rapidly change spanwise orientation and shape, possibly driven by the pressure gradient change-over caused by the hump and the development of CF reversal. The amplitude and shape deformation of the primary CF instabilities are found to contribute to a long-lasting suboptimal growth downstream of the hump, eventually leading to transition delay. ...
Journal article (2025) - Jordi Casacuberta, Sven Westerbeek, Juan Alberto Franco, Koen J. Groot, Stefan Hickel, Stefan Hein, Marios Kotsonis
Stationary velocity-perturbation streaks have recently been identified in laminar swept-wing boundary-layer flow interacting with a surface forward-facing step. Streaky structures at the step promote early laminar-turbulent transition under certain conditions. This work utilizes direct numerical simulations to explore the mechanisms of growth of stationary streaks at the step and provides insight into their origin, nature, and spatial organization. The analysis is mainly focused on, but not restricted to, incoming perturbations in the form of stationary crossflow instability. Stationary streaky structures are found to be universal to swept forward-facing-step flow subjected to three-dimensional perturbations in the incoming boundary layer. The streaks at the step are primarily ascribed to the lift-up effect. They appear as a linear perturbation response of the highly sheared step flow to the cross-stream pattern of incoming perturbations. A mechanism of base-flow deceleration additionally contributes to feeding growth to the streaks. Linear stability analysis carried out through the harmonic Navier-Stokes method confirms that the streaks are a linear perturbation phenomenon. Effects of spanwise perturbation wavelength and effective sweep angle on the mechanisms of the streaks are also assessed. ...

Mechanisms of interaction between crossflow instabilities and forward-facing steps

Doctoral thesis (2025) - J. Casacuberta Puig, M. Kotsonis, S. Hickel
This thesis presents a theoretical and numerical investigation of how a surface forward-facing step alters the stability and transition mechanisms of a laminar incompressible swept-wing boundary layer dominated by stationary crossflow instabilities (CFI). The results elucidate the processes that cause significant transition advancement due to the step and, in turn, challenge the classic paradigm in fluid mechanics that rapid, localised surface-geometry variations are universally detrimental to laminar flow. Through Direct Numerical Simulations (DNS) and new theoretical modelling, this thesis shows that, under specific conditions, a forwardfacing step stabilises pre-existing stationary crossflow vortices and thereby delays laminar-turbulent transition... ...
Book chapter (2024) - J. Casacuberta, S. Hickel, M. Kotsonis
Achieving and maintaining laminar flow on large swept lifting surfaces of subsonic aircraft poses a considerable challenge. ...
Journal article (2024) - Jordi Casacuberta, Stefan Hickel, Marios Kotsonis
A novel mechanism is identified, through which a spanwise-invariant surface feature (a two-dimensional forward-facing step) significantly stabilizes the stationary crossflow instability of a three-dimensional boundary layer. The mechanism is termed here as reverse lift-up effect, inasmuch as it acts reversely to the classic lift-up effect; that is, kinetic energy of an already-existing shear-flow instability is transferred to the underlying laminar flow through the action of cross-stream perturbations. To characterize corresponding energy-transfer mechanisms, a theoretical framework is presented, which is applicable to generic three-dimensional flows and surface features of arbitrary shape with one invariant spatial direction. The identification of a passive geometry-induced effect responsible for dampening stationary crossflow vortices is a promising finding for laminar flow control applications. ...
Although several previous studies have reported a potential drag-reducing effect of dimpled surfaces in turbulent boundary layers, there is a lack of replicability across experiments performed by different research groups. To contribute to the dialogue, we scrutinize one of the most studied dimple geometries reported in the literature, which has a dimple diameter of 20 mm and a depth of 0.5 mm. There is no general consensus in literature on the drag-reduction performance of this particular dimple geometry, with some studies suggesting a drag reduction, while others report a drag increase. The present combined experimental and numerical study comprises two sets of wind tunnel experiments and a well-resolved large-eddy simulation. The wind tunnel experiments and the large-eddy simulation both depict a total drag increase of around 1%–2% compared to the flat reference case. This finding agrees with a recent study by Spalart et al. (2019). Furthermore, the present wind tunnel experiments have shed light on a plausible reason behind the discrepancy between the study by Spalart et al. (2019) and earlier results from van Nesselrooij et al. (2016). Lastly, the large-eddy simulation results reveal that the pressure drag is the main contributor to the increase in the total drag of the dimpled surface. We believe that these results will contribute to a new consensus on the drag performance of this dimple geometry. ...
Conference paper (2022) - J. Casacuberta Puig, K.J. Groot, S. Hickel, M. Kotsonis
The evolution of secondary instabilities in a three-dimensional stationary-crossflow-domina- ted boundary layer is investigated by means of Direct Numerical Simulations (DNS) and linear spanwise BiGlobal stability analysis. Single-frequency unsteady disturbances and a critical stationary crossflow mode are considered. Unsteady perturbation content at 1 kHz manifests in the form of the type-III secondary instability mechanism in the lower portion of the boundary layer in the both the DNS and the stability approach. Considering disturbances at 6 kHz, the results from the stability analysis reveal the existence of largely amplified type-I and type-II secondary instability mechanisms. Strong growth displayed by the former is measured in the DNS, which potentially overshadows manifestations of the type-II mechanism. Laminar- turbulent transition primarily induced by the growth of type-I disturbances is captured in the 6 kHz case. Overall, we report good agreement between DNS and stability analysis in terms of perturbation organization and growth rate for all cases studied. ...
Journal article (2022) - J. Casacuberta, S. Hickel, S. Westerbeek, M. Kotsonis
The interaction between forward-facing steps of several heights and a pre-existing critical stationary crossflow instability of a swept-wing boundary layer is analysed. Direct numerical simulations (DNS) are performed of the incompressible three-dimensional laminar base flow and the stationary distorted flow that arise from the interaction between an imposed primary stationary crossflow perturbation and the steps. These DNS are complemented with solutions of the linear and the nonlinear parabolised stability equations, used towards identifying the influence of linearity and non-parallelism near the step. A fully stationary solution of the Navier-Stokes equations is enforced numerically, in order to isolate the mechanisms pertaining to the interaction of the stationary disturbance with the step. Results provide insight into the salient modifications of the base laminar boundary layer due to the step, and the response of the incoming crossflow instability to these changes. The fundamental spanwise Fourier mode of the disturbance field gradually lifts up as it approaches the step and passes over it. The flow environment around the step is characterised by a sudden spanwise modulation of the base-flow streamlines. Additional stationary perturbation structures are induced at the step, which manifest in the form of spanwise-aligned velocity streaks near the wall. Shortly downstream of the step, the fundamental component of the crossflow perturbation maintains a rather constant amplification for the smallest steps studied. For the largest step, however, the fundamental crossflow perturbation is stabilised significantly shortly downstream of the largest step. This surprising result is ascribed to a modulation of the kinetic energy transfer between the base flow and the fundamental perturbation field, which is brought forward as a new step interaction mechanism. Possible non-modal growth effects at the step are discussed. Furthermore, the results from DNS indicate significant amplification of the high-order harmonic crossflow components downstream of the step. ...
Conference paper (2021) - J. Casacuberta Puig, S. Hickel, M. Kotsonis
We study the interaction between a stationary crossflow instability and forward-facing steps in a swept-wing boundary layer using Direct Numerical Simulations (DNS). The station- ary primary crossflow mode is imposed at the inflow. Steps of several heights are modeled. Particular emphasis is placed on ensuring a fully stationary solution, in order to isolate the modulation of the primary instability at the step and the generation of stationary secondary perturbations. The main features of the base flow are step-induced flow reversal and secondary inflection points. When approaching the step, the incoming crossflow instability lifts up and passes over it. Additional perturbation streaks arise, which are accompanied by secondary stationary vortices that have the same spanwise wavenumber as the fundamental crossflow disturbance. For sufficiently high steps, secondary perturbation structures attain amplitude values comparable in magnitude to primary instability in the vicinity of the step. We propose metrics to quantify the impact of the step on the development of the crossflow disturbance upon interaction. Considering energy-based criteria, we find that the crossflow perturbation is amplified upstream of the step. Interestingly, the incoming primary crossflow instability undergoes strong stabilization immediately downstream of the largest step. The mechanisms revealed in this work provide a first insight into the possible causes of transition on swept wings due to forward-facing steps. ...
Journal article (2019) - J. Casacuberta, K. J. Groot, Q. Ye, S. Hickel
Micro-ramps are popular passive flow control devices which can delay flow separation by re-energising the lower portion of the boundary layer. We compute the laminar base flow, the instantaneous transitional flow, and the mean flow around a micro-ramp immersed in a quasi-incompressible boundary layer at supercritical roughness Reynolds number. Results of our Direct Numerical Simulations (DNS) are compared with results of BiLocal stability analysis on the DNS base flow and independent tomographic Particle Image Velocimetry (tomo-PIV) experiments. We analyse relevant flow structures developing in the micro-ramp wake and assess their role in the micro-ramp functionality, i.e., in increasing the near-wall momentum. The main flow feature of the base flow is a pair of streamwise counter-rotating vortices induced by the micro-ramp, the so-called primary vortex pair. In the instantaneous transitional flow, the primary vortex pair breaks up into large-scale hairpin vortices, which arise due to linear varicose instability of the base flow, and unsteady secondary vortices develop. Instantaneous vortical structures obtained by DNS and experiments are in good agreement. Matching linear disturbance growth rates from DNS and linear stability analysis are obtained until eight micro-ramp heights downstream of the micro-ramp. For the setup considered in this article, we show that the working principle of the micro-ramp is different from that of classical vortex generators; we find that transitional perturbations are more efficient in increasing the near-wall momentum in the mean flow than the laminar primary vortices in the base flow. ...
Micro-ramps are deployed to prevent boundary layer separation by creating a momentum excess close to the wall. Through Direct Numerical Simulations (DNS) of the base, instantaneous and mean flow, we identify that the perturbation dynamics in the wake of the microramp play an essential role in creating the near-wall momentum excess. To identify the origin of the perturbations, we deploy BiGlobal stability analysis on the laminar base flow. We demonstrate that the amplification of the most unstable linear mode is closely related to the time-averaged amplitude of the unsteady perturbations. The flow structure corresponding to this mode has a varicose symmetry with respect to the symmetry plane and matches with the early development of the hairpin vortices in the instantaneous flow field. It is concluded that the varicose instability supported by the laminar base flow represents the mechanism that generates the hairpins. ...
Journal article (2018) - Jordi Casacuberta, Koen J. Groot, Henry J. Tol, Stefan Hickel
Selective Frequency Damping (SFD) is a popular method for the computation of globally unstable steady-state solutions in fluid dynamics. The approach has two model parameters whose selection is generally unclear. In this article, a detailed analysis of the influence of these parameters is presented, answering several open questions with regard to the effectiveness, optimum efficiency and limitations of the method. In particular, we show that SFD is always capable of stabilising a globally unstable systems ruled by one unsteady unstable eigenmode and derive analytical formulas for optimum parameter values. We show that the numerical feasibility of the approach depends on the complex phase angle of the most unstable eigenvalue. A numerical technique for characterising the pertinent eigenmodes is presented. In combination with analytical expressions, this technique allows finding optimal parameters that minimise the spectral radius of a simulation, without having to perform an independent stability analysis. An extension to multiple unstable eigenmodes is derived. As computational example, a two-dimensional cylinder flow case is optimally stabilised using this method. We provide a physical interpretation of the stabilisation mechanism based on, but not limited to, this Navier–Stokes example. ...
We computed the base, instantaneous and mean flow around a micro-ramp immersed in an incompressible boundary layer. Results of our Direct Numerical Simulations (DNS) are compared with an independent stability analysis and experiments. We analyse flow structures and mechanisms that contribute to the micro-ramp functionality and find that transitional perturbations increase the near-wall momentum in the mean flow more efficiently than the primary vortices in the laminar base flow. ...