Numerical simulations of inviscid and viscous free‐surface flows with application to nonlinear water waves

The purpose of the present study is to establish a numerical model appropriate for solving inviscid/viscous free‐surface flows related to nonlinear water wave propagation. The viscous model presented herein is based on the Navier–Stokes equations, and the free‐surface is calculated through an arbitrary Lagrangian–Eulerian streamfunction‐vorticity formulation. The streamfunction field is governed by the Poisson equation, and the vorticity is obtained on the basis of the vorticity transport equation. For computing the inviscid flow the Laplace streamfunction equation is used. These equations together with the respective (appropriate) fully nonlinear free‐surface boundary conditions are solved using a finite difference method. To demonstrate the model feasibility, in the present study we first simulate collision processes of two solitary waves of different amplitudes, and compute the phenomenon of overtaking of such solitary waves. The developed model is subsequently applied to calculate (both inviscid and the viscous) flow field, as induced by passing of a solitary wave over submerged rectangular structures and rigid ripple beds. Our study provides a reasonably good understanding of the behavior of (inviscid/viscous) free‐surface flows, within the framework of streamfunction‐vorticity formulation. The successful simulation of the above‐mentioned test cases seems to suggest that the arbitrary Lagrangian–Eulerian/streamfunction‐vorticity formulation is a potentially powerful approach, capable of effectively solving the fully nonlinear inviscid/viscous free‐surface flow interactions. Copyright © 2009 John Wiley & Sons, Ltd.     

Towards a transparent boundary condition for compressible Navier–Stokes equations

A new artificial boundary condition for two‐dimensional subsonic flows governed by the compressible Navier–Stokes equations is derived. It is based on the hyperbolic part of the equations, according to the way of propagation of the characteristic waves. A reference flow, as well as a convection velocity, is used to properly discretize the terms corresponding to the entering waves. Numerical tests on various classical model problems, whose solutions are known, and comparisons with other boundary conditions (BCs), show the efficiency of the BC. Direct numerical simulations of more complex flows over a dihedral plate are simulated, without creation of acoustic waves going back in the flow. Copyright © 2001 John Wiley & Sons, Ltd.     

A study of the effect of step excrescences and free‐stream disturbances on boundary layer stability

In this work, a study of the mechanism by which free‐stream acoustic and vorticity disturbances interact with a boundary layer flow developing over a flat plate featuring a step excrescence located at a certain distance from a blunt leading edge is included. The numerical tool is a high‐fidelity implicit numerical algorithm solving for the unsteady, compressible form of the Navier–Stokes equations in a body‐fitted curvilinear coordinates and employing high‐accurate compact differencing schemes with Pade‐type filters. Acoustic and vorticity waves are generated using a source term in the momentum and energy equations, as opposed to using inflow boundary conditions, to avoid spurious waves that may propagate from boundaries. The results show that the receptivity to surface step excrescences is largely the result of an overall adverse pressure gradient posed by the step, and that the free‐stream disturbances accelerate the generation of instabilities in the downstream. As expected, it is found that the acoustic disturbance interacting with the surface imperfection is more efficient in exciting the Tollmien–Schlichting waves than the vorticity disturbance. The latter generates Tollmien–Schlichting waves that are grouped in wave packets consistent with the wavelength of the free‐stream disturbance. Copyright © 2015 John Wiley & Sons, Ltd.     

Non-homogeneous radiation and outflow boundary conditions simulating incoming acoustic and vorticity waves for exterior computational aeroacoustics problems

A set of non-homogeneous radiation and outflow boundary conditions which automatically generate prescribed incoming acoustic or vorticity waves and, at the same time, are almost transparent to outgoing sound waves produced internally in a finite computation domain is proposed. This type of boundary condition is needed for the numerical solution of many exterior aeroacoustics problems. In computational aeroacoustics, the computation scheme must be as non-dispersive and non-dissipative as possible. It must also support waves with wave speeds which are nearly the same as those of the original linearized Euler equations. To meet these requirements, a high-order/large-stencil scheme is often necessary. The proposed non-homogeneous radiation and outflow boundary conditions are designed primarily for use in conjunction with such high-order/large-stencil finite difference schemes. © 1998 John Wiley & Sons, Ltd.     

