Quantitative Analysis of Mechanisms That Govern Red Blood Cell Age Structure and Dynamics during Anaemia

Mathematical modelling has proven an important tool in elucidating and quantifying mechanisms that govern the age structure and population dynamics of red blood cells (RBCs). Here we synthesise ideas from previous experimental data and the mathematical modelling literature with new data in order to test hypotheses and generate new predictions about these mechanisms. The result is a set of competing hypotheses about three intrinsic mechanisms: the feedback from circulating RBC concentration to production rate of immature RBCs (reticulocytes) in bone marrow, the release of reticulocytes from bone marrow into the circulation, and their subsequent ageing and clearance. In addition we examine two mechanisms specific to our experimental system: the effect of phenylhydrazine (PHZ) and blood sampling on RBC dynamics. We performed a set of experiments to quantify the dynamics of reticulocyte proportion, RBC concentration, and erythropoietin concentration in PHZ-induced anaemic mice. By quantifying experimental error we are able to fit and assess each hypothesis against our data and recover parameter estimates using Markov chain Monte Carlo based Bayesian inference. We find that, under normal conditions, about 3% of reticulocytes are released early from bone marrow and upon maturation all cells are released immediately. In the circulation, RBCs undergo random clearance but have a maximum lifespan of about 50 days. Under anaemic conditions reticulocyte production rate is linearly correlated with the difference between normal and anaemic RBC concentrations, and their release rate is exponentially correlated with the same. PHZ appears to age rather than kill RBCs, and younger RBCs are affected more than older RBCs. Blood sampling caused short aperiodic spikes in the proportion of reticulocytes which appear to have a different developmental pathway than normal reticulocytes. We also provide evidence of large diurnal oscillations in serum erythropoietin levels during anaemia

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Modelling of the effect of ELMs on fuel retention at the bulk W divertor of JET

Effect of ELMs on fuel retention at the bulk W target of JET ITER-Like Wall was studied with multi-scale calculations. Plasma input parameters were taken from ELMy H-mode plasma experiment. The energetic intra-ELM fuel particles get implanted and create near-surface defects up to depths of few tens of nm, which act as the main fuel trapping sites during ELMs. Clustering of implantation-induced vacancies were found to take place. The incoming flux of inter-ELM plasma particles increases the different filling levels of trapped fuel in defects. The temperature increase of the W target during the pulse increases the fuel detrapping rate. The inter-ELM fuel particle flux refills the partially emptied trapping sites and fills new sites. This leads to a competing effect on the retention and release rates of the implanted particles. At high temperatures the main retention appeared in larger vacancy clusters due to increased clustering rate

Impact of fast ions on density peaking in JET : fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position rho(t) = 0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile

Impact of fast ions on density peaking in JET: fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/ trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position ρ$_{t}$=0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile