Individual-based modelling and simulation: towards a better understanding of growth dynamics from small inocula

Postprint (published version)

Prevalence and predictors of suicide attempt in an incidence cohort of 661 young people with first-episode psychosis

OBJECTIVES: Studies investigating suicidal behaviour in psychosis rarely focus on incidence cohorts of first-episode patients. This is important, because patients who refuse study participation have higher rates of comorbid substance use disorders and longer duration of untreated psychosis as well as worse course illness, variables potentially linked to higher prevalence of suicidal behaviour. The aims of the present study were therefore to examine the prevalence and predictors of suicide and suicide attempt before and during the first 18-24 months of treatment. METHOD: A retrospective file audit of 661 patients was carried out. RESULTS: Six patients (0.9%) died by suicide, 93 (14.3%) attempted suicide prior to entry, and 57 (8.7%) did so during treatment. Predictors of suicide attempt were: previous attempt (odds ratio (OR)=45.54, 95% confidence interval (CI)=9.46-219.15), sexual abuse (OR=8.46, 95%CI=1.88-38.03), comorbid polysubstance (OR=13.63, 95%CI=2.58-71.99), greater insight (OR=0.17, 95%CI=0.06-0.49), lower baseline Global Assessment of Functioning Scale and Scale of Occupational and Functional Assessment score (OR=0.96, 95%CI=0.62-0.91; OR=0.98, 95%CI=0.95-0.99), and longer time in treatment (OR=1.05, 95%CI=1.03-1.08). CONCLUSIONS: The prevalence of suicidal behaviour was high, indicating that suicidal behaviour in incidence populations is higher than in non-epidemiological cohorts of first-episode patients. The rate of repetition of suicide attempt among the sample, however, was lower than expected, suggesting that specialist services can play a role in reducing suicide risk.

Energy-based modelling and simulation of the interconnection of a back-to-back converter and a doubly-fed induction machine

