Modeling and simulation of population balances for particulate processes

This work focuses on the modeling and numerical approximations of population balance equations (PBEs) for the simulation of different phenomena occurring in process engineering. The population balance equation (PBE) is considered to be a statement of continuity. It tracks the change in particle size distribution as particles are born, die, grow or leave a given control volume. In the population balance models the one independent variable represents the time, the other(s) are property coordinate(s), e.g., the particle volume (size) in the present case. They typically describe the temporal evolution of the number density functions and have been used to model various processes such as granulation, crystallization, polymerization, emulsion and cell dynamics. The semi-discrete high resolution schemes are proposed for solving PBEs modeling one and two-dimensional batch crystallization models. The schemes are discrete in property coordinates but continuous in time. The resulting ordinary differential equations can be solved by any standard ODE solver. To improve the numerical accuracy of the schemes a moving mesh technique is introduced in both one and two-dimensional cases ...

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Obedience-disobedience dynamic and the role of responsibility

Obedience has been thoroughly studied in social psychology, both in its positive and negative aspects. Nevertheless, in these empirical studies disobedience has been considered to be the opposite of obedience and indeed its negation. Instead, some recent studies suggest that if obedience to authority is important in ensuring the continuity of social and group life, disobedience is crucial, under some circumstances, in stopping the authority relationship from degenerating into an authoritarian relationship. In this perspective, disobedience may be conceived of as a protest undermining the legitimacy of authority, or else it can represent an instrument of the community for controlling the legitimacy of the authority's demands, becoming a factor safeguarding against authoritarianism. The aim of the present study was to empirically verify the dynamics existing between disobedience and obedience. The results show that people who attach importance to both obedience and disobedience in the relationship between the individual and society recognize the importance of democratic values and consider themselves responsible for the defence of human rights. Instead, people who only recognize the value of obedience and consider disobedience as a threat to the status quo are more authoritarian, individualistic people.

Characterization and simulation of a CdTe detector for use in PET

The Voxel Imaging PET (VIP) Path nder project got the 4 year European Research Council FP7 grant in 2010 to prove the feasibility of using CdTe detectors in a novel conceptual design of PET scanner. The work presented in this thesis is a part of the VIP project and consists of, on the one hand, the characterization of a CdTe detector in terms of energy resolution and coincidence time resolution and, on the other hand, the simulation of the setup with the single detector in order to extend the results to the full PET scanner. An energy resolution of 0.98% at 511 keV with a bias voltage of 1000 V/mm has been measured at low temperature T=-8 ºC. The coincidence time distribution of two twin detectors has been found to be as low as 6 ns FWHM for events with energies above 500 keV under the same temperature and bias conditions. The measured energy and time resolution values are compatible with similar ndings available in the literature and prove the excellent potential of CdTe for PET applications. This results have been presented in form of a poster contribution at the IEEE NSS/MIC & RTSD 2011 conference in October 2011 in Valencia and at the iWoRID 2012 conference in July 2012 in Coimbra, Portugal. They have been also submitted for publication to "Journal of Instrumentation (JINST)" in September 2012.

Prediction of infant drug exposure through breastfeeding: population PK modeling and simulation of fluoxetine

