Flow behaviour in rearward-facing cavities

A group finite element method is described for the prediction of steady, incompressible, laminar flow in rearward-facing cavities with blowing and suction, The cavity is formed by a downstream lip extension to a reaward-facing step. The steady-state solution is obtained as the converged solution of a consistent pseudo-transient approximate factorization of the stream-function vorticity equations. Subsequently the global pressure distribution is obtained from a poisson equation for the Bernoulli energy variable. The flow over reaward-facing cavities with depth/breadth ratios of 0.56 and 1.18 is compared with the flow ever a reaward-facing step. The flowfield is dominated by a bubble of slowly recirculating fluid which grows in downstream extent with increasing Reynolds number, almost independently of lip length. The recirculation bubble extends slightly upstream of the lip and, for the longer lip, induces a very weak counter-rotating flow in the base of the cavity. The use of equal magnitude blowing and suction in the base of the cavity under the longer lip sets up a two cell structure in the cavity for clockwise blowing and suction and a single cell structure for anticlockwise blowing and suction. These additional cells appear to limit the influence of the blowing and suction on the strength and orientation of the main recirculation region downstream of the lip.     

Weakly 3D Transonic Flows in Narrow Channels

Steady transonic flows through channels so slender that the classical boundary layer approach fails are considered. The resulting viscous inviscid interaction problem for weakly three-dimensional laminar flows is formulated for perfect gases under the requirement that the channel is sufficiently narrow so that the flow outside the viscous wall layers becomes two-dimensional in the leading order approximation. The behavior of the flow upstream of a surface mounted three-dimensional obstacle will be demonstrated. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A control volume technique based on integrated RBFNs for the convection–diffusion equation

This article reports a new high‐order control‐volume discretization for the convection–diffusion equation in one and two dimensions. Diffusive fluxes at the faces of a control volume and other terms embracing the unknown field variable are all approximated using one‐dimensional integrated radial‐basis‐function networks; line integrals involving these fluxes and other integrals are evaluated using a high‐order numerical integration scheme. The accuracy of the proposed technique is investigated numerically through the solution of several linear and nonlinear test problems, including a benchmark thermally driven cavity flow. High‐order convergence solutions are obtained. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Large-Eddy Simulations of cavitation in a square surface cavity

We report on the development and application of a multiphase approach to the prediction of cavitation induced by high-speed flow over and within a square surface cavity. The approach entails employing a full cavitation model in conjunction with Large-Eddy Simulations in order to capture the initiation and development of bubble formations in turbulent-flow conditions. The incipient formation of the bubble cloud, and the flow processes of vortex shedding and shear-layer oscillations are tracked using the Volume of Fluid method. The validity of the computational approach was assessed by comparisons with experiments on cavitating flow over a hydrofoil. Application to the case of flow over and within a two-dimensional square cavity with cavitation clearly reveal the presence of traveling cavitation at the corner of the cavity trailing edge, and vortex cavitation within the cavity. It is shown that the collapse of cavitation bubbles results in an impact frequency that is higher than the frequency of the shear-layer oscillations. This implies that structural damage due to cavitation is likely to be most severe at the corner formed at the intersection of the cavity’s trailing edge and the flat surface upstream of it.     

Simplex finite element analysis of viscous incompressible flow with penalty function formulation

Viscous flow calculations are important for the determination of separated flows, recirculating flows, secondary flows and so on. This paper presents a penalty function approach for the finite element analysis of steady incompressible viscous flow. A simplex element is used with linear velocity and constant pressure in contrast to other works which usually employ higher order elements. Simplex elements yield analytical expressions for the element matrices which in turn lead to efficient solutions. Earlier works have partially indicated how constrain and lock-up problems might be avoided for simplex elements. This paper extends the earlier works by indicating the approach in detail and verifying that it is successful for several applications not discussed in the literature so far. Solution times and accuracy considerations are discussed for Couette flow, plane Poiseuille flow, a driven cavity problem, and laminar and turbulent flow over a step.     

