Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

Customer information system data pre-processing with feature selection techniques for non-technical losses prediction in an electricity market

Non-technical losses (NTL) identification and prediction are important tasks for many utilities. Data from customer information system (CIS) can be used for NTL analysis. However, in order to accurately and efficiently perform NTL analysis, the original data from CIS need to be pre-processed before any detailed NTL analysis can be carried out. In this paper, we propose a feature selection based method for CIS data pre-processing in order to extract the most relevant information for further analysis such as clustering and classifications. By removing irrelevant and redundant features, feature selection is an essential step in data mining process in finding optimal subset of features to improve the quality of result by giving faster time processing, higher accuracy and simpler results with fewer features. Detailed feature selection analysis is presented in the paper. Both time-domain and load shape data are compared based on the accuracy, consistency and statistical dependencies between features.

Tailoring magnetic field gradient design to magnet cryostat geometry

Eddy currents induced within a magnetic resonance imaging (MRI) cryostat bore during pulsing of gradient coils can be applied constructively together with the gradient currents that generate them, to obtain good quality gradient uniformities within a specified imaging volume over time. This can be achieved by simultaneously optimizing the spatial distribution and temporal pre-emphasis of the gradient coil current, to account for the spatial and temporal variation of the secondary magnetic fields due to the induced eddy currents. This method allows the tailored design of gradient coil/magnet configurations and consequent engineering trade-offs. To compute the transient eddy currents within a realistic cryostat vessel, a low-frequency finite-difference time-domain (FDTD) method using total-field scattered-field (TFSF) scheme has been performed and validated

FDTD study of an UWB radar technique for breast tumor detection

A finite difference time domain (FDTD) method is applied to investigate capabilities of an ultra-wide band (UWB) radar system to detect a breast tumor. The first part of the investigations concerns FDTD simulations of a phantom formed by a plastic container with liquid and a small reflecting target. The second part focuses on a three-dimensional numerical breast model with a small tumor. FDTD simulations are carried out assuming a planar incident wave. Various time snap shots of the electromagnetic field are recorded to learn about the physical phenomenon of reflection and scattering in different layers of the phantom.

Mathematical modelling and optimal multivariable control of chemical processes

This work reports the developnent of a mathenatical model and distributed, multi variable computer-control for a pilot plant double-effect climbing-film evaporator. A distributed-parameter model of the plant has been developed and the time-domain model transformed into the Laplace domain. The model has been further transformed into an integral domain conforming to an algebraic ring of polynomials, to eliminate the transcendental terms which arise in the Laplace domain due to the distributed nature of the plant model. This has made possible the application of linear control theories to a set of linear-partial differential equations. The models obtained have well tracked the experimental results of the plant. A distributed-computer network has been interfaced with the plant to implement digital controllers in a hierarchical structure. A modern rnultivariable Wiener-Hopf controller has been applled to the plant model. The application has revealed a limitation condition that the plant matrix should be positive-definite along the infinite frequency axis. A new multi variable control theory has emerged fram this study, which avoids the above limitation. The controller has the structure of the modern Wiener-Hopf controller, but with a unique feature enabling a designer to specify the closed-loop poles in advance and to shape the sensitivity matrix as required. In this way, the method treats directly the interaction problems found in the chemical processes with good tracking and regulation performances. Though the ability of the analytical design methods to determine once and for all whether a given set of specifications can be met is one of its chief advantages over the conventional trial-and-error design procedures. However, one disadvantage that offsets to some degree the enormous advantages is the relatively complicated algebra that must be employed in working out all but the simplest problem. Mathematical algorithms and computer software have been developed to treat some of the mathematical operations defined over the integral domain, such as matrix fraction description, spectral factorization, the Bezout identity, and the general manipulation of polynomial matrices. Hence, the design problems of Wiener-Hopf type of controllers and other similar algebraic design methods can be easily solved.

Cognition, ocular accommodation, and cardiovascular function in emmetropes and late-onset myopes

