Physics - Fluid Dynamics Publications (50)


Physics - Fluid Dynamics Publications

The main objective of this article is to study the nonlinear stability and dynamic transitions of the basic (zonal) shear flows for the three-dimensional continuously stratified rotating Boussinesq model. The model equations are fundamental equations in geophysical fluid dynamics, and dynamics associated with their basic zonal shear flows play a crucial role in understanding many important geophysical fluid dynamical processes, such as the meridional overturning oceanic circulation and the geophysical baroclinic instability. In this paper, first we derive a threshold for the energy stability of the basic shear flow, and obtain a criteria for nonlinear stability in terms of the critical horizontal wavenumbers and the system parameters such as the Froude number, the Rossby number, the Prandtl number and the strength of the shear flow. Read More

The flow in a shock tube is extremely complex with dynamic multi-scale structures of sharp fronts, flow separation, and vortices due to the interaction of the shock wave, the contact surface, and the boundary layer over the side wall of the tube. Prediction and understanding of the complex fluid dynamics is of theoretical and practical importance. It is also an extremely challenging problem for numerical simulation, especially at relatively high Reynolds numbers. Read More

Direct numerical simulation is performed to study compressible, viscous flow around a circular cylinder. The present study considers two-dimensional, shock-free continuum flow by varying the Reynolds number between 20 and 100 and the freestream Mach number between 0 and 0.5. Read More

Generating anatomically realistic three-dimensional (3D) models of the human sinonasal cavity for numerical investigations of sprayed drug transport presents a host of methodological ambiguities. For example, subject-specific radiographic images used for 3D reconstructions typically exclude spray bottles. Subtracting a bottle contour from the 3D airspace and dilating the anterior nasal vestibule for nozzle placement augment the complexity of model-building. Read More

We present a theoretical and numerical study of Fourier space triad phase dynamics in one-dimensional stochastically forced Burgers equation at Reynolds number $\mathrm{Re} \approx 2.7 \times 10^4$. We demonstrate that Fourier triad phases over the inertial range display a collective behaviour characterised by intermittent periods of synchronisation and alignment, reminiscent of Kuramoto model (Kuramoto 1984) and directly related to collisions of shocks in physical space. Read More

We perform numerical simulations of a turbulent channel flow over an hyper-elastic wall. In the fluid region the flow is governed by the incompressible Navier-Stokes (NS) equations, while the solid is a neo-Hookean material satisfying the incompressible Mooney-Rivlin law. The multiphase flow is solved with a one-continuum formulation, using a monolithic velocity field for both the fluid and solid phase, which allows the use of a fully Eulerian formulation. Read More

Wall modeling is the key technology for making industrial high-Reynolds-number flows accessible to computational analysis. In this work, we present an innovative multiscale approach to hybrid RANS/LES wall modeling, which overcomes the problem of the RANS/LES transition and enables coarse meshes near the boundary. In a layer near the wall, the Navier-Stokes equations are solved for an LES and a RANS component in one single equation. Read More

We study the inertial migration of finite-size neutrally buoyant spherical particles in dilute and semi-dilute suspensions in laminar square duct flow. We perform several direct numerical simulations using an immersed boundary method to investigate the effects of the bulk Reynolds number $Re_b$, particle Reynolds number $Re_p$ and duct to particle size ratio $h/a$ at different solid volume fractions $\phi$, from very dilute conditions to $20\%$. We show that the bulk Reynolds number $Re_b$ is the key parameter in inertial migration of particles in dilute suspensions. Read More

We have researched the motion of gas in the subnanochannel with functional surface which wettability has a gradient for the fluid by using molecular dynamics simulation. The results show that the gas is driven to flow under a single heat source and without any other work or energy applied to the system. The driving source is owed to the potential gradient of the functional face which keeps the fluid running in the subnanochannel. Read More

