Showing 1–27 of 27 results for author: Duval, B

Expulsion of runaway electrons using ECRH in the TCV tokamak

Authors: J. Decker, M. Hoppe, U. Sheikh, B. P. Duval, G. Papp, L. Simons, T. Wijkamp, J. Cazabonne, S. Coda, E. Devlaminck, O. Ficker, R. Hellinga, U. Kumar, Y. Savoye-Peysson, L. Porte, C. Reux, C. Sommariva, A. Tema Biwolé, B. Vincent, L. Votta, the TCV Team, the EUROfusion Tokamak Exploitation Team

Abstract: Runaway electrons (REs) are a concern for tokamak fusion reactors from discharge startup to termination. A sudden localized loss of a multi-megaampere RE beam can inflict severe damage to the first wall. Should a disruption occur, the existence of a RE seed may play a significant role in the formation of a RE beam and the magnitude of its current. The application of central electron cyclotron reso… ▽ More Runaway electrons (REs) are a concern for tokamak fusion reactors from discharge startup to termination. A sudden localized loss of a multi-megaampere RE beam can inflict severe damage to the first wall. Should a disruption occur, the existence of a RE seed may play a significant role in the formation of a RE beam and the magnitude of its current. The application of central electron cyclotron resonance heating (ECRH) in the Tokamak à Configuration Variable (TCV) reduces an existing RE seed population by up to three orders of magnitude within only a few hundred milliseconds. Applying ECRH before a disruption can also prevent the formation of a post-disruption RE beam in TCV where it would otherwise be expected. The RE expulsion rate and consequent RE current reduction are found to increase with applied ECRH power. Whereas central ECRH is effective in expelling REs, off-axis ECRH has a comparatively limited effect. A simple 0-D model for the evolution of the RE population is presented that explains the effective ECRH-induced RE expulsion results from the combined effects of increased electron temperature and enhanced RE transport. △ Less

Submitted 22 July, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Investigating the influence of divertor baffles on nitrogen-seeded detachment in TCV with SOLPS-ITER simulations and TCV experiments

Authors: G. Sun, H. Reimerdes, C. Theiler, B. P. Duval, M. Carpita, C. Colandrea, R. Ducker, O. Fevrier, S. Gorno, L. Simons, E. Tonello

Abstract: Plasma edge simulations with the SOLPS-ITER code are performed to study the influence of divertor baffles on nitrogen-seeded detachment in TCV single-null, L-mode discharges. Scans of nitrogen seeding rate are conducted in both baffled and unbaffled TCV divertors, where the nitrogen seeding with baffles is found to yield lower target temperatures and heat fluxes than with baffles-only and with see… ▽ More Plasma edge simulations with the SOLPS-ITER code are performed to study the influence of divertor baffles on nitrogen-seeded detachment in TCV single-null, L-mode discharges. Scans of nitrogen seeding rate are conducted in both baffled and unbaffled TCV divertors, where the nitrogen seeding with baffles is found to yield lower target temperatures and heat fluxes than with baffles-only and with seeding-only. The cumulative effects of baffles and seeding on target parameters are explained by the two-point model. The divertor neutral density and neutral compression increase with baffles, due to lower divertor to main chamber neutral conductance, as explained by a schematic neutral transport model with baffles. The nitrogen retention, defined as the ratio of average nitrogen nuclei density in divertor and main chamber, increases with the seeding rate if baffled, and remains constant if unbaffled. At the same outboard mid-plane separatrix plasma density, the nitrogen retention with baffles is lower than the unbaffled retention at low seeding levels and is higher at high seeding levels, which is explained by the changes of nitrogen ion and neutral transport with baffles and seeding. The baffled carbon retention is higher than the unbaffled retention due to lower divertor to main chamber carbon neutral conductance. Baffles increase the divertor radiation. The predicted trends of target parameters, the distribution of neutrals and radiations are well supported by TCV experiments, though discrepancies in the absolute values remain. The simulations yield an overall colder and denser divertor, consistent with previous SOLPS-ITER simulations of Ohmically heated L-modes in TCV. The successful comparison of simulation and experiment, together with the understanding gained from the neutral transport model, increases the confidence in the SOLPS simulations for the next TCV divertor upgrade. △ Less

Submitted 14 October, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comparison of detachment in Ohmic plasmas with positive and negative triangularity

Authors: O. Février, C. K. Tsui, G. Durr-Legoupil-Nicoud, C. Theiler, M. Carpita, S. Coda, C. Colandrea, B. P. Duval, S. Gorno, E. Huett, B. Linehan, A. Perek, L. Porte, H. Reimerdes, O. Sauter, E. Tonello, M. Zurita, T. Bolzonella, F. Sciortino, the TCV Team, the EUROfusion Tokamak Exploitation Team

Abstract: In recent years, negative triangularity (NT) has emerged as a potential high-confinement L-mode reactor solution. In this work, detachment is investigated using core density ramps in lower single null Ohmic L-mode plasmas across a wide range of upper, lower, and average triangularity (the mean of upper and lower triangularity: $δ$) in the TCV tokamak. It is universally found that detachment is mor… ▽ More In recent years, negative triangularity (NT) has emerged as a potential high-confinement L-mode reactor solution. In this work, detachment is investigated using core density ramps in lower single null Ohmic L-mode plasmas across a wide range of upper, lower, and average triangularity (the mean of upper and lower triangularity: $δ$) in the TCV tokamak. It is universally found that detachment is more difficult to access for NT shaping. The outer divertor leg of discharges with $δ\approx -0.3$ could not be cooled to below $5~\mathrm{eV}$ through core density ramps alone. The behavior of the upstream plasma and geometrical divertor effects (e.g. a reduced connection length with negative lower triangularity) do not fully explain the challenges in detaching NT plasmas. Langmuir probe measurements of the target heat flux widths ($λ_q$) were constant to within 30% across an upper triangularity scan, while the spreading factor $S$ was lower by up to 50% for NT, indicating a generally lower integral Scrape-Off Layer width, $λ_{int}$. The line-averaged core density was typically higher for NT discharges for a given fuelling rate, possibly linked to higher particle confinement in NT. Conversely, the divertor neutral pressure and integrated particle fluxes to the targets were typically lower for the same line-averaged density, indicating that NT configurations may be closer to the sheath-limited regime than their PT counterparts, which may explain why NT is more challenging to detach. △ Less

