Researcher profile

Kai Schneider

Kai Schneider contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

18 published item(s)

preprint2026arXiv

On the local anisotropy of quasi-two-dimensional forced shallow flow: An experimental study

We experimentally investigate quasi-two-dimensional (Q2D) forced shallow flows in the presence of solid boundaries and analyze the deviation from the Kolmogorov-Kraichnan (KK) theory. Complex motion is generated using a thin electrolyte subject to electromagnetic forces, and we employ particle tracking velocimetry to resolve the flow properties down to the Kolmogorov scale. Although the velocity probability distribution function closely resembles a Gaussian, deviations from Gaussianity emerge for velocity increments as scales decrease. The second-order structure function supports the onset of local anisotropy at small scales. The sign of the third-order structure function indicates the dominance of the inverse cascade in energy transfer, and the cross-correlation between longitudinal and transverse directions proves to be significant at large scales. The breakdown of local isotropy is consistent with the effect of bottom friction, which primarily affects the longitudinal motion, while leaving the perpendicular direction unaffected. This local anisotropy propagates to larger scales via the inverse energy cascade, with nonlinear interactions eventually influencing the perpendicular direction.

preprint2025arXiv

Nondiffusive transport of inertial heavy impurities in drift-wave turbulence

We investigate the transport behavior of tungsten impurities with finite inertia in drift-wave turbulence using the Hasegawa-Wakatani model. Unlike previous tracer-based models, our simulations reveal a transition to non-diffusive dynamics for a range of charge states. This novel mechanism offers a turbulence-driven route to core impurity accumulation. This finding underscores the nontrivial role of particle inertia in impurity dynamics and has strong implications for impurity control in future fusion devices.

preprint2024arXiv

Singularity formation of vortex sheets in 2D Euler equations using the characteristic mapping method

The goal of this numerical study is to get insight into singular solutions of the two-dimensional (2D) Euler equations for non-smooth initial data, in particular for vortex sheets. To this end high resolution computations of vortex layers in 2D incompressible Euler flows are performed using the characteristic mapping method (CMM). This semi-Lagrangian method evolves the flow map using the gradient-augmented level set method (GALS). The semi-group structure of the flow map allows its decomposition into sub-maps (each over a finite time interval), and thus the precision can be controlled by choosing appropriate remapping times. Composing the flow map yields exponential resolution in linear time, a unique feature of CMM, and thus fine scale flow structures can be resolved in great detail. Here the roll-up process of vortex layers is studied varying the thickness of the layer showing its impact on the growth of palinstrophy and possible blow up of absolute vorticity. The curvature of the vortex sheet shows a singular-like behavior. The self-similar structure of the vortex core is investigated in the vanishing thickness limit. Conclusions on the non-uniqueness of weak solutions of 2D Euler for non-smooth initial data are drawn and the presence of flow singularities is revealed tracking them in the complex plane.

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.

preprint2021arXiv

A wavelet-adaptive method for multiscale simulation of turbulent flows in flying insects

We present a wavelet-based adaptive method for computing 3D multiscale flows in complex, time-dependent geometries, implemented on massively parallel computers. While our focus is on simulations of flapping insects, it can be used for other flow problems, including turbulence, as well. The incompressible fluid is modeled with an artificial compressibility approach in order to avoid solving elliptical problems. No-slip and in/outflow boundary conditions are imposed using volume penalization. The governing equations are discretized on a locally uniform Cartesian grid with centered finite differences, and integrated in time with a Runge--Kutta scheme, both of 4th order. The domain is partitioned into cubic blocks with equidistant grids with different resolution and, for each block, biorthogonal interpolating wavelets are used as refinement indicators and prediction operators. Thresholding the wavelet coefficients allows to generate dynamically evolving grids, and an adaption strategy tracks the solution in both space and scale. Blocks are distributed among MPI processes and the global topology of the grid is encoded using a tree-like data structure. Analyzing the different physical and numerical parameters allows balancing their individual error contributions and thus ensures optimal convergence while minimizing computational effort. Different validation tests score accuracy and performance of our new open source code, WABBIT (Wavelet Adaptive Block-Based solver for Interactions with Turbulence), on massively parallel computers using fully adaptive grids. Flow simulations of flapping insects demonstrate its applicability to complex, bio-inspired problems.

