Time-domain simulation of wave-induced ship motions by a Rankine panel method

In this thesis, a numerical program has been developed to simulate the wave-induced ship motions in the time domain. Wave-body interactions have been studied for various ships and floating bodies through forced motion and free motion simulations in a wide range of wave frequencies. A three-dimensional Rankine panel method is applied to solve the boundary value problem for the wave-body interactions. The velocity potentials and normal velocities on the boundaries are obtained in the time domain by solving the mixed boundary integral equations in relation to the source and dipole distributions. The hydrodynamic forces are calculated by the integration of the instantaneous hydrodynamic pressures over the body surface. The equations of ship motion are solved simultaneously with the boundary value problem for each time step. The wave elevation is computed by applying the linear free surface conditions. A numerical damping zone is adopted to absorb the outgoing waves in order to satisfy the radiation condition for the truncated free surface. A numerical filter is applied on the free surface for the smoothing of the wave elevation. Good convergence has been reached for both forced motion simulations and free motion simulations. The computed added-mass and damping coefficients, wave exciting forces, and motion responses for ships and floating bodies are in good agreement with the numerical results from other programs and experimental data.

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.

Beryllium isotopes of sediments from the Portuguese margin

A high resolution study of authigenic Be isotopes (10Be and 9Be) combined with continuous relative paleointensity records has been performed along the same marine sedimentary sequences from the Portuguese margin (NE Atlantic) covering the past 300 kyr in order to assess relationships between geomagnetic moment variations and 10Be production rate variations. A careful examination of the various ways of taking into account environmental disturbing effects on the authigenic 10Be concentration leads to the conclusion that the most reliable proxy of cosmonuclide production rates is presently the authigenic 10Be/9Be ratio. Eight intervals of significant authigenic 10Be/9Be enhancement evidence geomagnetic moment drops related to global paleomagnetic excursions, some being already admitted, others being proposed as new geomagnetic features. Since, unlike sedimentary magnetic remanence, the authigenic 10Be/9Be records dipole moment variations without significant acquisition delay, it provides better constraints on their timing. Comparison of 10Be/9Be and benthic delta18O records from the same cores suggests that dipole moment lows preferentially occurred during or at the end of interglacial episodes, with a quasi-period of 100 kyr.

Resonance fluorescence spectrum of a \Lambda-type quantum emitter close to a metallic nanoparticle

We theoretically study the resonance fluorescence spectrum of a three-level quantum emitter coupled to a spherical metallic nanoparticle. We consider the case that the quantum emitter is driven by a single laser field along one of the optical transitions. We show that the development of the spectrum depends on the relative orientation of the dipole moments of the optical transitions in relation to the metal nanoparticle. In addition, we demonstrate that the location and width of the peaks in the spectrum are strongly modified by the exciton-plasmon coupling and the laser detuning, allowing to achieve controlled strongly subnatural spectral line. A strong antibunching of the fluorescent photons along the undriven transition is also obtained. Our results may be used for creating a tunable source of photons which could be used for a probabilistic entanglement scheme in the field of quantum information processing.

Localising the auditory N1m with event-related beamformers:localisation accuracy following bilateral and unilateral stimulation

The auditory evoked N1m-P2m response complex presents a challenging case for MEG source-modelling, because symmetrical, phase-locked activity occurs in the hemispheres both contralateral and ipsilateral to stimulation. Beamformer methods, in particular, can be susceptible to localisation bias and spurious sources under these conditions. This study explored the accuracy and efficiency of event-related beamformer source models for auditory MEG data under typical experimental conditions: monaural and diotic stimulation; and whole-head beamformer analysis compared to a half-head analysis using only sensors from the hemisphere contralateral to stimulation. Event-related beamformer localisations were also compared with more traditional single-dipole models. At the group level, the event-related beamformer performed equally well as the single-dipole models in terms of accuracy for both the N1m and the P2m, and in terms of efficiency (number of successful source models) for the N1m. The results yielded by the half-head analysis did not differ significantly from those produced by the traditional whole-head analysis. Any localisation bias caused by the presence of correlated sources is minimal in the context of the inter-individual variability in source localisations. In conclusion, event-related beamformers provide a useful alternative to equivalent-current dipole models in localisation of auditory evoked responses.

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.

The application of aryne chemistry and use of α-diazo carbonyl derivatives in indazole synthesis with biological evaluation