A fluid-like element for flow and heat convection problems

Application of mixed order finite element method to solve the laminar flow and heat convection problems is presented in this study. Solutions in primitive variables are obtained using the T4/C3 mini element on unstructured grids. The T4/C3 element is an element with the fewest degrees of freedom which also meets the LBB condition. Meanwhile, with one extra node in the triangle element, the T4/C3 element is more capable to represent the realistic phenomena of flow and heat transfer by taking the central physical behavior of triangle into consideration. This numerical approach has been validated by the applications of two test cases: natural convection in a differentially heated square cavity and forced convection in a sudden expansion. The approach is then applied to predict the flow over a heated channel with the opposite walls roughened by three pairs of repeated transverse ribs. Based on the numerical validations and experiments, the good numerical performance of applying the developed model in T4/C3 elements is verified.     

Influence of power-law index on transitional Reynolds numbers for flow over a semi-circular cylinder

In this work, the governing partial differential equations (continuity and Cauchy’s momentum equations) describing the flow of power-law type non-Newtonian fluids across a semi-circular cylinder (oriented with its curved surface in the upstream direction) have been solved numerically. In particular, consideration has been given to the delineation of the critical Reynolds numbers denoting the onset of flow separation from the surface of the cylinder and the onset of the laminar vortex shedding regime. This information is germane to establish the scaling of the macroscopic characteristics like drag coefficient and Strouhal number on the governing parameters, namely, Reynolds number and power-law index. The present results clearly suggest that the transitional Reynolds numbers show a strong dependence on the type (shear-thinning and shear-thickening) of fluid behavior as well as on the severity of the shear-dependence of the viscosity. With reference to the behavior seen in Newtonian fluids, the flow remains not only attached to the surface up to higher Reynolds numbers, but shear-thinning behavior also delays the onset of the laminar vortex shedding regime. As expected, shear-thickening fluids, of course, display the opposite characteristics.     

Algebraic multigrid,mixed‐order interpolation,and incompressible fluid flow

This paper presents the results of numerical experiments on the use of equal‐order and mixed‐order interpolations in algebraic multigrid (AMG) solvers for the fully coupled equations of incompressible fluid flow. Several standard test problems are addressed for Reynolds numbers spanning the laminar range. The range of unstructured meshes spans over two orders of problem size (over one order of mesh bandwidth). Deficiencies in performance are identified for AMG based on equal‐order interpolations (both zero‐order and first‐order). They take the form of poor, fragile, mesh‐dependent convergence rates. The evidence suggests that a degraded representation of the inter‐field coupling in the coarse‐grid approximation is the cause. Mixed‐order interpolation (first‐order for the vectors, zero‐order for the scalars) is shown to address these deficiencies. Convergence is then robust, independent of the number of coarse grids and (almost) of the mesh bandwidth. The AMG algorithms used are reviewed. Copyright © 2009 John Wiley & Sons, Ltd.     

Experiments on the transition to time‐dependence in highly strained vortex flow inside a rectangular cavity

The flows in a two‐sided lid‐driven cavity with a cross‐sectional aspect ratio Γ = 1.96 is experimentally investigated by laser‐Doppler anemometry and by hot‐film measurements. When the two facing walls move in opposite directions, a robust three‐dimensional flow consisting of four cells arises for Reynolds numbers Re ≥ 275 as a result of the elliptic instability. The flow within these steady cells is measured to high precision. For Re > 825 the four cell flow becomes oscillatory. The oscillations are due to a standing wave with the same wavelength as the underlying cellular flow.     

Low and moderate Reynolds number transient flow simulations using space filtered Navier-Stokes equations

