BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kai Schneider

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

ResearcherUnclaimed source record
physics.flu-dynmath.NAphysics.plasm-phphysics.comp-phNumerical AnalysisBiological Physicsmath.APnlin.AOnlin.CDphysics.class-phphysics.data-anQuantitative Methods

Catalog footprint

What is connected

31works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparison of adaptive multiresolution and adaptive mesh refinement applied to simulations of the compressible Euler equations

We present a detailed comparison between two adaptive numerical approaches to solve partial differential equations (PDEs), adaptive multiresolution (MR) and adaptive mesh refinement (AMR). Both discretizations are based on finite volumes in space with second order shock-capturing, and explicit time integration either with or without local time-stepping. The two methods are benchmarked for the compressible Euler equations in Cartesian geometry. As test cases a 2D Riemann problem, Lax-Liu 6, and a 3D ellipsoidally expanding shock wave have been chosen. We compare and assess their computational efficiency in terms of CPU time and memory requirements. We evaluate the accuracy by comparing the results of the adaptive computations with those obtained with the corresponding FV scheme using a regular fine mesh. We find that both approaches yield similar trends for CPU time compression for increasing number of refinement levels. MR exhibits more efficient memory compression than AMR and shows slightly enhanced convergence; however, a larger absolute overhead is measured for the tested codes.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An adaptive multiresolution method for ideal magnetohydrodynamics using divergence cleaning with parabolic-hyperbolic correction

We present an adaptive multiresolution method for the numerical simulation of ideal magnetohydrodynamics in two space dimensions. The discretization uses a finite volume scheme based on a Cartesian mesh and an explicit compact Rung-Kutta scheme for time integration. Harten's cell average multiresolution allows to introduce a locally refined spatial mesh while controlling the error. The incompressibility of the magnetic field is controlled by using a Generalized Lagrangian Multiplier (GLM) approach with a mixed hyperbolic-parabolic correction. Different applications to two-dimensional problems illustrate the properties of the method. For each application CPU time and memory savings are reported and numerical aspects of the method are discussed. The accuracy of the adaptive computations is assessed by comparison with reference solutions computed on a regular fine mesh.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Immersed boundary methods for numerical simulation of confined fluid and plasma turbulence in complex geometries: a review

Immersed boundary methods for computing confined fluid and plasma flows in complex geometries are reviewed. The mathematical principle of the volume penalization technique is described and simple examples for imposing Dirichlet and Neumann boundary conditions in one dimension are given. Applications for fluid and plasma turbulence in two and three space dimensions illustrate the applicability and the efficiency of the method in computing flows in complex geometries, for example in toroidal geometries with asymmetric poloidal cross-sections.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Wavelet transforms and their applications to MHD and plasma turbulence: a review

Wavelet analysis and compression tools are reviewed and different applications to study MHD and plasma turbulence are presented. We introduce the continuous and the orthogonal wavelet transform and detail several statistical diagnostics based on the wavelet coefficients. We then show how to extract coherent structures out of fully developed turbulent flows using wavelet-based denoising. Finally some multiscale numerical simulation schemes using wavelets are described. Several examples for analyzing, compressing and computing one, two and three dimensional turbulent MHD or plasma flows are presented.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Wavelet methods to eliminate resonances in the Galerkin-truncated Burgers and Euler equations

It is well known that solutions to the Fourier-Galerkin truncation of the inviscid Burgers equation (and other hyperbolic conservation laws) do not converge to the physically relevant entropy solution after the formation of the first shock. This loss of convergence was recently studied in detail in [S. S. Ray et al., Phys. Rev. E 84, 016301 (2011)], and traced back to the appearance of a spatially localized resonance phenomenon perturbing the solution. In this work, we propose a way to remove this resonance by filtering a wavelet representation of the Galerkin-truncated equations. A method previously developed with a complex-valued wavelet frame is applied and expanded to embrace the use of real-valued orthogonal wavelet basis, which we show to yield satisfactory results only under the condition of adding a safety zone in wavelet space. We also apply the complex-valued wavelet based method to the 2D Euler equation problem, showing that it is able to filter the resonances in this case as well.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Intrinsic rotation of toroidally confined magnetohydrodynamics

The spatiotemporal self-organization of viscoresistive magnetohydrodynamics (MHD) in a toroidal geometry is studied. Curl-free toroidal magnetic and electric fields are imposed. It is observed in our simulations that a flow is generated, which evolves from dominantly poloidal to toroidal when the Lundquist numbers are increased. It is shown that this toroidal organization of the flow is consistent with the tendency of the velocity field to align with the magnetic field. Up-down asymmetry of the geometry causes the generation of a non-zero toroidal angular momentum.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access