Resonance in rarefied gases

Dispersion and damping of ultrasound waves are a standard test for mathematical models of rarefied gas flows. Normally, one considers waves in semi-infinite systems in relatively large distance of the source. For a more complete picture, ultrasound propagation in finite closed systems of length L is studied by means of several models for rarefied gas flows: the Navier-Stokes-Fourier equations, Grad’s 13 moment equations, the regularized 13 moment equations, and the Burnett equations. All systems of equations are considered in simple 1-D geometry with their appropriate jump and slip boundary conditions. Damping and resonance are studied in dependence of frequency and length. For small L, all wave modes contribute to the solution.     

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.     

Comparison of outflow boundary conditions for subsonic aeroacoustic simulations

Aeroacoustics simulations require much more precise boundary conditions than classical aerodynamics. Two classes of non‐reflecting boundary conditions for aeroacoustics are compared in the present work: the characteristic analysis‐based methods and the Tam and Dong approach. In the characteristic methods, waves are identified and manipulated at the boundaries, whereas the Tam and Dong approach use modified linearized Euler equations in a buffer zone near outlets to mimic a non‐reflecting boundary. The principles of both approaches are recalled, and recent characteristic methods incorporating the treatment of transverse terms are discussed. Three characteristic techniques—the original Navier–Stokes characteristic boundary conditions (NSCBC) of Poinsot and Lele and two versions of the modified method of Yoo and Im—are compared with the Tam and Dong method for four typical aeroacoustics problems: vortex convection on a uniform flow, vortex convection on a shear flow, acoustic propagation from a monopole, and acoustic propagation from a dipole. Results demonstrate that the Tam and Dong method generally provides the best results and is a serious alternative solution to characteristic methods even though its implementation might require more care than the usual NSCBC approaches. Copyright © 2011 John Wiley & Sons, Ltd.     

A hybrid approach for nonlinear computational aeroacoustics predictions

In many aeroacoustics applications involving nonlinear waves and obstructions in the far-field, approaches based on the classical acoustic analogy theory or the linearised Euler equations are unable to fully characterise the acoustic field. Therefore, computational aeroacoustics hybrid methods that incorporate nonlinear wave propagation have to be constructed. In this study, a hybrid approach coupling Navier–Stokes equations in the acoustic source region with nonlinear Euler equations in the acoustic propagation region is introduced and tested. The full Navier–Stokes equations are solved in the source region to identify the acoustic sources. The flow variables of interest are then transferred from the source region to the acoustic propagation region, where the full nonlinear Euler equations with source terms are solved. The transition between the two regions is made through a buffer zone where the flow variables are penalised via a source term added to the Euler equations. Tests were conducted on simple acoustic and vorticity disturbances, two-dimensional jets (Mach 0.9 and 2), and a three-dimensional jet (Mach 1.5), impinging on a wall. The method is proven to be effective and accurate in predicting sound pressure levels associated with the propagation of linear and nonlinear waves in the near- and far-field regions.     

Multidimensional upwind schemes for the shallow water equations

A multidimensional discretisation of the shallow water equations governing unsteady free-surface flow is proposed. The method, based on a residual distribution discretisation, relies on a characteristic eigenvector decomposition of each cell residual, and the use of appropriate distribution schemes. For uncoupled equations, multidimensional convection schemes on compact stencils are used, while for coupled equations, either system distribution schemes such as the Lax–Wendroff scheme or scalar schemes may be used. For steady subcritical flows, the equations can be partially diagonalised into a purely convective equation of hyperbolic nature, and a set of coupled equations of elliptic nature. The multidimensional discretisation, which is second-order-accurate at steady state, is shown to be superior to the standard Lax–Wendroff discretisation. For steady supercritical flows, the equations can be fully diagonalised into a set of convective equations corresponding to the steady state characteristics. Discontinuities such as hydraulic jumps, are captured in a sharp and non-oscillatory way. For unsteady flows, the characteristic equations remain coupled. An appropriate treatment of the coupling terms allows the discretisation of these equations at the scalar level. Although presently only first-order-accurate in space and time, the classical dam-break problem demonstrates the validity of the approach. © 1998 John Wiley & Sons, Ltd.     