In this study, low and moderate Reynolds number flow problems in the laminar range are solved numerically with grids that do not resolve all the significant scales of motion. Spatial averaging or filtering of the Navier-Stokes equations and Taylor series approximations to the filtered advective terms are used in order to account for the effects of the unresolved or subgrid scales on the resolved scales. Numerical experiments with a transient 2-D lid driven cavity flow problem, using a penalty method Galerkin finite element code, show that this approach enhances the momentum transfer properties of the numerical solution, eliminates 2Δx type oscillations, and enables the use of coarser grids. The significance and order of the terms that describe the interaction between the resolved and the subgrid scales is studied and the success of the series approximations to these terms is demonstrated. © 1994 John Wiley & Sons, Inc.     

Method for solving the equations describing the interaction of a three-dimensional boundary layer with an outer inviscid flow

A new method is developed for solving the three-dimensional time-independent equations describing the interaction of a laminar boundary layer with an outer inviscid flow. The method also applies to the interaction of plane flows. By applying the method, the problem of the three-dimensional viscous supersonic gas flow over a roughness element (a hump and a cavity) is solved for the first time within the framework of the classical triple-deck theory. The asymptotic height of the roughness element corresponding to the nonseparated flow is determined, and separated flow patterns are constructed.     

Comparison of Some Preconditioners for the Incompressible Navier-Stokes Equations

In this paper we explore the performance of the SIMPLER, augmented Lagrangian, 'grad-div' preconditioners and their new variants for the two-by-two block systems arising in the incompressible Navier-Stokes equations. The lid-driven cavity and flow over a finite flat plate are chosen as the benchmark problems. For each problem the Reynolds number varies from a low to the limiting number for a laminar flow.     

时间空间均为二阶的新型NND差分格式

张涵信的研究表明,为了避免激波前后差分解的波动,在差分格式的改型方程中三阶导数的系数在激波上游必须是正的,而在激波下游则必须是负的.据此提出了一种新型的无波动、无自由参数耗散性的差分格式,它对时间和空间都是二阶的.证明了此格式是TVD的,而且是推广的二阶Годунов格式.在处理有激波的流场时,此格式是Lax-Wendroff格式的改进和推广.给出了若干算例,计算结果表明,此格式不仅无波动,而且具有形式紧凑、应用方便、分辨率高、稳定性准则中的Courant数较大的优点.     

Spatial decay bounds for time dependent magnetohydrodynamic geophysical flow

This article is concerned with spatial decay bounds for the time dependent magnetohydrodynamic geophysical flow in an infinite pipe when homogeneous lateral surface boundary conditions are applied. Assuming that the entrance velocity and magnetic field data are small enough and the fluid flow converges to laminar flow as the distance down the pipe tends to infinity, we derive a second order differential inequality that leads to an exponential decay estimate for the “energy” associated with the velocity and magnetic field represented by the difference between the entrance flow and fully developed laminar flow. We also show how to establish the explicit decay bounds for the total energy.     

Towards a better understanding of wall-driven square cavity flows using the lattice Boltzmann method

Wall-driven flow in square cavities has been studied extensively, yet it appears some main flow characteristics have not been fully investigated. Previous research on the classic lid-driven cavity (S1) flow has produced the critical Reynolds numbers separating the laminar steady and unsteady flows. Wall-driven cavities with two opposing walls moving at the same speed and the same (S2p) or opposite (S2a) directions have seldom been studied in the literature and no critical Reynolds numbers characterizing transitional flows have ever been investigated. After validating the LBM code for the three configurations studied, extensive numerical simulations have been undertaken to provide approximate ranges for the critical Hopf and Neimark-Sacker bifurcations for the classic and two two-sided cavity configurations. The threshold for transition to chaotic motion is also reported. The symmetries of the solutions are monitored across the various bifurcations for the two-sided wall driven cavities. The mirror-symmetry of the base solution for case S2p is lost at the Hopf bifurcation. The exact same scenario occurs with the pi-rotational symmetry of the base state for case S2a.     

