Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Data Envelopment Analysis and non-discretionary inputs : How to select the most suitable model using multi-criteria decision analysis

Within Data Envelopment Analysis, several alternative models allow for an environmental adjustment. The majority of them deliver divergent results. Decision makers face the difficult task of selecting the most suitable model. This study is performed to overcome this difficulty. By doing so, it fills a research gap. First, a two-step web-based survey is conducted. It aims (1) to identify the selection criteria, (2) to prioritize and weight the selection criteria with respect to the goal of selecting the most suitable model and (3) to collect the preferences about which model is preferable to fulfil each selection criterion. Second, Analytic Hierarchy Process is used to quantify the preferences expressed in the survey. Results show that the understandability, the applicability and the acceptability of the alternative models are valid selection criteria. The selection of the most suitable model depends on the preferences of the decision makers with regards to these criteria.

Non-Linear Journal Bearing Model for Analysis of Superharmonic Vibrations of Rotor Systems

A rotating machine usually consists of a rotor and bearings that supports it. The nonidealities in these components may excite vibration of the rotating system. The uncontrolled vibrations may lead to excessive wearing of the components of the rotating machine or reduce the process quality. Vibrations may be harmful even when amplitudes are seemingly low, as is usually the case in superharmonic vibration that takes place below the first critical speed of the rotating machine. Superharmonic vibration is excited when the rotational velocity of the machine is a fraction of the natural frequency of the system. In such a situation, a part of the machine’s rotational energy is transformed into vibration energy. The amount of vibration energy should be minimised in the design of rotating machines. The superharmonic vibration phenomena can be studied by analysing the coupled rotor-bearing system employing a multibody simulation approach. This research is focused on the modelling of hydrodynamic journal bearings and rotorbearing systems supported by journal bearings. In particular, the non-idealities affecting the rotor-bearing system and their effect on the superharmonic vibration of the rotating system are analysed. A comparison of computationally efficient journal bearing models is carried out in order to validate one model for further development. The selected bearing model is improved in order to take the waviness of the shaft journal into account. The improved model is implemented and analyzed in a multibody simulation code. A rotor-bearing system that consists of a flexible tube roll, two journal bearings and a supporting structure is analysed employing the multibody simulation technique. The modelled non-idealities are the shell thickness variation in the tube roll and the waviness of the shaft journal in the bearing assembly. Both modelled non-idealities may cause subharmonic resonance in the system. In multibody simulation, the coupled effect of the non-idealities can be captured in the analysis. Additionally one non-ideality is presented that does not excite the vibrations itself but affects the response of the rotorbearing system, namely the waviness of the bearing bushing which is the non-rotating part of the bearing system. The modelled system is verified with measurements performed on a test rig. In the measurements the waviness of bearing bushing was not measured and therefore it’s affect on the response was not verified. In conclusion, the selected modelling approach is an appropriate method when analysing the response of the rotor-bearing system. When comparing the simulated results to the measured ones, the overall agreement between the results is concluded to be good.

Non-Markovian Dynamics of the Spin-Boson Model

Avointen kvanttisysteemien teoria tutkii, miten vuorovaikutus laajan ympäristön kanssa muuttaa kvanttisysteemin dynamiikkaa. Kaikki kvanttisysteemit ovat pohjimmiltaan avoimia, joten tällaisten systeemien ymmärtäminen on välttämätöntä kvanttimaailman ymmärtämiseksi. Spin-bosoni -malli kuvaa avointa kaksitasosysteemiä, joka vuorovaikuttaa bosoneista koostuvan ympäristön kanssa. Useat fysikaaliset systeemit käyttäytyvät kaksitasosysteemin tavoin matalissa lämpötiloissa, jonka vuoksi spin-bosoni -mallilla on merkittävä rooli avointen kvanttisysteemien tutkimuksessa. Tutkielman ensimmäinen luku käsittelee yleistä avointen kvanttisysteemien teoriaa ja menetelmiä ei-Markovisen liikeyhtälön johtamiseksi. Ei-Markovinen liikeyhtälö kuvaa ympäristöstä systeemiin takaisin virtaavan informaation vaikutusta systeemin dynamiikkaan. Toisessa luvussa on johdettu spin-bosoni -mallille ei-Markovinen liikeyhtälö ja tarkasteltu spin-systeemin dynamiikkaa joissain erityistapauksissa. Kolmannessa luvussa on esitelty kvanttihyppymenetelmiä avointen systeemien dynamiikan ratkaisemiseksi ja sovellettu ei-Markovisten kvanttihyppyjen teoriaa spinbosoni -mallin tutkimiseksi. Ei-Markovisten kvanttihyppyjen teoriaa käyttämällä voidaan selittää spin-systeemin tiheysmatriisin populaatioissa ja koherensseissa havaittavien ei-Markovisten efektien alkuperä käännettyjen kvanttihyppyjen avulla. Tutkielman viimeisessä osassa esitellään muutamia spin-bosoni -mallin sovellutuksia fysiikan eri alueilta.

