Source author record

Dmitry Kolomenskiy

Dmitry Kolomenskiy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

12 published item(s)

preprint2022arXiv

An experimental data-driven mass-spring model of flexible Calliphora wings

Insect wings can undergo significant deformation during flapping motion owing to inertial, elastic and aerodynamic forces. Changes in shape then alter aerodynamic forces, resulting in a fully coupled Fluid-Structure Interaction (FSI) problem. Here, we present detailed three-dimensional FSI simulations of deformable blowfly (Calliphora vomitoria) wings in flapping flight. A wing model is proposed using a multi-parameter mass-spring approach, chosen for its implementation simplicity and computational efficiency. We train the model to reproduce static elasticity measurements by optimizing its parameters using a genetic algorithm with covariance matrix adaptation (CMA-ES). Wing models trained with experimental data are then coupled to a high-performance flow solver run on massively parallel supercomputers. Different features of the modeling approach and the intra-species variability of elastic properties are discussed. We found that individuals with different wing stiffness exhibit similar aerodynamic properties characterized by dimensionless forces and power at the same Reynolds number. We further study the influence of wing flexibility by comparing between the flexible wings and their rigid counterparts. Under equal prescribed kinematic conditions for rigid and flexible wings, wing flexibility improves lift-to-drag ratio as well as lift-to-power ratio and reduces peak force observed during wing rotation.

preprint2020arXiv

A mass-spring fluid-structure interaction solver: Application to flexible revolving wings

The secret to the spectacular flight capabilities of flapping insects lies in their wings, which are often approximated as flat, rigid plates. Real wings are however delicate structures, composed of veins and membranes, and can undergo significant deformation. In the present work, we present detailed numerical simulations of such deformable wings. Our results are obtained with a fluid-structure interaction solver, coupling a mass-spring model for the flexible wing with a pseudo-spectral code solving the incompressible Navier-Stokes equations. We impose the no-slip boundary condition through the volume penalization method; the time-dependent complex geometry is then completely described by a mask function. This allows solving the governing equations of the fluid on a regular Cartesian grid. Our implementation for massively parallel computers allows us to perform high resolution computations with up to 500 million grid points. The mass-spring model uses a functional approach, thus modeling the different mechanical behaviors of the veins and the membranes of the wing. We perform a series of numerical simulations of a flexible revolving bumblebee wing at a Reynolds number Re = 1800. In order to assess the influence of wing flexibility on the aerodynamics, we vary the elasticity parameters and study rigid, flexible and highly flexible wing models. Code validation is carried out by computing classical benchmarks.

preprint2020arXiv

Fluid-structure interaction using volume penalization and mass-spring models with application to flapping bumblebee flight

Wing flexibility plays an essential role in the aerodynamic performance of insects due to the considerable deformation of their wings during flight under the impact of inertial and aerodynamic forces. These forces come from the complex wing kinematics of insects. In this study, both wing structural dynamics and flapping wing motion are taken into account to investigate the effect of wing deformation on the aerodynamic efficiency of a bumblebee in tethered flight. A fluid-structure interaction solver, coupling a mass-spring model for the flexible wing with a pseudo-spectral code solving the incompressible Navier-Stokes equations, is implemented for this purpose. We first consider a tethered bumblebee flying in laminar flow with flexible wings. Compared to the rigid model, flexible wings generate smaller aerodynamic forces but require much less power. Finally, the bumblebee model is put into a turbulent flow to investigate its influence on the force production of flexible wings.

preprint2019arXiv

WaveRange: Wavelet-based data compression for three-dimensional numerical simulations on regular grids

A wavelet-based method for compression of three-dimensional simulation data is presented and its software framework is described. It uses wavelet decomposition and subsequent range coding with quantization suitable for floating-point data. The effectiveness of this method is demonstrated by applying it to example numerical tests, ranging from idealized configurations to realistic global-scale simulations.

preprint2016arXiv

Added costs of insect-scale flapping flight in unsteady airflows

