Researcher profile

Kunihiko Taira

Kunihiko Taira contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
6topics
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

14 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.

preprint2025arXiv

Phase autoencoder for rapid data-driven synchronization of rhythmic spatiotemporal patterns

We present a machine-learning method for data-driven synchronization of rhythmic spatiotemporal patterns in reaction-diffusion systems. Based on the phase autoencoder [Yawata {\it et al.}, Chaos {\bf 34}, 063111 (2024)], we map high-dimensional field variables of the reaction-diffusion system to low-dimensional latent variables characterizing the asymptotic phase and amplitudes of the field variables. This yields a reduced phase description of the limit cycle underlying the rhythmic spatiotemporal dynamics in a data-driven manner. We propose a method to drive the system along the tangential direction of the limit cycle, enabling phase control without inducing amplitude deviations. With examples of 1D oscillating spots and 2D spiral waves in the FitzHugh-Nagumo reaction-diffusion system, we show that the method achieves rapid synchronization in both reference-based and coupling-based settings. These results demonstrate the potential of data-driven phase description based on the phase autoencoder for synchronization of high-dimensional spatiotemporal dynamics.

preprint2022arXiv

A Database for Reduced-Complexity Modeling of Fluid Flows

We present a publicly accessible database designed to aid in the conception, training, demonstration, evaluation, and comparison of reduced-complexity models for fluid mechanics. Availability of high-quality flow data is essential for all of these aspects of model development for both data-driven and physics-based methods. The database contains time-resolved data for six distinct datasets: a large eddy simulation of a turbulent jet, direct numerical simulations of a zero-pressure-gradient turbulent boundary layer, particle-image-velocimetry measurements for the same boundary layer at several Reynolds numbers, direct numerical simulations of laminar stationary and pitching flat-plate airfoils, particle-image-velocimetry and force measurements of an airfoil encountering a gust, and a large eddy simulation of the separated, turbulent flow over an airfoil. These six cases span several key flow categories: laminar and turbulent, statistically stationary and transient, tonal and broadband spectral content, canonical and application-oriented, wall-bounded and free-shear flow, and simulation and experimental measurements. For each dataset, we describe the flow setup and computational/experimental methods, catalog the data available in the database, and provide examples of how these data can be used for reduced-complexity modeling. All data can be downloaded using a browser interface or Globus. Our vision is that the common testbed provided by this database will aid the fluid mechanics community in clarifying the distinct capabilities of new and existing methods.

preprint2022arXiv

Transition, intermittency and phase interference effects in airfoil secondary tones and acoustic feedback loop

A large eddy simulation is performed to study secondary tones generated by a NACA0012 airfoil at angle of attack of $α= 3^{\circ}$ with freestream Mach number of $M_{\infty} = 0.3$ and Reynolds number of $Re = 5 \times 10^4$. Laminar separation bubbles are observed over the suction side and near the trailing edge, on the pressure side. Flow visualization and spectral analysis are employed to investigate vortex shedding aft of the suction side separation bubble. Vortex interaction results in merging or bursting such that coherent structures or turbulent packets are advected towards the trailing edge leading to different levels of noise emission. Despite the intermittent occurrence of laminar-turbulent transition, the noise spectrum depicts a main tone with multiple equidistant secondary tones. To understand the role of flow instabilities on the tones, the linearized Navier-Stokes equations are examined in its operator form through bi-global stability and resolvent analyses, and by time evolution of disturbances using a matrix-free method. These linear global analyses reveal amplification of disturbances over the suction side separation bubble. Non-normality of the linear operator leads to further transient amplification due to modal interaction among eigenvectors. Two-point, one time autocovariance calculations of pressure along the spanwise direction elucidate aspects of the acoustic feedback loop mechanism in the non-linear solutions. This feedback process is self-sustained by acoustic waves radiated from the trailing edge, which reach the most sensitive flow location between the leading edge and the separation bubble, as identified by the resolvent analysis. Leading edge disturbances arising from secondary diffraction and phase interference among the most unstable frequencies computed in the eigenspectrum are also shown to have an important role in the feedback loop.

preprint2022arXiv

Wing sweep effects on laminar separated flows

