Low-loss directional cloaks without superluminal velocity or magnetic response.

The possibility of making an optically large (many wavelengths in diameter) object appear invisible has been a subject of many recent studies. Exact invisibility scenarios for large (relative to the wavelength) objects involve (meta)materials with superluminal phase velocity [refractive index (RI) less than unity] and/or magnetic response. We introduce a new approximation applicable to certain device geometries in the eikonal limit: piecewise-uniform scaling of the RI. This transformation preserves the ray trajectories but leads to a uniform phase delay. We show how to take advantage of phase delays to achieve a limited (directional and wavelength-dependent) form of invisibility that does not require loss-ridden (meta)materials with superluminal phase velocities.

Tissue-engineered cardiac patch for advanced functional maturation of human ESC-derived cardiomyocytes.

Human embryonic stem cell-derived cardiomyocytes (hESC-CMs) provide a promising source for cell therapy and drug screening. Several high-yield protocols exist for hESC-CM production; however, methods to significantly advance hESC-CM maturation are still lacking. Building on our previous experience with mouse ESC-CMs, we investigated the effects of 3-dimensional (3D) tissue-engineered culture environment and cardiomyocyte purity on structural and functional maturation of hESC-CMs. 2D monolayer and 3D fibrin-based cardiac patch cultures were generated using dissociated cells from differentiated Hes2 embryoid bodies containing varying percentage (48-90%) of CD172a (SIRPA)-positive cardiomyocytes. hESC-CMs within the patch were aligned uniformly by locally controlling the direction of passive tension. Compared to hESC-CMs in age (2 weeks) and purity (48-65%) matched 2D monolayers, hESC-CMs in 3D patches exhibited significantly higher conduction velocities (CVs), longer sarcomeres (2.09 ± 0.02 vs. 1.77 ± 0.01 μm), and enhanced expression of genes involved in cardiac contractile function, including cTnT, αMHC, CASQ2 and SERCA2. The CVs in cardiac patches increased with cardiomyocyte purity, reaching 25.1 cm/s in patches constructed with 90% hESC-CMs. Maximum contractile force amplitudes and active stresses of cardiac patches averaged to 3.0 ± 1.1 mN and 11.8 ± 4.5 mN/mm(2), respectively. Moreover, contractile force per input cardiomyocyte averaged to 5.7 ± 1.1 nN/cell and showed a negative correlation with hESC-CM purity. Finally, patches exhibited significant positive inotropy with isoproterenol administration (1.7 ± 0.3-fold force increase, EC50 = 95.1 nm). These results demonstrate highly advanced levels of hESC-CM maturation after 2 weeks of 3D cardiac patch culture and carry important implications for future drug development and cell therapy studies.

Cyclic Pursuit: Symmetry, Reduction and Nonlinear Dynamics

In this dissertation, we explore the use of pursuit interactions as a building block for collective behavior, primarily in the context of constant bearing (CB) cyclic pursuit. Pursuit phenomena are observed throughout the natural environment and also play an important role in technological contexts, such as missile-aircraft encounters and interactions between unmanned vehicles. While pursuit is typically regarded as adversarial, we demonstrate that pursuit interactions within a cyclic pursuit framework give rise to seemingly coordinated group maneuvers. We model a system of agents (e.g. birds, vehicles) as particles tracing out curves in the plane, and illustrate reduction to the shape space of relative positions and velocities. Introducing the CB pursuit strategy and associated pursuit law, we consider the case for which agent i pursues agent i+1 (modulo n) with the CB pursuit law. After deriving closed-loop cyclic pursuit dynamics, we demonstrate asymptotic convergence to an invariant submanifold (corresponding to each agent attaining the CB pursuit strategy), and proceed by analysis of the reduced dynamics restricted to the submanifold. For the general setting, we derive existence conditions for relative equilibria (circling and rectilinear) as well as for system trajectories which preserve the shape of the collective (up to similarity), which we refer to as pure shape equilibria. For two illustrative low-dimensional cases, we provide a more comprehensive analysis, deriving explicit trajectory solutions for the two-particle "mutual pursuit" case, and detailing the stability properties of three-particle relative equilibria and pure shape equilibria. For the three-particle case, we show that a particular choice of CB pursuit parameters gives rise to remarkable almost-periodic trajectories in the physical space. We also extend our study to consider CB pursuit in three dimensions, deriving a feedback law for executing the CB pursuit strategy, and providing a detailed analysis of the two-particle mutual pursuit case. We complete the work by considering evasive strategies to counter the motion camouflage (MC) pursuit law. After demonstrating that a stochastically steering evader is unable to thwart the MC pursuit strategy, we propose a (deterministic) feedback law for the evader and demonstrate the existence of circling equilibria for the closed-loop pursuer-evader dynamics.

