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

Rainer Grauer

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

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-phcond-mat.stat-mechphysics.plasm-phastro-ph.SRmath-phmath.APmath.MPmath.NAmath.PRNumerical Analysisphysics.data-anphysics.space-ph

Catalog footprint

What is connected

18works
13topics
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

18 published item(s)

preprint2022arXiv

Spontaneous Symmetry Breaking for Extreme Vorticity and Strain in the 3D Navier-Stokes Equations

We investigate the spatio-temporal structure of the most likely configurations realising extremely high vorticity or strain in the stochastically forced 3D incompressible Navier-Stokes equations. Most likely configurations are computed by numerically finding the highest probability velocity field realising an extreme constraint as solution of a large optimisation problem. High-vorticity configurations are identified as pinched vortex filaments with swirl, while high-strain configurations correspond to counter-rotating vortex rings. We additionally observe that the most likely configurations for vorticity and strain spontaneously break their rotational symmetry for extremely high observable values. Instanton calculus and large deviation theory allow us to show that these maximum likelihood realisations determine the tail probabilities of the observed quantities. In particular, we are able to demonstrate that artificially enforcing rotational symmetry for large strain configurations leads to a severe underestimate of their probability, as it is dominated in likelihood by an exponentially more likely symmetry broken vortex-sheet configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stochastic interpolation of sparsely sampled time series by a superstatistical random process and its synthesis in Fourier and wavelet space

We present a novel method for stochastic interpolation of sparsely sampled time signals based on a superstatistical random process generated from a multivariate Gaussian scale mixture. In comparison to other stochastic interpolation methods such as Gaussian process regression, our method possesses strong multifractal properties and is thus applicable to a broad range of real-world time series, e.g. from solar wind or atmospheric turbulence. Furthermore, we provide a sampling algorithm in terms of a mixing procedure that consists of generating a 1 + 1-dimensional field u(t, ξ), where each Gaussian component uξ(t) is synthesized with identical underlying noise but different covariance function Cξ(t,s) parameterized by a log-normally distributed parameter ξ. Due to the Gaussianity of each component uξ(t), we can exploit standard sampling alogrithms such as Fourier or wavelet methods and, most importantly, methods to constrain the process on the sparse measurement points. The scale mixture u(t) is then initialized by assigning each point in time t a ξ(t) and therefore a specific value from u(t, ξ), where the time-dependent parameter ξ(t) follows a log-normal process with a large correlation time scale compared to the correlation time of u(t, ξ). We juxtapose Fourier and wavelet methods and show that a multiwavelet-based hierarchical approximation of the interpolating paths, which produce a sparse covariance structure, provide an adequate method to locally interpolate large and sparse datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Symmetries and zero modes in sample path large deviations

Sharp large deviation estimates for stochastic differential equations with small noise, based on minimizing the Freidlin-Wentzell action functional under appropriate boundary conditions, can be obtained by integrating certain matrix Riccati differential equations along the large deviation minimizers or instantons, either forward or backward in time. Previous works in this direction often rely on the existence of isolated minimizers with positive definite second variation. By adopting techniques from field theory and explicitly evaluating the large deviation prefactors as functional determinant ratios using Forman's theorem, we extend the approach to general systems where degenerate submanifolds of minimizers exist. The key technique for this is a boundary-type regularization of the second variation operator. This extension is particularly relevant if the system possesses continuous symmetries that are broken by the instantons. We find that removing the vanishing eigenvalues associated with the zero modes is possible within the Riccati formulation and amounts to modifying the initial or final conditions and evaluation of the Riccati matrices. We apply our results in multiple examples including a dynamical phase transition for the average surface height in short-time large deviations of the one-dimensional Kardar-Parisi-Zhang equation with flat initial profile.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events

Vlasov solvers that operate on a phase-space grid are highly accurate but also numerically demanding. Coarse velocity space resolutions, which are unproblematic in particle-in-cell (PIC) simulations, can lead to numerical heating or oscillations in standard continuum Vlasov methods. We present a new dual Vlasov solver which is based on an established positivity preserving advection scheme for the update of the distribution function and an energy conserving partial differential equation solver for the kinetic update of mean velocity and temperature. The solvers work together via moment fitting during which the maximum entropy part of the distribution function is replaced by the solution from the partial differential equation solver. This numerical scheme makes continuum Vlasov methods competitive with PIC methods concerning computational cost and enables us to model large scale reconnection in Earth's magnetosphere with a fully kinetic continuum method. The simulation results agree well with measurements by the MMS spacecraft.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficient Computation of Instantons for Multi-Dimensional Turbulent Flows with Large Scale Forcing