Submitted 23 January, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Parallel flows as a key component to interpret Super-X divertor experiments

Authors: M. Carpita, O. Février, H. Reimerdes, C. Theiler, B. P. Duval, C. Colandrea, G. Durr-Legoupil-Nicoud, D. Galassi, S. Gorno, E. Huett, J. Loizu, L. Martinelli, A. Perek, L. Simons, G. Sun, E. Tonello, C. Wüthrich, the TCV team

Abstract: The Super-X Divertor (SXD) is an alternative divertor configuration leveraging total flux expansion at the Outer Strike Point (OSP). While the extended 2-Point Model (2PM) predicts facilitated detachment access and control in the SXD configuration, these attractive features are not always retrieved experimentally. These discrepancies are at least partially explained by the effect of parallel flows… ▽ More The Super-X Divertor (SXD) is an alternative divertor configuration leveraging total flux expansion at the Outer Strike Point (OSP). While the extended 2-Point Model (2PM) predicts facilitated detachment access and control in the SXD configuration, these attractive features are not always retrieved experimentally. These discrepancies are at least partially explained by the effect of parallel flows which, when self-consistently included in the 2PM, reveal the role of total flux expansion on the pressure balance and weaken the total flux expansion effect on detachment access and control, compared to the original predictions. This new model can partially explain the discrepancies between the 2PM and experiments performed on tokamak à configuration variable (TCV), in ohmic L-mode scenarios, which are particularly apparent when scanning the OSP major radius Rt. In core density ramps in lower Single-Null (SN) configuration, the impact of Rt on the CIII emission front movement in the divertor outer leg - used as a proxy for the plasma temperature in the divertor - is substantially weaker than 2PM predictions. Furthermore, in OSP radial sweeps in lower and upper SN configurations, in ohmic L-mode scenarios with a constant core density, the peak parallel particle flux density at the OSP is almost independent of Rt, while the 2PM predicts a linear dependence. Finally, analytical and numerical modeling of parallel flows in the divertor is presented. It is shown that an increase in total flux expansion can favour supersonic flows at the OSP. Parallel flows are also shown to be relevant by analysing SOLPS-ITER simulations of TCV. △ Less

Submitted 27 February, 2024; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: This version of the paper is outdated. The most recent, and published, version of the paper can be found at the corresponding Journal / DOI reference

Journal ref: Nucl. Fusion 64 046019 (2024)

Performance assessment of a tightly baffled, long-legged divertor configuration in TCV with SOLPS-ITER

Authors: G. Sun, H. Reimerdes, C. Theiler, B. P. Duval, M. Carpita, C. Colandrea, O. Février, the TCV team

Abstract: Numerical simulations explore the possibility to test the tightly baffled, long-legged divertor (TBLLD) concept in a future upgrade of the Tokamak à configuration variable (TCV). The SOLPS-ITER code package is used to compare the exhaust performance of several TBLLD configurations with existing unbaffled and baffled TCV configurations. The TBLLDs feature a range of radial gaps between the separatr… ▽ More Numerical simulations explore the possibility to test the tightly baffled, long-legged divertor (TBLLD) concept in a future upgrade of the Tokamak à configuration variable (TCV). The SOLPS-ITER code package is used to compare the exhaust performance of several TBLLD configurations with existing unbaffled and baffled TCV configurations. The TBLLDs feature a range of radial gaps between the separatrix and the outer leg side walls. All considered TBLLDs are predicted to lead to a denser and colder plasma in front of the targets and improve the power handling by factors of 2-3 compared to the present, baffled divertor and by up to a factor of 12 compared to the original, unbaffled configuration. The improved TBLLD performance is mainly due to a better neutral confinement with improved plasma-neutral interactions in the divertor region. Both power handling capability and neutral confinement increases when reducing the radial gap. The core compatibility of TBLLDs with nitrogen seeding is also evaluated and the detachment window with acceptable core pollution for the proposed TBLLDs is explored, showing a reduction of required upstream impurity concentration up to 18% to achieve the detachment with thinner radial gap. △ Less

Submitted 20 April, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Investigating the impact of the molecular charge-exchange rate on detached SOLPS-ITER simulations

Authors: Kevin Verhaegh, Aelwyn C Williams, David Moulton, Bruce Lipschultz, Basil P. Duval, Olivier Fevrier, Alexandre Fil, Nick Osborne, Holger Reimerdes, Christian Theiler

Abstract: Plasma-molecular interactions generate molecular ions which react with the plasma and contribute to detachment through molecular activated recombination (MAR), reducing the ion target flux, and molecular activated dissociation (MAD), both of which create excited atoms. Hydrogenic emission from these atoms have been detected experimentally in detached TCV, JET and MAST-U deuterium plasmas. The TCV… ▽ More Plasma-molecular interactions generate molecular ions which react with the plasma and contribute to detachment through molecular activated recombination (MAR), reducing the ion target flux, and molecular activated dissociation (MAD), both of which create excited atoms. Hydrogenic emission from these atoms have been detected experimentally in detached TCV, JET and MAST-U deuterium plasmas. The TCV findings, however, were in disagreement with SOLPS-ITER simulations for deuterium indicating a molecular ion density ($D_2^+$) that was insufficient to lead to significant hydrogenic emission, which was attributed to underestimates of the molecular charge exchange rate ($D_2 + D^+ \rightarrow D_2^+ + D$) for deuterium (obtained by rescaling the hydrogen rates by their isotope mass). In this work, we have performed new SOLPS-ITER simulations with the default rate setup and a modified rate setup where ion isotope mass rescaling was disabled. This increased the $D_2^+$ content by $> \times 100$. By disabling ion isotope mass rescaling: 1) the total ion sinks are more than doubled due to the inclusion of MAR; 2) the additional MAR causes the ion target flux to roll-over during detachment; 3) the total $Dα$ emission in the divertor increases during deep detachment by roughly a factor four; 4) the neutral atom density in the divertor is doubled due to MAD, leading to a 50\% increase in neutral pressure; 5) total hydrogenic power loss is increased by up to 60\% due to MAD. These differences result in an improved agreement between the experiment and the simulations in terms of spectroscopic measurements, ion source/sink inferences and the occurrence of an ion target flux roll-over. △ Less