This paper describes the port interconnection of two subsystems: a power electronics subsystem (a back-to-back AC/CA converter (B2B), coupled to a phase of the power grid), and an electromechanical subsystem (a doubly-fed induction machine (DFIM). The B2B is a variable structure system (VSS), due to presence of control-actuated switches: however, from a modelling simulation, as well as a control-design, point of view, it is sensible to consider modulated transformers (MTF in the bond graph language) instead of the pairs of complementary switches. The port-Hamiltonian models of both subsystems are presented and, using a power-preserving interconnection, the Hamiltonian description of the whole system is obtained; detailed bond graphs of all subsystems and the complete system are also provided. Using passivity-based controllers computed in the Hamiltonian formalism for both subsystems, the whole model is simulated; simulations are run to rest the correctness and efficiency of the Hamiltonian network modelling approach used in this work.

General electromagnetic simulation tool to predict the microwave nonlinear response of planar, arbitrarily-shaped HTS structures

This work describes a simulation tool being developed at UPC to predict the microwave nonlinear behavior of planar superconducting structures with very few restrictions on the geometry of the planar layout. The software is intended to be applicable to most structures used in planar HTS circuits, including line, patch, and quasi-lumped microstrip resonators. The tool combines Method of Moments (MoM) algorithms for general electromagnetic simulation with Harmonic Balance algorithms to take into account the nonlinearities in the HTS material. The Method of Moments code is based on discretization of the Electric Field Integral Equation in Rao, Wilton and Glisson Basis Functions. The multilayer dyadic Green's function is used with Sommerfeld integral formulation. The Harmonic Balance algorithm has been adapted to this application where the nonlinearity is distributed and where compatibility with the MoM algorithm is required. Tests of the algorithm in TM010 disk resonators agree with closed-form equations for both the fundamental and third-order intermodulation currents. Simulations of hairpin resonators show good qualitative agreement with previously published results, but it is found that a finer meshing would be necessary to get correct quantitative results. Possible improvements are suggested.

Effectiveness of anti-psychotics and related drugs in the Huntington French-speaking group cohort.

PURPOSE: Huntington's disease is a rare condition. Patients are commonly treated with antipsychotics and tetrabenazine. The evidence of their effect on disease progression is limited and no comparative study between these drugs has been conducted. We therefore compared the effectiveness of antipsychotics on disease progression. METHODS: 956 patients from the Huntington French Speaking Group were followed for up to 8 years between 2002 and 2010. The effectiveness of treatments was assessed using Unified Huntington's Disease Rating Scale (UHDRS) scores and then compared using a mixed model adjusted on a multiple propensity score. RESULTS: 63% of patients were treated with antipsychotics during the survey period. The most commonly prescribed medications were dibenzodiazepines (38%), risperidone (13%), tetrabenazine (12%) and benzamides (12%). There was no difference between treatments on the motor and behavioural declines observed, after taking the patient profiles at the start of the drug prescription into account. In contrast, the functional decline was lower in the dibenzodiazepine group than the other antipsychotic groups (Total Functional Capacity: 0.41 ± 0.17 units per year vs. risperidone and 0.54 ± 0.19 vs. tetrabenazine, both p<0.05). Benzamides were less effective than other antipsychotics on cognitive evolution (Stroop interference, Stroop color and Literal fluency: p<0.05). CONCLUSIONS: Antipsychotics are widely used to treat patients with Huntington's disease. Although differences in motor or behavioural profiles between patients according to the antipsychotics used were small, there were differences in drug effectiveness on the evolution of functional and cognitive scores.

Modeling and Simulation ofReal-Time System as a Tool to Improve Capacity Management - Case Fenix

Tämän tutkimustyön kohteena on TietoEnator Oy:n kehittämän Fenix-tietojärjestelmän kapasiteettitarpeen ennustaminen. Työn tavoitteena on tutustua Fenix-järjestelmän eri osa-alueisiin, löytää tapa eritellä ja mallintaa eri osa-alueiden vaikutus järjestelmän kuormitukseen ja selvittää alustavasti mitkä parametrit vaikuttavat kyseisten osa-alueiden luomaan kuormitukseen. Osa tätä työtä on tutkia eri vaihtoehtoja simuloinnille ja selvittää eri vaihtoehtojen soveltuvuus monimutkaisten järjestelmien mallintamiseen. Kerätyn tiedon pohjaltaluodaan järjestelmäntietovaraston kuormitusta kuvaava simulaatiomalli. Hyödyntämällä mallista saatua tietoa ja tuotantojärjestelmästä mitattua tietoa mallia kehitetään vastaamaan yhä lähemmin todellisen järjestelmän toimintaa. Mallista tarkastellaan esimerkiksi simuloitua järjestelmäkuormaa ja jonojen käyttäytymistä. Tuotantojärjestelmästä mitataan eri kuormalähteiden käytösmuutoksia esimerkiksi käyttäjämäärän ja kellonajan suhteessa. Tämän työn tulosten on tarkoitus toimia pohjana myöhemmin tehtävälle jatkotutkimukselle, jossa osa-alueiden parametrisointia tarkennetaan lisää, mallin kykyä kuvata todellista järjestelmää tehostetaanja mallin laajuutta kasvatetaan.

Genetic polymorphism of CCR5 gene and HIV disease: the heterozygous (CCR5/delta ccr5) genotype is neither essential nor sufficient for protection against disease progression. Swiss HIV Cohort.

Homozygous (delta ccr5/delta ccr5) and heterozygous (CCR5/delta ccr5) deletions in the beta-chemokine receptor 5 (CCR5) gene, which encodes for the major co-receptor for macrophage-tropic HIV-1 entry, have been implicated in resistance to HIV infection and in protection against disease progression, respectively. The CCR5/delta ccr5 genotype was found more frequently in long-term nonprogressors (LTNP) (31.0%) than in progressors (10.6%, p < 0.0001), in agreement with previous studies. Kaplan-Meier survival analyses showed that a slower progression of disease, i.e. higher proportion of subjects with CD4+ T cell counts > 500/microl (p = 0.0006) and a trend toward a slower progression to AIDS (p = 0.077), was associated with the CCR5/delta ccr5 genotype. However, when LTNP were analyzed separately, no significant differences in CD4+ T cell counts (p = 0.12) and viremia levels (p = 0.65) were observed between the wild-type (69% of LTNP) and the heterozygous (31.0%) genotypes. Therefore, there are other factors which play a major role in determining the status of nonprogression in the majority of LTNP. Furthermore, there was no evidence that the CCR5/delta ccr5 genotype was associated with different rates of disease progression in the group of progressors. Taken together, these results indicate that the CCR5/delta ccr5 genotype is neither essential nor sufficient for protection against the progression of HIV disease.

Childhood cancer and nuclear power plants in Switzerland: a census-based cohort study.

BACKGROUND: Previous studies on childhood cancer and nuclear power plants (NPPs) produced conflicting results. We used a cohort approach to examine whether residence near NPPs was associated with leukaemia or any childhood cancer in Switzerland. METHODS: We computed person-years at risk for children aged 0-15 years born in Switzerland from 1985 to 2009, based on the Swiss censuses 1990 and 2000 and identified cancer cases from the Swiss Childhood Cancer Registry. We geo-coded place of residence at birth and calculated incidence rate ratios (IRRs) with 95% confidence intervals (CIs) comparing the risk of cancer in children born <5 km, 5-10 km and 10-15 km from the nearest NPP with children born >15 km away, using Poisson regression models. RESULTS: We included 2925 children diagnosed with cancer during 21 117 524 person-years of follow-up; 953 (32.6%) had leukaemia. Eight and 12 children diagnosed with leukaemia at ages 0-4 and 0-15 years, and 18 and 31 children diagnosed with any cancer were born <5 km from a NPP. Compared with children born >15 km away, the IRRs (95% CI) for leukaemia in 0-4 and 0-15 year olds were 1.20 (0.60-2.41) and 1.05 (0.60-1.86), respectively. For any cancer, corresponding IRRs were 0.97 (0.61-1.54) and 0.89 (0.63-1.27). There was no evidence of a dose-response relationship with distance (P > 0.30). Results were similar for residence at diagnosis and at birth, and when adjusted for potential confounders. Results from sensitivity analyses were consistent with main results. CONCLUSIONS: This nationwide cohort study found little evidence of an association between residence near NPPs and the risk of leukaemia or any childhood cancer.

The Dynamic Analysis of a Human Skeleton Using the Flexile Multibody Simulation Approach

The human motion study, which relies on mathematical and computational models ingeneral, and multibody dynamic biomechanical models in particular, has become asubject of many recent researches. The human body model can be applied to different physical exercises and many important results such as muscle forces, which are difficult to be measured through practical experiments, can be obtained easily. In the work, human skeletal lower limb model consisting of three bodies in build using the flexible multibody dynamics simulation approach. The floating frame of reference formulation is used to account for the flexibility in the bones of the human lower limb model. The main reason of considering the flexibility inthe human bones is to measure the strains in the bone result from different physical exercises. It has been perceived the bone under strain will become stronger in order to cope with the exercise. On the other hand, the bone strength is considered and important factors in reducing the bone fractures. The simulation approach and model developed in this work are used to measure the bone strain results from applying raising the sole of the foot exercise. The simulation results are compared to the results available in literature. The comparison shows goof agreement. This study sheds the light on the importance of using the flexible multibody dynamic simulation approach to build human biomechanical models, which can be used in developing some exercises to achieve the optimalbone strength.

Development of a simulation tool for torsional vibration of branched power transmission system

In this master's thesis a mechanical model that is driven with variable speed synchronous machine was developed. The developed mechanical model simulates the mechanics of power transmission and its torsional vibrations. The mechanical model was developed for the need of the branched mechanics of a rolling mill and the propulsion system of a tanker. First, the scope of the thesis was to clarify the concepts connected to the mechanical model. The clarified concepts are the variable speed drive, the mechanics of power transmission and the vibrationsin the power transmission. Next, the mechanical model with straight shaft line and twelve moments of inertia that existed in the beginning was developed to be branched considering the case of parallel machines and the case of parallel rolls. Additionally, the model was expanded for the need of moreaccurate simulation to up to thirty moments of inertia. The model was also enhanced to enable three phase short circuit situation of the simulated machine. After that the mechanical model was validated by comparing the results of the developed simulation tool to results of other simulation tools. The compared results are the natural frequencies and mode shapes of torsional vibration, the response of the load torque step and the stress in the mechanical system occurred by the permutation of the magnetic field that is arisen from the three phase short circuit situation. The comparisons were accomplished well and the mechanical model was validated for the compared cases. Further development to be made is to develop the load torque to be time-dependent and to install two frequency converters and two FEM modeled machines to be simulated parallel.

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Simulation of stromatolite growth using diffusion-limited aggregation and a discrete model of biogeochemical exchanges in microbial mat

Resume : Mieux comprendre les stromatolithes et les tapis microbiens est un sujet important en biogéosciences puisque cela aide à l'étude des premières formes de vie sur Terre, a mieux cerner l'écologie des communautés microbiennes et la contribution des microorganismes a la biominéralisation, et même à poser certains fondements dans les recherches en exobiologie. D'autre part, la modélisation est un outil puissant utilisé dans les sciences naturelles pour appréhender différents phénomènes de façon théorique. Les modèles sont généralement construits sur un système d'équations différentielles et les résultats sont obtenus en résolvant ce système. Les logiciels disponibles pour implémenter les modèles incluent les logiciels mathématiques et les logiciels généraux de simulation. L'objectif principal de cette thèse est de développer des modèles et des logiciels pour aider a comprendre, via la simulation, le fonctionnement des stromatolithes et des tapis microbiens. Ces logiciels ont été développés en C++ en ne partant d'aucun pré-requis de façon a privilégier performance et flexibilité maximales. Cette démarche permet de construire des modèles bien plus spécifiques et plus appropriés aux phénomènes a modéliser. Premièrement, nous avons étudié la croissance et la morphologie des stromatolithes. Nous avons construit un modèle tridimensionnel fondé sur l'agrégation par diffusion limitée. Le modèle a été implémenté en deux applications C++: un moteur de simulation capable d'exécuter un batch de simulations et de produire des fichiers de résultats, et un outil de visualisation qui permet d'analyser les résultats en trois dimensions. Après avoir vérifié que ce modèle peut en effet reproduire la croissance et la morphologie de plusieurs types de stromatolithes, nous avons introduit un processus de sédimentation comme facteur externe. Ceci nous a mené a des résultats intéressants, et permis de soutenir l'hypothèse que la morphologie des stromatolithes pourrait être le résultat de facteurs externes autant que de facteurs internes. Ceci est important car la classification des stromatolithes est généralement fondée sur leur morphologie, imposant que la forme d'un stromatolithe est dépendante de facteurs internes uniquement (c'est-à-dire les tapis microbiens). Les résultats avancés dans ce mémoire contredisent donc ces assertions communément admises. Ensuite, nous avons décidé de mener des recherches plus en profondeur sur les aspects fonctionnels des tapis microbiens. Nous avons construit un modèle bidimensionnel de réaction-diffusion fondé sur la simulation discrète. Ce modèle a été implémenté dans une application C++ qui permet de paramétrer et exécuter des simulations. Nous avons ensuite pu comparer les résultats de simulation avec des données du monde réel et vérifier que le modèle peut en effet imiter le comportement de certains tapis microbiens. Ainsi, nous avons pu émettre et vérifier des hypothèses sur le fonctionnement de certains tapis microbiens pour nous aider à mieux en comprendre certains aspects, comme la dynamique des éléments, en particulier le soufre et l'oxygène. En conclusion, ce travail a abouti à l'écriture de logiciels dédiés à la simulation de tapis microbiens d'un point de vue tant morphologique que fonctionnel, suivant deux approches différentes, l'une holistique, l'autre plus analytique. Ces logiciels sont gratuits et diffusés sous licence GPL (General Public License). Abstract : Better understanding of stromatolites and microbial mats is an important topic in biogeosciences as it helps studying the early forms of life on Earth, provides clues re- garding the ecology of microbial ecosystems and their contribution to biomineralization, and gives basis to a new science, exobiology. On the other hand, modelling is a powerful tool used in natural sciences for the theoretical approach of various phenomena. Models are usually built on a system of differential equations and results are obtained by solving that system. Available software to implement models includes mathematical solvers and general simulation software. The main objective of this thesis is to develop models and software able to help to understand the functioning of stromatolites and microbial mats. Software was developed in C++ from scratch for maximum performance and flexibility. This allows to build models much more specific to a phenomenon rather than general software. First, we studied stromatolite growth and morphology. We built a three-dimensional model based on diffusion-limited aggregation. The model was implemented in two C++ applications: a simulator engine, which can run a batch of simulations and produce result files, and a Visualization tool, which allows results to be analysed in three dimensions. After verifying that our model can indeed reproduce the growth and morphology of several types of stromatolites, we introduced a sedimentation process as an external factor. This lead to interesting results, and allowed to emit the hypothesis that stromatolite morphology may be the result of external factors as much as internal factors. This is important as stromatolite classification is usually based on their morphology, imposing that a stromatolite shape is dependant on internal factors only (i.e. the microbial mat). This statement is contradicted by our findings, Second, we decided to investigate deeper the functioning of microbial mats, We built a two-dimensional reaction-diffusion model based on discrete simulation, The model was implemented in a C++ application that allows setting and running simulations. We could then compare simulation results with real world data and verify that our model can indeed mimic the behaviour of some microbial mats. Thus, we have proposed and verified hypotheses regarding microbial mats functioning in order to help to better understand them, e.g. the cycle of some elements such as oxygen or sulfur. ln conclusion, this PhD provides a simulation software, dealing with two different approaches. This software is free and available under a GPL licence.

