Predicting turbulent flow in a staggered tube bundle

This paper presents the results of calculations performed for the turbulent, incompressible flow around a staggered array of tubes for which carefully obtained experimental results are available as part of an established ERCOFTAC-IAHR test case. The Reynolds-averaged Navier–Stokes equations are solved using a pressure-based finite volume algorithm, using collocated cell vertex store on an unstructured and adaptive mesh of tetrahedra. Turbulence closure is obtained with a truncated form of a low-Reynolds number k– model developed by Yang and Shih. The computational domain covers all seven rows of tubes used in the experimental study and periodic flow is allowed to develop naturally. The results of the computations are surprisingly good and compare favourably with results obtained by others using a wide range of alternative k– models for a single cylinder with periodic inflow and outflow boundaries on structured meshes.     

Computational and experimental study of annular photo-reactor hydrodynamics

The performance of ultraviolet (UV) reactors used for water treatment is greatly influenced by the reactor hydrodynamics due to the non-homogeneity of the radiation field. Reliable modeling of the reactor flow structures is therefore crucial for the design process. In this study, the turbulent flow through two characteristic annular UV-reactor configurations, with inlets concentric (L-shape) and normal (U-shape) to the reactor axis, was investigated through computational fluid dynamic (CFD). The modeling results were evaluated with the velocity profiles from particle image velocimetry (PIV) experiments. The influence of mesh structure and density, as well as three turbulence models: Standard κ–, Realizable κ–, and Reynolds stress model (RSM), on the simulation results were evaluated. Mesh-independent solutions were achieved at mean cell volumes of 5 × 10⁻⁹ m³. The Realizable κ– displayed the best overall match to the experimental PIV measurements. In general, the CFD models showed a close agreement with the experimental data for the majority of the reactor domain and captured the influences of reactor configuration and internal reactor structures on the flow distribution. The validated CFD hydrodynamic models could be integrated with kinetic and radiation distribution models for UV-reactor performance simulation.     

Experimental study of surface waves with Faraday resonance excitation

The results of an investigation of standing two-dimensional gravity waves on the free surface of a homogeneous liquid, induced by the vertical oscillations of a rectangular vessel under Faraday resonance conditions, are presented. The frequency ranges of excitation are determined and resonance relationships for the second and third modes are obtained and analyzed. Nonlinearities of the waves generated, such as wave profile asymmetry and node oscillations, are evaluated. Wave breakdown and the onset of unstable oscillation modes are considered. Experimental results are compared with the theoretical data.The experimental studies [1–4], devoted to Faraday resonance, deal mainly with the conditions under which resonance arises and the frequency response.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–129, January–February, 1995.     

A coupled-mode model for water wave scattering by vertically sheared currents in variable bathymetry regions

A coupled-mode model is developed for treating the wave–current–seabed interaction problem, with application to wave scattering by non-homogeneous, sheared current with linear vertical velocity profile, over general bottom topography. The wave potential is represented by a series of local vertical modes containing the propagating and evanescent modes, plus additional terms accounting for the satisfaction of the boundary conditions. Using the above representation, in conjunction with a variational principle, a coupled system of differential equations on the horizontal plane is derived, with respect to the unknown modal amplitudes. In the case of small-amplitude waves, a linearized version of the above coupled-mode system is obtained, extending previous analysis by Belibassakis et al. (2011) to the propagation of water waves over variable bathymetry regions in the presence of vertically sheared currents. Keeping only the propagating mode in the vertical expansion of the wave potential, the present system reduces to a one-equation model, that is shown to extend known mild-slope mild vertical shear equation concerning wave–current interaction over slowly varying topography. After additional simplifications, the latter model is shown to be compatible with the extended mild-slope mild-shear equation by Touboul et al. (2016). Results are presented for various representative test cases demonstrating the usefulness of the present coupled mode system and the importance of various terms in the modal expansion, and compared against experimental data collected in wave flume validating the present method. The analytical structure of the present system facilitates extensions to model non-linear effects and applications concerning wave scattering by inhomogeneous currents in coastal regions with general 3D bottom topography.     

Hybrid finite‐volume finite‐difference scheme for the solution of Boussinesq equations