Submitted 16 April, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Power exhaust and core-divertor compatibility of the baffled snowflake divertor in TCV

Authors: Sophie Gorno, Claudia Colandrea, Olivier Février, Holger Reimerdes, Christian Theiler, Basil P Duval, Tilmann Lunt, Harshita Raj, Umar A Sheikh, Luke Simons, Andrew Thornton

Abstract: A baffled Snowflake Minus Low-Field Side (SF-LFS) is geometrically-optimised in TCV, increasing divertor neutral pressure, to evaluate the roles of divertor closure (comparing with an unbaffled SF-LFS) and magnetic geometry (comparing with a baffled Single Null, SN) in power exhaust and core-divertor compatibility. Ohmically-heated L-mode discharges in deuterium, with a line-averaged core density… ▽ More A baffled Snowflake Minus Low-Field Side (SF-LFS) is geometrically-optimised in TCV, increasing divertor neutral pressure, to evaluate the roles of divertor closure (comparing with an unbaffled SF-LFS) and magnetic geometry (comparing with a baffled Single Null, SN) in power exhaust and core-divertor compatibility. Ohmically-heated L-mode discharges in deuterium, with a line-averaged core density of approximately 4x10(19) m-3, are seeded with nitrogen to approach detached conditions. Baffles in the SF-LFS configuration are found to reduce the peak outer target heat flux by up to 23%, without significantly affecting the location of the inter-null radiation region or the core-divertor compatibility. When compared to the baffled SN, the baffled SF-LFS exhibits a reduction in outer target heat flux by up to 66% and the ability to balance the strike-point distribution of heat flux. These benefits are less significant with N2 seeding, with similar peak target quantities (such as heat flux, electron temperature and ion flux) and divertor radiated power. Despite a radiating region located farther from the confined plasma for the SF-LFS than the baffled SN, no change in core confinement is observed. Core effective charge even indicates an increase in core impurity penetration for the SF-LFS. These experiments constitute a good reference for detailed model validations and extrapolations, exploring important physics such as core impurity shielding and the dependence of divertor cross-field transport on magnetic geometry. △ Less

Submitted 23 November, 2022; v1 submitted 7 September, 2022; originally announced September 2022.

Improved Heat and Particle Flux Mitigation in High Core Confinement, Baffled, Alternate Divertor Configurations in the TCV tokamak

Authors: Harshita Raj, C. Theiler, A. Thornton, O. Fevrier, S. Gorno, F. Bagnato, P. Blanchard, C. Colandrea, H. de Oliveira, B. P. Duval, B. Labit, A. Perek, H. Reimerdes, U. Sheikh, M. Vallar, B. Vincent

Abstract: Nitrogen seeded detachment has been achieved in the Tokamak a Configuration Variable (TCV) in advanced divertor configurations (ADCs), namely X-divertor and X-point target, with and without baffles in H-mode plasmas with high core confinement. Both ADCs show a remarkable reduction in the inter-ELM particle and heat fluxes to the target compared to the standard divertor configuration. 95-98% of the… ▽ More Nitrogen seeded detachment has been achieved in the Tokamak a Configuration Variable (TCV) in advanced divertor configurations (ADCs), namely X-divertor and X-point target, with and without baffles in H-mode plasmas with high core confinement. Both ADCs show a remarkable reduction in the inter-ELM particle and heat fluxes to the target compared to the standard divertor configuration. 95-98% of the peak heat flux to the target is mitigated as a synergetic effect of ADCs, baffling, and nitrogen seeded detachment. The effect of divertor geometry and baffles on core-divertor compatibility is investigated in detail. The power balance in these experiments is also investigated to explore the physics behind the observed reduction in heat fluxes in the ADCs. △ Less

Submitted 27 July, 2022; originally announced July 2022.

New insights on divertor parallel flows, ExB drifts, and fluctuations from in situ, two-dimensional probe measurement in the Tokamak à Configuration Variable

Authors: H. De Oliveira, C. Theiler, O. Février, H. Reimerdes, B. P. Duval, C. K. Tsui, S. Gorno, D. S. Oliveira, A. Perek

Abstract: In-situ, two-dimensional (2D) Langmuir probe measurements across a large part of the TCV divertor are reported in L-mode discharges with and without divertor baffles. This provides detailed insights into time averaged profiles, particle fluxes, and fluctuations behavior in different divertor regimes. The presence of the baffles is shown to substantially increase the divertor neutral pressure for a… ▽ More In-situ, two-dimensional (2D) Langmuir probe measurements across a large part of the TCV divertor are reported in L-mode discharges with and without divertor baffles. This provides detailed insights into time averaged profiles, particle fluxes, and fluctuations behavior in different divertor regimes. The presence of the baffles is shown to substantially increase the divertor neutral pressure for a given upstream density and to facilitate the access to detachment, an effect that increases with plasma current. The detailed, 2D probe measurements allow for a divertor particle balance, including ion flux contributions from parallel flows and ExB drifts. The poloidal flux contribution from the latter is often comparable or even larger than the former, such that the divertor parallel flow direction reverses in some conditions, pointing away from the target. In most conditions, the integrated particle flux at the outer target can be predominantly ascribed to ionization along the outer divertor leg, consistent with a closed-box approximation of the divertor. The exception is a strongly detached divertor, achieved here only with baffles, where the total poloidal ion flux even decreases towards the outer target, indicative of significant plasma recombination. The most striking observation from relative density fluctuation measurements along the outer divertor leg is the transition from poloidally uniform fluctuation levels in attached conditions to fluctuations strongly peaking near the X-point when approaching detachment. △ Less

Submitted 18 May, 2022; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: Added reference

Spectroscopic investigations of detachment on the MAST Upgrade Super-X divertor