Extreme events play a crucial role in fluid turbulence. Inspired by methods from field theory, these extreme events, their evolution and probability can be computed with help of the instanton formalism as minimizers of a suitable action functional. Due to the high number of degrees of freedom in multi-dimensional fluid flows, traditional global minimization techniques quickly become prohibitive in their memory requirements. We outline a novel method for finding the minimizing trajectory in a wide class of problems that typically occurs in turbulence setups, where the underlying dynamical system is a non-gradient, non-linear partial differential equation, and the forcing is restricted to a limited length scale. We demonstrate the efficiency of the algorithm in terms of performance and memory by computing high resolution instanton field configurations corresponding to viscous shocks for 1D and 2D compressible flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Irreversibility and small-scale generation in 3D turbulent flows

In three-dimensional turbulent flows energy is supplied at large scales and cascades down to the smallest scales where viscosity dominates. The flux of energy through scales implies the generation of small scales from larger ones, which is the fundamental reason for the irreversibility of the dynamics of turbulent flows. As we showed recently, this irreversibility manifests itself by an asymmetry of the probability distribution of the instantaneous power $p$ of the forces acting on fluid elements. In particular, the third moment of $p$ was found to be negative. Yet, a physical connection between the irreversibility manifested in the distribution of $p$ and the energy flux or small-scale generation in turbulence has not been established. Here, with analytical calculations and support from numerical simulations of fully developed turbulence, we connect the asymmetry in the power distribution, {\it i.e.}, the negative value of $\langle p^3 \rangle$, to the generation of small scales, or more precisely, to the amplification (stretching) of vorticity in turbulent flows. Our result is the first step towards a quantitative understanding of the origin of the irreversibility observed at the level of individual Lagrangian trajectories in turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The instanton method and its numerical implementation in fluid mechanics

A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional fluid dynamical problems. We illustrate these ideas using the two-dimensional Burgers equation and the three-dimensional Navier-Stokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relevance of instantons in Burgers turbulence

Instanton calculations are performed in the context of stationary Burgers turbulence to estimate the tails of the probability density function (PDF) of velocity gradients. These results are then compared to those obtained from massive direct numerical simulations (DNS) of the randomly forced Burgers equation. The instanton predictions are shown to agree with the DNS in a wide range of regimes, including those that are far from the limiting cases previously considered in the literature. These results settle the controversy of the relevance of the instanton approach for the prediction of the velocity gradient PDF tail exponents. They also demonstrate the usefulness of the instanton formalism in Burgers turbulence, and suggest that this approach may be applicable in other contexts, such as 2D and 3D turbulence in compressible and incompressible flows.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structures and Lagrangian statistics of the Taylor-Green Dynamo

The evolution of a Taylor-Green forced magnetohydrodynamic (MHD) system showing dynamo activity is analyzed via direct numerical simulations. The statistical properties of the velocity and magnetic field in Eulerian coordinates and along trajectories of fluid elements (Lagrangian coordinates) are found to change between the kinematic, non-linear and saturated regime. Fluid element (tracer) trajectories change from chaotic quasi-isotropic (kinematic phase) to mean magnetic field aligned (saturated phase). The probability density functions (PDFs) of the magnetic field change from strongly non-Gaussian in the kinematic to quasi-Gaussian PDFs in the saturated regime so that their flatness give a precise handle on the definition of the limiting points of the three regimes. Also the statistics of the fluctuations of the kinetic and magnetic energy along fluid trajectories change. All this goes along with a dramatic increase of the correlation time of velocity and magnetic field fluctuations experienced by tracers significantly exceeding one turbulent large-eddy turn-over time. A remarkable consequence is an intermittent scaling regime of the Lagrangian magnetic field structure functions at unusually long time scales.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effect of turbulent fluctuations on the drag and lift forces on a towed sphere and its boundary layer

