Numerical evaluation of e-fields induced by body motion near high-field MRI scanner

In modern magnetic resonance imaging (MRI), both patients and radiologists are exposed to strong, nonuniform static magnetic fields inside or outside of the scanner, in which the body movement may be able to induce electric currents in tissues which could be possibly harmful. This paper presents theoretical investigations into the spatial distribution of induced E-fields in the human model when moving at various positions around the magnet. The numerical calculations are based on an efficient, quasistatic, finite-difference scheme and an anatomically realistic, full-body, male model. 3D field profiles from an actively-shielded 4 T magnet system are used and the body model projected through the field profile with normalized velocity. The simulation shows that it is possible to induce E-fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The results are easy to extrapolate to very high field strengths for the safety evaluation at a variety of field strengths and motion velocities.

Velocity distributions in a plate heat exchanger channel

The development of a Laser Doppler Anemometer technique to measure the velocity distribution in a commercial plate heat exchanger is described. Detailed velocity profiles are presented and a preliminary investigation is reported on flow behaviour through a single cell in the channel matrix. The objective of the study was to extend previous investigations of plate heat exchanger flow patterns in the laminar range with the eventual aim of establishing the effect of flow patterns on heat transfer performance, thus leading to improved plate heat exchanger design and design methods. Accurate point velocities were obtained by Laser Anemometry in a perspex replica of the metal channel. Oil was used as a circulating liquid with a refractive index matched to that of the perspex so that the laser beams were not distorted. Cell-by-cell velocity measurements over a range of Reynolds number up to ten showed significant liquid mal-distribution. Local cell velocities were found to be as high as twenty seven times average velocity, contrary to the previously held belief of four times. The degree of mal-distribution varied across the channel as well as in the vertical direction, and depended on the upward or downward direction of flow. At Reynolds numbers less than one, flow zig-zagged from one side of the channel to the other in wave form, but increases in Reynolds number improved liquid distribution. A detailed examination of selected cells showed velocity variations in different directions, together with variation within individual cells. Experimental results are also reported on the flow split when passing through a single cell in a section of a channel . These observations were used to explain mal-distribution in the perspex channel itself.

The application of H controller synthesis to high speed independent drive systems

This thesis describes work completed on the application of H controller synthesis to the design of controllers for single axis high speed independent drive design examples. H controller synthesis was used in a single controller format and in a self-tuning regulator, a type of adaptive controller. Three types of industrial design examples were attempted using H controller synthesis, both in simulation and on a Drives Test Facility at Aston University. The results were benchmarked against a Proportional, Integral and Derivative (PID) with velocity feedforward controller (VFF), the industrial standard for this application. An analysis of the differences between a H and PID with VFF controller was completed. A direct-form H controller was determined for a limited class of weighting function and plants which shows the relationship between the weighting function, nominal plant and the controller parameters. The direct-form controller was utilised in two ways. Firstly it allowed the production of simple guidelines for the industrial design of H controllers. Secondly it was used as the controller modifier in a self-tuning regulator (STR). The STR had a controller modification time (including nominal model parameter estimation) of 8ms. A Set-Point Gain Scheduling (SPGS) controller was developed and applied to an industrial design example. The applicability of each control strategy, PID with VFF, H, SPGS and STR, was investigated and a set of general guidelines for their use was determined. All controllers developed were implemented using standard industrial equipment.

Inscription and characterization of waveguides written into borosilicate glass by a high-repetition-rate femtosecond laser at 800 nm

A series of waveguides was inscribed in a borosilicate glass (BK7) by an 11 MHz repetition rate femtosecond laser operating with pulse energies from 16 to 30 nJ and focused at various depths within the bulk material. The index modification was measured using a quantitative phase microscopy technique that revealed central index changes ranging from 5×10-3 to 10-2, leading to waveguides that exhibited propagation losses of 0.2 dB/cm at a wavelength of 633 nm and 0.6 dB/cm at a wavelength of 1550 nm with efficient mode matching, less than 0.2 dB, to standard optical fibers. Analysis of the experimental data shows that, for a given inscription energy, the index modification has a strong dependence on inscription scanning velocity. At higher energies, the index modification increases with increasing inscription scanning velocity with other fabrication parameters constant.

A compartmental CFD-PBM model of high shear wet granulation

The conventional, geometrically lumped description of the physical processes inside a high shear granulator is not reliable for process design and scale-up. In this study, a compartmental Population Balance Model (PBM) with spatial dependence is developed and validated in two lab-scale high shear granulation processes using a 1.9L MiPro granulator and 4L DIOSNA granulator. The compartmental structure is built using a heuristic approach based on computational fluid dynamics (CFD) analysis, which includes the overall flow pattern, velocity and solids concentration. The constant volume Monte Carlo approach is implemented to solve the multi-compartment population balance equations. Different spatial dependent mechanisms are included in the compartmental PBM to describe granule growth. It is concluded that for both cases (low and high liquid content), the adjustment of parameters (e.g. layering, coalescence and breakage rate) can provide a quantitative prediction of the granulation process.

