High‐order boundary conditions for linearized shallow water equations with stratification,dispersion and advection

The two‐dimensional linearized shallow water equations are considered in unbounded domains with density stratification. Wave dispersion and advection effects are also taken into account. The infinite domain is truncated via a rectangular artificial boundary ??, and a high‐order open boundary condition (OBC) is imposed on ??. Then the problem is solved numerically in the finite domain bounded by ??. A recently developed boundary scheme is employed, which is based on a reformulation of the sequence of OBCs originally proposed by Higdon. The OBCs can easily be used up to any desired order. They are incorporated here in a finite difference scheme. Numerical examples are used to demonstrate the performance and advantages of the computational method, with an emphasis is on the effect of stratification. Published in 2004 by John Wiley & Sons, Ltd.     

Poiseuille-Bénard流的出口边界条件及其谱元法计算

研究二维矩形管道中底部加热的不可压缩Ｐｏｉｓｅｕｉｌｌｅ－Ｂｅｎａｒｄ流的谱元法数值计算问题．讨论各种不同的出口边界条件的处理及其对谱元法数值模拟的影响．通过干扰区、干扰幅度和计算时间的比较,确定比较理想的出口边界条件．     

Comparison of several open boundary numerical treatments for laminar recirculating flows

Several open boundary conditions (OBCs) are compared and evaluated in the framework of the SIMPLE algorithm using staggered and non-staggered grid systems. The benchmark laminar flow test cases used for the OBC evaluation are Poiseuille-Benard flow in a channel and stratified backward-facing step flow. The investigated OBCs are linear explicit step space extrapolation, Orlanski's monochromatic wave, and pressure extrapolation. Orlanski's and pressure extrapolation open boundary treatment for unsteady and steady flows, respectively, yield little reflection and has proved to be adequate for engineering calculations.     

Parameter estimation for a three‐dimensional numerical barotropic tidal model with adjoint method

The parameters of a three‐dimensional (3‐D) barotropic tidal model are estimated using the adjoint method. The mode splitting technique is employed in both forward and adjoint models. In the external mode, the alternating direction implicit method is used to discretize the two‐dimensional depth‐averaged equations and a semi‐implicit scheme is used for the 3‐D internal mode computations. In this model the bottom friction is expressed in terms of bottom velocity which is different from the previous works. Besides, the bottom friction coefficients (BFCs) are supposed to be spatially varying, i.e. the BFC at some grid points are selected as the independent BFC, while the BFC at the other grid points can be obtained through linear interpolation with these independent BFCs. On the basis of the simulation of M₂ tide in the Bohai and North Yellow Seas (BNYS), twin experiments are carried out to invert the prescribed distributions of model parameters. The parameters inverted are the Fourier coefficients of open boundary conditions (OBCs), the BFC and the vertical eddy viscosity profiles. In these twin experiments, the real topography of BNYS is installed. The ‘observations’ are produced by the tidal model and recorded at the position of TOPEX/Poseidon altimeter data, tidal gauge data and current data. The experiments discuss the influence of initial guesses, model errors and data number. The inversion has obtained satisfactory results and the prescribed distributions have been successfully inverted. The results indicate that the inversion of BFC is more sensitive to data error than that of OBC and the vertical eddy viscosity profiles. Copyright © 2007 John Wiley & Sons, Ltd.     

A three‐dimensional numerical internal tidal model involving adjoint method

