Motion correction algorithm of lung tumors for respiratory gated PET images

Respiratory gating in lung PET imaging to compensate for respiratory motion artifacts is a current research issue with broad potential impact on quantitation, diagnosis and clinical management of lung tumors. However, PET images collected at discrete bins can be significantly affected by noise as there are lower activity counts in each gated bin unless the total PET acquisition time is prolonged, so that gating methods should be combined with imaging-based motion correction and registration methods. The aim of this study was to develop and validate a fast and practical solution to the problem of respiratory motion for the detection and accurate quantitation of lung tumors in PET images. This included: (1) developing a computer-assisted algorithm for PET/CT images that automatically segments lung regions in CT images, identifies and localizes lung tumors of PET images; (2) developing and comparing different registration algorithms which processes all the information within the entire respiratory cycle and integrate all the tumor in different gated bins into a single reference bin. Four registration/integration algorithms: Centroid Based, Intensity Based, Rigid Body and Optical Flow registration were compared as well as two registration schemes: Direct Scheme and Successive Scheme. Validation was demonstrated by conducting experiments with the computerized 4D NCAT phantom and with a dynamic lung-chest phantom imaged using a GE PET/CT System. Iterations were conducted on different size simulated tumors and different noise levels. Static tumors without respiratory motion were used as gold standard; quantitative results were compared with respect to tumor activity concentration, cross-correlation coefficient, relative noise level and computation time. Comparing the results of the tumors before and after correction, the tumor activity values and tumor volumes were closer to the static tumors (gold standard). Higher correlation values and lower noise were also achieved after applying the correction algorithms. With this method the compromise between short PET scan time and reduced image noise can be achieved, while quantification and clinical analysis become fast and precise.

Mitigating congestion by integrating time forecasting and realtime information aggregation in cellular networks

An iterative travel time forecasting scheme, named the Advanced Multilane Prediction based Real-time Fastest Path (AMPRFP) algorithm, is presented in this dissertation. This scheme is derived from the conventional kernel estimator based prediction model by the association of real-time nonlinear impacts that caused by neighboring arcs’ traffic patterns with the historical traffic behaviors. The AMPRFP algorithm is evaluated by prediction of the travel time of congested arcs in the urban area of Jacksonville City. Experiment results illustrate that the proposed scheme is able to significantly reduce both the relative mean error (RME) and the root-mean-squared error (RMSE) of the predicted travel time. To obtain high quality real-time traffic information, which is essential to the performance of the AMPRFP algorithm, a data clean scheme enhanced empirical learning (DCSEEL) algorithm is also introduced. This novel method investigates the correlation between distance and direction in the geometrical map, which is not considered in existing fingerprint localization methods. Specifically, empirical learning methods are applied to minimize the error that exists in the estimated distance. A direction filter is developed to clean joints that have negative influence to the localization accuracy. Synthetic experiments in urban, suburban and rural environments are designed to evaluate the performance of DCSEEL algorithm in determining the cellular probe’s position. The results show that the cellular probe’s localization accuracy can be notably improved by the DCSEEL algorithm. Additionally, a new fast correlation technique for overcoming the time efficiency problem of the existing correlation algorithm based floating car data (FCD) technique is developed. The matching process is transformed into a 1-dimensional (1-D) curve matching problem and the Fast Normalized Cross-Correlation (FNCC) algorithm is introduced to supersede the Pearson product Moment Correlation Co-efficient (PMCC) algorithm in order to achieve the real-time requirement of the FCD method. The fast correlation technique shows a significant improvement in reducing the computational cost without affecting the accuracy of the matching process.

Peak -to -average power ratio reduced parallel interference cancellation Multicarrier-Code Division Multiple Access system with anti -interference property

