A FINITE VOLUME SOLVER FOR THE SHALLOW WATER EQUATIONS IN VARYING BED ENVIRONMENTS

Abstract

A numerical scheme for solving the shallow-water equations is presented. An analogy is made between flows governed by shallow-water equations and the Euler system of equations used in gas dynamics. An emphasis is placed on the difference presented by the bathymetry in hydraulic systems. The discretization of the governing equations is based on Roe's flux difference-splitting solver, initially developed for solving inviscid compressible flows. The spatial discretization is handled within a finite-volume context by using triangles or quadrilaterals as the basic control-volume cells. This approach enables an easy and flexible treatment of general geometries. A development of the boundary conditions tailored for the current scheme is given. Fundamental validation tests are presented.     

Hybrid mesh generation using advancing reduction technique

This paper presents an application of the advancing reduction technique for 2D hybrid mesh generation (triangles + quadrilaterals). Based on an initial rectangle mesh (RM) covering the whole domain, the advancing reduction technique coarsens the base RM in a marching way from the boundary to the interior of the domain so that different zones of sub‐RMs with different edge lengths are recognized. These sub‐RMs are connected to each other with the so‐called transition layers which consist of the transition triangles and quadrilaterals. As demonstrated by examples, the proposed method is simple, efficient, and easy to implement. Copyright © 2015 John Wiley & Sons, Ltd.     

On the Green's Functions and Boundary Integral Formulation of Elastic Planes with Cutouts

ABSTRACT

The problem of an infinite elastic plane that contains a hole of arbitrary shape and is subjected to a concentrated unit load is considered. The Green's function (influence function) for the problem is formulated by means of two complex potential functions. This is accomplished by mapping the region that is exterior to the hole onto a unit circle. A class of closed contour hole shapes is analyzed. Green's functions for an elliptical hole and a class of triangular holes are determined. Green's functions for a class of rectangular holes are also discussed. In order to determine stress and displacement fields for the finite plane problem, Green's function is employed and an indirect boundary integral equation is formulated, with the integrand of the integral equation incorporating the effect of the hole. The contour of the hole is no longer considered a part of the boundary and only the contour of the region that is exterior to the hole is subdivided into boundary elements. Examples for elliptical and triangular holes are solved.     

Complex variable function method fob hole shape optimization in an elastic plane

In this paper, a complex variable function method for solving the hole shape optimization problem in an elastic plane is presented. In this method, the stresses in hole problems are analysed by taking advantage of the efficiency of the complex variable function method. To optimize the hole shape, the coeffecients in conformal mapping functions are taken as design variables, and the sensitivity analysis and gradient methods are used to reduce the largest circumferential stress in absolute value and at the same time to make the second largest circumferential stress in absolute value not to exceed the largest one (in fact, these two stresses are the stationary values of the circumferential stresses). The coefficients in conformal mapping function are revised by iteration step by step until the largest circumferential stress in absolute value is reduced to the second largest stress. This method guarantees the continuity, differentiability and accuracy of the stress solution along the boundary, and it is evident that this method is better than either the difference method or the finite element method.     

The degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary

This paper investigates the degenerate scale problem for the Laplace equation and plane elasticity in a multiply connected region with an outer circular boundary. Inside the boundary, there are many voids with arbitrary configurations. The problem is analyzed with a relevant homogenous BIE (boundary integral equation). It is assumed that all the inner void boundary tractions are equal to zero, and tractions on the outer circular boundary are constant. Therefore, all the integrations in BIE are performed on the outer circular boundary only. By using the relation z * conjg(z) = a * a, or conjg(z) = a * a/z on the circular boundary with radius a, all integrals can be reduced to an integral for complex variable and they can be integrated in closed form. The degenerate scale a = 1 is found in the Laplace equation and in plane elasticity regardless of the void configuration.     

The Linearized Crocco Equation

In this paper, we study the existence and uniqueness of a degenerate parabolic equation, with nonhomogeneous boundary conditions, coming from the linearization of the Crocco equation [12]. The Crocco equation is a nonlinear degenerate parabolic equation obtained from the Prandtl equations with the so-called Crocco transformation. The linearized Crocco equation plays a major role in stabilization problems of fluid flows described by the Prandtl equations [5]. To study the infinitesimal generator associated with the adjoint linearized Crocco equation – with homogeneous boundary conditions – we first study degenerate parabolic equations in which the x-variable plays the role of a time variable. This equation is doubly degenerate: the coefficient in front of ∂_x vanishes on a part of the boundary, and the coefficient of the elliptic operator vanishes in another part of the boundary. This makes very delicate the proof of uniqueness of solution. To overcome this difficulty, a uniqueness result is first obtained for an equation in which the elliptic operator is symmetric, and it is next extended to the original equation by combining an iterative process and a fixed point argument (see Th. 4.9). This kind of argument is also used to prove estimates, which cannot be obtained in a classical way.     

