Numerical Simulations and Coherence Properties of Supercontinuum Generation in Photonic Crystal and Tapered Optical Fibers

Numerical simulations have been used to study broad-band supercontinuum generation in optical fibers with dispersion and nonlinearity characteristics typical and photonic crystal or tapered fibers structures. The simulations include optical shock and Raman nonlinearity terms, with quantum noise taken into account phenomenologically by including in the input field a noise seed of one photon per mode with random phase. For input pulses of 150-fs duration injected in the anomalous dispersion regime, the effect of modulational instability is shown to lead to severe temporal jitter in the output, and associated fluctuations in the spectral amplitude and phase across the generated supercontinuum. The spectral phase fluctuations are quantified by performing multiple simulations and calculating both the standard deviation of the phase and, more rigorously, the degree of first-order coherence as a function of wavelength across the spectrum. By performing simulations over a range of input pulse durations and wavelengths, we can identify the conditions under which coherent supercontinua with a well-defined spectral phase are generated.

A volume of fluid numerical model for wave impacts at coastal structures

The first stages in the development of a new design tool, to be used by coastal engineers to improve the efficiency, analysis, design, management and operation of a wide range of coastal and harbour structures, are described. The tool is based on a two-dimensional numerical model, NEWMOTICS-2D, using the volume of fluid (VOF) method, which permits the rapid calculation of wave hydrodynamics at impermeable natural and man-made structures. The critical hydrodynamic flow processes and forces are identified together with the equations that describe these key processes. The different possible numerical approaches for the solution of these equations, and the types of numerical models currently available, are examined and assessed. Preliminary tests of the model, using comparisons with results from a series of hydraulic model test cases, are described. The results of these tests demonstrate that the VOF approach is particularly appropriate for the simulation of the dynamics of waves at coastal structures because of its flexibility in representing the complex free surfaces encountered during wave impact and breaking. The further programme of work, required to develop the existing model into a tool for use in routine engineering design, is outlined.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Experimental and numerical study of the cold crucible melting process

The cold crucible, or induction skull melting process as is otherwise known, has the potential to produce high purity melts of a range of difficult to melt materials, including Ti–Al and Ti6Al4V alloys for Aerospace, Ti–Ta and other biocompatible materials for surgical implants, silicon for photovoltaic and electronic applications, etc. A water cooled AC coil surrounds the crucible causing induction currents to melt the alloy and partially suspend it against gravity away from water-cooled surfaces. Strong stirring takes place in the melt due to the induced electromagnetic Lorentz forces and very high temperatures are attainable under the right conditions (i.e., provided contact with water cooled walls is minimised). In a joint numerical and experimental research programme, various aspects of the design and operation of this process are investigated to increase our understanding of the physical mechanisms involved and to maximise process efficiency. A combination of FV and Spectral CFD techniques are used at Greenwich to tackle this problem numerically, with the experimental work taking place at Birmingham University. Results of this study, presented here, highlight the influence of turbulence and free surface behaviour on attained superheat and also discuss coil design variations and dual frequency options that may lead to winning crucible designs.

Numerical modelling for electromagnetic processing of materials

Electromagnetic processing of materials (EPM) is one of the most widely practiced and fast growing applications of magnetic and electric forces to fluid flow. EPM is encountered in both industrial processes and laboratory investigations. Applications range in scale from nano-particle manipulation to tonnes of liquid metal treated in the presence of various configurations of magnetic fields. Some of these processes are specifically designed and made possible by the use of the electromagnetic force, like the magnetic levitation of liquid droplets, whilst others involve electric currents essential for electrothermal or electrochemical reasons, for instance, in electrolytic metal production and in induction melting. An insight for the range of established and novel EPM applications can be found in the review presented by Asai [1] in the EPM-2003 conference proceedings.

Numerical analysis of bodyworn UHF antenna systems

Bodyworn antennas are found in a wide range of medical, military and personal communication applications, yet reliable communication from the surface of the human body still presents a range of engineering challenges. At UHF and microwave frequencies, bodyworn antennas can suffer from reduced efficiency due to electromagnetic absorption in tissue, radiation pattern fragmentation and variations in feed-point impedance. The significance and nature of these effects are system specific and depend on the operating frequency, propagation environment and physical constraints on the antenna itself. This paper describes how numerical electromagnetic modelling techniques such as FDTD (finite-difference time-domain) can be used in the design of bodyworn antennas. Examples are presented for 418 MHz, 916 .5 MHz and 2 . 45 GHz, in the context of both biomedical signalling and wireless personal-area networking applications such as the Bluetooth(TM)* wireless technology.