Orthogonal Frequency-Division Multiplexing (OFDM) has been proved to be a promising technology that enables the transmission of higher data rate. Multicarrier Code-Division Multiple Access (MC-CDMA) is a transmission technique which combines the advantages of both OFDM and Code-Division Multiplexing Access (CDMA), so as to allow high transmission rates over severe time-dispersive multi-path channels without the need of a complex receiver implementation. Also MC-CDMA exploits frequency diversity via the different subcarriers, and therefore allows the high code rates systems to achieve good Bit Error Rate (BER) performances. Furthermore, the spreading in the frequency domain makes the time synchronization requirement much lower than traditional direct sequence CDMA schemes. There are still some problems when we use MC-CDMA. One is the high Peak-to-Average Power Ratio (PAPR) of the transmit signal. High PAPR leads to nonlinear distortion of the amplifier and results in inter-carrier self-interference plus out-of-band radiation. On the other hand, suppressing the Multiple Access Interference (MAI) is another crucial problem in the MC-CDMA system. Imperfect cross-correlation characteristics of the spreading codes and the multipath fading destroy the orthogonality among the users, and then cause MAI, which produces serious BER degradation in the system. Moreover, in uplink system the received signals at a base station are always asynchronous. This also destroys the orthogonality among the users, and hence, generates MAI which degrades the system performance. Besides those two problems, the interference should always be considered seriously for any communication system. In this dissertation, we design a novel MC-CDMA system, which has low PAPR and mitigated MAI. The new Semi-blind channel estimation and multi-user data detection based on Parallel Interference Cancellation (PIC) have been applied in the system. The Low Density Parity Codes (LDPC) has also been introduced into the system to improve the performance. Different interference models are analyzed in multi-carrier communication systems and then the effective interference suppression for MC-CDMA systems is employed in this dissertation. The experimental results indicate that our system not only significantly reduces the PAPR and MAI but also effectively suppresses the outside interference with low complexity. Finally, we present a practical cognitive application of the proposed system over the software defined radio platform.

Peak-to-Average Power Ratio Reduced Parallel Interference Cancellation Multicarrier-Code Division Multiple Access System with Anti-Interference Property

Orthogonal Frequency-Division Multiplexing (OFDM) has been proved to be a promising technology that enables the transmission of higher data rate. Multicarrier Code-Division Multiple Access (MC-CDMA) is a transmission technique which combines the advantages of both OFDM and Code-Division Multiplexing Access (CDMA), so as to allow high transmission rates over severe time-dispersive multi-path channels without the need of a complex receiver implementation. Also MC-CDMA exploits frequency diversity via the different subcarriers, and therefore allows the high code rates systems to achieve good Bit Error Rate (BER) performances. Furthermore, the spreading in the frequency domain makes the time synchronization requirement much lower than traditional direct sequence CDMA schemes. There are still some problems when we use MC-CDMA. One is the high Peak-to-Average Power Ratio (PAPR) of the transmit signal. High PAPR leads to nonlinear distortion of the amplifier and results in inter-carrier self-interference plus out-of-band radiation. On the other hand, suppressing the Multiple Access Interference (MAI) is another crucial problem in the MC-CDMA system. Imperfect cross-correlation characteristics of the spreading codes and the multipath fading destroy the orthogonality among the users, and then cause MAI, which produces serious BER degradation in the system. Moreover, in uplink system the received signals at a base station are always asynchronous. This also destroys the orthogonality among the users, and hence, generates MAI which degrades the system performance. Besides those two problems, the interference should always be considered seriously for any communication system. In this dissertation, we design a novel MC-CDMA system, which has low PAPR and mitigated MAI. The new Semi-blind channel estimation and multi-user data detection based on Parallel Interference Cancellation (PIC) have been applied in the system. The Low Density Parity Codes (LDPC) has also been introduced into the system to improve the performance. Different interference models are analyzed in multi-carrier communication systems and then the effective interference suppression for MC-CDMA systems is employed in this dissertation. The experimental results indicate that our system not only significantly reduces the PAPR and MAI but also effectively suppresses the outside interference with low complexity. Finally, we present a practical cognitive application of the proposed system over the software defined radio platform.

Motion Correction Algorithm of Lung Tumors for Respiratory Gated PET Images

Respiratory gating in lung PET imaging to compensate for respiratory motion artifacts is a current research issue with broad potential impact on quantitation, diagnosis and clinical management of lung tumors. However, PET images collected at discrete bins can be significantly affected by noise as there are lower activity counts in each gated bin unless the total PET acquisition time is prolonged, so that gating methods should be combined with imaging-based motion correction and registration methods. The aim of this study was to develop and validate a fast and practical solution to the problem of respiratory motion for the detection and accurate quantitation of lung tumors in PET images. This included: (1) developing a computer-assisted algorithm for PET/CT images that automatically segments lung regions in CT images, identifies and localizes lung tumors of PET images; (2) developing and comparing different registration algorithms which processes all the information within the entire respiratory cycle and integrate all the tumor in different gated bins into a single reference bin. Four registration/integration algorithms: Centroid Based, Intensity Based, Rigid Body and Optical Flow registration were compared as well as two registration schemes: Direct Scheme and Successive Scheme. Validation was demonstrated by conducting experiments with the computerized 4D NCAT phantom and with a dynamic lung-chest phantom imaged using a GE PET/CT System. Iterations were conducted on different size simulated tumors and different noise levels. Static tumors without respiratory motion were used as gold standard; quantitative results were compared with respect to tumor activity concentration, cross-correlation coefficient, relative noise level and computation time. Comparing the results of the tumors before and after correction, the tumor activity values and tumor volumes were closer to the static tumors (gold standard). Higher correlation values and lower noise were also achieved after applying the correction algorithms. With this method the compromise between short PET scan time and reduced image noise can be achieved, while quantification and clinical analysis become fast and precise.