A hybrid scheme composed of finite‐volume and finite‐difference methods is introduced for the solution of the Boussinesq equations. While the finite‐volume method with a Riemann solver is applied to the conservative part of the equations, the higher‐order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy in space for the finite‐volume solution is achieved using the MUSCL‐TVD scheme. Within this, four limiters have been tested, of which van‐Leer limiter is found to be the most suitable. The Adams–Basforth third‐order predictor and Adams–Moulton fourth‐order corrector methods are used to obtain fourth‐order accuracy in time. A recently introduced surface gradient technique is employed for the treatment of the bottom slope. A new model ‘HYWAVE’, based on this hybrid solution, has been applied to a number of wave propagation examples, most of which are taken from previous studies. Examples include sinusoidal waves and bi‐chromatic wave propagation in deep water, sinusoidal wave propagation in shallow water and sinusoidal wave propagation from deep to shallow water demonstrating the linear shoaling properties of the model. Finally, sinusoidal wave propagation over a bar is simulated. The results are in good agreement with the theoretical expectations and published experimental results. Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental investigation of cylindrically diverging solitons in an electric lattice

The results of an experimental investigation of cylindrical solitons in a two-dimensional electric LC-lattice are given. It is shown that in the continuum limit, propagation of cylindrical waves far from the center of symmetry in such a lattice may be described for each ray tube by a known modification of the Korteweg-de Vries equation which takes account of the cylindrical divergence. The dispersion term in this equation depends on the direction of wave propagation relative to the direction of the main axes of the lattice. Formation of solitons from non-soliton-shaped pulses was observed. The variations of soliton amplitude and duration with distance have been determined. They agree well with the numerical calculations by Maxon & Viecelli [2] and Dorfman [9]. Comparison of the obtained experimental data with the known theoretical laws of amplitude attenuation for diverging solitons [2, 12, 14] seems to favor the validity of the law A r^−2/3 rather than A r^−1/2.     

The Green’s function of the mild-slope equation:: The case of a monotonic bed profile

In the present work the Green’s function of the mild-slope and the modified mild-slope equations is studied. An effective numerical Fourier inversion scheme has been developed and applied to the construction and study of the source-generated water-wave potential over an uneven bottom profile with different depths at infinity. In this sense, the present work is a prerequisite to the study of the diffraction of water waves by localized bed irregularities superimposed over an uneven bottom. In the case of a monotonic bed profile, the main characteristics of the far-field are: (i) the formation of a shadow zone with an ever expanding width, which is located along the bottom irregularity, and (ii) in each of the two sectors not including the bottom irregularity the asymptotic behavior of the wave field approaches the form of an outgoing cylindrical wave, propagating with an amplitude of order O(R^−1/2), where R is the distance from the source, and wavelength corresponding to the sector-depth at infinity. Moreover, the weak wave system propagating in the shadow zone is of order O(R^−3/2), and along the bottom irregularity consists of the superposition of two outgoing waves with wavelengths corresponding to the two depths at infinity.     

Wavelet Analysis of the Evolution of a Solitary Wave in a Composite Material

The evolution of a solitary wave propagating through a microstructural material (composite) is studied on the basis of wavelet analysis. A specific feature of the solution technique proposed is the use of Mexican hat (MH) wavelets, which are elastic wavelets, i.e., they are solutions of the basic system of wave equations for an elastic material with a microstructure. The initial wave profile is also chosen in the form of the MH-wavelet. Primary attention is given to the relationship among the profile behavior, wave bottom length, and characteristic microstructure length. A computer analysis conducted demonstrates that the approach proposed allows us to detect the basic wave effects: splitting of the wave into two modes with different phase velocities, simultaneous propagation of both modes in the components of the composite, and strong dependence of the evolution rate on the characteristic lengths of the wave and microstructure     

On one-dimensional shock waves

In this paper uniform asymptotic expansions for the solutions of a system of differential equations are obtained in the domain containing a shock wave. It is shown, in particular, that the function θ(t,x)/ε contained in the expansions and describing the behavior of the solution in the neighborhood of the wave front has, generally speaking, a discontinuity of derivatives at the front. The results are applicable to one-dimensional problems in gas dynamics with low viscosity and heat-conductivity.     

