Vortex penalization method for bluff body flows

In this work, a penalization method is discussed in the context of vortex methods for incompressible flows around complex geometries. In particular, we illustrate the method in two cases: the flow around a rotating blade for Reynolds numbers 1000 and 10,000 and the flow past a semi‐circular body consisting of a porous layer surrounding a rigid body at Reynolds numbers 550 and 3000. In the latter example, the results are interpreted in terms of control strategy. Copyright © 2015 John Wiley & Sons, Ltd.     

Passive control of the flow around a square cylinder using porous media

The passive control of bluff body flows using porous media is investigated by means of the penalization method. This method is used to create intermediate porous media between solid obstacles and the fluid in order to modify the boundary layer behaviour. The study covers a wide range of two‐dimensional flows from low transitional flow to fully established turbulence by direct numerical simulation of incompressible Navier–Stokes equations. A parametric study is performed to illustrate the effect of the porous layer permeability and thickness on the passive control. The numerical results reveal the ability of porous media to both regularize the flow and to reduce the drag forces up to 30%. Copyright © 2004 John Wiley & Sons, Ltd.     

Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Incompressible Viscous Fluid Flows in a Thin Spherical Shell

Linearized stability of incompressible viscous fluid flows in a thin spherical shell is studied by using the two-dimensional Navier–Stokes equations on a sphere. The stationary flow on the sphere has two singularities (a sink and a source) at the North and South poles of the sphere. We prove analytically for the linearized Navier–Stokes equations that the stationary flow is asymptotically stable. When the spherical layer is truncated between two symmetrical rings, we study eigenvalues of the linearized equations numerically by using power series solutions and show that the stationary flow remains asymptotically stable for all Reynolds numbers.      

A Two-Layer Mathematical Model of Blood Flow in Porous Constricted Blood Vessels

In this paper, we discussed a mathematical model for two-layered non-Newtonian blood flow through porous constricted blood vessels. The core region of blood flow contains the suspension of erythrocytes as non-Newtonian Casson fluid and the peripheral region contains the plasma flow as Newtonian fluid. The wall of porous constricted blood vessel configured as thin transition Brinkman layer over layered by Darcy region. The boundary of fluid layer is defined as stress jump condition of Ocha-Tapiya and Beavers–Joseph. In this paper, we obtained an analytic expression for velocity, flow rate, wall shear stress. The effect of permeability, plasma layer thickness, yield stress and shape of the constriction on velocity in core & peripheral region, wall shear stress and flow rate is discussed graphically. This is found throughout the discussion that permeability and plasma layer thickness have accountable effect on various flow parameters which gives an important observation for diseased blood vessels.     

Existence and Stability of Asymmetric Burgers Vortices

Burgers vortices are stationary solutions of the three-dimensional Navier–Stokes equations in the presence of a background straining flow. These solutions are given by explicit formulas only when the strain is axisymmetric. In this paper we consider a weakly asymmetric strain and prove in that case that non-axisymmetric vortices exist for all values of the Reynolds number. In the limit of large Reynolds numbers, we recover the asymptotic results of Moffatt, Kida & Ohkitani [11]. We also show that the asymmetric vortices are stable with respect to localized two-dimensional perturbations.     

Influence of Vertical Vibrations on the Stability of a Binary Mixture in a Horizontal Porous Layer Subjected to a Vertical Heat Flux

Macroscale three-dimensional modeling of fluid flow in a thin porous layer under unsaturated conditions is a challenging task. One major issue is that such layers do not satisfy the representative elementary volume length-scale requirement. Recently, a new approach, called reduced continua model (RCM), has been developed to describe multiphase fluid flow in a stack of thin porous layers. In that approach, flow equations are formulated in terms of thickness-averaged variables and properties. In this work, we have performed a set of experiments, where a wet \(260\hbox {-}\upmu \hbox {m}\)-thin porous layer was placed on top of a dry layer of the same material. We measured the change of average saturation with time using a single-sided low-field nuclear magnetic resonance device known as NMR-MOUSE. We have employed both RCM and the traditional Richards equation-based models to simulate our experimental results. We found that the traditional unsaturated flow model cannot simulate experimental results satisfactorily. Very close agreement was obtained by including the dynamic capillary term as postulated by Hassanizadeh and Gray in the traditional equations. The reduced continua model was found to be in good agreement with the experimental result without adding dynamic capillarity term. Moreover, the computational effort needed for RCM simulations was one order of magnitude less than that of traditional models.     