Investigation and Development of a Fully 3D Tilt Capable Hybrid SPECT - CT System for Dedicated Breast Imaging

X-ray mammography has been the gold standard for breast imaging for decades, despite the significant limitations posed by the two dimensional (2D) image acquisitions. Difficulty in diagnosing lesions close to the chest wall and axilla, high amount of structural overlap and patient discomfort due to compression are only some of these limitations. To overcome these drawbacks, three dimensional (3D) breast imaging modalities have been developed including dual modality single photon emission computed tomography (SPECT) and computed tomography (CT) systems. This thesis focuses on the development and integration of the next generation of such a device for dedicated breast imaging. The goals of this dissertation work are to: [1] understand and characterize any effects of fully 3-D trajectories on reconstructed image scatter correction, absorbed dose and Hounsifeld Unit accuracy, and [2] design, develop and implement the fully flexible, third generation hybrid SPECT-CT system capable of traversing complex 3D orbits about a pendant breast volume, without interference from the other. Such a system would overcome artifacts resulting from incompletely sampled divergent cone beam imaging schemes and allow imaging closer to the chest wall, which other systems currently under research and development elsewhere cannot achieve.

The dependence of x-ray scatter radiation on object shape, size, material composition and the CT acquisition trajectory, was investigated with a well-established beam stop array (BSA) scatter correction method. While the 2D scatter to primary ratio (SPR) was the main metric used to characterize total system scatter, a new metric called ‘normalized scatter contribution’ was developed to compare the results of scatter correction on 3D reconstructed volumes. Scatter estimation studies were undertaken with a sinusoidal saddle (±15° polar tilt) orbit and a traditional circular (AZOR) orbit. Clinical studies to acquire data for scatter correction were used to evaluate the 2D SPR on a small set of patients scanned with the AZOR orbit. Clinical SPR results showed clear dependence of scatter on breast composition and glandular tissue distribution, otherwise consistent with the overall phantom-based size and density measurements. Additionally, SPR dependence was also observed on the acquisition trajectory where 2D scatter increased with an increase in the polar tilt angle of the system.

The dose delivered by any imaging system is of primary importance from the patient’s point of view, and therefore trajectory related differences in the dose distribution in a target volume were evaluated. Monte Carlo simulations as well as physical measurements using radiochromic film were undertaken using saddle and AZOR orbits. Results illustrated that both orbits deliver comparable dose to the target volume, and only slightly differ in distribution within the volume. Simulations and measurements showed similar results, and all measured dose values were within the standard screening mammography-specific, 6 mGy dose limit, which is used as a benchmark for dose comparisons.

Hounsfield Units (HU) are used clinically in differentiating tissue types in a reconstructed CT image, and therefore the HU accuracy of a system is very important, especially when using non-traditional trajectories. Uniform phantoms filled with various uniform density fluids were used to investigate differences in HU accuracy between saddle and AZOR orbits. Results illustrate the considerably better performance of the saddle orbit, especially close to the chest and nipple region of what would clinically be a pedant breast volume. The AZOR orbit causes shading artifacts near the nipple, due to insufficient sampling, rendering a major portion of the scanned phantom unusable, whereas the saddle orbit performs exceptionally well and provides a tighter distribution of HU values in reconstructed volumes.