We study turbulent channel flows of monodisperse and polydisperse suspensions of finite-size spheres by means of Direct Numerical Simulations using an immersed boundary method to account for the dispersed phase. Suspensions with 3 different Gaussian distributions of particle radii are considered (i.e. Read More

Many reduced order models are neither robust with respect to the parameter changes nor cost-effective enough for handling the nonlinear dependence of complex dynamical systems. In this study, we put forth a robust machine learning framework for projection based reduced order modeling of such nonlinear and nonstationary systems. As a demonstration, we focus on a nonlinear advection-diffusion system given by the viscous Burgers equation, which is a prototype setting of more realistic fluid dynamics applications with the same quadratic nonlinearity. Read More

We report development and application of a fluid-structure interaction (FSI) solver for compressible flows with large-scale flow-induced deformation of the structure. The FSI solver utilizes partitioned approach to strongly couple a sharp-interface immersed boundary method based flow solver with an open-source finite-element structure dynamics solver. The flow solver is based on a higher-order finite-difference method on Cartesian grid and employs ghost-cell methodology to impose boundary conditions on the immersed boundary. Read More

The project A2 of the LIMTECH Alliance aimed at a better understanding of those magnetohydrodynamic instabilities that are relevant for the generation and the action of cosmic magnetic fields. These comprise the hydromagnetic dynamo effect and various magnetically triggered flow instabilities, such as the magnetorotational instability and the Tayler instability. The project was intended to support the experimental capabilities to become available in the framework of the DREsden Sodium facility for DYNamo and thermohydraulic studies (DRESDYN). Read More

In this work, a high-order discontinuous Galerkin (DG) method is used to perform a large-eddy simulation (LES) of a subsonic isothermal jet at high Reynolds number Re D = 10^6 on a fully un-structured mesh. Its radiated acoustic field is computed using the Ffowcs Williams and Hawkings formulation. In order to assess the accuracy of the DG method, the simulation results are compared to experimental measurements and a reference simulation based on a finite volume method. Read More

We examine discrete vortex dynamics in two-dimensional flow through a network-theoretic approach. The interaction of the vortices is represented with a graph, which allows the use of network-theoretic approaches to identify key vortex-to-vortex interactions. We employ sparsification techniques on these graph representations based on spectral theory for constructing sparsified models and evaluating the dynamics of vortices in the sparsified setup. Read More

The present paper reports on our effort to characterize vortical interactions in complex fluid flows through the use of network analysis. In particular, we examine the vortex interactions in two-dimensional decaying isotropic turbulence and find that the vortical interaction network can be characterized by a weighted scale-free network. It is found that the turbulent flow network retains its scale-free behavior until the characteristic value of circulation reaches a critical value. Read More

We test the recently developed hierarchy of diffusive moment closures for gas dynamics together with the near-wall viscosity scaling on the Poiseuille flow of argon and nitrogen in a one micrometer wide channel, and compare it against the corresponding Direct Simulation Monte Carlo computations. We find that the diffusive regularized Grad equations with viscosity scaling provide the most accurate approximation to the benchmark DSMC results. At the same time, the conventional Navier-Stokes equations without the near-wall viscosity scaling are found to be the least accurate among the tested closures. Read More

The need for improved engine efficiencies has motivated the development of high-pressure combustion systems, in which operating conditions achieve and exceed critical conditions. Associated with these conditions are large thermodynamic gradients and strong variations in transport properties as the fluid undergoes mixing and phase transition. Accurately simulating these real-fluid environments remains a main challenge. Read More

We analyze the angular dynamics of small triaxial ellipsoids in a viscous shear flow subject to weak thermal noise. By numerically integrating an overdamped angular Langevin equation, we find the steady angular probability distribution for a range of triaxial particle shapes. From this distribution we compute the intrinsic viscosity of a dilute suspension of triaxial particles. Read More

Undesired wave reflections, which occur at domain boundaries in flow simulations with free-surface waves, can be minimized by applying source terms in the vicinity of the boundary to damp the waves. Examples of such approaches are absorbing layers, damping zones, forcing zones, relaxation zones and sponge layers. A problem with these approaches is that the effectivity of the wave damping depends on the parameters in the source term functions, which are case-dependent and must be adjusted to the wave. Read More

A theory is presented which analytically predicts the reflection coefficients for absorbing layers. In this manner, the absorbing layer parameters can be tuned before running the actual simulations, to ensure optimum absorption for every case. The validity of the theory is demonstrated in finite-volume-based flow simulations for pressure waves in liquid water and in an ideal gas. Read More

We report that the lift force on a single large particle segregating in a monodisperse dense granular flow is correlated with a downstream velocity lag. This correlation suggests a viscous-inertial origin for the lift force, similar to the Saffman lift force in (micro) fluids. This insight is relevant for modelling of particle-size segregation and our approach opens up a new avenue for the numerical study of granular flows. Read More

We give the first correction to the suspension viscosity due to fluid elasticity for a dilute suspension of spheres in a viscoelastic medium. Our perturbation theory is valid to $O(\mathrm{Wi}^2)$ in the Weissenberg number $\mathrm{Wi}=\dot\gamma \lambda$, where $\dot\gamma$ is the typical magnitude of the suspension velocity gradient, and $\lambda$ is the relaxation time of the viscoelastic fluid. For shear flow we find that the suspension shear-thickens due to elastic stretching in strain 'hot spots' near the particle, despite the fact that the stress inside the particles decreases relative to the Newtonian case. Read More

We examine knotted solutions, the most simple of which is the "Hopfion", from the point of view of relations between electromagnetism and ideal fluid dynamics. A map between fluid dynamics and electromagnetism works for initial conditions or for linear perturbations, allowing us to find new knotted fluid solutions. Knotted solutions are also found to to be solutions of nonlinear generalizations of electromagnetism, and of quantum-corrected actions for electromagnetism coupled to other modes. Read More

Hammack & Segur (1978) conducted a series of surface water-wave experiments in which the evolution of long waves of depression was measured and studied. This present work compares time series from these experiments with predictions from numerical simulations of the KdV, Serre, and several unidirectional and bidirectional Whitham-type equations. These comparisons show that accurate predictions come from models that contain an accurate reproduction of the Euler phase velocity, sufficient nonlinearity, and surface tension effects. Read More

The velocity profile in a water bridge is reanalyzed. Assuming hypothetically that the bulk charge has a radial distribution, a surface potential is formed that is analogous to the Zeta potential. The Navier Stokes equation is solved, neglecting the convective term; then, analytically and for special field and potential ranges, a sign change of the total mass flow is reported caused by the radial charge distribution. Read More

We develop a theory of three-dimensional slow Rossby waves in rotating spherical density stratified convection. The excited by a non-axisymmetric instability, slow Rossby waves with frequency that is much smaller than the rotating frequency, interact with the density stratified convection and the inertial waves. The density stratification is taken into account using the anelastic approximation for very low-Mach-number flows. Read More

In the first part of this paper we establish a uniqueness result for continuity equations with velocity field whose derivative can be represented by a singular integral operator of an $L^1$ function, extending the Lagrangian theory in \cite{BouchutCrippa13}. The proof is based on a combination of a stability estimate via optimal transport techniques developed in \cite{Seis16a} and some tools from harmonic analysis introduced in \cite{BouchutCrippa13}. In the second part of the paper, we address a question that arose in \cite{FilhoMazzucatoNussenzveig06}, namely whether 2D Euler solutions obtained via vanishing viscosity are renormalized (in the sense of DiPerna and Lions) when the initial data has low integrability. Read More


This paper presents an experimental study of different instability scenarios in a parallelogram-shaped internal wave attractor in a trapezoidal domain filled with a uniformly stratified fluid. Energy is injected into the system via the oscillatory motion of a vertical wall of the trapezoidal domain. Whole-field velocity measurements are performed with the conventional PIV technique. Read More

The surface shear viscosity of an insoluble surfactant monolayer often depends strongly on its surface pressure. Here, we show that a particle moving within a bounded monolayer breaks the kinematic reversibility of low-Reynolds-number flows. The Lorentz reciprocal theorem allows such irreversibilities to be computed without solving the full nonlinear equations, giving the leading-order contribution of surface-pressure-dependent surface viscosity. Read More

Small drops impinging angularly on thin flowing soap films frequently demonstrate the rare emergence of bulk elastic effects working in-tandem with the more common-place hydrodynamic interactions. Three collision regimes are observable: (a) drop piercing through the film, (b) it coalescing with the flow, and (c) it bouncing off the film surface. During impact, the drop deforms along with a bulk elastic deformation of the film. Read More

We present an extensive numerical comparison of a family of balance models appropriate to the semi-geostrophic limit of the rotating shallow water equations, and derived by variational asymptotics in Oliver (2006) for small Rossby numbers ${\mathrm{Ro}}$. This family of generalized large-scale semi-geostrophic (GLSG) models contains the $L_1$-model introduced by Salmon (1983) as a special case. We use these models to produce balanced initial states for the full shallow water equations. Read More

For the purpose of Uncertainty Quantification (UQ) of Reynolds-Averaged Navier-Stokes closures, we introduce a framework in which perturbations in the eigenvalues of the anisotropy tensor are made in order to bound a Quantity-of-Interest based on limiting states of turbulence. To make the perturbations representative of local flow features, we introduce two additional transport equations for linear combinations of these aforementioned eigenvalues. The location, magnitude and direction of the eigenvalue perturbations are now governed by the model transport equations. Read More

Solutions of the forced Navier-Stokes equation have been conjectured to thermalize at scales larger than the forcing scale, similar to an absolute equilibrium obtained for the spectrally-truncated Euler equation. Using direct numeric simulations of Taylor-Green flows and general-periodic helical flows, we present results on the probability density function, energy spectrum, auto-correlation function and correlation time that compare the two systems. In the case of highly helical flows, we derive an analytic expression describing the correlation time for the absolute equilibrium of helical flows that is different from the $E^{-1/2}k^{-1}$-scaling law of weakly helical flows. Read More

This paper presents the development of an Adaptive Algebraic Multiscale Solver for Compressible flow (C-AMS) in heterogeneous porous media. Similar to the recently developed AMS for incompressible (linear) flows [Wang et al., JCP, 2014], C-AMS operates by defining primal and dual-coarse blocks on top of the fine-scale grid. Read More

This paper introduces an Algebraic MultiScale method for simulation of flow in heterogeneous porous media with embedded discrete Fractures (F-AMS). First, multiscale coarse grids are independently constructed for both porous matrix and fracture networks. Then, a map between coarse- and fine-scale is obtained by algebraically computing basis functions with local support. Read More

The turbulent Rayleigh--Taylor system in a rotating reference frame is investigated by direct numerical simulations within the Oberbeck-Boussinesq approximation. On the basis of theoretical arguments, supported by our simulations, we show that the Rossby number decreases in time, and therefore the Coriolis force becomes more important as the system evolves and produces many effects on Rayleigh--Taylor turbulence. We find that rotation reduces the intensity of turbulent velocity fluctuations and therefore the growth rate of the temperature mixing layer. Read More

Two well-known turbulence models to describe the inertial and dissipative ranges simultaneously are by Pao~[Phys. Fluids {\bf 8}, 1063 (1965)] and Pope~[{\em Turbulent Flows.} Cambridge University Press, 2000]. Read More

Sound waves on a fluid stream, in a de Laval nozzle, are shown to correspond to quasinormal modes emitted by black holes that are physical solutions in a quadratic curvature gravity with cosmological constant. Sound waves patterns in transsonic regimes at a laboratory are employed here to provide experimental data regarding generalized theories of gravity, comprised by the exact de Sitter-like solution and a perturbative solution around the Schwarzschild-de Sitter standard solution. Using the classical tests of General Relativity to bound free parameters in these solutions, acoustic perturbations on fluid flows in nozzles are then regarded to study quasinormal modes of these black holes solutions, providing deviations of the de Laval nozzle cross-sectional area, when compared to the Schwarzschild solution. Read More

Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. The objective of this work is to develop and disseminate an open source LES tool for low-Mach number turbulent combustion using the OpenFOAM framework. Read More

We study the Rolie-Poly model for entangled polymers, using a singular perturbation analysis for the limit of large relaxation time. In this limit, it is shown that the model displays the characteristic features of thixotropic yield stress fluids, including yield stress hysteresis, delayed yielding and long term persistence of a decreased viscosity after cessation of flow. We focus on the startup and cessation of shear flow. Read More

We present an analytical study of oscillatory laminar shear flow over a compliant viscoelastic layer on a rigid base. This problem relates to oscillating blood flow in viscoelastic vessels. The deeper motivation for this study, however, is the possible use of compliant coatings for turbulent drag reduction. Read More

We give a survey of recent results on weak-strong uniqueness for compressible and incompressible Euler and Navier-Stokes equations, and also make some new observations. The importance of the weak-strong uniqueness principle stems, on the one hand, from the instances of non-uniqueness for the Euler equations exhibited in the past years; and on the other hand from the question of convergence of singular limits, for which weak-strong uniqueness represents an elegant tool. Read More

We find steady channel flows that are locally optimal for transferring heat from fixed-temperature walls, under the constraint of a fixed rate of viscous dissipation (enstrophy = $Pe^2$), also the power needed to pump the fluid through the channel. We generate the optima with net flux as a continuation parameter, starting from parabolic (Poiseuille) flow, the unique optimum at maximum net flux. Decreasing the flux, we eventually reach optimal flows that concentrate the enstrophy in boundary layers of thickness $\sim Pe^{-2/5}$ at the channel walls, and have a uniform flow with speed $\sim Pe^{4/5}$ outside the boundary layers. Read More

The unsteady characteristics of the flow over thick flatback airfoils have been investigated by means of CFD calculations. Sandia airfoils which have 35% maximum thickness with three different trailing edge thicknesses were selected. The calculations provided good results compared with available experimental data with regard to the lift curve and the impact of trailing edge thickness. Read More

CFD studies have been performed to demonstrate the effects of rotation on large wind turbine blades. The 10 MW blade of the Advanced Aerodynamic Tools for Large Rotors (AVATAR) project was selected in the present investigation. In order to demonstrate the rotational augmentation, lift coefficients of the blade sections were compared with the two-dimensional simulations at consistent inflow conditions. Read More

It is shown that the addition of small amounts of microscopic rods in a viscous fluid at low Reynolds number causes a significant increase of the flow resistance. Numerical simulations of the dynamics of the solution reveal that this phenomenon is associated to a transition from laminar to chaotic flow. Polymer stresses give rise to flow instabilities which, in turn, perturb the alignment of the rods. Read More

We propose a multi-layer approach to simulate hyperpycnal and hypopycnal plumes in flows with free surface. The model allows to compute the vertical profile of the horizontal and the vertical components of the velocity of the fluid flow. The model can describe as well the vertical profile of the sediment concentration and the velocity components of each one of the sediment species that form the turbidity current. Read More

In this paper, we investigate both experimentally and theoretically the dynamics of a liquid plug driven by a cyclic forcing inside a cylindrical rigid capillary tube. First, it is shown that depending on the type of forcing (flow rate or pressure cycle), the dynamics of the liquid plug can either be stable and periodic, or conversely accelerative and eventually leading to the plug rupture. In the latter case, we identify the source of the instability to be a combination of a flow memory resulting from liquid film deposition on the walls and a lubrication effect, i. Read More

It is shown that a channel flow of a dilute polymer solution between two widely spaced cylinders hindering the flow is an important paradigm of an unbounded flow in the case in which the channel wall is located sufficiently far from the cylinders. The quantitative characterization of instabilities in a creeping viscoelastic channel flow between two widely spaced cylinders reveals two elastically driven transitions, which are associated with the breaking of time-reversal and mirror symmetries: Hopf and forward bifurcations described by two order parameters $\mbox{v}_{rms}$ and $\bar{\omega}$, respectively. We suggest that a decrease of the normalized distance between the obstacles leads to a collapse of the two bifurcations into a codimension-2 point, a situation general for many non-equilibrium systems. Read More