Showing papers on "Incompressible flow published in 2005"

Spectral/hp Element Methods for Computational Fluid Dynamics

Abstract: Introduction Fundamental concepts in one dimension Multi-dimensional expansion bases Multi-dimensional formulations Diffusion equation Advection and advection-diffusion Non-conforming elements Algorithms for incompressible flows Incompressible flow simulations:verification and validation Hyperbolic conservation laws Appendices Jacobi polynomials Gauss-Type integration Collocation differentiation Co discontinuous expansion bases Characteristic flux decomposition References Index

Locomotion of articulated bodies in a perfect fluid

Abstract: This paper is concerned with modeling the dynamics of N articulated solid bodies submerged in an ideal fluid. The model is used to analyze the locomotion of aquatic animals due to the coupling between their shape changes and the fluid dynamics in their environment. The equations of motion are obtained by making use of a two-stage reduction process which leads to significant mathematical and computational simplifications. The first reduction exploits particle relabeling symmetry: that is, the symmetry associated with the conservation of circulation for ideal, incompressible fluids. As a result, the equations of motion for the submerged solid bodies can be formulated without explicitly incorporating the fluid variables. This reduction by the fluid variables is a key difference with earlier methods, and it is appropriate since one is mainly interested in the location of the bodies, not the fluid particles. The second reduction is associated with the invariance of the dynamics under superimposed rigid motions. This invariance corresponds to the conservation of total momentum of the solid-fluid system. Due to this symmetry, the net locomotion of the solid system is realized as the sum of geometric and dynamic phases over the shape space consisting of allowable relative motions, or deformations, of the solids. In particular, reconstruction equations that govern the net locomotion at zero momentum, that is, the geometric phases, are obtained. As an illustrative example, a planar three-link mechanism is shown to propel and steer itself at zero momentum by periodically changing its shape. Two solutions are presented: one corresponds to a hydrodynamically decoupled mechanism and one is based on accurately computing the added inertias using a boundary element method. The hydrodynamically decoupled model produces smaller net motion than the more accurate model, indicating that it is important to consider the hydrodynamic interaction of the links.

Flow and Forced Convection Heat Transfer in Crossflow of Non-Newtonian Fluids over a Circular Cylinder

Abstract: The steady and incompressible flow of power-law type non-Newtonian fluids across an unconfined, heated circular cylinder is investigated numerically to determine the dependence of the individual drag components and of the heat transfer characteristics on power-law index (0.5 ≤ n ≤ 1.4), Prandtl number (1 ≤ Pr ≤ 100), and Reynolds number (5 ≤ Re ≤ 40). The momentum and energy equations are expressed in the stream function/vorticity formulation and are solved using a second-order accurate finite difference method to determine the pressure drag and frictional drag as well as the local and surface-averaged Nusselt numbers and to map the temperature field near the cylinder. The accuracy of the numerical procedure is established using previously available numerical and analytical results for momentum and heat transfer in Newtonian and power-law fluids. The results reported herein provide fundamental knowledge of the flow and heat transfer behavior for the flow of non-Newtonian fluids over a circular cylinder; t...

Mathematical Issues Concerning the Navier–Stokes Equations and Some of Its Generalizations

Abstract: This chapter primarily deals with internal, isothermal, unsteady flows of a class of incompressible fluids with both constant and shear or pressure dependent viscosity that includes the Navier–Stokes fluid as a special subclass. We begin with a description of fluids within the framework of a continuum. We then discuss various ways in which the response of a fluid can depart from that of a Navier–Stokes fluid. Next, we introduce a general thermodynamic framework that has been successful in describing the disparate response of continua that includes those of inelasticity, solid-to-solid transformation, viscoelasticity, granular materials, blood and asphalt rheology, etc. Here, it leads to a novel derivation of the constitutive equation for the Cauchy stress for fluids with constant, or shear and/or pressure, or density dependent viscosity within a full thermomechanical setting. One advantage of this approach consists in a transparent treatment of the constraint of incompressibility. We then concentrate on the mathematical analysis of three-dimensional unsteady flows of fluids with shear dependent viscosity that includes the Navier–Stokes model and Ladyzhenskaya's model as special cases. We are interested in the issues connected with mathematical self-consistency of the models, i.e., we are interested in knowing whether (1) flows exist for reasonable, but arbitrary initial data and all instants of time, (2) flows are uniquely determined, (3) the velocity is bounded and (4) the long-time behavior of all possible flows can be captured by a finite-dimensional, small (compact) set attracting all flow trajectories exponentially. For simplicity, we eliminate the choice of boundary conditions and their influence on the flows by assuming that all functions are spatially periodic with zero mean value over a periodic cell. All these results can however be extended to internal flows wherein the tangential component of the velocity satisfies Navier's slip at the boundary. Most of the results also hold for the no-slip boundary condition. While the mathematical consistency understood in the above sense for the Navier–Stokes model in three dimensions has not been established as yet, we will show that Ladyzhenskaya's model and some of its generalization enjoy all above characteristics for a certain range of parameters. We also discuss briefly further results related to generalizations of the Navier–Stokes equations.