BEM/FEM non-linear model applied to transient analysis with viscous damping

In this article a two-dimensional transient boundary element formulation based on the mass matrix approach is discussed. The implicit formulation of the method to deal with elastoplastic analysis is considered, as well as the way to deal with viscous damping effects. The time integration processes are based on the Newmark rhoand Houbolt methods, while the domain integrals for mass, elastoplastic and damping effects are carried out by the well known cell approximation technique. The boundary element algebraic relations are also coupled with finite element frame relations to solve stiffened domains. Some examples to illustrate the accuracy and efficiency of the proposed formulation are also presented.

Non-linear analysis of laminated metal matrix composites by an integrated micro/macro-mechanical model

The present paper describes an integrated micro/macro mechanical study of the elastic-viscoplastic behavior of unidirectional metal matrix composites (MMC). The micromechanical analysis of the elastic moduli is based on the Composites Cylinder Assemblage model (CCA) with comparisons also draw with a Representative Unit Cell (RUC) technique. These "homogenization" techniques are later incorporated into the Vanishing Fiber Diameter (VFD) model and a new formulation is proposed. The concept of a smeared element procedure is employed in conjunction with two different versions of the Bodner and Partom elastic-viscoplastic constitutive model for the associated macroscopic analysis. The formulations developed are also compared against experimental and analytical results available in the literature.

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Is pancreas development abnormal in the non-obese diabetic mouse, a spontaneous model of type I diabetes?

Despite extensive genetic and immunological research, the complex etiology and pathogenesis of type I diabetes remains unresolved. During the last few years, our attention has been focused on factors such as abnormalities of islet function and/or microenvironment, that could interact with immune partners in the spontaneous model of the disease, the non-obese diabetic (NOD) mouse. Intriguingly, the first anomalies that we noted in NOD mice, compared to control strains, are already present at birth and consist of 1) higher numbers of paradoxically hyperactive ß cells, assessed by in situ preproinsulin II expression; 2) high percentages of immature islets, representing islet neogenesis related to neonatal ß-cell hyperactivity and suggestive of in utero ß-cell stimulation; 3) elevated levels of some types of antigen-presenting cells and FasL+ cells, and 4) abnormalities of extracellular matrix (ECM) protein expression. However, the colocalization in all control mouse strains studied of fibroblast-like cells (anti-TR-7 labeling), some ECM proteins (particularly, fibronectin and collagen I), antigen-presenting cells and a few FasL+ cells at the periphery of islets undergoing neogenesis suggests that remodeling phenomena that normally take place during postnatal pancreas development could be disturbed in NOD mice. These data show that from birth onwards there is an intricate relationship between endocrine and immune events in the NOD mouse. They also suggest that tissue-specific autoimmune reactions could arise from developmental phenomena taking place during fetal life in which ECM-immune cell interaction(s) may play a key role.

Sleep deprivation reduces the lymphocyte count in a non-obese mouse model of type 1 diabetes mellitus

