Development of a numerical wave tank for analysis of nonlinear and irregular wave field

A numerical wave-absorption filter has been developed for an open boundary condition in the analysis of nonlinear and irregular wave evolution. The filter is composed of a simulated sponge layer and Sommerfeld's radiation condition at the outer edge of the layer. The wave-absorption characteristics of the filter have been investigated by applying the linear potential theory and a two-dimensional nonlinear boundary element model. In both cases, the filter is found to he applicable for a wide range of wave parameters. In order to realize an idealized “numerical wave tank”, the present model also incorporates a nonreflective wave generator in the computational domain composed of a series of vertically aligned point sources. Numerous numerical experiments demonstrate that the present approach is effective in generating an arbitrary wave profile without reflection not only at the open boundaries but also at the wave generator.     

The application of the Lattice Boltzmann method to the one-dimensional modeling of pulse waves in elastic vessels

The one-dimensional nonlinear equations for the blood flow motion in distensible vessels are considered using the kinetic approach. It is shown that the Lattice Boltzmann (LB) model for non-ideal gas is asymptotically equivalent to the blood flow equations for compliant vessels at the limit of low Knudsen numbers. The equations of state for non-ideal gas are transformed to the pressure-luminal area response. This property allows to model arbitrary pressure-luminal area relations. Several test problems are considered: the propagation of a sole nonlinear wave in an elastic vessel, the propagation of a pulse wave in a vessel with varying mechanical properties (artery stiffening) and in an artery bifurcation, in the last problem Resistor–Capacitor–Resistor (RCR) boundary conditions are considered. The comparison with the previous results shows a good precision.     

Tests of a new numerical scheme on a “non-reflection and free-transmission” open-sea boundary for longwaves

A new numerical scheme of a “non-reflection and free-transmission” boundary for longwave equations proposed by Hino (1987) has been tested for a variety of cases. The test results verify the effectiveness of the method for (a) a single progressive wave train on a horizontal bottom, (b) two wave trains each propagating in opposite directions on a horizontal bottom, (c) a single wave train propagating on a sloping bottom with friction, (d) oscillatory flood waves in an open channel flow, (e) two-dimensional waves travelling obliquely to open boundaries and (f) water surface oscillation in a harbor by waves incident through an opening.     

The cumulative interaction of finite amplitude waves in strings

A “two time scale” asymptotic expansion procedure describing the modulation of a propagating simple wave governed by a system of non-linear partial differential equations is applied to the deflection waves of non-linear elastic strings. Rapid deflection signals propagating into a general slowly varying disturbance are modulated. In addition, they themselves affect the equations for that disturbance. The two effects are separated naturally when, to prevent the cumulative growth inherent in most “high frequency” procedures, an averaging technique is introduced. The interaction of two deflection waves is given as a specific example.     

Damping of large-amplitude solitary waves

We consider the damping of large-amplitude solitary waves in the framework of the extended Korteweg-de Vries equation (that is, the usual Korteweg-de Vries equation supplemented with a cubic nonlinear term) modified by the inclusion of a small damping term. The damping of a solitary wave is studied for several different forms of friction, using both the analytical adiabatic asymptotic theory and numerical simulations. When the coefficient of the cubic nonlinear term has the opposite sign to the coefficient of the linear dispersive term, the extended Kortweg-de Vries equation can support large-amplitude “thick” solitary waves. Under the influence of friction, these “thick” solitary waves decay and may produce one or more secondary solitary waves in this process. On the other hand, when the coefficient of the cubic nonlinear term has the same sign as the coefficient of the linear dispersive term, but the opposite sign to the coefficient of the quadratic nonlinear term, the action of friction may cause a solitary wave to decay into a wave packet.     

Stability analysis of annular flow structure, using a discrete form of the “two-fluid model”

The differential form of the “two-fluid model” for annular flow, neglecting surface tension, is ill-posed, and it is not suited for examining the stability of the steady-state solutions with respect to the average film thickness. It is shown here that a discrete (difference) representation of the two-fluid model may lead to an appropriate criterion for the stability of the steady-state solutions. Exactly the same criterion is obtained from the requirement that the kinematic waves will propagate in the downstream direction. The suggested discrete form of the “two-fluid model” is used to perform transient simulation and for examining the system response to finite disturbances.     

Weakly nonlocal envelope solitary waves: numerical calculations for the Klein-Gordon (φ) equation