Diffusion and Mixing in Fluid Flow

Abstract: We study enhancement of diffusive mixing on a compact Riemannian manifold by a fast incompressible flow. Our main result is a sharp description of the class of flows that make the deviation of the solution from its average arbitrarily small in an arbitrarily short time, provided that the flow amplitude is large enough. The necessary and sufficient condition on such flows is expressed naturally in terms of the spectral properties of the dynamical system associated with the flow. In particular, we find that weakly mixing flows always enhance dissipation in this sense. The proofs are based on a general criterion for the decay of the semigroup generated by an operator of the form G+iAL with a negative unbounded self-adjoint operator G, a self-adjoint operator L, and parameter A >> 1. In particular, they employ the RAGE theorem describing evolution of a quantum state belonging to the continuous spectral subspace of the hamiltonian (related to a classical theorem of Wiener on Fourier transforms of measures). Applications to quenching in reaction-diffusion equations are also considered.

On the inviscid limit for two-dimensional incompressible flow with navier friction condition

Abstract: In [Nonlinearity, 11 (1998), pp. 1625--1636], Clopeau, Mikelic, and Robert studied the inviscid limit of the two-dimensional incompressible Navier--Stokes equations in a bounded domain subject to Navier friction--type boundary conditions. They proved that the inviscid limit satisfies the incompressible Euler equations, and their result ultimately includes flows generated by bounded initial vorticities. Our purpose in this article is to adapt and, to some extent, simplify their argument in order to include pth power integrable initial vorticities, with p > 2.

Vortex-induced vibrations at subcritical Re

Abstract: Flow past a stationary cylinder becomes unstable at Re. Flow-induced vibrations of an elastically mounted cylinder, of low non-dimensional mass, is investigated at subcritical Reynolds numbers. A stabilized finite-element formulation is used to solve the incompressible flow equations and the cylinder motion in two dimensions. The cylinder is free to vibrate in both the transverse and in-line directions. It is found that, for certain natural frequencies of the spring–mass system, vortex shedding and self-excited vibrations of the cylinder are possible for Re as low as 20. Lock-in is observed in all cases. However, the mass of the oscillator plays a major role in determining the proximity of the vortex-shedding frequency to the natural frequency of the oscillator. A global linear stability analysis (LSA) for the combined flow and oscillator is carried out. The results from the LSA are in good agreement with the two-dimensional direct numerical simulations.

Hybridized globally divergence-free LDG methods. Part I: The Stokes problem

Abstract: We devise and analyze a new local discontinuous Galerkin (LDG) method for the Stokes equations of incompressible fluid flow. This optimally convergent method is obtained by using an LDG method to discretize a vorticity-velocity formulation of the Stokes equations and by applying a new hybridization to the resulting discretization. One of the main features of the hybridized method is that it provides a globally divergence-free approximate velocity without having to construct globally divergence-free finite-dimensional spaces; only elementwise divergence-free basis functions are used. Another important feature is that it has significantly less degrees of freedom than all other LDG methods in the current literature; in particular, the approximation to the pressure is only defined on the faces of the elements. On the other hand, we show that, as expected, the condition number of the Schur-complement matrix for this approximate pressure is of order h -2 in the mesh size h. Finally, we present numerical experiments that confirm the sharpness of our theoretical a priori error estimates.

Domain-decomposition method for parallel lattice Boltzmann simulation of incompressible flow in porous media.

Abstract: The lattice Boltzmann method has proven to be a promising method to simulate flow in porous media. Its practical application often relies on parallel computation because of the demand for a large domain and fine grid resolution to adequately resolve pore heterogeneity. The existing domain-decomposition methods for parallel computation usually decompose a domain into a number of subdomains first and then recover the interfaces and perform the load balance. Normally, the interface recovery and the load balance have to be performed iteratively until an acceptable load balance is achieved; this costs time. In this paper we propose a cell-based domain-decomposition method for parallel lattice Boltzmann simulation of flow in porous media. Unlike the existing methods, the cell-based method performs the load balance first to divide the total number of fluid cells into a number of groups (or subdomains), in which the difference of fluid cells in each group is either 0 or 1, depending on if the total number of fluid cells is a multiple of the processor numbers; the interfaces between the subdomains are recovered at last. The cell-based method is to recover the interfaces rather than the load balance; it does not need iteration and gives an exact load balance. The performance of the proposed method is analyzed and compared using different computer systems; the results indicate that it reaches the theoretical parallel efficiency. The method is then applied to simulate flow in a three-dimensional porous medium obtained with microfocus x-ray computed tomography to calculate the permeability, and the result shows good agreement with the experimental data.

Systematic Derivation of Jump Conditions for the Immersed Interface Method in Three-Dimensional Flow Simulation

