Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

A new method for analysing the equilibrium and time-dependent behaviour of Markovian models

Many large-scale stochastic systems, such as telecommunications networks, can be modelled using a continuous-time Markov chain. However, it is frequently the case that a satisfactory analysis of their time-dependent, or even equilibrium, behaviour is impossible. In this paper, we propose a new method of analyzing Markovian models, whereby the existing transition structure is replaced by a more amenable one. Using rates of transition given by the equilibrium expected rates of the corresponding transitions of the original chain, we are able to approximate its behaviour. We present two formulations of the idea of expected rates. The first provides a method for analysing time-dependent behaviour, while the second provides a highly accurate means of analysing equilibrium behaviour. We shall illustrate our approach with reference to a variety of models, giving particular attention to queueing and loss networks. (C) 2003 Elsevier Ltd. All rights reserved.

A type system for value-dependent information flow analysis

Information systems are widespread and used by anyone with computing devices as well as corporations and governments. It is often the case that security leaks are introduced during the development of an application. Reasons for these security bugs are multiple but among them one can easily identify that it is very hard to define and enforce relevant security policies in modern software. This is because modern applications often rely on container sharing and multi-tenancy where, for instance, data can be stored in the same physical space but is logically mapped into different security compartments or data structures. In turn, these security compartments, to which data is classified into in security policies, can also be dynamic and depend on runtime data. In this thesis we introduce and develop the novel notion of dependent information flow types, and focus on the problem of ensuring data confidentiality in data-centric software. Dependent information flow types fit within the standard framework of dependent type theory, but, unlike usual dependent types, crucially allow the security level of a type, rather than just the structural data type itself, to depend on runtime values. Our dependent function and dependent sum information flow types provide a direct, natural and elegant way to express and enforce fine grained security policies on programs. Namely programs that manipulate structured data types in which the security level of a structure field may depend on values dynamically stored in other fields The main contribution of this work is an efficient analysis that allows programmers to verify, during the development phase, whether programs have information leaks, that is, it verifies whether programs protect the confidentiality of the information they manipulate. As such, we also implemented a prototype typechecker that can be found at http://ctp.di.fct.unl.pt/DIFTprototype/.

Nonparametric estimation of time-dependent ROC curves conditional on a continuous covariate

The receiver-operating characteristic (ROC) curve is the most widely used measure for evaluating the performance of a diagnostic biomarker when predicting a binary disease outcome. The ROC curve displays the true positive rate (or sensitivity) and the false positive rate (or 1-specificity) for different cut-off values used to classify an individual as healthy or diseased. In time-to-event studies, however, the disease status (e.g. death or alive) of an individual is not a fixed characteristic, and it varies along the study. In such cases, when evaluating the performance of the biomarker, several issues should be taken into account: first, the time-dependent nature of the disease status; and second, the presence of incomplete data (e.g. censored data typically present in survival studies). Accordingly, to assess the discrimination power of continuous biomarkers for time-dependent disease outcomes, time-dependent extensions of true positive rate, false positive rate, and ROC curve have been recently proposed. In this work, we present new nonparametric estimators of the cumulative/dynamic time-dependent ROC curve that allow accounting for the possible modifying effect of current or past covariate measures on the discriminatory power of the biomarker. The proposed estimators can accommodate right-censored data, as well as covariate-dependent censoring. The behavior of the estimators proposed in this study will be explored through simulations and illustrated using data from a cohort of patients who suffered from acute coronary syndrome.

Lagrangian methods for visualization and analysis of time-dependent vector fields

Magdeburg, Univ., Fak. für Informatik, Diss., 2013

A stable and 40Ar/39Ar isotope study of a major thrust in the Helvetic nappes (Swiss Alps) - Evidence for fluid-flow and constraints on nappe kinematics

