Phase autoencoder for limit-cycle oscillators.

We present a phase autoencoder that encodes the asymptotic phase of a limit-cycle oscillator, a fundamental quantity characterizing its synchronization dynamics. This autoencoder is trained in such a way that its latent variables directly represent the asymptotic phase of the oscillator. The trained autoencoder can perform two functions without relying on the mathematical model of the oscillator: first, it can evaluate the asymptotic phase and the phase sensitivity function of the oscillator; second, it can reconstruct the oscillator state on the limit cycle in the original space from the phase value as an input. Using several examples of limit-cycle oscillators, we demonstrate that the asymptotic phase and the phase sensitivity function can be estimated only from time-series data by the trained autoencoder. We also present a simple method for globally synchronizing two oscillators as an application of the trained autoencoder.

Grasping extreme aerodynamics on a low-dimensional manifold.

Modern air vehicles perform a wide range of operations, including transportation, defense, surveillance, and rescue. These aircraft can fly in calm conditions but avoid operations in gusty environments, encountered in urban canyons, over mountainous terrains, and in ship wakes. With extreme weather becoming ever more frequent due to global warming, it is anticipated that aircraft, especially those that are smaller in size, will encounter sizeable atmospheric disturbances and still be expected to achieve stable flight. However, there exists virtually no theoretical fluid-dynamic foundation to describe the influence of extreme vortical gusts on wings. To compound this difficulty, there is a large parameter space for gust-wing interactions. While such interactions are seemingly complex and different for each combination of gust parameters, we show that the fundamental physics behind extreme aerodynamics is far simpler and lower-rank than traditionally expected. We reveal that the nonlinear vortical flow field over time and parameter space can be compressed to only three variables with a lift-augmented autoencoder while holding the essence of the original high-dimensional physics. Extreme aerodynamic flows can be compressed through machine learning into a low-dimensional manifold, which can enable real-time sparse reconstruction, dynamical modeling, and control of extremely unsteady gusty flows. The present findings offer support for the stable flight of next-generation small air vehicles in atmosphere conditions traditionally considered unflyable.

Data-driven time-dependent state estimation for interfacial fluid mechanics in evaporating droplets.

Droplet evaporation plays crucial roles in biodiagnostics, microfabrication, and inkjet printing. Experimentally studying the evolution of a sessile droplet consisting of two or more components needs sophisticated equipment to control the vast parameter space affecting the physical process. On the other hand, the non-axisymmetric nature of the problem, attributed to compositional perturbations, introduces challenges to numerical methods. In this work, droplet evaporation problem is studied from a new perspective. We analyze a sessile methanol droplet evolution through data-driven classification and regression techniques. The models are trained using experimental data of methanol droplet evolution under various environmental humidity levels and substrate temperatures. At higher humidity levels, the interfacial tension and subsequently contact angle increase due to higher water uptake into droplet. Therefore, different regimes of evolution are observed due to adsorption-absorption and possible condensation of water which turns the droplet from a single component into a binary system. In this work, machine learning and data-driven techniques are utilized to estimate the regime of droplet evaporation, the time evolution of droplet base diameter and contact angle, and level of surrounding humidity. Droplet regime is estimated by classification algorithms through point-by-point analysis of droplet profile. Decision tree demonstrates a better performance compared to Naïve Bayes (NB) classifier. Additionally, the level of surrounding humidity, as well as the time evolution of droplet base diameter and contact angle, are estimated by regression algorithms. The estimation results show promising performance for four cases of methanol droplet evolution under conditions unseen by the model, demonstrating the model's capability to capture the complex physics underlying binary droplet evolution.

Randomized methods to characterize large-scale vortical flow networks.

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.

Sparsification of long range force networks for molecular dynamics simulations.

Atomic interactions in solid materials are described using network theory. The tools of network theory focus on understanding the properties of a system based upon the underlying interactions which govern their dynamics. While the full atomistic network is dense, we apply a spectral sparsification technique to construct a sparse interaction network model that reduces the computational complexity while preserving macroscopic conservation properties. This sparse network is compared to a reduced network created using a cut-off radius (threshold method) that is commonly used to speed-up computations while approximating interatomic forces. The approximations used to estimate the total forces on each atom are quantified to assess how local interatomic force errors propagate errors at the global or continuum scale by comparing spectral sparsification to thresholding. In particular, we quantify the performance of the spectral sparsification algorithm for the short-range Lennard-Jones potential and the long-range Coulomb potential. Spectral sparsification of the Lennard-Jones potential yields comparable results to thresholding while spectral sparsification yields improvements when considering a long-range Coulomb potential. The present network-theoretic formulation is implemented on two sample problems: relaxation of atoms near a surface and a tensile test of a solid with a circular hole.

Network community-based model reduction for vortical flows.

A network community-based reduced-order model is developed to capture key interactions among coherent structures in high-dimensional unsteady vortical flows. The present approach is data-inspired and founded on network-theoretic techniques to identify important vortical communities that are comprised of vortical elements that share similar dynamical behavior. The overall interaction-based physics of the high-dimensional flow field is distilled into the vortical community centroids, considerably reducing the system dimension. Taking advantage of these vortical interactions, the proposed methodology is applied to formulate reduced-order models for the inter-community dynamics of vortical flows, and predict lift and drag forces on bodies in wake flows. We demonstrate the capabilities of these models by accurately capturing the macroscopic dynamics of a collection of discrete point vortices, and the complex unsteady aerodynamic forces on a circular cylinder and an airfoil with a Gurney flap. The present formulation is found to be robust against simulated experimental noise and turbulence due to its integrating nature of the system reduction.

Networked-oscillator-based modeling and control of unsteady wake flows.

A networked-oscillator-based analysis is performed to examine and control the transfer of kinetic energy for periodic bluff body flows. The dynamics of energy fluctuations in the flow field are described by a set of oscillators defined by conjugate pairs of spatial proper orthogonal decomposition (POD) modes. To extract the network of interactions among oscillators, impulse responses of the oscillators to amplitude and phase perturbations are tracked. Tracking small energy inputs and using linear regression, a networked-oscillator model is constructed that reveals energy exchange among the modes. The model captures the nonlinear interactions among the modal oscillators through a linear approximation. A large collection of system responses is aggregated to capture the general network structure of oscillator interactions. The present networked-oscillator model describes the modal perturbation dynamics more accurately than the empirical Galerkin reduced-order model. The linear network model for nonlinear dynamics is subsequently utilized to design a model-based feedback controller. The controller suppresses the modal amplitudes that result in wake unsteadiness leading to drag reduction. The strength of the proposed approach is demonstrated for a canonical example of two-dimensional unsteady flow over a circular cylinder. The present formulation enables the characterization of modal interactions to control fundamental energy transfers in unsteady bluff body flows.