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.     

Finite analytic boundary condition for a higher‐order finite analytic scheme

A higher‐order finite analytic scheme based on one‐dimensional finite analytic solutions is used to discretize three‐dimensional equations governing turbulent incompressible free surface flow. In order to preserve the accuracy of the numerical scheme, a new, finite analytic boundary condition is proposed for an accurate numerical solution of the partial differential equation. This condition has higher‐order accuracy. Thus, the same order of accuracy is used for the boundary. Boundary conditions were formulated and derived for fluid inflow, outflow, impermeable surfaces and symmetry planes. The derived boundary conditions are treated implicitly and updated with the solution of the problem. The basic idea for the derivation of boundary conditions was to use the discretized form of the governing equations for the fluid flow simplified on the boundaries and flow information. To illustrate the influence of the higher‐order effects at the boundaries, another, lower‐order finite analytic boundary condition, is suggested. The simulations are performed to demonstrate the validity of the present scheme and boundary conditions for a Wigley hull advancing in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

AUGMENTED LAGRANGIAN METHOD AND OPEN BOUNDARY CONDITIONS IN 2D SIMULATION OF POISEUILLE–BÉNARD CHANNEL FLOW

The main objective of this study is to compare the influence of different boundary conditions upon the incompressible Poiseuille –Bénard channel flow (PBCF) in a 2D rectangular duct heated from below. In a first technical part the algorithm used to carry out this work, based on the augmented Lagrangian method, is presented. The implementation details of the five different open boundary conditions (OBCs) and the periodic boundary conditions (PBCs) tested in the present paper are also given. The study is then carried out for 1800<Ra≤ 10,000, 0<Re≤10 and 0·67≤Pr≤6·4. The five selected OBCs, applied at the outlet of the computational domain, respectively express the following conditions: a square profile for the velocity (OBC1), mass conservation (OBC2), zero second derivative of the horizontal velocity component (OBC3), a mixed boundary condition combining Dirichlet and Neumann conditions (OBC4) and an Orlanski-type boundary condition (OBC5). A good estimation of the perturbation amplitude and of the length of the perturbed zone at the outlet boundary is proposed. It is shown that OBC5 causes very little perturbation in the recirculating flow compared with the other OBCs. © 1997 John Wiley & Sons, Ltd.     

Boundary element analysis of driven cavity flow for low and moderate Reynolds numbers

A boundary element method for steady two‐dimensional low‐to‐moderate‐Reynolds number flows of incompressible fluids, using primitive variables, is presented. The velocity gradients in the Navier–Stokes equations are evaluated using the alternatives of upwind and central finite difference approximations, and derivatives of finite element shape functions. A direct iterative scheme is used to cope with the non‐linear character of the integral equations. In order to achieve convergence, an underrelaxation technique is employed at relatively high Reynolds numbers. Driven cavity flow in a square domain is considered to validate the proposed method by comparison with other published data. Copyright © 2001 John Wiley & Sons, Ltd.     

Coupled Navier–Stokes/molecular dynamics simulations in nonperiodic domains based on particle forcing

This paper focuses on coupling methods for hybrid Navier–Stokes/molecular dynamics (MD) simulations. The computational domain is split in a continuum flow region, where a finite‐volume discretisation of the Navier–Stokes equations is used, and one or more particle domains, where molecular level modelling of the flow is employed. The domains are defined with a partial overlap, in which the flow states are coupled through an exchange of the velocity components. For the steady flows considered, an under‐relaxed Newton iteration method is used to drive the coupled system to convergence. The main focus of the present work is on methods to impose nonperiodic boundary conditions on the particle domain(s). A particle forcing is applied in the direction normal to the particle domain boundary to impose the boundary normal velocity component. A novel aspect of the present work is the extension of this method to more general nonplanar particle domain boundaries. The main contribution of the paper is the development of a particle forcing method in the direction tangential to the domain boundary, which is based on the equivalent continuum‐flow boundary shear stresses along with an iterative forcing strength adjustment based on the extrapolated particle boundary velocity. Furthermore, an adaptation scheme is presented, which uses the finite‐volume flux residuals of the particle bin averaged velocity field as a truncation criterion for the iterative force‐update scheme. It is demonstrated that by comparing the residual reduction for the momentum equation in the nonhomogeneous directions during the molecular dynamics simulations with that for a homogeneous direction, the forcing iteration at which the statistical noise in the velocity field dominates the uncertainty in the forcing strength can be determined. At this point the iteration can be truncated. It is shown that with adaptive schemes of this type, the total number of MD evaluations required in a coupled Navier–Stokes/MD simulation can be reduced relative to a hybrid scheme with a fixed number of forcing‐strength updates. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite volume solution of the Navier–Stokes equations in velocity–vorticity formulation

For the incompressible Navier–Stokes equations, vorticity‐based formulations have many attractive features over primitive‐variable velocity–pressure formulations. However, some features interfere with the use of the numerical methods based on the vorticity formulations, one of them being the lack of a boundary conditions on vorticity. In this paper, a novel approach is presented to solve the velocity–vorticity integro‐differential formulations. The general numerical method is based on standard finite volume scheme. The velocities needed at the vertexes of each control volume are calculated by a so‐called generalized Biot–Savart formula combined with a fast summation algorithm, which makes the velocity boundary conditions implicitly satisfied by maintaining the kinematic compatibility of the velocity and vorticity fields. The well‐known fractional step approaches are used to solve the vorticity transport equation. The paper describes in detail how we accurately impose no normal‐flow and no tangential‐flow boundary conditions. We impose a no‐flux boundary condition on solid objects by the introduction of a proper amount of vorticity at wall. The diffusion term in the transport equation is treated implicitly using a conservative finite update. The diffusive fluxes of vorticity into flow domain from solid boundaries are determined by an iterative process in order to satisfy the no tangential‐flow boundary condition. As application examples, the impulsively started flows through a flat plate and a circular cylinder are computed using the method. The present results are compared with the analytical solution and other numerical results and show good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Global volume conservation in unsteady free surface flows with energy absorbing far‐end boundaries

A wave absorption filter for the far‐end boundary of semi‐infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier–Stokes equations and depth‐integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far‐end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic energy of the outgoing waves. The computational scheme is tested to verify the conservation of total fluid volume in the domain for long simulation durations. Combination of the Sommerfeld boundary and dissipation zone can effectively minimize reflections and prevent cumulative changes in total fluid volume in the domain. Solitary wave, nonlinear periodic waves and irregular waves are simulated to illustrate the numerical developments. Earthquake excited surface waves and nonlinear hydrodynamic pressures in a dam–reservoir are computed. Copyright © 2009 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).     

Accuracy and stability of a set of free‐surface time‐domain boundary element models based on B‐splines