The impact of turbulent fluctuations on the forces exerted by a fluid on a towed spherical particle is investigated by means of high-resolution direct numerical simulations. The measurements are carried out using a novel scheme to integrate the two-way coupling between the particle and the incompressible surrounding fluid flow maintained in a high-Reynolds-number turbulent regime. The main idea consists in combining a Fourier pseudo-spectral method for the fluid with an immersed-boundary technique to impose the no-slip boundary condition on the surface of the particle. Benchmarking of the code shows a good agreement with experimental and numerical measurements from other groups. A study of the turbulent wake downstream the sphere is also reported. The mean velocity deficit is shown to behave as the inverse of the distance from the particle, as predicted from classical similarity analysis. This law is reinterpreted in terms of the principle of "permanence of large eddies" that relates infrared asymptotic self-similarity to the law of decay of energy in homogeneous turbulence. The developed method is then used to attack the problem of an upstream flow that is in a developed turbulent regime. It is shown that the average drag force increases as a function of the turbulent intensity and the particle Reynolds number. This increase is significantly larger than predicted by standard drag correlations based on laminar upstream flows. It is found that the relevant parameter is the ratio of the viscous boundary layer thickness to the dissipation scale of the ambient turbulent flow. The drag enhancement can be motivated by the modification of the mean velocity and pressure profile around the sphere by small scale turbulent fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A novel code for numerical 3-D MHD studies of CME expansion

A recent third-order, essentially non-oscillatory central scheme to advance the equations of single-fluid magnetohydrodynamics (MHD) in time has been implemented into a new numerical code. This code operates on a 3-D Cartesian, non-staggered grid, and is able to handle shock-like gradients without producing spurious oscillations. To demonstrate the suitability of our code for the simulation of coronal mass ejections (CMEs) and similar heliospheric transients, we present selected results from test cases and perform studies of the solar wind expansion during phases of minimum solar activity. We can demonstrate convergence of the system into a stable Parker-like steady state for both hydrodynamic and MHD winds. The model is subsequently applied to expansion studies of CME-like plasma bubbles, and their evolution is monitored until a stationary state similar to the initial one is achieved. In spite of the model's (current) simplicity, we can confirm the CME's nearly self-similar evolution close to the Sun, thus highlighting the importance of detailed modelling especially at small heliospheric radii. Additionally, alternative methods to implement boundary conditions at the coronal base, as well as strategies to ensure a solenoidal magnetic field, are discussed and evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-Time Euler singularities: A Lagrangian perspective

We address the question whether a singularity in a three-dimensional incompressible inviscid fluid flow can occur in finite time. Analytical considerations and numerical simulations suggest high-symmetry flows being a promising candidate for a finite-time blowup. Utilizing Lagrangian and geometric non-blowup criteria, we present numerical evidence against the formation of a finite-time singularity for the high-symmetry vortex dodecapole initial condition. We use data obtained from high resolution adaptively refined numerical simulations and inject Lagrangian tracer particles to monitor geometric properties of vortex line segments. We then verify the assumptions made by analytical non-blowup criteria introduced by Deng et. al [Commun. PDE 31 (2006)] connecting vortex line geometry (curvature, spreading) to velocity increase to rule out singular behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Instanton filtering for the stochastic Burgers equation

We address the question whether one can identify instantons in direct numerical simulations of the stochastically driven Burgers equation. For this purpose, we first solve the instanton equations using the Chernykh-Stepanov method [Phys. Rev. E 64, 026306 (2001)]. These results are then compared to direct numerical simulations by introducing a filtering technique to extract prescribed rare events from massive data sets of realizations. Using this approach we can extract the entire time history of the instanton evolution which allows us to identify the different phases predicted by the direct method of Chernykh and Stepanov with remarkable agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lagrangian and geometric analysis of finite-time Euler singularities

We present a numerical method of analyzing possibly singular incompressible 3D Euler flows using massively parallel high-resolution adaptively refined numerical simulations up to 8192^3 mesh points. Geometrical properties of Lagrangian vortex line segments are used in combination with analytical non-blowup criteria by Deng et al [Commun. PDE 31 (2006)] to reliably distinguish between singular and near-singular flow evolution. We then apply the presented technique to a class of high-symmetry initial conditions and present numerical evidence against the formation of a finite-time singularity in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Longitudinal and Transverse structure functions in high Reynolds-number turbulence

