Identifying and Imposing Partial Equilibrium in Chemically Reacting Systems

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76236/1/AIAA-2003-667-389.pd

An implementation of a grid-independent upwind scheme for the Euler equations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/77364/1/AIAA-1989-1931-789.pd

A Magnetohydrodynamic Model of the M87 Jet I: Superluminal Knot Ejections from HST-1 as Trails of Quad Relativistic MHD Shocks

This is the first in a series of papers that introduces a new paradigm for understanding the jet in M87: a collimated relativistic flow in which strong magnetic fields play a dominant dynamical role. Here wefocus on the flow downstream of HST-1 - an essentially stationary flaring feature that ejects trails of superluminal components. We propose that these components are quad relativistic magnetohydrodynamic shock fronts (forward/reverse fast and slow modes) in a narrow jet with a helically twisted magnetic structure. And we demonstrate the properties of such shocks with simple one-dimensional numerical simulations. Quasi-periodic ejections of similar component trails may be responsible for the M87 jet substructures observed further downstream on 100 - 1,000 pc scales. This new paradigm requires the assimilation of some new concepts into the astrophysical jet community, particularly the behavior of slow/fast-mode waves/shocks and of current-driven helical kink instabilities. However, the prospects of these ideas applying to a large number of other jet systems may make this worth the effort.Comment: 7 pages, 4 figures, Accepted for Publication in Ap

Three-Dimensional Simulations of Jets from Keplerian Disks: Self--Regulatory Stability

We present the extension of previous two-dimensional simulations of the time-dependent evolution of non-relativistic outflows from the surface of Keplerian accretion disks, to three dimensions. The accretion disk itself is taken to provide a set of fixed boundary conditions for the problem. The 3-D results are consistent with the theory of steady, axisymmetric, centrifugally driven disk winds up to the Alfv\'en surface of the outflow. Beyond the Alfv\'en surface however, the jet in 3-D becomes unstable to non-axisymmetric, Kelvin-Helmholtz instabilities. We show that jets maintain their long-term stability through a self-limiting process wherein the average Alfv\'enic Mach number within the jet is maintained to order unity. This is accomplished in at least two ways. First, poloidal magnetic field is concentrated along the central axis of the jet forming a ``backbone'' in which the Alfv\'en speed is sufficiently high to reduce the average jet Alfv\'enic Mach number to unity. Second, the onset of higher order Kelvin-Helmholtz ``flute'' modes (m \ge 2) reduce the efficiency with which the jet material is accelerated, and transfer kinetic energy of the outflow into the stretched, poloidal field lines of the distorted jet. This too has the effect of increasing the Alfv\'en speed, and thus reducing the Alfv\'enic Mach number. The jet is able to survive the onset of the more destructive m=1 mode in this way. Our simulations also show that jets can acquire corkscrew, or wobbling types of geometries in this relatively stable end-state, depending on the nature of the perturbations upon them. Finally, we suggest that jets go into alternating periods of low and high activity as the disappearance of unstable modes in the sub-Alfv\'enic regime enables another cycle of acceleration to super-Alfv\'enic speeds.Comment: 57 pages, 22 figures, submitted to Ap

A comparison of Eulerian and Lagrangian transport and non-linear reaction algorithms

