Nuclear Fusion最新文献

{"title":"Electron-scale turbulence characteristics with varying electron temperature gradient in LHD","authors":"Tatsuhiro Nasu, T. Tokuzawa, Motoki Nakata, K. Ida, S. Inagaki, Masaki Nishiura, Yasuo Yoshimura, Ryoma Yanai, Kenji Tanaka, M. Yoshinuma, Tatsuya Kobayashi, Akira Ejiri, K. Y. Watanabe, Ichihiro Yamada","doi":"10.1088/1741-4326/ad5d7c","DOIUrl":"https://doi.org/10.1088/1741-4326/ad5d7c","url":null,"abstract":"\u0000 The electron-scale turbulence, whose wavelength is about the electron Larmor radius, is thought to have the potential to cause the stiffness of the electron temperature gradient and degrade the confinement of future burning plasma in which the electron heating by alpha particles is dominant. The dependence of electron-scale turbulence and electron heat flux on the electron temperature inverse gradient length Rax/LTe were investigated. The electron temperature gradient is successfully varied in the range of −3 < Rax/LTe < 12 by controlling the injection power of the on/off-axis electron cyclotron heating. The results show a significant increase in electron-scale turbulence with increasing Rax/LTe, especially in conditions where the Electron Temperature Gradient (ETG) instability is linearly unstable, suggesting the presence of ETG turbulence at high Rax/LTe. The electron heat flux also increases steeply with increasing Rax/LTe. In addition, electron-scale turbulence is observed even at Rax/LTe ∼ 0, which is stable in linear GKV calculations. Finding the cause of this phenomenon is an interesting task for the future.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"43 5","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141713091","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Proof-of-principle of parametric stellarator neutronics modeling using Serpent2","authors":"Tommi Tapani Lyytinen, Antti Snicker, Joonas Virtanen, Iole Palermo, Javier Alguacil, Timo Jos Bogaarts, Felix Warmer","doi":"10.1088/1741-4326/ad4f9f","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4f9f","url":null,"abstract":"\u0000 This contribution presents neutron transport studies for the 5-period HELIAS stellarator using the Serpent2 code. These studies utilize a parametric geometry model, enabling scans in neutronics modeling by varying the thickness of the reactor layers. For example, the tritium breeding ratio (TBR) can be determined by exploring various blanket material options and thicknesses to identify the threshold configuration that meets the TBR design criterion of 1.15. We found out that with the helium-cooled pebble ped (HCPB) candidate option, the TBR criterion is met with a breeding zone thickness of 26 cm, while with the dual-coolant lithium lead (DCLL) the threshold is exceeded at a thickness of 46 cm. Furthermore, the geometry includes non-planar field coils, allowing to study the fast neutron flux in these superconducting coils with a technological limit of 1e9 1/cm2s. It is shown that the neutron fast flux is not constant at the coils, necessitating a neutron transport simulation to determine the distribution of the fast-flux at the coils. We show that the peak flux can be more than a factor of 2 higher than the average flux, and that the peak flux location rotates helically.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"21 8","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141105507","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: Actuator design for particle and divertor condition control","authors":"Jae-Sun Park, Jeremy D Lore, Matthew L Reinke, Adam Q. Kuang, Sebastian De Pascuale, Alexander J Creely","doi":"10.1088/1741-4326/ad4f9d","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4f9d","url":null,"abstract":"\u0000 This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1·10-3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"17 6","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141107139","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Tokamak divertor plasma emulation with machine learning","authors":"G. Holt, A. Keats, Stanislas Pamela, Mike Kryjak, A. Agnello, N. Amorisco, Benjamin Dudson, M. Smyrnakis","doi":"10.1088/1741-4326/ad4f9e","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4f9e","url":null,"abstract":"\u0000 Future tokamak devices that aim to create conditions relevant to power plant operations must consider strategies for mitigating damage to plasma facing components in the divertor. One of the goals of MAST-U tokamak operations is to inform these considerations by researching advanced divertor configurations that aid stable plasma detachment. Machine design, scenario planning and detachment control would all greatly benefit from tools that enable rapid calculation of scenario-relevant quantities given some input parameters. This paper presents a method for generating large, simulated scrape-off layer data sets, which was applied to generate a data set of steady-state Hermes-3 simulations of the MAST-U tokamak. A machine learning model was constructed using a Bayesian approach to hyperparameter optimisation to predict diagnosable output quantities given control-relevant input features. The resulting best-performing model, which is based on a feedfoward neural network, achieves high accuracy when predicting electron temperature at the divertor target and carbon impurity radiation front position and runs in around 1 ms in inference mode. Techniques for interpreting the predictions made by the model were applied, and a high-resolution parameter scan of upstream conditions was performed to demonstrate the utility of rapidly generating accurate predictions using the emulator. This work represents a step forward in the design of machine learning-driven emulators of tokamak exhaust simulation codes in operational modes relevant to divertor detachment control and plasma scenario design.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"36 46","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141104070","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Super-X and Conventional divertor configurations in MAST-U ohmic L-mode; a comparison facilitated by interpretative modelling","authors":"David Moulton, J. R. Harrison, L. Xiang, Peter Ryan, Andrew Kirk, Kevin Verhaegh, T. Wijkamp, F. Federici, James G Clark, B. Lipschultz","doi":"10.1088/1741-4326/ad4f9c","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4f9c","url":null,"abstract":"\u0000 Measurements are presented, alongside corresponding interpretative SOLPS-ITER simulations, of the first MAST-U experiments comparing ohmically heated L-mode fuelling scans in Conventional divertor (CD) and Super-X divertor (SXD) configurations. In experiment, at comparable outer mid-plane separatrix electron density, $n_{e,rm{sep,OMP}}$, the maximum lower outer target heat load was found to be a factor 16$pm7$ lower in SXD compared to CD. In simulation, a factor 26.8 reduction was found (slightly higher than the experimental range), suggesting an additional reduction in SXD compared to the factor 9.3 expected from geometric considerations alone. According to the simulations, this additional reduction in the SXD is due to a net radial transport of the energy remaining downstream of the $T_e=5$ eV location. This energy is carried out of the critical (highest heat load) flux tube by deuterium atoms, demonstrating the importance of a longer legged divertor which provides space for this to occur. Importantly, in both simulation and experiment, the SXD has minimal impact on the upstream $n_e$ and $T_e$ profiles. Spectral inferences of detachment front movement in SXD compare well between simulation and experiment. In regions of high magnetic field gradient, the parallel movement of the front towards the X-point becomes less sensitive to increasing $n_{e,rm{sep,OMP}}$, in qualitative agreement with simplified models and previous predictive simulations. Additional aspects, regarding the target ion flux rollover, upstream separatrix temperature and drift effects, are also presented and discussed.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"47 3","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141103466","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Experimental investigation on PDI bifurcation of lower hybrid wave during electron density ramp-up in EAST","authors":"Guanghou Yan, Bojiang Ding, Chenbin Wu, Miaohui Li, Seung Gyou Baek, Francesco Napoli, Taotao Zhou, K. Geng, N. Yan, Mao Wang, Xinjun Zhang, Handong Xu, J. H. Yang, W. Ma","doi":"10.1088/1741-4326/ad4fa0","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4fa0","url":null,"abstract":"\u0000 Effect of parametric decay instability (PDI) on current drive efficiency of 4.6 GHz lower hybrid wave in EAST is investigated experimentally, showing the PDI channel bifurcation of lower hybrid (LH) wave for the first time in EAST. Firstly, experiments with 3 platforms of LH power were performed, achieving the LH power required for the PDI occurrence. Secondly, PDI bifurcation experiments were further carried out by ramping up plasma electron density. Loop voltage increases with density increase, implying the decrease of current driven by LH wave. PDI bifurcation during electron density ramp-up was studied by analyzing the parallel refractive index (n_∥) and the frequency spectrum broadening, which is measured with a radio frequency (RF) magnetic probe array recently installed close to the LH antenna. It is observed for the first time that both of them firstly increase with density, then there is not much variation and a clear sideband in the frequency spectrum is also observed when the density is up to 4×10^19 m^(-3), suggesting the change of the PDI channel. Calculation of the mode growth rate driven by PDI shows that when the edge electron density is up to 1.9×10^18 m^(-3), the growth rate of ion cyclotron quasi-mode (ICQM) will exceed that of ion sound quasi-mode (ISQM), quantitively explaining that with a density increase, PDI channel partly transits from ISQM to ICQM channel. Studies provide a possible way to reduce the power deposition in the edge region and improve LHCD capability by means of mitigating PDI behavior.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"39 13","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141103594","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Research on the efficient process-oriented structural optimization method of the large-scale vacuum cryostat for fusion reactor","authors":"Qingzhou Yu, Hao Xu, Zhaoxi Chen, Qingxi Yang","doi":"10.1088/1741-4326/ad4ef2","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4ef2","url":null,"abstract":"\u0000 An efficient optimization design for large and complex components of fusion reactor is crucial to address the engineering design requirements and further promote technical standardization. Based on research status, current engineering designs for fusion reactors have some deficiencies such as time and energy wastes, inefficiency and the difficulties in covering the typical “multi-variable, multi-objective” design requirements. They are pressing and common problems that urgently need to be overcome. To deal with the aforementioned technical challenges, it is vitally important to design an efficient, precise, and normalized approach tailored for the development of future fusion reactors. Therefore, this paper proposes a process-oriented optimization design method, which involves Coupled external parameterized modeling, Experimental points design, Response surface optimization and Structural integrity validation (CERS), to improve currently inefficient design methods. And the vacuum cryostat, the largest and complex component of tokamak, is taken as an example to present the basic procedures of CERS. Firstly, the functions, basic structures, load types, analysis methods and verification criteria of the cryostat are presented in detail. And then real-time data interaction between external global parametric variables and ANSYS via coupling is established by CERS, which achieves parametric modeling of cryostat and the efficient experiment point design and optimization analysis with multi-variables and multi-objectives in an automatic way. Subsequently, this study demonstrates the significance and sensitivity of various structural parameters of cryostat from such objectives as maximum deformation, maximum equivalent stress, and total mass. And the optimal set of its structural parameters is obtained by establishing mathematical optimization model. Finally, the structural integrity is verified. The result indicates that the optimized cryostat maintains a minimum safety margin of 23% and will not suffer fatigue damage under various load events during its service. Besides, the nonlinear buckling load multiplier ∅ is 5.4, obtained by analyzing the load-displacement curve of cryostat according to the zero-curvature criterion. This shows that the designed cryostat is stable enough. The proposed method is simple, efficient and reliable which can be applied to both cryostat and other complex components of fusion reactors in the engineering design fields. It has great value of practically technical reference and can further promote the standardization of engineering design technology for future fusion reactor.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"57 8","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141110480","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Power and isotope effects in the ITER Baseline Scenario with tungsten and tungsten-equivalent radiators in DIII-D","authors":"A. C. C. Sips, Francesca Turco, Charles Mark Greenfield, Lothar Schmitz, Timothy C Luce, Tomas Odstrcil, Adam G McLean, Igor Bykov, A W Hyatt, Tom H Osborne","doi":"10.1088/1741-4326/ad4ef3","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4ef3","url":null,"abstract":"\u0000 FEC 2023 contribution: Experiments in DIII-D document the ITER Baseline Scenario (IBS) at q95 ~ 3 and PIN/PLH ~ 1-2, in both deuterium and hydrogen utilizing Kr and Xe as Tungsten-equivalent radiators. The power threshold for H-mode operation (PLH) was found to be about a factor of two higher than the scaling law. In recent IBS experiments in deuterium, intrinsic levels of metals such as Tungsten (W) or molybdenum and inconel are present that reduce the pedestal pressure by 20-25%. A complete radiative collapse of deuterium IBS plasmas occurs at W core concentrations CW = 10-5. Simulations show that for core temperatures expected for ITER, the plasmas would not have a radiative collapse at CW = 1x10-5, moreover Q = 8-10 would still be achieved for CW up to 3x10-5. In contrast to deuterium, the IBS in hydrogen is not affected by intrinsic high-Z impurities. Krypton was injected in a matrix scan of input power and impurity flow in IBS hydrogen discharges. Krypton impurity density profiles in hydrogen are similar to deuterium plasmas, but at Kr flows that are 2-3 times higher for the same input power. Krypton is transported into the core and affects the whole radius; at the highest injection rates a radiative collapse occurs at core radiation fractions of 0.3-0.35, consistent with the expected maximum W radiation fraction for ITER core plasmas. Comparing the results with previous ITPA database studies of the IBS confirms that at higher radiation fraction due to high-Z impurities, a drop in H98 of >10% is observed. On the other hand, the results using Kr as a W-equivalent radiator indicate that metal (W) devices at lower core temperatures than ITER may provide overly pessimistic performance extrapolations to ITER for deuterium-tritium operation.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"8 6","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141110170","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stability analysis of plasma waves driven by runaway electrons in tokamak hot plasmas","authors":"Carmine Castaldo, Lorenzo Della Volpe, Renato Fedele, William Bin, Paolo Buratti, Alessandro Cardinali, Francesco Napoli, Massimo Marinucci, Gerarda Apruzzese, C. Cianfarani, Edmondo Giovannozzi, O. Tudisco","doi":"10.1088/1741-4326/ad4ef6","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4ef6","url":null,"abstract":"\u0000 The local stability analysis of plasma waves driven by runaway electrons (RE) has been performed considering hot plasma Maxwellian background, with electron and ion temperatures of the order of 1 keV. It is shown that hot plasma waves, namely electron plasma waves (EPW) and ion Bernstein waves (IBW) can be driven unstable by RE at their coalescence frequency via Cherenkov resonance by RE with energy distribution peaked at about 8 MeV. A skew-normal distribution is used as a model of the RE energy distribution. The EPW and IBW couples of waves occur between any successive ion-cyclotron harmonics frequencies nfci, above the lower hybrid resonance. At their confluence, the perpendicular group velocity vanishes and significant RF emissions are expected. The frequency gap between two successive confluences is ~fci. Groups of RF line emissions, separated by almost constant frequency gap ~fci/2 are detected during various quiescent runaway plasma discharges in the FTU tokamak. The analysis of a specific discharge suggests that the frequencies of the line emissions observed and the frequencies occurring at the EPW-IBW confluences are in reasonable agreement. A possible explanation of the line emissions with ~fci/2 gap in terms of nonlinear mode coupling is proposed.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"3 12","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141111232","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"The role of divertor pumping combined with full drifts in particle exhaust and divertor plasma","authors":"Xuele Zhao, C. Sang, Yilin Wang, Chen Zhang, Dezhen Wang","doi":"10.1088/1741-4326/ad4ef7","DOIUrl":"https://doi.org/10.1088/1741-4326/ad4ef7","url":null,"abstract":"\u0000 The effect of drifts combined with pumping on particle exhaust is assessed using the SOLPS-ITER code package, considering full drifts. Both drifts and pumping speed S can affect particle exhaust. Drifts change the neutral density by influencing plasma flow and the resulting particle recycling. This leads to the accumulation of neutral particles either far away or close to the pump opening location. While the particle exhaust is enhanced as S raises. When the pump opening is positioned at the common flux region (CFR) of the outer divertor (referred to as Pump CFR/OD), particle exhaust is suppressed by drifts in forward Bt, while it is enhanced by drifts in reversed Bt, with fixed S. On the other hand, when the pump is situated in the private flux region (PFR) of the outer divertor (referred to as Pump PFR/OD), particle exhaust is enhanced by drifts in both reversed and forward Bt compared to the case without drifts. Moreover, the effective pumping in reversed Bt is stronger than in forward Bt. In the same Bt direction, Pump PFR/OD has a higher effective pumping than Pump CFR/OD. Increased S results in higher particle exhaust in all Bt direction and pump location cases. The plasma detachment is affected by drifts, S and pump opening location, respectively. For the specified Bt direction and pump opening location case, higher S suppresses plasma detachment. For an identical particle exhaust rate, stronger pumping capacity can promote plasma detachment. Therefore, Pump PFR/OD can more easily achieve outer divertor detachment than Pump CFR/OD in the same Bt direction. Overall, placing the pump at the PFR side of the outer divertor while running in reversed Bt is the best option from the divertor particle exhaust and plasma detachment point of view.","PeriodicalId":503481,"journal":{"name":"Nuclear Fusion","volume":"62 36","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141110740","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}