Derivation of Wave Loads for the Design of Oscillating Wave Surge Converters

This paper deals with the problem of estimating wave pressure loads acting on Oscillating Wave Surge Converters (OWSC) for assessment of fatigue on their components. Recent wave loading data issued from experimental testing of a 25th scale model of a box-shaped OWSC are here used to review the accuracy of the predictions made by an engineering method previously developed to derive wave pressure loads on OWSCs from experimental data. Predictions are shown underestimate wave pressure loads, and other methods subsequently developed are presented. A simplistic experimental method taking in consideration variations of the wetted surface area of the flap is shown to lead to relatively good estimates of wave pressure loads that could be used for fatigue calculations.

PHOTOSIGNAL ENHANCEMENT IN AL-GAAS DIODES DUE TO SURFACE-PLASMONS AND GUIDED-WAVE MODES

Light of wavelength 632.8 nm and p-polarization is incident on a prism-air gap (varied from 0.7 to 7 mum)-Al-GaAs arrangement. Both the photosignal generated by the Schottky diode and the reflectance are measured as a function of the internal angle of incidence in the prism. There is significant, well-defined enhancement of the photosignal, up to a factor of approximately 7.5, associated with two different types of enhanced absorption modes. For air gaps <1.5 mum there is photosignal enhancement due to an enhanced absorption feature (reflectance dip) that occurs at an angle of incidence just above critical angle in the prism; this feature corresponds to the excitation of a surface plasmon polariton at the Al-air interface. For air gaps > 1 mum there are between one and ten photoresponse peaks at input angles less than the critical angle. The corresponding enhanced absorption features are due to leaky guided wave modes set up in the air gap.

The Vertical Distribution and Evolution of Slam Pressure on an Oscillating Wave Surge Converter

The accurate definition of the extreme wave loads which act on offshore structures represents a significant challenge for design engineers and even with decades of empirical data to base designs upon there are still failures attributed to wave loading. The environmental conditions which cause these loads are infrequent and highly non-linear which means that they are not well understood or simple to describe. If the structure is large enough to affect the incident wave significantly further non-linear effects can influence the loading. Moreover if the structure is floating and excited by the wave field then its responses, which are also likely to be highly non-linear, must be included in the analysis. This makes the description of the loading on such a structure difficult to determine and the design codes will often suggest employing various tools including small scale experiments, numerical and analytical methods, as well as empirical data if available.
Wave Energy Converters (WECs) are a new class of offshore structure which pose new design challenges, lacking the design codes and empirical data found in other industries. These machines are located in highly exposed and energetic sites, designed to be excited by the waves and will be expected to withstand extreme conditions over their 25 year design life. One such WEC is being developed by Aquamarine Power Ltd and is called Oyster. Oyster is a buoyant flap which is hinged close to the seabed, in water depths of 10 to 15m, piercing the water surface. The flap is driven back and forth by the action of the waves and this mechanical energy is then converted to electricity.
It has been identified in previous experiments that Oyster is not only subject to wave impacts but it occasionally slams into the water surface with high angular velocity. This slamming effect has been identified as an extreme load case and work is ongoing to describe it in terms of the pressure exerted on the outer skin and the transfer of this short duration impulsive load through various parts of the structure.
This paper describes a series of 40th scale experiments undertaken to investigate the pressure on the face of the flap during the slamming event. A vertical array of pressure sensors are used to measure the pressure exerted on the flap. Characteristics of the slam pressure such as the rise time, magnitude, spatial distribution and temporal evolution are revealed. Similarities are drawn between this slamming phenomenon and the classical water entry problems, such as ship hull slamming. With this similitude identified, common analytical tools are used to predict the slam pressure which is compared to that measured in the experiment.

Investigation of Utilizing Single Surface Depressions in Optimization of NLF Aerofoil Design

An environment has been created for the optimisation of aerofoil profiles with inclusion of small surface features. For TS wave dominated flows, the paper examines the consequences of the addition of a depression on the aerodynamic optimisation of an NLF aerofoil, and describes the geometry definition fidelity and optimisation algorithm employed in the development process. The variables that define the depression for this optimisation investigation have been fixed, however a preliminary study is presented demonstrating the sensitivity of the flow to the depression characteristics. Solutions to the optimisation problem are then presented using both gradient-based and genetic algorithm techniques, and for accurate representation of the inclusion of small surface perturbations it is concluded that a global optimisation method is required for this type of aerofoil optimisation task due to the nature of the response surface generated. When dealing with surface features, changes in the transition onset are likely to be of a non-linear nature so it is highly critical to have an optimisation algorithm that is robust, suggesting that for this framework, gradient-based methods alone are not suited.