Analytical approximations to numerical solutions of theoretical emission measure distributions

Emission line fluxes from cool stars are widely used to establish an apparent emission measure distribution, EmdApp(Te), between temperatures characteristic of the low transition region and the low corona. The true emission measure distribution, EmdTrue(Te), is determined by the energy balance and geometry adopted and, with a numerical model, can be used to predict EmdApp(Te), to guide further modelling. The scaling laws that exist between coronal parameters arise from the dimensions of the terms in the energy balance equation. Here, analytical approximations to numerical solutions for EmdTrue(Te) are presented, which show how the constants in the coronal scaling laws are determined. The apparent emission measure distributions show a minimum value at some T0 and a maximum at the mean coronal temperature Tc (although in some stars, emission from active regions can contribute). It is shown that, for the energy balance and geometry adopted, the analytical values of the emission measure and electron pressure at T0 and Tc depend on only three parameters: the stellar surface gravity and the values of T0 and Tc. The results are tested against full numerical solutions for e Eri (K2 V) and are applied to Procyon (a CMi, F5 IV/V). The analytical approximations can be used to restrict the required range of full numerical solutions, to check the assumed geometry and to show where the adopted energy balance may not be appropriate. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

A Numerical Study of Inlet Geometry for a Low Inertia Mixed Flow Turbocharger Turbine

Mixed flow turbines can offer improvements over typical radial turbines used in automotive turbochargers, with regards to transient performance and low velocity ratio efficiency. Turbine rotor mass dominates the rotating inertia of the turbocharger, and any reductions of mass in the outer radii of the wheel, including the rotor back-disk, can significantly reduce this inertia and improve the acceleration of the assembly. Off-design, low velocity ratio conditions are typified by highly tangential flow at the rotor inlet and a non-zero inlet blade angle is preferred for such operating conditions. This is achievable in a Mixed Flow Turbine without increasing bending stresses within the rotor blade, which is beneficial in high speed and high inlet temperature turbine design. A range of mixed flow turbine rotors was designed with varying cone angle and inlet blade angle and each was assessed at a number of operating points. These rotors were based on an existing radial flow turbine, and both the hub and shroud contours and exducer geometry were maintained. The inertia of each rotor was also considered. The results indicated that there was a trade-off between efficiency and inertia for the rotors and certain designs may be beneficial for the transient performance of downsized, turbocharged engines.

A Review of Numerical Modelling of Wave Energy Converter Arrays

Large-scale commercial exploitation of wave energy is certain to require the deployment of wave energy converters (WECs) in arrays, creating ‘WEC farms’. An understanding of the hydrodynamic interactions in such arrays is essential for determining optimum layouts of WECs, as well as calculating the area of ocean that the farms will require. It is equally important to consider the potential impact of wave farms on the local and distal wave climates and coastal processes; a poor understanding of the resulting environmental impact may hamper progress, as it would make planning consents more difficult to obtain. It is therefore clear that an understanding the interactions between WECs within a farm is vital for the continued development of the wave energy industry.To support WEC farm design, a range of different numerical models have been developed, with both wave phase-resolving and wave phase-averaging models now available. Phase-resolving methods are primarily based on potential flow models and include semi-analytical techniques, boundary element methods and methods involving the mild-slope equations. Phase-averaging methods are all based around spectral wave models, with supra-grid and sub-grid wave farm models available as alternative implementations.The aims, underlying principles, strengths, weaknesses and obtained results of the main numerical methods currently used for modelling wave energy converter arrays are described in this paper, using a common framework. This allows a qualitative comparative analysis of the different methods to be performed at the end of the paper. This includes consideration of the conditions under which the models may be applied, the output of the models and the relationship between array size and computational effort. Guidance for developers is also presented on the most suitable numerical method to use for given aspects of WEC farm design. For instance, certain models are more suitable for studying near-field effects, whilst others are preferable for investigating far-field effects of the WEC farms. Furthermore, the analysis presented in this paper identifies areas in which the numerical modelling of WEC arrays is relatively weak and thus highlights those in which future developments are required.

Integration of numerical models in geographic databases: the case of Sado estuary management

Geographic information systems (GIS) are now widely applied in coastal resource management. Their ability to organise and interface information from a large range of public and private data sources, and their ability to combine this information, using management criteria, to develop a comprehensive picture of the system explains the success of GIS in this area. The use of numerical models as a tool to improve coastal management is also widespread. Less usual is a GIS-based management to ol implementing a comprehensive management model and integrating a numerical modelling system into itself. In this paper such a methodology is proposed. A GIS-based management tool based on the DPSIR model is presented. An overview of the MOHID numerical modelling system is given and the method of integrating this model in the management tool is described. This system is applied to the Sado Estuary (Portugal). Some preliminary results of the integration are presented, demonstrating the capabilities of the management system.

