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.

Ocean waves across the Arctic: attenuation due to dissipation dominates over scattering for periods longer than 19 s

The poorly understood attenuation of surface waves in sea ice is generally attributed to the combination of scattering and dissipation. Scattering and dissipation have very different effects on the directional and temporal distribution of wave energy, making it possible to better understand their relative importance by analysis of swell directional spreading and arrival times. Here we compare results of a spectral wave model – using adjustable scattering and dissipation attenuation formulations – with wave measurements far inside the ice pack. In this case, scattering plays a negligible role in the attenuation of long swells. Specifically, scattering-dominated attenuation would produce directional wave spectra much broader than the ones recorded, and swell events arriving later and lasting much longer than observed. Details of the dissipation process remain uncertain. Average dissipation rates are consistent with creep effects but are 12 times those expected for a laminar boundary layer under a smooth solid ice plate.

Benefits of using a spectral hindcast database for wave power extraction assessment

When considering deployment of wave energy converters at a given site, it is of prime importance from both a technical and an economical point of view to accurately assess the total yearly energy that can be extracted by the given device. Especially, to be considered is the assessment of the efficiency of the device over the widest span of the sea-states spectral bandwidth. Hence, the aim of this study is to assess the biases and errors introduced on extracted power classically computed using spectral data derived from analytical functions such as a JONSWAP spectrum, compared to the power derived using actual wave spectra obtained from a spectral hindcast database.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

A energia das ondas e o seu aproveitamento

As crises energéticas surgidas no decorrer do último século, incluindo a crise do petróleo, obrigaram o Homem a procurar cada vez mais fontes de energia alternativas e preferencialmente inesgotáveis. Desta situação, resultou uma forte aposta na exploração das fontes de energias renováveis, que são uma das principais alternativas para responder a um aumento de procura, e também, face às exigências de consumos actuais, beneficiando de ao se apostar numa energia limpa e renovável existir uma forte redução nos impactes ambientais que outras fontes de energia não apresentam. O aproveitamento dos recursos provenientes de fontes de energia renováveis para a produção de energia já existe há vários anos, e, em alguns casos, atingiram já um estado de maturidade considerável, como é caso da energia eólica. Em comparação, o mesmo já não acontece com a energia das ondas. Embora o oceano apresente um recurso com enorme potencial para ser explorado, incluindo as ondas e correntes oceânicas, os dispositivos tecnológicos necessários para a exploração deste recurso encontram-se maioritariamente ainda em fase experimental, havendo casos pontuais que atingiram a fase pré-comercial. Assim, não existe até à data um dispositivo padrão para a exploração da energia das ondas em grande escala, contrariamente ao que acontece com a energia eólica. Para esta situação, contribuiu o elevado número de dispositivos patenteados para a exploração da energia das ondas, nenhum deles com vantagens significativas relativamente a outros, e também, devido ao facto de a exploração deste tipo de energia não poder ser feito de igual modo na costa ou a muitos quilómetros dela. Na presente dissertação são apresentados alguns dos principais dispositivos existentes para a extracção de energia proveniente das ondas oceânicas, com especial atenção para os dispositivos de coluna de água oscilante.

Longshore Transport Estimation on Ofir Beach in Northwest Portugal: Sand-Tracer Experiment

This work aims to shed some light on longshore sediment transport (LST) in the highly energetic northwest coast of Portugal. Data achieved through a sand-tracer experiment are compared with data obtained from the original and the new re-evaluated longshore sediment transport formulas (USACE Waterways Experiment Station’s Coastal Engineering and Research Center, Kamphuis, and Bayram bulk formulas) to assess their performance. The field experiment with dyed sand was held at Ofir Beach during one tidal cycle under medium wave-energy conditions. Local hydrodynamic conditions and beach topography were recorded. The tracer was driven southward in response to the local swell and wind- and wave-induced currents (Hsb=0.75mHsb=0.75m, Tp=11.5sTp=11.5s, θb=8−12°θb=8−12°). The LST was estimated by using a linear sediment transport flux approach. The obtained value (2.3×10−3m3⋅s−12.3×10−3m3⋅s−1) approached the estimation provided by the original Bayram formula (2.5×10−3m3⋅s−12.5×10−3m3⋅s−1). The other formulas overestimated the transport, but the estimations resulting from the new re-evaluated formulas also yield approximate results. Therefore, the results of this work indicated that the Bayram formula may give satisfactory results for predicting the longshore sediment transport on Ofir Beach.

A maritimidade portuguesa : a extensão da plataforma continental

