240,775 research outputs found
Electrical Flows, Laplacian Systems, and Faster Approximation of Maximum Flow in Undirected Graphs
We introduce a new approach to computing an approximately maximum s-t flow in
a capacitated, undirected graph. This flow is computed by solving a sequence of
electrical flow problems. Each electrical flow is given by the solution of a
system of linear equations in a Laplacian matrix, and thus may be approximately
computed in nearly-linear time.
Using this approach, we develop the fastest known algorithm for computing
approximately maximum s-t flows. For a graph having n vertices and m edges, our
algorithm computes a (1-\epsilon)-approximately maximum s-t flow in time
\tilde{O}(mn^{1/3} \epsilon^{-11/3}). A dual version of our approach computes a
(1+\epsilon)-approximately minimum s-t cut in time
\tilde{O}(m+n^{4/3}\eps^{-8/3}), which is the fastest known algorithm for this
problem as well. Previously, the best dependence on m and n was achieved by the
algorithm of Goldberg and Rao (J. ACM 1998), which can be used to compute
approximately maximum s-t flows in time \tilde{O}(m\sqrt{n}\epsilon^{-1}), and
approximately minimum s-t cuts in time \tilde{O}(m+n^{3/2}\epsilon^{-3})
The combined Lagrangian advection method
We present and test a new hybrid numerical method for simulating layerwise-two-dimensional geophysical flows. The method radically extends the original Contour-Advective Semi-Lagrangian (CASL) algorithm by combining three computational elements for the advection of general tracers (e.g. potential vorticity, water vapor, etc.): (1) a pseudospectral method for large scales, (2) Lagrangian contours for intermediate to small scales, and (3) Lagrangian particles for the representation of general forcing and dissipation. The pseudo-spectral method is both efficient and highly accurate at large scales, while contour advection is efficient and accurate at small scales, allowing one to simulate extremely finescale structure well below the basic grid scale used to represent the velocity field. The particles allow one to efficiently incorporate general forcing and dissipation
The rheology of a suspension of nearly spherical particles subject to Brownian rotations
A set of constitutive equations, valid for arbitrary linear bulk flows, is derived for a dilute suspension of nearly spherical, rigid particles which are subject to rotary Brownian couples. These constitutive equations are subsequently applied to find the resulting stress patterns for a variety of time-dependent bulk flow fields. The rheological responses are found to exhibit many of the same qualitative features as have been observed in recent experimental investigations of polymeric solutions and other complex materials
Dynamics of viscoelastic pipe flow in the maximum drag reduction limit
Polymer additives can substantially reduce the drag of turbulent flows and
the upper limit, the so called "maximum drag reduction" (MDR) asymptote is
universal, i.e. independent of the type of polymer and solvent used. Until
recently, the consensus was that, in this limit, flows are in a marginal state
where only a minimal level of turbulence activity persists. Observations in
direct numerical simulations using minimal sized channels appeared to support
this view and reported long "hibernation" periods where turbulence is
marginalized. In simulations of pipe flow we find that, indeed, with increasing
Weissenberg number (Wi), turbulence expresses long periods of hibernation if
the domain size is small. However, with increasing pipe length, the temporal
hibernation continuously alters to spatio-temporal intermittency and here the
flow consists of turbulent puffs surrounded by laminar flow. Moreover, upon an
increase in Wi, the flow fully relaminarises, in agreement with recent
experiments. At even larger Wi, a different instability is encountered causing
a drag increase towards MDR. Our findings hence link earlier minimal flow unit
simulations with recent experiments and confirm that the addition of polymers
initially suppresses Newtonian turbulence and leads to a reverse transition.
The MDR state on the other hand results from a separate instability and the
underlying dynamics corresponds to the recently proposed state of
elasto-inertial-turbulence (EIT).Comment: 18 pages, 5 figure
On self-sustained oscillations in two-dimensional compressible flow over rectangular cavities
Numerical simulations are used to investigate the resonant instabilities in two-dimensional flow past an open cavity. The compressible Navier–Stokes equations are solved directly (no turbulence model) for cavities with laminar boundary layers upstream. The computational domain is large enough to directly resolve a portion of the radiated acoustic field, which is shown to be in good visual agreement with schlieren photographs from experiments at several different Mach numbers. The results show a transition from a shear-layer mode, primarily for shorter cavities and lower Mach numbers, to a wake mode for longer cavities and higher Mach numbers. The shear-layer mode is characterized well by the acoustic feedback process described by Rossiter (1964), and disturbances in the shear layer compare well with predictions based on linear stability analysis of the Kelvin–Helmholtz mode. The wake mode is characterized instead by a large-scale vortex shedding with Strouhal number independent of Mach number. The wake mode oscillation is similar in many ways to that reported by Gharib & Roshko (1987) for incompressible flow with a laminar upstream boundary layer. Transition to wake mode occurs as the length and/or depth of the cavity becomes large compared to the upstream boundary-layer thickness, or as the Mach and/or Reynolds numbers are raised. Under these conditions, it is shown that the Kelvin–Helmholtz instability grows to sufficient strength that a strong recirculating flow is induced in the cavity. The resulting mean flow is similar to wake profiles that are absolutely unstable, and absolute instability may provide an explanation of the hydrodynamic feedback mechanism that leads to wake mode. Predictive criteria for the onset of shear-layer oscillations (from steady flow) and for the transition to wake mode are developed based on linear theory for amplification rates in the shear layer, and a simple model for the acoustic efficiency of edge scattering
- …