A three‐dimensional internal tidal model involving the adjoint method is constructed based on the nonlinear, time‐dependent, free‐surface hydrodynamic equations in spherical coordinates horizontally, and isopycnic coordinates vertically, subject to the hydrostatic approximations. This model consists of two submodels: the forward model is used for the simulation of internal tides, while the adjoint model is used for optimization of modal parameters. Mode splitting technique is employed in both forward and adjoint models. In this model, the adjoint method is employed to estimate model parameters by assimilating the interior observations. As a preliminary feasibility study, a set of ideal experiments with the model‐generated pseudo‐observations of surface currents are performed to invert the open boundary conditions (OBCs). In the ideal experiments, 14 kinds of bottom topographies and six kinds of predetermined distributions of OBCs are considered to examine their influence on experiment results. The inversion obtained satisfying results and all the predetermined distributions were successfully inverted. Analysis of results suggests the following: in the case where the spatial variation of the OBC distribution is great or the open boundary is close to a rough topography, the results will be comparatively poor, but still satisfactory; both the tidal elevations and currents can be simulated very accurately with the surface currents at several observation points; the assimilation precision could be reliable and able to reflect both of the inversion and simulation results in the whole field. The performance and results of ideal experiments give a preliminary indication that the construction of this model is successful. Copyright © 2011 John Wiley & Sons, Ltd.     

Résumé and remarks on the open boundary condition minisymposium

The incompressible Navier-Stokes equations—and their thermal convection and stratified flow analogue, the Boussinesq equations—possess solutions in bounded domains only when appropriate/legitimate boundary conditions (BCs) are appended at all points on the domain boundary. When the boundary—or, more commonly, a portion of it—is not endowed with a Dirichlet BC, we are faced with selecting what are called open boundary conditions (OBCs), because the fluid may presumably enter or leave the domain through such boundaries. The two minisymposia on OBCs that are summarized in this paper had the objective of finding the best OBCs for a small subset of two-dimensional test problems. This objective, which of course is not really well-defined, was not met (we believe), but the contributions obtained probably raised many more questions/issues than were resolved—notable among them being the advent of a new class of OBCs that we call FBCs (fuzzy boundary conditions).     

On the critical problem of F. D. pressure treatment for laminar flows confined by permeable walls

In the present work the viscous (low Reynolds) flow in plane ducts confined by permeable walls has been studied. A simple model of the filtrating walls has been used, with the normal velocity component proportional to the pressure jump across the wall, resulting in a non-standard boundary value Navier-Stokes problem. A critical analysis of the appropriate boundary condition and pressure problem has led to the conclusions of employing a simple explicit finite volume approach, and of avoiding the use of higher order finite difference schemes. In this paper a special emphasis on the structure of the involved computational matrices has been given to illustrate the chosen algorithm. The latter yields a steady state solution that is second order accurate in space, and it has an accuracy in time of order ≤ Δt (the time step), due to the explicit treatment of the velocity boundary conditions along the membrane. The model has been tested to study the effects of the inlet/outlet conditions, Reynolds number and filtrating wall constant.     

The solution of viscous incompressible jet flows using non‐staggered boundary fitted co‐ordinate methods

A new approach for the solution of the steady incompressible Navier–Stokes equations in a domain bounded in part by a free surface is presented. The procedure is based on the finite difference technique, with the non‐staggered grid fractional step method used to solve the flow equations written in terms of primitive variables. The physical domain is transformed to a rectangle by means of a numerical mapping technique. In order to design an effective free solution scheme, we distinguish between flows dominated by surface tension and those dominated by inertia and viscosity. When the surface tension effect is insignificant we used the kinematic condition to update the surface; whereas, in the opposite case, we used the normal stress condition to obtain the free surface boundary. Results obtained with the improved boundary conditions for a plane Newtonian jet are found to compare well with the available two‐dimensional numerical solutions for Reynolds numbers, up to Re=100, and Capillary numbers in the range of 0≤Ca<1000. Copyright © 2001 John Wiley & Sons, Ltd.     

FLOW VISUALIZATION AROUND A CIRCULAR CYLINDER NEAR TO A PLANE WALL