Treatment patterns and health care resource utilization in a 1-year observational cohort study of outpatients with schizophrenia at risk of nonadherence treated with long-acting injectable antipsychotics

Purpose: To describe (1) the clinical profiles and the patterns of use of long-acting injectable (LAI) antipsychotics in patients with schizophrenia at risk of nonadherence with oral antipsychotics, and in those who started treatment with LAI antipsychotics, (2) health care resource utilization and associated costs. Patients and methods: A total of 597 outpatients with schizophrenia at risk of nonadherence, according to the psychiatrist's clinical judgment, were recruited at 59 centers in a noninterventional prospective observational study of 1-year follow-up when their treatment was modified. In a post hoc analysis, the profiles of patients starting LAI or continuing with oral antipsychotics were described, and descriptive analyses of treatments, health resource utilization, and direct costs were performed in those who started an LAI antipsychotic. Results: Therapy modifications involved the antipsychotic medications in 84.8% of patients, mostly because of insufficient efficacy of prior regimen. Ninety-two (15.4%) patients started an LAI antipsychotic at recruitment. Of these, only 13 (14.1%) were prescribed with first-generation antipsychotics. During 1 year, 16.3% of patients who started and 14.9% of patients who did not start an LAI antipsychotic at recruitment relapsed, contrasting with the 20.9% who had been hospitalized only within the prior 6 months. After 1 year, 74.3% of patients who started an LAI antipsychotic continued concomitant treatment with oral antipsychotics. The mean (median) total direct health care cost per patient per month during the study year among the patients starting any LAI antipsychotic at baseline was 1,407 ( 897.7). Medication costs (including oral and LAI antipsychotics and concomitant medication) represented almost 44%, whereas nonmedication costs accounted for more than 55% of the mean total direct health care costs. Conclusion: LAI antipsychotics were infrequently prescribed in spite of a psychiatrist-perceived risk of nonadherence to oral antipsychotics. Mean medication costs were lower than nonmedication costs.