Hopf Bifurcation in Viscous Incompressible Flow Down an Inclined Plane

We provide the Hopf bifurcation theorem, which guarantees the existence of time periodic solution bifurcating from the stationary flow down an inclined plane under certain assumptions on the eigenvalues of the problem obtained by linearization around the stationary flow. Since we reduce the problem to the fixed domain, the inhomogeneous terms of reduced equations and reduced boundary conditions contain the highest derivatives. To deal with these we apply the Lyapunov–Schmidt decomposition directly.     

Solutions similaires pour un problème de couches limites en milieux poreux

We consider heat transfer from a surface embedded in an unbounded porous medium, saturated with a fluid at rest. A thermal boundary-layer approximation, based on the assumption that convection takes place in a thin layer around the heating surface, is done. The use of boundary layer technics show that we can find similarity solutions by solving a one-dimensional boundary value problem, involving a third-order nonlinear differential equation depending on a parameter. We prove existence and uniqueness results for some values of the parameter, non-existence for the other ones, and when it is possible, we construct explicit solutions.     

Flow through a thin porous layer (grid)

This paper gives an exact hydrodynamic formulation of the problem of flow though a grid, regarded as a limiting case of a thin porous layer. The formulated coupling conditions on the grid contain two characteristics: the resistance coefficient and the refractive index, which are connected by a universal relation, independent of the grid resistance law and the nature of the flow. Problems of a weakly disturbed flow of fluid passing through the grid in a plane channel are considered as an illustration. A qualitative explanation of the calculated effects is given. Smoothing grids, distributing and supporting grids, in, for instance, apparatuses with a stationary or fluidized bed, and thin porous heat-producing layers are widely used. The interaction of a flow with a smoothing grid has been investigated in the hydraulic formulation [1]. In this case the grid acted only as a source of pressure difference and was characterized by the resistance coefficient alone. Yet a grid (screen), which is a thin structured layer, will, generally speaking, also have a refractive effect on the flow. In fact, in the case of inclined incidence the flow will tend to turn in the direction of least resistance, i.e., perpendicular to the layer. It is significant that such acceleration of the flow within the layer cannot be described within the framework of the Darcy model, which assumes that the motion is inertialess. In view of this, the exact hydrodynamic formulation of the problem of flow through a grid, regarded as a surface of discontinuity, is based on a previously proposed model^* that takes inertial effects into account.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 74–80, July–August, 1978.     

Computing the force distribution on the surface of complex,deforming geometries using vortex methods and Brinkman penalization

The distribution of forces on the surface of complex, deforming geometries is an invaluable output of flow simulations. One particular example of such geometries involves self‐propelled swimmers. Surface forces can provide significant information about the flow field sensed by the swimmers and are difficult to obtain experimentally. At the same time, simulations of flow around complex, deforming shapes can be computationally prohibitive when body‐fitted grids are used. Alternatively, such simulations may use penalization techniques. Penalization methods rely on simple Cartesian grids to discretize the governing equations, which are enhanced by a penalty term to account for the boundary conditions. They have been shown to provide a robust estimation of mean quantities, such as drag and propulsion velocity, but the computation of surface force distribution remains a challenge. We present a method for determining flow‐induced forces on the surface of both rigid and deforming bodies, in simulations using remeshed vortex methods and Brinkman penalization. The pressure field is recovered from the velocity by solving a Poisson's equation using the Green's function approach, augmented with a fast multipole expansion and a tree‐code algorithm. The viscous forces are determined by evaluating the strain‐rate tensor on the surface of deforming bodies, and on a “lifted” surface in simulations involving rigid objects. We present results for benchmark flows demonstrating that we can obtain an accurate distribution of flow‐induced surface forces. The capabilities of our method are demonstrated using simulations of self‐propelled swimmers, where we obtain the pressure and shear distribution on their deforming surfaces.     

Parametric modulation in the Taylor-Couette ferrofluid flow

A parametric instability of the Taylor-Couette ferrofluid flow excited by a periodically oscillating magnetic field, has been investigated numerically. The Floquet analysis has been employed. It has been found that the modulation of the applied magnetic field affects the stability of the basic flow. The instability response has been found to be synchronous with respect to the frequency of periodically oscillating magnetic field.     

Computation of turbulent flow past an array of cylinders using a spectral method with Brinkman penalization

