Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks

This paper studies the fracturing process in low-porous rocks during uniaxial compressive tests considering the original defects and the new mechanical cracks in the material. For this purpose, five different kinds of rocks have been chosen with carbonate mineralogy and low porosity (lower than 2%). The characterization of the fracture damage is carried out using three different techniques: ultrasounds, mercury porosimetry and X-ray computed tomography. The proposed methodology allows quantifying the evolution of the porous system as well as studying the location of new cracks in the rock samples. Intercrystalline porosity (the smallest pores with pore radius < 1 μm) shows a limited development during loading, disappearing rapidly from the porosimetry curves and it is directly related to the initial plastic behaviour in the stress–strain patterns. However, the biggest pores (corresponding to the cracks) suffer a continuous enlargement until the unstable propagation of fractures. The measured crack initiation stress varies between 0.25 σp and 0.50 σp for marbles and between 0.50 σp and 0.85 σp for micrite limestone. The unstable propagation of cracks is assumed to occur very close to the peak strength. Crack propagation through the sample is completely independent of pre-existing defects (porous bands, stylolites, fractures and veins). The ultrasonic response in the time-domain is less sensitive to the fracture damage than the frequency-domain. P-wave velocity increases during loading test until the beginning of the unstable crack propagation. This increase is higher for marbles (between 15% and 30% from initial vp values) and lower for micrite limestones (between 5% and 10%). When the mechanical cracks propagate unstably, the velocity stops to increase and decreases only when rock damage is very high. Frequency analysis of the ultrasonic signals shows clear changes during the loading process. The spectrum of treated waveforms shows two main frequency peaks centred at low (~ 20 kHz) and high (~ 35 kHz) values. When new fractures appear and grow the amplitude of the high-frequency peak decreases, while that of the low-frequency peak increases. Besides, a slight frequency shift is observed towards higher frequencies.

Parameters of the Earths orbit for the last 5 Million years in 1 kyr resolution

New values for the astronomical parameters of the Earth's orbit and rotation (eccentricity, obliquity and precession) are proposed for paleoclimatic research related to the Late Miocene, the Pliocene and the Quaternary. They have been obtained from a numerical solution of the Lagrangian system of the planetary point masses and from an analytical solution of the Poisson equations of the Earth-Moon system. The analytical expansion developed in this paper allows the direct determination of the main frequencies with their phase and amplitude. Numerical and analytical comparisons with the former astronomical solution BER78 are performed so that the accuracy and the interval of time over which the new solution is valid can be estimated. The corresponding insolation values have also been computed and compared to the former ones. This analysis leads to the conclusion that the new values are expected to be reliable over the last 5 Ma in the time domain and at least over the last 10 Ma in the frequency domain.

Experimental muscle pain changes feedforward postural responses of the trunk muscles

Many studies have identified changes in trunk muscle recruitment in clinical low back pain (LBP). However, due to the heterogeneity of the LBP population these changes have been variable and it has been impossible to identify a cause-effect relationship. Several studies have identified a consistent change in the feed-forward postural response of transversus abdominis (TrA), the deepest abdominal muscle, in association with arm movements in chronic LBP. This study aimed to determine whether the feedforward recruitment of the trunk muscles in a postural task could be altered by acute experimentally induced LBP. Electromyographic (EMG) recordings of the abdominal and paraspinal muscles were made during arm movements in a control trial, following the injection of isotonic (non-painful) and hypertonic (painful) saline into the longissimus muscle at L4, and during a 1-h follow-up. Movements included rapid arm flexion in response to a light and repetitive arm flexion-extension. Temporal and spatial EMG parameters were measured. The onset and amplitude of EMG of most muscles was changed in a variable manner during the period of experimentally induced pain. However, across movement trials and subjects the activation of TrA was consistently reduced in amplitude or delayed. Analyses in the time and frequency domain were used to confirm these findings. The results suggest that acute experimentally induced pain may affect feedforward postural activity of the trunk muscles. Although the response was variable, pain produced differential changes in the motor control of the trunk muscles, with consistent impairment of TrA activity.

