Two‐dimensional dispersion analyses of finite element approximations to the shallow water equations

Dispersion analysis of discrete solutions to the shallow water equations has been extensively used as a tool to define the relationships between frequency and wave number and to determine if an algorithm leads to a dual wave number response and near 2Δx oscillations. In this paper, we explore the application of two‐dimensional dispersion analysis to cluster based and Galerkin finite element‐based discretizations of the primitive shallow water equations and the generalized wave continuity equation (GWCE) reformulation of the harmonic shallow water equations on a number of grid configurations. It is demonstrated that for various algorithms and grid configurations, contradictions exist between the results of one‐dimensional and two‐dimensional dispersion analysis as a result of subtle changes in the mass matrix. Numerical experiments indicate that the two‐dimensional dispersion analysis correctly predicts the existence and onset of near 2Δx noise in the solution. Copyright © 2004 John Wiley & Sons, Ltd.     

The criterion of the existence or inexistence of transverse shock wave at wedge supported oblique detonation wave

A simplified theoretic method and numerical simulations were carried out to investigate the characterization of propagation of transverse shock wave at wedge supported oblique detonation wave.After solution validation,a criterion which is associated with the ratio Φ (u 2 /u CJ) of existence or inexistence of the transverse shock wave at the region of the primary triple was deduced systematically by 38 cases.It is observed that for abrupt oblique shock wave (OSW)/oblique detonation wave (ODW) transition,a transverse shock wave is generated at the region of the primary triple when Φ < 1,however,such a transverse shock wave does not take place for the smooth OSW/ODW transition when Φ > 1.The parameter Φ can be expressed as the Mach number behind the ODW front for stable CJ detonation.When 0.9 < Φ < 1.0,the reflected shock wave can pass across the contact discontinuity and interact with transverse waves which are originating from the ODW front.When 0.8 < Φ < 0.9,the reflected shock wave can not pass across the contact discontinuity and only reflects at the contact discontinuity.The condition (0.8 < Φ < 0.9) agrees well with the ratio (D ave /D CJ) in the critical detonation.     

A fibre sensor based frequency analysis of surface waves at hollow cone nozzles

The characterization of the disintegration process caused by sheet break up, e.g. at hollow cone nozzles (HCN) is of fundamental importance for a better understanding of the drop formation process. This paper introduces our new measuring device capable of resolving surface wave formation and propagation on the sheet along main streamlines of the flow. Furthermore, it reveals the development of adapted algorithms for the data and image evaluation. Two different frequency ranges are found for the oscillations on the sheet at higher viscosities, showing a superposition of different spectra in each case. The wave bands in the lower frequency range 10<f _low<80 s⁻¹ are considered to be irrelevant for the disintegration process. Surface waves with higher energy content at given conditions are found within a range 300<f _high<1,500 s⁻¹ probably leading to sheet disintegration. Results obtained from a photo-based evaluation show good agreement compared to the frequency analysis.     

The effects of boundary to the ECRH and the energy absorption of ICRH

Electron-cyclotron resonant heating (ECRH) of Tokamak plasma is examined. When plasma is heated by waves, we must consider the distribution of incident wave energy toO andX modes as the wave is incident from vacuum to the surface of plasma as well as the absorption efficiency ofO mode andX mode. Numerical calculation shows that for small incident angle, the incident energy transfers principally intoO mode when the electric fieldE ⁱ of incident wave is parallel to the incident plane, therefore it is efficient to heat the plasma byO mode. WhenE ⁱ is perpendicular to the incident plane, the energy transfers principally intoX mode and heating the plasma byX mode is efficient. Ion-cyclotron resonant heating (ICRH) is also considered, the formula of the energy of ion-cyclotron wave absorbed by plasma is obtained.     

A compact numerical algorithm for solving the time‐dependent mild slope equation

The mild slope equation has been widely used to describe combined wave refraction and diffraction. In this study, a new numerical algorithm is developed to solve the time‐dependent mild slope equation in a second‐order hyperbolic form. The numerical algorithm is based on a compact and explicit finite difference method that is second‐order accurate in both time and space. The algorithm has the similar structure to the leap‐frog method but is constructed on three time levels for the second‐order time derivative term. The numerical model has the capability of simulating transient wave motion by correctly predicting the speed of wave energy propagation, which is important for the real‐time forecast of the arrival time of storm waves generated in the far field. The model is validated against analytical solution for wave shoaling and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope (Coastal Eng. 1982; 6 :255). Lastly, the realistic scale Homma's island (Geophys. Mag. 1950; 21 :199) is studied with the use of various wave periods of T = 720s, T = 120 s, and T = 24 s. These wave periods correspond to long, intermediate, and short waves for the given topography, respectively. Comparisons are made between numerical results and existing analytical solutions in terms of the wave amplification around the island, which serves as the indicator for the potential wave runup. Excellent agreements are obtained. The model runs on a PC (Pentium IV 1.8GHz) and the computer capacity allows the computation of a mesh system up to 3000 × 3000, which is equivalent to about 150 × 150 waves or a large area of 540km × 540km for a wave train with the period of T = 60 s. Copyright 2004 John Wiley & Sons, Ltd.     