The aerial environment in the operating domain of small-scale natural and artificial flapping wing fliers is highly complex, unsteady and generally turbulent. Considering flapping flight in an unsteady wind environment with a periodically varying lateral velocity component, we show that body rotations experienced by flapping wing fliers result in the reorientation of the aerodynamic force vector that can render a substantial cumulative deficit in the vertical force. We derive quantitative estimates of the body roll amplitude and the related energetic requirements to maintain the weight support in free flight under such conditions. We conduct force measurements of a miniature hummingbird-inspired robotic flapper and numerical simulations of a bumblebee. In both cases, we demonstrate the loss of weight support due to body roll rotations. Using semi-restrained flight measurements, we demonstrate the increased power requirements to maintain altitude in unsteady winds, achieved by increasing the flapping frequency. Flapping fliers may increase their flapping frequency as well as the stroke amplitude to produce the required increase in aerodynamic force, both of these two types of compensatory control requiring additional energetic cost. We analyze the existing data from experiments on animals flying in von Kármán streets and find reasonable agreement with the proposed theoretical model.

preprint2016arXiv

Aerodynamic ground effect in fruitfly sized insect takeoff

Aerodynamic ground effect in flapping-wing insect flight is of importance to comparative morphologies and of interest to the micro-air-vehicle (MAV) community. Recent studies, however, show apparently contradictory results of either some significant extra lift or power savings, or zero ground effect. Here we present a numerical study of fruitfly sized insect takeoff with a specific focus on the significance of leg thrust and wing kinematics. Flapping-wing takeoff is studied using numerical modelling and high performance computing. The aerodynamic forces are calculated using a three-dimensional Navier--Stokes solver based on a pseudo-spectral method with volume penalization. It is coupled with a flight dynamics solver that accounts for the body weight, inertia and the leg thrust, while only having two degrees of freedom: the vertical and the longitudinal horizontal displacement. The natural voluntary takeoff of a fruitfly is considered as reference. The parameters of the model are then varied to explore possible effects of interaction between the flapping-wing model and the ground plane. These modified takeoffs include cases with decreased leg thrust parameter, and/or with periodic wing kinematics, constant body pitch angle. The results show that the ground effect during natural voluntary takeoff is negligible. In the modified takeoffs, when the rate of climb is slow, the difference in the aerodynamic forces due to the interaction with the ground is up to 6%. Surprisingly, depending on the kinematics, the difference is either positive or negative, in contrast to the intuition based on the helicopter theory, which suggests positive excess lift. This effect is attributed to unsteady wing-wake interactions. A similar effect is found during hovering.

preprint2016arXiv

Bumblebees minimize control challenges by combining active and passive modes in unsteady winds

The natural wind environment that volant insects encounter is unsteady and highly complex, posing significant flight control and stability challenges. Unsteady airflows can range from structured chains of discrete vortices shed in the wake of an object to fully developed chaotic turbulence. It is critical to understand the flight control strategies insects employ to safely navigate in natural environments. We combined experiments on free flying bumblebees with high fidelity numerical simulations and lower order modeling to identify the salient mechanics that mediate insect flight in unsteady winds. We trained bumblebees to fly upwind towards an artificial flower in a wind tunnel under steady wind and in a von Karman street (23Hz) formed in the wake of a cylinder. The bees displayed significantly higher movement in the unsteady vortex street compared to steady winds. Correlation analysis revealed that at lower frequencies, less than 10 Hz, in both steady and unsteady winds the bees mediated lateral movement with body roll, typical casting motion. At higher frequencies in unsteady winds there was a negative correlation between body roll and lateral accelerations. Numerical simulations of a bumblebee in similar conditions permitted the separation of the passive and active components of the flight trajectories. Comparison between the free-flying and numerical bees revealed a novel mechanism that enables bees to passively ride out high frequency perturbations while performing active maneuvers and corrections at lower frequencies. The capacity of maintaining stability by combining passive and active modes at different timescales provides a viable means for volant animals and machines to tackle the control challenges posed by complex airflows.