Pseudo‐characteristic formulation and dynamic boundary conditions for computational aeroacoustics

The present paper deals with the use of the pseudo‐characteristic formulation of the Navier–Stokes and Euler equations recently introduced by Sesterhenn (Comput. Fluid. 2001; 30 :37–67) for the simulation of acoustic wave propagation. The emphasis is put on the formulation of an efficient method on structured curvilinear grids, along with the definition and implementation of efficient boundary conditions. The cases of inflow, outflow, rigid/compliant walls and walls with prescribed impedance are addressed. The proposed boundary conditions are assessed on generic cases. The pseudo‐characteristic formulation enables a straightforward and optimal use of high‐order upwind dispersion‐relation‐preserving schemes, yielding an efficient method. Copyright © 2006 John Wiley & Sons, Ltd.     

Unsteady interaction of uniformly vortex flows

The problem of the decay of an arbitrary discontinuity (the Riemann problem) for the system of equations describing vortex plane-parallel flows of an ideal incompressible liquid with a free boundary is studied in a long-wave approximation. A class of particular solutions that correspond to flows with piecewise-constant vorticity is considered. Under certain restrictions on the initial data of the problem, it is proved that this class contains self-similar solutions that describe the propagation of strong and weak discontinuities and the simple waves resulting from the nonlinear interaction of the specified vortex flows. An algorithm for determining the type of resulting wave configurations from initial data is proposed. It extends the known approaches of the theory of one-dimensional gas flows to the case of substantially two-dimensional flows. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 5, pp. 55–66, September–October, 1998.     

Sound and vorticity interactions: transmission and scattering

We review several aspects of the propagation of sound in vortical flows. We restrict ourselves to isothermal, humidity-free flows at low Mach number M. Since vorticity plays a major role in vortex-flow interactions we focus on vortical flows. We consider two main canonical situations. The first concerns the transmission of sound. We analyze the evolution of acoustic wavefronts as they propagate across a single vortex. The second situation addresses the scattering of sound waves by nonstationary vortices. We study the evolution of the acoustic pressure emitted in the far field, at an angle with the initial direction of propagation. In this geometry one performs direct spectroscopy of the flow vorticity field. In each case, we review theoretical results and compare with experimental measurements and numerical simulations when available. We also briefly report how the following new acoustic techniques have recently been used to study complex or turbulent flows: time-resolved acoustic spectroscopy, speckle interferometry and Lagrangian particle tracking. PACS 43.25, 43.28, 47.32, 67.40     

位势流动和非均匀介质中声波方程的36种形式