Finally, the third generation, fully-suspended SPECT-CT system was designed in and developed in our lab. A novel mechanical method using a linear motor was developed for tilting the CT system. A new x-ray source and a custom made 40 x 30 cm2 detector were integrated on to this system. The SPECT system was nested, in the center of the gantry, orthogonal to the CT source-detector pair. The SPECT system tilts on a goniometer, and the newly developed CT tilting mechanism allows ±15° maximum polar tilting of the CT system. The entire gantry is mounted on a rotation stage, allowing complex arbitrary trajectories for each system, without interference from the other, while having a common field of view. This hybrid system shows potential to be used clinically as a diagnostic tool for dedicated breast imaging.

Proper-motion binaries in the Hipparcos catalogue: comparison with radial velocity data

Context. This paper is the last in a series devoted to the analysis of the binary content of the Hipparcos Catalogue. Aims. The comparison of the proper motions constructed from positions spanning a short (Hipparcos) or long time (Tycho-2) makes it possible to uncover binaries with periods of the order of or somewhat larger than the short time span (in this case, the 3 yr duration of the Hipparcos mission), since the unrecognised orbital motion will then add to the proper motion. Methods. A list of candidate proper motion binaries is constructed from a carefully designed χ2 test evaluating the statistical significance of the difference between the Tycho-2 and Hipparcos proper motions for 103 134 stars in common between the two catalogues (excluding components of visual systems). Since similar lists of proper-motion binaries have already been constructed, the present paper focuses on the evaluation of the detection efficiency of proper-motion binaries, using different kinds of control data (mostly radial velocities). The detection rate for entries from the Ninth Catalogue of Spectroscopic Binary Orbits (SB9) is evaluated, as well as for stars like barium stars, which are known to be all binaries, and finally for spectroscopic binaries identified from radial velocity data in the Geneva-Copenhagen survey of F and G dwarfs in the solar neighbourhood. Results. Proper motion binaries are efficiently detected for systems with parallaxes in excess of ∼20 mas, and periods in the range 1000-30 000 d. The shortest periods in this range (1000-2000 d, i.e. once to twice the duration of the Hipparcos mission) may appear only as DMSA/G binaries (accelerated proper motion in the Hipparcos Double and Multiple System Annex). Proper motion binaries detected among SB9 systems having periods shorter than about 400 d hint at triple systems, the proper-motion binary involving a component with a longer orbital period. A list of 19 candidate triple systems is provided. Binaries suspected of having low-mass (brown-dwarf-like) companions are listed as well. Among the 37 barium stars with parallaxes larger than 5 mas, only 7 exhibit no evidence for duplicity whatsoever (be it spectroscopic or astrometric). Finally, the fraction of proper-motion binaries shows no significant variation among the various (regular) spectral classes, when due account is taken for the detection biases. © ESO 2007.

Nanoscopic modelling of the adhesion, indentation and fracture characteristics of metallic systems via molecular dynamics simulations

Large-scale molecular dynamics simulations have been performed on canonical ensembles to model the adhesion and indentation characteristics of 3-D metallic nano-scale junctions in tip-substrate geometries, and the crack propagation in 2-D metallic lattices. It is shown that irreversible flows in nano-volumes of materials control the behaviour of the 3-D nano-contacts, and that local diffusional flow constitutes the atomistic mechanism underlying these plastic flows. These simulations show that the force of adhesion in metallic nano-contacts is reduced when adsorbate monolayers are present at the metal—metal junctions. Our results are in agreement with the conclusions of very accurate point-contact experiments carried out in this field. Our fracture simulations reveal that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures, the nucleation of dislocations is shown to cause a brittle-to-ductile transition. Limiting crack propagation velocities are computed for different strain rates and a dynamic instability is shown to control the crack movement beyond this limiting velocity, in line with the recent experimental results.

An experimental and numerical CFD study of turbulence in a tundish container