preprint2021arXiv

Adaptive two- and three-dimensional multiresolution computations of resistive magnetohydrodynamics

Fully adaptive computations of the resistive magnetohydrodynamic (MHD) equations are presented in two and three space dimensions using a finite volume discretization on locally refined dyadic grids. Divergence cleaning is used to control the incompressibility constraint of the magnetic field. For automatic grid adaptation a cell-averaged multiresolution analysis is applied which guarantees the precision of the adaptive computations, while reducing CPU time and memory requirements. Implementation issues of the open source code CARMEN-MHD are discussed. To illustrate its precision and efficiency different benchmark computations including shock-cloud interaction and magnetic reconnection are presented.

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.

preprint2018arXiv

Wavelet Methods for Studying the Onset of Strong Plasma Turbulence

Wavelet basis functions are a natural tool for analyzing turbulent flows containing localized coherent structures of different spatial scales. Here, wavelets are used to study the onset and subsequent transition to fully developed turbulence from a laminar state. Originally applied to neutral fluid turbulence, an iterative wavelet technique decomposes the field into coherent and incoherent contributions. In contrast to Fourier power spectra, finite time Lyapunov exponents (FTLE), and simple measures of intermittency such as non-Gaussian statistics of field increments, the wavelet technique is found to provide a quantitative measure for the onset of turbulence and to track the transition to fully developed turbulence. The wavelet method makes no assumptions about the structure of the coherent current sheets or the underlying plasma model. Temporal evolution of the coherent and incoherent wavelet fluctuations is found to be highly correlated with the magnetic field energy and plasma thermal energy, respectively. The onset of turbulence is identified with the rapid growth of a background of incoherent fluctuations spreading across a range of scales and a corresponding drop in the coherent components. This is suggestive of the interpretation of the coherent and incoherent wavelet fluctuations as measures of coherent structures (e.g., current sheets) and dissipation, respectively. The ratio of the incoherent to coherent fluctuations $R_{ic}$ is found to be fairly uniform across different plasma models and provides an empirical threshold for turbulence onset. The technique is illustrated through examples. First, it is applied to the Kelvin--Helmholtz instability from different simulation models including fully kinetic, hybrid (kinetic ion/fluid electron), and Hall MHD simulations. Second, it is applied to the development of turbulence downstream of the bowshock in a magnetosphere simulation.

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.

preprint2012arXiv

Conditional vorticity budget of coherent and incoherent flow contributions in fully developed homogeneous isotropic turbulence

We investigate the conditional vorticity budget of fully developed three-dimensional homogeneous isotropic turbulence with respect to coherent and incoherent flow contributions. The Coherent Vorticity Extraction based on orthogonal wavelets allows to decompose the vorticity field into coherent and incoherent contributions, of which the latter are noise-like. The impact of the vortex structures observed in fully developed turbulence on statistical balance equations is quantified considering the conditional vorticity budget. The connection between the basic structures present in the flow and their statistical implications is thereby assessed. The results are compared to those obtained for large- and small-scale contributions using a Fourier decomposition, which reveals pronounced differences.

preprint2012arXiv

Fully adaptive multiresolution schemes for strongly degenerate parabolic equations in one space dimension

We present a fully adaptive multiresolution scheme for spatially one-dimensional quasilinear strongly degenerate parabolic equations with zero-flux and periodic boundary conditions. The numerical scheme is based on a finite volume discretization using the Engquist-Osher numerical flux and explicit time stepping. An adaptive multiresolution scheme based on cell averages is then used to speed up the CPU time and the memory requirements of the underlying finite volume scheme, whose first-order version is known to converge to an entropy solution of the problem. A particular feature of the method is the storage of the multiresolution representation of the solution in a graded tree, whose leaves are the non-uniform finite volumes on which the numerical divergence is eventually evaluated. Moreover using the $L^1$ contraction of the discrete time evolution operator we derive the optimal choice of the threshold in the adaptive multiresolution method. Numerical examples illustrate the computational efficiency together with the convergence properties.