PURPOSE. To investigate objectively and noninvasively the role of cognitive demand on autonomic control of systemic cardiovascular and ocular accommodative responses in emmetropes and myopes of late-onset. METHODS. Sixteen subjects (10 men, 6 women) aged between 18 and 34 years (mean ± SD: 22.6 ± 4.4 years), eight emmetropes (EMMs; mean spherical equivalent [MSE] refractive error ± SD: 0.05 ± 0.24 D) and eight with late-onset myopia (LOMs; MSE ± SD: -3.66 ± 2.31 D) participated in the study. Subjects viewed stationary numerical digits monocularly within a Badal optical system (at both 0.0 and -3.0 D) while performing a two-alternative, forced-choice paradigm that matched cognitive loading across subjects. Five individually matched cognitive levels of increasing difficulty were used in random order for each subject. Five 20-second, continuous-objective recordings of the accommodative response measured with an open-view infrared autorefractor were obtained for each cognitive level, whereas simultaneous measurement of heart rate was continuously recorded with a finger-mounted piezoelectric pulse transducer for 5 minutes. Fast Fourier transformation of cardiovascular function allowed the relative power of the autonomic components to be assessed in the frequency domain, whereas heart period gave an indication of the time-domain response. RESULTS. Increasing the cognitive demand led to a significant reduction in the accommodative response in all subjects (0.0 D: by -0.35 ± 0.33 D; -3.0 D: by -0.31 ± 0.40 D, P < 0.001). The greater lag of LOMs compared with EMMs was not significant (P = 0.07) at both distance (0.38 ± 0.35 D) and near (0.14 ± 0.42 D). Mean heart period reduced with increasing levels of workload (P < 0.0005). LOMs exhibited a relative elevation in sympathetic system activity compared to EMMs. Within refractive groups, however, accommodative shifts with increasing cognition correlated with parasympathetic activity (r = 0.99, P < 0.001), more than with sympathetic activity (r = 0.62, P > 0.05). CONCLUSIONS. In an equivalent workload paradigm, increasing cognitive demand caused a reduction in accommodative response that was attributable principally to a concurrent reduction in the relative power of the parasympathetic component of the autonomic nervous system (ANS). The disparity in accommodative response between EMMs and LOMs, however, appears to be augmented by changes in the sympathetic nervous component of the systemic ANS. Copyright © Association for Research in Vision and Ophthalmology.

Autonomic correlates of ocular accommodation and cardiovascular function

PURPOSE: To evaluate the hypothesis that objective measures of open- and closed-loop ocular accommodation are related to systemic cardiovascular function, and ipso facto autonomic nervous system activity. METHODS: Sixty subjects (29 male; 31 female) varying in age from 18 to 33 years (average: 20.3 +/- 2.9 years) with a range of refractive errors [mean spherical equivalent (MSE): -7.12 to +1.82 D] participated in the study. Five 20-s continuous objective recordings of the accommodative response, measured with an open-view IR autorefractor (Shin-Nippon SRW-5000), were obtained for a variety of open- and closed-loop accommodative demands while simultaneous continuous measurement of heart rate was recorded with a finger-mounted piezo-electric pulse transducer for 5 min. Fast Fourier Transformation of cardiovascular function allowed the absolute and relative power of the autonomic components to be assessed in the frequency-domain, whereas heart period gave an indication of the time-domain response. RESULTS: Increasing closed-loop accommodative demand led to a concurrent increase in heart rate of approximately 2 beats/min for a 4.0 D increase in accommodative demand. The increase was attributable to a reduction in the absolute (p < 0.05) and normalised (p < 0.001) input of the systemic parasympathetic nervous system, and was unaffected by refractive group. The interaction with refractive group failed to reach significance. CONCLUSIONS: For sustained accommodation effort, the data demonstrate covariation between the oculomotor and cardiovascular systems which implies that a near visual task can significantly influence cardiovascular behaviour. Accommodative effort alone, however, is not a sufficient stimulus to induce autonomic differences between refractive groups. The data suggest that both the oculomotor and cardiovascular systems are predominantly attributable to changes in the systemic parasympathetic nervous system.

Temporal finite elements in the dynamic analysis of structures

This thesis demonstrates that the use of finite elements need not be confined to space alone, but that they may also be used in the time domain, It is shown that finite element methods may be used successfully to obtain the response of systems to applied forces, including, for example, the accelerations in a tall structure subjected to an earthquake shock. It is further demonstrated that at least one of these methods may be considered to be a practical alternative to more usual methods of solution. A detailed investigation of the accuracy and stability of finite element solutions is included, and methods of applications to both single- and multi-degree of freedom systems are described. Solutions using two different temporal finite elements are compared with those obtained by conventional methods, and a comparison of computation times for the different methods is given. The application of finite element methods to distributed systems is described, using both separate discretizations in space and time, and a combined space-time discretization. The inclusion of both viscous and hysteretic damping is shown to add little to the difficulty of the solution. Temporal finite elements are also seen to be of considerable interest when applied to non-linear systems, both when the system parameters are time-dependent and also when they are functions of displacement. Solutions are given for many different examples, and the computer programs used for the finite element methods are included in an Appendix.

Force sensing for measuring human body movement