As potencialidades da plataforma continental estendida são enormes e variadas, desde as que podem ser obtidas através de sectores tradicionais – como portos e marinas ou turismo náutico – até as que advirão de novos sectores como a exploração dos fundos do mar ou a energia das ondas, entre outras. Com efeito, devido ao alargamento resultante das negociações nas Nações Unidas, é praticamente garantido que Portugal passe a controlar um espaço marítimo, acrescido de 2.1 milhões de km2, isto é, vinte vezes a extensão da sua superfície terrestre. A maritimidade de Portugal constitui um traço profundamente marcado na sua história, bem como nas suas tradicionais práticas económicas, sociais e simbólicas. Na verdade, os últimos anos testemunharam a criação de um amplo consenso quanto ao papel determinante que o desenvolvimento costeiro, a superfície e o comprimento dos limites do “Mar Português” terá no futuro do país. Da articulação entre o Território Nacional, o Oceano Atlântico e o Continente Europeu, associada à confluência das vias de comunicação marítimas que ligam a Europa à América do Norte e do Sul configurar-e-ão os aspectos essenciais do designado “Espaço Estratégico de Interesse Nacional Permanente”. A área geográfica prioritária para o desenvolvimento das acções militares de Defesa Nacional e dos interesses vitais permanentes inclui assim – para além do Território Nacional – o Espaço Interterritorial, o Mar Territorial (MT), o espaço aéreo sob responsabilidade nacional e a Zona Económica Exclusiva (ZEE). O País é titular de soberania no território continental e nos arquipélagos da Madeira dos Açores, para além de estar comprometido com alguns Estados da CPLP, que têm no Atlântico Sul uma importante e poderosa presença. Pela geografia e pela definição dos interesses dos Estados e dos grandes espaços em que estão incluídos, Portugal está na articulação da segurança do Atlântico Norte com a segurança do Atlântico Sul. Esta realidade acarreta novos desafios para a Defesa Nacional, assumindo como preocupações acrescidas o terrorismo, o tráfico e a pirataria marítimas, entre outras.

Seismic Wave Velocity as a Means of In-Place Density Measurement, Part 1, HR-114

Velocity-density tests conducted in the laboratory involved small 4-inch diameter by 4.58-inch-long compacted soil cylinders made up of 3 differing soil types and for varying degrees of density and moisture content, the latter being varied well beyond optimum moisture values. Seventeen specimens were tested, 9 with velocity determinations made along two elements of the cylinder, 180 degrees apart, and 8 along three elements, 120 degrees apart. Seismic energy was developed by blows of a small tack hammer on a 5/8-inch diameter steel ball placed at the center of the top of the cylinder, with the detector placed successively at four points spaced 1/2-inch apart on the side of the specimen involving wave travel paths varying from 3.36 inches to 4.66 inches in length. Time intervals were measured using a model 217 micro-seismic timer in both laboratory and field measurements. Forty blows of the hammer were required for each velocity determination, which amounted to 80 blows on 9 laboratory specimens and 120 blows on the remaining 8 cylinders. Thirty-five field tests were made over the three selected soil types, all fine-grained, using a 2-foot seismic line with hammer-impact points at 6-inch intervals. The small tack hammer and 5/8-inch steel ball was, again, used to develop seismic wave energy. Generally, the densities obtained from the velocity measurements were lower than those measured in the conventional field testing. Conclusions were reached that: (1) the method does not appear to be usable for measurement of density of essentially fine-grained soils when the moisture content greatly exceeds the optimum for compaction, and (2) due to a gradual reduction in velocity upon aging, apparently because of gradual absorption of pore water into the expandable interlayer region of the clay, the seismic test should be conducted immediately after soil compaction to obtain a meaningful velocity value.

The impact of a seasonally ice free Arctic Ocean on the temperature, precipitation and surface mass balance of Svalbard

The observed decline in summer sea ice extent since the 1970s is predicted to continue until the Arctic Ocean is seasonally ice free during the 21st Century. This will lead to a much perturbed Arctic climate with large changes in ocean surface energy flux. Svalbard, located on the present day sea ice edge, contains many low lying ice caps and glaciers and is expected to experience rapid warming over the 21st Century. The total sea level rise if all the land ice on Svalbard were to melt completely is 0.02 m. The purpose of this study is to quantify the impact of climate change on Svalbard’s surface mass balance (SMB) and to determine, in particular, what proportion of the projected changes in precipitation and SMB are a result of changes to the Arctic sea ice cover. To investigate this a regional climate model was forced with monthly mean climatologies of sea surface temperature (SST) and sea ice concentration for the periods 1961–1990 and 2061–2090 under two emission scenarios. In a novel forcing experiment, 20th Century SSTs and 21st Century sea ice were used to force one simulation to investigate the role of sea ice forcing. This experiment results in a 3.5 m water equivalent increase in Svalbard’s SMB compared to the present day. This is because over 50 % of the projected increase in winter precipitation over Svalbard under the A1B emissions scenario is due to an increase in lower atmosphere moisture content associated with evaporation from the ice free ocean. These results indicate that increases in precipitation due to sea ice decline may act to moderate mass loss from Svalbard’s glaciers due to future Arctic warming.