An analysis is given for the accuracy and stability of some perturbation‐based time‐domain boundary element models (BEMs) with B‐spline basis functions, solving hydrodynamic free‐surface problems, including forward speed effects. The spatial convergence rate is found as a function of the order of the B‐spline basis. It is shown that for all the models examined the mixed implicit–explicit Euler time integration scheme is correct to second order. Stability diagrams are found for models based on B‐splines of orders third through to sixth for two different time integration schemes. The stability analysis can be regarded as an extension of the analysis by Vada and Nakos [Vada T, Nakos DE. Time marching schemes for ship motion simulations. In Proceedings of the 8th International Workshop on Water Waves and Floating Bodies, St. John's, Newfoundland, Canada, 1993; 155–158] to include B‐splines of orders higher than three (piecewise quadratic polynomials) and to include finite water depth and a current at an oblique angle to the model grid. Copyright © 2000 John Wiley & Sons, Ltd.     

Extension of domain‐free discretization method to simulate compressible flows over fixed and moving bodies

This paper is the first endeavour to present the local domain‐free discretization (DFD) method for the solution of compressible Navier–Stokes/Euler equations in conservative form. The discretization strategy of DFD is that for any complex geometry, there is no need to introduce coordinate transformation and the discrete form of governing equations at an interior point may involve some points outside the solution domain. The functional values at the exterior dependent points are updated at each time step to impose the wall boundary condition by the approximate form of solution near the boundary. Some points inside the solution domain are constructed for the approximate form of solution, and the flow variables at constructed points are evaluated by the linear interpolation on triangles. The numerical schemes used in DFD are the finite element Galerkin method for spatial discretization and the dual‐time scheme for temporal discretization. Some numerical results of compressible flows over fixed and moving bodies are presented to validate the local DFD method. Copyright © 2006 John Wiley & Sons, Ltd.     

Artificial boundary method for the exterior Stokes flow in three dimensions

In this paper, the artificial boundary method is considered for the numerical simulation of the exterior Stokes flow in three dimensions. First, an exact relation between the normal stress and the velocity field is obtained on a spherical artificial boundary. With the relation specified on the artificial boundary, the original problem is reduced to a new one only defined on a finite domain. After that, an variational problem equivalent to the reduced problem is derived. By truncating the series term in the formulation, a sequence of approximate variational problems are obtained, which can then be solved with a suitable finite‐element scheme. Finally, a numerical example is presented to show the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

Towards simulation of flapping wings using immersed boundary method

In this work the immersed boundary method is applied to simulate incompressible turbulent flows around stationary and moving objects. The goal is to demonstrate that the immersed boundary technique along with a large eddy simulation approach is capable of simulating the effect of the so‐called leading edge vortex (LEV), which can be found in flapping wing aerodynamics. A Lagrangian method is used to approximate the solutions in the freshly cleared cells that lay within solid objects at one time step and emerge into fluid domain at the next time step. Flow around a stationary cylinder at Re_D = 20, 40, and 3900 (based on cylinder diameter D) is first studied to validate the immersed boundary solver based on the finite volume scheme using a staggered grid. Then, a harmonically oscillating cylinder at Re_D = 10 000 is considered to test the solver after the Lagrangian method is implemented to interpolate the solution in the freshly cleared cells. Finally, this approach is used to study flows around a stationary flat‐plate at several angles of attack and fast pitching flat‐plate. The rapidly pitching plate creates a dynamic LEV that can be used to improve the efficiency of flapping wings of micro air vehicle and to determine the optimum flapping frequency. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical study of heat island flows: Stationary solutions

We present two‐dimensional numerical simulations of a natural convection problem in an unbounded domain. The flow circulation is induced by a heat island located on the ground and thermal stratification is applied in the vertical direction. The main effect of this stably stratified environment is to induce the propagation of thermal perturbations in the horizontal direction far from the local thermal source. Numerical stationary solutions at Ra?10⁵ are computed in large elongated computational domains: convergence with respect to the domain sizes is investigated at different resolutions. On fine grids, with mesh size , a thermal sponge layer is added at the vertical boundaries: this local damping technique improves the convergence with respect to the domain length. Boussinesq equations are discretized with a second‐order finite volume scheme on a staggered grid combined with a second‐order projection method for the time integration. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Numerical simulation of fully nonlinear irregular wave tank in three dimension

A fully nonlinear irregular wave tank has been developed using a three‐dimensional higher‐order boundary element method (HOBEM) in the time domain. The Laplace equation is solved at each time step by an integral equation method. Based on image theory, a new Green function is applied in the whole fluid domain so that only the incident surface and free surface are discretized for the integral equation. The fully nonlinear free surface boundary conditions are integrated with time to update the wave profile and boundary values on it by a semi‐mixed Eulerian–Lagrangian time marching scheme. The incident waves are generated by feeding analytic forms on the input boundary and a ramp function is introduced at the start of simulation to avoid the initial transient disturbance. The outgoing waves are sufficiently dissipated by using a spatially varying artificial damping on the free surface before they reach the downstream boundary. Numerous numerical simulations of linear and nonlinear waves are performed and the simulated results are compared with the theoretical input waves. Copyright © 2006 John Wiley & Sons, Ltd.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.     

A boundary element method for the solution of finite mobility ratio immiscible displacement in a Hele‐Shaw cell

In this paper, the interaction between two immiscible fluids with a finite mobility ratio is investigated numerically within a Hele‐Shaw cell. Fingering instabilities initiated at the interface between a low‐viscosity fluid and a high‐viscosity fluid are analysed at varying capillary numbers and mobility ratios using a finite mobility ratio model. The present work is motivated by the possible development of interfacial instabilities that can occur in porous media during the process of CO₂ sequestration but does not pretend to analyse this complex problem. Instead, we present a detailed study of the analogous problem occurring in a Hele‐Shaw cell, giving indications of possible plume patterns that can develop during the CO₂ injection. The numerical scheme utilises a boundary element method in which the normal velocity at the interface of the two fluids is directly computed through the evaluation of a hypersingular integral. The boundary integral equation is solved using a Neumann convergent series with cubic B‐Spline boundary discretisation, exhibiting sixth‐order spatial convergence. The convergent series allows the long‐term nonlinear dynamics of growing viscous fingers to be explored accurately and efficiently. Simulations in low‐mobility ratio regimes reveal large differences in fingering patterns compared with those predicted by previous high‐mobility ratio models. Most significantly, classical finger shielding between competing fingers is inhibited. Secondary fingers can possess significant velocity, allowing greater interaction with primary fingers compared with high‐mobility ratio flows. Eventually, this interaction can lead to base thinning and the breaking of fingers into separate bubbles. Copyright © 2015 John Wiley & Sons, Ltd.     

Neural networks for BEM analysis of steady viscous flows

This paper presents a new neural network‐boundary integral approach for analysis of steady viscous fluid flows. Indirect radial basis function networks (IRBFNs) which perform better than element‐based methods for function interpolation, are introduced into the BEM scheme to represent the variations of velocity and traction along the boundary from the nodal values. In order to assess the effect of IRBFNs, the other features used in the present work remain the same as those used in the standard BEM. For example, Picard‐type scheme is utilized in the iterative procedure to deal with the non‐linear convective terms while the calculation of volume integrals and velocity gradients are based on the linear finite element‐based method. The proposed IRBFN‐BEM is verified on the driven cavity viscous flow problem and can achieve a moderate Reynolds number of 1400 using a relatively coarse uniform mesh. The results obtained such as the velocity profiles along the horizontal and vertical centrelines as well as the properties of the primary vortex are in very good agreement with the benchmark solution. Furthermore, the secondary vortices are also captured by the present method. Thus, it appears that an ability to represent the boundary solution accurately can significantly improve the overall solution accuracy of the BEM. Copyright © 2003 John Wiley & Sons, Ltd.