考虑桩土相互作用效应的桩顶纵向振动时域响应分析

在考虑桩土耦合作用以及土竖向波动效应条件下，对均质滞回材料阻尼土中弹性支承桩桩顶纵向振动时域响应进行了理论研究。首先建立了桩与滞回阻尼土在谐和振动情况下的定解问题，然后先对土层动力平衡方程进行求解并得到土体振动位移形式解，接着依据平衡条件将该形式解祸合进桩身动力平衡方程，并通过对桩动力平衡方程的求解，最终得到桩顶位移和速度频域响应解析解和半正弦脉冲激励作用下桩顶时域响应的半解析解。通过与其它相关理论解的对比验证了本文解的正确性和适用性，并基于所得解对桩顶时域响应特性进行了分析，最后将理论曲线与现场工程实测曲线进行了拟合对比，结果表明两者符合较好。     

A second order stationary solution of the three dimensional wave diffraction potential

A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solution of an initial value problem, and using Laplace transformation method and Tauberian theorem, the boundary value problems of stationary solution of the first and second order diffraction potential have been derived in this paper. Furthermore, the explicit solution of the second order stationary diffraction potential has been obtained with the method of the double Fourier transformation. It is found that the asymptotic behaviour of the second order stationary solution at far field is dependent on two wave systems, the first is “free wave”, travelling independently of the first order wave system, the other is “phase locked waves”, which accompany the first order waves. At the same time, the radiation conditions of the second order diffraction problems are derived. We also find that one can still pursue a steady state formulation with the inclusion of Rayleigh damping. Finally, as an example, the second order wave forces upon a fixed vertical circular cylinder have been calculated, and the numerical results agree well with the experimental data.     

Uni-directional waves over slowly varying bottom. Part I: Derivation of a KdV-type of equation

The exact equations for surface waves over an uneven bottom can be formulated as a Hamiltonian system, with the total energy of the fluid as Hamiltonian. If the bottom variations are over a length scale L that is longer than the characteristic wavelength ℓ, approximating the kinetic energy for the case of “rather long and rather low” waves gives Boussinesq type of equations. If in the case of an even bottom one restricts further to uni-directional waves, the Korteweg-de Vries (KdV) is obtained. For slowly varying bottom this uni-directionalization will be studied in detail in this part I, in a very direct way which is simpler than other derivations found in the literature. The surface elevation is shown to be described by a forced KdV-type of equation. The modification of the obtained KdV-equation shares the property of the standard KdV-equation that it has a Hamiltonian structure, but now the structure map depends explicitly on the spatial variable through the bottom topography. The forcing is derived explicitly, and the order of the forcing, compared to the first order contributions of dispersion and nonlinearity in KdV, is shown to depend on the ratio between ℓ and L; for very mild bottom variations, the forcing is negligible. For localized topography the effect of this forcing is investigated. In part II the distortion of solitary waves will be studied.     

Exact analysis of Burgers equation on semiline with flux condition at the origin

The initial boundary value problem for the Burgers equation in the domain x 0, t > 0 with flux boundary condition at x = 0 has been solved exactly. The behaviour of the solution as t tends to infinity is studied and the “asymptotic profile at infinity” is obtained. In addition, the uniqueness of the solution of the initial boundary value problem is proved and its inviscid limit as → 0 is obtained.     

Propagation and localization of Rossby waves over random topography (two-layer model)

We study propagation of Rossby waves over randomly stratified bottom topography in a two-layer system. The problem is reduced to coupled stochastic wave equations in latitudinal variable with suitable boundary conditions. When the two-dimensional system is broken into its basic barotropic and baroclinic modes, each one is scattered by the medium, and generates the like-wise and the cross-wise components. We study statistics of the reflection and transmission coefficients for both types of generated modes, and show localization for the like-wise components, and propagation for the cross-wise components. The localization lengths and the transmission rates are estimated in terms of the basic parameters of the system and the correlation-function of topographic fluctuations.     

SINGLE-LAYER FINITE ELEMENT APPROXIMATION OF 3D NAVIER–STOKES EQUATIONS WITH FREE BOUNDARIES FOR SHALLOW WATER FLOWS