BACKGROUND: Risks of significant infant drug exposurethrough breastmilk are poorly defined for many drugs, and largescalepopulation data are lacking. We used population pharmacokinetics(PK) modeling to predict fluoxetine exposure levels ofinfants via mother's milk in a simulated population of 1000 motherinfantpairs.METHODS: Using our original data on fluoxetine PK of 25breastfeeding women, a population PK model was developed withNONMEM and parameters, including milk concentrations, wereestimated. An exponential distribution model was used to account forindividual variation. Simulation random and distribution-constrainedassignment of doses, dosing time, feeding intervals and milk volumewas conducted to generate 1000 mother-infant pairs with characteristicssuch as the steady-state serum concentrations (Css) and infantdose relative to the maternal weight-adjusted dose (relative infantdose: RID). Full bioavailability and a conservative point estimate of1-month-old infant CYP2D6 activity to be 20% of the adult value(adjusted by weigth) according to a recent study, were assumed forinfant Css calculations.RESULTS: A linear 2-compartment model was selected as thebest model. Derived parameters, including milk-to-plasma ratios(mean: 0.66; SD: 0.34; range, 0 - 1.1) were consistent with the valuesreported in the literature. The estimated RID was below 10% in >95%of infants. The model predicted median infant-mother Css ratio was0.096 (range 0.035 - 0.25); literature reported mean was 0.07 (range0-0.59). Moreover, the predicted incidence of infant-mother Css ratioof >0.2 was less than 1%.CONCLUSION: Our in silico model prediction is consistent withclinical observations, suggesting that substantial systemic fluoxetineexposure in infants through human milk is rare, but further analysisshould include active metabolites. Our approach may be valid forother drugs. [supported by CIHR and Swiss National Science Foundation(SNSF)]

Aqueous dynamic and histological findings after deep sclerectomy with collagen implant in an animal model.

AIM: The use of an animal model to study the aqueous dynamic and the histological findings after deep sclerectomy with (DSCI) and without collagen implant. METHODS: Deep sclerectomy was performed on rabbits' eyes. Eyes were randomly assigned to receive collagen implants. Measurements of intraocular pressure (IOP) and aqueous outflow facility using the constant pressure method through cannulation of the anterior chamber were performed. The system was filled with BSS and cationised ferritin. Histological assessment of the operative site was performed. Sections were stained with haematoxylin and eosin and with Prussian blue. Aqueous drainage vessels were identified by the reaction between ferritin and Prussian blue. All eyes were coded so that the investigator was blind to the type of surgery until the evaluation was completed. RESULTS: A significant decrease in IOP (p<0.05) was observed during the first 6 weeks after DSCI (mean IOP was 13.07 (2.95) mm Hg preoperatively and 9.08 (2.25) mm Hg at 6 weeks); DS without collagen implant revealed a significant decrease in IOP at weeks 4 and 8 after surgery (mean IOP 12.57 (3.52) mm Hg preoperatively, 9.45 (3.38) mm Hg at 4 weeks, and 9.22 (3.39) mm Hg at 8 weeks). Outflow facility was significantly increased throughout the 9 months of follow up in both DSCI and DS groups (p<0.05). The preoperative outflow facility (OF) was 0.15 (0.02) micro l/min/mm Hg. At 9 months, OF was 0.52 (0.28) microl/min/mm Hg and 0.46 (0.07) micro l/min/mm Hg for DSCI and DS respectively. Light microscopy studies showed the appearance of new aqueous drainage vessels in the sclera adjacent to the dissection site in DSCI and DS and the apparition of spindle cells lining the collagen implant in DSCI after 2 months. CONCLUSION: A significant IOP decrease was observed during the first weeks after DSCI and DS. DS with or without collagen implant provided a significant increase in outflow facility throughout the 9 months of follow up. This might be partly explained by new drainage vessels in the sclera surrounding the operated site. Microscopic studies revealed the appearance of spindle cells lining the collagen implant in DSCI after 2 months.

Predicting the potential for lightning activity in mediterranean storms based on the WRF model dynamic and microphysical fields

A new parameter is introduced: the lightning potential index (LPI), which is a measure of the potential for charge generation and separation that leads to lightning flashes in convective thunderstorms. The LPI is calculated within the charge separation region of clouds between 0 C and 20 C, where the noninductive mechanism involving collisions of ice and graupel particles in the presence of supercooled water is most effective. As shown in several case studies using the Weather Research and Forecasting (WRF) model with explicit microphysics, the LPI is highly correlated with observed lightning. It is suggested that the LPI may be a useful parameter for predicting lightning as well as a tool for improving weather forecasting of convective storms and heavy rainfall.