Low current and nadir CD4+ T-cell counts are associated with higher hepatitis C virus RNA levels in the Swiss HIV cohort study.

BACKGROUND: The aim of this study was to evaluate the effect of CD4+ T-cell counts and other characteristics of HIV-infected individuals on hepatitis C virus (HCV) RNA levels. METHODS: All HIV-HCV-coinfected Swiss HIV Cohort Study participants with available HCV RNA levels and concurrent CD4+ T-cell counts before starting HCV therapy were included. Potential predictors of HCV RNA levels were assessed by multivariate censored linear regression models that adjust for censored values. RESULTS: The study included 1,031 individuals. Low current and nadir CD4+ T-cell counts were significantly associated with higher HCV RNA levels (P = 0.004 and 0.001, respectively). In individuals with current CD4+ T-cell counts < 200/microl, median HCV RNA levels (6.22 log10 IU/ml) were +0.14 and +0.24 log10 IU/ml higher than those with CD4+ T-cell counts of 200-500/microl and > 500/microl. Based on nadir CD4+ T-cell counts, median HCV RNA levels (6.12 log10 IU/ml) in individuals with < 200/microl CD4+ T-cells were +0.06 and +0.44 log10 IU/ml higher than those with nadir T-cell counts of 200-500/microl and > 500/microl. Median HCV RNA levels were also significantly associated with HCV genotype: lower values were associated with genotype 4 and higher values with genotype 2, as compared with genotype 1. Additional significant predictors of lower HCV RNA levels were female gender and HIV transmission through male homosexual contacts. In multivariate analyses, only CD4+ T-cell counts and HCV genotype remained significant predictors of HCV RNA levels. Conclusions: Higher HCV RNA levels were associated with CD4+ T-cell depletion. This finding is in line with the crucial role of CD4+ T-cells in the control of HCV infection.

Simulation of Boundary Layers and Heat Transfer in the Entrance Region of Pipes

The study of fluid flow in pipes is one of the main topic of interest for engineers in industries. In this thesis, an effort is made to study the boundary layers formed near the wall of the pipe and how it behaves as a resistance to heat transfer. Before few decades, the scientists used to derive the analytical and empirical results by hand as there were limited means available to solve the complex fluid flow phenomena. Due to the increase in technology, now it has been practically possible to understand and analyze the actual fluid flow in any type of geometry. Several methodologies have been used in the past to analyze the boundary layer equations and to derive the expression for heat transfer. An integral relation approach is used for the analytical solution of the boundary layer equations and is compared with the FLUENT simulations for the laminar case. Law of the wall approach is used to derive the empirical correlation between dimensionless numbers and is then compared with the results from FLUENT for the turbulent case. In this thesis, different approaches like analytical, empirical and numerical are compared for the same set of fluid flow equations.