Fluctuations and single molecules

The fluorescence of single molecules coupled to a thermal bath is studied both experimentally and theoretically. The effect of different fluctuations on the coherence properties of resonance fluorescence is considered first. Coherence is measured in an interference experiment where a single molecule is used as a light source. A standard approach based on the optical Bloch equations apparently provides quite an accurate description of the interference experiment. Systems with long correlation times (where spectra are time dependent on any timescale) are considered next. It is shown that intensity-time-frequency correlation spectroscopy, which provides both high signal-to-noise ratio and high time resolution, is very suitable for such a case. The Bloch equations are further tested in an experiment where the shape of an excitation spectral line of a single molecule is accurately measured over six orders of magnitude of the exciting laser power. Significant deviations from the predictions of the Bloch equations are found. The role of critical parameters-the correlation time of the bath, the Rabi oscillation period, and the coupling constant between the bath and the molecule-is discussed. The paper also includes a short general introduction to the methodology of single-molecule studies.

Rainfall variability at decadal and longer time scales: signal or noise?

Rainfall variability occurs over a wide range of temporal scales. Knowledge and understanding of such variability can lead to improved risk management practices in agricultural and other industries. Analyses of temporal patterns in 100 yr of observed monthly global sea surface temperature and sea level pressure data show that the single most important cause of explainable, terrestrial rainfall variability resides within the El Nino-Southern Oscillation (ENSO) frequency domain (2.5-8.0 yr), followed by a slightly weaker but highly significant decadal signal (9-13 yr), with some evidence of lesser but significant rainfall variability at interclecadal time scales (15-18 yr). Most of the rainfall variability significantly linked to frequencies tower than ENSO occurs in the Australasian region, with smaller effects in North and South America, central and southern Africa, and western Europe. While low-frequency (LF) signals at a decadal frequency are dominant, the variability evident was ENSO-like in all the frequency domains considered. The extent to which such LF variability is (i) predictable and (ii) either part of the overall ENSO variability or caused by independent processes remains an as yet unanswered question. Further progress can only be made through mechanistic studies using a variety of models.

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.

Sviluppo e applicazione di nuovi criteri per la verifica a fatica multiassiale di componenti strutturali in leghe di acciaio e alluminio