Using exact relations between velocity structure functions (Hill, Hill and Boratav, and Yakhot) and neglecting pressure contributions in a first approximation, we obtain a closed system and derive simple order-dependent rescaling relationships between longitudinal and transverse structure functions. By means of numerical data with turbulent Reynolds numbers ranging from $\Re_λ=320$ to $\Re_λ=730$, we establish a clear correspondence between their respective scaling range, while confirming that their scaling exponents do differ. This difference does not seem to depend on Reynolds number. Making use of the Mellin transform, we further map longitudinal to (rescaled) transverse probability density functions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Turbulence properties and global regularity of a modified Navier-Stokes equation

We introduce a modification of the Navier-Stokes equation that has the remarkable property of possessing an infinite number of conserved quantities in the inviscid limit. This new equation is studied numerically and turbulence properties are analyzed concerning energy spectra and scaling of structure functions. The dissipative structures arising in this new equation are curled vortex sheets contrary to vortex tubes arising in Navier-Stokes turbulence. The numerically calculated scaling of structure functions is compared with a phenomenological model based on the She-Lévêque approach. Finally, for this equation we demonstrate global well-posedness for sufficiently smooth initial conditions in the periodic case and in $\mathbb R^3$. The key feature is the availability of an additional estimate which shows that the $L^4$-norm of the velocity field remains finite.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Simulation of Current Sheet Formation in a Quasi-Separatrix Layer using Adaptive Mesh Refinement

The formation of a thin current sheet in a magnetic quasi-separatrix layer (QSL) is investigated by means of numerical simulation using a simplified ideal, low-$β$, MHD model. The initial configuration and driving boundary conditions are relevant to phenomena observed in the solar corona and were studied earlier by Aulanier et al., A&A 444, 961 (2005). In extension to that work, we use the technique of adaptive mesh refinement (AMR) to significantly enhance the local spatial resolution of the current sheet during its formation, which enables us to follow the evolution into a later stage. Our simulations are in good agreement with the results of Aulanier et al. up to the calculated time in that work. In a later phase, we observe a basically unarrested collapse of the sheet to length scales that are more than one order of magnitude smaller than those reported earlier. The current density attains correspondingly larger maximum values within the sheet. During this thinning process, which is finally limited by lack of resolution even in the AMR studies, the current sheet moves upward, following a global expansion of the magnetic structure during the quasi-static evolution. The sheet is locally one-dimensional and the plasma flow in its vicinity, when transformed into a co-moving frame, qualitatively resembles a stagnation point flow. In conclusion, our simulations support the idea that extremely high current densities are generated in the vicinities of QSLs as a response to external perturbations, with no sign of saturation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Conditional Eulerian and Lagrangian velocity increment statistics of fully developed turbulent flow

Conditional statistics of homogeneous isotropic turbulent flow is investigated by means of high-Reynolds number direct numerical simulations performed with $2048^3$ collocation points. Eulerian as well as Lagrangian velocity increment statistics under several conditions are analyzed and compared. In agreement with experimental data longitudinal probability density functions $P(δ^{||}_l u|ε_l)$ conditioned on a scale-averaged energy dissipation rate are close to Gaussian distributions over all scales within the inertial range of scales. Also transverse increments conditioned on either the dissipation rate or the square of the vorticity have quasi-Gaussian probability distribution functions (PDFs). Concerning Lagrangian statistics we found that conditioning on a trajectory averaged energy-dissipation rate $ε_τ$ significantly reduces the scale dependence of the increment PDFs $P(δ_τu_i|ε_τ)$. By means of dimensional arguments we propose a novel condition for Lagrangian increments which is shown to reduce even more the flatness of the corresponding PDFs and thus intermittency in the inertial range of scales. The conditioned Lagrangian PDF corresponding to the smallest increment considered is reasonably well described by the K41-prediction of the PDF of acceleration. Conditioned structure functions show approximately K41-scaling with a larger scaling range than the unconditioned ones.

0 citations0 reviewsNo code linkedOpen access