The research developed in this thesis explores the sensing and inference of human movement in a dynamic way, as opposed to conventional measurement systems, that are only concerned with discrete evaluations of stimuli in sequential time. Typically, conventional approaches are used to infer the dynamic movement of the body; such as vision and motion tracking devices, with either a human diagnosis or complex image processing algorithm to classify the movement. This research is therefore the first of its kind to attempt and provide a movement classifying algorithm through the use of minimal sensing points, with the application for this novel system, to classify human movement during a golf swing. There are two main categories of force sensing. Firstly, array-type systems consisting of many sensing elements, and are the most commonly researched and commercially available. Secondly, reduced force sensing element systems (RFSES) also known as distributive systems have only been recently exploited in the academic world. The fundamental difference between these systems is that array systems handle the data captured from each sensor as unique outputs and suffer the effects of resolution. The effect of resolution, is the error in the load position measurement between sensing elements, as the output is quantized in terms of position. This can be compared to a reduced sensor element system that maximises that data received through the coupling of data from a distribution of sensing points to describe the output in discrete time. Also this can be extended to a coupling of transients in the time domain to describe an activity or dynamic movement. It is the RFSES that is to be examined and exploited in the commercial sector due to its advantages over array-based approaches such as reduced design, computational complexity and cost.

The application of parallel processing methods to vibrational analysis

The trend in modal extraction algorithms is to use all the available frequency response functions data to obtain a global estimate of the natural frequencies, damping ratio and mode shapes. Improvements in transducer and signal processing technology allow the simultaneous measurement of many hundreds of channels of response data. The quantity of data available and the complexity of the extraction algorithms make considerable demands on the available computer power and require a powerful computer or dedicated workstation to perform satisfactorily. An alternative to waiting for faster sequential processors is to implement the algorithm in parallel, for example on a network of Transputers. Parallel architectures are a cost effective means of increasing computational power, and a larger number of response channels would simply require more processors. This thesis considers how two typical modal extraction algorithms, the Rational Fraction Polynomial method and the Ibrahim Time Domain method, may be implemented on a network of transputers. The Rational Fraction Polynomial Method is a well known and robust frequency domain 'curve fitting' algorithm. The Ibrahim Time Domain method is an efficient algorithm that 'curve fits' in the time domain. This thesis reviews the algorithms, considers the problems involved in a parallel implementation, and shows how they were implemented on a real Transputer network.

All-optical signal regeneration by temporal slicing of nonlinearly flattened optical waveform

A novel all-optical time domain regeneration technique using nonlinear pulse broadening and flattening in normal dispersion fiber and subsequent temporal slicing by an amplitude modulator (or a device performing a similar function) is proposed. Substantial suppression of the timing jitter of jitter-degraded optical signals is demonstrated using the proposed approach.

All-optical signal regeneration by temporal slicing of nonlinearly flattened optical waveform

A novel all-optical time domain regeneration technique using nonlinear pulse broadening and flattening in normal dispersion fiber and subsequent temporal slicing by an amplitude modulator (or a device performing a similar function) is proposed. Substantial suppression of the timing jitter of jitter-degraded optical signals is demonstrated using the proposed approach.

All-optical signal regeneration by temporal slicing of nonlinearly flattened optical waveform

A novel all-optical time domain regeneration technique using nonlinear pulse broadening and flattening in normal dispersion fiber and subsequent temporal slicing by an amplitude modulator (or a device performing a similar function) is proposed. Substantial suppression of the timing jitter of jitter-degraded optical signals is demonstrated using the proposed approach.

All-optical signal regeneration by temporal slicing of nonlinearly flattened optical waveform

A novel all-optical time domain regeneration technique using nonlinear pulse broadening and flattening in normal dispersion fiber and subsequent temporal slicing by an amplitude modulator (or a device performing a similar function) is proposed. Substantial suppression of the timing jitter of jitter-degraded optical signals is demonstrated using the proposed approach. © 2005 IEEE.

Full-vectorial modeling of femtosecond pulses for laser inscription of photonic structures

During the last decade, microfabrication of photonic devices by means of intense femtosecond (fs) laser pulses has emerged as a novel technology. A common requirement for the production of these devices is that the refractive index modification pitch size should be smaller than the inscribing wavelength. This can be achieved by making use of the nonlinear propagation of intense fs laser pulses. Nonlinear propagation of intense fs laser pulses is an extremely complicated phenomenon featuring complex multiscale spatiotemporal dynamics of the laser pulses. We have utilized a principal approach based on finite difference time domain (FDTD) modeling of the full set of Maxwell's equations coupled to the conventional Drude model for generated plasma. Nonlinear effects are included, such as self-phase modulation and multiphoton absorption. Such an approach resolves most problems related to the inscription of subwavelength structures, when the paraxial approximation is not applicable to correctly describe the creation of and scattering on the structures. In a representative simulation of the inscription process, the signature of degenerate four wave mixing has been found. © 2012 Optical Society of America.