La presente Tesi ha per oggetto lo sviluppo e la validazione di nuovi criteri per la verifica a fatica multiassiale di componenti strutturali metallici . In particolare, i nuovi criteri formulati risultano applicabili a componenti metallici, soggetti ad un’ampia gamma di configurazioni di carico: carichi multiassiali variabili nel tempo, in modo ciclico e random, per alto e basso/medio numero di cicli di carico. Tali criteri costituiscono un utile strumento nell’ambito della valutazione della resistenza/vita a fatica di elementi strutturali metallici, essendo di semplice implementazione, e richiedendo tempi di calcolo piuttosto modesti. Nel primo Capitolo vengono presentate le problematiche relative alla fatica multiassiale, introducendo alcuni aspetti teorici utili a descrivere il meccanismo di danneggiamento a fatica (propagazione della fessura e frattura finale) di componenti strutturali metallici soggetti a carichi variabili nel tempo. Vengono poi presentati i diversi approcci disponibili in letteratura per la verifica a fatica multiassiale di tali componenti, con particolare attenzione all'approccio del piano critico. Infine, vengono definite le grandezze ingegneristiche correlate al piano critico, utilizzate nella progettazione a fatica in presenza di carichi multiassiali ciclici per alto e basso/medio numero di cicli di carico. Il secondo Capitolo è dedicato allo sviluppo di un nuovo criterio per la valutazione della resistenza a fatica di elementi strutturali metallici soggetti a carichi multiassiali ciclici e alto numero di cicli. Il criterio risulta basato sull'approccio del piano critico ed è formulato in termini di tensioni. Lo sviluppo del criterio viene affrontato intervenendo in modo significativo su una precedente formulazione proposta da Carpinteri e collaboratori nel 2011. In particolare, il primo intervento riguarda la determinazione della giacitura del piano critico: nuove espressioni dell'angolo che lega la giacitura del piano critico a quella del piano di frattura vengono implementate nell'algoritmo del criterio. Il secondo intervento è relativo alla definizione dell'ampiezza della tensione tangenziale e un nuovo metodo, noto come Prismatic Hull (PH) method (di Araújo e collaboratori), viene implementato nell'algoritmo. L'affidabilità del criterio viene poi verificata impiegando numerosi dati di prove sperimentali disponibili in letteratura. Nel terzo Capitolo viene proposto un criterio di nuova formulazione per la valutazione della vita a fatica di elementi strutturali metallici soggetti a carichi multiassiali ciclici e basso/medio numero di cicli. Il criterio risulta basato sull'approccio del piano critico, ed è formulato in termini di deformazioni. In particolare, la formulazione proposta trae spunto, come impostazione generale, dal criterio di fatica multiassiale in regime di alto numero di cicli discusso nel secondo Capitolo. Poiché in presenza di deformazioni plastiche significative (come quelle caratterizzanti la fatica per basso/medio numero di cicli di carico) è necessario conoscere il valore del coefficiente efficace di Poisson del materiale, vengono impiegate tre differenti strategie. In particolare, tale coefficiente viene calcolato sia per via analitica, che per via numerica, che impiegando un valore costante frequentemente adottato in letteratura. Successivamente, per validarne l'affidabilità vengono impiegati numerosi dati di prove sperimentali disponibili in letteratura; i risultati numerici sono ottenuti al variare del valore del coefficiente efficace di Poisson. Inoltre, al fine di considerare i significativi gradienti tensionali che si verificano in presenza di discontinuità geometriche, come gli intagli, il criterio viene anche esteso al caso dei componenti strutturali intagliati. Il criterio, riformulato implementando il concetto del volume di controllo proposto da Lazzarin e collaboratori, viene utilizzato per stimare la vita a fatica di provini con un severo intaglio a V, realizzati in lega di titanio grado 5. Il quarto Capitolo è rivolto allo sviluppo di un nuovo criterio per la valutazione del danno a fatica di elementi strutturali metallici soggetti a carichi multiassiali random e alto numero di cicli. Il criterio risulta basato sull'approccio del piano critico ed è formulato nel dominio della frequenza. Lo sviluppo del criterio viene affrontato intervenendo in modo significativo su una precedente formulazione proposta da Carpinteri e collaboratori nel 2014. In particolare, l’intervento riguarda la determinazione della giacitura del piano critico, e nuove espressioni dell'angolo che lega la giacitura del piano critico con quella del piano di frattura vengono implementate nell'algoritmo del criterio. Infine, l’affidabilità del criterio viene verificata impiegando numerosi dati di prove sperimentali disponibili in letteratura.

O EFEITO DO CONTROLE RESPIRATÓRIO EM VARIÁVEIS ELETROFISIOLÓGICAS DA ATENÇÃO