Estudo da estrutura subsuperficial da província Borborema com correlação de ruído sísmico

Ambient seismic noise has traditionally been considered as an unwanted perturbation in seismic data acquisition that "contaminates" the clean recording of earthquakes. Over the last decade, however, it has been demonstrated that consistent information about the subsurface structure can be extracted from cross-correlation of ambient seismic noise. In this context, the rules are reversed: the ambient seismic noise becomes the desired seismic signal, while earthquakes become the unwanted perturbation that needs to be removed. At periods lower than 30 s, the spectrum of ambient seismic noise is dominated by microseism, which originates from distant atmospheric perturbations over the oceans. The microsseism is the most continuous seismic signal and can be classified as primary – when observed in the range 10-20 s – and secondary – when observed in the range 5-10 s. The Green‘s function of the propagating medium between two receivers (seismic stations) can be reconstructed by cross-correlating seismic noise simultaneously recorded at the receivers. The reconstruction of the Green‘s function is generally proportional to the surface-wave portion of the seismic wavefield, as microsseismic energy travels mostly as surface-waves. In this work, 194 Green‘s functions obtained from stacking of one month of daily cross-correlations of ambient seismic noise recorded in the vertical component of several pairs of broadband seismic stations in Northeast Brazil are presented. The daily cross-correlations were stacked using a timefrequency, phase-weighted scheme that enhances weak coherent signals by reducing incoherent noise. The cross-correlations show that, as expected, the emerged signal is dominated by Rayleigh waves, with dispersion velocities being reliably measured for periods ranging between 5 and 20 s. Both permanent stations from a monitoring seismic network and temporary stations from past passive experiments in the region are considered, resulting in a combined network of 33 stations separated by distances between 60 and 1311 km, approximately. The Rayleigh-wave, dispersion velocity measurements are then used to develop tomographic images of group velocity variation for the Borborema Province of Northeast Brazil. The tomographic maps allow to satisfactorily map buried structural features in the region. At short periods (~5 s) the images reflect shallow crustal structure, clearly delineating intra-continental and marginal sedimentary basins, as well as portions of important shear zones traversing the Borborema Province. At longer periods (10 – 20 s) the images are sensitive to deeper structure in the upper crust, and most of the shallower anomalies fade away. Interestingly, some of them do persist. The deep anomalies do not correlate with either the location of Cenozoic volcanism and uplift - which marked the evolution of the Borborema Province in the Cenozoic - or available maps of surface heat-flow, and the origin of the deep anomalies remains enigmatic.

Caracterização do ritmo circadiano de atividade e repouso em saguis idosos (Callithrix Jacchus)

