Development of a high-frequency surface-wave radar cross section model for icebergs

In this thesis, the first-order radar cross section (RCS) of an iceberg is derived and simulated. This analysis takes place in the context of a monostatic high frequency surface wave radar with a vertical dipole source that is driven by a pulsed waveform. The starting point of this work is a general electric field equation derived previ- ously for an arbitrarily shaped iceberg region surrounded by an ocean surface. The condition of monostatic backscatter is applied to this general field equation and the resulting expression is inverse Fourier transformed. In the time domain the excitation current of the transmit antenna is specified to be a pulsed sinusoid signal. The result- ing electric field equation is simplified and its physical significance is assessed. The field equation is then further simplified by restricting the iceberg's size to fit within a single radar patch width. The power received by the radar is calculated using this electric field equation. Comparing the received power with the radar range equation gives a general expression for the iceberg RCS. The iceberg RCS equation is found to depend on several parameters including the geometry of the iceberg, the radar frequency, and the electrical parameters of both the iceberg and the ocean surface. The RCS is rewritten in a form suitable for simulations and simulations are carried out for rectangularly shaped icebergs. Simulation results are discussed and are found to be consistent with existing research.

Analisi del catalogo Europeo-Mediterraneo degli RCMT: 18 anni di dati.

L'RCMT (Regional Centroid Moment Tensor), realizzato e gestito dai ricercatori dell'INGV (Istituto Nazionale di Geofisica e Vulcanologia), è dal 1997 il catalogo di riferimento per gli eventi sismici avvenuti nell'area Europea-Mediterranea, ossia nella regione avente longitudine compresa tra 10° W e 40° E e latitudine compresa tra 25° N e 60° N. Tale regione è caratterizzata da un'attività tettonica complessa, legata non soltanto alla convergenza delle placche Euroasiatica ed Africana, ma anche al movimento di altre placche minori (ad esempio, la placca Arabica), che tutte insieme danno origine ad una vasta gamma di regimi tettonici. Col termine RCMT si indica un particolare tipo di tensore momento sismico, la cui determinazione avviene su scala regionale, per eventi sismici aventi M_w >= 4.5 (M_w >= 4.0 per gli eventi che avvengono nella penisola italica). Il tensore momento sismico è uno strumento fondamentale per caratterizzare natura ed entità di un terremoto. Da esso, infatti, oltre alla magnitudo momento M_w, si ricava anche il meccanismo focale. Comunemente rappresentato sotto forma di beach ball, consente di individuare il tipo di movimento (distensivo, compressivo o trascorrente, o anche una combinazione del primo o del secondo con il terzo) avvenuto sulla faglia che ha provocato il terremoto. I tensori momento sismico permettono, quindi, di identificare le faglie che si attivano durante una sequenza sismica, di comprendere la loro cinematica e di ipotizzare la successiva evoluzione a breve termine. Scopo di questa relazione di laurea è stato derivare le relazioni che intercorrono fra le M_w dell'RCMT e le M_w del CMT (Centroid Moment Tensor della Columbia University), del GFZ (Deutsches GeoForschungsZentrum di Postdam) e del TDMT (Time Domain Moment Tensor). Le relazioni sono state ottenute applicando il metodo dei minimi quadrati agli eventi comuni, che sono stati selezionati utilizzando alcuni semplici programmi scritti in Fortran.

Modelagem e análise dinâmica de um sistema de conversão de energia eólica dotado de gerador síncrono de imã permanente utilizando a plataforma ATP

