BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tim Colonius

Tim Colonius contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynphysics.comp-phmath.NANumerical AnalysisMachine LearningComputational Engineering, Finance, and ScienceNeural and Evolutionary Computingphysics.data-an

Trust snapshot

Quick read

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

preprint2026arXiv

Resolvent4py: a parallel Python package for analysis, model reduction and control of large-scale linear systems

In this paper, we present resolvent4py, a parallel Python package for the analysis, model reduction and control of large-scale linear systems with millions or billions of degrees of freedom. This package provides the user with a friendly Python-like experience (akin to that of well-established libraries such as numpy and scipy), while enabling MPI-based parallelism through mpi4py, petsc4py and slepc4py. In turn, this allows for the development of streamlined and efficient Python code that can be used to solve several problems in fluid mechanics, solid mechanics, graph theory, molecular dynamics and several other fields.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Normality-based analysis of multiscale velocity gradients and energy transfer in direct and large-eddy simulations of isotropic turbulence

Symmetry-based analyses of multiscale velocity gradients highlight that strain self-amplification (SS) and vortex stretching (VS) drive forward energy transfer in turbulent flows. By contrast, a strain-vorticity covariance mechanism produces backscatter that contributes to the bottleneck effect in the subinertial range of the energy cascade. We extend these analyses by using a normality-based decomposition of filtered velocity gradients in forced isotropic turbulence to distinguish contributions from normal straining, pure shearing and rigid rotation at a given scale. Our analysis of direct numerical simulation (DNS) data illuminates the importance of shear layers in the inertial range and (especially) the subinertial range of the cascade. Shear layers contribute significantly to SS and VS and play a dominant role in the backscatter mechanism responsible for the bottleneck effect. Our concurrent analysis of large-eddy simulation (LES) data characterizes how different closure models affect the flow structure and energy transfer throughout the resolved scales. We thoroughly demonstrate that the multiscale flow features produced by a mixed model closely resemble those in a filtered DNS, whereas the features produced by an eddy viscosity model resemble those in an unfiltered DNS at a lower Reynolds number. This analysis helps explain how small-scale shear layers, whose imprint is mitigated upon filtering, amplify the artificial bottleneck effect produced by the eddy viscosity model in the inertial range of the cascade. Altogether, the present results provide a refined interpretation of the flow structures and mechanisms underlying the energy cascade and insight for designing and evaluating LES closure models.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Velocity gradient partitioning in turbulent flows

The velocity gradient tensor can be decomposed into normal straining, pure shearing and rigid rotation tensors, each with distinct symmetry and normality properties. We partition the strength of turbulent velocity gradients based on the relative contributions of these constituents in several canonical flows. These flows include forced isotropic turbulence, turbulent channels and turbulent boundary layers. For forced isotropic turbulence, the partitioning is in excellent agreement with previous results. For wall-bounded turbulence, the partitioning collapses onto the isotropic partitioning far from the wall, where the mean shearing is relatively weak. By contrast, the near-wall partitioning is dominated by shearing. Between these two regimes, the partitioning collapses well at sufficiently high friction Reynolds numbers and its variations in the buffer layer and the log-law region can be reasonably modelled as a function of the mean shearing strength. Altogether, our results highlight the expressivity and broad applicability of the velocity gradient partitioning as advantages for turbulence modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hybrid quadrature moment method for accurate and stable representation of non-Gaussian processes and their dynamics

Solving the population balance equation (PBE) for the dynamics of a dispersed phase coupled to a continuous fluid is expensive. Still, one can reduce the cost by representing the evolving particle density function in terms of its moments. In particular, quadrature-based moment methods (QBMMs) invert these moments with a quadrature rule, approximating the required statistics. QBMMs have been shown to accurately model sprays and soot with a relatively compact set of moments. However, significantly non-Gaussian processes such as bubble dynamics lead to numerical instabilities when extending their moment sets accordingly. We solve this problem by training a recurrent neural network (RNN) that adjusts the QBMM quadrature to evaluate unclosed moments with higher accuracy. The proposed method is tested on a simple model of bubbles oscillating in response to a temporally fluctuating pressure field. The approach decreases model-form error by a factor of 10 when compared to traditional QBMMs. It is both numerically stable and computationally efficient since it does not expand the baseline moment set. Additional quadrature points are also assessed, optimally placed and weighted according to an additional RNN. These points further decrease the error at low cost since the moment set is again unchanged.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterizing viscoelastic materials via ensemble-based data assimilation of bubble collapse observations

