Computation by the method of characteristics of disturbance in an elastic half-space caused by blast-like surface loading

Mode of access: Internet.

Measurement of water surface roughness and wind shear stress by the use of a pitot tube in a laboratory wave channel /

Mode of access: Internet.

HARPA : a versatile three-dimensional Hamiltonian ray-tracing program for acoustic waves in the atmosphere above irregular terrain /

Mode of access: Internet.

Exponentially decaying pressure pulse moving with superseismic velocity on the surface of a half space of Von Mises elastoplastic material

Mode of access: Internet.

Investigation of small surface protuberances upstream of turbulent boundary layer separation produced by a skewed shock wave at Mach 3 /

"Interim report for period January 1976-February 1976."

Experimental observations of watertable waves in an unconfined aquifer with a sloping boundary

Observations of horizontal and vertical variations in piezometric head in a homogeneous, laboratory aquifer are presented and discussed. The observed fluctuations are induced by a simple harmonic oscillation in the clear water reservoir acting across a sloping boundary. The data qualitatively supports existing theories in that higher harmonics are generated in the active forcing zone and that a significant increase in the inland, asymptotic watertable over height (relative to that found for the vertical boundary case) is observed. The observed overheight is shown to be accurately reproduced by existing small-amplitude perturbation theory. Detailed measurements in the vicinity of the sloping boundary reveal that the signal of generated higher harmonics is strongest near the sand surface and that vertical flows are significant in this region. The aquifer is of finite-depth and is influenced by capillary effects, the experimental data therefore exposes limitations of theories which are based on the assumption of a shallow aquifer free of capillary effects. The dispersive properties of the measured pressure wave in the aquifer are comparable to those found from field observations and likewise do not agree with those predicted by the capillary free, shallow aquifer theory. Although some improvement is obtained, discrepancies between the data and theory persist even when a finite-depth aquifer and capillary effects are considered in the theoretical model. Further sand column experiments eliminate a truncated capillary fringe as a possible contributor to these discrepancies. However, the neglect of horizontal flows in the fringe may have caused the discrepancies. (C) 2004 Elsevier Ltd. All rights reserved.

Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

Surface wave techniques for the study of engineering structures and materials

A wire drive pulse echo method of measuring the spectrum of solid bodies described. Using an 's' plane representation, a general analysis of the transient response of such solids has been carried out. This was used for the study of the stepped amplitude transient of high order modes of disks and for the case where there are two adjacent resonant frequencies. The techniques developed have been applied to the measurenent of the elasticities of refractory materials at high temperatures. In the experimental study of the high order in-plane resonances of thin disks it was found that the energy travelled at the edge of the disk and this initiated the work on one dimensional Rayleigh waves.Their properties were established for the straight edge condition by following an analysis similar to that of the two dimensional case. Experiments were then carried out on the velocity dispersion of various circuits including the disk and a hole in a large plate - the negative curvature condition.Theoretical analysis established the phase and group velocities for these cases and experimental tests on aluminium and glass gave good agreement with theory. At high frequencies all velocities approach that of the one dimensional Rayleigh waves. When applied to crack detection it was observed that a signal burst travelling round a disk showed an anomalous amplitude effect. In certain cases the signal which travelled the greater distance had the greater amplitude.An experiment was designed to investigate the phenanenon and it was established that the energy travelled in two nodes with different velocities.It was found by analysis that as well as the Rayleigh surface wave on the edge, a seoond node travelling at about the shear velocity was excited and the calculated results gave reasonable agreement with the experiments.

Stokes waves revisited:exact solutions in the asymptotic limit

The Stokes perturbative solution of the nonlinear (boundary value dependent) surface gravity wave problem is known to provide results of reasonable accuracy to engineers in estimating the phase speed and amplitudes of such nonlinear waves. The weakling in this structure though is the presence of aperiodic “secular variation” in the solution that does not agree with the known periodic propagation of surface waves. This has historically necessitated increasingly higher-ordered (perturbative) approximations in the representation of the velocity profile. The present article ameliorates this long-standing theoretical insufficiency by invoking a compact exact n-ordered solution in the asymptotic infinite depth limit, primarily based on a representation structured around the third-ordered perturbative solution, that leads to a seamless extension to higher-order (e.g., fifth-order) forms existing in the literature. The result from this study is expected to improve phenomenological engineering estimates, now that any desired higher-ordered expansion may be compacted within the same representation, but without any aperiodicity in the spectral pattern of the wave guides.

Amplitude hysteresis of the surface reactance of a layered superconductor

A new nonlinear electrodynamic phenomenon in layered superconducting slabs irradiated on one side by plane electromagnetic waves in the terahertz range is predicted and studied theoretically. It is shown that the surface reactance of a sample and its reflection coefficient have hysteresis behavior when the amplitude of the incident wave is changed. The analogy between the electrodynamic problem of the electromagnetic field distribution in a superconductor and the mechanical problem of particle motion in a central field is also discussed.

Sub-wavelength resonators as seismic Rayleigh waves shield