“Weakly nonlocal” solitary waves differ from ordinary solitary waves by possessing small amplitude, oscillatory “wings” that extend indefinitely from the large amplitude “core”. Such generalized solitary waves have been discovered in capillarygravity water waves, particle physics models, and geophysical Rossby waves. In this work, we present explicit calculations of weakly nonlocal envelope solitary waves. Each is a sine wave modulated by a slowly-varying “envelope” that itself propagates at the group velocity. Our example is the cubically nonlinear Klein-Gordon equation, which is a model in particle physics (φ⁴ theory) and in electrical engineering (with a different sign). Both cases have weakly nonlocal“breather” solitons. Via the Lorentz invariance, each breather generates a one-parameter family of nonlocal envelope solitary waves. The φ⁴ breather was described and calculated in earlier work. This generates envelope solitons which have “wings” that are (mostly) proportional to the second harmonic of the sinusoidal factor. In this article, we calculate breathers and envelope solitary waves for the second, electrical engineering case. Since these, unlike the φ⁴ waves, contain only odd harmonics, the envelope solitary waves are nonlocal only via the third harmonic.     

Mach Wave Effect on Laminar-Turbulent Transition in Supersonic Flow over a Flat Plate

The effect of a Mach wave (N wave) on laminar-turbulent transition induced by the first instability mode (Tollmien–Schlichting wave) in the flat-plate boundary layer is investigated on the basis of the numerical solution of Navier–Stokes equations at the freestream Mach number of 2.5. In accordance with the experiment, the N wave is generated by a two-dimensional roughness at the computation domain boundary corresponding to the side wall of the test section of a wind tunnel. It is shown that the disturbance induced by the backward front of the N wave in the boundary layer has no effect on the beginning of transition but displaces downstream the nonlinear stage of the first mode development. The disturbance induced by the forward front of the N wave displaces the beginning of transition upstream.     

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.     

The moving-load problem in soil dynamics—the vertical displacement approximation

The moving point load problem in soil dynamics is analyzed in the vertical particle displacement approximation. Prior to its motion, the load is stationary. From the instant at which it is set into motion it moves, with constant speed, along a straight path on the (horizontal) planar surface of a semi-infinite elastic medium. The modified Cagniard method for solving transient wave problems is used to determine closed-form expressions for the vertical component of the particle displacement from the elastodynamic wave equation of which only the vertical component is taken into account. The relevant approximation is standard in soil dynamics. Both the cases of “subsonic” and “supersonic” surface load speeds are considered. Methods to include losses in the model are briefly discussed. The study has been initiated with a view to the application of the results to the analysis of the ground motion generated by high-speed trains traveling on a poorly consolidated soil.     

Investigation of a combustion driven oscillation in a refinery flare: Part B: Visualisation of a periodic flow instability in a bifurcating duct following a contraction

A flow visualisation study was performed to investigate a periodic flow instability in a bifurcating duct within the tip of the flares at the Shell refinery in Clyde, NSW, to verify the trigger of a combustion-driven oscillation proposed in Part A of this study, and to identify its features. The model study assessed only the flow instability, uncoupled from the acoustic resonance and the combustion that are also present in the actual flare. Three strong, coupled flow oscillations were found to be present in three regions of the fuel line in the flare tip model. A periodic flow separation was found to occur within the contraction at the inlet to the tip, a coupled, periodic flow oscillation was found in the two transverse “cross-over ducts” from the central pipe to the outer annulus and an oscillating flow recirculation was present in the “end-cap” region of the central pipe. The dimensionless frequency of these oscillations in the model was found to match that measured in the full-scale plant for high fuel flow rates. This, and the strength of these flow oscillations, gives confidence that they are integral to the full-scale combustion-driven oscillation and most likely the primary trigger. The evidence indicates that the periodic flow instability is initiated by the separation and roll-up of the annular boundary layer at the start of the contraction in the fuel section of the flare tip. The separation generates an annular vortex which interacts with the blind-ended pipe downstream, leading to a pressure wave which propagates back upstream, initiating the next separation event and repeating the cycle. The study also investigated flow control devices with a view to finding a practical approach to mitigate the oscillations. The shape of these devices was constrained to allow installation without removing the tip of the flare. This aspect of the study highlighted the strength and nature of the coupled oscillation, since it proved to be very difficult to mitigate the oscillation in this way. An effective configuration is presented, comprising of three individual components, all three of which were found to be necessary to eliminate the oscillation completely.     

Elastic wave scattering by a three-dimensional inhomogeneity in an elastic half space

In this paper we will consider scattering of elastic waves in a half space. The half space is an isotropic, linear and homogeneous medium except for a finite inhomogeneity. The T-matrix method (also called the “extended boundary condition method” or “null field approach”) is extended to derive expressions for the elastic field inside the half space and the surface field on the interface. The assumptions on the source that excites the half space are fairly weak. In the numerical applications found in this paper we assume a Rayleigh surface wave to be the incoming field, and we only compute the surface displacements. We make illustrations on some simple types of scatterers (spheres and spheroids; the latter ones can be arbitrarily oriented).     

Permanent Fronts in Two-Phase Flows in a Porous Medium