Advanced age may become a limiting factor for the maintenance of rhythms in organisms, reducing the capacity of generation and synchronization of biological rhythms. In this study, the influence of aging on the expression of endogenous periodicity and synchronization (photic and social) of the circadian activity rhythm (CAR) was evaluated in a diurnal primate, the marmoset (Callithrix jacchus). This study had two approaches: one with longitudinal design, performed with a male marmoset in two different phases: adult (three years) and older (9 y.o.) (study 1) and the second, a transversal approach, with 6 old (♂: 9.7 ± 2.0 y.o.) and 11 adults animals (♂: 4.2 ± 0.8 y.o.) (study 2). The evaluation of the photic synchronization involved two conditions in LD (natural and artificial illuminations). In study 1, the animal was subjected to the following stages: LD (12:12 ~ 350: ~ 2 lx), LL (~ 350 lx) and LD resynchronization. In the second study, the animals were initially evaluated in natural LD, and then the same sequence stages of study 1. During the LL stage in study 2, the vocalizations of conspecifics kept in natural LD on the outside of the colony were considered temporal cue to the social synchronization. The record of the activity was performed automatically at intervals of five minutes through infrared sensor and actimeters, in studies 1 and 2, respectively. In general, the aged showed a more fragmented activity pattern (> IV < H and > PSD, ANOVA, p < 0.05), lower levels of activity (ANOVA, p < 0.05) and shorter duration of active phase (ANOVA, p < 0.05) in LD conditions, when compared to adults. In natural LD, the aged presented phase delay pronounced for onset and offset of active phase (ANOVA, p < 0.05), while the adults had the active phase more adjusted to light phase. Under artificial LD, there was phase advance and greater adjustment of onset and offset of activity in relation to the LD in the aged (ANOVA, p < 0.05). In LL, there was a positive correlation between age and the endogenous period () in the first 20 days (Spearman correlation, p < 0.05), with prolonged  held in two aged animals. In this condition, most adults showed free-running period of the circadian activity rhythm with  < 24 h for the first 30 days and later on relative coordination mediated by auditory cues. In study 2, the cross-correlation analysis between the activity profiles of the animals in LL with control animals kept under natural LD, found that there was less social synchronization in the aged. With the resubmission to the LD, the resynchronization rate was slower in the aged (t-test; p < 0.05) and in just one aged animal there was a loss of resynchronization capability. According to the data set, it is suggested that the aging in marmosets may be related to: 1) lower amplitude and greater fragmentation of the activity, accompanied to phase delay with extension of period, caused by changes in a photic input, in the generation and behavioral expression of the CAR; 2) lower capacity of the circadian activity rhythm to photic synchronization, that can become more robust in artificial lighting conditions, possibly due to the higher light intensities at the beginning of the active phase due to the abrupt transitions between the light and dark phases; and 3) smaller capacity of non-photic synchronization for auditory cues from conspecifics, possibly due to reducing sensory inputs and responsiveness of the circadian oscillators to auditory cues, what can make the aged marmoset most vulnerable, as these social cues may act as an important supporting factor for the photic synchronization.

Caracterização do ritmo circadiano de atividade e repouso em saguis idosos (Callithrix Jacchus)

Advanced age may become a limiting factor for the maintenance of rhythms in organisms, reducing the capacity of generation and synchronization of biological rhythms. In this study, the influence of aging on the expression of endogenous periodicity and synchronization (photic and social) of the circadian activity rhythm (CAR) was evaluated in a diurnal primate, the marmoset (Callithrix jacchus). This study had two approaches: one with longitudinal design, performed with a male marmoset in two different phases: adult (three years) and older (9 y.o.) (study 1) and the second, a transversal approach, with 6 old (♂: 9.7 ± 2.0 y.o.) and 11 adults animals (♂: 4.2 ± 0.8 y.o.) (study 2). The evaluation of the photic synchronization involved two conditions in LD (natural and artificial illuminations). In study 1, the animal was subjected to the following stages: LD (12:12 ~ 350: ~ 2 lx), LL (~ 350 lx) and LD resynchronization. In the second study, the animals were initially evaluated in natural LD, and then the same sequence stages of study 1. During the LL stage in study 2, the vocalizations of conspecifics kept in natural LD on the outside of the colony were considered temporal cue to the social synchronization. The record of the activity was performed automatically at intervals of five minutes through infrared sensor and actimeters, in studies 1 and 2, respectively. In general, the aged showed a more fragmented activity pattern (> IV < H and > PSD, ANOVA, p < 0.05), lower levels of activity (ANOVA, p < 0.05) and shorter duration of active phase (ANOVA, p < 0.05) in LD conditions, when compared to adults. In natural LD, the aged presented phase delay pronounced for onset and offset of active phase (ANOVA, p < 0.05), while the adults had the active phase more adjusted to light phase. Under artificial LD, there was phase advance and greater adjustment of onset and offset of activity in relation to the LD in the aged (ANOVA, p < 0.05). In LL, there was a positive correlation between age and the endogenous period () in the first 20 days (Spearman correlation, p < 0.05), with prolonged  held in two aged animals. In this condition, most adults showed free-running period of the circadian activity rhythm with  < 24 h for the first 30 days and later on relative coordination mediated by auditory cues. In study 2, the cross-correlation analysis between the activity profiles of the animals in LL with control animals kept under natural LD, found that there was less social synchronization in the aged. With the resubmission to the LD, the resynchronization rate was slower in the aged (t-test; p < 0.05) and in just one aged animal there was a loss of resynchronization capability. According to the data set, it is suggested that the aging in marmosets may be related to: 1) lower amplitude and greater fragmentation of the activity, accompanied to phase delay with extension of period, caused by changes in a photic input, in the generation and behavioral expression of the CAR; 2) lower capacity of the circadian activity rhythm to photic synchronization, that can become more robust in artificial lighting conditions, possibly due to the higher light intensities at the beginning of the active phase due to the abrupt transitions between the light and dark phases; and 3) smaller capacity of non-photic synchronization for auditory cues from conspecifics, possibly due to reducing sensory inputs and responsiveness of the circadian oscillators to auditory cues, what can make the aged marmoset most vulnerable, as these social cues may act as an important supporting factor for the photic synchronization.