Well-Posedness in Smooth Function Spaces for the Moving-Boundary Three-Dimensional Compressible Euler Equations in Physical Vacuum

We prove well-posedness for the three-dimensional compressible Euler equations with moving physical vacuum boundary, with an equation of state given by p(ρ) =  C _γ ρ ^γ for γ > 1. The physical vacuum singularity requires the sound speed c to go to zero as the square-root of the distance to the moving boundary, and thus creates a degenerate and characteristic hyperbolic free-boundary system wherein the density vanishes on the free-boundary, the uniform Kreiss–Lopatinskii condition is violated, and manifest derivative loss ensues. Nevertheless, we are able to establish the existence of unique solutions to this system on a short time-interval, which are smooth (in Sobolev spaces) all the way to the moving boundary, and our estimates have no derivative loss with respect to initial data. Our proof is founded on an approximation of the Euler equations by a degenerate parabolic regularization obtained from a specific choice of a degenerate artificial viscosity term, chosen to preserve as much of the geometric structure of the Euler equations as possible. We first construct solutions to this degenerate parabolic regularization using a higher-order version of Hardy’s inequality; we then establish estimates for solutions to this degenerate parabolic system which are independent of the artificial viscosity parameter. Solutions to the compressible Euler equations are found in the limit as the artificial viscosity tends to zero. Our regular solutions can be viewed as degenerate viscosity solutions. Our methodology can be applied to many other systems of degenerate and characteristic hyperbolic systems of conservation laws.     

Analytical and numerical investigation for true and spurious eigensolutions of an elliptical membrane using the real-part dual BIEM/BEM

The paper performs analytical and numerical investigation of the true and spurious eigensolutions of an elliptical membrane using the real-part boundary integral equation method (BIEM) following the successful work on a circular case by using the dual boundary element method (BEM) (Kuo et al. in Int. J. Numer. Methods Eng. 48:1401–1422, 2000). We extend to the elliptical case in this paper. To analytically study the eigenproblems of an elliptical membrane, the elliptical coordinates and Mathieu functions are adopted. The fundamental solution is expanded into the degenerate kernel by using the elliptical coordinates and the boundary densities are expanded by using the eigenfunction expansion. The Jacobian terms may exist in the degenerate kernel, boundary density and boundary contour integration but they can cancel each other out. Therefore, the orthogonal relations are reserved in the boundary contour integral. It is interesting to find that the BIEM using the real or the imaginary-part kernel to deal with an elliptical membrane yields spurious eigensolutions. This finding agrees with those corresponding to the circular case. The spurious eigenvalues in the real-part BIEM are found to be the zeros of the mth-order (even or odd) modified Mathieu functions of the second kind or their derivatives. To verify this finding, the BEM is implemented. Furthermore, the commercial finite-element code ABAQUS is also utilized to provide eigensolutions for comparisons. It is found that good agreement is obtained.     

Electrophoresis of a 2-particle cluster near a plane boundaryElectrophorèse de deux particules en présence d'une paroi plane

Particle–boundary and particle–particle interactions in Electrophoresis are examined by considering a 2-particle cluster near a plane boundary. The advocated treatment holds for two insulating particles of arbitrary shapes and zeta potential functions and resorts to 13 boundary-integral equations. Preliminary results reveal that, depending upon the addressed velocity nature (translational or angular), wall–particle may be stronger or weaker than particle–particle interactions. To cite this article: A. Sellier, C. R. Mecanique 331 (2003).     

A finite volume method for fluid flow in polar cylindrical grids

In the numerical simulation of fluid flows using a polar cylindrical grid, grid lines meet at a single point on the axis of the polar cylindrical grid system; this makes the grids around the axis degenerate from being general quadrilaterals into triangles. Therefore, a special treatment must be performed when the axis has to be included in the computational domain in order to solve a non-axisymmetrical fluid flow problem. In this paper a new numerical method has been developed to deal with the difficulty of the axis when the control volume technique is used with a non-staggered grid arrangement. Two illustrative examples of the proposed method are presented for simulating the fluid flows on the axis and all the numerical results obtained for the two examples are shown to be in good agreement with the available analytical solutions. © 1998 John Wiley & Sons, Ltd.     

Existence of a Global Smooth Solution¶for a Degenerate Goursat Problem¶of Gas Dynamics

The Goursat problem of a mixed type equation , P≥ 0, is considered. At the ends of its supports we have P=0, which means it is degenerate hyperbolic. We prove the global existence of a smooth solution to the degenerate Goursat problem up to a boundary where P=0. This problem comes from the expansion of a wedge of gas with constant velocity into vacuum, in two-dimensional pressure-gradient equations in gas dynamics, where P is the pressure and P=0 means vacuum. Accepted June 16, 2000?Published online December 6, 2000     