Derivation of lake mixing and stratification indices from high-resolution lake buoy data

Lake Analyzer is a numerical code coupled with supporting visualization tools for determining indices of mixing and stratification that are critical to the biogeochemical cycles of lakes and reservoirs. Stability indices, including Lake Number, Wedderburn Number, Schmidt Stability, and thermocline depth are calculated according to established literature definitions and returned to the user in a time series format. The program was created for the analysis of high-frequency data collected from instrumented lake buoys, in support of the emerging field of aquatic sensor network science. Available outputs for the Lake Analyzer program are: water temperature (error-checked and/or down-sampled), wind speed (error-checked and/or down-sampled), metalimnion extent (top and bottom), thermocline depth, friction velocity, Lake Number, Wedderburn Number, Schmidt Stability, mode-1 vertical seiche period, and Brunt-Väisälä buoyancy frequency. Secondary outputs for several of these indices delineate the parent thermocline depth (seasonal thermocline) from the shallower secondary or diurnal thermocline. Lake Analyzer provides a program suite and best practices for the comparison of mixing and stratification indices in lakes across gradients of climate, hydro-physiography, and time, and enables a more detailed understanding of the resulting biogeochemical transformations at different spatial and temporal scales.

Physical properties measured on sediment cores from high latitudes

This collection of 359 data sets represents raw data of physical properties measurements on Polarstern sediment cores from both polar oceans, sampled and measured between 1985 and 1995.

High-speed label-free functional photoacoustic microscopy of mouse brain in action.

We present fast functional photoacoustic microscopy (PAM) for three-dimensional high-resolution, high-speed imaging of the mouse brain, complementary to other imaging modalities. We implemented a single-wavelength pulse-width-based method with a one-dimensional imaging rate of 100 kHz to image blood oxygenation with capillary-level resolution. We applied PAM to image the vascular morphology, blood oxygenation, blood flow and oxygen metabolism in both resting and stimulated states in the mouse brain.

Near real time surface ocean velocity, Tyrrhenian Sea and Ligurian Sea in September 2016

The data set consists of maps of total velocity of the surface current in the North-Western Tyrrhenian Sea and Ligurian Sea averaged over a time interval of 1 hour around the cardinal hour. Surface ocean velocities estimated by HF Radar are representative of the upper 0.3-2.5 meters of the ocean. Total velocities are derived using least square fit that maps radial velocities measured from individual sites onto a cartesian grid. The final product is a map of the horizontal components of the ocean currents on a regular grid in the area of overlap of two or more radar stations.

Near real time surface ocean velocity in the Bay of Biscay - Basque Country from Mai to September 2016

The data set consists of maps of total velocity of the surface current in the Southeastern Bay of Biscay averaged over a time interval of 1 hour around the cardinal hour. Surface ocean velocities estimated by this HF Radar(4.65 MHz) are representative of the upper 2-3 meters of the ocean. The main objective of near real time processing is to produce the best product from available data at the time of processing. Total velocities are derived using least square fit that maps radial velocities measured from individual sites onto a cartesian grid. The final product is a map of the horizontal components of the ocean currents on a regular grid in the area of overlap of two or more radar stations.

Seawater carbonate chemistry and its effects on properties and sinking velocity of aggregates of the coccolithophore Emiliania huxleyi, 2010

Coccolithophores play an important role in organic matter export due to their production of the mineral calcite that can act as ballast. Recent studies indicated that calcification in coccolithophores may be affected by changes in seawater carbonate chemistry. We investigated the influence of CO2 on the aggregation and sinking behaviour of the coccolithophore Emiliania huxleyi (PML B92/11) during a laboratory experiment. The coccolithophores were grown under low (~180 µatm), medium (~380 µatm), and high (~750 µatm) CO2 conditions. Aggregation of the cells was promoted using roller tables. Size and settling velocity of aggregates were determined during the incubation using video image analysis. Our results indicate that aggregate properties are sensitive to changes in the degree of ballasting, as evoked by ocean acidification. Average sinking velocity was highest for low CO2 aggregates (~1292 m d-1) that also had the highest particulate inorganic to particulate organic carbon (PIC/POC) ratio. Lowest PIC/POC ratios and lowest sinking velocity (~366 m d-1) at comparable sizes were observed for aggregates of the high CO2 treatment. Aggregates of the high CO2 treatment showed a 4-fold lower excess density (~4.2*10**-4 g cm**-3) when compared to aggregates from the medium and low CO2 treatments (~1.7 g*10**-3 cm**-3). We also observed that more aggregates formed in the high CO2 treatment, and that those aggregates contained more bacteria than aggregates in the medium and low CO2 treatment. If applicable to the future ocean, our findings suggest that a CO2 induced reduction of the calcite content of aggregates could weaken the deep export of organic matter in the ocean, particularly in areas dominated by coccolithophores.