Viscoelastic material properties at high strain rates are needed to model many biological and medical systems. Bubble cavitation can induce such strain rates, and the resulting bubble dynamics are sensitive to the material properties. Thus, in principle, these properties can be inferred via measurements of the bubble dynamics. Estrada et al. (2018) demonstrated such bubble-dynamic high-strain-rate rheometry by using least-squares shooting to minimize the difference between simulated and experimental bubble radius histories. We generalize their technique to account for additional uncertainties in the model, initial conditions, and material properties needed to uniquely simulate the bubble dynamics. Ensemble-based data assimilation minimizes the computational expense associated with the bubble cavitation model. We test an ensemble Kalman filter (EnKF), an iterative ensemble Kalman smoother (IEnKS), and a hybrid ensemble-based 4D--Var method (En4D--Var) on synthetic data, assessing their estimations of the viscosity and shear modulus of a Kelvin--Voigt material. Results show that En4D--Var and IEnKS provide better moduli estimates than EnKF. Applying these methods to the experimental data of Estrada et al. (2018) yields similar material property estimates to those they obtained, but provides additional information about uncertainties. In particular, the En4D--Var yields lower viscosity estimates for some experiments, and the dynamic estimators reveal a potential mechanism that is unaccounted for in the model, whereby the viscosity is reduced in some cases due to material damage occurring at bubble collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enhancement of shock-capturing methods via machine learning

In recent years, machine learning has been used to create data-driven solutions to problems for which an algorithmic solution is intractable, as well as fine-tuning existing algorithms. This research applies machine learning to the development of an improved finite-volume method for simulating PDEs with discontinuous solutions. Shock capturing methods make use of nonlinear switching functions that are not guaranteed to be optimal. Because data can be used to learn nonlinear relationships, we train a neural network to improve the results of a fifth-order WENO method. We post-process the outputs of the neural network to guarantee that the method is consistent. The training data consists of the exact mapping between cell averages and interpolated values for a set of integrable functions that represent waveforms we would expect to see while simulating a PDE. We demonstrate our method on linear advection of a discontinuous function, the inviscid Burgers' equation, and the 1-D Euler equations. For the latter, we examine the Shu-Osher model problem for turbulence-shockwave interactions. We find that our method outperforms WENO in simulations where the numerical solution becomes overly diffused due to numerical viscosity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

FiniteNet: A Fully Convolutional LSTM Network Architecture for Time-Dependent Partial Differential Equations

In this work, we present a machine learning approach for reducing the error when numerically solving time-dependent partial differential equations (PDE). We use a fully convolutional LSTM network to exploit the spatiotemporal dynamics of PDEs. The neural network serves to enhance finite-difference and finite-volume methods (FDM/FVM) that are commonly used to solve PDEs, allowing us to maintain guarantees on the order of convergence of our method. We train the network on simulation data, and show that our network can reduce error by a factor of 2 to 3 compared to the baseline algorithms. We demonstrate our method on three PDEs that each feature qualitatively different dynamics. We look at the linear advection equation, which propagates its initial conditions at a constant speed, the inviscid Burgers' equation, which develops shockwaves, and the Kuramoto-Sivashinsky (KS) equation, which is chaotic.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lift Coefficient Estimation for a Rapidly Pitching Airfoil

We develop a method for estimating the instantaneous lift coefficient on a rapidly pitching airfoil that uses a small number of pressure sensors and a measurement of the angle of attack. The approach assimilates four surface pressure measurements with a modified nonlinear state space model (Goman-Khrabrov model) through a Kalman filter. The error of lift coefficient estimates based only on a weighted-sum of the measured pressures are found to be noisy and biased, which leads to inaccurate estimates. The estimate is improved by including the predictive model in an conventional Kalman filter. The Goman-Khrabrov model is shown to be a linear parameter-varying system and can therefore be used in the Kalman filter without the need for linearization. Additional improvement is realized by modifying the algorithm to provide more accurate estimate of the lift coefficient. The improved Kalman filtering approach results in a bias-free lift coefficient estimate that is more precise than either the pressure-based estimate or the Goman-Khrabrov model on their own. The new method will enable performance enhancements in aerodynamic systems whose performance relies on lift.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lift-up, Kelvin-Helmholtz and Orr mechanisms in turbulent jets

Three amplification mechanisms present in turbulent jets, namely lift-up, Kelvin-Helmholtz, and Orr, are characterized via global resolvent analysis and spectral proper orthogonal decomposition (SPOD) over a range of Mach numbers. The lift-up mechanism was recently identified in turbulent jets via local analysis by Nogueira et al. (J. Fluid Mech., vol. 873, 2019, pp. 211-237) at low Strouhal number ($St$) and non-zero azimuthal wavenumbers ($m$). In these limits, a global SPOD analysis of data from high-fidelity simulations reveals streamwise vortices and streaks similar to those found in turbulent wall-bounded flows. These structures are in qualitative agreement with the global resolvent analysis, which shows that they are a response to upstream forcing of streamwise vorticity near the nozzle exit. Analysis of mode shapes, component-wise amplitudes, and sensitivity analysis distinguishes the three mechanisms and the regions of frequency-wavenumber space where each dominates, finding lift-up to be dominant as $St/m \rightarrow 0$. Finally, SPOD and resolvent analyses of localized regions show that the lift-up mechanism is present throughout the jet, with a dominant azimuthal wavenumber inversely proportional to streamwise distance from the nozzle, with streaks of azimuthal wavenumber exceeding five near the nozzle, and wavenumbers one and two most energetic far downstream of the potential core.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-linear input/output analysis: application to boundary layer transition

