A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model. (C) 2004 Elsevier Inc. All rights reserved.

Analysis of chaotic instabilities in a rotating body with internal energy dissipation

Melnikov's method is used to analytically predict the onset of chaotic instability in a rotating body with internal energy dissipation. The model has been found to exhibit chaotic instability when a harmonic disturbance torque is applied to the system for a range of forcing amplitude and frequency. Such a model may be considered to be representative of the dynamical behavior of a number of physical systems such as a spinning spacecraft. In spacecraft, disturbance torques may arise under malfunction of the control system, from an unbalanced rotor, from vibrations in appendages or from orbital variations. Chaotic instabilities arising from such disturbances could introduce uncertainties and irregularities into the motion of the multibody system and consequently could have disastrous effects on its intended operation. A comprehensive stability analysis is performed and regions of nonlinear behavior are identified. Subsequently, the closed form analytical solution for the unperturbed system is obtained in order to identify homoclinic orbits. Melnikov's method is then applied on the system once transformed into Hamiltonian form. The resulting analytical criterion for the onset of chaotic instability is obtained in terms of critical system parameters. The sufficient criterion is shown to be a useful predictor of the phenomenon via comparisons with numerical results. Finally, for the purposes of providing a complete, self-contained investigation of this fundamental system, the control of chaotic instability is demonstated using Lyapunov's method.

Quantitative scattering of melanin solutions

The optical scattering coefficient of a dilute, well-solubilized eumelanin solution has been accurately measured as a function of incident wavelength, and found to contribute < 6% of the total optical attenuation between 210 and 325 nm. At longer wavelengths (325-800 nm), the scattering was less than the minimum sensitivity of our instrument. This indicates that ultraviolet and visible optical density spectra can be interpreted as true absorption with a high degree of confidence. The scattering coefficient versus wavelength was found to be consistent with Rayleigh theory for a particle radius of 38 6 1 nm. Our results shed important light on the role of melanins as photoprotectants.

An equivalent distributed magnetic current based FDTD method for the calculation of e-fields induced by gradient coils in MRI

This paper evaluates a low-frequency FDTD method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current (DEMC) is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretizing of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modeling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multilayered spherical phantom model and a complete body model.

Canonical solution groups of a homomorphism : manipulator kinematics, nonparametric regression and distributed object systems

The kinematic mapping of a rigid open-link manipulator is a homomorphism between Lie groups. The homomorphisrn has solution groups that act on an inverse kinematic solution element. A canonical representation of solution group operators that act on a solution element of three and seven degree-of-freedom (do!) dextrous manipulators is determined by geometric analysis. Seven canonical solution groups are determined for the seven do! Robotics Research K-1207 and Hollerbach arms. The solution element of a dextrous manipulator is a collection of trivial fibre bundles with solution fibres homotopic to the Torus. If fibre solutions are parameterised by a scalar, a direct inverse funct.ion that maps the scalar and Cartesian base space coordinates to solution element fibre coordinates may be defined. A direct inverse pararneterisation of a solution element may be approximated by a local linear map generated by an inverse augmented Jacobian correction of a linear interpolation. The action of canonical solution group operators on a local linear approximation of the solution element of inverse kinematics of dextrous manipulators generates cyclical solutions. The solution representation is proposed as a model of inverse kinematic transformations in primate nervous systems. Simultaneous calibration of a composition of stereo-camera and manipulator kinematic models is under-determined by equi-output parameter groups in the composition of stereo-camera and Denavit Hartenberg (DH) rnodels. An error measure for simultaneous calibration of a composition of models is derived and parameter subsets with no equi-output groups are determined by numerical experiments to simultaneously calibrate the composition of homogeneous or pan-tilt stereo-camera with DH models. For acceleration of exact Newton second-order re-calibration of DH parameters after a sequential calibration of stereo-camera and DH parameters, an optimal numerical evaluation of DH matrix first order and second order error derivatives with respect to a re-calibration error function is derived, implemented and tested. A distributed object environment for point and click image-based tele-command of manipulators and stereo-cameras is specified and implemented that supports rapid prototyping of numerical experiments in distributed system control. The environment is validated by a hierarchical k-fold cross validated calibration to Cartesian space of a radial basis function regression correction of an affine stereo model. Basic design and performance requirements are defined for scalable virtual micro-kernels that broker inter-Java-virtual-machine remote method invocations between components of secure manageable fault-tolerant open distributed agile Total Quality Managed ISO 9000+ conformant Just in Time manufacturing systems.