Flow visualization, particle image velocimetry and hot-film anemometry have been employed to study the fluid flow around a circular cylinder near to a plane wall for Reynolds numbers, based on cylinder diameter, between 1200 and 4960. The effect of changing the gap between the cylinder and the wall, G, from G=0 (cylinder touching the wall) to G/D=2, was investigated. It is shown that the flow may be characterized by four distinct regions. (a) For very small gaps, G/D≤0·125, the gap flow is suppressed or extremely weak, and separation of the boundary layer occurs both upstream and downstream of the cylinder. Although there is no regular vortex shedding, there is a periodicity associated with the outer shear-layer. (b) In the “small gap ratio” region, 0·125<G/D<0·5, the flow is very similar to that for very small gaps, except that there is now a pronounced pairing between the inner shear-layer shed from the cylinder and the wall boundary layer. (c) Intermediate gap ratios, 0·5<G/D<0·75, are characterized by the onset of vortex shedding from the cylinder. (d) For the fourth region, characterized by the largest gap ratios considered, G/D>1·0, there is no separation of the wall boundary layer, either upstream or downstream of the cylinder.     

Adaptive wall functions for the v2‐f turbulence model

Adaptive wall functions for the v²‐f turbulence model have been derived for the flow over a flat plate at zero pressure gradient. These wall functions were implemented via tables for the turbulence quantities and the friction velocity u_τ. A special treatment for the ε and f boundary conditions is proposed. On fine grids (y⁺<1) this approach yields results consistent with the wall integration solution. Detailed numerical results are presented for a zero pressure gradient boundary layer and separated flow over a ramp. Copyright © 2006 John Wiley & Sons, Ltd.     

An evaluation of a 3D free‐energy‐based lattice Boltzmann model for multiphase flows with large density ratio

In this paper, the 3D Navier–Stokes (N–S) equation and Cahn–Hilliard (C–H) equations were solved using a free‐energy‐based lattice Boltzmann (LB) model. In this model, a LB equation with a D3Q19 velocity model is used to recover continuity and N–S equations while another LB equation with D3Q7 velocity model for solving C–H equation (Int. J. Numer. Meth. Fluids, 2008; 56 :1653–1671) is applied to solve the 3D C–H equation. To avoid the excessive use of computational resources, a moving reference frame is adopted to allow long‐time simulation of a bubble rising. How to handle the inlet/outlet and moving‐wall boundary conditions are suggested. These boundary conditions are simple and easy for implementation. This model's performance on two‐phase flows was investigated and the mass conservation of this model was evaluated. The model is validated by its application to simulate the 3D air bubble rising in viscous liquid (density ratio is 1000). Good agreement was obtained between the present numerical results and experimental results when Re is small. However, for high‐Re cases, the mass conservation seems not so good as the low‐Re case. Copyright © 2009 John Wiley & Sons, Ltd.     

Centreline velocity decay characterisation in low-velocity jets downstream from an extended conical diffuser

The centreline velocity decay of round airflow jets issuing from extended conical diffusers with length-to-diameter ratio 1.2≤L _t /d≤20 is studied for moderate bulk Reynolds numbers 1131≤Re _b ≤9054. The centreline velocity decay varies as a function of the initial conditions. The functional correlation between the centreline velocity decay coefficient and the initial centreline turbulence level observed on convergent nozzles (Malmström et al. in J. Fluid Mech. 246:363–377, 1997) breaks down as the initial centreline turbulence level exceeds 20 %. In addition, the centreline velocity decay coefficient expressed as function of the bulk velocity U _b decreases for U _b <3 m/s instead of initial mean velocity U ₀<6 m/s as reported for convergent nozzles (Malmström et al. in J. Fluid Mech. 246:363–377, 1997). The asymptotic values of the decay coefficient for U _b >3 m/s decrease linearly when expressed as function of the initial centreline turbulence intensity u ₀/U ₀. Studied flow and geometrical conditions are relevant to flow through the human upper airways.     

Quadratic programming algorithm for wall slip and free boundary pressure condition

