Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

A novel fractional sampling offset estimation in an OFDM system

The Orthogonal Frequency Division Multiplexing (OFDM) is a form of Multi-Carrier Modulation where the data stream is transmitted over a number of carriers which are orthogonal to each other i.e. the carrier spacing is selected such that each carrier is located at the zeroes of all other carriers in the spectral domain. This paper proposes a new novel sampling offset estimation algorithm for an OFDM system in order to receive the OFDM data symbols error-free over the noisy channel at the receiver and to achieve fine timing synchronization between the transmitter and the receiver. The performance of this algorithm has been studied in AWGN, ADSL and SUI channels successfully.

Volume change behaviour of fine grained solis

Soils in and and semi-arid zones undergoes volume changes due to wetting. Depending upon the type of clay minerals present, degree of saturation, externally applied load and bonding, the fine grained soils either swells or compresses. One of the parameter that affects the volume change behaviour is the primary clay mineral present in their clay size fraction. A simple method of identifying the same has been presented. It has been brought out that in an expansive unsaturated undisturbed soil, the diffuse double layer repulsion, the stress state and the bonding play significant role in their volume change behaviour. In non-expansive fine grained unsaturated undisturbed soils, the shearing resistance at particle level (including the matrix suction and bonding) and fabric play a significant role in influencing the volume change behaviour. While both the mechanism co-exist, one of them play a dominant role depending upon the primary clay mineral is swelling or non swelling.

A novel fractional sampling offset estimation in an OFDM system

The Orthogonal Frequency Division Multiplexing (OFDM) is a form of Multi-Carrier Modulation where the data stream is transmitted over a number of carriers which are orthogonal to each other i.e. the carrier spacing is selected such that each carrier is located at the zeroes of all other carriers in the spectral domain. This paper proposes a new novel sampling offset estimation algorithm for an OFDM system in order to receive the OFDM data symbols error-free over the noisy channel at the receiver and to achieve fine timing synchronization between the transmitter and the receiver. The performance of this algorithm has been studied in AWGN, ADSL and SUI channels successfully.

Influence of Oxide Particle Network Morphology on Ion Solvation and Transport in ``Soggy Sand''Electrolytes

The role of oxide surface chemical composition and solvent on ion solvation and ion transport of ``soggy sand'' electrolytes are discussed here. A ``soggy sand'' electrolyte system comprising dispersions of hydrophilic/hydrophobic functionalized aerosil silica in lithium perchlorate methoxy polyethylene glycol solution was employed for the study. Static and dynamic rheology measurements show formation of an attractive particle network in the case of the composite with unmodified aerosil silica (i.e., with surface silanol groups) as well as composites with hydrophobic alkane groups. While particle network in the composite with hydrophilic aerosil silica (unmodified) were due to hydrogen bonding, hydrophobic aerosil silica particles were held together via van der Waals forces. The network strength in the latter case (i.e., for hydrophobic composites) were weaker compared with the composite with unmodified aerosil silica. Both unmodified silica as well as hydrophobic silica composites displayed solid-like mechanical strength. No enhancement in ionic conductivity compared to the liquid electrolyte was observed in the case of the unmodified silica. This was attributed to the existence of a very strong particle network, which led to the ``expulsion'' of all conducting entities from the interfacial region between adjacent particles. The ionic conductivity for composites with hydrophobic aerosil particles displayed ionic conductivity dependent on the size of the hydrophobic chemical moiety. No spanning attractive particle network was observed for aerosil particles with surfaces modified with stronger hydrophilic groups (than silanol). The composite resembled a sol, and no percolation in ionic conductivity was observed.

Dual-Drug Delivery System Based on In Situ Gel-Forming Nanosuspension of Forskolin to Enhance Antiglaucoma Efficacy

The present study was designed to improve the bioavailability of forskolin by the influence of precorneal residence time and dissolution characteristics. Nanosizing is an advanced approach to overcome the issue of poor aqueous solubility of active pharmaceutical ingredients. Forskolin nanocrystals have been successfully manufactured and stabilized by poloxamer 407. These nanocrystals have been characterized in terms of particle size by scanning electron microscopy and dynamic light scattering. By formulating Noveon AA-1 polycarbophil/poloxamer 407 platforms, at specific concentrations, it was possible to obtain a pH and thermoreversible gel with a pH(gel)/T-gel close to eye pH/temperature. The addition of forskolin nanocrystals did not alter the gelation properties of Noveon AA-1 polycarbophil/poloxamer 407 and nanocrystal properties of forskolin. The formulation was stable over a period of 6 months at room temperature. In vitro release experiments indicated that the optimized platform was able to prolong and control forskolin release for more than 5 h. The in vivo studies on dexamethasone-induced glaucomatous rabbits indicated that the intraocular pressure lowering efficacy for nanosuspension/hydrogel systems was 31% and lasted for 12 h, which is significantly better than the effect of traditional eye suspension (18%, 4-6 h). Hence, our investigations successfully prove that the pH and thermoreversible polymeric in situ gel-forming nanosuspension with ability of controlled drug release exhibits a greater potential for glaucoma therapy.