An experimental investigation on the dynamic axial buckling of cylindrical shells using a Kolsky Bar

Several experiments were performed with a Kolsky Bar (Split Hopkinson Pressure Bar) device to investigate the dynamic axial buckling of cylindrical shells. The Kolsky Bar is a loading as well as a measuring device which can subject the shells to a fairly good square pulse. An attempt is made to understand the interaction between the stress wave and the dynamic buckling of cylindrical shells. It is suggested that the dynamic axial buckling of the shells, elastic or elasto-plastic, is mainly due to the compressive wave rather than the flexural or bending wave. The experimental results seem to support the two critical velocity theory for plastic buckling, withV _c1 corresponding to an axisymmetric buckling mode andV _c2 corresponding to a non-symmetric buckling mode. The project supported by National Natural Science Foundation of China     

Genetic algorithm optimization of phononic bandgap structures

This paper describes the use of genetic algorithms (GAs) for the optimal design of phononic bandgaps in periodic elastic two-phase media. In particular, we link a GA with a computational finite element method for solving the acoustic wave equation, and find optimal designs for both metal–matrix composite systems consisting of Ti/SiC, and H₂O-filled porous ceramic media, by maximizing the relative acoustic bandgap for these media. The term acoustic here implies that, for simplicity, only dilatational wave propagation is considered, although this is not an essential limitation of the method. The inclusion material is found to have a lower longitudinal modulus (and lower wave speed) than the surrounding matrix material, a result consistent with observations that stronger scattering is observed if the inclusion material has a lower wave velocity than the matrix material.     

Two new predictor-corrector algorithms for second-order cone programming

Based on the ideas of infeasible interior-point methods and predictor-corrector algorithms, two interior-point predictor-corrector algorithms for the second-order cone programming (SOCP) are presented. The two algorithms use the Newton direction and the Euler direction as the predictor directions, respectively. The corrector directions belong to the category of the Alizadeh-Haeberly-Overton (AHO) directions. These algorithms are suitable to the cases of feasible and infeasible interior iterative points. A simpler neighborhood of the central path for the SOCP is proposed, which is the pivotal difference from other interior-point predictor-corrector algorithms. Under some assumptions, the algorithms possess the global, linear, and quadratic convergence. The complexity bound O(rln(ɛ ₀/ɛ)) is obtained, where r denotes the number of the second-order cones in the SOCP problem. The numerical results show that the proposed algorithms are effective.     

An experimental investigation on the airside performance of fin-and-tube heat exchangers having sinusoidal wave fins

The heat transfer and friction characteristics of the heat exchangers having sinusoidal wave fins were experimentally investigated. Twenty-nine samples having different waffle heights (1.5 and 2.0 mm), fin pitches (1.3–1.7 mm) and tube rows (1–3) were tested. Focus was given to the effect of waffle configuration (herringbone or sinusoidal) on the heat transfer and friction characteristics. Results show that the sinusoidal wave geometry provides higher heat transfer coefficients and friction factors than the herringbone wave geometry, and the difference increases as the number of row increases. The j/f ratios of the herringbone wave geometry, however, are larger than those of the sinusoidal wave geometry. Compared with the herringbone wave geometry, the sinusoidal wave geometry yielded a weak row effect, which suggests a superior heat transfer performance at the fully developed flow region for the sinusoidal wave geometry. Possible reasoning is provided considering the flow characteristics in wavy channels. Within the present geometric variations, the effect of waffle height on the heat transfer coefficient was not prominent. The effect of fin pitch was also negligible. Existing correlations highly overpredicted both the heat transfer coefficients and friction factors. A new correlation was developed based on the present data.     

Effect of acceleration on the wavy Taylor vortex flow