Abstract: In this paper, we systematically derive jump conditions for the immersed interface method [SIAM J. Numer. Anal., 31 (1994), pp. 1019-1044; SIAM J. Sci. Comput., 18 (1997), pp. 709-735] to simulate three-dimensional incompressible viscous flows subject to moving surfaces. The surfaces are represented as singular forces in the Navier--Stokes equations, which give rise to discontinuities of flow quantities. The principal jump conditions across a closed surface of the velocity, the pressure, and their normal derivatives have been derived by Lai and Li [Appl. Math. Lett., 14 (2001), pp. 149-154]. In this paper, we first extend their derivation to generalized surface parametrization. Starting from the principal jump conditions, we then derive the jump conditions of all first-, second-, and third-order spatial derivatives of the velocity and the pressure. We also derive the jump conditions of first- and second-order temporal derivatives of the velocity. Using these jump conditions, the immersed interface method is applicable to the simulation of three-dimensional incompressible viscous flows subject to moving surfaces, where near the surfaces the first- and second-order spatial derivatives of the velocity and the pressure can be discretized with, respectively, third- and second-order accuracy, and the first-order temporal derivatives of the velocity can be discretized with second-order accuracy.

Numerical Simulation of Swirling Flow in Complex Hydroturbine Draft Tube Using Unsteady Statistical Turbulence Models

Abstract: A numerical method is developed for carrying out unsteady Reynolds-averaged Navier-Stokes (URANS) simulations and detached-eddy simulations (DESs) in complex 3D geometries. The method is applied to simulate incompressible swirling flow in a typical hydroturbine draft tube, which consists of a strongly curved 90° elbow and two piers. The governing equations are solved with a second-order-accurate, finite-volume, dual-time-stepping artificial compressibility approach for a Reynolds number of 1.1 million on a mesh with 1.8 million nodes. The geometrical complexities of the draft tube are handled using domain decomposition with overset (chimera) grids. Numerical simulations show that unsteady statistical turbulence models can capture very complex 3D flow phenomena dominated by geometry-induced, large-scale instabilities and unsteady coherent structures such as the onset of vortex breakdown and the formation of the unsteady rope vortex downstream of the turbine runner. Both URANS and DES appear to yield the general shape and magnitude of mean velocity profiles in reasonable agreement with measurements. Significant discrepancies among the DES and URANS predictions of the turbulence statistics are also observed in the straight downstream diffuser.

Moving Mesh Finite Element Methods for the Incompressible Navier--Stokes Equations

Abstract: This work presents the first effort in designing a moving mesh algorithm to solve the incompressible Navier--Stokes equations in the primitive variables formulation. The main difficulty in developing this moving mesh scheme is how to keep it divergence-free for the velocity field at each time level. The proposed numerical scheme extends a recent moving grid method based on harmonic mapping [R. Li, T. Tang, and P. W. Zhang, J. Comput. Phys., 170 (2001), pp. 562--588], which decouples the PDE solver and the mesh-moving algorithm. This approach requires interpolating the solution on the newly generated mesh. Designing a divergence-free-preserving interpolation algorithm is the first goal of this work. Selecting suitable monitor functions is important and is found challenging for the incompressible flow simulations, which is the second goal of this study. The performance of the moving mesh scheme is tested on the standard periodic double shear layer problem. No spurious vorticity patterns appear when even fairly coarse grids are used.

Extending the validity of squeezed-film damper models with elongations of surface dimensions

Abstract: Compact models for micromechanical squeezed-film dampers with gap sizes comparable to the surface dimensions are presented. Two different models considering both the border flow and non-uniform pressure distribution effects are first derived for small squeeze numbers. In the first 'surface extension' model the border effects are considered simply by calculating the damping with extended surface dimensions, and in the second 'border flow channel' model an additional short fictitious flow channel is placed at the damper borders. Utilizing a large amount of two-dimensional (2D) FEM simulation results by varying the damper dimensions, mainly the ratio a/h between the surface length and the air gap height, surface elongations are extracted using both elongation models. Both linear and torsional modes of motion are considered at the continuum flow regime. These results show that the 'surface extension' model is superior, since the extracted elongation Δa is almost constant (Δa = 1.3h), leading to a very simple model. Next, the rare gas effects are included in the 'surface extension' model in the slip flow regime (Knudsen number 0 4 in the linear motion and for a/h > 10 in the torsional motion. The model assumes incompressible flow and thus the maximum frequency where the models are valid is limited. In typical MEMS topologies where the elongations must be considered, this means that the models are valid below frequencies of 500 kHz. To also model rectangular 2D squeezed-film dampers, these elongations are applied directly in the surface length and width used in the compact models. Comparison with three-dimensional (3D) FEM simulations shows that the new model gives excellent results, and it extends the validity range of existing compact models. The maximum relative error of the models is smaller than 10% for a/h > 16 in the linear motion and for a/h > 16 in the torsional motion. The new surface extension model is useful in simulating both the circuit level and the system level behavior of gas-damped microelectromechanical devices with aspect ratios greater than 2 in the time and frequency domains.