A prática do ioga tem se tornado cada vez mais popular, não apenas pelos benefícios físicos, mas principalmente pelo bem-estar psicológico trazido pela sua prática. Um dos componentes do ioga é o Prãnãyama, ou controle da respiração. A atenção e a respiração são dois mecanismos fisiológicos e involuntários requeridos para a execução do Prãnãyama. O principal objetivo desse estudo foi verificar se variáveis contínuas do EEG (potência de diferentes faixas que o compõem) seriam moduladas pelo controle respiratório, comparando-se separadamente as duas fases do ciclo respiratório (inspiração e expiração), na situação de respiração espontânea e controlada. Fizeram parte do estudo 19 sujeitos (7 homens/12 mulheres, idade média de 36,89 e DP = ± 14,46) que foram convidados a participar da pesquisa nas dependências da Faculdade de Saúde da Universidade Metodista de São Paulo. Para o registro do eletroencefalograma foi utilizado um sistema de posicionamento de cinco eletrodos Ag AgCl (FPz, Fz, Cz, Pz e Oz) fixados a uma touca de posicionamento rápido (Quick-Cap, Neuromedical Supplies®), em sistema 10-20. Foram obtidos valores de máxima amplitude de potência (espectro de potência no domínio da frequência) nas frequências teta, alfa e beta e delta e calculada a razão teta/beta nas diferentes fases do ciclo respiratório (inspiração e expiração), separadamente, nas condições de respiração espontânea e de controle respiratório. Para o registro do ciclo respiratório, foi utilizada uma cinta de esforço respiratório M01 (Pletismógrafo). Os resultados mostram diferenças significativas entre as condições de respiração espontânea e de controle com valores das médias da razão teta/beta menores na respiração controlada do que na respiração espontânea e valores de média da potência alfa sempre maiores no controle respiratório. Diferenças significativas foram encontradas na comparação entre inspiração e expiração da respiração controlada com diminuição dos valores das médias da razão teta/beta na inspiração e aumento nos valores das médias da potência alfa, sobretudo na expiração. Os achados deste estudo trazem evidências de que o controle respiratório modula variáveis eletrofisiológicas relativas à atenção refletindo um estado de alerta, porém mais relaxado do que na situação de respiração espontânea.

Autonomic dysfunction in unselected and untreated primary open angle glaucoma patients:A pilot study

Purpose: To investigate the presence of silent cardiac ischaemic episodes and the status of autonomic function in consecutive, newly diagnosed and untreated primary open-angle glaucoma patients. Methods: Twenty-four consecutively diagnosed glaucoma patients and 22 age-matched controls were subjected to ambulatory 24-h blood pressure (BP) and electrocardiogram (ECG) monitoring by using Cardiotens-01 (Meditech Ltd). Based on the ECG recordings, heart rate variability (HRV) frequency domain parameters [low-frequency (LF), high-frequency (HF) and LF/HF ratio] were calculated and analysed in the two study groups. Results: Glaucoma patients demonstrated higher LF and LF/HF values than normal subjects for both the active period (p = 0.020 and 0.029) and the passive period (p = 0.044 and 0.049 respectively). HRV parameters were similar in patients and controls suffering from silent cardiac ischaemia (p > 0.05); however, glaucoma patients with normal ECG demonstrated higher LF and LF/HF values during the active period of the 24-h measurement period than control subjects characterized by the same cardiac activity (p = 0.010 and 0.021 respectively). Conclusion: Independent of a history and/or clinical signs of cardiovascular disease, glaucoma patients exhibit abnormal autonomic function. © 2007 The Authors.

Writing and characterisation methods for UV-photosensitivity based devices and their applications

The underlying work to this thesis focused on the exploitation and investigation of photosensitivity mechanisms in optical fibres and planar waveguides for the fabrication of advanced integrated optical devices for telecoms and sensing applications. One major scope is the improvement of grating fabrication specifications by introducing new writing techniques and the use of advanced characterisation methods for grating testing. For the first time the polarisation control method for advanced grating fabrication has successfully been converted to apodised planar waveguide fabrication and the development of a holographic method for the inscription of chirped gratings at arbitrary wavelength is presented. The latter resulted in the fabrication of gratings for pulse-width suppression and wavelength selection in diode lasers. In co-operation with research partners a number of samples were tested using optical frequency domain and optical low coherence reflectometry for a better insight into the limitations of grating writing techniques. Using a variety of different fabrication methods, custom apodised and chirped fibre Bragg gratings were written for the use as filter elements for multiplexer-demultiplexer devices, as well as for short pulse generation and wavelength selection in telecommunication transmission systems. Long period grating based devices in standard, speciality and tapered fibres are presented, showing great potential for multi-parameter sensing. One particular scope is the development of vectorial curvature and refractive index sensors with potential for medical, chemical and biological sensing. In addition the design of an optically tunable Mach-Zehnder based multiwavelength filter is introduced. The discovery of a Type IA grating type through overexposure of hydrogen loaded standard and Boron-Germanium co-doped fibres strengthened the assumption of UV-photosensitivity being a highly non-linear process. Gratings of this type show a significantly lower thermal sensitivity compared to standard gratings, which makes them useful for sensing applications. An Oxford Lasers copper-vapour laser operating at 255 nm in pulsed mode was used for their inscription, in contrast to previous work using CW-Argon-Ion lasers and contributing to differences in the processes of the photorefractive index change