声波方程是对大多数声学问题进行数学描述的出发点. 那些得到 广泛应用的经典波动方程及对流波动方程都存在苛刻的适用条件, 即仅适用于描述处于静态或匀速运动状态的定常 均匀介质中的线性无耗散声波. 然而, 很多实际场合并不满足这些严格的适用条件. 本文对经典声波方程和对流声波 方程进行推广, 导出了编号为W1$\sim$W36的36种不同形式的声波方程, 涵盖了处于静止、势流或旋涡流状态下的非均匀 和/或非定常介质中的声波传播问题. 所考虑的声波传播情形包括: (1) 线性波, 即具有小梯度(小振幅)性质; (2)非线性波, 即具有陡峭梯度性质, 包括``波纹'(小振幅大梯度)或者大振幅波. 本文仅考虑非耗散声波, 即排除了由剪切、体积黏度及热传导所引起的耗散. 对具有匀熵或等熵(熵沿流线守恒)性质的均匀介质和非均匀介质中的声传播进行了研究但非等熵(即耗散)情况除外; 另外, 对非定常介质中的 声波问题也进行了分析. 所涉及的介质可以处于静止、匀速运动状态, 或者是非匀速的和/或非定常的平均流动, 包括: (1)低Mach数的势平均流(即不可压缩的平均态), 或高速势平均流(即非均匀可压缩的平均流); ② 变截面管 道中的准一维传播, 包括无平均流的号管和具有低或高Mach数平均流的喷管; 或③平面的、空间的、或轴对称的单 向剪切平均流. 本文没有探讨其他类型的旋涡平均流(将与耗散及其他情形一起留待下一步研究), 例如, 可能与剪切效应相结合的轴对称旋转平均流. 通过对流体力学的一般方程进行消元处理或根据声学变分原理, 导出了36种波动方程, 对一些波动方程还采用这两种方法进行相互校验. 尽管声波方程的36种形式没有涵盖非线性、非均匀与非定常及非匀速运动介质 这3个效应的所有可能的组合情形, 但它们的确包括了孤立状态下的各种效应, 并包括了多种多重效应组合的 情形. 虽然经典波动方程和对流波动方程仅适用于处于静止(或匀速运动)的均匀定常介质中的线性无耗散声波, 但它们在 相关文献中已被广泛采用; 本文给出的36种声波方程提供了它们多种有用的推广形式. 在许多实际应用中, 经典波动方 程和对流波动方程仅是粗略的近似, 声波方程的更一般形式可提供更令人满意的理论模型. 本文每节末尾给出了这些应用 的众多范例. 在这篇评论文章中引用了240篇参考文献.     

The effect of virtual mass on the basic equations for unsteady one-dimensional heterogeneous flows

A separated two-component flow model is presented which includes virtual mass forces coupling the momentum equations of the two components. It is shown that for physically realistic situations four real roots of the characteristics determinant can exist. These are associated with the acoustic propagation velocities and the fiow velocities of the constitutive phases. Direct analytical solution of the full characteristic determinant is difficult. However, for low Mach number flows an acoustic propagation velocity is obtained which falls between the well-known true separated and homogeneous wave speeds, and compares favorably with experimental data for glass/water and air/water mixtures.     

Bleustein–Gulyaev waves in 6 mm piezoelectric materials loaded with a viscous liquid layer of finite thickness

In this paper, we analyze the propagation of Bleustein–Gulyaev waves in an unbounded piezoelectric half-space loaded with a viscous liquid layer of finite thickness within the linear elastic theories. Exact solutions of the phase velocity equations are obtained in the cases of both electrically open circuit and short circuit by solving the equilibrium equations of piezoelectric materials and the diffusion equation of viscous liquid. A PZT-5H/Glycerin system is selected to perform the numerical calculation. The results show that the mass density and the viscous coefficient have different effects on the propagation attenuation and phase velocity under different electrical boundary conditions. In particular, the penetration depth of the waves is of the same order as the wavelength in the case of electrically short circuit. These effects can be used to manipulate the behavior of the waves and have implications in the application of acoustic wave devices.     

Toward a reduction of mesh imprinting

This paper is the initial investigation into a new Lagrangian cell‐centered hydrodynamic scheme that is motivated by the desire for an algorithm that resists mesh imprinting and has reduced complexity. Key attributes of the new approach include multidimensional construction, the use of flux‐corrected transport (FCT) to achieve second order accuracy, automatic determination of the mesh motion through vertex fluxes, and vorticity control. Toward this end, vorticity preserving Lax–Wendroff (VPLW) type schemes for the two‐dimensional acoustic equations were analyzed and then implemented in a FCT algorithm. Here, mesh imprinting takes the form of anisotropic dispersion relationships. If the stencil for the LW methods is limited to nine points, four free parameters exist. Two parameters were fixed by insisting that no spurious vorticity be created. Dispersion analysis was used to understand how the remaining two parameters could be chosen to increase isotropy. This led to new VPLW schemes that suffer less mesh imprinting than the rotated Richtmyer method. A multidimensional, vorticity preserving FCT implementation was then sought using the most promising VPLW scheme to address the problem of spurious extrema. A well‐behaved first order scheme and a new flux limiter were devised in the process. The flux limiter is unique in that it acts on temporal changes and does not place a priori bounds on the solution. Numerical results have demonstrated that the vorticity preserving FCT scheme has comparable performance to an unsplit MUSCL‐H algorithm at high Courant numbers but with reduced mesh imprinting and superior symmetry preservation. Copyright © 2014 John Wiley & Sons, Ltd.     

