Flow stabilization with active hydrodynamic cloaks

We demonstrate that fluid flow cloaking solutions, based on active hydrodynamic metamaterials, exist for two-dimensional flows past a cylinder in a wide range of Reynolds numbers (Re's), up to approximately 200. Within the framework of the classical Brinkman equation for homogenized porous flow, we demonstrate using two different methods that such cloaked flows can be dynamically stable for Re's in the range of 5-119. The first highly efficient method is based on a linearization of the Brinkman-Navier-Stokes equation and finding the eigenfrequencies of the least stable eigenperturbations; the second method is a direct numerical integration in the time domain. We show that, by suppressing the von Kármán vortex street in the weakly turbulent wake, porous flow cloaks can raise the critical Reynolds number up to about 120 or five times greater than for a bare uncloaked cylinder. © 2012 American Physical Society.

A theory of hypoellipticity and unique ergodicity for semilinear stochastic PDEs

We present a theory of hypoellipticity and unique ergodicity for semilinear parabolic stochastic PDEs with "polynomial" nonlinearities and additive noise, considered as abstract evolution equations in some Hilbert space. It is shown that if Hörmander's bracket condition holds at every point of this Hilbert space, then a lower bound on the Malliavin covariance operatorμt can be obtained. Informally, this bound can be read as "Fix any finite-dimensional projection on a subspace of sufficiently regular functions. Then the eigenfunctions of μt with small eigenvalues have only a very small component in the image of Π." We also show how to use a priori bounds on the solutions to the equation to obtain good control on the dependency of the bounds on the Malliavin matrix on the initial condition. These bounds are sufficient in many cases to obtain the asymptotic strong Feller property introduced in [HM06]. One of the main novel technical tools is an almost sure bound from below on the size of "Wiener polynomials," where the coefficients are possibly non-adapted stochastic processes satisfying a Lips chitz condition. By exploiting the polynomial structure of the equations, this result can be used to replace Norris' lemma, which is unavailable in the present context. We conclude by showing that the two-dimensional stochastic Navier-Stokes equations and a large class of reaction-diffusion equations fit the framework of our theory.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Inline holographic coherent anti-Stokes Raman microscopy.

We demonstrate a simple approach for inline holographic coherent anti-Stokes Raman scattering (CARS) microscopy, in which a layer of uniform nonlinear medium is placed in front of a specimen to be imaged. The reference wave created by four-wave mixing in the nonlinear medium can interfere with the CARS signal generated in the specimen to result in an inline hologram. We experimentally and theoretically investigate the inline CARS holography and show that it has chemical selectivity and can allow for three-dimensional imaging.

Dimensional dependence of the Stokes-Einstein relation and its violation

We generalize to higher spatial dimensions the Stokes-Einstein relation (SER) as well as the leading correction to diffusivity in finite systems with periodic boundary conditions, and validate these results with numerical simulations. We then investigate the evolution of the high-density SER violation with dimension in simple hard sphere glass formers. The analysis suggests that this SER violation disappears around dimension du = 8, above which it is not observed. The critical exponent associated with the violation appears to evolve linearly in 8 - d, below d = 8, as predicted by Biroli and Bouchaud [J. Phys.: Condens. Matter 19, 205101 (2007)], but the linear coefficient is not consistent with the prediction. The SER violation with d establishes a new benchmark for theory, and its complete description remains an open problem. © 2013 AIP Publishing LLC.