Dynamic and Regulated Association of Caveolin with Lipid Bodies: Modulation of Lipid Body Motility and Function by a Dominant Negative Mutant

Caveolins are a crucial component of caveolae but have also been localized to the Golgi complex, and, under some experimental conditions, to lipid bodies (LBs). The physiological relevance and dynamics of LB association remain unclear. We now show that endogenous caveolin-1 and caveolin-2 redistribute to LBs in lipid loaded A431 and FRT cells. Association with LBs is regulated and reversible; removal of fatty acids causes caveolin to rapidly leave the lipid body. We also show by subcellular fractionation, light and electron microscopy that during the first hours of liver regeneration, caveolins show a dramatic redistribution from the cell surface to the newly formed LBs. At later stages of the regeneration process (when LBs are still abundant), the levels of caveolins in LBs decrease dramatically. As a model system to study association of caveolins with LBs we have used brefeldin A (BFA). BFA causes rapid redistribution of endogenous caveolins to LBs and this association was reversed upon BFA washout. Finally, we have used a dominant negative LB-associated caveolin mutant (cavDGV) to study LB formation and to examine its effect on LB function. We now show that the cavDGV mutant inhibits microtubule-dependent LB motility and blocks the reversal of lipid accumulation in LBs.

Mobile phone use and brain tumors in children and adolescents: a multicenter case-control study.

Background It has been hypothesized that children and adolescents might be more vulnerable to possible health effects from mobile phone exposure than adults. We investigated whether mobile phone use is associated with brain tumor risk among children and adolescents. Methods CEFALO is a multicenter case-control study conducted in Denmark, Sweden, Norway, and Switzerland that includes all children and adolescents aged 7-19 years who were diagnosed with a brain tumor between 2004 and 2008. We conducted interviews, in person, with 352 case patients (participation rate: 83%) and 646 control subjects (participation rate: 71%) and their parents. Control subjects were randomly selected from population registries and matched by age, sex, and geographical region. We asked about mobile phone use and included mobile phone operator records when available. Odds ratios (ORs) for brain tumor risk and 95% confidence intervals (CIs) were calculated using conditional logistic regression models. Results Regular users of mobile phones were not statistically significantly more likely to have been diagnosed with brain tumors compared with nonusers (OR = 1.36; 95% CI = 0.92 to 2.02). Children who started to use mobile phones at least 5 years ago were not at increased risk compared with those who had never regularly used mobile phones (OR = 1.26, 95% CI = 0.70 to 2.28). In a subset of study participants for whom operator recorded data were available, brain tumor risk was related to the time elapsed since the mobile phone subscription was started but not to amount of use. No increased risk of brain tumors was observed for brain areas receiving the highest amount of exposure. Conclusion The absence of an exposure-response relationship either in terms of the amount of mobile phone use or by localization of the brain tumor argues against a causal association.

Measurement and simulation of anisotropy in the Infrared and Raman spectra of beta-FeSi2 single crystals

In this paper we show that the orthorhombic phase of FeSi2 (stable at room temperature) displays a sizable anisotropy in the infrared spectra, with minor effects in the Raman data too. This fact is not trivial at all, since the crystal structure corresponds to a moderate distortion of the fluorite symmetry. Our analysis is carried out on small single crystals grown by flux transport, through polarization-resolved far-infrared reflectivity and Raman measurements. Their interpretation has been obtained by means of the simulated spectra with tight-binding molecular dynamics.

Application of Visualization and Simulation Program to Improve Work Zone Safety and Mobility, January 2010

A previous study sponsored by the Smart Work Zone Deployment Initiative, “Feasibility of Visualization and Simulation Applications to Improve Work Zone Safety and Mobility,” demonstrated the feasibility of combining readily available, inexpensive software programs, such as SketchUp and Google Earth, with standard two-dimensional civil engineering design programs, such as MicroStation, to create animations of construction work zones. The animations reflect changes in work zone configurations as the project progresses, representing an opportunity to visually present complex information to drivers, construction workers, agency personnel, and the general public. The purpose of this study is to continue the work from the previous study to determine the added value and resource demands created by including more complex data, specifically traffic volume, movement, and vehicle type. This report describes the changes that were made to the simulation, including incorporating additional data and converting the simulation from a desktop application to a web application.