Human development requires a broad balance between ecological, social and economic factors in order to ensure its own sustainability. In this sense, the search for new sources of energy generation, with low deployment and operation costs, which cause the least possible impact to the environment, has been the focus of attention of all society segments. To do so, the reduction in exploration of fossil fuels and the encouragement of using renewable energy resources for distributed generation have proved interesting alternatives to the expansion of the energy matrix of various countries in the world. In this sense, the wind energy has acquired an increasingly significant role, presenting increasing rates of power grid penetration and highlighting technological innovations such as the use of permanent magnet synchronous generators (PMSG). In Brazil, this fact has also been noted and, as a result, the impact of the inclusion of this source in the distribution and sub-transmission power grid has been a major concern of utilities and agents connected to Brazilian electrical sector. Thus, it is relevant the development of appropriate computational tools that allow detailed predictive studies about the dynamic behavior of wind farms, either operating with isolated load, either connected to the main grid, taking also into account the implementation of control strategies for active/reactive power generation and the keeping of adequate levels of voltage and frequency. This work fits in this context since it comprises mathematical and computational developments of a complete wind energy conversion system (WECS) endowed with PMSG using time domain techniques of Alternative Transients Program (ATP), which prides itself a recognized reputation by scientific and academic communities as well as by electricity professionals in Brazil and elsewhere. The modeling procedures performed allowed the elaboration of blocks representing each of the elements of a real WECS, comprising the primary source (the wind), the wind turbine, the PMSG, the frequency converter, the step up transformer, the load composition and the power grid equivalent. Special attention is also given to the implementation of wind turbine control techniques, mainly the pitch control responsible for keeping the generator under the maximum power operation point, and the vector theory that aims at adjusting the active/reactive power flow between the wind turbine and the power grid. Several simulations are performed to investigate the dynamic behavior of the wind farm when subjected to different operating conditions and/or on the occurrence of wind intensity variations. The results have shown the effectiveness of both mathematical and computational modeling developed for the wind turbine and the associated controls.

Influence of the sputtering flow regime on the structural properties and magnetic behavior of Fe-Ga thin films (Ga ∼ 30 at.%)

In this paper we analyze the structure of Fe-Ga layers with a Ga content of ∼30 at.% deposited by the sputtering technique under two different regimes. We also studied the correlation between the structure and magnetic behavior of the samples. Keeping the Ar pressure fixed, we modified the flow regime from ballistic to diffusive by increasing the distance between the target and the substrate. X-ray diffraction measurements have shown a lower structural quality when growing in the diffusive flow. We investigated the impact of the growth regime by means of x-ray absorption fine structure (XAFS) measurements and obtained signs of its influence on the local atomic order. Full multiple scattering and finite difference calculations based on XAFS measurements point to a more relevant presence of a disordered A2 phase and of orthorhombic Ga clusters on the Fe-Ga alloy deposited under a diffusive regime; however, in the ballistic sample, a higher presence of D0_3/B2 phases is evidenced. Structural characteristics, from local to long range, seem to determine the magnetic behavior of the layers. Whereas a clear in-plane magnetic anisotropy is observed in the film deposited under ballistic flow, the diffusive sample is magnetically isotropic. Therefore, our experimental results provide evidence of a correlation between flow regime and structural properties and its impact on the magnetic behavior of a rather unexplored compositional region of Fe-Ga compounds.

Stochastically forced dislocation density distribution in plastic deformation

The dynamical evolution of dislocations in plastically deformed metals is controlled by both deterministic factors arising out of applied loads and stochastic effects appearing due to fluctuations of internal stress. Such type of stochastic dislocation processes and the associated spatially inhomogeneous modes lead to randomness in the observed deformation structure. Previous studies have analyzed the role of randomness in such textural evolution but none of these models have considered the impact of a finite decay time (all previous models assumed instantaneous relaxation which is "unphysical") of the stochastic perturbations in the overall dynamics of the system. The present article bridges this knowledge gap by introducing a colored noise in the form of an Ornstein-Uhlenbeck noise in the analysis of a class of linear and nonlinear Wiener and Ornstein-Uhlenbeck processes that these structural dislocation dynamics could be mapped on to. Based on an analysis of the relevant Fokker-Planck model, our results show that linear Wiener processes remain unaffected by the second time scale in the problem but all nonlinear processes, both Wiener type and Ornstein-Uhlenbeck type, scale as a function of the noise decay time τ. The results are expected to ramify existing experimental observations and inspire new numerical and laboratory tests to gain further insight into the competition between deterministic and random effects in modeling plastically deformed samples.