Asymptotic Behavior of Global Classical Solutions to the Mixed Initial-Boundary Value Problem for Quasilinear Hyperbolic Systems with Small BV Data

In this paper, we investigate the asymptotic behavior of global classical solutions to the mixed initial-boundary value problem with small BV data for linearly degenerate quasilinear hyperbolic systems with general nonlinear boundary conditions in the half space {(t,x)|t≥0,x≥0}. Based on the existence result on the global classical solution, we prove that when t tends to the infinity, the solution approaches a combination of C ¹ traveling wave solutions, provided that the C ¹ norm of the initial and boundary data is bounded and the BV norm of the initial and boundary data is sufficiently small. Applications to quasilinear hyperbolic systems arising in physics and mechanics, particularly to the system describing the motion of the relativistic string in the Minkowski space-time R ¹⁺ⁿ , are also given.     

A boundary element formulation using the simple method

A new boundary element method is described for calculation of the steady incompressible laminar flows. The method is based on the well-known SIMPLE algorithm. The new boundary element method allows one to find the fields of the pressure and velocity corrections without inner iterations, thus reducing the computational time drastically. This makes it different from the method developed by Patankar and Spalding.³² However, the new method demands a much larger computer strorage. The boundary integral equations are discretized with the help of constant boundary elements and constant cells. The values of the integrals along the boundary elements and the cells for the two-dimensional domain are found analytically. To preserve the stability in the iteration process, under-relaxation for the convection terms is used. This paper gives the results of calculations of the flows between two plane parallel plates at Re = 20 and Re = 200, the flows in a square cavity with a moving upper lid at Re = 1 and Re = 100 and the flow in a plane channel with sudden symmetric expansion at Re =46·6.     

The radar cross section of the semi-infinite elliptic cone

In order to determine the radar cross section of a semi-infinite elliptic cone first the solution of the pertinent canonical boundary value problem in sphero-conal coordinates is derived. For this purpose one has to solve a two-parameter eigenvalue problem with two coupled Lamé differential equations. The exact nose-on radar cross section of the semi-infinite elliptic cone has been evaluated for two orthogonal polarizations of the incident plane wave. In the degenerate case of a circular cone the known formula first deduced by Hansen and Schiff is obtained. The theoretical results also hold for the case of a plane angular sector which is another degenerate case of the elliptic cone.     

Mixed Boundary Value Problem in a Non-homogeneous Medium Under Steady Distribution of Temperature

This paper is concerned with a mixed boundary value problem of a non-homogeneous medium under steady distribution of temperature whose elastic constants are exponential functions of y. By using Fourier cosine transforms the mixed boundary value problem of heat conduction is reduced to a Fredholm integral equation of the second kind. Then the elastic problem of the non-homogeneous semi-infinite half-plane under distribution of load over a plane face is solved. The influence of the non-homogeneity of the material on the thermal stress distribution is illustrated graphically.     

Transient Growth in Exactly Counter-Rotating Couette–Taylor Flow

Transient growth due to non-normality is investigated for the Couette-Taylor problem with counter-rotating cylinders as a function of aspect ratio η and Reynolds number Re. For all Re≤500, transient growth is enhanced by curvature, i.e. is greater for η<1 than for η=1, the plane Couette limit. For fixed Re>130, it is found that the greatest transient growth is achieved for η on the linear stability boundary. Transient growth is approximately 20% higher near the Couette-Taylor linear stability boundary at Re=310, η=0.986 than at Re=310, η=1, near the threshold observed for transition in plane Couette flow. For 106<Re<130, the greatest transient growth occurs for a value of η between the linear stability boundary and one. For Re<106, the flow is linearly stable and the greatest transient growth occurs for a value of η less than one. The energy in the optimal inputs is primarily meridional; that in the optimal outputs is primarily azimuthal. Pseudospectra are calculated for two contrasting cases. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 Communicated by H.J.S. Fernando     

Energy release or absorption due to simultaneous expansion of many interacting nanoholes in elastic materials