Stable isotope and Ar-40/Ar-39 measurements,were made on samples associated with a major tectonic discontinuity in the Helvetic Alps, the basal thrust of the Diablerets nappe (external zone of the Alpine Belt) in order to determine both the importance of fluids in this thrust zone and the timing of thrusting. A systematic decrease in the delta(18)O values (up to 6 parts per thousand) of calcite, quartz, and white mica exists within a 10- to 70-m-wide zone over a distance of 37 km along the thrust, and they become more pronounced toward the root of the nappe. A similar decrease in the delta(13)C values of calcite is observed only in the deepest sections (up to 3 parts per thousand). The delta D-SMOW (SMOW = standard mean ocean water) values of white mica are -54 parts per thousand +/- 8 parts per thousand (n = 22) and are independent of the distance from the thrust. These variations are interpreted to reflect syntectonic solution reprecipitation during fluid passage along the thrust. The calculated delta(18)O and delta D values (versus SMOW) for the fluid in equilibrium with the analyzed minerals is 12 parts per thousand to 16 parts per thousand and -30 parts per thousand to +5 parts per thousand, respectively, for assumed temperatures of 250 to 450 degrees C. The isotopic and structural data are consistent with fluids derived from the deep-seated roots of the Helvetic nappes where large volumes of Mesozoic sediments were metamorphosed to the amphibolite facies, It is suggested that connate and metamorphic waters, overpressured by rapid tectonic burial in a subductive system escaped by upward infiltration along moderately dipping pathways until they reached the main shear zone at the base of the moving pile, where they were channeled toward the surface, This model also explains the mechanism by which large amounts of fluid were removed from the Mesozoic sediments during Alpine metamorphism. White mica Ar-49/Ar-39 ages vary from 27 Ma far from the Diablerets thrust to 15 Ma along the thrust. An older component is observed in micas far from the thrust, interpreted as a detrital signature, and indicates that regional metamorphic temperatures were less than about 350 degrees C. The;plateau and near plateau ages nearest the thrust are consistent with either neocrystallization of white mica or argon loss by recrystallization during thrusting, which may have been enhanced in the zones of highest fluid flow. The 15 Ma Ar-40/Ar-39 age plateau measured on white mica sampled exactly on the thrust surface dates the end of both fluid flow and tectonic transport.

Improving the accuracy of heat balance integral methods applied to thermal problems with time dependent boundary conditions

In this paper the two main drawbacks of the heat balance integral methods are examined. Firstly we investigate the choice of approximating function. For a standard polynomial form it is shown that combining the Heat Balance and Refined Integral methods to determine the power of the highest order term will either lead to the same, or more often, greatly improved accuracy on standard methods. Secondly we examine thermal problems with a time-dependent boundary condition. In doing so we develop a logarithmic approximating function. This new function allows us to model moving peaks in the temperature profile, a feature that previous heat balance methods cannot capture. If the boundary temperature varies so that at some time t & 0 it equals the far-field temperature, then standard methods predict that the temperature is everywhere at this constant value. The new method predicts the correct behaviour. It is also shown that this function provides even more accurate results, when coupled with the new CIM, than the polynomial profile. Analysis primarily focuses on a specified constant boundary temperature and is then extended to constant flux, Newton cooling and time dependent boundary conditions.

Monitoring time-dependent degradation of phospholipids in sectioned tissues by MALDI imaging mass spectrometry.

Imaging mass spectrometry (IMS) is useful for visualizing the localization of phospholipids on biological tissue surfaces creating great opportunities for IMS in lipidomic investigations. With advancements in IMS of lipids, there is a demand for large-scale tissue studies necessitating stable, efficient and well-defined sample handling procedures. Our work within this article shows the effects of different storage conditions on the phospholipid composition of sectioned tissues from mouse organs. We have taken serial sections from mouse brain, kidney and liver thaw mounted unto ITO-coated glass slides and stored them under various conditions later analyzing them at fixed time points. A global decrease in phospholipid signal intensity is shown to occur and to be a function of time and temperature. Contrary to the global decrease, oxidized phospholipid and lysophospholipid species are found to increase within 2 h and 24 h, respectively, when mounted sections are kept at ambient room conditions. Imaging experiments reveal that degradation products increase globally across the tissue. Degradation is shown to be inhibited by cold temperatures, with sample integrity maintained up to a week after storage in −80 °C freezer under N2 atmosphere. Overall, the results demonstrate a timeline of the effects of lipid degradation specific to sectioned tissues and provide several lipid species which can serve as markers of degradation. Importantly, the timeline demonstrates oxidative sample degradation begins appearing within the normal timescale of IMS sample preparation of lipids (i.e. 1-2 h) and that long-term degradation is global. Taken together, these results strengthen the notion that standardized procedures are required for phospholipid IMS of large sample sets, or in studies where many serial sections are prepared together but analyzed over time such as in 3-D IMS reconstruction experiments.