A major difficulty in computing engineering flows at high Reynolds number is the need for non-uniform grids adapted to solid boundaries that may be moving or changing shape. These non-uniform grids are expensive to calculate and cannot be used with the most accurate or efficient numerical schemes. We present one solution to this problem: a Brinkman (volume) penalization of the obstacle which allows an efficient pseudo-spectral method to be used to solve the Navier–Stokes equations on a Cartesian grid. Although this is the most severe test of the penalization (due to the global support of the Fourier basis), it is shown that the method still yields reasonable results. We also present an analytical solution of Stokes flow calculated using the penalization which illustrates the error and continuity properties of the approach. Work is currently underway to implement the penalization approach in a wavelet basis.     

T-RANS Simulations of Subcritical Flow with Heat Transfer Past a Circular Cylinder Surrounded by a Thin Porous Layer

We study flow and heat transfer to a cylinder in cross flow at Re = 3,900–80,000 by means of three-dimensional transient RANS (T-RANS) simulations, employing an RNG k − ε turbulence model. Both the case of a bare solid cylinder and that of a solid cylinder surrounded at some fixed distance by a thin porous layer have been studied. The latter configuration is a standard test geometry for measuring the insulating and protective performance of garments. In this geometry, the flow in the space between the solid cylinder and the porous layer is laminar but periodic, whereas the outer flow is transitional and characterized by vortex shedding in the wake of the cylinder. The results from the T-RANS simulations are validated against data from Direct Numerical Simulations and experiments. It is found that T-RANS is very well suited for simulating this type of flow. The transient nature of the flow underneath the porous layer is well reproduced, as well as the influence of vortex shedding on the heat transfer in the downstream stagnation zone. T-RANS results are found to be in much better agreement with DNS and experimental data than results from steady-state RANS.     

Analysis of Thermal Flow in a Rotating Porous U-Turn Duct Using Lattice Boltzmann Method

The lattice Boltzmann method is carried out to investigate the heat transfer enhancement in a U-turn duct which is partially filled with a porous media. The porous layer is inserted at the core of the duct and is modeled using the Brinkman–Forchheimer assumptions. In order to validate the results, first a channel flow problem without any porous layer is compared with available data. Second, the porous Couette flow and partially porous channel flow are successfully compared with the studies of other researchers. Then, fluid flow in a clear U-turn duct is studied looking carefully at the velocity, curvature and rotation effects. Finally, the effects of porous layer thickness on the rate of heat transfer and pressure drop are investigated. Parametric studies are conducted to evaluate the effects of various parameters (i.e., Reynolds number, Darcy number, rotation number), highlighting their influences on the thermo-hydrodynamics behavior of the flow. The optimum values of porous layer thickness are presented for specific flow parameters.     

Effects of Hall current and slip condition on unsteady flow of a viscous fluid due to non-coaxial rotation of a porous disk and a fluid at infinity

The unsteady hydromagnetic flow due to non-coaxial rotations of a porous disk with slip condition and a fluid at infinity has been studied on taking Hall currents into account. An exact solution of the governing equation has been obtained by the Laplace transform technique. Asymptotic solution is obtained for large time. It is found that for large time there exists a thin boundary layer near the disk. The thickness of this layer decreases with increase in either suction or magnetic parameter.     

A New Derivation of Jeffery’s Equation

In this article, we present a modern derivation of Jeffery’s equation for the motion of a small rigid body immersed in a Navier–Stokes flow, using methods of asymptotic analysis. While Jeffery’s result represents the leading order equations of a singularly perturbed flow problem involving ellipsoidal bodies, our formulation is for bodies of general shape and we also derive the equations of the next relevant order.      

A Remark on the Existence of 2-D Steady Navier–Stokes Flow in Bounded Symmetric Domain Under General Outflow Condition

Let Ω be a 2-dimensional bounded domain, symmetric with respect to the x₂-axis. The boundary has several connected components, intersecting the x₂-axis. The boundary value is symmetric with respect to the x₂-axis satisfying the general outflow condition. The existence of the symmetric solution to the steady Navier–Stokes equations was established by Amick [2] and Fujita [4]. Fujita [4] proved a key lemma concerning the solenoidal extension of the boundary value by virtual drain method. In this note, we give a different proof via elementary approach by means of the stream function.     

Convection in Hele-Shaw cell with parametric excitation

The linear stability of a thin vertical fluid layer heated from below is considered. Here the gravitational field consists of two parts: a constant part and a time-dependent part varying periodically. The time-dependent part has been expressed in Fourier series. The effect of gravity modulation on the fluid layer is examined. Using an asymptotic analysis the convective threshold has been determined. Some comparisons are made with the known results.