Surface properties and catalytic activity of ferrospinels of nickel, cobalt and copper, prepared by soft chemical methods

Ferrospinels of nickel, cobalt and copper and their sulphated analogues were prepared by the room temperature coprecipitation route to yield samples with high surface areas. The intrinsic acidity among the ferrites was found to decrease in the order: cobalt> nickel> copper. Sulphation caused an increase in the number of weak and medium strong acid sites, whereas the strong acid sites were left unaffected. Electron donor studies revealed that copper ferrite has both the highest proportion of strong sites and the lowest proportion of weak basic sites. All the ferrite samples proved to be good catalysts for the benzoy lation of toluene with benzoyl chloride. copper and cobalt ferrites being much more active than nickel ferrite. The catalytic activity for benzoylation was not much influenced by sulphation, but it increased remarkably with calcination temperature of the catalyst. Surface Lewis acid sites, provided by the octahedral cations on the spinel surface, are suggested to be responsible for the catalytic activity for the benzoylation reaction.

Surface studies of some metallic thin films by optical methods

The thesis aims to present the results of experimental investigations on the changes of optical properties of metallic thin films due to heating. The parameters which are measured are reflectivity, refractive indices and the ellipsometric quantities V and A . The materials used in the studies are metals like Silver, Aluminium and Copper. By applying the optical method the interdiffusion taking place in multilayer ‘films of Aluminium and Silver has also been studied. Special interest has been taken to reveal the mechanisms of the hillock growth and surface roughness caused by heating and their relation with the stress in the film

Efficient evaluation of highly oscillatory acoustic scattering surface integrals

We consider the approximation of some highly oscillatory weakly singular surface integrals, arising from boundary integral methods with smooth global basis functions for solving problems of high frequency acoustic scattering by three-dimensional convex obstacles, described globally in spherical coordinates. As the frequency of the incident wave increases, the performance of standard quadrature schemes deteriorates. Naive application of asymptotic schemes also fails due to the weak singularity. We propose here a new scheme based on a combination of an asymptotic approach and exact treatment of singularities in an appropriate coordinate system. For the case of a spherical scatterer we demonstrate via error analysis and numerical results that, provided the observation point is sufficiently far from the shadow boundary, a high level of accuracy can be achieved with a minimal computational cost.

On integral equation and least squares methods for scattering by diffraction gratings

In this paper we consider the scattering of a plane acoustic or electromagnetic wave by a one-dimensional, periodic rough surface. We restrict the discussion to the case when the boundary is sound soft in the acoustic case, perfectly reflecting with TE polarization in the EM case, so that the total field vanishes on the boundary. We propose a uniquely solvable first kind integral equation formulation of the problem, which amounts to a requirement that the normal derivative of the Green's representation formula for the total field vanish on a horizontal line below the scattering surface. We then discuss the numerical solution by Galerkin's method of this (ill-posed) integral equation. We point out that, with two particular choices of the trial and test spaces, we recover the so-called SC (spectral-coordinate) and SS (spectral-spectral) numerical schemes of DeSanto et al., Waves Random Media, 8, 315-414 1998. We next propose a new Galerkin scheme, a modification of the SS method that we term the SS* method, which is an instance of the well-known dual least squares Galerkin method. We show that the SS* method is always well-defined and is optimally convergent as the size of the approximation space increases. Moreover, we make a connection with the classical least squares method, in which the coefficients in the Rayleigh expansion of the solution are determined by enforcing the boundary condition in a least squares sense, pointing out that the linear system to be solved in the SS* method is identical to that in the least squares method. Using this connection we show that (reflecting the ill-posed nature of the integral equation solved) the condition number of the linear system in the SS* and least squares methods approaches infinity as the approximation space increases in size. We also provide theoretical error bounds on the condition number and on the errors induced in the numerical solution computed as a result of ill-conditioning. Numerical results confirm the convergence of the SS* method and illustrate the ill-conditioning that arises.

Solving surface structures from normal incidence X-ray standing wave data