Authors: Kevin Verhaegh, Bruce Lipschultz, James Harrison, Nick Osborne, Aelwyn Williams, Peter Ryan, James Clark, Fabio Federici, Bob Kool, Tijs Wijkamp, Alexandre Fil, David Moulton, Omkar Myatra, Andrew Thornton, Thomas Bosman, Geof Cunningham, Basil Duval, Stuart Henderson, Rory Scannell, the MAST Upgrade team

Abstract: We present the first analysis of the atomic and molecular processes at play during detachment in the MAST-U Super-X divertor using divertor spectroscopy data. Our analysis indicates detachment in the MAST-U Super-X divertor can be separated into four sequential phases: First, the ionisation region detaches from the target at detachment onset leaving a region of increased molecular densities downst… ▽ More We present the first analysis of the atomic and molecular processes at play during detachment in the MAST-U Super-X divertor using divertor spectroscopy data. Our analysis indicates detachment in the MAST-U Super-X divertor can be separated into four sequential phases: First, the ionisation region detaches from the target at detachment onset leaving a region of increased molecular densities downstream. The plasma interacts with these molecules, resulting in molecular ions ($D_2^+$ and/or $D_2^- \rightarrow D + D^-$) that further react with the plasma leading to Molecular Activated Recombination and Dissociation (MAR and MAD), which results in excited atoms and significant Balmer line emission. Second, the MAR region detaches from the target leaving a sub-eV temperature region downstream. Third, an onset of strong emission from electron-ion recombination (EIR) ensues. Finally, the electron density decays near the target, resulting in a density front moving upstream. The analysis in this paper indicates that plasma-molecule interactions have a larger impact than previously reported and play a critical role in the intensity and interpretation of hydrogen atomic line emission characteristics on MAST-U. Furthermore, we find that the Fulcher band emission profile in the divertor can be used as a proxy for the ionisation region and may also be employed as a plasma temperature diagnostic for improving the separation of hydrogenic emission arising from electron-impact excitation and that from plasma-molecular interactions. We provide evidences for the presence of low electron temperatures ($<0.5$ eV) during detachment phases III-IV based on quantitative spectroscopy analysis, a Boltzmann relation of the high-n Balmer line transitions together with an analysis of the brightness of high-n Balmer lines. △ Less

Submitted 18 October, 2022; v1 submitted 5 April, 2022; originally announced April 2022.

X-point and divertor filament dynamics from Gas Puff Imaging on TCV

Authors: Curdin Wüthrich, Christian Theiler, Nicola Offeddu, Davide Galassi, Diego Sales de Oliveira, Basil Duval, Olivier Février, Theodore Golfinopoulos, Woonghee Han, Earl Marmar, Jim Terry, Cedric Tsui

Abstract: A new Gas Puff Imaging (GPI) diagnostic has been installed on the TCV tokamak, providing two-dimensional insights into Scrape-Off-Layer (SOL) turbulence dynamics above, at and below the magnetic X-point. A detailed study in L-mode, attached, lower single-null discharges shows that statistical properties have little poloidal variations, while vast differences are present in the 2D behaviour of inte… ▽ More A new Gas Puff Imaging (GPI) diagnostic has been installed on the TCV tokamak, providing two-dimensional insights into Scrape-Off-Layer (SOL) turbulence dynamics above, at and below the magnetic X-point. A detailed study in L-mode, attached, lower single-null discharges shows that statistical properties have little poloidal variations, while vast differences are present in the 2D behaviour of intermittent filaments. Strongly elongated filaments, just above the X-point and in the divertor far-SOL, show a good consistency in shape and dynamics with field-line tracing from filaments at the outboard midplane, highlighting their connection. In the near-SOL of the outer divertor leg, shortlived, high frequency and more circular (diameter $\sim$15 sound Larmour radii) filaments are observed. These divertor-localised filaments appear born radially at the position of maximum density and display a radially outward motion with velocity $\approx$400\,m/s that is comparable to radial velocities of upstream-connected filaments. Conversely, in these discharges ($B\times\nabla B$ pointing away from the divertor), these divertor filaments' poloidal velocities differ strongly from those of upstream-connected filaments. The importance of divertor-localised filaments upon radial transport and profile broadening is explored using filament statistics and in-situ kinetic profile measurements along the divertor leg. This provides evidence that these filaments contribute significantly to electron density profile broadening in the divertor. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: This is the version of the article before peer review or editing, as submitted by an author to IOPScience Nuclear Fusion. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. 15 pages, 15 figures

Surrogate-Based Black-Box Optimization Method for Costly Molecular Properties

Authors: Jules Leguy, Thomas Cauchy, Beatrice Duval, Benoit Da Mota

Abstract: AI-assisted molecular optimization is a very active research field as it is expected to provide the next-generation drugs and molecular materials. An important difficulty is that the properties to be optimized rely on costly evaluations. Machine learning methods are investigated with success to predict these properties, but show generalization issues on less known areas of the chemical space. We p… ▽ More AI-assisted molecular optimization is a very active research field as it is expected to provide the next-generation drugs and molecular materials. An important difficulty is that the properties to be optimized rely on costly evaluations. Machine learning methods are investigated with success to predict these properties, but show generalization issues on less known areas of the chemical space. We propose here a surrogate-based black box optimization method, to tackle jointly the optimization and machine learning problems. It consists in optimizing the expected improvement of the surrogate of a molecular property using an evolutionary algorithm. The surrogate is defined as a Gaussian Process Regression (GPR) model, learned on a relevant area of the search space with respect to the property to be optimized. We show that our approach can successfully optimize a costly property of interest much faster than a purely metaheuristic approach. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Submitted to ICTAI 2021

The role of plasma-molecule interactions on power and particle balance during detachment on the TCV tokamak

Authors: Kevin Verhaegh, Bruce Lipschultz, James Harrison, Basil Duval, Alex Fil, Mirko Wensing, Chris Bowman, Daljeet Gahle, Andrei Kukushkin, David Moulton, Artur Perek, Andrey Pshenov, Fabio Federici, Olivier Février, Omkar Myatra, Andreas Smolders, Christian Theiler