In this paper, we use a laser optical technique to investigate the characteristics of a wavy Taylor vortex flow between two concentric cylinders, with the inner cylinder subjected to a wide range of predetermined acceleration and the outer one at rest. We focus on the inner/outer radius ratio of 0.894, with an acceleration (dRe/dt*) from 0.1123 to 2,247, and Reynolds number from Re/Re _c =1.0 to 36. The results show that, with increasing Reynolds number, there is an initial increase in the wavelength of the wavy vortex flow (λ), and a decrease in the wave speed (c) before they asymptote to a constant value, which is a function of the acceleration. As for the wave amplitude (A), it is found that the effect of acceleration is significant only in a very narrow range of Reynolds numbers. Received: 21 August 2001 / Accepted: 22 November 2001     

Study of Iterative Algorithms for Solution of Dynamic Contact Problems for Elastic Cracked Bodies

A plane problem is solved for the contact interaction between the faces of a rectilinear crack under the action of a normally incident harmonic tension–compression wave. Iterative algorithms are presented to solve the problem for both given initial distribution of contact forces and given initial discontinuity in the displacements of the crack faces. The convergence rates of the algorithms, the maximum contact forces, and displacement discontinuities are compared.     

Contribution to the dynamics of a horizontally unstretchable half-space of the westergaard type

Summary The transfer matrix of a layer is derived from the dynamic equations of a continuum reinforced with very thin unstretchable membranes. For the case of a periodically spaced reinforcement, the finite difference method enables to find closed solutions. The elastic wave propagation in a half-space is investigated, and it is shown that the only possible nonattenuating wave has the velocityc ₂ of the shear wave.     

Interaction of a gasdynamic shock with small disturbances

The stability of shock wave based on the definition of Landau and Lifschitz^[1] is treated in this paper. This is tantamount to solving the problem of interaction of small disturbances with a shock wave. Small disturbances are introduced on both sides of a steady, non-dissipative, plane shock wave. Landau et al.^[1] obtained the stability criterionM ₁>1,M ₂<1 for small disturbances which are travelling in the direction perpendicular to the shock wave. In the present paper, we assume that the small disturbances may be two dimensional, i.e. they may be propagating in the direction inclined to the shock wave. The conclusions obtained are: regardless of whether the incident wave and diverging wave are defined according to the direction of the phase velocity or the group velocity, the shock wave is unstable for some frequencies and longitudinal wave lengths of the disturbances, even if the conditionsM ₁>1,M ₂<1 are fulfilled. Then several experiments are proposed, and the problem of ways to define the incident wave and diverging wave is discussed. The meaning of this problem is illustrated. The same results can be obtained for the steady shock wave in a tube.     

Spectral Analysis of the Balance Equation of Ground Water Hydrology

The spectral analysis of the balance equation of ground water flow, associated to an asymptotic expansion of the conductivity (K) and head (h) fields, permits to show that the high wave number components of the source terms, F, and of the conductivity, K, are attenuated when h is computed as solution to the balance equation. This has important consequences on the inverse mapping (h, F) → K: in fact it is not possible to recover in a reliable way the high wave number components of K, because small errors on the corresponding components of h are amplified so that they can hidden the true signal.     

Mass,momentum and total excess energy transported by a weak planeN wave

A weakly nonlinear plane acoustic wave is emitted into an ideal gas of semi-infinite extent from an infinite plate by its sinusoidal motion of single period. The wave develops into anN wave in the far field, as long as the energy dissipation is negligible everywhere except for discontinuous shock fronts. The third-order effects at shock fronts are evaluated, i.e., the generation of reflected acoustic wave as a result of the interaction of shock and expansion wave and the production of entropy by the energy dissipation at shock fronts. Consideration of these effects enables one to estimate the whole mass, momentum and total excess energy (sum of the kinetic energy and excess of internal energy over an initial undisturbed value) transported by theN wave to the accuracy of third order of wave amplitude. It is shown that the mass and total excess energy transported by theN wave increase and the momentum decreases to asymptotic limits as the wave propagates. The result shows good agreement with a numerical result obtained by solving the Euler equations with a high-resolution TVD upwind scheme.     

Estimates of azimuthal numbers associated with elementary elliptic cylinder wave functions