Time scale analysis for fluidizedbedmeltgranulation-II: binder spreading rate

The spreading time of liquid binder droplet on the surface a primary particle is analyzed for Fluidized Bed Melt Granulation (FBMG). As discussed in the first paper of this series (Chua et al., in press) the droplet spreading rate has been identified as one of the important parameters affecting the probability of particles aggregation in FBMG. In this paper, the binder droplet spreading time has been estimated using Computational Fluid Dynamic modeling (CFD) based on Volume of Fluid approach (VOF). A simplified analytical solution has been developed and tested to explore its validity for predicting the spreading time. For the purpose of models validation, the droplet spreading evolution was recorded using a high speed video camera. Based on the validated model, a generalized correlative equation for binder spreading time is proposed. For the operating conditions considered here, the spreading time for Polyethylene Glycol (PEG1500) binder was found to fall within the range of 10-2 to 10-5 s. The study also included a number of other common binders used in FBMG. The results obtained here will be further used in paper III, where the binder solidification rate is discussed.

Time scale analysis for fluidizedbedmeltgranulation III: binder solidification rate

In series I and II of this study ([Chua et al., 2010a] and [Chua et al., 2010b]), we discussed the time scale of granule–granule collision, droplet–granule collision and droplet spreading in Fluidized Bed Melt Granulation (FBMG). In this third one, we consider the rate at which binder solidifies. Simple analytical solution, based on classical formulation for conduction across a semi-infinite slab, was used to obtain a generalized equation for binder solidification time. A multi-physics simulation package (Comsol) was used to predict the binder solidification time for various operating conditions usually considered in FBMG. The simulation results were validated with experimental temperature data obtained with a high speed infrared camera during solidification of ‘macroscopic’ (mm scale) droplets. For the range of microscopic droplet size and operating conditions considered for a FBMG process, the binder solidification time was found to fall approximately between 10-3 and 10-1 s. This is the slowest compared to the other three major FBMG microscopic events discussed in this series (granule–granule collision, granule–droplet collision and droplet spreading).

Optical activity in photorefractive Bi12TiO20

The influence of optical activity on two-wave mixing (TWM) in photorefractive BTO and BSO crystals in the absence of an applied field is studied both theoretically and experimentally. For the conventinal orientations of the grating vector, K [001] and K[001], the piezoelectric and photoelastic effects are either zero or negligible. This makes an analytical treatment of the TWM problem possible. We obtain an analytical solution for the coupled wave equations of TWM valid for arbitrary optical activity. This result is of special importance for BTO crystals. In these crystals under the condition of maximum energy transfer (|K|rD=1, where rD is the Debye radius) neither the approximation of small optical activity nor the one of dominating optical activity is applicable and our analytical solution becomes essential. Our experimental setup uses beams with a trapezoidal overlap that allows us to study the thickness-dependence of the gain in a single measurement. Experimental and theoretical results for a BTO crystal are compared with those for a BSO crystal and are explained in the framework of the model used.

Vibration analysis of submerged rectangular microplates with distributed mass loading

This paper investigates the vibration characteristics of the coupling system of a microscale fluid-loaded rectangular isotropic plate attached to a uniformly distributed mass. Previous literature has, respectively, studied the changes in the plate vibration induced by an acoustic field or by the attached mass loading. This paper investigates the issue of involving these two types of loading simultaneously. Based on Lamb's assumption of the fluid-loaded structure and the Rayleigh–Ritz energy method, this paper presents an analytical solution for the natural frequencies and mode shapes of the coupling system. Numerical results for microplates with different types of boundary conditions have also been obtained and compared with experimental and numerical results from previous literature. The theoretical model and novel analytical solution are of particular interest in the design of microplate-based biosensing devices.