Abstract: This paper shows experimental results from the TCV tokamak that indicate plasma-molecule interactions involving $D_2^+$ and possibly $D^-$ play an important role as sinks of energy (through hydrogenic radiation as well as dissociation) and particles during divertor detachment if low target temperatures ($< 3$ eV) are achieved. Both molecular activated recombination (MAR) and ion source reduction d… ▽ More This paper shows experimental results from the TCV tokamak that indicate plasma-molecule interactions involving $D_2^+$ and possibly $D^-$ play an important role as sinks of energy (through hydrogenic radiation as well as dissociation) and particles during divertor detachment if low target temperatures ($< 3$ eV) are achieved. Both molecular activated recombination (MAR) and ion source reduction due to a power limitation effect are shown to be important in reducing the ion target flux during a density ramp. In contrast, the electron-ion recombination (EIR) ion sink is too small to play an important role in reducing the ion target flux. MAR or EIR do not occur during $N_2$ seeding induced detachment as the target temperatures are not sufficiently low. The impact of $D_2^+$ is shown to be underestimated in present (vibrationally unresolved) SOLPS-ITER simulations, which could result from an underestimated $D_2 + D^+ \rightarrow D_2^+ + D$ rate. The converged SOLPS-ITER simulations are post-processed with alternative reaction rates, resulting in considerable contributions of $D_2^+$ to particle and power losses as well as dissociation below the $D_2$ dissociation area. Those findings are in quantitative agreement with the experimental results. △ Less

Submitted 30 July, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Latest results on quiescent and post-disruption runaway electron mitigation experiments at Frascati Tokamak Upgrade

Authors: D. Carnevale, P. Buratti, M. Baruzzo, W. Bin, F. Bombarda, L. Boncagni, C. Paz-Soldan, L. Calacci, M. Cappelli, C. Castaldo, S. Ceccuzzi, C. Centioli, C. Cianfarani, S. Coda, F. Cordella, O. D Arcangelo, J. Decker, B. Duval, B. Esposito, L. Gabellieri, S. Galeani, S. Garavaglia, C. Galperti, G. Ghillardi, G. Granucci , et al. (16 additional authors not shown)

Abstract: Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissi… ▽ More Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissipative effects in terms of critical electrical field increase due to fan-like instabilities. Analysis of large fan-like instabilities on post-disruption RE beams, that seem to be correlated with low electrical field and background density drops, reveal their strong RE energy suppression capability suggesting a new strategy for RE energy suppression controlling large fan instabilities. We demonstrate how such density drops can be induced using modulated ECRH power on post-disruption beams. △ Less

Submitted 25 May, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: DRAFT that will be submitted to NF (the conference version is on FEC/IAEA conference 2020)

Tomographic reconstruction of the runaway distribution function in TCV using multispectral synchrotron images

Authors: T. A. Wijkamp, A. Perek, J. Decker, B. Duval, M. Hoppe, G. Papp, U. A. Sheikh, I. G. J. Classen, R. J. E. Jaspers, the TCV team, the EUROfusion MST1 team

Abstract: Synchrotron radiation observed in a quiescent TCV runaway discharge is studied using filtered camera images targeting three distinct wavelength intervals. Through the tomographic SART procedure the high momentum, high pitch angle part of the spatial and momentum distribution of these relativistic particles is reconstructed. Experimental estimates of the distribution are important for verification… ▽ More Synchrotron radiation observed in a quiescent TCV runaway discharge is studied using filtered camera images targeting three distinct wavelength intervals. Through the tomographic SART procedure the high momentum, high pitch angle part of the spatial and momentum distribution of these relativistic particles is reconstructed. Experimental estimates of the distribution are important for verification and refinement of formation-, decay- and transport-models underlying runaway avoidance and mitigation strategy design. Using a test distribution it is demonstrated that the inversion procedure provides estimates accurate to within a few tens of percent in the region of phase-space contributing most to the synchrotron image. We find that combining images filtered around different parts of the emission spectrum widens the probed part of momentum-space and reduces reconstruction errors. Next, the SART algorithm is used to obtain information on the spatiotemporal runaway momentum distribution in a selected TCV discharge. The momentum distribution is found to relax towards an avalanche-like exponentially decaying profile. Anomalously high pitch angles and a radial profile increasing towards the edge are found for the most strongly emitting particles in the distribution. Pitch angle scattering by toroidal magnetic field ripple is consistent with this picture. An alternative explanation is the presence of high frequency instabilities in combination with the formation of a runaway shell at the edge of the plasma. △ Less

Submitted 5 February, 2021; v1 submitted 30 October, 2020; originally announced November 2020.

Comments: 18 pages, 13 figures, submitted to Nuclear Fusion

A novel hydrogenic spectroscopic technique for inferring the role of plasma-molecule interaction on power and particle balance during detached conditions

Authors: K Verhaegh, B Lipschultz, C Bowman, B P Duval, U Fantz, A Fil, J R Harrison, D Moulton, O Myatra, D Wünderlich, F Federici, D S Gahle, A Perek, M Wensing, the TCV team, the EuroFusion MST1 team

Abstract: Detachment, an important mechanism for reducing target heat deposition, is achieved through reductions in power, particle and momentum; which are induced through plasma-atom and plasma-molecule interactions. Experimental research in how those reactions precisely contribute to detachment is limited. In this work, we investigate a new spectroscopic technique to utilise Hydrogen Balmer line measureme… ▽ More Detachment, an important mechanism for reducing target heat deposition, is achieved through reductions in power, particle and momentum; which are induced through plasma-atom and plasma-molecule interactions. Experimental research in how those reactions precisely contribute to detachment is limited. In this work, we investigate a new spectroscopic technique to utilise Hydrogen Balmer line measurements to 1) disentangle the Balmer line emission from the various plasma-atom and plasma-molecule interactions; and 2) quantify their contributions to ionisation, recombination and radiative power losses. During detachment, the observed $Hα$ emission often strongly increases, which could be an indicator for plasma-molecule interactions involving $H_2^+$ and/or $H^-$. Our analysis technique quantifies the $Hα$ emission due to plasma-molecule interactions and uses this to 1) quantify the Balmer line emission contribution due to $H_2^+$ and/or $H^-$; 2) subsequently estimate its resulting particle sinks/sources and radiative power losses. Its performance is verified using synthetic diagnostic techniques of both detached TCV and MAST-U SOLPS-ITER simulations. Experimental results of this technique on TCV data show a bifurcation occurs between the measured total $Hα$ and the atomic estimate of $Hα$ emission, indicative of the presence of additional $Hα$ due to plasma-molecule interactions with $H_2^+$ (and/or $H^-$). An example analysis shows that the hydrogenic line series, even $Lyα$ as well as the medium-n Balmer lines can be significantly influenced by plasma-molecule interactions by tens of percent during which significant Molecular Activated Recombination (MAR) is expected. △ Less