This thesis outlines the synthetic chemistry involved in the preparation of a range of novel indazole compounds and details the subsequent investigation into their potential as biologically active agents. The synthetic route utilised in this research to form the indazole structure was the [3+2] dipolar cycloaddition of diazo carbonyl compounds with reactive aryne intermediates generated in situ. The preparation of further novel indazole derivatives containing different functional groups and substituents was performed by synthesising alternative 1,3- dipole and dipolarophile analogues and provided additionally diverse compounds. Further derivatisation of the indazole product was made possible by deacylation and alkylation methods. Transformation reactions were performed on alkenecontaining ester side chains to provide novel epoxide, aldehyde and tertiary amine derivatives. The first chapter is a review of the literature beginning with a short overview on the structure, reactivity and common synthetic routes to diazo carbonyl derivatives. More attention is given to the use of diazo compounds as 1,3-dipoles in cycloaddition reactions or where the diazo group is incorporated into the final product. A review of the interesting background, structure and reactivity of aryne intermediates is also presented. In addition, some common syntheses of indazole compounds are presented as well as a brief discussion on the importance of indazole compounds as therapeutic agents. The second chapter discusses the synthetic routes employed towards the synthesis of the range of indazoles. Initially, the syntheses of the diazo carbonyl and aryne precursors are described. Next, the synthetic methods to prepare the indazole compounds are provided followed by discussion on derivatisation of the indazole compounds including N-deacylation, N-benzylation and ester side-chain transformation of some alkene-containing indazoles. A series of novel indazole derivatives were submitted for anti-cancer screening at the U.S National Cancer Institute (NCI). A number of these derivatives were identified as hit compounds, with excellent growth inhibition. The results obtained from biological evaluation from the NCI are provided with further results pending from the Community for Open Antimicrobial Drug Discovery. The third chapter details the full experimental procedures, including spectroscopic and analytical data for all the compounds prepared during this research.

Radiocarbon age, Mg/Ca and d18O measurements on planktonic foraminifera of sediment core GeoB12605-3

The sea surface temperature (SST) of the tropical Indian Ocean is a major component of global climate teleconnections. While the Holocene SST history is documented for regions affected by the Indian and Arabian monsoons, data from the near-equatorial western Indian Ocean are sparse. Reconstructing past zonal and meridional SST gradients requires additional information on past temperatures from the western boundary current region. We present a unique record of Holocene SST and thermocline depth variations in the tropical western Indian Ocean as documented in foraminiferal Mg/Ca ratios and d18O from a sediment core off northern Tanzania. For Mg/Ca and thermocline d18O, most variance is concentrated in the centennial to bicentennial periodicity band. On the millennial time scale, an early to mid-Holocene (~7.8-5.6 ka) warm phase is followed by a temperature drop by up to 2°C, leading to a mid-Holocene cool interval (5.6-4.2 ka). The shift is accompanied by an initial reduction in the difference between surface and thermocline foraminiferal d18O, consistent with the thickening of the mixed layer and suggestions of a strengthened Walker circulation. However, we cannot confirm the expected enhanced zonal SST gradient, as the cooling of similar magnitude had previously been found in SSTs from the upwelling region off Sumatra and in Flores air temperatures. The SST pattern probably reflects the tropical Indian Ocean expression of a large-scale climate anomaly rather than a positive Indian Ocean Dipole-like mean state.

Planktonic foraminiferal flux and shell geochemistry data in the sediment trap time series JAM 1-3 off south Java

Results from sediment trap experiments conducted in the seasonal upwelling area off south Java from November 2000 until July 2003 revealed significant monsoon-, El Niño-Southern Oscillation-, and Indian Ocean Dipole-induced seasonal and interannual variations in flux and shell geochemistry of planktonic foraminifera. Surface net primary production rates together with total and species-specific planktonic foraminiferal flux rates were highest during the SE monsoon-induced coastal upwelling period from July to October, with three species Globigerina bulloides, Neogloboquadrina pachyderma dex., and Globigerinita glutinata contributing to 40% of the total foraminiferal flux. Shell stable oxygen isotopes (d18O) and Mg/Ca data of Globigerinoides ruber sensu stricto (s.s.), G. ruber sensu lato (s.l.), Neogloboquadrina dutertrei, Pulleniatina obliquiloculata, and Globorotalia menardii in the sediment trap time series recorded surface and subsurface conditions. We infer habitats of 0-30 m for G. ruber at the mixed layer depth, 60-80 m (60-90 m) for P. obliquiloculata (N. dutertrei) at the upper thermocline depth, and 90-110 m (100-150 m) for G. menardii in the 355-500 mm (>500 µm) size fraction corresponding to the (lower) thermocline depth in the study area. Shell Mg/Ca ratio of G. ruber (s.l. and s.s.) reveals an exponential relationship with temperature that agrees with published relationships particularly with the Anand et al. (2003) equations. Flux-weighted foraminiferal data in sediment trap are consistent with average values in surface sediment samples off SW Indonesia. This consistency confirms the excellent potential of these proxies for reconstructing past environmental conditions in this part of the ocean realm.