Age determination and sedimentology on profile Lz1020 from Lake Hoare, Antarctica

Up to 2.3 m long sediment sequences were recovered from the deepest part of Lake Hoare in Taylor Valley, southern Victoria Land, Antarctica. Sedimentological, biogeochemical, and mineralogical analyses revealed a high spatial variability of these parameters in Lake Hoare. Five distinct lithological units were recognized. Radiocarbon dating of bulk organic carbon samples from the sediment sequences yielded apparently too old ages and significant age reversals, which prevented the establishment of reliable age-depth models. However, cross correlation of the sedimentary characteristics with those of sediment records from neighbouring Lake Fryxell indicates that the lowermost two units of the Lake Hoare sediment sequences were probably deposited during the final phase of proglacial Lake Washburn, which occupied Taylor Valley during the late Pleistocene and early Holocene. High amounts of angular gravel and the absence of fine-grained material imply a complete desiccation with subaerial conditions in the Lake Hoare basin in the middle of the Holocene. The late Holocene (< c. 3300 calendar yr BP) is characterized by the establishment of environmental conditions similar to those existing today. A late Holocene desiccation event, such as proposed in former studies, is not indicated in the sediment sequences recovered.

Physical properties and sedimentology of 4 profiles fom the South Atlantic

Ultrasonic P wavc transmission seismograms recorded on sediment cores have been analyzed to study the acoustic and estimate the clastic properties of marine sediments from different provinces dominated by terrigenous, calcareous, amI diatomaceous sedimentation. Instantaneous frequencies computed from the transmission seismograms are displayed as gray-shaded images to give an acoustic overview of the lithology of each core. Ccntirneter-scale variations in the ultrasonic waveforms associated with lithological changes are illustrated by wiggle traces in detail. Cross-correlation, multiple-filter, and spectral ratio techniques are applied to derive P wave velocities and attenuation coefficients. S wave velocities and attenuation coefficients, elastic moduli, and permeabilities are calculated by an inversion scheme based on the Biot-Stoll viscoelastic model. Together wilh porosity measurements, P and S wave scatter diagrams are constructed to characterize different sediment types by their velocity- and attenuation-porosity relationships. They demonstrate that terrigenous, calcareous, and diatomaceous sediments cover different velocity- and attenuation-porosity ranges. In terrigcnous sediments, P wave vclocities and attenuation coefficients decrease rapidly with increasing porosity, whereas S wave velocities and shear moduli are very low. Calcareous sediments behave similarly at relatively higher porosities. Foraminifera skeletons in compositions of terrigenous mud and calcareous ooze cause a stiffening of the frame accompanied by higher shear moduli, P wave velocities, and attenuation coefficients. In diatomaceous ooze the contribution of the shear modulus becomes increasingly important and is controlled by the opal content, whereas attenuation is very low. This leads to the opportunity to predict the opal content from nondestructive P wave velocity measurements at centimeter-scale resolution.

RELIBILITY ANALYSIS OF SIMPLE SLOPES AND SOIL-STRUCTURES WITH LINEAR LIMIT STATES