A family of exact solutions for a model of a one-dimensional horizontal flow of two immiscible, incompressible fluids in a porous medium, including the effects of capillary pressure, is obtained analytically by solving the governing singular parabolic nonlinear diffusion equation. Each solution has the form of a permanent front propagating with a constant velocity. It is shown that, for every propagation velocity, there exists a set of permanent fronts all of which are moving with this velocity in an inflowing wetting–outflowing non-wetting flow configuration. Global bifurcations of this set, with the front velocity as a bifurcation parameter, are investigated analytically and numerically in detail in the case when the permeabilities and the capillary pressure are linear functions of the wetting phase saturation. Main results for the nonlinear Brooks–Corey model are also presented. In both models three global bifurcations occur. By using a geometric dynamical system approach, the nonlinear stability of the permanent fronts is established analytically. Based on the permanent front solutions, an interpretation of the dynamics of an arbitrary front of finite extent in the model is given as follows. The instantaneous upstream (downstream) velocity of an arbitrary non-quasistationary front is equal to the velocity of a permanent front whose shape coincides up to two leading orders with the instantaneous shape of the non-quasistationary front at the upstream (respectively, downstream) location. The upstream and downstream locations of the front undergo instantaneous translations governed by modified nonsingular hyperbolic equations. The portion of the front in between these locations undergoes a diffusive redistribution governed by a nonsingular nonlinear parabolic diffusion equation. We have proposed a numerical approach based on a parabolic–hyperbolic domain decomposition for computing non-quasistationary fronts.     

Numerical simulation of fully nonlinear irregular wave tank in three dimension

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

Propagation characteristics of pressure disturbances originated by gas jets in fluidized beds

An experimental investigation has been carried out on velocities and amplitudes of pressure disturbances in fluidized beds made of 100–200 μm glass ballotini. Disturbances were originated by gas jetting in a 0.35 m i.d. fluidized bed. A fluidization tube 0.10 m i.d. has also been used. Different types of disturbances have been induced in the bed contained in this tube: injection of a freely rising bubble and of a captive bubble; injection of a bubble chain; and compression of the bed free surface. The dynamic wave character of the disturbances has been shown. Velocities and amplitudes of waves moving through the beds have been measured. In particular, wave velocities have been compared with theoretical results obtained by the application of “pseudo-homogeneous” and “separated phase flow” models.     

A new Finite Element Method for magneto-dynamical problems: two-dimensional results

A mixed Lagrange finite element technique is used to solve the Maxwell equations in the magneto-hydrodynamic (MHD) limit in an hybrid domain composed of vacuum and conducting regions. The originality of the approach is that no artificial boundary condition is enforced at the interface between the conducting and the insulating regions and the non-conducting medium is not approximated by a weakly conducting medium as is frequently done in the literature. As a first evaluation of the performance of the method, we study two-dimensional (2D) configurations, where the flow streamlines of the conducting fluid are planar, i.e., invariant in one direction, and either the magnetic field (“magnetic scalar” case) or the electric field (“electric scalar” case) is parallel to the invariant direction. Induction heating, eddy current generation, and magnetic field stretching are investigated showing the usefulness of finite element methods to solve magneto-dynamical problems with complex insulating boundaries.     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

Subharmonic generation in physical systems: An interaction-box approach

A so-called “interaction-box” formalism, which has recently been introduced to describe hysteresis in dynamical systems in the case of higher harmonic generation, is further discussed and generalized to describe the phenomenon of subharmonic generation. In this case, the increase in the periodicity of the response is reflected in the formation of multiple loops in the Effect (output) vs. Cause (input) diagrams. Conversely, we show how this type of response represents a sort of “signature” of the system, and can thus be employed to draw general conclusions about the features of the latter. A specific example of a nonlinear system is chosen to illustrate the approach, namely a vibrating cantilever beam with a breathing crack. Effect vs. Cause curves are calculated for this system in the presence of higher harmonics and subharmonics.     

DYNAMICS AND STABILITY OF PINNED–CLAMPED AND CLAMPED–PINNED CYLINDRICAL SHELLS CONVEYING FLUID

The paper examines the dynamics and stability of fluid-conveying cylindrical shells having pinned–clamped or clamped–pinned boundary conditions, where “pinned” is an abbreviation for “simply supported”. Flügge's equations are used to describe the shell motion, while the fluid-dynamic perturbation pressure is obtained utilizing the linearized potential flow theory. The solution is obtained using two methods — the travelling wave method and the Fourier-transform approach. The results obtained by both methods suggest that the negative damping of the clamped–pinned systems and positive damping of the pinned–clamped systems, observed by previous investigators for any arbitrarily small flow velocity, are simply numerical artefacts; this is reinforced by energy considerations, in which the work done by the fluid on the shell is shown to be zero. Hence, it is concluded that both systems are conservative.