Direct detection of atmospheric particle formation using the Neutral cluster and Air Ion Spectrometer

Aerosol particles play an important role in the Earth s atmosphere and in the climate system: they scatter and absorb solar radiation, facilitate chemical processes, and serve as seeds for cloud formation. Secondary new particle formation (NPF) is a globally important source of these particles. Currently, the mechanisms of particle formation and the vapors participating in this process are, however, not truly understood. In order to fully explain atmospheric NPF and subsequent growth, we need to measure directly the very initial steps of the formation processes. This thesis investigates the possibility to study atmospheric particle formation using a recently developed Neutral cluster and Air Ion Spectrometer (NAIS). First, the NAIS was calibrated and intercompared, and found to be in good agreement with the reference instruments both in the laboratory and in the field. It was concluded that NAIS can be reliably used to measure small atmospheric ions and particles directly at the sizes where NPF begins. Second, several NAIS systems were deployed simultaneously at 12 European measurement sites to quantify the spatial and temporal distribution of particle formation events. The sites represented a variety of geographical and atmospheric conditions. The NPF events were detected using NAIS systems at all of the sites during the year-long measurement period. Various particle formation characteristics, such as formation and growth rates, were used as indicators of the relevant processes and participating compounds in the initial formation. In a case of parallel ion and neutral cluster measurements, we also estimated the relative contribution of ion-induced and neutral nucleation to the total particle formation. At most sites, the particle growth rate increased with the increasing particle size indicating that different condensing vapors are participating in the growth of different-sized particles. The results suggest that, in addition to sulfuric acid, organic vapors contribute to the initial steps of NPF and to the subsequent growth, not just later steps of the particle growth. As a significant new result, we found out that the total particle formation rate varied much more between the different sites than the formation rate of charged particles. The results infer that the ion-induced nucleation has a minor contribution to particle formation in the boundary layer in most of the environments. These results give tools to better quantify the aerosol source provided by secondary NPF in various environments. The particle formation characteristics determined in this thesis can be used in global models to assess NPF s climatic effects.

Position-sensitive silicon strip detector characterization using particle beams

Silicon strip detectors are fast, cost-effective and have an excellent spatial resolution. They are widely used in many high-energy physics experiments. Modern high energy physics experiments impose harsh operation conditions on the detectors, e.g., of LHC experiments. The high radiation doses cause the detectors to eventually fail as a result of excessive radiation damage. This has led to a need to study radiation tolerance using various techniques. At the same time, a need to operate sensors approaching the end their lifetimes has arisen. The goal of this work is to demonstrate that novel detectors can survive the environment that is foreseen for future high-energy physics experiments. To reach this goal, measurement apparatuses are built. The devices are then used to measure the properties of irradiated detectors. The measurement data are analyzed, and conclusions are drawn. Three measurement apparatuses built as a part of this work are described: two telescopes measuring the tracks of the beam of a particle accelerator and one telescope measuring the tracks of cosmic particles. The telescopes comprise layers of reference detectors providing the reference track, slots for the devices under test, the supporting mechanics, electronics, software, and the trigger system. All three devices work. The differences between these devices are discussed. The reconstruction of the reference tracks and analysis of the device under test are presented. Traditionally, silicon detectors have produced a very clear response to the particles being measured. In the case of detectors nearing the end of their lifefimes, this is no longer true. A new method benefitting from the reference tracks to form clusters is presented. The method provides less biased results compared to the traditional analysis, especially when studying the response of heavily irradiated detectors. Means to avoid false results in demonstrating the particle-finding capabilities of a detector are also discussed. The devices and analysis methods are primarily used to study strip detectors made of Magnetic Czochralski silicon. The detectors studied were irradiated to various fluences prior to measurement. The results show that Magnetic Czochralski silicon has a good radiation tolerance and is suitable for future high-energy physics experiments.

Root system traits of Norway spruce, Scots pine, and silver birch in mixed boreal forests: an analysis of root architecture, morphology, and anatomy

The aim of this thesis was to unravel the functional-structural characteristics of root systems of Betula pendula Roth., Picea abies (L.) Karst., and Pinus sylvestris L. in mixed boreal forest stands differing in their developmental stage and site fertility. The root systems of these species had similar structural regularities: horizontally-oriented shallow roots defined the horizontal area of influence, and within this area, each species placed fine roots in the uppermost soil layers, while sinker roots defined the maximum rooting depth. Large radial spread and high ramification of coarse roots, and the high specific root length (SRL) and root length density (RLD) of fine roots indicated the high belowground competitiveness and root plasticity of B. pendula. Smaller radial root spread and sparser branching of coarse roots, and low SRL and RLD of fine roots of the conifers could indicate their more conservative resource use and high association with and dependence on ectomycorrhiza-forming fungi. The vertical fine root distributions of the species were mostly overlapping, implying the possibility for intense belowground competition for nutrients. In each species, conduits tapered and their frequency increased from distal roots to the stem, from the stem to the branches, and to leaf petioles in B. pendula. Conduit tapering was organ-specific in each species violating the assumptions of the general vascular scaling model (WBE). This reflects the hierarchical organization of a tree and differences between organs in the relative importance of transport, safety, and mechanical demands. The applied root model was capable of depicting the mass, length and spread of coarse roots of B. pendula and P. abies, and to the lesser extent in P. sylvestris. The roots did not follow self-similar fractal branching, because the parameter values varied within the root systems. Model parameters indicate differences in rooting behavior, and therefore different ecophysiological adaptations between species.