The first objective of this research was to develop closed-form and numerical probabilistic methods of analysis that can be applied to otherwise conventional methods of unreinforced and geosynthetic reinforced slopes and walls. These probabilistic methods explicitly include random variability of soil and reinforcement, spatial variability of the soil, and cross-correlation between soil input parameters on probability of failure. The quantitative impact of simultaneously considering the influence of random and/or spatial variability in soil properties in combination with cross-correlation in soil properties is investigated for the first time in the research literature. Depending on the magnitude of these statistical descriptors, margins of safety based on conventional notions of safety may be very different from margins of safety expressed in terms of probability of failure (or reliability index). The thesis work also shows that intuitive notions of margin of safety using conventional factor of safety and probability of failure can be brought into alignment when cross-correlation between soil properties is considered in a rigorous manner. The second objective of this thesis work was to develop a general closed-form solution to compute the true probability of failure (or reliability index) of a simple linear limit state function with one load term and one resistance term expressed first in general probabilistic terms and then migrated to a LRFD format for the purpose of LRFD calibration. The formulation considers contributions to probability of failure due to model type, uncertainty in bias values, bias dependencies, uncertainty in estimates of nominal values for correlated and uncorrelated load and resistance terms, and average margin of safety expressed as the operational factor of safety (OFS). Bias is defined as the ratio of measured to predicted value. Parametric analyses were carried out to show that ignoring possible correlations between random variables can lead to conservative (safe) values of resistance factor in some cases and in other cases to non-conservative (unsafe) values. Example LRFD calibrations were carried out using different load and resistance models for the pullout internal stability limit state of steel strip and geosynthetic reinforced soil walls together with matching bias data reported in the literature.

Italian Ragweed Pollen Inventory

This study provides the first spatially detailed and complete inventory of Ambrosia pollen sources in Italy – the third largest centre of ragweed in Europe. The inventory relies on a well tested top-down approach that combines local knowledge, detailed land cover, pollen observations and a digital elevation model that assumes permanent ragweed populations mainly grow below 745m. The pollen data were obtained from 92 volumetric pollen traps located throughout Italy during 2004-2013. Land cover is derived from Corine Land cover information with 100m resolution. The digital elevation model is based on the NASA shuttle radar mission with 90m resolution. The inventory is produced using a combination of ArcGIS and Python for automation and validated using cross-correlation and has a final resolution of 5km x 5km. The method includes a harmonization of the inventory with other European inventories for the Pannonian Plain, France and Austria in order to provide a coherent picture of all major ragweed sources. The results show that the mean annual pollen index varies from 0 in South Italy to 6779 in the Po Valley. The results also show that very large pollen indexes are observed in the Milan region, but this region has smaller amounts of ragweed habitats compared to other parts of the Po Valley and known ragweed areas in France and the Pannonian Plain. A significant decrease in Ambrosia pollen concentrations was recorded in 2013 by pollen monitoring stations located in the Po Valley, particularly in the Northwest of Milan. This was the same year as the appearance of the Ophraella communa leaf beetle in Northern Italy. These results suggest that ragweed habitats near to the Milan region have very high densities of Ambrosia plants compared to other known ragweed habitats in Europe. The Milan region therefore appears to contain habitats with the largest ragweed infestation in Europe, but a smaller amount of habitats is a likely cause the pollen index to be lower compared to central parts of the Pannonian Plain. A low number of densely packed habitats may have increased the impact of the Ophraella beetle and might account for the documented decrease in airborne Ambrosia pollen levels, an event that cannot be explained by meteorology alone. Further investigations that model atmospheric pollen before and after the appearance of the beetle in this part of Northern Italy are needed to assess the influence of the beetle on airborne Ambrosia pollen concentrations. Future work will focus on short distance transport episodes for stations located in the Po Valley, and long distance transport events for stations in Central Italy that exhibit peaks in daily airborne Ambrosia pollen levels.

Hydrometallurgical process modeling using advanced mathematical tools

Hydrometallurgical process modeling is the main objective of this Master’s thesis work. Three different leaching processes namely, high pressure pyrite oxidation, direct oxidation zinc concentrate (sphalerite) leaching and gold chloride leaching using rotating disc electrode (RDE) are modeled and simulated using gPROMS process simulation program in order to evaluate its model building capabilities. The leaching mechanism in each case is described in terms of a shrinking core model. The mathematical modeling carried out included process model development based on available literature, estimation of reaction kinetic parameters and assessment of the model reliability by checking the goodness fit and checking the cross correlation between the estimated parameters through the use of correlation matrices. The estimated parameter values in each case were compared with those obtained using the Modest simulation program. Further, based on the estimated reaction kinetic parameters, reactor simulation and modeling for direct oxidation zinc concentrate (sphalerite) leaching is carried out in Aspen Plus V8.6. The zinc leaching autoclave is based on Cominco reactor configuration and is modeled as a series of continuous stirred reactors (CSTRs). The sphalerite conversion is calculated and a sensitivity analysis is carried out so to determine the optimum reactor operation temperature and optimum oxygen mass flow rate. In this way, the implementation of reaction kinetic models into the process flowsheet simulation environment has been demonstrated.

Seismic waveform inversion for core-mantle boundary topography

International audience

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience