Heat transfer analysis of viscous incompressible fluid by combined natural convection and radiation in an open cavity

The effect of radiation on natural convection of Newtonian fluid contained in an open cavity is investigated in this study. The governing partial differential equations are solved numerically using the Alternate Direct Implicit method together with the Successive Over Relaxation method. The study is focused on studying the flow pattern and the convective and radiative heat transfer rates are studied for different values of radiation parameters namely, the optical thickness of the fluid, scattering albedo, and the Planck number. It was found that in the optically thin limit, an increase in the optical thickness of the fluid raises the temperature and radiation heat transfer of the fluid. However, a further increase in the optical thickness decreases the radiative heat transfer rate due to increase in the energy level of the fluid, which ultimately reduces the total heat transfer rate within the fluid.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Formulation of the theory of turbulence in an incompressible fluid.

"This report was prepared in connection with Space Technology Laboratories Program of General Research."

Adaptive discontinuous Galerkin methods for eigenvalue problems arising in incompressible fluid flows

In this article we consider the a posteriori error estimation and adaptive mesh refinement of discontinuous Galerkin finite element approximations of the hydrodynamic stability problem associated with the incompressible Navier-Stokes equations. Particular attention is given to the reliable error estimation of the eigenvalue problem in channel and pipe geometries. Here, computable a posteriori error bounds are derived based on employing the generalization of the standard Dual-Weighted-Residual approach, originally developed for the estimation of target functionals of the solution, to eigenvalue/stability problems. The underlying analysis consists of constructing both a dual eigenvalue problem and a dual problem for the original base solution. In this way, errors stemming from both the numerical approximation of the original nonlinear flow problem, as well as the underlying linear eigenvalue problem are correctly controlled. Numerical experiments highlighting the practical performance of the proposed a posteriori error indicator on adaptively refined computational meshes are presented.

Adaptivity and a posteriori error control for bifurcation problems II: Incompressible fluid flow in open systems with Z_2 symmetry

In this article we consider the a posteriori error estimation and adaptive mesh refinement of discontinuous Galerkin finite element approximations of the bifurcation problem associated with the steady incompressible Navier-Stokes equations. Particular attention is given to the reliable error estimation of the critical Reynolds number at which a steady pitchfork or Hopf bifurcation occurs when the underlying physical system possesses reflectional or Z_2 symmetry. Here, computable a posteriori error bounds are derived based on employing the generalization of the standard Dual-Weighted-Residual approach, originally developed for the estimation of target functionals of the solution, to bifurcation problems. Numerical experiments highlighting the practical performance of the proposed a posteriori error indicator on adaptively refined computational meshes are presented.

Double diffusive magneto-convection fluid flow in a strong cross magnetic field with uniform surface heat and mass flux

In this study, magnetohydrodynamic natural convection boundary layer flow of an electrically conducting and viscous incompressible fluid along a heated vertical flat plate with uniform heat and mass flux in the presence of strong cross magnetic field has been investigated. For smooth integrations the boundary layer equations are transformed in to a convenient dimensionless form by using stream function formulation as well as the free variable formulation. The nonsimilar parabolic partial differential equations are integrated numerically for Pr ≪1 that is appropriate for liquid metals against the local Hartmann parameter ξ . Further, asymptotic solutions are obtained near the leading edge using regular perturbation method for smaller values of ξ . Solutions for values of ξ ≫ 1 are also obtained by employing the matched asymptotic technique. The results obtained for small, large and all ξ regimes are examined in terms of shear stress, τw, rate of heat transfer, qw, and rate of mass transfer, mw, for important physical parameter. Attention has been given to the influence of Schmidt number, Sc, buoyancy ratio parameter, N and local Hartmann parameter, ξ on velocity, temperature and concentration distributions and noted that velocity and temperature of the fluid achieve their asymptotic profiles for Sc ≥ 10:0.

Uncertainty analysis of vibrational frequencies of an incompressible liquid in a rectangular tank with and without a baffle using polynomial chaos expansion

Polynomial chaos expansion (PCE) with Latin hypercube sampling (LHS) is employed for calculating the vibrational frequencies of an inviscid incompressible fluid partially filled in a rectangular tank with and without a baffle. Vibration frequencies of the coupled system are described through their projections on the PCE which uses orthogonal basis functions. PCE coefficients are evaluated using LHS. Convergence on the coefficient of variation is used to find the orthogonal polynomial basis function order which is employed in PCE. It is observed that the dispersion in the eigenvalues is more in the case of a rectangular tank with a baffle. The accuracy of the PCE method is verified with standard MCS results and is found to be more efficient.