A program is provided to determine structural parameters of atoms in or adsorbed on surfaces by refinement of atomistic models towards experimentally determined data generated by the normal incidence X-ray standing wave (NIXSW) technique. The method employs a combination of Differential Evolution Genetic Algorithms and Steepest Descent Line Minimisations to provide a fast, reliable and user friendly tool for experimentalists to interpret complex multidimensional NIXSW data sets.

Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the $p$-version

Plane wave discontinuous Galerkin (PWDG) methods are a class of Trefftz-type methods for the spatial discretization of boundary value problems for the Helmholtz operator $-\Delta-\omega^2$, $\omega>0$. They include the so-called ultra weak variational formulation from [O. Cessenat and B. Després, SIAM J. Numer. Anal., 35 (1998), pp. 255–299]. This paper is concerned with the a priori convergence analysis of PWDG in the case of $p$-refinement, that is, the study of the asymptotic behavior of relevant error norms as the number of plane wave directions in the local trial spaces is increased. For convex domains in two space dimensions, we derive convergence rates, employing mesh skeleton-based norms, duality techniques from [P. Monk and D. Wang, Comput. Methods Appl. Mech. Engrg., 175 (1999), pp. 121–136], and plane wave approximation theory.

The domain derivative in rough-surface scattering and rigorous estimates for first-order perturbation theory

We investigate Fréchet differentiability of the scattered field with respect to variation in the boundary in the case of time–harmonic acoustic scattering by an unbounded, sound–soft, one–dimensional rough surface. We rigorously prove the differentiability of the scattered field and derive a characterization of the Fréchet derivative as the solution to a Dirichlet boundary value problem. As an application of these results we give rigorous error estimates for first–order perturbation theory, justifying small perturbation methods that have a long history in the engineering literature. As an application of our rigorous estimates we show that a plane acoustic wave incident on a sound–soft rough surface can produce an unbounded scattered field.

Plane wave discontinuous Galerkin methods: exponential convergence of the hp-version

We consider the two-dimensional Helmholtz equation with constant coefficients on a domain with piecewise analytic boundary, modelling the scattering of acoustic waves at a sound-soft obstacle. Our discretisation relies on the Trefftz-discontinuous Galerkin approach with plane wave basis functions on meshes with very general element shapes, geometrically graded towards domain corners. We prove exponential convergence of the discrete solution in terms of number of unknowns.

A priori error analysis of space–time Trefftz discontinuous Galerkin methods for wave problems

We present and analyse a space–time discontinuous Galerkin method for wave propagation problems. The special feature of the scheme is that it is a Trefftz method, namely that trial and test functions are solution of the partial differential equation to be discretised in each element of the (space–time) mesh. The method considered is a modification of the discontinuous Galerkin schemes of Kretzschmar et al. (2014) and of Monk & Richter (2005). For Maxwell’s equations in one space dimension, we prove stability of the method, quasi-optimality, best approximation estimates for polynomial Trefftz spaces and (fully explicit) error bounds with high order in the meshwidth and in the polynomial degree. The analysis framework also applies to scalar wave problems and Maxwell’s equations in higher space dimensions. Some numerical experiments demonstrate the theoretical results proved and the faster convergence compared to the non-Trefftz version of the scheme.

Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki

The performance of three urban land surface models, run in offline mode, with their default external parameters, is evaluated for two distinctly different sites in Helsinki: Torni and Kumpula. The former is a dense city centre site with 22% vegetation, while the latter is a suburban site with over 50% vegetation. At both locations the models are compared against sensible and latent heat fluxes measured using the eddy covariance technique, along with snow depth observations. The cold climate experienced by the city causes strong seasonal variations that include snow cover and stable atmospheric conditions. Most of the time the three models are able to account for the differences between the study areas as well as the seasonal and diurnal variability of the energy balance components. However, the performances are not systematic across the modelled components, season and surface type. The net all-wave radiation is well simulated, with the greatest uncertainties related to snowmelt timing, when the fraction of snow cover has a key role, particularly in determining the surface albedo. For the turbulent fluxes, more variation between the models is seen which can partly be explained by the different methods in their calculation and partly by surface parameter values. For the sensible heat flux, simulation of wintertime values was the main problem, which also leads to issues in predicting near-surface stabilities particularly at the dense city centre site. All models have the most difficulties in simulating latent heat flux. This study particularly emphasizes that improvements are needed in the parameterization of anthropogenic heat flux and thermal parameters in winter, snow cover in spring and evapotranspiration in order to improve the surface energy balance modelling in cold climate cities.