Simulation of non‐hydrostatic,density‐stratified flow in irregular domains

A numerical model has been developed for simulating density‐stratified flow in domains with irregular but simple topography. The model was designed for simulating strong interactions between internal gravity waves and topography, e.g. exchange flows in contracting channels, tidally or convectively driven flow over two‐dimensional sills or waves propagating onto a shoaling bed. The model is based on the non‐hydrostatic, Boussinesq equations of motion for a continuously stratified fluid in a rotating frame, subject to user‐configurable boundary conditions. An orthogonal boundary fitting co‐ordinate system is used for the numerical computations, which rely on a fourth‐order compact differentiation scheme, a third‐order explicit time stepping and a multi‐grid based pressure projection algorithm. The numerical techniques are described and a suite of validation studies are presented. The validation studies include a pointwise comparison of numerical simulations with both analytical solutions and laboratory measurements of non‐linear solitary wave propagation. Simulation results for flows lacking analytical or laboratory data are analysed a posteriori to demonstrate satisfaction of the potential energy balance. Computational results are compared with two‐layer hydraulic predictions in the case of exchange flow through a contracting channel. Finally, a simulation of circulation driven by spatially non‐uniform surface buoyancy flux in an irregular basin is discussed. Copyright © 2000 John Wiley & Sons, Ltd.     

Over-reflection and instabilities in shear flows of shallow water with bottom friction

In a two-dimensional shear flow of shallow water, the bottom friction relates uniquely the spanwise profile of the depth-averaged velocity to the bottom topography. If the basic flow varies weakly in the spanwise direction, the local analysis of stability at every spanwise position gives the region of the flow parameters for which the classic hydraulic instability due to the bottom friction cannot occur. In this region, the linear analyses of the waves scattering and instability due to the lateral shear can be performed effectively by means of the frictionless linearized equations if both the bottom slope and friction are equally small.The energy of the total perturbed flow can be split into three main parts that correspond to the basic flow, small amplitude wave motion and induced mean flow. The waves can be either amplified or damped near the critical layers, where their streamwise phase velocity equals the velocity of the basic flow. Two physical mechanisms of this amplification exist. The first one is similar to that suggested by Takehiro and Hayashi for a linear frictionless shallow water flow. The incident and transmitted waves carry energy of opposite signs, which results in an increase in the amplitude of the reflected wave compared to that of the incident one. This mechanism of over-reflection operates for any combination of the flow parameters. The other mechanism is similar to Landau damping in plasma flows; it is related to the energy exchange between the waves and fluid particles at the critical layers due to the velocity synchronism. It may lead to either additional amplification or damping of the waves for different flow conditions. In particular, its significance can be reduced by stronger bottom friction. If the basic flow has uniform potential vorticity, Landau damping is negligible, and over-reflection always occurs. If the feed-back is provided by another critical layer, the net over-reflection results in the formation of trapped modes.     

Transonic viscous-inviscid interaction by a finite element method

A method is outlined for solving two-dimensional transonic viscous flow problems, in which the velocity vector is split into the gradient of a potential and a rotational component. The approach takes advantage of the fact that for high-Reynolds-number flows the viscous terms of the Navier-Stokes equations are important only in a thin shear layer and therefore solution of the full equations may not be needed everywhere. Most of the flow can be considered inviscid and, neglecting the entropy and vorticity effects, a potential model is a good approximation in the flow core. The rotational part of the flow can then be calculated by solution of the potential, streamfunction and vorticity transport equations. Implementation of the no-slip and no-penetration boundary conditions at the walls provides a simple mechanism for the interaction between the viscous and inviscid solutions and no extra coupling procedures are needed. Results are presented for turbulent transonic internal choked flows.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.