Preconditioned Minimal Residual Methods for Chebyshev Spectral Caluclations

The problem of preconditioning the pseudospectral Chebyshev approximation of an elliptic operator is considered. The numerical sensitiveness to variations of the coefficients of the operator are investigated for two classes of preconditioning matrices: one arising from finite differences, the other from finite elements. The preconditioned system is solved by a conjugate gradient type method, and by a DuFort-Frankel method with dynamical parameters. The methods are compared on some test problems with the Richardson method and with the minimal residual Richardson method

On the boundary treatment in spectral methods for hyperbolic systems

Spectral methods were successfully applied to the simulation of slow transients in gas transportation networks. Implicit time advancing techniques are naturally suggested by the nature of the problem. The correct treatment of the boundary conditions are clarified in order to avoid any stability restriction originated by the boundaries. The Beam and Warming and the Lerat schemes are unconditionally linearly stable when used with a Chebyshev pseudospectral method. Engineering accuracy for a gas transportation problem is achieved at Courant numbers up to 100

Non-axisymmetric instability of shear-banded Taylor-Couette flow

Recent experiments show that shear-banded flows of semi-dilute worm-like micelles in Taylor-Couette geometry exhibit a flow instability in the form of Taylor-like vortices. Here we perform the non-axisymmetric linear stability analysis of the diffusive Johnson-Segalman model of shear banding and show that the nature of this instability depends on the applied shear rate. For the experimentally relevant parameters, we find that at the beginning of the stress plateau the instability is driven by the interface between the bands, while most of the stress plateau is occupied by the bulk instability of the high-shear-rate band. Our work significantly alters the recently proposed stability diagram of shear-banded flows based on axisymmetric analysis.Comment: 6 pages, 5 figures, main text and supplementary material; accepted to Phys. Rev. Let

An Eulerian Approach to the Analysis of Krause's Consensus Models

Abstract. In this paper we analyze a class of multi-agent consensus dynamical systems inspired by Krause’s original model. As in Krause’s, the basic assumption is the so-called bounded confidence: two agents can influence each other only when their state values are below a given distance threshold R. We study the system under an Eulerian point of view considering (possibly continuous) probability distributions of agents and we present original convergence results. The limit distribution is always necessarily a convex combination of delta functions at least R far apart from each other: in other terms these models are locally aggregating. The Eulerian perspective provides the natural framework for designing a numerical algorithm, by which we obtain several simulations in 1 and 2 dimensions

Modification of three-dimensional transition in the wake of a rotationally oscillating cylinder

A study of the flow past an oscillatory rotating cylinder has been conducted, where the frequency of oscillation has been matched to the natural frequency of the vortex street generated in the wake of a stationary cylinder, at Reynolds number 300. The focus is on the wake transition to three-dimensional flow and, in particular, the changes induced in this transition by the addition of the oscillatory rotation. Using Floquet stability analysis, it is found that the fine-scale three-dimensional mode that typically dominates the wake at a Reynolds number beyond that at the second transition to three-dimensional flow (referred to as mode B) is suppressed for amplitudes of rotation beyond a critical amplitude, in agreement with past studies. However, the rotation does not suppress the development of three-dimensionality completely, as other modes are discovered that would lead to three-dimensional flow. In particular, the longer-wavelength mode that leads the three-dimensional transition in the wake of a stationary cylinder (referred to as mode A) is left essentially unaffected at low amplitudes of rotation. At higher amplitudes of oscillation, mode A is also suppressed as the two-dimensional near wake changes in character from a single- to a double- row wake; however, another mode is predicted to render the flow three-dimensional, dubbed mode D (for double row). This mode has the same spatio-temporal symmetries as mode A

Spectral methods for exterior elliptic problems

Spectral approximations for exterior elliptic problems in two dimensions are discussed. As in the conventional finite difference or finite element methods, the accuracy of the numerical solutions is limited by the order of the numerical farfield conditions. A spectral boundary treatment is introduced at infinity which is compatible with the infinite order interior spectral scheme. Computational results are presented to demonstrate the spectral accuracy attainable. Although a simple Laplace problem is examined, the analysis covers more complex and general cases