Projection of global wave climate change towards the end of the 21st century

Wind generated waves at the sea surface are of outstanding importance for both their practical relevance in many aspects, such as coastal erosion, protection, or safety of navigation, and for their scientific relevance in modifying fluxes at the air-sea interface. So far long-term changes in ocean wave climate have been studied mostly from a regional perspective with global dynamical studies emerging only recently. Here a global wave climate study is presented, in which a global wave model (WAM) is driven by atmospheric forcing from a global climate model (ECHAM5) for present day and potential future climate conditions represented by the IPCC (Intergovernmental Panel for Climate Change) A1B emission scenario. It is found that changes in mean and extreme wave climate towards the end of the twenty-first century are small to moderate, with the largest signals being a poleward shift in the annual mean and extreme significant wave heights in the mid-latitudes of both hemispheres, more pronounced in the Southern Hemisphere, and most likely associated with a corresponding shift in mid-latitude storm tracks. These changes are broadly consistent with results from the few studies available so far. The projected changes in the mean wave periods, associated with the changes in the wave climate in the mid to high latitudes, are also shown, revealing a moderate increase in the equatorial eastern side of the ocean basins. This study presents a step forward towards a larger ensemble of global wave climate projections required to better assess robustness and uncertainty of potential future wave climate change.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Self-consistent wave-particle interactions in dispersive scale long-period field-line-resonances

Using 1D Vlasov drift-kinetic computer simulations, it is shown that electron trapping in long period standing shear Alfven waves (SAWs) provides an efficient energy sink for wave energy that is much more effective than Landau damping. It is also suggested that the plasma environment of low altitude auroral-zone geomagnetic field lines is more suited to electron acceleration by inertial or kinetic scale Alfven waves. This is due to the self-consistent response of the electron distribution function to SAWs, which must accommodate the low altitude large-scale current system in standing waves. We characterize these effects in terms of the relative magnitude of the wave phase and electron thermal velocities. While particle trapping is shown to be significant across a wide range of plasma temperatures and wave frequencies, we find that electron beam formation in long period waves is more effective in relatively cold plasma.

Wave spectra assimilation in typhoon wave modeling for the East China Sea

The East China Sea is a hot area for typhoon waves to occur. A wave spectra assimilation model has been developed to predict the typhoon wave more accurately and operationally. This is the first time where wave data from Taiwan have been used to predict typhoon wave along the mainland China coast. The two-dimensional spectra observed in Taiwan northeast coast modify the wave field output by SWAN model through the technology of optimal interpolation (OI) scheme. The wind field correction is not involved as it contributes less than a quarter of the correction achieved by assimilation of waves. The initialization issue for assimilation is discussed. A linear evolution law for noise in the wave field is derived from the SWAN governing equations. A two-dimensional digital low-pass filter is used to obtain the initialized wave fields. The data assimilation model is optimized during the typhoon Sinlaku. During typhoons Krosa and Morakot, data assimilation significantly improves the low frequency wave energy and wave propagation direction in Taiwan coast. For the far-field region, the assimilation model shows an expected ability of improving typhoon wave forecast as well, as data assimilation enhances the low frequency wave energy. The proportion of positive assimilation indexes is over 81% for all the periods of comparison. The paper also finds that the impact of data assimilation on the far-field region depends on the state of the typhoon developing and the swell propagation direction.

SURFACE ALFVEN WAVE DAMPING IN A THREE-DIMENSIONAL SIMULATION OF THE SOLAR WIND

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in the regions of strong inhomogeneities in density or magnetic field (e.g., the border between open and closed magnetic field lines). Using a three-dimensional (3D) magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1 and 4 R(circle dot) (solar radii), the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere: 3 x 10(-6) to 10(-1) Hz. In the region between open and closed field lines, within a few R(circle dot) of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We demonstrate that projection effects from the plane of the sky to 3D are significant in the calculation of field line expansion. We determine that waves with frequencies >2.8 x 10(-4) Hz are damped between 1 and 4 R(circle dot). In quiet-Sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find it as an order of magnitude larger than needed for quiet-Sun regions. For active regions, the contribution to the heating is 20%. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self-consistently as an energy driver for the wind in global MHD models.