We extend linear input/output (resolvent) analysis to take into account nonlinear triadic interactions by considering a finite number of harmonics in the frequency domain using the harmonic balance method. Forcing mechanisms that maximize the drag are calculated using a gradient-based ascent algorithm. By including nonlinearity in the analysis, the proposed frequency-domain framework identifies the worst-case disturbances for laminar-turbulent transition. We demonstrate the framework on a flat-plate boundary layer by considering three-dimensional spanwise-periodic perturbations triggered by a few optimal forcing modes of finite amplitude. Two types of volumetric forcing are considered, one corresponding to a single frequency/spanwise-wavenumber pair, and a multi-harmonic where a harmonic frequency and wavenumber are also added. Depending on the forcing strategy, we recover a range of transition scenarios associated with K-type and H-type mechanisms, including oblique and Tollmien-Schlichting waves, streaks and their breakdown. We show that nonlinearity plays a critical role in optimizing growth by combining and redistributing energy between the linear mechanisms and the higher perturbation harmonics. With a very limited range of frequencies and wavenumbers, the calculations appear to reach the early stages of the turbulent regime through the generation and breakdown of hairpin and quasi-streamwise staggered vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A fast multi-resolution lattice Green's function method for elliptic difference equations

We propose a mesh refinement technique for solving elliptic difference equations on unbounded domains based on the fast lattice Green's function (FLGF) method. The FLGF method exploits the regularity of the Cartesian mesh and uses the fast multipole method in conjunction with fast Fourier transforms to yield linear complexity and decrease time-to-solution. We extend this method to a multi-resolution scheme and allow for locally refined Cartesian blocks embedded in the computational domain. Appropriately chosen interpolation and regularization operators retain consistency between the discrete Laplace operator and its inverse on the unbounded domain. Second-order accuracy and linear complexity are maintained, while significantly reducing the number of degrees of freedom and hence the computational cost.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Gaussian moment method and its augmentation via LSTM recurrent neural networks for the statistics of cavitating bubble populations

Phase-averaged dilute bubbly flow models require high-order statistical moments of the bubble population. The method of classes, which directly evolve bins of bubbles in the probability space, are accurate but computationally expensive. Moment-based methods based upon a Gaussian closure present an opportunity to accelerate this approach, particularly when the bubble size distributions are broad (polydisperse). For linear bubble dynamics a Gaussian closure is exact, but for bubbles undergoing large and nonlinear oscillations, it results in a large error from misrepresented higher-order moments. Long short-term memory recurrent neural networks, trained on Monte Carlo truth data, are proposed to improve these model predictions. The networks are used to correct the low-order moment evolution equations and improve prediction of higher-order moments based upon the low-order ones. Results show that the networks can reduce model errors to less than $1\%$ of their unaugmented values.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock-bubble interaction, and gas bubble cavitation. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock-bubble, shock-droplet, and shock-water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas-liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock-bubble-vessel-wall and acoustic-bubble-net interactions are used to demonstrate the full capabilities of MFC.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulation of humpback whale bubble-net feeding models

Humpback whales can generate intricate bubbly regions, called bubble nets, via their blowholes. They appear to exploit these bubble nets for feeding via loud vocalizations. A fully-coupled phase-averaging approach is used to model the flow, bubble dynamics, and corresponding acoustics. A previously hypothesized waveguiding mechanism is assessed for varying acoustic frequencies and net void fractions. Reflections within the bubbly region result in observable waveguiding for only a small range of flow parameters. A configuration of multiple whales surrounding and vocalizing towards an annular bubble net is also analyzed. For a range of flow parameters the bubble net keeps its core region substantially quieter than the exterior. This approach appears more viable, though it relies upon the cooperation of multiple whales. A spiral bubble net configuration that circumvents this requirement is also investigated. The acoustic wave behaviors in the spiral interior vary qualitatively with the vocalization frequency and net void fraction. The competing effects of vocalization guiding and acoustic attenuation are quantified. Low void fraction cases allow low-frequency waves to partially escape the spiral region, with the remaining vocalizations still exciting the net interior. Higher void fraction nets appear preferable, guiding even low-frequency vocalizations while still maintaining a quiet net interior.

0 citations0 reviewsNo code linkedOpen access