DEVELOPMENT OF A CONDENSED PHASE VELOCITY MAP IMAGING APPARATUS: DISSOCIATION PHOTODYNAMICS OF NITROGEN DIOXIDE AT 355 NM

The photochemistry of the polar regions of Earth, as well as the interstellar medium, is driven by the effect of ultraviolet radiation on ice surfaces and on the materials trapped within them. While the area of ice photochemistry is vast and much research has been completed, it has only recently been possible to study the dynamics of these processes on a microscopic level. One of the leading techniques for studying photoreaction dynamics is Velocity Map Imaging (VMI). This technique has been used extensively to study several types of reaction dynamics processes. Although the majority of these studies have utilized molecular beams as the main medium for reactants, new studies showed the versatility of the technique when applied to molecular dynamics of molecules adsorbed on metal surfaces. Herein the development of a velocity map imaging apparatus capable of studying the photochemistry of condensed phase materials is described. The apparatus is used to study of the photo-reactivity of NO2 condensed within argon matrices to illustrate its capabilities. A doped ice surface is formed by condensing Ar and NO2 gas onto a sapphire rod which is cooled using a helium compressor to 20 K. The matrix is irradiated using an Nd:YAG laser at 355 nm, and the resulting NO fragment is state-selectively ionized using an excimer-pumped dye laser. In all, we are able to detect transient photochemically generated species and can collect information on their quantum state and kinetic energy distribution. It is found that the REMPI spectra changes as different sections of the dissociating cloud are probed. The rotational and translational energy populations are found to be bimodal with a low temperature component roughly at the temperature of the matrix, and a second component with much higher temperature, the rotational temperature showing a possible population inversion, and the translational temperature of 100-200 K. The low temperature translational component is found to dominate at long delay times between dissociation and ionization, while at short time delays the high temperature component plays a larger role. The velocity map imaging technique allows for the detection of both the axial and radial components of the translational energy. The distribution of excess energy over the rotational, electronic and translational states of the NO photofragments provides evidence for collisional quenching of the fragments in the Ar-matrix prior to their desorption.

Optimising power take-off of an oscillating wave surge converter using high fidelity numerical simulations

Oscillating wave surge converters are a promising technology to harvest ocean wave energy in the near shore region. Although research has been going on for many years, the characteristics of the wave action on the structure and especially the phase relation between the driving force and wave quantities like velocity or surface elevation have not been investigated in detail. The main reason for this is the lack of suitable methods. Experimental investigations using tank tests do not give direct access to overall hydrodynamic loads, only damping torque of a power take off system can be measured directly. Non-linear computational fluid dynamics methods have only recently been applied in the research of this type of devices. This paper presents a new metric named wave torque, which is the total hydrodynamic torque minus the still water pitch stiffness at any given angle of rotation. Changes in characteristics of that metric over a wave cycle and for different power take off settings are investigated using computational fluid dynamics methods. Firstly, it is shown that linearised methods cannot predict optimum damping in typical operating states of OWSCs. We then present phase relationships between main kinetic parameters for different damping levels. Although the flap seems to operate close to resonance, as predicted by linear theory, no obvious condition defining optimum damping is found.

Multiphase modeling of melting / solidification with high density variations using XFEM

La modélisation de la cryolite, utilisée dans la fabrication de l’aluminium, implique plusieurs défis, notament la présence de discontinuités dans la solution et l’inclusion de la difference de densité entre les phases solide et liquide. Pour surmonter ces défis, plusieurs éléments novateurs ont été développés dans cette thèse. En premier lieu, le problème du changement de phase, communément appelé problème de Stefan, a été résolu en deux dimensions en utilisant la méthode des éléments finis étendue. Une formulation utilisant un multiplicateur de Lagrange stable spécialement développée et une interpolation enrichie a été utilisée pour imposer la température de fusion à l’interface. La vitesse de l’interface est déterminée par le saut dans le flux de chaleur à travers l’interface et a été calculée en utilisant la solution du multiplicateur de Lagrange. En second lieu, les effets convectifs ont été inclus par la résolution des équations de Stokes dans la phase liquide en utilisant la méthode des éléments finis étendue aussi. Troisièmement, le changement de densité entre les phases solide et liquide, généralement négligé dans la littérature, a été pris en compte par l’ajout d’une condition aux limites de vitesse non nulle à l’interface solide-liquide pour respecter la conservation de la masse dans le système. Des problèmes analytiques et numériques ont été résolus pour valider les divers composants du modèle et le système d’équations couplés. Les solutions aux problèmes numériques ont été comparées aux solutions obtenues avec l’algorithme de déplacement de maillage de Comsol. Ces comparaisons démontrent que le modèle par éléments finis étendue reproduit correctement le problème de changement phase avec densités variables.