Polarization studies in electromagnetic scattering by small Solar-system particles

In remote-sensing studies, particles that are comparable to the wavelength exhibit characteristic features in electromagnetic scattering, especially in the degree of linear polarization. These features vary with the physical properties of the particles, such as shape, size, refractive index, and orientation. In the thesis, the direct problem of computing the unknown scattered quantities using the known properties of the particles and the incident radiation is solved at both optical and radar spectral regions in a unique way. The internal electromagnetic fields of wavelength-scale particles are analyzed by using both novel and established methods to show how the internal fields are related to the scattered fields in the far zone. This is achieved by using the tools and methods that were developed specifically to reveal the internal field structure of particles and to study the mechanisms that relate the structure to the scattering characteristics of those particles. It is shown that, for spherical particles, the internal field is a combination of a forward propagating wave with the apparent wavelength determined by the refractive index of the particle, and a standing wave pattern with the apparent wavelength the same as for the incident wave. Due to the surface curvature and dielectric nature of the particle, the incident wave front undergoes a phase shift, and the resulting internal wave is focused mostly at the forward part of the particle similar to an optical lens. This focusing is also seen for irregular particles. It is concluded that, for both spherical and nonspherical particles, the interference at the far field between the partial waves that originate from these concentrated areas in the particle interior, is responsible for the specific polarization features that are common for wavelength-scale particles, such as negative values and local extrema in the degree of linear polarization, asymmetry of the phase function, and enhancement of intensity near the backscattering direction. The papers presented in this thesis solve the direct problem for particles with both simple and irregular shapes to demonstrate that these interference mechanisms are common for all dielectric wavelength-scale particles. Furthermore, it is shown that these mechanisms can be applied to both regolith particles in the optical wavelengths and hydrometeors at microwave frequencies. An advantage from this kind of study is that it does not matter whether the observation is active (e.g., polarimetric radar) or passive (e.g., optical telescope). In both cases, the internal field is computed for two mutually perpendicular incident polarizations, so that the polarization characteristics can then be analyzed according to the relation between these fields and the scattered far field.

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Chromium-manganese iron alloy system design cast in metal and sand moulds for erosion resistance: a positron lifetime study

Erosion characteristics of high chromium (Cr, 16-19%) alloy cast iron with 5% and 10% manganese (Mn) prepared in metal and sand moulds through induction melting are investigated using jet erosion test setup in both as-cast and heat-treated conditions. The samples were characterised for hardness and microstructural properties. A new and novel non-destructive evaluation technique namely positron lifetime spectroscopy has also been used for the first time to characterise the microstructure of the material in terms of defects and their concentration. We found that the hardness decreases irrespective of the sample condition when the mould type is changed from metal to sand, On the other hand, the erosion volume loss shows an increasing trend. Since the macroscopic properties have a bearing on the microstructure, good credence is obtained from the microstructural features as seen from light and scanning electron micrographs. Faster cooling in the metal mould yielded fine carbide precipitation on the surface. The defect size and their concentration derived from positron method are higher for sand mould compared to metal mould. Lower erosion loss corresponds to smaller size defects in metal mould are the results of quicker heat transfer in the metal mould compared to the sand mould. Heat treatment effects are clearly seen as the reduced concentration of defects and spherodisation of carbides points to this. The erosion loss with respect to the defects size and concentration correlate very well.

Symmetry-breaking transition in a two-level system coupled to an Ohmic bath: A squeezed-state approach

Displaced squeezed states are proposed as variational ground states for phonons (Bose fields) coupled to two-level systems (spin systems). We have investigated the zero-temperature phase diagram for the localization-delocalization transition of a tunneling particle interacting with an Ohmic heat bath. Our results are compared with known existing approximate treatments. A modified phase diagram using the displaced squeezed state is presented.

On the spurious increase in particle volume fraction during coarsening simulations

The problem of spurious increase in volume fraction of second-phase particles during computer simulations of coarsening is examined. The origin of this problem is traced to the use of too long a time step (used for numerical integration of growth rates with respect to time) which leads to small particles with large negative growth rates shrinking to negative radii at the end of the time step. Such a shrinkage to negative sizes has the effect of pumping solute into the system. It is therefore suggested that the length of the time step be chosen in accordance with the size of the smallest particle present in the system. It is shown that spurious increase in particle Volume has a significant effect on the particle size distributions in the scaling regime (making them broader and more skewed in the Lifshitz-Slyozov-Wagner model). Its effect on coarsening kinetics, however, is found to be small.