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

Divakar Viswanath

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

ResearcherUnclaimed source record
math.NAphysics.flu-dynnlin.CDmath.APmath.DSphysics.ao-ph

Catalog footprint

What is connected

11works
6topics
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

11 published item(s)

preprint2015arXiv

Error Analysis of Finite Differences and the Mapping Parameter in Spectral Differentiation

The Chebyshev points are commonly used for spectral differentiation in non-periodic domains. The rounding error in the Chebyshev approximation to the $n$-the derivative increases at a rate greater than $n^{2m}$ for the $m$-th derivative. The mapping technique of Kosloff and Tal-Ezer (\emph{J. Comp. Physics}, vol. 104 (1993), p. 457-469) ameliorates this increase in rounding error. We show that the argument used to justify the choice of the mapping parameter is substantially incomplete. We analyze rounding error as well as discretization error and give a more complete argument for the choice of the mapping parameter. If the discrete cosine transform is used to compute derivatives, we show that a different choice of the mapping parameter yields greater accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intermittency at Fine Scales and Complex Singularities of Turbulent Couette Flow

Fine scales of turbulent velocity fields, beyond the inertial range and well into the dissipative range, are highly intermittent. It has been hypothesized that complex plane singularities are the principal mechanism behind fine scale intermittency. In this article, we view the velocity field of a turbulent flow as an analytic function of time. Although the function is only available for real values of time, we present a numerical technique to analytically continue the function to complex values of time, and with sufficient fidelity to locate and visualize the singularity closest to the real axis. Using this technique, we demonstrate a robust connection between temporal intermittency and the location of singularities in the complex plane.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Navier-Stokes solver using Green's functions II: spectral integration of channel flow and plane Couette flow

The Kleiser-Schumann algorithm has been widely used for the direct numerical simulation of turbulence in rectangular geometries. At the heart of the algorithm is the solution of linear systems which are tridiagonal except for one row. This note shows how to solve the Kleiser-Schumann problem using perfectly triangular matrices. An advantage is the ability to use functions in the LAPACK library. The method is used to simulate turbulence in channel flow at $Re=80,000$ (and $Re_τ=2400$) using $10^{9}$ grid points. An assessment of the length of time necessary to eliminate transient effects in the initial state is included.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accuracy and stability of inversion of power series

This article considers the numerical inversion of the power series $p(x)=1+b_{1}x+b_{2}x^{2}+\cdots$ to compute the inverse series $q(x)$ satisfying $p(x)q(x)=1$. Numerical inversion is a special case of triangular back-substitution, which has been known for its beguiling numerical stability since the classic work of Wilkinson (1961). We prove the numerical stability of inversion of power series and obtain bounds on numerical error. A range of examples show these bounds to be quite good. When $p(x)$ is a polynomial and $x=a$ is a root with $p(a)=0$, we show that root deflation via the simple division $p(x)/(x-a)$ can trigger instabilities relevant to polynomial root finding and computation of finite-difference weights. When $p(x)$ is a polynomial, the accuracy of the computed inverse $q(x)$ is connected to the pseudozeros of $p(x)$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Barycentric Hermite Interpolation

Let $z_{1},\ldots,z_{K}$ be distinct grid points. If $f_{k,0}$ is the prescribed value of a function at the grid point $z_{k}$, and $f_{k,r}$ the prescribed value of the $r$\foreignlanguage{american}{-th} derivative, for $1\leq r\leq n_{k}-1$, the Hermite interpolant is the unique polynomial of degree $N-1$ ($N=n_{1}+\cdots+n_{K}$) which interpolates the prescribed function values and function derivatives. We obtain another derivation of a method for Hermite interpolation recently proposed by Butcher et al. {[}\emph{Numerical Algorithms, vol. 56 (2011), p. 319-347}{]}. One advantage of our derivation is that it leads to an efficient method for updating the barycentric weights. If an additional derivative is prescribed at one of the interpolation points, we show how to update the barycentric coefficients using only $\mathcal{O}\left(N\right)$ operations. Even in the context of confluent Newton series, a comparably efficient and general method to update the coefficients appears not to be known. If the method is properly implemented, it computes the barycentric weights with fewer operations than other methods and has very good numerical stability even when derivatives of high order are involved. We give a partial explanation of its numerical stability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Metric Entropy and the Optimal Prediction of Chaotic Signals

Suppose we are given a time series or a signal $x(t)$ for $0\leq t\leq T$. We consider the problem of predicting the signal in the interval $T<t\leq T+t_{f}$ from a knowledge of its history and nothing more. We ask the following question: what is the largest value of $t_{f}$ for which a prediction can be made? We show that the answer to this question is contained in a fundamental result of information theory due to Wyner, Ziv, Ornstein, and Weiss. In particular, for the class of chaotic signals, the upper bound is $t_{f}\leq\log_{2}T/H$ in the limit $T\rightarrow\infty$, with $H$ being entropy in a sense that is explained in the text. If $\bigl|x(T-s)-x(t^{\ast}-s)\bigr|$ is small for $0\leq s\leqτ$, where $τ$ is of the order of a characteristic time scale, the pattern of events leading up to $t=T$ is similar to the pattern of events leading up to $t=t^{\ast}$. It is reasonable to expect $x(t^{\ast}+t_{f})$ to be a good predictor of $x(T+t_{f}).$ All existing methods for prediction use this idea in some way or the other. Unfortunately, this intuitively reasonable idea is fundamentally deficient and all existing methods fall well short of the Wyner-Ziv entropy bound on $t_{f}$. An optimal predictor should decompose the distance between the pattern of events leading up to $t=T$ and the pattern leading up to $t=t^{\ast}$ into stable and unstable components. A good match should have suitably small unstable components but will in general allow stable components which are as large as the tolerance for correct prediction. For the special case of hyperbolic toral automorphisms, we derive an optimal predictor using Pade approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Navier-Stokes solver using Green's functions I: channel flow and plane Couette flow