We reveal the effects of sweep on the wake dynamics around NACA 0015 wings at high angles of attack using direct numerical simulations and resolvent analysis. The influence of sweep on the wake dynamics is considered for sweep angles from $0^\circ$ to $45^\circ$ and angles of attack from $16^\circ$ to $30^\circ$ for a spanwise periodic wing at a chord-based Reynolds number of $400$ and a Mach number of $0.1$. Wing sweep affects the wake dynamics, especially in terms of stability and spanwise fluctuations with implications on the development of three-dimensional wakes. We observe that wing sweep attenuates spanwise fluctuations. Even as the sweep angle influences the wake, force and pressure coefficients can be collapsed for low angles of attack when examined in wall-normal and wingspan-normal independent flow components. Some small deviations at high sweep and incidence angles are attributed to vortical wake structures that impose secondary aerodynamic loads, revealed through the force element analysis. Furthermore, we conduct global resolvent analysis to uncover oblique modes with high disturbance amplification. The resolvent analysis also reveals the presence of wavemakers in the shear-dominated region associated with the emergence of three-dimensional wakes at high angles of attack. For flows at high sweep angles, the optimal convection speed of the response modes is shown to be faster than the optimal wavemakers speed suggesting a mechanism for the attenuation of perturbations. The present findings serve as a fundamental stepping stone to understanding separated flows at higher Reynolds numbers.

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.

preprint2020arXiv

Active flow control of a pump-induced wall-normal vortex with steady blowing

The emergence of a submerged vortex upstream of a pump can reduce pump intake efficiency and cause structural damage. In this study, we consider the use of active flow control with steady blowing to increase the pressure distribution within a single-phase pump-induced wall-normal vortex model, which is based on the Burgers vortex with a no-slip boundary condition prescribed along its symmetry plane. The goal of our control is to modify the vortex core velocity profile. These changes are sought to increase the core pressure such that detrimental effects on the pump are alleviated. Three-dimensional direct numerical simulations (DNS) are performed to examine the dynamics of the vortex with the application of axial momentum injection at and around the root of the vortex. We find that the active flow control approach can effectively modify the wall-normal vortical structure and significantly increase the low-core pressure by up to 81% compared to that of the uncontrolled case. The result shows that the control setup is also effective when it is introduced in an off-centered manner. Compared to the unsteady blowing and suction based actuation from our previous work (Liu et al. 2018), the current steady control technique offers an effective and simple flow control setup that can support robust operations of pumps.

preprint2020arXiv

Assessment of supervised machine learning methods for fluid flows

We apply supervised machine learning techniques to a number of regression problems in fluid dynamics. Four machine learning architectures are examined in terms of their characteristics, accuracy, computational cost, and robustness for canonical flow problems. We consider the estimation of force coefficients and wakes from a limited number of sensors on the surface for flows over a cylinder and NACA0012 airfoil with a Gurney flap. The influence of the temporal density of the training data is also examined. Furthermore, we consider the use of convolutional neural network in the context of super-resolution analysis of two-dimensional cylinder wake, two-dimensional decaying isotropic turbulence, and three-dimensional turbulent channel flow. In the concluding remarks, we summarize on findings from a range of regression type problems considered herein.

preprint2020arXiv

Machine learning based spatio-temporal super resolution reconstruction of turbulent flows

We present a new turbulent data reconstruction method with supervised machine learning techniques inspired by super resolution and inbetweening, which can recover high-resolution turbulent flows from grossly coarse flow data in space and time. For the present machine learning based data reconstruction, we use the downsampled skip-connection/multi-scale model based on a convolutional neural network to incorporate the multi-scale nature of fluid flows into its network structure. As an initial example, the model is applied to a two-dimensional cylinder wake at $Re_D$ = 100. The reconstructed flow fields by the proposed method show great agreement with the reference data obtained by direct numerical simulation. Next, we examine the capability of the proposed model for a two-dimensional decaying homogeneous isotropic turbulence. The machine-learned models can follow the decaying evolution from coarse input data in space and time, according to the assessment with the turbulence statistics. The proposed concept is further investigated for a complex turbulent channel flow over a three-dimensional domain at $Re_τ$ =180. The present model can reconstruct high-resolved turbulent flows from very coarse input data in space, and it can also reproduce the temporal evolution when the time interval is appropriately chosen. The dependence on the amount of training snapshots and duration between the first and last frames based on a temporal two-point correlation coefficient are also assessed to reveal the capability and robustness of spatio-temporal super resolution reconstruction. These results suggest that the present method can meet a range of flow reconstructions for supporting computational and experimental efforts.