This paper describes work performed at IRSID/USINOR in France and the University of Greenwich, UK, to investigate flow structures and turbulence in a water-model container, simulating aspects typical of metal tundish operation. Extensive mean and fluctuating velocity measurements were performed at IRSID using LDA to determine the flow field and these form the basis for a numerical model validation. This apparently simple problem poses several difficulties for the CFD modelling. The flow is driven by the strong impinging jet at the inlet. Accurate description of the jet is most important and requires a localized fine grid, but also a turbulence model that predicts the correct spreading rates of jet and impinging wall boundary layers. The velocities in the bulk of the tundish tend to be (indeed need to be) much smaller than those of the jet, leading to damping of turbulence, or even laminar flow. The authors have developed several low-Reynolds number (low-Re) k–var epsilon model variants to compute this flow and compare against measurements. Best agreement is obtained when turbulence damping is introduced to account not only for walls, but also for low-Re regions in the bulk – the k–var epsilon model otherwise allows turbulence to accumulate in the container due to the restricted outlet. Several damping functions are tested and the results reported here. The k–ω model, which is more suited to transitional flow, also seems to perform well in this problem.

Characterisation of weldpool shapes in the laser welding of thick plates

We consider the problem of finding the heat distribution and the shape of the liquid fraction during laser welding of a thick steel plate using the finite volume CFD package PHYSICA. Since the shape of the keyhole is not known in advance, the following two-step approach to handling this problem has been employed. In the first stage, we determine the geometry of the keyhole for the steady-state case and form an appropriate mesh that includes both the workpiece and the keyhole. In the second stage, we impose the boundary conditions by assigning temperature to the walls of the keyhole and find the heat distribution and the shape of the liquid fraction for a given welding speed and material properties. We construct a fairly accurate approximation of the keyhole as a sequence of include sliced cones. A formula for finding the initial radius of the keyhole is derived by determining the radius of the vaporisation isotherm for the line heat source. We report on the results of a series of computational experiments for various heat input values and welding velocities.

AC and DC magnetic levitation: melting, fluid flow and oscillations

The presented numerical modelling for the magnetic levitation involves coupling of the electromagnetic field, liquid shape change, fluid velocities and the temperature field at every time step during the simulation in time evolution. Combination of the AC and DC magnetic fields can be used to achieve high temperature, stable levitation conditions. The oscillation frequency spectra are analysed for droplets levitated in AC and DC magnetic fields at various combinations. An electrically poorly conducting, diamagnetic droplet (e.g. water) can be stably levitated using the dia- and para-magnetic properties of the sample material in a high intensity, gradient DC field.

Computational model for prediction of particle degradation during dilute phase pneumatic conveying: the use of a laboratory scale degradation tester for the determination of degradation propensity

The overall objective of this work is to develop a computational model of particle degradation during dilute-phasepneumatic conveying. A key feature of such a model is the prediction of particle breakage due to particle–wall collisions in pipeline bends. This paper presents a method for calculating particle impact degradation propensity under a range of particle velocities and particle sizes. It is based on interpolation on impact data obtained in a new laboratory-scale degradation tester. The method is tested and validated against experimental results for degradation at 90± impact angle of a full-size distribution sample of granulated sugar. In a subsequent work, the calculation of degradation propensity is coupled with a ow model of the solids and gas phases in the pipeline.

Investigation into the performance of turbulence models for fluid flow and heat transfer phenomena in electronic applications

Computational Fluid Dynamics (CFD) is gradually becoming a powerful and almost essential tool for the design, development and optimization of engineering applications. However the mathematical modelling of the erratic turbulent motion remains the key issue when tackling such flow phenomena. The reliability of CFD analysis depends heavily on the turbulence model employed together with the wall functions implemented. In order to resolve the abrupt changes in the turbulent energy and other parameters situated at near wall regions a particularly fine mesh is necessary which inevitably increases the computer storage and run-time requirements. Turbulence modelling can be considered to be one of the three key elements in CFD. Precise mathematical theories have evolved for the other two key elements, grid generation and algorithm development. The principal objective of turbulence modelling is to enhance computational procedures of efficient accuracy to reproduce the main structures of three dimensional fluid flows. The flow within an electronic system can be characterized as being in a transitional state due to the low velocities and relatively small dimensions encountered. This paper presents simulated CFD results for an investigation into the predictive capability of turbulence models when considering both fluid flow and heat transfer phenomena. Also a new two-layer hybrid kε / kl turbulence model for electronic application areas will be presented which holds the advantages of being cheap in terms of the computational mesh required and is also economical with regards to run-time.