Numerical solvers of the incompressible Navier-Stokes equations have reproduced turbulence phenomena such as the law of the wall, the dependence of turbulence intensities on the Reynolds number, and experimentally observed properties of turbulence energy production. In this article, we begin a sequence of investigations whose eventual aim is to derive and implement numerical solvers that can reach higher Reynolds numbers than is currently possible. Every time step of a Navier-Stokes solver in effect solves a linear boundary value problem. The use of Green's functions leads to numerical solvers which are highly accurate in resolving the boundary layer, which is a source of delicate but exceedingly important physical effects at high Reynolds numbers. The use of Green's functions brings with it a need for careful quadrature rules and a reconsideration of time steppers. We derive and implement Green's function based solvers for the channel flow and plane Couette flow geometries. The solvers are validated by reproducing turbulent signals which are in good qualitative and quantitative agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectral integration of linear boundary value problems

Spectral integration was deployed by Orszag and co-workers (1977, 1980, 1981) to obtain stable and efficient solvers for the incompressible Navier-Stokes equation in rectangular geometries. Two methods in current use for channel flow and plane Couette flow, namely, Kleiser-Schumann (1980) and Kim-Moin-Moser (1977), rely on the same technique. In its current form, the technique of spectral integration, as applied to the Navier-Stokes equations, is dominated by rounding errors at higher Reynolds numbers which would otherwise be within reach. In this article, we derive a number of versions of spectral integration and explicate their properties, with a view to extending the Kleiser-Schumann and Kim-Moin-Moser algorithms to higher Reynolds numbers. More specifically, we show how spectral integration matrices that are banded, but bordered by dense rows, can be reduced to purely banded matrices. Key properties, such as the accuracy of spectral integration even when Green's functions are not resolved by the underlying grid, the accuracy of spectral integration in spite of ill-conditioning of underlying linear systems, and the accuracy of derivatives, are thoroughly explained.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite Difference Weights, Spectral Differentiation, and Superconvergence

Let $z_{1},z_{2},...,z_{N}$ be a sequence of distinct grid points. A finite difference formula approximates the $m$-th derivative $f^{(m)}(0)$ as $\sum w_{k}f(z_{k})$, with $w_{k}$ being the weights. We derive an algorithm for finding the weights $w_{k}$ which is an improvement of an algorithm of Fornberg (\emph{Mathematics of Computation}, vol. 51 (1988), p. 699-706). This algorithm uses fewer arithmetic operations than that of Fornberg by a factor of $4/(5m+5)$ while being equally accurate. The algorithm that we derive computes finite difference weights accurately even when $m$, the order of the derivative, is as high as 16. In addition, the algorithm generalizes easily to the efficient computation of spectral differentiation matrices. The order of accuracy of the finite difference formula for $f^{(m)}(0)$ with grid points $hz_{k}$, $1\leq k\leq N$, is typically $\mathcal{O}(h^{N-m})$. However, the most commonly used finite difference formulas have an order of accuracy that is higher than the typical. For instance, the centered difference approximation $(f(h)-2f(0)+f(-h))/h^{2}$ to $f"(0)$ has an order of accuracy equal to 2 not 1. Even unsymmetric finite difference formulas can exhibit such superconvergence or boosted order of accuracy, as shown by the explicit algebraic condition that we derive. If the grid points are real, we prove a basic result stating that the order of accuracy can never be boosted by more than 1.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Complex Singularities and the Lorenz Attractor

The Lorenz attractor is one of the best known examples of applied mathematics. However, much of what is known about it is a result of numerical calculations and not of mathematical analysis. As a step toward mathematical analysis, we allow the time variable in the three dimensional Lorenz system to be complex, hoping that solutions that have resisted analysis on the real line will give up their secrets in the complex plane. Knowledge of singularities being fundamental to any investigation in the complex plane, we build upon earlier work and give a complete and consistent formal development of complex singularities of the Lorenz system using {\it psi series}. The psi series contain two undetermined constants. In addition, the location of the singularity is undetermined as a consequence of the autonomous nature of the Lorenz system. We prove that the psi series converge, using a technique that is simpler and more powerful than that of Hille, thus implying a two parameter family of singular solutions of the Lorenz system. We pose three questions, answers to which may bring us closer to understanding the connection of complex singularities to Lorenz dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Exact and Asymptotic Conditions on Traveling Wave Solutions of the Navier-Stokes Equations

We derive necessary conditions that traveling wave solutions of the Navier-Stokes equations must satisfy in the pipe, Couette, and channel flow geometries. Some conditions are exact and must hold for any traveling wave solution irrespective of the Reynolds number ($Re$). Other conditions are asymptotic in the limit $Re\to\infty$. The exact conditions are likely to be useful tools in the study of transitional structures. For the pipe flow geometry, we give computations up to $Re=100000$ showing the connection of our asymptotic conditions to critical layers that accompany vortex structures at high $Re$.

0 citations0 reviewsNo code linkedOpen access