Electromechanical resonators as measuring transducers

With the increasing use of digital computers for data acquisition and digital process control, frequency domain transducers have become very attractive due to their virtual digital output. Essentially they are electrically maintained oscillators where the sensor is the controlling resonator.They are designed to make the frequency a function of the physical parameter being measured. Because of their high quality factor, mechanical resonators give very good frequency stability and are widely used as sensors. For this work symmetrical mechanical resonators such as the tuning fork were considered, to be the most promising. These are dynamically clamped and can be designed to have extensive regions where no vibrations occur.This enables the resonators to be robustly mounted in a way convenient for various applications. Designs for the measurement of fluid density and tension have been produced. The principle of the design of the resonator for fluid density measurement is a thin gap (trapping a lamina of fluid) between its two members which vibrate in antiphase.An analysis of the inter­ action between this resonator and the fluid lamina has carried out.In gases narrow gaps are needed for a good sensitivity and the use of the material fused quartz, because of its low density and very low temperature coefficient, is ideally suitable. In liquids an adequate sensitivity is achieved even with a wide lamina gap. Practical designs of such transducers have been evolved. The accuracy for liquid measurements is better than 1%. For gases it was found that, in air, a change of atmospheric pressure of 0.3% could be detected. In constructing a tension transducer using such a mechanical sensor as a wire or a beam, major difficulties are encountered in making an efficient clamping arrangement for the sensor. The use of dynamically clamped beams has been found to overcome the problem and this is the basis of the transducer investigated.

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.

Parallel implementation of image segmentation algorithms on a network of transputers

Image segmentation is one of the most computationally intensive operations in image processing and computer vision. This is because a large volume of data is involved and many different features have to be extracted from the image data. This thesis is concerned with the investigation of practical issues related to the implementation of several classes of image segmentation algorithms on parallel architectures. The Transputer is used as the basic building block of hardware architectures and Occam is used as the programming language. The segmentation methods chosen for implementation are convolution, for edge-based segmentation; the Split and Merge algorithm for segmenting non-textured regions; and the Granlund method for segmentation of textured images. Three different convolution methods have been implemented. The direct method of convolution, carried out in the spatial domain, uses the array architecture. The other two methods, based on convolution in the frequency domain, require the use of the two-dimensional Fourier transform. Parallel implementations of two different Fast Fourier Transform algorithms have been developed, incorporating original solutions. For the Row-Column method the array architecture has been adopted, and for the Vector-Radix method, the pyramid architecture. The texture segmentation algorithm, for which a system-level design is given, demonstrates a further application of the Vector-Radix Fourier transform. A novel concurrent version of the quad-tree based Split and Merge algorithm has been implemented on the pyramid architecture. The performance of the developed parallel implementations is analysed. Many of the obtained speed-up and efficiency measures show values close to their respective theoretical maxima. Where appropriate comparisons are drawn between different implementations. The thesis concludes with comments on general issues related to the use of the Transputer system as a development tool for image processing applications; and on the issues related to the engineering of concurrent image processing applications.

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.