(Table 2) U-channel paleomagnetic record of ODP Site 178-1098

Thick Holocene sedimentary sections (>45 m) cored in the Palmer Deep by the United States Antarctic Program (USAP) and during Ocean Drilling Program (ODP) Leg 178 provide the first opportunity to examine past geomagnetic field behavior at high southern latitudes. After removal of a low-coercivity drilling overprint the sediments display a stable, single-component remanent magnetization. Two short cores that recovered the uppermost 2.6 m of sediment have inclinations that fluctuate about the present day inclination (-57°) measured at Faraday Station, and several features with wavelengths of 10 to 20 cm appear to be correlative. However, shipboard measurements of inclination fluctuations on split-core samples from three holes drilled at ODP Site 1098 do not correlate well with each other, even though the intensity and susceptibility data correlate very well and the overall mean inclination for cores from each hole is consistent with the expected geocentric axial dipole (GAD) inclination. The correlation is improved dramatically by using inclinations measured on u-channels taken from the pristine center of a split core. Consequently, the anomalous directions and the resulting poor between-hole correlation of inclinations obtained from shipboard data can be attributed to coring-induced deformation, which is common on the outer edge of ODP piston cores, and/or measurement artifacts in the split-core data. Our preferred inclination record is thus derived from u-channel results. The upper ~25 m represents continuous sedimentation over the past 9000 yr, with an average sedimentation rate exceeding 250 cm/kyr (0.25 cm/yr). Given that remanence measurements on u-channels average over an interval <7 cm long, we obtained independent measurements of the paleo-geomagnetic field that average over only ~30 yr. This high-resolution record is characterized by an inclination that fluctuates within +/-15° of the current GAD inclination.

Investigating the potential magnetic origin of wind variability in OB stars

In this thesis, the origin of large-scale structures in hot star winds, believed to be responsible for the presence of discrete absorption components (DACs) in the absorption troughs of ultraviolet resonance lines, is constrained using both observations and numerical simulations. These structures are understood as arising from bright regions on the stellar surface, although their physical cause remains unknown. First, we use high quality circular spectropolarimetric observations of 13 well-studied OB stars to evaluate the potential role of dipolar magnetic fields in producing DACs. We perform longitudinal field measurements and place limits on the field strength using Bayesian inference, assuming that it is dipolar. No magnetic field was detected within this sample. The derived constraints statistically refute any significant dynamical influence from a magnetic dipole on the wind for all of these stars, ruling out such fields as a cause for DACs. Second, we perform numerical simulations using bright spots constrained by broadband optical photometric observations. We calculate hydrodynamical wind models using three sets of spot sizes and strengths. Co-rotating interaction regions are yielded in each model, and radiative transfer shows that the properties of the variations in the UV resonance lines synthesized from these models are consistent with those found in observed UV spectra, establishing the first consistent link between UV spectroscopic line profile variability and photometric variations and thus supporting the bright spot paradigm (BSP). Finally, we develop and apply a phenomenological model to quantify the measurable effects co-rotating bright spots would have on broadband optical photometry and on the profiles of photopheric lines in optical spectra. This model can be used to evaluate the existence of these spots, and, in the event of their detection, characterize them. Furthermore, a tentative spot evolution model is presented. A preliminary analysis of its output, compared to the observed photometric variations of xi Persei, suggests the possible existence of “active longitudes” on the surface of this star. Future work will expand the range of observational diagnostics that can be interpreted within the BSP, and link phenomenology (bright spots) to physical processes (magnetic spots or non-radial pulsations).

The Rotational Evolution and Magnetospheric Emission of the Magnetic Early B-type Stars