The paper deals with issues related to the construction of solutions, 2 π-periodic in the angular variable, of the Mathieu differential equation for the circular elliptic cylinder harmonics, the associated characteristic values, and the azimuthal numbers needed to form the elementary elliptic cylinder wave functions. A superposition of the latter is one possible form for representing the analytic solution of the thermoelastic wave propagation problem in long waveguides with elliptic cross-section contour. The classical Sturm-Liouville problem for the Mathieu equation is reduced to a spectral problem for a linear self-adjoint operator in the Hilbert space of infinite square summable two-sided sequences. An approach is proposed that permits one to derive rather simple algorithms for computing the characteristic values of the angular Mathieu equation with real parameters and the corresponding eigenfunctions. Priority is given to the application of the most symmetric forms and equations that have not yet been used in the theory of the Mathieu equation. These algorithms amount to constructing a matrix diagonalizing an infinite symmetric pentadiagonal matrix. The problem of generalizing the notion of azimuthal number of a wave propagating in a cylindrical waveguide to the case of elliptic geometry is considered. Two-sided mutually refining estimates are constructed for the spectral values of the Mathieu differential operator with periodic and half-periodic (antiperiodic) boundary conditions.     

Efficient Algorithms for Modeling the Transport and Biodegradation of Chlorinated Ethenes in Groundwater

Predicting the fate of chlorinated ethenes in groundwater requires the solution of equations that describe both the transport and the biodegradation of the contaminants. Here, we present a model that accounts for (1) transport of chlorinated ethenes in flowing groundwater, (2) mass transfer of contaminants between mobile groundwater and stationary biofilms, and (3) diffusion and biodegradation within the biofilms. Equations for biodegradation kinetics account for biomass growth within the biofilms, the effect of hydrogen on dechlorination, and competitive inhibition between vinyl chloride and cis–dichloroethene. The overall model consists of coupled, non-linear, partial differential equations; solution of such a model is challenging and requires innovative numerical algorithms. We developed and tested two new numerical algorithms to solve the equations in the model; these are called system splitting with operator splitting (SSOS) and system splitting with Picard iteration (SSPI). We discuss the conditions under which one of these algorithms is superior to the other. The contributions of this paper are as follows: first, we believe that the mathematical model presented here is the first transport model that also accounts for diffusion and non-linear biodegradation of chlorinated ethenes in biofilms; second, the SSOS and SSPI are new computational algorithms developed specifically for problems of transport, mass transfer, and non-linear reaction; third, we have identified which of the two new algorithms is computationally more efficient for the case of chlorinated ethenes; and finally, we applied the model to compare the biodegradation behavior under diffusion-limited, metabolism-limited, and hydrogen-limited (donor-limited) conditions.     

Variation in the wave resistance of an amphibian air-cushion vehicle moving over a broken-ice field

The nonstationary rectilinear motion of an amphibian air-cushion vehicle (AACV) on a water surface covered with finely broken ice is considered for various modes of velocity variation. The influence of the water depth, flotation parameters, and mode of motion on the wave resistance of the vehicle is analyzed. Maneuvering methods for increasing or decreasing the wave resistance of AACVs are proposed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 97–102, January–February, 2007.     

The Benjamin–Lighthill Conjecture for Near-Critical Values of Bernoulli’s Constant

In 1954, Benjamin and Lighthill made a conjecture concerning the classical nonlinear problem of steady gravity waves on water of finite depth. According to this conjecture, a point of some cusped region on the (r, s)-plane (r and s are the non-dimensional Bernoulli’s constant and the flow force, respectively), corresponds to every steady wave motion described by the problem. Conversely, at least one steady flow corresponds to every point of the region. In the present paper, this conjecture is proved for near-critical flows (when r attains values close to one), under the assumption that the slopes of wave profiles are bounded. Another question studied here concerns the uniqueness of solutions, and it is proved that for every near-critical value of r only the following waves do exist: (i) a unique (up to translations) solitary wave; (ii) a family of Stokes waves (unique up to translations), which is parametrised by the distance from the bottom to the wave crest. The latter parameter belongs to the interval bounded below by the depth of the subcritical uniform stream and above by the distance from the bottom to the crest of solitary wave corresponding to the chosen value of r.     

Special Cases in Optimization Problems for Stationary Linear Closed-Loop Systems

Consideration is given to problems of solving the algebraic Riccati equation (ARE)—J-factorization of matrix polynomials and J-factorization of rational matrices—to which traditional solution algorithms are not applicable. In this connection, solution algorithms for these problems are discussed where the eigenvalues of the Hamiltonian matrix corresponding to the ARE and the zeros of matrix polynomials are located on the imaginary axis. Moreover, a procedure is set forth for asymptotic expansion of a stabilizing solution of the ARE in the neighborhood of a point at which the ARE has no stabilizing solution. It is shown how this expansion can be used for constructing canonical J-factorization of matrix polynomials that is nearly a noncanonical J-factorization. It is pointed out that the algorithms described can be implemented with the help of MATLAB routines