Assessment of primary energy conversions of oscillating water columns. I. Hydrodynamic analysis

This is an investigation on the development of a numerical assessment method for the hydrodynamic performance of an oscillating water column (OWC) wave energy converter. In the research work, a systematic study has been carried out on how the hydrodynamic problem can be solved and represented reliably, focusing on the phenomena of the interactions of the wave-structure and the wave-internal water surface. These phenomena are extensively examined numerically to show how the hydrodynamic parameters can be reliably obtained and used for the OWC performance assessment. In studying the dynamic system, a two-body system is used for the OWC wave energy converter. The first body is the device itself, and the second body is an imaginary “piston,” which replaces part of the water at the internal water surface in the water column. One advantage of the two-body system for an OWC wave energy converter is its physical representations, and therefore, the relevant mathematical expressions and the numerical simulation can be straightforward. That is, the main hydrodynamic parameters can be assessed using the boundary element method of the potential flow in frequency domain, and the relevant parameters are transformed directly from frequency domain to time domain for the two-body system. However, as it is shown in the research, an appropriate representation of the “imaginary” piston is very important, especially when the relevant parameters have to be transformed from frequency-domain to time domain for a further analysis. The examples given in the research have shown that the correct parameters transformed from frequency domain to time domain can be a vital factor for a successful numerical simulation.

Hydrodynamics of oscillating water column wave energy converters

This work deals with the numerical studies on hydrodynamics of oscillating water column (OWC) wave energy converters and its damping optimization on maximizing wave energy conversion by the OWC device. As a fundamental step, the hydrodynamic problems have been systematically studied by considering the interactions of the wave-structure and of the wave-internal water surface. Our first attention is on how the hydrodynamic performance can be reliably assessed, especially when it comes to the time-domain analysis, and what the physics behind the considerations is. Further on, a damping optimization for the OWC wave energy converter is also present based on the dynamics of the linear system, and a study on how we can optimize the damping for the given sea states so that the power conversion from irregular waves from irregular waves can be maximized.

Primary wave energy conversions of oscillating water columns

This paper presents a study on the numerical simulation of the primary wave energy conversion in the oscillating water column (OWC) wave energy converters (WECs). The new proposed numerical approach consists of three major components: potential flow analysis for the conventional hydrodynamic parameters, such as added mass, damping coefficients, restoring force coefficients and wave excitations; the thermodynamic analysis of the air in the air chamber, which is under the assumptions of the given power take-off characteristics and an isentropic process of air flow. In the formulation, the air compressibility and its effects have been included; and a time-domain analysis by combining the linear potential flow and the thermodynamics of the air flow in the chamber, in which the hydrodynamics and thermodynamics/aerodynamics have been coupled together by the force generated by the pressurised and de-pressurised air in the air chamber, which in turn has effects on the motions of the structure and the internal water surface. As an example, the new developed approach has been applied to a fixed OWC device. The comparisons of the measured data and the simulation results show the new method is very capable of predicting the performance of the OWC devices.

Wave energy conversion of oscillating water column devices including air compressibility