Submitted 17 February, 2021; v1 submitted 1 August, 2020; originally announced August 2020.

Journal ref: Plasma Phys. Control. Fusion 63 035018 (2021)

A study of the influence of plasma-molecule interactions on particle balance during detachment

Authors: K Verhaegh, B Lipschultz, J R Harrison, B P Duval, C Bowman, A Fil, D S Gahle, D Moulton, O Myatra, A Perek, C Theiler, M Wensing

Abstract: In this work we provide experimental insights into the impact of plasma-molecule interactions on the target ion flux decrease during divertor detachment achieved through a core density ramp in the TCV tokamak. Our improved analysis of the hydrogen Balmer series shows that plasma-molecule processes are strongly contributing to the Balmer series intensities and substantially alter the divertor detac… ▽ More In this work we provide experimental insights into the impact of plasma-molecule interactions on the target ion flux decrease during divertor detachment achieved through a core density ramp in the TCV tokamak. Our improved analysis of the hydrogen Balmer series shows that plasma-molecule processes are strongly contributing to the Balmer series intensities and substantially alter the divertor detachment particle balance. We find that Molecular Activated Recombination (MAR) ion sinks from $H_2^+$ and/or $H^-$ are a factor $\sim$ 5 larger than Electron-Ion Recombination (EIR) and are a significant contributor to the observed reduction in the outer divertor ion target flux. Molecular Activated Ionisation (MAI) may also be significant during detachment. Plasma-molecule interactions enhance the Balmer line series emission strongly near the target as detachment proceeds. This indicates enhancements of the Lyman series, potentially affecting power balance in the divertor. As those enhancements vary spatially in the divertor and are different for different transitions, they are expected to result in a separation of the $Lyβ$ and $Lyα$ emission regions. This may have implications for the treatment and diagnosis of divertor opacity. The demonstrated enhancement of the Balmer series through plasma-molecule processes potentially poses a challenge to using the Balmer series for understanding and diagnosing detachment based only on atom-plasma processes. △ Less

Submitted 17 February, 2021; v1 submitted 1 August, 2020; originally announced August 2020.

Journal ref: Nuclear Materials and Energy, 100922, 2021

Runaway electron synchrotron radiation in a vertically translated plasma

Authors: M. Hoppe, G. Papp, T. Wikjamp, A. Perek, J. Decker, B. Duval, O. Embreus, T. Fülöp, U. A. Sheikh

Abstract: Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron di… ▽ More Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron diagnostic Soft. This experimental validation lends confidence to the theory underlying the synthetic diagnostics which are used for benchmarking theoretical models of and probing runaway dynamics. We present a comparison of synchrotron measurements in TCV with predictions of kinetic theory for runaway dynamics in uniform magnetic fields. We find that to explain the detected synchrotron emission, significant non-collisional pitch angle scattering as well as radial transport of REs would be needed. Such effects could be caused by the presence of magnetic perturbations, which should be further investigated in future TCV experiments. △ Less

Submitted 18 August, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures. Accepted for publication in Nuclear Fusion

Journal ref: Nucl. Fusion 60 094002 (2020)

Population-based Gradient Descent Weight Learning for Graph Coloring Problems

Authors: Olivier Goudet, Béatrice Duval, Jin-Kao Hao

Abstract: Graph coloring involves assigning colors to the vertices of a graph such that two vertices linked by an edge receive different colors. Graph coloring problems are general models that are very useful to formulate many relevant applications and, however, are computationally difficult. In this work, a general population-based weight learning framework for solving graph coloring problems is presented.… ▽ More Graph coloring involves assigning colors to the vertices of a graph such that two vertices linked by an edge receive different colors. Graph coloring problems are general models that are very useful to formulate many relevant applications and, however, are computationally difficult. In this work, a general population-based weight learning framework for solving graph coloring problems is presented. Unlike existing methods for graph coloring that are specific to the considered problem, the presented work targets a generic objective by introducing a unified method that can be applied to different graph coloring problems. This work distinguishes itself by its solving approach that formulates the search of a solution as a continuous weight tensor optimization problem and takes advantage of a gradient descent method computed in parallel on graphics processing units. The proposed approach is also characterized by its general global loss function that can easily be adapted to different graph coloring problems. The usefulness of the proposed approach is demonstrated by applying it to solve two typical graph coloring problems and performing large computational studies on popular benchmarks. Improved best-known results (new upper bounds) are reported for several large graphs. △ Less

Submitted 26 October, 2020; v1 submitted 5 September, 2019; originally announced September 2019.

Novel inferences of ionisation & recombination for particle/power balance during detached discharges using deuterium Balmer line spectroscopy

Authors: K. Verhaegh, B. Lipschultz, B. P. Duval, A. Fil, M. Wensing, C. Bowman, D. S. Gahle