Developing a numerical inverse-theory-based extraction of orientation-dependent relaxation rates from partially- relaxed spectra

Second-rank tensor interactions, such as quadrupolar interactions between the spin- 1 deuterium nuclei and the electric field gradients created by chemical bonds, are affected by rapid random molecular motions that modulate the orientation of the molecule with respect to the external magnetic field. In biological and model membrane systems, where a distribution of dynamically averaged anisotropies (quadrupolar splittings, chemical shift anisotropies, etc.) is present and where, in addition, various parts of the sample may undergo a partial magnetic alignment, the numerical analysis of the resulting Nuclear Magnetic Resonance (NMR) spectra is a mathematically ill-posed problem. However, numerical methods (de-Pakeing, Tikhonov regularization) exist that allow for a simultaneous determination of both the anisotropy and orientational distributions. An additional complication arises when relaxation is taken into account. This work presents a method of obtaining the orientation dependence of the relaxation rates that can be used for the analysis of the molecular motions on a broad range of time scales. An arbitrary set of exponential decay rates is described by a three-term truncated Legendre polynomial expansion in the orientation dependence, as appropriate for a second-rank tensor interaction, and a linear approximation to the individual decay rates is made. Thus a severe numerical instability caused by the presence of noise in the experimental data is avoided. At the same time, enough flexibility in the inversion algorithm is retained to achieve a meaningful mapping from raw experimental data to a set of intermediate, model-free

Numerical modelling of the impact of climate change on the morphology of Saint-Lawrence tributaries

Cette thèse examine les impacts sur la morphologie des tributaires du fleuve Saint-Laurent des changements dans leur débit et leur niveau de base engendrés par les changements climatiques prévus pour la période 2010–2099. Les tributaires sélectionnés (rivières Batiscan, Richelieu, Saint-Maurice, Saint-François et Yamachiche) ont été choisis en raison de leurs différences de taille, de débit et de contexte morphologique. Non seulement ces tributaires subissent-ils un régime hydrologique modifié en raison des changements climatiques, mais leur niveau de base (niveau d’eau du fleuve Saint-Laurent) sera aussi affecté. Le modèle morphodynamique en une dimension (1D) SEDROUT, à l’origine développé pour des rivières graveleuses en mode d’aggradation, a été adapté pour le contexte spécifique des tributaires des basses-terres du Saint-Laurent afin de simuler des rivières sablonneuses avec un débit quotidien variable et des fluctuations du niveau d’eau à l’aval. Un module pour simuler le partage des sédiments autour d’îles a aussi été ajouté au modèle. Le modèle ainsi amélioré (SEDROUT4-M), qui a été testé à l’aide de simulations à petite échelle et avec les conditions actuelles d’écoulement et de transport de sédiments dans quatre tributaires du fleuve Saint-Laurent, peut maintenant simuler une gamme de problèmes morphodynamiques de rivières. Les changements d’élévation du lit et d’apport en sédiments au fleuve Saint-Laurent pour la période 2010–2099 ont été simulés avec SEDROUT4-M pour les rivières Batiscan, Richelieu et Saint-François pour toutes les combinaisons de sept régimes hydrologiques (conditions actuelles et celles prédites par trois modèles de climat globaux (MCG) et deux scénarios de gaz à effet de serre) et de trois scénarios de changements du niveau de base du fleuve Saint-Laurent (aucun changement, baisse graduelle, baisse abrupte). Les impacts sur l’apport de sédiments et l’élévation du lit diffèrent entre les MCG et semblent reliés au statut des cours d’eau (selon qu’ils soient en état d’aggradation, de dégradation ou d’équilibre), ce qui illustre l’importance d’examiner plusieurs rivières avec différents modèles climatiques afin d’établir des tendances dans les effets des changements climatiques. Malgré le fait que le débit journalier moyen et le débit annuel moyen demeurent près de leur valeur actuelle dans les trois scénarios de MCG, des changements importants dans les taux de transport de sédiments simulés pour chaque tributaire sont observés. Ceci est dû à l’impact important de fortes crues plus fréquentes dans un climat futur de même qu’à l’arrivée plus hâtive de la crue printanière, ce qui résulte en une variabilité accrue dans les taux de transport en charge de fond. Certaines complications avec l’approche de modélisation en 1D pour représenter la géométrie complexe des rivières Saint-Maurice et Saint-François suggèrent qu’une approche bi-dimensionnelle (2D) devrait être sérieusement considérée afin de simuler de façon plus exacte la répartition des débits aux bifurcations autour des îles. La rivière Saint-François est utilisée comme étude de cas pour le modèle 2D H2D2, qui performe bien d’un point de vue hydraulique, mais qui requiert des ajustements pour être en mesure de pleinement simuler les ajustements morphologiques des cours d’eau.