This paper presents an investigation on air compressibility in the air chamber and its effects on the power conversion of oscillating water column (OWC) devices. As it is well known that for practical OWC plants, their air chambers may be large enough for accommodating significant air compressibility, the “spring effect,” an effect that is frequently and simply regarded to store and release energy during the reciprocating process of a wave cycle. Its insight effects on the device’s performance and power conversion, however, have not been studied in detail. This research will investigate the phenomena with a special focus on the effects of air compressibility on wave energy conversion. Air compressibility itself is a complicated nonlinear process in nature, but it can be linearised for numerical simulations under certain assumptions for frequency domain analysis. In this research work, air compressibility in the OWC devices is first linearised and further coupled with the hydrodynamics of the OWC. It is able to show mathematically that in frequency-domain, air compressibility can increase the spring coefficients of both the water body motion and the device motion (if it is a floating device), and enhance the coupling effects between the water body and the structure. Corresponding to these changes, the OWC performance, the capture power, and the optimised Power Take-off (PTO) damping coefficient in the wave energy conversion can be all modified due to air compressibility. To validate the frequency-domain results and understand the problems better, the more accurate time-domain simulations with fewer assumptions have been used for comparison. It is shown that air compressibility may significantly change the dynamic responses and the capacity of converting wave energy of the OWC devices if the air chamber is very large.

Dynamic effects of anchor positional tolerance on tension moored floating wind turbine

For water depths greater than 60m floating wind turbines will become the most economical option for generating offshore wind energy. Tension mooring stabilised units are one type of platform being considered by the offshore wind energy industry. The complex mooring arrangement used by this type of platform means that the dynamics are greatly effected by offsets in the positioning of the anchors. This paper examines the issue of tendon anchor position tolerances. The dynamic effects of three positional tolerances are analysed in survival state using the time domain FASTLink. The severe impact of worst case anchor positional offsets on platform and turbine survivability is shown. The worst anchor misposition combinations are highlighted and should be strongly avoided. Novel methods to mitigate this issue are presented.

Particle positions of numerical experiments 1 to 4

In mixed sediment beds, erosion resistance can change relative to that of beds composed of a uniform sediment because of varying textural and/or other grain-size parameters, with effects on pore water flow that are difficult to quantify by means of analogue techniques. To overcome this difficulty, a three-dimensional numerical model was developed using a finite difference method (FDM) flow model coupled with a distinct element method (DEM) particle model. The main aim was to investigate, at a high spatial resolution, the physical processes occurring during the initiation of motion of single grains at the sediment-water interface and in the shallow subsurface of simplified sediment beds under different flow velocities. Increasing proportions of very fine sand (D50=0.08 mm) were mixed into a coarse sand matrix (D50=0.6 mm) to simulate mixed sediment beds, starting with a pure coarse sand bed in experiment 1 (0 wt% fines), and proceeding through experiment 2 (6.5 wt% fines), experiment 3 (10.5 wt% fines), and experiment 4 (28.7 wt% fines). All mixed beds were tested for their erosion behavior at predefined flow velocities varying in the range of U 1-5=10-30 cm/s. The experiments show that, with increasing fine content, the smaller particles increasingly fill the spaces between the larger particles. As a consequence, pore water inflow into the sediment is increasingly blocked, i.e., there is a decrease in pore water flow velocity and, hence, in the flow momentum available to entrain particles. These findings are portrayed in a new conceptual model of enhanced sediment bed stabilization.

Stress Inversion Using LiDAR Tunnel Deformation Measurements

Far-field stresses are those present in a volume of rock prior to excavations being created. Estimates of the orientation and magnitude of far-field stresses, often used in mine design, are generally obtained by single-point measurements of stress, or large-scale, regional trends. Point measurements can be a poor representation of far-field stresses as a result of excavation-induced stresses and geological structures. For these reasons, far-field stress estimates can be associated with high levels of uncertainty. The purpose of this thesis is to investigate the practical feasibility, applications, and limitations of calibrating far-field stress estimates through tunnel deformation measurements captured using LiDAR imaging. A method that estimates the orientation and magnitude of excavation-induced principal stress changes through back-analysis of deformation measurements from LiDAR imaged tunnels was developed and tested using synthetic data. If excavation-induced stress change orientations and magnitudes can be accurately estimated, they can be used in the calibration of far-field stress input to numerical models. LiDAR point clouds have been proven to have a number of underground applications, thus it is desired to explore their use in numerical model calibration. The back-analysis method is founded on the superposition of stresses and requires a two-dimensional numerical model of the deforming tunnel. Principal stress changes of known orientation and magnitude are applied to the model to create calibration curves. Estimation can then be performed by minimizing squared differences between the measured tunnel and sets of calibration curve deformations. In addition to the back-analysis estimation method, a procedure consisting of previously existing techniques to measure tunnel deformation using LiDAR imaging was documented. Under ideal conditions, the back-analysis method estimated principal stress change orientations within ±5° and magnitudes within ±2 MPa. Results were comparable for four different tunnel profile shapes. Preliminary testing using plastic deformation, a rough tunnel profile, and profile occlusions suggests that the method can work under more realistic conditions. The results from this thesis set the groundwork for the continued development of a new, inexpensive, and efficient far-field stress estimate calibration method.