The energy release or absorption due to simultaneous expansion of many interacting nanoholes in elastic materials under plane strain deformation is studied as influenced by the surface effect along rims of nanoholes. The M-integral classically used in macro mechanics with defects is extended to treat the problem with many interacting nanoholes. After some manipulations, the energy change due to the simultaneous expansion of many nanoholes represented by the M-integral is evaluated. Four different arrays of many nanoholes under a monotonically increasing tensile loading are considered. Attention is focused on the influence induced from the surface tension, the surface Lamé constants, and the interaction among many nanoholes on the M-integral. It is concluded that the surface tension yields significant influence on the M-integral, whereas the surface Lamé constants offer much smaller influence, which could be neglected with some relative errors less than 2%. It is found that, unlike those in macro mechanics with defects, the simultaneous expansion could either release energy (the positive value of the M-integral) or absorb energy (the negative value of the M-integral), depending on the loading levels. There is a neutral loading point, at which the M-integral transforms from a negative value to a positive value in all arrays of nanovoids under consideration. It is also found that the interaction among multiple nanoholes influences the value of the neutral loading point significantly because the mutual influence induced from both the interacting effect and the surface effect yields a quite different feature from those induced from the interacting effect only. That is, the surface effect always inhibits the influence of the interacting effect on the M-integral.     

Using a piecewise linear bottom to fit the bed variation in a laterally averaged,z‐co‐ordinate hydrodynamic model

In developing a 3D or laterally averaged 2D model for free‐surface flows using the finite difference method, the water depth is generally discretized either with the z‐co‐ordinate (z‐levels) or a transformed co‐ordinate (e.g. the so‐called σ‐co‐ordinate or σ‐levels). In a z‐level model, the water depth is discretized without any transformation, while in a σ‐level model, the water depth is discretized after a so‐called σ‐transformation that converts the water column to a unit, so that the free surface will be 0 (or 1) and the bottom will be ‐1 (or 0) in the stretched co‐ordinate system. Both discretization methods have their own advantages and drawbacks. It is generally not conclusive that one discretization method always works better than the other. The biggest problem for the z‐level model normally stems from the fact that it cannot fit the topography properly, while a σ‐level model does not have this kind of a topography‐fitting problem. To solve the topography‐fitting problem in a laterally averaged, 2D model using z‐levels, a piecewise linear bottom is proposed in this paper. Since the resulting computational cells are not necessarily rectangular looking at the x–z plane, flux‐based finite difference equations are used in the model to solve the governing equations. In addition to the piecewise linear bottom, the model can also be run with full cells or partial cells (both full cell and partial cell options yield a staircase bottom that does not fit the real bottom topography). Two frictionless wave cases were chosen to evaluate the responses of the model to different treatments of the topography. One wave case is a boundary value problem, while the other is an initial value problem. To verify that the piecewise linear bottom does not cause increased diffusions for areas with steep bottom slopes, a barotropic case in a symmetric triangular basin was tested. The model was also applied to a real estuary using various topography treatments. The model results demonstrate that fitting the topography is important for the initial value problem. For the boundary value problem, topography‐fitting may not be very critical if the vertical spacing is appropriate. Copyright © 2004 John Wiley & Sons, Ltd.     

Exact calculation of the flow of a staggered tandem cascade with moving second blade row

Summary  The plane flow around a tandem cascade of flat plates is calculated by means of conformal mapping. The blades of the two rows are perpendicular to each other. The first row is stationary, the second row moves with constant velocity. The conformal mapping will be constructed by a “mapping flow”. The blades of one row are stream lines and those of the other row are potential lines of the flow. By conformal mapping, the physical flow around the tandem cascade of the physical ζ-plane is converted into a flow between infinitely long straight walls in the z-plane, each wall corresponding to one of the blades. The conditions far upstream and far downstream of the cascade are represented by source-vortices. In the z-plane, the boundary conditions may be easily fulfilled by reflection and repetition of the source-vortices, and the flow may be calculated by well-known methods. The physical flow searched for is then obtained by inverse mapping. Received 24 July 2000; accepted for publication 6 December 2000     

Numerical prediction of turbulent boundary layer noise from a sharp-edged flat plate

An efficient hybrid uncorrelated wall plane waves–boundary element method (UWPW-BEM) technique is proposed to predict the flow-induced noise from a structure in low Mach number turbulent flow. Reynolds-averaged Navier-Stokes equations are used to estimate the turbulent boundary layer parameters such as convective velocity, boundary layer thickness, and wall shear stress over the surface of the structure. The spectrum of the wall pressure fluctuations is evaluated from the turbulent boundary layer parameters and by using semi-empirical models from literature. The wall pressure field underneath the turbulent boundary layer is synthesized by realizations of uncorrelated wall plane waves (UWPW). An acoustic BEM solver is then employed to compute the acoustic pressure scattered by the structure from the synthesized wall pressure field. Finally, the acoustic response of the structure in turbulent flow is obtained as an ensemble average of the acoustic pressures due to all realizations of uncorrelated plane waves. To demonstrate the hybrid UWPW-BEM approach, the self-noise generated by a flat plate in turbulent flow with Reynolds number based on chord Re_c = 4.9 × 10⁵ is predicted. The results are compared with those obtained from a large eddy simulation (LES)-BEM technique as well as with experimental data from literature.