Abstract: The physics of divertor detachment is determined by divertor power, particle and momentum balance. This work provides a novel analysis technique of the Balmer line series to obtain a full particle/power balance measurement of the divertor. This supplies new information to understand what controls the divertor target ion flux during detachment. Atomic deuterium excitation emission is separated fr… ▽ More The physics of divertor detachment is determined by divertor power, particle and momentum balance. This work provides a novel analysis technique of the Balmer line series to obtain a full particle/power balance measurement of the divertor. This supplies new information to understand what controls the divertor target ion flux during detachment. Atomic deuterium excitation emission is separated from recombination quantitatively using Balmer series line ratios. This enables analysing those two components individually, providing ionisation/recombination source/sinks and hydrogenic power loss measurements. Probabilistic Monte Carlo techniques were employed to obtain full error propagation - eventually resulting in probability density functions for each output variable. Both local and overall particle and power balance in the divertor are then obtained. These techniques and their assumptions have been verified by comparing the analysed synthetic diagnostic 'measurements' obtained from SOLPS simulation results for the same discharge. Power/particle balance measurements have been obtained during attached and detached conditions on the TCV tokamak. △ Less

Submitted 27 November, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: The analysis results of this paper were formerly in arXiv:1810.04969

Journal ref: 2019 Plasma Phys. Control. Fusion 61 125018

140 Gbaud On-Off Keying Links in C-Band for Short-Reach Optical Interconnects

Authors: Oskars Ozolins, José Manuel Estaran, Aleksejs Udalcovs, Filipe Jorge, Haik Mardoyan, Agnieszka Konczykowska, Muriel Riet, Bernadette Duval, Virginie Nodjiadjim, Jean-Yves Dupuy, Xiaodan Pang, Urban Westergren, Jiajia Chen, Sergei Popov, Sébastien Bigo

Abstract: We demonstrate 140 Gbaud intensity modulated direct detection dispersion-uncompensated links with Mach Zehnder modulator and distributed feedback travelling-wave electro-absorption modulator over 5500 and 960 meters of standard single mode fibre, respectively, enabled by compact packaged ultra-high speed InP-based 2:1-Selector. We demonstrate 140 Gbaud intensity modulated direct detection dispersion-uncompensated links with Mach Zehnder modulator and distributed feedback travelling-wave electro-absorption modulator over 5500 and 960 meters of standard single mode fibre, respectively, enabled by compact packaged ultra-high speed InP-based 2:1-Selector. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Comments: 3 pages, 44th European Conference on Optical Communication (ECOC 2018), Rome, Italy, 2018

An improved understanding of the roles of atomic processes and power balance in divertor target ion current loss during detachment

Authors: Kevin Verhaegh, Bruce Lipschultz, Basil Duval, Olivier Février, Alexandre Fil, Christian Theiler, Mirko Wensing, Christopher Bowman, Daljeet Gahle, James Harrison, Benoit Labit, Claudio Marini, Roberto Maurizio, Hugo de Oliveira, Holger Reimerdes, Umar Sheikh, Cedric Tsui, Nicola Vianello, Wouter Vijvers

Abstract: The process of divertor detachment, whereby heat and particle fluxes to divertor surfaces are strongly diminished, is required to reduce heat loading and erosion in a magnetic fusion reactor to acceptable levels. In this paper the physics leading to the decrease of the total divertor ion current (It), or 'roll-over', is experimentally explored on the TCV tokamak through characterization of the loc… ▽ More The process of divertor detachment, whereby heat and particle fluxes to divertor surfaces are strongly diminished, is required to reduce heat loading and erosion in a magnetic fusion reactor to acceptable levels. In this paper the physics leading to the decrease of the total divertor ion current (It), or 'roll-over', is experimentally explored on the TCV tokamak through characterization of the location, magnitude and role of the various divertor ion sinks and sources including a complete analysis of particle and power balance. These first measurements of the profiles of divertor ionisation and hydrogenic radiation along the divertor leg are enabled through novel spectroscopic techniques. Over a range in TCV plasma conditions (plasma current and electron density, with/without impurity-seeding) the $I_t$ roll-over is ascribed to a drop in the divertor ion source; recombination remains small or negligible farther into the detachment process. The ion source reduction is driven by both a reduction in the power available for ionization, Precl, and concurrent increase in the energy required per ionisation, $E_{ion}$: often described as 'power starvation' (or 'power limitation'). The detachment threshold is found experimentally (in agreement with analytic model predictions) to be $\sim P_{recl}/I_t {E_{ion}} \sim 2$, corresponding to a target electron temperature, $T_t \sim E_{ion}/γ$ where $γ$ is the sheath transmission coefficient. The target pressure reduction, required to reduce the target ion current, is driven both by volumetric momentum loss as well as upstream pressure loss. The measured evolution through detachment of the divertor profile of various ion sources/sinks as well as power losses are quantitatively reproduced through full 2D SOLPS modelling through the detachment process as the core density is varied. △ Less

Submitted 20 August, 2019; v1 submitted 11 October, 2018; originally announced October 2018.

Comments: The analysis part of this paper is now presented arxiv:1903.08157

Journal ref: 2019 Nucl. Fusion 59 126038

When data mining meets optimization: A case study on the quadratic assignment problem

Authors: Yangming Zhou, Jin-Kao Hao, Béatrice Duval

Abstract: This paper presents a hybrid approach called frequent pattern based search that combines data mining and optimization. The proposed method uses a data mining procedure to mine frequent patterns from a set of high-quality solutions collected from previous search, and the mined frequent patterns are then employed to build starting solutions that are improved by an optimization procedure. After prese… ▽ More This paper presents a hybrid approach called frequent pattern based search that combines data mining and optimization. The proposed method uses a data mining procedure to mine frequent patterns from a set of high-quality solutions collected from previous search, and the mined frequent patterns are then employed to build starting solutions that are improved by an optimization procedure. After presenting the general approach and its composing ingredients, we illustrate its application to solve the well-known and challenging quadratic assignment problem. Computational results on the 21 hardest benchmark instances show that the proposed approach competes favorably with state-of-the-art algorithms both in terms of solution quality and computing time. △ Less

Submitted 7 October, 2017; v1 submitted 17 August, 2017; originally announced August 2017.

Experimental observations and modelling of intrinsic rotation reversals in tokamaks

Authors: Y. Camenen, C. Angioni, A. Bortolon, B. P. Duval, E. Fable, W. A. Hornsby, R. M. Mcdermott, D. H. Na, Y-S. Na, A. G. Peeters, J. E. Rice