When laboratory-measured chemical reaction rates are used in simulations at the field-scale, the models typically overpredict the apparent reaction rates. The discrepancy is primarily due to poorer mixing of chemically distinct waters at the larger scale. As a result, realistic field-scale predictions require accurate simulation of the degree of mixing between fluids. The Lagrangian particle-tracking (PT) method is a now-standard way to simulate the transport of conservative or sorbing solutes. The method’s main advantage is the absence of numerical dispersion (and its artificial mixing) when simulating advection. New algorithms allow particles of different species to interact in nonlinear (e.g., bimolecular) reactions. Therefore, the PT methods hold a promise of more accurate field-scale simulation of reactive transport because they eliminate the masking effects of spurious mixing due to advection errors inherent in grid-based methods. A hypothetical field-scale reaction scenario is constructed and run in PT and Eulerian (finite-volume/finite-difference) simulators. Grid-based advection schemes considered here include 1st- to 3rd-order spatially accurate total-variation-diminishing flux-limiting schemes, both of which are widely used in current transport/reaction codes. A homogeneous velocity field in which the Courant number is everywhere unity, so that the chosen Eulerian methods incur no error when simulating advection, shows that both the Eulerian and PT methods can achieve convergence in the L1 (integrated concentration) norm, but neither shows stricter pointwise convergence. In this specific case with a constant dispersion coefficient and bimolecular reaction A+B¿P, the correct total amount of product is 0.221MA0, where MA0 is the original mass of reactant A. When the Courant number drops, the grid-based simulations can show remarkable errors due to spurious over- and under-mixing. In a heterogeneous velocity field (keeping the same constant and isotropic dispersion), the PT simulations show an increased reaction total from 0.221MA0 to 0.372MA0 due to fluid deformation, while the 1st-order Eulerian simulations using ˜ 106 cells (with a classical grid Peclet number ¿x/aL of 10) have total product of 0.53MA0, or approximately twice as much additional reaction due to advection error. The 3rd-order TVD algorithm fares better, with total product of 0.394MA0, or about 1.14 times the increased reaction total. A very strict requirement on grid Peclet numbers for Eulerian simulations will be required for realistic reactions because of their nonlinear nature. We analytically estimate the magnitude of the effect for the end-member cases of very fast and very slow reactions and show that in either case, the mass produced is proportional to View the MathML source where Pe is the Peclet number. Therefore, extra mass is produced according to View the MathML source where the dispersion includes any numerical dispersion error. We test two PT methods, one that kills particles upon reaction and another that decrements a particle’s mass. For the bimolecular reaction studied here, the computational demands of the particle-killing methods are much smaller than, and the particle-number-preserving algorithm are on par with, the fastest Eulerian methods.Peer ReviewedPostprint (author's final draft

Process-evaluation of tropospheric humidity simulated by general circulation models using water vapor isotopic observations: 2. Using isotopic diagnostics to understand the mid and upper tropospheric moist bias in the tropics and subtropics

Evaluating the representation of processes controlling tropical and subtropical tropospheric relative humidity (RH) in atmospheric general circulation models (GCMs) is crucial to assess the credibility of predicted climate changes. GCMs have long exhibited a moist bias in the tropical and subtropical mid and upper troposphere, which could be due to the mis-representation of cloud processes or of the large-scale circulation, or to excessive diffusion during water vapor transport. The goal of this study is to use observations of the water vapor isotopic ratio to understand the cause of this bias. We compare the three-dimensional distribution of the water vapor isotopic ratio measured from space and ground to that simulated by several versions of the isotopic GCM LMDZ. We show that the combined evaluation of RH and of the water vapor isotopic composition makes it possible to discriminate the most likely cause of RH biases. Models characterized either by an excessive vertical diffusion, an excessive convective detrainment or an underestimated in situ cloud condensation will all produce a moist bias in the free troposphere. However, only an excessive vertical diffusion can lead to a reversed seasonality of the free tropospheric isotopic composition in the subtropics compared to observations. Comparing seven isotopic GCMs suggests that the moist bias found in many GCMs in the mid and upper troposphere most frequently results from an excessive diffusion during vertical water vapor transport. This study demonstrates the added value of water vapor isotopic measurements for interpreting shortcomings in the simulation of RH by climate models

Recommendations for enterovirus diagnostics and characterisation within and beyond Europe.

Enteroviruses (EV) can cause severe neurological and respiratory infections, and occasionally lead to devastating outbreaks as previously demonstrated with EV-A71 and EV-D68 in Europe. However, these infections are still often underdiagnosed and EV typing data is not currently collected at European level. In order to improve EV diagnostics, collate data on severe EV infections and monitor the circulation of EV types, we have established European non-polio enterovirus network (ENPEN). First task of this cross-border network has been to ensure prompt and adequate diagnosis of these infections in Europe, and hence we present recommendations for non-polio EV detection and typing based on the consensus view of this multidisciplinary team including experts from over 20 European countries. We recommend that respiratory and stool samples in addition to cerebrospinal fluid (CSF) and blood samples are submitted for EV testing from patients with suspected neurological infections. This is vital since viruses like EV-D68 are rarely detectable in CSF or stool samples. Furthermore, reverse transcriptase PCR (RT-PCR) targeting the 5'noncoding regions (5'NCR) should be used for diagnosis of EVs due to their sensitivity, specificity and short turnaround time. Sequencing of the VP1 capsid protein gene is recommended for EV typing; EV typing cannot be based on the 5'NCR sequences due to frequent recombination events and should not rely on virus isolation. Effective and standardized laboratory diagnostics and characterisation of circulating virus strains are the first step towards effective and continuous surveillance activities, which in turn will be used to provide better estimation on EV disease burden

Time to Switch to Second-line Antiretroviral Therapy in Children With Human Immunodeficiency Virus in Europe and Thailand.

Background: Data on durability of first-line antiretroviral therapy (ART) in children with human immunodeficiency virus (HIV) are limited. We assessed time to switch to second-line therapy in 16 European countries and Thailand. Methods: Children aged <18 years initiating combination ART (≥2 nucleoside reverse transcriptase inhibitors [NRTIs] plus nonnucleoside reverse transcriptase inhibitor [NNRTI] or boosted protease inhibitor [PI]) were included. Switch to second-line was defined as (i) change across drug class (PI to NNRTI or vice versa) or within PI class plus change of ≥1 NRTI; (ii) change from single to dual PI; or (iii) addition of a new drug class. Cumulative incidence of switch was calculated with death and loss to follow-up as competing risks. Results: Of 3668 children included, median age at ART initiation was 6.1 (interquartile range (IQR), 1.7-10.5) years. Initial regimens were 32% PI based, 34% nevirapine (NVP) based, and 33% efavirenz based. Median duration of follow-up was 5.4 (IQR, 2.9-8.3) years. Cumulative incidence of switch at 5 years was 21% (95% confidence interval, 20%-23%), with significant regional variations. Median time to switch was 30 (IQR, 16-58) months; two-thirds of switches were related to treatment failure. In multivariable analysis, older age, severe immunosuppression and higher viral load (VL) at ART start, and NVP-based initial regimens were associated with increased risk of switch. Conclusions: One in 5 children switched to a second-line regimen by 5 years of ART, with two-thirds failure related. Advanced HIV, older age, and NVP-based regimens were associated with increased risk of switch

Development and Validation of a Risk Score for Chronic Kidney Disease in HIV Infection Using Prospective Cohort Data from the D:A:D Study

Ristola M. on työryhmien DAD Study Grp ; Royal Free Hosp Clin Cohort ; INSIGHT Study Grp ; SMART Study Grp ; ESPRIT Study Grp jäsen.Background Chronic kidney disease (CKD) is a major health issue for HIV-positive individuals, associated with increased morbidity and mortality. Development and implementation of a risk score model for CKD would allow comparison of the risks and benefits of adding potentially nephrotoxic antiretrovirals to a treatment regimen and would identify those at greatest risk of CKD. The aims of this study were to develop a simple, externally validated, and widely applicable long-term risk score model for CKD in HIV-positive individuals that can guide decision making in clinical practice. Methods and Findings A total of 17,954 HIV-positive individuals from the Data Collection on Adverse Events of Anti-HIV Drugs (D:A:D) study with >= 3 estimated glomerular filtration rate (eGFR) values after 1 January 2004 were included. Baseline was defined as the first eGFR > 60 ml/min/1.73 m2 after 1 January 2004; individuals with exposure to tenofovir, atazanavir, atazanavir/ritonavir, lopinavir/ritonavir, other boosted protease inhibitors before baseline were excluded. CKD was defined as confirmed (>3 mo apart) eGFR In the D:A:D study, 641 individuals developed CKD during 103,185 person-years of follow-up (PYFU; incidence 6.2/1,000 PYFU, 95% CI 5.7-6.7; median follow-up 6.1 y, range 0.3-9.1 y). Older age, intravenous drug use, hepatitis C coinfection, lower baseline eGFR, female gender, lower CD4 count nadir, hypertension, diabetes, and cardiovascular disease (CVD) predicted CKD. The adjusted incidence rate ratios of these nine categorical variables were scaled and summed to create the risk score. The median risk score at baseline was -2 (interquartile range -4 to 2). There was a 1: 393 chance of developing CKD in the next 5 y in the low risk group (risk score = 5, 505 events), respectively. Number needed to harm (NNTH) at 5 y when starting unboosted atazanavir or lopinavir/ritonavir among those with a low risk score was 1,702 (95% CI 1,166-3,367); NNTH was 202 (95% CI 159-278) and 21 (95% CI 19-23), respectively, for those with a medium and high risk score. NNTH was 739 (95% CI 506-1462), 88 (95% CI 69-121), and 9 (95% CI 8-10) for those with a low, medium, and high risk score, respectively, starting tenofovir, atazanavir/ritonavir, or another boosted protease inhibitor. The Royal Free Hospital Clinic Cohort included 2,548 individuals, of whom 94 individuals developed CKD (3.7%) during 18,376 PYFU (median follow-up 7.4 y, range 0.3-12.7 y). Of 2,013 individuals included from the SMART/ESPRIT control arms, 32 individuals developed CKD (1.6%) during 8,452 PYFU (median follow-up 4.1 y, range 0.6-8.1 y). External validation showed that the risk score predicted well in these cohorts. Limitations of this study included limited data on race and no information on proteinuria. Conclusions Both traditional and HIV-related risk factors were predictive of CKD. These factors were used to develop a risk score for CKD in HIV infection, externally validated, that has direct clinical relevance for patients and clinicians to weigh the benefits of certain antiretrovirals against the risk of CKD and to identify those at greatest risk of CKD.Peer reviewe