The objective of the present study was to determine whether sleep deprivation (SD) would promote changes in lymphocyte numbers in a type 1 diabetes model (non-obese diabetic, NOD, mouse strain) and to determine whether SD would affect female and male NOD compared to Swiss mice. The number of lymphocytes in peripheral blood after 24 and 96 h of SD (by multiple platform method) or equivalent period of time in home-cage controls was examined prior to the onset of diabetes. SD for 96 h significantly reduced lymphocytes in male Swiss mice compared to control (8.6 ± 2.1 vs 4.1 ± 0.7 10³/µL; P < 0.02). In male NOD animals, 24- and 96-h SD caused a significant decrease of lymphocytes compared to control (4.4 ± 0.3 vs 1.6 ± 0.5; P < 0.001 and 4.4 ± 0.3 vs 0.9 ± 0.1 10³/µL; P < 0.00001, respectively). Both 24- and 96-h SD induced a reduction in the number of lymphocytes in female Swiss (7.5 ± 0.5 vs 4.5 ± 0.5, 4.4 ± 0.6 10³/µL; P < 0.001, respectively) and NOD mice (4 ± 0.6 vs 1.8 ± 0.2, 1.2 ± 0.4 10³/µL; P < 0.01, respectively) compared to the respective controls. Loss of sleep induced lymphopenia in peripheral blood in both genders and strains used. Since many cases of autoimmunity present reduced numbers of lymphocytes and, in this study, it was more evident in the NOD strain, our results suggest that SD should be considered a risk factor in the onset of autoimmune disorders.

Sorafenib prevents liver fibrosis in a non-alcoholic steatohepatitis (NASH) rodent model

Liver fibrosis occurring as an outcome of non-alcoholic steatohepatitis (NASH) can precede the development of cirrhosis. We investigated the effects of sorafenib in preventing liver fibrosis in a rodent model of NASH. Adult Sprague-Dawley rats were fed a choline-deficient high-fat diet and exposed to diethylnitrosamine for 6 weeks. The NASH group (n=10) received vehicle and the sorafenib group (n=10) received 2.5 mg·kg-1·day-1 by gavage. A control group (n=4) received only standard diet and vehicle. Following treatment, animals were sacrificed and liver tissue was collected for histologic examination, mRNA isolation, and analysis of mitochondrial function. Genes related to fibrosis (MMP9, TIMP1, TIMP2), oxidative stress (HSP60, HSP90, GST), and mitochondrial biogenesis (PGC1α) were evaluated by real-time quantitative polymerase chain reaction (RT-qPCR). Liver mitochondrial oxidation activity was measured by a polarographic method, and cytokines by enzyme-linked immunosorbent assay (ELISA). Sorafenib treatment restored mitochondrial function and reduced collagen deposition by nearly 63% compared to the NASH group. Sorafenib upregulated PGC1α and MMP9 and reduced TIMP1 and TIMP2 mRNA and IL-6 and IL-10 protein expression. There were no differences in HSP60, HSP90 and GST expression. Sorafenib modulated PGC1α expression, improved mitochondrial respiration and prevented collagen deposition. It may, therefore, be useful in the treatment of liver fibrosis in NASH.

Evaluation of a new mathematical model to describe Clostridium perfringens growth during the cooling of cooked ground beef

A mathematical model previously developed to study microbial growth in food products under an isothermal environment was adapted to a time-varying temperature regime. The resulting model was applied to study the growth of Clostridium perfringens in meat products. This micro-organism is of particular relevance to public health and economy due to the loss of productivity caused by it. Results showed a similar performance of the model used compared to the Baranyi model under an isothermal situation and a slightly better performance under a non-isothermal temperature profile.

Non-offensive defence : a criteria model of military credibility

Finnish Defence Studies is published under the auspices of the War College, and the contributions reflect the fields of research and teaching of the College. Finnish Defence Studies will occasionally feature documentation on Finnish Security Policy. Views expressed are those of the authors and do not necessarily imply endorsement by the War College.

Asymptotic Approximations in the Near-Integrated Model with a Non-Zero Initial Condition

