Combinatorial Continuous Maximal Flows

Maximum flow (and minimum cut) algorithms have had a strong impact on computer vision. In particular, graph cuts algorithms provide a mechanism for the discrete optimization of an energy functional which has been used in a variety of applications such as image segmentation, stereo, image stitching and texture synthesis. Algorithms based on the classical formulation of max-flow defined on a graph are known to exhibit metrication artefacts in the solution. Therefore, a recent trend has been to instead employ a spatially continuous maximum flow (or the dual min-cut problem) in these same applications to produce solutions with no metrication errors. However, known fast continuous max-flow algorithms have no stopping criteria or have not been proved to converge. In this work, we revisit the continuous max-flow problem and show that the analogous discrete formulation is different from the classical max-flow problem. We then apply an appropriate combinatorial optimization technique to this combinatorial continuous max-flow CCMF problem to find a null-divergence solution that exhibits no metrication artefacts and may be solved exactly by a fast, efficient algorithm with provable convergence. Finally, by exhibiting the dual problem of our CCMF formulation, we clarify the fact, already proved by Nozawa in the continuous setting, that the max-flow and the total variation problems are not always equivalent.Comment: 26 page

Optimized multi-objective design of herringbone micromixers

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.A design method which systematically integrates Computational Fluids Dynamics (CFD) with an optimization scheme based on the use of the techniques Design of Experiments (DOE), Function Approximation technique (FA) and Multi-Objective Genetic Algorithm (MOGA), has been applied to the shape optimization of the staggered herringbone micromixer (SHM) at different Reynolds numbers. To quantify the mixing intensity in the mixer a Mixing index is defined on the basis of the intensity of segregation of the mass concentration on the outlet section. Four geometric parameters, i.e., aspect ratio of the mixing channel, ratio of groove depth to channel height, ratio of groove width to groove pitch and the asymmetry factor (offset) of groove, are the design variables selected for optimization. The mixing index at the outlet section and the pressure drop in the mixing channel are the performance criteria used as objective functions. The Pareto front with the optimum trade-offs, maximum mixing index with minimum pressure drop, is obtained. Experiments for qualitative and quantitative validation have been implemented.This study is supported by the Dorothy Hodgkin Postgraduate Award (DHPA) of the Engineering and Physical Sciences Research Council (EPSRC) of United Kingdom and Ebara Research Co. Ltd. of Japan

Localizing Region-Based Active Contours

©2008 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.DOI: 10.1109/TIP.2008.2004611In this paper, we propose a natural framework that allows any region-based segmentation energy to be re-formulated in a local way. We consider local rather than global image statistics and evolve a contour based on local information. Localized contours are capable of segmenting objects with heterogeneous feature profiles that would be difficult to capture correctly using a standard global method. The presented technique is versatile enough to be used with any global region-based active contour energy and instill in it the benefits of localization. We describe this framework and demonstrate the localization of three well-known energies in order to illustrate how our framework can be applied to any energy. We then compare each localized energy to its global counterpart to show the improvements that can be achieved. Next, an in-depth study of the behaviors of these energies in response to the degree of localization is given. Finally, we show results on challenging images to illustrate the robust and accurate segmentations that are possible with this new class of active contour models

Cluster-based feedback control of turbulent post-stall separated flows

We propose a novel model-free self-learning cluster-based control strategy for general nonlinear feedback flow control technique, benchmarked for high-fidelity simulations of post-stall separated flows over an airfoil. The present approach partitions the flow trajectories (force measurements) into clusters, which correspond to characteristic coarse-grained phases in a low-dimensional feature space. A feedback control law is then sought for each cluster state through iterative evaluation and downhill simplex search to minimize power consumption in flight. Unsupervised clustering of the flow trajectories for in-situ learning and optimization of coarse-grained control laws are implemented in an automated manner as key enablers. Re-routing the flow trajectories, the optimized control laws shift the cluster populations to the aerodynamically favorable states. Utilizing limited number of sensor measurements for both clustering and optimization, these feedback laws were determined in only $O(10)$ iterations. The objective of the present work is not necessarily to suppress flow separation but to minimize the desired cost function to achieve enhanced aerodynamic performance. The present control approach is applied to the control of two and three-dimensional separated flows over a NACA 0012 airfoil with large-eddy simulations at an angle of attack of $9^\circ$, Reynolds number $Re = 23,000$ and free-stream Mach number $M_\infty = 0.3$. The optimized control laws effectively minimize the flight power consumption enabling the flows to reach a low-drag state. The present work aims to address the challenges associated with adaptive feedback control design for turbulent separated flows at moderate Reynolds number.Comment: 32 pages, 18 figure

Bypass to Turbulence in Hydrodynamic Accretion Disks: An Eigenvalue Approach

Cold accretion disks such as those in star-forming systems, quiescent cataclysmic variables, and some active galactic nuclei, are expected to have neutral gas which does not couple well to magnetic fields. The turbulent viscosity in such disks must be hydrodynamic in origin, not magnetohydrodynamic. We investigate the growth of hydrodynamic perturbations in a linear shear flow sandwiched between two parallel walls. The unperturbed flow is similar to plane Couette flow but with a Coriolis force included. Although there are no exponentially growing eigenmodes in this system, nevertheless, because of the non-normal nature of the eigenmodes, it is possible to have a large transient growth in the energy of perturbations. For a constant angular momentum disk, we find that the perturbation with maximum growth has a wave-vector in the vertical direction. The energy grows by more than a factor of 100 for a Reynolds number R=300 and more than a factor of 1000 for R=1000. Turbulence can be easily excited in such a disk, as found in previous numerical simulations. For a Keplerian disk, on the other hand, similar vertical perturbations grow by no more than a factor of 4, explaining why the same simulations did not find turbulence in this system. However, certain other two-dimensional perturbations with no vertical structure do exhibit modest growth. For the optimum two-dimensional perturbation, the energy grows by a factor of ~100 for R~10^4.5 and by a factor of 1000 for R~10^6. It is conceivable that these two-dimensional disturbances might lead to self-sustained turbulence. The Reynolds numbers of cold astrophysical disks are much larger even than 10^6, therefore, hydrodynamic turbulence may be possible in disks.Comment: 39 pages including 9 figures; Final version to appear in The Astrophysical Journa

Convex optimization of contour deformations

We discuss various formal aspects of contour deformations used to alleviate sign problems; most importantly, relating these contour deformations to a certain convex optimization problem. As a consequence of this connection we describe a general method for proving upper bounds on the average phase achievable by the contour deformation method. Using this method we show that Abelian lattice Yang-Mills in two spacetime dimensions possesses, for many values of the complex coupling, an exponential sign problem that cannot be removed via any contour deformation.Comment: 13 pages; comments welcom