preprint2020arXiv

Network broadcast analysis and control of turbulent flows

We present a network-based modal analysis technique that identifies key dynamical paths along which perturbations amplify over a time-varying base flow. This analysis is built upon the Katz centrality, which reveals the flow structures that can effectively spread perturbations over a time-evolving network of vortical interactions on the base flow. Motivated by the resolvent form of the Katz function, we take the singular value decomposition of the resulting communicability matrix, complementing the resolvent analysis for fluid flows. The right-singular vectors, referred to as the broadcast modes, give insights into the sensitive regions where introduced perturbations can be effectively spread and amplified over the entire fluid-flow network that evolves in time. We apply this analysis to a two-dimensional decaying isotropic turbulence. The broadcast mode reveals that vortex dipoles are important structures in spreading perturbations. By perturbing the flow with the principal broadcast mode, we demonstrate the utility of the insights gained from the present analysis to effectively modify the evolution of turbulent flows. The current network-inspired work presents a novel use of network analysis to guide flow control efforts, in particular for time-varying base flows.

preprint2020arXiv

Randomized resolvent analysis

Performing global resolvent analysis for high-Reynolds-number turbulent flow calls for the handling of a large discrete operator. Even though such large operator is required in the analysis, most applications of resolvent analysis extracts only a few dominant resolvent response and forcing modes. Here, we consider the use of randomized numerical linear algebra to reduce the dimension of the resolvent operator for achieving computational speed up and memory saving compared to the standard resolvent analysis. To accomplish this goal, we utilize sketching of the linear operator with random test matrices with a Gaussian distribution and with insights from the base flow incorporated to perform singular value decomposition on a low-rank matrix holding dominant characteristics of the full resolvent operator. The strength of the randomized resolvent analysis is demonstrated on a turbulent separated flow over an airfoil. This randomized approach clears the path towards tackling resolvent analysis for higher-Reynolds number bi- and tri-global base flows.

preprint2020arXiv

Resolvent analysis of an airfoil laminar separation bubble at Re = 500,000

We perform resolvent analysis to examine the perturbation dynamics over the laminar separation bubble (LSB) that forms near the leading edge of a NACA 0012 airfoil at a chord-based Reynolds number of 500,000 and an angle of attack of 8 degrees. Randomized SVD is adopted in the present analysis to relieve the computational cost associated with the high-Re global base flow. To examine the local physics over the LSB, we consider the use of exponential discounting to limit the time horizon that allows for the instability to develop with respect to the base flow. With discounting, the gain distribution over frequency accurately captures the spectral content over the LSB obtained from flow simulation. The peak-gain frequency also agrees with previous flow control results on suppressing dynamic stall over a pitching airfoil. According to the gain distribution and the modal structures, we conclude that the dominant energy-amplification mechanism is the Kelvin-Helmholtz instability. In addition to discounting, we also examine the use of spatial windows for both the forcing and response. From the response-windowed analysis, we find that the LSB serves the main role of energy amplifier, with the amplification saturating at the reattachment point. The input window imposes the constraint of surface forcing, and the results show that the optimal actuator location is slightly upstream of the separation point. The surface-forcing mode also suggest the optimal momentum forcing in the surface-tangent direction, with strong uni-directionality that is ideal for synthetic-jet-type actuators. This study demonstrates the strength of randomized resolvent analysis in tackling high-Reynolds-number base flows and calls attention to the care needed for base-flow instabilities.

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.

preprint2019arXiv

Randomized methods to characterize large-scale vortical flow network

We demonstrate the effective use of randomized methods for linear algebra to perform network-based analysis of complex vortical flows. Network theoretic approaches can reveal the connectivity structures among a set of vortical elements and analyze their collective dynamics. These approaches have recently been generalized to analyze high-dimensional turbulent flows, for which network computations can become prohibitively expensive. In this work, we propose efficient methods to approximate network quantities, such as the leading eigendecomposition of the adjacency matrix, using randomized methods. Specifically, we use the Nyström method to approximate the leading eigenvalues and eigenvectors, achieving significant computational savings and reduced memory requirements. The effectiveness of the proposed technique is demonstrated on two high-dimensional flow fields: two-dimensional flow past an airfoil and two-dimensional turbulence. We find that quasi-uniform column sampling outperforms uniform column sampling, while both feature the same computational complexity.