A streamwise-constant model of turbulent pipe flow

A streamwise-constant model is presented to investigate the basic mechanisms responsible for the change in mean flow occuring during pipe flow transition. Using a single forced momentum balance equation, we show that the shape of the velocity profile is robust to changes in the forcing profile and that both linear non-normal and nonlinear effects are required to capture the change in mean flow associated with transition to turbulence. The particularly simple form of the model allows for the study of the momentum transfer directly by inspection of the equations. The distribution of the high- and low-speed streaks over the cross-section of the pipe produced by our model is remarkably similar to one observed in the velocity field near the trailing edge of the puff structures present in pipe flow transition. Under stochastic forcing, the model exhibits a quasi-periodic self-sustaining cycle characterized by the creation and subsequent decay of "streamwise-constant puffs", so-called due to the good agreement between the temporal evolution of their velocity field and the projection of the velocity field associated with three-dimensional puffs in a frame of reference moving at the bulk velocity. We establish that the flow dynamics are relatively insensitive to the regeneration mechanisms invoked to produce near-wall streamwise vortices and that using small, unstructured background disturbances to regenerate the streamwise vortices is sufficient to capture the formation of the high- and low-speed streaks and their segregation leading to the blunting of the velocity profile characteristic of turbulent pipe flow

Travelling-waves consistent with turbulence-driven secondary flow in a square duct

We present numerically determined travelling-wave solutions for pressure-driven flow through a straight duct with a square cross-section. This family of solutions represents typical coherent structures (a staggered array of counter-rotating streamwise vortices and an associated low-speed streak) on each wall. Their streamwise average flow in the cross-sectional plane corresponds to an eight vortex pattern much alike the secondary flow found in the turbulent regime

Turbulence characteristics of the B\"{o}dewadt layer in a large enclosed rotor-stator system

A three-dimensional (3D) direct numerical simulation is combined with a laboratory study to describe the turbulent flow in an enclosed annular rotor-stator cavity characterized by a large aspect ratio G=(b-a)/h=18.32 and a small radius ratio a/b=0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate Omega under consideration is equivalent to the rotational Reynolds number Re=Omegab2/nu=9.5 x 104, where nu is the kinematic viscosity of the fluid. This corresponds to a value at which an experiment carried out at the laboratory has shown that the stator boundary layer is turbulent, whereas the rotor boundary layer is still laminar. Comparisons of the 3D computed solution with velocity measurements have given good agreement for the mean and turbulent fields. The results enhance evidence of weak turbulence at this Reynolds number, by comparing the turbulence properties with available data in the literature. An approximately self-similar boundary layer behavior is observed along the stator side. The reduction of the structural parameter a1 under the typical value 0.15 and the variation in the wall-normal direction of the different characteristic angles show that this boundary layer is three-dimensional. A quadrant analysis of conditionally averaged velocities is performed to identify the contributions of different events (ejections and sweeps) on the Reynolds shear stress producing vortical structures. The asymmetries observed in the conditionally averaged quadrant analysis are dominated by Reynolds stress-producing events in this B\"{o}dewadt layer. Moreover, case 1 vortices (with a positive wall induced velocity) are found to be the major source of generation of special strong events, in agreement with the conclusions of Lygren and Andersson.Comment: 16 page

Robust control stability using the error loop

The paper briefly formulates the error loop as a tool for designing robust stability control systems in front of structured and unstructured uncertainties. The error loop indicates that a tool for accommodating such uncertainties is the noise estimator, which is the unique feedback channel from plant to control. It is shown that the causality constraint preventing perfect cancellation of causal uncertainties (unknown disturbance), makes also control law to play a role, offering a further degree of freedom. Employing asymptotic expansions of the closed-loop transfer functions, simple, explicit design formulae derive from stability inequalities: they relate closed-loop eigenvalues to model parameter and requirements. A simple example is provided from a ball and beam plan