Wall slip is often observed in a highly sheared fluid film in a solid gap. This makes a difficulty in mathematical analysis for the hydrodynamic effect because fluid velocity at the liquid–solid interfaces is not known a priori. If the gap has a convergent–divergent wedge, a free boundary pressure condition, i.e. Reynolds pressure boundary condition, is usually used in the outlet zone in numerical solution. This paper, based on finite element method and parametric quadratic programming technique, gives a numerical solution technique for a coupled boundary non‐linearity of wall slip and free boundary pressure condition. It is found that the numerical error decreases with the number of elements in a negative power law having an index larger than 2. Our method does not need an iterative process and can simultaneously gives rise to fluid film pressure distribution, wall slip velocity and surface shear stress. Wall slip always decreases the hydrodynamic pressure. Large wall slip even causes a null hydrodynamic pressure in a pure sliding solid gap. Copyright © 2005 John Wiley & Sons, Ltd.     

Vorticity and Regularity for Viscous Incompressible Flows under the Dirichlet Boundary Condition. Results and Related Open Problems

In reference [7] it is proved that the solution of the evolution Navier–Stokes equations in the whole of R ³ must be smooth if the direction of the vorticity is Lipschitz continuous with respect to the space variables. In reference [5] the authors improve the above result by showing that Lipschitz continuity may be replaced by 1/2-H?lder continuity. A central point in the proofs is to estimate the integral of the term (ω · ∇)u · ω, where u is the velocity and ω = ∇ × u is the vorticity. In reference [4] we extend the main estimates on the above integral term to solutions under the slip boundary condition in the half-space R ₊³. This allows an immediate extension to this problem of the 1/2-H?lder sufficient condition. The aim of these notes is to show that under the non-slip boundary condition the above integral term may be estimated as well in a similar, even simpler, way. Nevertheless, without further hypotheses, we are not able now to extend to the non slip (or adherence) boundary condition the 1/2-H?lder sufficient condition. This is not due to the “nonlinear" term (ω · ∇)u · ω but to a boundary integral which is due to the combination of viscosity and adherence to the boundary. On the other hand, by appealing to the properties of Green functions, we are able to consider here a regular, arbitrary open set Ω.      

Pressure boundary condition for the time‐dependent incompressible Navier–Stokes equations

In Gresho and Sani (Int. J. Numer. Methods Fluids 1987; 7 :1111–1145; Incompressible Flow and the Finite Element Method, vol. 2. Wiley: New York, 2000) was proposed an important hypothesis regarding the pressure Poisson equation (PPE) for incompressible flow: Stated there but not proven was a so‐called equivalence theorem (assertion) that stated/asserted that if the Navier–Stokes momentum equation is solved simultaneously with the PPE whose boundary condition (BC) is the Neumann condition obtained by applying the normal component of the momentum equation on the boundary on which the normal component of velocity is specified as a Dirichlet BC, the solution (u, p) would be exactly the same as if the ‘primitive’ equations, in which the PPE plus Neumann BC is replaced by the usual divergence‐free constraint (? · u = 0), were solved instead. This issue is explored in sufficient detail in this paper so as to actually prove the theorem for at least some situations. Additionally, like the original/primitive equations that require no BC for the pressure, the new results establish the same thing when the PPE approach is employed. Copyright © 2005 John Wiley & Sons, Ltd.     

The prediction of spatially periodic flows using a finite‐volume model

A periodic boundary condition has been developed that can be used in conjunction with a specified flow rate to produce accurate results in spatially periodic geometries. This condition is useful in situations where the flow rate is known, or more importantly, in cases where the pressure gradient is not known a priori, such as in countercurrent flows. Using the present condition, the flow rate is imposed at the inlet in terms of a bulk velocity, but the velocity field evolves as part of the solution. The condition is formulated to be suitable for both fixed and moving periodic domains. For the case of a moving domain, a correction is introduced to account for changes in the instantaneous velocity through the periodic edges. Under periodic conditions, these corrections integrate to zero over a complete (temporal) period. The new periodic condition is shown to produce accurate results for flat and wavy‐walled channels under both induced flow and countercurrent conditions. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical heat transfer in a cavity with a solar control coating deposited to a vertical semitransparent wall