We explore the thesis that tall structures can be protected by means of seismic metamaterials. Seismic metamaterials can be built as some elements are created over soil layer with different shapes, dimensions, patterns and from different materials. Resonances in these elements are acting as locally resonant metamaterials for Rayleigh surface waves in the geophysics context. Analytically we proved that if we put infinite chain of SDOF resonator over the soil layer as an elastic, homogeneous and isotropic material, vertical component of Rayleigh wave, longitudinal resonance of oscillators will couple with each other, they would create a Rayleigh bandgap frequency, and wave will experience attenuation before it reaches the structure. As it is impossible to use infinite chain of resonators over soil layer, we considered finite number of resonators throughout our simulations. Analytical work is interpreted using finite element simulations that demonstrates the observed attenuation is due to bandgaps when oscillators are arranged at sub-wavelength scale with respect to the incident Rayleigh wave. For wavelength less than 5 meters, the resulting bandgaps are remarkably large and strongly attenuating when impedance of oscillators matches impedance of soil. Since longitudinal resonance of SDOF resonator are proportional to its length inversely, a formed array of resonators that attenuates Rayleigh waves at frequency ≤10 Hz could be designed starting from vertical pillars coupled to the ground. Optimum number of vertical pillars and their interval spacing called effective area of resonators are investigated. For 10 pillars with effective area of 1 meter and resonance frequency of 4.9 Hz, bandgap frequency causes attenuation and a sinusoidal impulsive force illustrate wave steering down phenomena. Simulation results proved analytical findings of this work.

The effect of ultrasonic waves on the phase behavior of a surfactant-brine-oil system

Acknowledgements The authors would like to gratefully acknowledge and appreciate the School of Engineering, University of Aberdeen, Aberdeen, Scotland, UK, for the provision of the laboratory facilities necessary for completing this work.

The effect of ultrasonic waves on the phase behavior of a surfactant-brine-oil system

Acknowledgements The authors would like to gratefully acknowledge and appreciate the School of Engineering, University of Aberdeen, Aberdeen, Scotland, UK, for the provision of the laboratory facilities necessary for completing this work.

Efficient Computation of Electromagnetic Waves in Hydrocarbon Exploration Using the Improved Numerical Mode Matching (NMM) Method

In this study, we developed and improved the numerical mode matching (NMM) method which has previously been shown to be a fast and robust semi-analytical solver to investigate the propagation of electromagnetic (EM) waves in an isotropic layered medium. The applicable models, such as cylindrical waveguide, optical fiber, and borehole with earth geological formation, are generally modeled as an axisymmetric structure which is an orthogonal-plano-cylindrically layered (OPCL) medium consisting of materials stratified planarly and layered concentrically in the orthogonal directions.

In this report, several important improvements have been made to extend applications of this efficient solver to the anisotropic OCPL medium. The formulas for anisotropic media with three different diagonal elements in the cylindrical coordinate system are deduced to expand its application to more general materials. The perfectly matched layer (PML) is incorporated along the radial direction as an absorbing boundary condition (ABC) to make the NMM method more accurate and efficient for wave diffusion problems in unbounded media and applicable to scattering problems with lossless media. We manipulate the weak form of Maxwell's equations and impose the correct boundary conditions at the cylindrical axis to solve the singularity problem which is ignored by all previous researchers. The spectral element method (SEM) is introduced to more efficiently compute the eigenmodes of higher accuracy with less unknowns, achieving a faster mode matching procedure between different horizontal layers. We also prove the relationship of the field between opposite mode indices for different types of excitations, which can reduce the computational time by half. The formulas for computing EM fields excited by an electric or magnetic dipole located at any position with an arbitrary orientation are deduced. And the excitation are generalized to line and surface current sources which can extend the application of NMM to the simulations of controlled source electromagnetic techniques. Numerical simulations have demonstrated the efficiency and accuracy of this method.

Finally, the improved numerical mode matching (NMM) method is introduced to efficiently compute the electromagnetic response of the induction tool from orthogonal transverse hydraulic fractures in open or cased boreholes in hydrocarbon exploration. The hydraulic fracture is modeled as a slim circular disk which is symmetric with respect to the borehole axis and filled with electrically conductive or magnetic proppant. The NMM solver is first validated by comparing the normalized secondary field with experimental measurements and a commercial software. Then we analyze quantitatively the induction response sensitivity of the fracture with different parameters, such as length, conductivity and permeability of the filled proppant, to evaluate the effectiveness of the induction logging tool for fracture detection and mapping. Casings with different thicknesses, conductivities and permeabilities are modeled together with the fractures in boreholes to investigate their effects for fracture detection. It reveals that the normalized secondary field will not be weakened at low frequencies, ensuring the induction tool is still applicable for fracture detection, though the attenuation of electromagnetic field through the casing is significant. A hybrid approach combining the NMM method and BCGS-FFT solver based integral equation has been proposed to efficiently simulate the open or cased borehole with tilted fractures which is a non-axisymmetric model.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.