The study and measurements of the behaviour and characteristics of electromagnetic relays

The thesis presents a theoretical and practical study of the dynamic behaviour of electromagnetic relays. After discussing the problem of solving the dynamicc equations analytically and presenting a historical survey of the earlier works in the relay and its dynamics, the simulation of a relay on the analogue computer is discussed. It is shown that the simulation may be used to obtain specific solutions to the dynamic equations. The computer analysis provides the dynamic characteristics for design purposes and may be used in the study of bouncing, rebound oscillations and stability of the armature motion. An approximate analytical solution to the two dynamic equations is given based on the assumption that the dynamic variation of the pull with the position of the armature is linear. The assumption is supported by the Computer-aided analysis and experimental results. The solution is intended to provide a basis for a rational design. A rigorous method of analysing the dynamic performance by using Ahlberg's theory is also presented. This method may be justified to be the extension of Ahlberg's theory by taking the mass and frictional damping forces into account. While calculating the armature motion mathematically, Ahlberg considers the equilibrium of two kinds of forces, namely pull and load, and disregards the mass and friction forces, whereas the present method deals with the equilibrium of all four kinds of forces. It is shown how this can be utilised to calculate the dynamic characteristics for a specific design. The utility of this method also extends to the study of stability, contact bounce and armature rebound. The magnetic circuit and other related topics which are essential to the study of relay dynamics are discussed and some necessary experimental results are given.

Optical activity in photorefractive Bi12TiO20

The influence of optical activity on two-wave mixing (TWM) in photorefractive BTO and BSO crystals in the absence of an applied field is studied both theoretically and experimentally. For the conventinal orientations of the grating vector, K [001] and K[001], the piezoelectric and photoelastic effects are either zero or negligible. This makes an analytical treatment of the TWM problem possible. We obtain an analytical solution for the coupled wave equations of TWM valid for arbitrary optical activity. This result is of special importance for BTO crystals. In these crystals under the condition of maximum energy transfer (|K|rD=1, where rD is the Debye radius) neither the approximation of small optical activity nor the one of dominating optical activity is applicable and our analytical solution becomes essential. Our experimental setup uses beams with a trapezoidal overlap that allows us to study the thickness-dependence of the gain in a single measurement. Experimental and theoretical results for a BTO crystal are compared with those for a BSO crystal and are explained in the framework of the model used.

Nonlinearity management in fibre links with distributed amplification

An approach to nonlinearity management in optical transmission lines with periodic dispersion compensation and distributed Raman amplification was presented. The optimization of a three-step dispersion map with forward and backward pumped distributed amplification was examined. The optimization was performed using the analytical solution obtained under the assumption of undepleted pumps and without inclusion of double Rayleigh Scattering (DRS), and by means of a full numerical approach accounting for all important effects. It was found that both procedures led to the same final solution.

Върху два подхода за преподаване на понятията в курса по компютърна топология

Kalin Georgiev, Dimitrina Stavrova - We consider two approaches for introducing topological notions in a course of Computational Topology. One is an intuitive inductive introduction of a notion; the other is theoretical-analytical one. As an example we treat the notions of interior, closure and boundary of a set in a topological space. We analyze several visual representations as well as analytical ones. Examples of tests and quiz problems are considered. Comparisons of student’s achievement on different type of problems are presented.

Monte Carlo Algorithm for Solving Integral Equations with Polynomial Non-Linearity. Parallel Implementation

An iterative Monte Carlo algorithm for evaluating linear functionals of the solution of integral equations with polynomial non-linearity is proposed and studied. The method uses a simulation of branching stochastic processes. It is proved that the mathematical expectation of the introduced random variable is equal to a linear functional of the solution. The algorithm uses the so-called almost optimal density function. Numerical examples are considered. Parallel implementation of the algorithm is also realized using the package ATHAPASCAN as an environment for parallel realization.The computational results demonstrate high parallel efficiency of the presented algorithm and give a good solution when almost optimal density function is used as a transition density.

Conflict-Controlled Processes Involving Fractional Differential Equations with Impulses

MSC 2010: 34A08, 34A37, 49N70