Block mesh refinement for incompressible flows in curvilinear domains

A method for the solution of the Navier–Stokes equation for the prediction of flows inside domains of arbitrary shaped bounds by the use of Cartesian grids with block-refinement in space is presented. In order to avoid the complexity of the body fitted numerical grid generation procedure, we use a saw tooth method for the curvilinear geometry approximation. By using block-nested refinement, we achieved the desired geometry Cartesian approximation in order to find an accurate solution of the N–S equations. The method is applied to incompressible laminar flows and is based on a cell-centred approximation. We present the numerical simulation of the flow field for two geometries, driven cavity and stenosed tubes. The utility of the algorithm is tested by comparing the convergence characteristics and accuracy to those of the standard single grid algorithm. The Cartesian block refinement algorithm can be used in any complex curvilinear geometry simulation, to accomplish a reduction in memory requirements and the computational time effort.     

An element-by-element BICGSTAB iterative method for three-dimensional steady Navier-Stokes equations

Construction of a stabilized Galerkin upwind finite element model for steady and incompressible Navier-Stokes equations in three dimensions is the main theme of this study. In the time-independent context, the weighted residuals statement is kept biased in favor of the upstream flow direction by adding an artificial damping term of physical plausibility to the Galerkin framework. This upwind approach has significant advantage of seeking solutions free from cross-stream diffusion error. Finite element solutions have been found by mixed formulation, implemented in quadratic cubic elements which are characterized as possessing the so-called LBB (Ladyzhenskaya-Babuška-Brezzi) condition. An element-by-element BICGSTAB solution solver is intended to alleviate difficulties regarding the asymmetry and indefiniteness arising from the use of a mixed formulation for incompressible fluid flows. The developed three-dimensional finite element code is first rectified by solving a problem amenable to analytic solution. A well-known lid-driven cavity flow problem in a cubical cavity is also studied.     

Mixed convection turbulent film condensation on a sphere

The present paper reports a research on condensation heat transfer of an isothermal sphere with an external flow of vapor. The high tangential velocity of the vapor flow is determined from potential flow theory. The transition criterion of the onset turbulence has been given in the local film Reynolds number (Re_Γ). An eddy diffusivity model along with an expression by [H. Kato, N.N. Shiwaki, M. Hirota, On the turbulent heat transfer by free convection from a vertical plate, Int. J. Heat Mass Transfer, 11(1968) 1117–1125] is used to model turbulence. And the local liquid–vapor interfacial shear which occurs for high velocity vapor flow across a sphere surface is defined by the Colburn analogy. The paper then presents analytical analysis for the local dimensionless film thickness and heat transfer characteristics for the film condensation. And a comparison with those generated by previous theoretical of laminar condensation is discussed. The comparison shows the heat transfer coefficient of turbulent film condensation is higher than laminar film condensation under the high vapor velocity.     

Numerical Investigation of Particle Focusing Patterns in Laminar Pipe Flow With High Reynolds Numbers北大核心CSCD

The inertial focusing characteristics of particles in laminar flow pipes with high Re numbers were studied based on the “relative motion model”. In order to solve the problem of long pipes with high Re number flow, periodic boundary conditions were imposed on the inlet and outlet of the pipe. The research results show that the use of periodic boundary conditions can effectively reduce the computational, and the mechanical properties of particles in high Re flow can be calculated by using L=4D pipe. The difference from the low Re number is that as the Re number continues to increase,the lift force of the particles in the radial direction is no longer distributed as a parabola. The lift curve has a concave area between r+ =0.5 ~ 0.7, and there is a tendency for a new inertial focus point to appear in this section. By means of particles of a+ =1/17 for Re > 1 000, this new focus point position is solvable. In addition, in the analysis of the flow field, a secondary flow occurs around the particle, and its intensity gradually increases with the Re number and the closeness of the particle to the wall. The generation of the secondary flow affects the spatial distribution of the particle lift. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.