Computational modelling of bubbles, droplets and particles in metals reduction and refining

A multi-phase framework is typically required for the CFD modelling of metals reduction processes. Such processes typically involve the interaction of liquid metals, a gas (often air) top space, liquid droplets in the top space and injection of both solid particles and gaseous bubbles into the bath. The exchange of mass, momentum and energy between the phases is fundamental to these processes. Multi-phase algorithms are complex and can be unreliable in terms of either or both convergence behaviour or in the extent to which the physics is captured. In this contribution, we discuss these multi-phase flow issues and describe an example of each of the main “single phase” approaches to modelling this class of problems (i.e., Eulerian–Lagrangian and Eulerian–Eulerian). Their utility is illustrated in the context of two problems – one involving the injection of sparging gases into a steel continuous slab caster and the other based on the development of a novel process for aluminium electrolysis. In the steel caster, the coupling of the Lagrangian tracking of the gas phase with the continuum enables the simulation of the transient motion of the metal–flux interface. The model of the electrolysis process employs a novel method for the calculation of slip velocities of oxygen bubbles, resulting from the dissolution of alumina, which allows the efficiency of the process to be predicted.

Multiscale modelling of dendrite growth during vacuum arc remelting

Vacuum Arc Remelting (VAR) is the accepted method for producing homogeneous, fine microstructures that are free of inclusions required for rotating grade applications. However, as ingot sizes are increasing INCONEL 718 becomes increasingly susceptible to defects such as freckles, tree rings, and white spots increases for large diameter billets. Therefore, predictive models of these defects are required to allow optimization of process parameters. In this paper, a multiscale and multi-physics model is presented to predict the development of microstructures in the VAR ingot during solidification. At the microscale, a combined stochastic nucleation approach and finite difference solution of the solute diffusion is applied in the semi-solid zone of the VAR ingot. The micromodel is coupled with a solution of the macroscale heat transfer, fluid flow and electromagnetism in the VAR process through the temperature, pressure and fluid flow fields. The main objective of this study is to achieve a better understanding of the formation of the defects in VAR by quantifying the influence of VAR processing parameters on grain nucleation and dendrite growth. In particular, the effect of different ingot growth velocities on the microstructure formation was investigated. It was found that reducing the velocity produces significantly more coarse grains.

Levitated liquid droplets in AC and DC magnetic field

Electromagnetic levitation of liquid metal droplets can be used to measure the properties of highly reactive liquid materials. Two independent numerical models, the commercial COMSOL and the spectral-collocation based free surface code SPHINX, have been applied to solve the transient electromagnetic, fluid flow and thermodynamic equations, which describe the levitated liquid motion and heating processes. The SPHINX model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the electromagnetic and gravity forces, temperature dependent surface tension, magnetically controlled turbulent momentum transport. The models are adapted to incorporate periodic laser heating at the top of the droplet, which is used to measure the thermal conductivity of the material. Novel effects in the levitated droplet of magnetically damped turbulence and nonlinear growth of velocities in high DC magnetic field are analysed.

Parametrical modeling and design optimization of blood plasma separation device with microchannel mechanism

This paper presents an analysis of biofluid behavior in a T-shaped microchannel device and a design optimization for improved biofluid performance in terms of particle liquid separation. The biofluid is modeled with single phase shear rate non-Newtonian flow with blood property. The separation of red blood cell from plasma is evident based on biofluid distribution in the microchannels against various relevant effects and findings, including Zweifach-Fung bifurcation law, Fahraeus effect, Fahraeus-Lindqvist effect and cell free phenomenon. The modeling with the initial device shows that this T-microchannel device can separate red blood cell from plasma but the separation efficiency among different bifurcations varies largely. In accordance with the imbalanced performance, a design optimization is conducted. This includes implementing a series of simulations to investigate the effect of the lengths of the main and branch channels to biofluid behavior and searching an improved design with optimal separation performance. It is found that changing relative lengths of branch channels is effective to both uniformity of flow rate ratio among bifurcations and reduction of difference of the flow velocities between the branch channels, whereas extending the length of the main channel from bifurcation region is only effective for uniformity of flow rate ratio.