preprint2011arXiv

The role of coherent vorticity in turbulent transport in resistive drift-wave turbulence

The coherent vortex extraction method, a wavelet technique for extracting coherent vortices out of turbulent flows, is applied to simulations of resistive drift-wave turbulence in magnetized plasma (Hasegawa-Wakatani system). The aim is to retain only the essential degrees of freedom, responsible for the transport. It is shown that the radial density flux is carried by these coherent modes. In the quasi-hydrodynamic regime, coherent vortices exhibit depletion of the polarization-drift nonlinearity and vorticity strongly dominates strain, in contrast to the quasiadiabatic regime.

preprint2010arXiv

Inertial range scaling of the scalar flux spectrum in two-dimensional turbulence

Two-dimensional statistically stationary isotropic turbulence with an imposed uniform scalar gradient is investigated. Dimensional arguments are presented to predict the inertial range scaling of the turbulent scalar flux spectrum in both the inverse cascade range and the enstrophy cascade range for small and unity Schmidt numbers. The scaling predictions are checked by direct numerical simulations and good agreement is observed.

preprint2010arXiv

Origin of Lagrangian Intermittency in Drift-Wave Turbulence

The Lagrangian velocity statistics of dissipative drift-wave turbulence are investigated. For large values of the adiabaticity (or small collisionality), the probability density function of the Lagrangian acceleration shows exponential tails, as opposed to the stretched exponential or algebraic tails, generally observed for the highly intermittent acceleration of Navier-Stokes turbulence. This exponential distribution is shown to be a robust feature independent of the Reynolds number. For small adiabaticity, algebraic tails are observed, suggesting the strong influence of point-vortex-like dynamics on the acceleration. A causal connection is found between the shape of the probability density function and the autocorrelation of the norm of the acceleration.

preprint2010arXiv

Rapid generation of angular momentum in bounded magnetized plasma

Direct numerical simulations of two-dimensional decaying MHD turbulence in bounded domains show the rapid generation of angular momentum in nonaxisymmetric geometries. It is found that magnetic fluctuations enhance this mechanism. On a larger time scale, the generation of a magnetic angular momentum, or angular field, is observed. For axisymmetric geometries, the generation of angular momentum is absent; nevertheless, a weak magnetic field can be observed. The derived evolution equations for both the angular momentum and angular field yield possible explanations for the observed behavior.

preprint2010arXiv

Self-organization and symmetry-breaking in two-dimensional plasma turbulence

The spontaneous self-organization of two-dimensional magnetized plasma is investigated within the framework of magnetohydrodynamics with a particular emphasis on the symmetry-breaking induced by the shape of the confining boundaries. This symmetry-breaking is quantified by the angular momentum, which is shown to be generated rapidly and spontaneously from initial conditions free from angular momentum as soon as the geometry lacks axisymmetry. This effect is illustrated by considering circular, square, and elliptical boundaries. It is shown that the generation of angular momentum in nonaxisymmetric geometries can be enhanced by increasing the magnetic pressure. The effect becomes stronger at higher Reynolds numbers. The generation of magnetic angular momentum (or angular field), previously observed at low Reynolds numbers, becomes weaker at larger Reynolds numbers.

preprint2008arXiv

Extreme Lagrangian acceleration in confined turbulent flow

A Lagrangian study of two-dimensional turbulence for two different geometries, a periodic and a confined circular geometry, is presented to investigate the influence of solid boundaries on the Lagrangian dynamics. It is found that the Lagrangian acceleration is even more intermittent in the confined domain than in the periodic domain. The flatness of the Lagrangian acceleration as a function of the radius shows that the influence of the wall on the Lagrangian dynamics becomes negligible in the center of the domain and it also reveals that the wall is responsible for the increased intermittency. The transition in the Lagrangian statistics between this region, not directly influenced by the walls, and a critical radius which defines a Lagrangian boundary layer, is shown to be very sharp with a sudden increase of the acceleration flatness from about 5 to about 20.