This paper considers various asymptotic approximations in the near-integrated firstorder autoregressive model with a non-zero initial condition. We first extend the work of Knight and Satchell (1993), who considered the random walk case with a zero initial condition, to derive the expansion of the relevant joint moment generating function in this more general framework. We also consider, as alternative approximations, the stochastic expansion of Phillips (1987c) and the continuous time approximation of Perron (1991). We assess how these alternative methods provide or not an adequate approximation to the finite-sample distribution of the least-squares estimator in a first-order autoregressive model. The results show that, when the initial condition is non-zero, Perron's (1991) continuous time approximation performs very well while the others only offer improvements when the initial condition is zero.

Abdominal aortic aneurysm follow-up after endovascular repair in a canine model with non-invasive vascular elastography

Le traitement chirurgical des anévrismes de l'aorte abdominale est de plus en plus remplacé par la réparation endovasculaire de l’anévrisme (« endovascular aneurysm repair », EVAR) en utilisant des endoprothèses (« stent-grafts », SGs). Cependant, l'efficacité de cette approche moins invasive est compromise par l'incidence de l'écoulement persistant dans l'anévrisme, appelé endofuites menant à une rupture d'anévrisme si elle n'est pas détectée. Par conséquent, une surveillance de longue durée par tomodensitométrie sur une base annuelle est nécessaire ce qui augmente le coût de la procédure EVAR, exposant le patient à un rayonnement ionisants et un agent de contraste néphrotoxique. Le mécanisme de rupture d'anévrisme secondaire à l'endofuite est lié à une pression du sac de l'anévrisme proche de la pression systémique. Il existe une relation entre la contraction ou l'expansion du sac et la pressurisation du sac. La pressurisation résiduelle de l'anévrisme aortique abdominale va induire une pulsation et une circulation sanguine à l'intérieur du sac empêchant ainsi la thrombose du sac et la guérison de l'anévrisme. L'élastographie vasculaire non-invasive (« non-invasive vascular elastography », NIVE) utilisant le « Lagrangian Speckle Model Estimator » (LSME) peut devenir une technique d'imagerie complémentaire pour le suivi des anévrismes après réparation endovasculaire. NIVE a la capacité de fournir des informations importantes sur l'organisation d'un thrombus dans le sac de l'anévrisme et sur la détection des endofuites. La caractérisation de l'organisation d'un thrombus n'a pas été possible dans une étude NIVE précédente. Une limitation de cette étude était l'absence d'examen tomodensitométrique comme étalon-or pour le diagnostic d'endofuites. Nous avons cherché à appliquer et optimiser la technique NIVE pour le suivi des anévrismes de l'aorte abdominale (AAA) après EVAR avec endoprothèse dans un modèle canin dans le but de détecter et caractériser les endofuites et l'organisation du thrombus. Des SGs ont été implantés dans un groupe de 18 chiens avec un anévrisme créé dans l'aorte abdominale. Des endofuites de type I ont été créés dans 4 anévrismes, de type II dans 13 anévrismes tandis qu’un anévrisme n’avait aucune endofuite. L'échographie Doppler (« Doppler ultrasound », DUS) et les examens NIVE ont été réalisés avant puis à 1 semaine, 1 mois, 3 mois et 6 mois après l’EVAR. Une angiographie, une tomodensitométrie et des coupes macroscopiques ont été réalisées au moment du sacrifice. Les valeurs de contrainte ont été calculées en utilisant l`algorithme LSME. Les régions d'endofuite, de thrombus frais (non organisé) et de thrombus solide (organisé) ont été identifiées et segmentées en comparant les résultats de la tomodensitométrie et de l’étude macroscopique. Les valeurs de contrainte dans les zones avec endofuite, thrombus frais et organisé ont été comparées. Les valeurs de contrainte étaient significativement différentes entre les zones d'endofuites, les zones de thrombus frais ou organisé et entre les zones de thrombus frais et organisé. Toutes les endofuites ont été clairement caractérisées par les examens d'élastographie. Aucune corrélation n'a été trouvée entre les valeurs de contrainte et le type d'endofuite, la pression de sac, la taille des endofuites et la taille de l'anévrisme.