Visuo-spatial processing in a dynamic and a static working memory paradigm in schizophrenia.

Recent findings suggest that the visuo-spatial sketchpad (VSSP) may be divided into two sub-components processing dynamic or static visual information. This model may be useful to elucidate the confusion of data concerning the functioning of the VSSP in schizophrenia. The present study examined patients with schizophrenia and matched controls in a new working memory paradigm involving dynamic (the Ball Flight Task - BFT) or static (the Static Pattern Task - SPT) visual stimuli. In the BFT, the responses of the patients were apparently based on the retention of the last set of segments of the perceived trajectory, whereas control subjects relied on a more global strategy. We assume that the patients' performances are the result of a reduced capacity in chunking visual information since they relied mainly on the retention of the last set of segments. This assumption is confirmed by the poor performance of the patients in the static task (SPT), which requires a combination of stimulus components into object representations. We assume that the static/dynamic distinction may help us to understand the VSSP deficits in schizophrenia. This distinction also raises questions about the hypothesis that visuo-spatial working memory can simply be dissociated into visual and spatial sub-components.

Prediction of Infant Drug Exposure Through Breastfeeding: Population PK Modeling and Simulation of Fluoxetine Exposure.

The likelihood of significant exposure to drugs in infants through breast milk is poorly defined, given the difficulties of conducting pharmacokinetics (PK) studies. Using fluoxetine (FX) as an example, we conducted a proof-of-principle study applying population PK (popPK) modeling and simulation to estimate drug exposure in infants through breast milk. We simulated data for 1,000 mother-infant pairs, assuming conservatively that the FX clearance in an infant is 20% of the allometrically adjusted value in adults. The model-generated estimate of the milk-to-plasma ratio for FX (mean: 0.59) was consistent with those reported in other studies. The median infant-to-mother ratio of FX steady-state plasma concentrations predicted by the simulation was 8.5%. Although the disposition of the active metabolite, norfluoxetine, could not be modeled, popPK-informed simulation may be valid for other drugs, particularly those without active metabolites, thereby providing a practical alternative to conventional PK studies for exposure risk assessment in this population.

Dosage optimization of treatments using population pharmacokinetic modeling and simulation.

Pharmacokinetic variability in drug levels represent for some drugs a major determinant of treatment success, since sub-therapeutic concentrations might lead to toxic reactions, treatment discontinuation or inefficacy. This is true for most antiretroviral drugs, which exhibit high inter-patient variability in their pharmacokinetics that has been partially explained by some genetic and non-genetic factors. The population pharmacokinetic approach represents a very useful tool for the description of the dose-concentration relationship, the quantification of variability in the target population of patients and the identification of influencing factors. It can thus be used to make predictions and dosage adjustment optimization based on Bayesian therapeutic drug monitoring (TDM). This approach has been used to characterize the pharmacokinetics of nevirapine (NVP) in 137 HIV-positive patients followed within the frame of a TDM program. Among tested covariates, body weight, co-administration of a cytochrome (CYP) 3A4 inducer or boosted atazanavir as well as elevated aspartate transaminases showed an effect on NVP elimination. In addition, genetic polymorphism in the CYP2B6 was associated with reduced NVP clearance. Altogether, these factors could explain 26% in NVP variability. Model-based simulations were used to compare the adequacy of different dosage regimens in relation to the therapeutic target associated with treatment efficacy. In conclusion, the population approach is very useful to characterize the pharmacokinetic profile of drugs in a population of interest. The quantification and the identification of the sources of variability is a rational approach to making optimal dosage decision for certain drugs administered chronically.