Abstract: The progress made in understanding spontaneous toroidal rotation reversals in tokamaks is reviewed and current ideas to solve this ten-year-old puzzle are explored. The paper includes a summarial synthesis of the experimental observations in AUG, C-Mod, KSTAR, MAST and TCV tokamaks, reasons why turbulent momentum transport is thought to be responsible for the reversals, a review of the theory of t… ▽ More The progress made in understanding spontaneous toroidal rotation reversals in tokamaks is reviewed and current ideas to solve this ten-year-old puzzle are explored. The paper includes a summarial synthesis of the experimental observations in AUG, C-Mod, KSTAR, MAST and TCV tokamaks, reasons why turbulent momentum transport is thought to be responsible for the reversals, a review of the theory of turbulent momentum transport and suggestions for future investigations. △ Less

Submitted 27 January, 2017; originally announced January 2017.

Journal ref: Plasma Physics and Controlled Fusion, IOP Publishing, 2017, 59, pp.34001 - 34001

Spectroscopic investigations of divertor detachment in TCV

Authors: K. Verhaegh, B. Lipschultz, B. P. Duval, J. R. Harrison, H. Reimerdes, C. Theiler, B. Labit, R. Maurizio, C. Marini, F. Nespoli, U. Sheikh, C. K. Tsui, N. Vianello, W. A. J. Vijvers, TCV team, MST1 team

Abstract: The aim of this work is to provide an understanding of detachment at TCV with emphasis on analysis of the Balmer line emission. A new Divertor Spectroscopy System has been developed for this purpose. Further development of Balmer line analysis techniques has allowed detailed information to be extracted from the three-body recombination contribution to the n=7 Balmer line intensity. During densit… ▽ More The aim of this work is to provide an understanding of detachment at TCV with emphasis on analysis of the Balmer line emission. A new Divertor Spectroscopy System has been developed for this purpose. Further development of Balmer line analysis techniques has allowed detailed information to be extracted from the three-body recombination contribution to the n=7 Balmer line intensity. During density ramps, the plasma at the target detaches as inferred from a drop in ion current to the target. At the same time the Balmer $6\rightarrow2$ and $7\rightarrow2$ line emission near the target is dominated by recombination. As the core density increases further, the density and recombination rate are rising all along the outer leg to the x-point while remaining highest at the target. Even at the highest core densities accessed (Greenwald fraction 0.7) the peaks in recombination and density may have moved not more than a few cm poloidally away from the target which is different to other, higher density tokamaks, where both the peak in recombination and density continue to move towards the x-point as the core density is increased. The inferred magnitude of recombination is small compared to the target ion current at the time detachment (particle flux drop) starts at the target. However, recombination may be having more localized effects (to a flux tube) which we cannot discern at this time. Later, at the highest densities achieved, the total recombination does reach levels similar to the particle flux. △ Less

Submitted 16 January, 2017; v1 submitted 15 July, 2016; originally announced July 2016.

Comments: Article accepted for publication in Journal of Nuclear Materials and Energy

Journal ref: Nuclear Materials and Energy Volume 12, August 2017, Pages 1112-1117

Reinforcement learning based local search for grouping problems: A case study on graph coloring

Authors: Yangming Zhou, Jin-Kao Hao, Béatrice Duval

Abstract: Grouping problems aim to partition a set of items into multiple mutually disjoint subsets according to some specific criterion and constraints. Grouping problems cover a large class of important combinatorial optimization problems that are generally computationally difficult. In this paper, we propose a general solution approach for grouping problems, i.e., reinforcement learning based local searc… ▽ More Grouping problems aim to partition a set of items into multiple mutually disjoint subsets according to some specific criterion and constraints. Grouping problems cover a large class of important combinatorial optimization problems that are generally computationally difficult. In this paper, we propose a general solution approach for grouping problems, i.e., reinforcement learning based local search (RLS), which combines reinforcement learning techniques with descent-based local search. The viability of the proposed approach is verified on a well-known representative grouping problem (graph coloring) where a very simple descent-based coloring algorithm is applied. Experimental studies on popular DIMACS and COLOR02 benchmark graphs indicate that RLS achieves competitive performances compared to a number of well-known coloring algorithms. △ Less

Submitted 1 April, 2016; originally announced April 2016.

An optimal real-time controller for vertical plasma stabilization

Authors: N. Cruz, J-M. Moret, S. Coda, B. P. Duval, H. B. Le, A. P. Rodrigues, C. A. F. Varandas, C. M. B. A. Correia, B. S. Goncalves

Abstract: Modern Tokamaks have evolved from the initial axisymmetric circular plasma shape to an elongated axisymmetric plasma shape that improves the energy confinement time and the triple product, which is a generally used figure of merit for the conditions needed for fusion reactor performance. However, the elongated plasma cross section introduces a vertical instability that demands a real-time feedback… ▽ More Modern Tokamaks have evolved from the initial axisymmetric circular plasma shape to an elongated axisymmetric plasma shape that improves the energy confinement time and the triple product, which is a generally used figure of merit for the conditions needed for fusion reactor performance. However, the elongated plasma cross section introduces a vertical instability that demands a real-time feedback control loop to stabilize the plasma vertical position and velocity. At the Tokamak à Configuration Variable (TCV) in-vessel poloidal field coils driven by fast switching power supplies are used to stabilize highly elongated plasmas. TCV plasma experiments have used a PID algorithm based controller to correct the plasma vertical position. In late 2013 experiments a new optimal real-time controller was tested improving the stability of the plasma. This contribution describes the new optimal real-time controller developed. The choice of the model that describes the plasma response to the actuators is discussed. The high order model that is initially implemented demands the application of a mathematical order reduction and the validation of the new reduced model. The lower order model is used to derive the time optimal control law. A new method for the construction of the switching curves of a bang-bang controller is presented that is based on the state-space trajectories that optimize the time to target of the system. A closed loop controller simulation tool was developed to test different possible algorithms and the results were used to improve the controller parameters. The final control algorithm and its implementation are described and preliminary experimental results are discussed. △ Less

Submitted 18 June, 2014; v1 submitted 13 June, 2014; originally announced June 2014.

Comments: Preprint submitted to IEEE Transactions on Nuclear Science