On the transport properties of fluid-particle flow

The hydrodynamic forces acting on a solid particle in a viscous, incompressible fluid medium at low Reynolds number flow is investigated mathematically as a prerequisite to the understanding of transport processes in two-phase flow involving solid particles and fluid. Viscous interaction between a small number of spherical particles and continuous solid boundaries as well as fluid interface are analyzed under a “point-force” approximation. Non-spherical and elastic spherical particles in a simple shear flow area are then considered. Non-steady motion of a spherical particle is briefly touched upon to illustrate the transient effect of particle motion.

A macroscopic continuum description of particle-fluid flow is formulated in terms of spatial averages yielding a set of particle continuum and bulk fluid equations. Phenomenological formulas describing the transport processes in a fluid medium are extended to cases where the volume concentration of solid particles is sufficiently high to exert an important influence. Hydrodynamic forces acting on a spherical solid particle in such a system, e.g. drag, torque, rotational coupling force, and viscous collision force between streams of different sized particles moving relative to each other are obtained. Phenomenological constants, such as the shear viscosity coefficient, and the diffusion coefficient of the bulk fluid, are found as a function of the material properties of the constituents of the two-phase system and the volume concentration of solid. For transient heat conduction phenomena, it is found that the introduction of a complex conductivity for the bulk fluid permits a simple mathematical description of this otherwise complicated process. The rate of heat transfer between particle continuum and bulk fluid is also investigated by means of an Oseen-type approximation to the energy equation.

Non-relativistic conformal symmetries in fluid mechanics

The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schrodinger group, which also involves, in addition, Schrodinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally invariant relativistic theory, the recently discussed conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.

La méthode IIM pour une membrane immergée dans un fluide incompressible

La méthode IIM (Immersed Interface Method) permet d'étendre certaines méthodes numériques à des problèmes présentant des discontinuités. Elle est utilisée ici pour étudier un fluide incompressible régi par les équations de Navier-Stokes, dans lequel est immergée une membrane exerçant une force singulière. Nous utilisons une méthode de projection dans une grille de différences finies de type MAC. Une dérivation très complète des conditions de saut dans le cas où la viscosité est continue est présentée en annexe. Deux exemples numériques sont présentés : l'un sans membrane, et l'un où la membrane est immobile. Le cas général d'une membrane mobile est aussi étudié en profondeur.

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

Criteria for mixed grids in computational fluid dynamics

The finite volume method is used as a numerical method for solving the fluid flow equations. This method is appropriate to employ under structured and unstructured meshes. Mixed grids, combining both types of grids, are investigated. The coupling of different grids is done by overlapping strategy. The computational effort for the mixed grid is evaluated by the CPU-time, with different percentage of covering area of the unstructured mesh. The present scheme is tested for the driven cavity problem, where the incompressible fluid is integrated by calculating the velocity fields and computing the pressure field in each time step. Several schemes for unstructured grid are examined, and the compatibility condition is applied to check their consistency. A scheme to verify the compatibility condition for the unstructured grids is presented. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.

Fluid dynamic analysis and modelling of industrial chemical equipment

The research is aimed at contributing to the identification of reliable fully predictive Computational Fluid Dynamics (CFD) methods for the numerical simulation of equipment typically adopted in the chemical and process industries. The apparatuses selected for the investigation, specifically membrane modules, stirred vessels and fluidized beds, were characterized by a different and often complex fluid dynamic behaviour and in some cases the momentum transfer phenomena were coupled with mass transfer or multiphase interactions. Firs of all, a novel modelling approach based on CFD for the prediction of the gas separation process in membrane modules for hydrogen purification is developed. The reliability of the gas velocity field calculated numerically is assessed by comparison of the predictions with experimental velocity data collected by Particle Image Velocimetry, while the applicability of the model to properly predict the separation process under a wide range of operating conditions is assessed through a strict comparison with permeation experimental data. Then, the effect of numerical issues on the RANS-based predictions of single phase stirred tanks is analysed. The homogenisation process of a scalar tracer is also investigated and simulation results are compared to original passive tracer homogenisation curves determined with Planar Laser Induced Fluorescence. The capability of a CFD approach based on the solution of RANS equations is also investigated for describing the fluid dynamic characteristics of the dispersion of organics in water. Finally, an Eulerian-Eulerian fluid-dynamic model is used to simulate mono-disperse suspensions of Geldart A Group particles fluidized by a Newtonian incompressible fluid as well as binary segregating fluidized beds of particles differing in size and density. The results obtained under a number of different operating conditions are compared with literature experimental data and the effect of numerical uncertainties on axial segregation is also discussed.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).