Optimization of the Mechanical and Optical Properties of Tunable Optical Sensor Arrays (TOSA) for a Nanospectrometer in the Visible and Near Infrared Spectral Range

Tunable Optical Sensor Arrays (TOSA) based on Fabry-Pérot (FP) filters, for high quality spectroscopic applications in the visible and near infrared spectral range are investigated within this work. The optical performance of the FP filters is improved by using ion beam sputtered niobium pentoxide (Nb2O5) and silicon dioxide (SiO2) Distributed Bragg Reflectors (DBRs) as mirrors. Due to their high refractive index contrast, only a few alternating pairs of Nb2O5 and SiO2 films can achieve DBRs with high reflectivity in a wide spectral range, while ion beam sputter deposition (IBSD) is utilized due to its ability to produce films with high optical purity. However, IBSD films are highly stressed; resulting in stress induced mirror curvature and suspension bending in the free standing filter suspensions of the MEMS (Micro-Electro-Mechanical Systems) FP filters. Stress induced mirror curvature results in filter transmission line degradation, while suspension bending results in high required filter tuning voltages. Moreover, stress induced suspension bending results in higher order mode filter operation which in turn degrades the optical resolution of the filter. Therefore, the deposition process is optimized to achieve both near zero absorption and low residual stress. High energy ion bombardment during film deposition is utilized to reduce the film density, and hence the film compressive stress. Utilizing this technique, the compressive stress of Nb2O5 is reduced by ~43%, while that for SiO2 is reduced by ~40%. Filters fabricated with stress reduced films show curvatures as low as 100 nm for 70 μm mirrors. To reduce the stress induced bending in the free standing filter suspensions, a stress optimized multi-layer suspension design is presented; with a tensile stressed metal sandwiched between two compressively stressed films. The stress in Physical Vapor Deposited (PVD) metals is therefore characterized for use as filter top-electrode and stress compensating layer. Surface micromachining is used to fabricate tunable FP filters in the visible spectral range using the above mentioned design. The upward bending of the suspensions is reduced from several micrometers to less than 100 nm and 250 nm for two different suspension layer combinations. Mechanical tuning of up to 188 nm is obtained by applying 40 V of actuation voltage. Alternatively, a filter line with transmission of 65.5%, Full Width at Half Maximum (FWHM) of 10.5 nm and a stopband of 170 nm (at an output wavelength of 594 nm) is achieved. Numerical model simulations are also performed to study the validity of the stress optimized suspension design for the near infrared spectral range, wherein membrane displacement and suspension deformation due to material residual stress is studied. Two bandpass filter designs based on quarter-wave and non-quarter-wave layers are presented as integral components of the TOSA. With a filter passband of 135 nm and a broad stopband of over 650 nm, high average filter transmission of 88% is achieved inside the passband, while maximum filter transmission of less than 1.6% outside the passband is achieved.

Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

A multiwave range test for obstacle reconstructions with unknown physical properties

We develop a new multiwave version of the range test for shape reconstruction in inverse scattering theory. The range test [R. Potthast, et al., A ‘range test’ for determining scatterers with unknown physical properties, Inverse Problems 19(3) (2003) 533–547] has originally been proposed to obtain knowledge about an unknown scatterer when the far field pattern for only one plane wave is given. Here, we extend the method to the case of multiple waves and show that the full shape of the unknown scatterer can be reconstructed. We further will clarify the relation between the range test methods, the potential method [A. Kirsch, R. Kress, On an integral equation of the first kind in inverse acoustic scattering, in: Inverse Problems (Oberwolfach, 1986), Internationale Schriftenreihe zur Numerischen Mathematik, vol. 77, Birkhäuser, Basel, 1986, pp. 93–102] and the singular sources method [R. Potthast, Point sources and multipoles in inverse scattering theory, Habilitation Thesis, Göttingen, 1999]. In particular, we propose a new version of the Kirsch–Kress method using the range test and a new approach to the singular sources method based on the range test and potential method. Numerical examples of reconstructions for all four methods are provided.