preprint2015arXiv

Adaptive gradient-augmented level set method with multiresolution error estimation

A space-time adaptive scheme is presented for solving advection equations in two space dimensions. The gradient-augmented level set method using a semi-Lagrangian formulation with backward time integration is coupled with a point value multiresolution analysis using Hermite interpolation. Thus locally refined dyadic spatial grids are introduced which are efficiently implemented with dynamic quadtree data structures. For adaptive time integration, an embedded Runge-Kutta method is employed. The precision of the new fully adaptive method is analysed and speed up of CPU time and memory compression with respect to the uniform grid discretization are reported.

preprint2015arXiv

FluSI: A novel parallel simulation tool for flapping insect flight using a Fourier method with volume penalization

FluSI, a fully parallel open source software for pseudo-spectral simulations of three-dimensional flapping flight in viscous flows, is presented. It is freely available for non-commercial use under [https://github.com/pseudospectators/FLUSI]. The computational framework runs on high performance computers with distributed memory architectures. The discretization of the three-dimensional incompressible Navier--Stokes equations is based on a Fourier pseudo-spectral method with adaptive time stepping. The complex time varying geometry of insects with rigid flapping wings is handled with the volume penalization method. The modules characterizing the insect geometry, the flight mechanics and the wing kinematics are described. Validation tests for different benchmarks illustrate the efficiency and precision of the approach. Finally, computations of a model insect in the turbulent regime demonstrate the versatility of the software.

preprint2014arXiv

Analysis and discretization of the volume penalized Laplace operator with Neumann boundary conditions

We study the properties of an approximation of the Laplace operator with Neumann boundary conditions using volume penalization. For the one-dimensional Poisson equation we compute explicitly the exact solution of the penalized equation and quantify the penalization error. Numerical simulations using finite differences allow then to assess the discretisation and penalization errors. The eigenvalue problem of the penalized Laplace operator with Neumann boundary conditions is also studied. As examples in two space dimensions, we consider a Poisson equation with Neumann boundary conditions in rectangular and circular domains.

preprint2014arXiv

Leading-edge vortex shedding from rotating wings

The paper presents a numerical investigation of the leading-edge vortices generated by rotating triangular wings at Reynolds number $Re=250$. A series of three-dimensional numerical simulations have been carried out using a Fourier pseudo-spectral method with volume penalization. The transition from stable attachment of the leading-edge vortex to periodic vortex shedding is explored, as a function of the wing aspect ratio and the angle of attack. It is found that, in a stable configuration, the spanwise flow in the recirculation bubble past the wing is due to the centrifugal force, incompressibility and viscous stresses. For the flow outside of the bubble, an inviscid model of spanwise flow is presented.

preprint2012arXiv

Approximation of the Laplace and Stokes operators with Dirichlet boundary conditions through volume penalization: a spectral viewpoint

We report the results of a detailed study of the spectral properties of Laplace and Stokes operators, modified with a volume penalization term designed to approximate Dirichlet conditions in the limit when a penalization parameter, $η$, tends to zero. The eigenvalues and eigenfunctions are determined either analytically or numerically as functions of $η$, both in the continuous case and after applying Fourier or finite difference discretization schemes. For fixed $η$, we find that only the part of the spectrum corresponding to eigenvalues $λ\lesssim η^{-1}$ approaches Dirichlet boundary conditions, while the remainder of the spectrum is made of uncontrolled, spurious wall modes. The penalization error for the controlled eigenfunctions is estimated as a function of $η$ and $λ$. Surprisingly, in the Stokes case, we show that the eigenfunctions approximately satisfy, with a precision $O(η)$, Navier slip boundary conditions with slip length equal to $\sqrtη$. Moreover, for a given discretization, we show that there exists a value of $η$, corresponding to a balance between penalization and discretization errors, below which no further gain in precision is achieved. These results shed light on the behavior of volume penalization schemes when solving the Navier-Stokes equations, outline the limitations of the method, and give indications on how to choose the penalization parameter in practical cases.