Fluid flow regulates stromal cell organization and CCL21 expression in a tissue-engineered lymph node microenvironment.

In the paracortex of the lymph node (LN), T zone fibroblastic reticular cells (TRCs) orchestrate an immune response by guiding lymphocyte migration both physically, by creating three-dimensional (3D) cell networks, and chemically, by secreting the chemokines CCL19 and CCL21 that direct interactions between CCR7-expressing cells, including mature dendritic cells and naive T cells. TRCs also enwrap matrix-based conduits that transport fluid from the subcapsular sinus to high endothelial venules, and fluid flow through the draining LN rapidly increases upon tissue injury or inflammation. To determine whether fluid flow affects TRC organization or function within a 3D network, we regenerated the 3D LN T zone stromal network by culturing murine TRC clones within a macroporous polyurethane scaffold containing type I collagen and Matrigel and applying slow interstitial flow (1-23 microm/min). We show that the 3D environment and slow interstitial flow are important regulators of TRC morphology, organization, and CCL21 secretion. Without flow, CCL21 expression could not be detected. Furthermore, when flow through the LN was blocked in mice in vivo, CCL21 gene expression was down-regulated within 2 h. These results highlight the importance of lymph flow as a homeostatic regulator of constitutive TRC activity and introduce the concept that increased lymph flow may act as an early inflammatory cue to enhance CCL21 expression by TRCs, thereby ensuring efficient immune cell trafficking, lymph sampling, and immune response induction.

Preferential transcription of Paracoccidioides brasiliensis genes: host niche and time-dependent expression

Paracoccidioides brasiliensis causes infection through inhalation by the host of airborne propagules from the mycelium phase of the fungus. This fungus reaches the lungs, differentiates into the yeast form and is then disseminated to virtually all parts of the body. Here we review the identification of differentially-expressed genes in host-interaction conditions. These genes were identified by analyzing expressed sequence tags (ESTs) from P. brasiliensis cDNA libraries. The P. brasiliensis was recovered from infected mouse liver as well as from fungal yeast cells incubated in human blood and plasma, mimicking fungal dissemination to organs and tissues and sites of infection with inflammation, respectively. In addition, ESTs from a cDNA library of P. brasiliensis mycelium undergoing the transition to yeast were previously analyzed. Together, these studies reveal significant changes in the expression of a number of genes of potential importance in the host-fungus interaction. In addition, the unique and divergent representation of transcripts when the cDNA libraries are compared suggests differential gene expression in response to specific niches in the host. This analysis of gene expression patterns provides details about host-pathogen interactions and peculiarities of sites within the host.

Profile order and time-dependent artifacts in contrast-enhanced coronary MR angiography at 3T: origin and prevention.

To enhance the clinical value of coronary magnetic resonance angiography (MRA), high-relaxivity contrast agents have recently been used at 3T. Here we examine a uniform bilateral shadowing artifact observed along the coronary arteries in MRA images collected using such a contrast agent. Simulations were performed to characterize this artifact, including its origin, to determine how best to mitigate this effect, and to optimize a data acquisition/injection scheme. An intraluminal contrast agent concentration model was used to simulate various acquisition strategies with two profile orders for a slow-infusion of a high-relaxivity contrast agent. Filtering effects from temporally variable weighting in k-space are prominent when a centric, radial (CR) profile order is applied during contrast infusion, resulting in decreased signal enhancement and underestimation of vessel width, while both pre- and postinfusion steady-state acquisitions result in overestimation of the vessel width. Acquisition during the brief postinfusion steady-state produces the greatest signal enhancement and minimizes k-space filtering artifacts.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.