A transient two‐dimensional computational model of combined natural convection, conduction, and radiation in a cavity with an aspect ratio of one, containing air as a laminar and non‐participating fluid, is presented. The cavity has two opaque adiabatic horizontal walls, one opaque isothermal vertical wall, and an opposite semitransparent wall, which consists of a 6‐mm glass sheet with a solar control coating of SnS–Cu_xS facing the cavity. The semitransparent wall also exchanges heat by convection and radiation from its external surface to the surroundings and allows solar radiation pass through into the interior of the cavity. The momentum and energy equations in the transient state were solved by finite differences using the alternating direction implicit (ADI) technique. The transient conduction equation and the radiative energy flux boundary conditions are coupled to these equations. The results in this paper are limited to the following conditions: 10⁴≤Gr≤10⁶, an isothermal vertical cold wall of 21°C, outside air temperatures in the range 30°C≤T₀≤40°C and incident solar radiation of AM2 (750 W m⁻²) normal to the semitransparent wall. The model allows calculation of the redistribution of the absorbed component of solar radiation to the inside and outside of the cavity. The influences of the time step and mesh size were considered. Using arguments of energy balance in the cavity, it was found that the percentage difference was less than 4 per cent, showing a possible total numerical error less than this number. For Gr=10⁶ a wave appeared in the upper side of the cavity, suggesting the influence of the boundary walls over the air flow inside the cavity. A Nusselt number correlation as a function of the Rayleigh number is presented. Copyright © 2000 John Wiley & Sons, Ltd.     

Free (open) boundary condition: some experiences with viscous flow simulations

The free (or open) boundary condition (FBC, OBC) was proposed by Papanastasiou et al. (A new outflow boundary condition, International Journal for Numerical Methods in Fluids, 1992; 14:587–608) to handle truncated domains with synthetic boundaries where the outflow conditions are unknown. In the present work, implementation of the FBC has been tested in several benchmark problems of viscous flow in fluid mechanics. The FEM is used to provide numerical results for both cases of planar and axisymmetric domains under laminar, isothermal or non‐isothermal, steady‐state conditions, for Newtonian fluids. The effects of inertia, gravity, compressibility, pressure dependence of the viscosity, slip at the wall, and surface tension are all considered individually in the extrudate‐swell benchmark problem for a wide range of the relevant parameters. The present results extend previous ones regarding the applicability of the FBC and show cases where the FBC is inappropriate, namely in the extrudate‐swell problem with gravity or surface‐tension effects. Particular emphasis has been given to the pressure at the outflow, which is the most sensitive quantity of the computations. In all cases where FBC is appropriate, excellent agreement has been found in comparisons with results from very long domains. The formulation for Picard‐type iterations is given in some detail, and the differences with the Newton–Raphson formulation are highlighted regarding some computational aspects. Copyright © 2011 John Wiley & Sons, Ltd.     

Consistent inlet and outlet boundary conditions for particle methods

In this paper, simple and consistent open boundary conditions are presented for the numerical simulation of viscous incompressible laminar flows. The present approach is based on an arbitrary Lagrangian-Eulerian particle method using upwind interpolation. Three kinds of inlet/outlet boundary conditions are proposed for particle methods, a pressure specified inlet/outlet condition, a velocity profile specified inlet/outlet condition, and a fully developed flow outlet condition. These inlet/outlet conditions are realized by using boundary particles and modification to the physical value such as velocity. Poiseuille flows, flows over a backward-facing step, and flows in a T-shape branch are calculated. The results are compared with those of mesh-based methods such as the finite volume method. The method presented herein exhibits accuracy and numerical stability.