In order to simulate flows in the shallow water limit, the full incompressible Navier–Stokes equations with free boundaries are solved using a single layer of finite elements. This implies a polynomial approximation of the velocity profile in the vertical direction, which in turn distorts the wave speed. This fact is verified by numerical results: the wave speed depends on the vertical discretization. When at least two layers of finite elements are used, the boundary layer at the bottom can be simulated and the correct solution for the shallow water limit is recovered. Then this algorithm is applied to the prediction of Tsunami event.     

Reflection of a beam of elastic waves by a periodic surface profile

Reflection of elastic waves from a traction-free solid-air boundary of periodic saw-tooth profile is investigated analytically and experimentally. For an incident plane wave the surface displacements on the profile are computed as the solution of a singular integral equation. The reflected field is subsequently obtained by using an integral representation. Incident beams of finite width are represented by Fourier superpositions of plane waves. The dependence of the reflected signal spectra on the incident beam width is examined closely near the fundamental surface resonance frequency. Experimental spectra which were obtained using two different diameter transducers, are compared to the corresponding theoretical spectra. It is found that the depth of the spectral minima depends on the incident beam width. Both analytical and experimental results exhibit the splitting of an incident beam of elastic waves into two reflected beams. The beam splitting is more pronounced for a narrower incident beam and for frequencies close to a resonance frequency of the profile.     

Aspects of self-similar diffraction gas flows and their numerical simulation

A large number of papers has been devoted to the investigation of the interaction of a plane shock wave with bodies of various geometric shapes, and they have been generalized and classified for a stationary body in [1, 2]. Separate results of experimental and theoretical investigations of the interaction of a shock wave with a wedge, cone, sphere, and cylinder moving with supersonic velocities are contained in [3–9]. Analysis of the available results shows that the features of the unsteady gas flows formed in this case largely depend on the nature of the boundary-value problem that arises for the system of differential gas dynamic equations. The question of the wave structure of the unsteady gas flow and the accuracy of the obtained solution is central to the numerical investigation of the present class of problems. The most characteristic types of unsteady self-similar gas flows that arise on the interaction of a plane shock wave with bodies of a wedge or convex corner type are calculated on the basis of an explicit numerical continuous calculation method of the second order of accuracy. The accuracy of the numerical solutions is discussed on the basis of a comparison with the experimental data. The case of the interaction of a shock wave with the rarefaction wave that arises in a supersonic flow past a convex corner is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 146–152, July–August, 1986.     

Gas-dynamic approach in the nonlinear theory of ion acoustic waves in a plasma: An exact solution

An exact solution of the problem of the acoustic wave structure in a plasma is obtained. Both plasma component are treated as gases with specified initial temperatures and adiabatic exponents. The system of equations describing the wave profile is solved using an original method consisting of reducing the system to the Bernoulli equation. A numerical example of the obtained general solution of the problem of the wave profile for arbitrary parameters is given. Curves are constructed that bound the region of existence of a stationary solitary ion acoustic wave in the parameter space. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 5, pp. 3–11, September–October, 2007.     

数值模拟在钱塘江涌潮分析中的应用──Ⅱ.计算结果和分析

关于钱塘江涌潮分析中用数值模拟的计算结果和分析,对于一维情况,采用有精确解问题的计算、水跃的数值模拟和溃坝的数值模拟来检验数值方法的正确性,并将该方法应用于钱塘江涌潮的计算．对于二维情况,对存在理论解的直道水激波与有实验数据的湾道溃坝和非线性水波爬坡问题进行了计算,并将该方法应用于钱塘江涌潮的计算,数值计算结果表明了方法的有效性．     

Wave flow regimes of a thin layer of viscous fluid subject to gravity

The studies of Kapitsa initiated the detailed experimental and theoretical study of the flow of a thin layer of viscous liquid (liquid film) over a solid surface [1–2]. Extensive experimental data on this question have now been accumulated. As a rule, the existing theories are based on linearization of the problem and diverge considerably from the experimental results. The present paper is also addressed to the theoretical solution of this problem. The solution method used enables consideration of the wave flow of the liquid as a nonlinear problem and on this basis permits determining all the parameters of the wave regime-amplitude, wavelength, wave propagation speed, frequency.