How do the magnetic fields of massive stars evolve over time? Are their gyrochronological ages consistent with ages inferred from evolutionary tracks? Why do most stars predicted to host Centrifugal Magnetospheres (CMs) display no H$\alpha$ emission? Does plasma escape from CMs via centrifugal breakout events, or by a steady-state leakage mechanism? This thesis investigates these questions via a population study with a sample of 51 magnetic early B-type stars. The longitudinal magnetic field \bz~was measured from Least Squares Deconvolution profiles extracted from high-resolution spectropolarimetric data. New rotational periods $P_{\rm rot}$ were determined for 15 stars from \bz, leaving only 3 stars for which $P_{\rm rot}$ is unknown. Projected rotational velocities \vsini~were measured from multiple spectral lines. Effective temperatures and surface gravities were measured via ionization balances and line profile fitting of H Balmer lines. Fundamental physical parameters, \bz, \vsini, and $P_{\rm rot}$ were then used to determine radii, masses, ages, dipole oblique rotator model, stellar wind, magnetospheric, and spindown parameters using a Monte Carlo approach that self-consistently calculates all parameters while accounting for all available constraints on stellar properties. Dipole magnetic field strengths $B_{\rm d}$ follow a log-normal distribution similar to that of Ap stars, and decline over time in a fashion consistent with the expected conservation of fossil magnetic flux. $P_{\rm rot}$ increases with fractional main sequence age, mass, and $B_{\rm d}$, as expected from magnetospheric braking. However, comparison of evolutionary track ages to maximum spindown ages $t_{\rm S,max}$ shows that initial rotation fractions may be far below critical for stars with $M_*>10 M_\odot$. Computing $t_{\rm S,max}$ with different mass-loss prescriptions indicates that the mass-loss rates of B-type stars are likely much lower than expected from extrapolation from O-type stars. Stars with H$\alpha$ in emission and absorption occupy distinct regions in the updated rotation-magnetic confinement diagram: H$\alpha$-bright stars are found to be younger, more rapidly rotating, and more strongly magnetized than the general population. Emission strength is sensitive both to the volume of the CM and to the mass-loss rate, favouring leakage over centrifugal breakout.

Optical properties of chiral thin films and microparticles

In this work I study the optical properties of helical particles and chiral sculptured thin films, using computational modeling (discrete dipole approximation, Berreman calculus), and experimental techniques (glancing angle deposition, ellipsometry, scatterometry, and non-linear optical measurements). The first part of this work focuses on linear optics, namely light scattering from helical microparticles. I study the influence of structural parameters and orientation on the optical properties of particles: circular dichroism (CD) and optical rotation (OR), and show that as a consequence of random orientation, CD and OR can have the opposite sign, compared to that of the oriented particle, potentially resulting in ambiguity of measurement interpretation. Additionally, particles in random orientation scatter light with circular and elliptical polarization states, which implies that in order to study multiple scattering from randomly oriented chiral particles, the polarization state of light cannot be disregarded. To perform experiments and attempt to produce particles, a newly constructed multi stage thin film coating chamber is calibrated. It enables the simultaneous fabrication of multiple sculptured thin film coatings, each with different structure. With it I successfully produce helical thin film coatings with Ti and TiO_{2}. The second part of this work focuses on non-linear optics, with special emphasis on second-harmonic generation. The scientific literature shows extensive experimental and theoretical work on second harmonic generation from chiral thin films. Such films are expected to always show this non-linear effect, due to their lack of inversion symmetry. However no experimental studies report non-linear response of chiral sculptured thin films. In this work I grow films suitable for a second harmonic generation experiment, and report the first measurements of non-linear response.

van der Waals coefficients for positronium interactions with atoms

The random-phase approximation with exchange (RPAE) is used with a B-spline basis to compute dynamic dipole polarizabilities of noble-gas atoms and several other closed-shell atoms (Be, Mg, Ca, Zn, Sr, Cd, and Ba). From these, values of the van der Waals C₆ constants for positronium interactions with these atoms are determined and compared with existing data. After correcting the RPAE polarizabilities to fit the most accurate static polarizability data, our best predictions of C₆ for Ps–noble-gas pairs are expected to be accurate to within 1%, and to within a few percent for the alkaline-earth metals. We also used accurate dynamic dipole polarizabilities from the literature to compute the C₆ coefficients for the alkali-metal atoms. Implications of increased C₆ values for Ps scattering from more polarizable atoms are discussed.

Ultracompact Retrodirective Antenna Arrays With Superdirective Radiation Patterns

It is shown that the direction-of-arrival (DoA) information carried by an incident electromagnetic (EM) wave can be encoded into the evanescent near field of an electrically small resonance antenna array with a spatial rate higher than that of the incident field oscillation rate in free space. Phase conjugation of the received signal leads to the retrodirection of the near field in the antenna array environment, which in turn generates a retrodirected far-field beam toward the original DoA. This EM phenomenon enables electrically small retrodirective antenna arrays with superdirective, angular super-resolution, auto-pointing properties for an arbitrary DoA. A theoretical explanation of the phenomenon based on first principal observations is given and full-wave simulations demonstrate a realizability route for the proposed retrodirective terminal that is comprised of resonance dipole antenna elements. Specifically, it is shown that a three-element disk-loaded retrodirective dipole array with 0.15\lambda spacings can achieve a 3.4-dBi maximal gain, 3-dBi front-to-back ratio, and 13% return loss fractional bandwidth (at the 10-dB level). Then, it is demonstrated that the radiation gain of a three-element array can be improved to approximately 6 dBi at the expense of the return loss fractional bandwidth reduction (2%).