Linking Parametric CAD with Adjoint Surface Sensitivities

The goal of this work is to present an efficient CAD-based adjoint process chain for calculating parametric sensitivities (derivatives of the objective function with respect to the CAD parameters) in timescales acceptable for industrial design processes. The idea is based on linking parametric design velocities (geometric sensitivities computed from the CAD model) with adjoint surface sensitivities. A CAD-based design velocity computation method has been implemented based on distances between discrete representations of perturbed geometries. This approach differs from other methods due to the fact that it works with existing commercial CAD packages (unlike most analytical approaches) and it can cope with the changes in CAD model topology and face labeling. Use of the proposed method allows computation of parametric sensitivities using adjoint data at a computational cost which scales with the number of objective functions being considered, while it is essentially independent of the number of design variables. The gradient computation is demonstrated on test cases for a Nozzle Guide Vane (NGV) model and a Turbine Rotor Blade model. The results are validated against finite difference values and good agreement is shown. This gradient information can be passed to an optimization algorithm, which will use it to update the CAD model parameters.

Modeling, Analysis and Performance Comparison of Two Direct Sampling DCSK Receivers under Frequency Nonselective Fading Channels

Two direct sampling correlator-type receivers for differential chaos shift keying (DCSK) communication systems under frequency non-selective fading channels are proposed. These receivers operate based on the same hardware platform with different architectures. In the first scheme, namely sum-delay-sum (SDS) receiver, the sum of all samples in a chip period is correlated with its delayed version. The correlation value obtained in each bit period is then compared with a fixed threshold to decide the binary value of recovered bit at the output. On the other hand, the second scheme, namely delay-sum-sum (DSS) receiver, calculates the correlation value of all samples with its delayed version in a chip period. The sum of correlation values in each bit period is then compared with the threshold to recover the data. The conventional DCSK transmitter, frequency non-selective Rayleigh fading channel, and two proposed receivers are mathematically modelled in discrete-time domain. The authors evaluated the bit error rate performance of the receivers by means of both theoretical analysis and numerical simulation. The performance comparison shows that the two proposed receivers can perform well under the studied channel, where the performances get better when the number of paths increases and the DSS receiver outperforms the SDS one.

A Bias Compensated Cross-Relation approach to Thermocouple Characterisation

The measurement of fast changing temperature fluctuations is a challenging problem due to the inherent limited bandwidth of temperature sensors. This results in a measured signal that is a lagged and attenuated version of the input. Compensation can be performed provided an accurate, parameterised sensor model is available. However, to account for the in influence of the measurement environment and changing conditions such as gas velocity, the model must be estimated in-situ. The cross-relation method of blind deconvolution is one approach for in-situ characterisation of sensors. However, a drawback with the method is that it becomes positively biased and unstable at high noise levels. In this paper, the cross-relation method is cast in the discrete-time domain and a bias compensation approach is developed. It is shown that the proposed compensation scheme is robust and yields unbiased estimates with lower estimation variance than the uncompensated version. All results are verified using Monte-Carlo simulations.