Researcher profile

Vassilios Theofilis

Vassilios Theofilis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2025arXiv

Biglobal resolvent analysis of separated flow over a NACA0012 airfoil

The effects of Reynolds number across $Re=1000$, $2500$, $5000$, and $10000$ on separated flow over a two-dimensional NACA0012 airfoil at an angle of attack of $α=14^\circ$ are investigated through the biglobal resolvent analysis. We identify modal structures and energy amplifications over a range of frequency, spanwise wavenumber, and discount parameter, providing insights across various timescales. Using temporal discounting, we find that the shear layer dynamics dominates over short time horizons, while the wake dynamics becomes the primary amplification mechanism over long time horizons. Spanwise effects also appear over long time horizon, sustained by low frequencies. At a fixed timescale, we investigate the influence of Reynolds number on response and forcing mode structures, as well as the energy gain over different frequencies. Across all Reynolds numbers, the response modes shift from wake-dominated structures at low frequencies to shear layer-dominated structures at higher frequencies. The frequency at which the dominant mechanism changes is independent of the Reynolds number. The response mode structures show similarities across different Reynolds numbers, with local streamwise wavelengths only depending on frequency. Comparisons at a different angle of attack ($α=9^\circ$) show that the transition from wake to shear layer dynamics with increasing frequency only occurs if the unsteady flow is three-dimensional. We also study the dominant frequencies associated with wake and shear layer dynamics across the angles of attack and Reynolds numbers, and present the characteristic scaling for each mechanism.

preprint2021arXiv

Analytical Prediction of Low-Frequency Fluctuations Inside a One-Dimensional Shock

Linear instability of high-speed boundary layers is routinely examined assuming quiescent edge conditions, without reference to the internal structure of shocks or to instabilities potentially generated in them. Our recent work has shown that the kinetically modeled internal nonequilibrium zone of straight shocks away from solid boundaries exhibits low-frequency molecular fluctuations. The presence of the dominant low frequencies observed using the Direct Simulation Monte Carlo (DSMC) method has been explained as a consequence of the well-known bimodal probability density function (PDF) of the energy of particles inside a shock. Here, PDFs of particle energies are derived in the upstream and downstream equilibrium regions, as well as inside shocks, and it is shown for the first time that they have the form of the non-central chi-squared (NCCS) distributions. A linear correlation is proposed to relate the change in the shape of the analytical PDFs as a function of Mach number, within the range $3 \le M \le 10$, with the DSMC-derived average characteristic low-frequency of shocks, as computed in our earlier work. At a given Mach number $M=7.2$, varying the input translational temperature in the range $89 \le T_{tr,1}/(K) \le 1420$, it is shown that the variation in DSMC-derived low-frequencies is correlated with the change in most-probable-speed inside shocks at the location of maximum bulk velocity gradient. Using the proposed linear functions, average low-frequencies are estimated within the examined ranges of Mach number and input temperature and a semi-empirical relationship is derived to predict low-frequency oscillations in shocks. Our model can be used to provide realistic physics-based boundary conditions in receptivity and linear stability analysis studies of laminar-turbulent transition in high-speed flows.

preprint2021arXiv

Kinetic modelling of three-dimensional shock/laminar separation bubble instabilities in hypersonic flows over a double wedge

Linear global instability of the three-dimensional (3-D), spanwise-homogeneous laminar separation bubble (LSB) induced by shock-wave/boundary-layer interaction (SBLI) in a Mach 7 flow of nitrogen over a $30^{\circ}-55^{\circ}$ double wedge is studied. At these conditions corresponding to a freestream unit Reynolds number, $Re_1=5.2\times 10^{4}$ m$^{-1}$, the flow exhibits rarefaction effects and comparable shock-thicknesses to the size of the boundary-layer at separation. This, in turn, requires the use of the high-fidelity Direct Simulation Monte Carlo (DSMC) method to accurately resolve unsteady flow features. We show for the first time that the LSB sustains self-excited, small-amplitude, 3-D perturbations that lead to spanwise-periodic flow structures not only in and downstream of the separated region, as seen in a multitude of experiments and numerical simulations, but also in the internal structure of the separation and detached shock layers. The spanwise-periodicity length and growth rate of the structures in the two zones are found to be identical. It is shown that the linear global instability leads to low-frequency unsteadiness of the triple point formed by the intersection of separation and detached shocks, corresponding to a Strouhal number of $St\sim0.02$. Linear superposition of the spanwise-homogeneous base flow and the leading 3-D flow eigenmode provides further evidence of the strong coupling between linear instability in the LSB and the shock layer.

preprint2021arXiv

Linear modal instabilities around post-stall swept finite-aspect ratio wings at low Reynolds numbers

Linear modal instabilities of flow over finite-span untapered wings have been investigated numerically at Reynolds number 400, at a range of angles of attack and sweep on two wings having aspect ratios 4 and 8. Base flows have been generated by direct numerical simulation, marching the unsteady incompressible three-dimensional Navier-Stokes equations to a steady state, or using selective frequency damping to obtain stationary linearly unstable flows. Unstable three-dimensional linear global modes of swept wings have been identified for the first time using spectral-element time-stepping solvers. The effect of the wing geometry and flow parameters on these modes has been examined in detail. An increase of the angle of attack was found to destabilize the flow, while an increase of the sweep angle had the opposite effect. On unswept wings, TriGlobal analysis revealed that the most unstable global mode peaks in the midspan region of the wake; the peak of the mode structure moves towards the tip as sweep is increased. Data-driven analysis was then employed to study the effects of wing geometry and flow conditions on the nonlinear wake. On unswept wings, the dominant mode at low angles of attack is a Kelvin-Helmholtz-like instability, qualitatively analogous with global modes of infinite-span wings under same conditions. At higher angles of attack and moderate sweep angles, the dominant mode is a structure denominated the interaction mode. At high sweep angles, this mode evolves into elongated streamwise vortices on higher aspect ratio wings, while on shorter wings it becomes indistinguishable from tip-vortex instability.

preprint2019arXiv

On the formation of three-dimensional flows over finite-aspect-ratio wings under tip effects

We perform DNS of flow over finite-aspect-ratio NACA 0015 wings to characterize the tip effects on the wake dynamics. This study focuses on the development of three-dimensional separated flow over the wing, and discuss flow structures formed on the wing surface as well as in the far field wake. Vorticity is introduced into the flow in a predominantly two-dimensional manner. The vortex sheet from the wing tip rolls up around the free end to form the tip vortex. At its inception, the tip vortex is weak and its effect is spatially confined. As the flow around the tip separates, the tip effects extend farther in the spanwise direction, generating noticeable three dimensionality in the wake. For low-aspect-ratio wings, the wake remains stable due to the strong downwash over the entire span. Increasing the aspect ratio allows unsteady vortical flow to emerge away from the tip at sufficiently high angles of attack. These unsteady vortices shed and form closed vortical loops. For higher-aspect-ratio wings, the tip effects retard the near-tip shedding process, which desynchronizes from the two-dimensional shedding over the mid-span region, yielding vortex dislocation. At high angles of attack, the tip vortex exhibits noticeable undulations due to the strong interaction with the unsteady shedding vortices. The spanwise distribution of force coefficients is related to the three-dimensional wake dynamics and the tip effects. Vortical elements in the wake that are responsible for the generation of lift and drag forces are identified through the force element analysis. We note that at high angles of attack, a stationary vortical structure forms at the leading edge near the tip, giving rise to locally high lift and drag forces. The analysis performed here reveals how the vortical flow around the tip influences the separation physics, the global wake dynamics, and the spanwise force distributions.