Soft X-ray Tomography and Sawtooth Oscillations during IBW in HT-7 Tokamak

Tomographic reconstructions of soft x-ray emissivity in the HT-7 tokamak have been obtained using 86 miniature soft-x-ray detectors. Because there are two-detector arrays which are distributed in horizontal and vertical directions on one poloidal plane, two-dimensional images of x-ray emissivity are obtained without having resort to a rotational model. The algorithm of Fourier-Bessel expansion was employed. The data analysed so far have been used principally to study MHD instabilities, and the illustrative examples of the sawtooth crash in IBW-heated plasma are described.

Soft X-ray Tomography on HT-7

The soft X-ray imaging system on HT-7 is introduced. Fourier-Bessel expansion method is used to reconstruct the soft X-ray emission distribution. Simulation results show that the inversion accuracy is associated with the angular expansion limit and the hot core displacement. In addition, tomographic results of sawtooth oscillations in a typical low hybrid current drive （LHCD） plasma are presented. It is shown that the reconstructions are reliable to analyze the perturbative mode and the evolution.

Study of sawtooth oscillations on the HT-7 tokamak using 2D tomography of soft x-ray signal

It is the first time so far as we know that two arrays of multi-channel soft x-ray detectors are used to generate twodimensional （2D） images of sawtooth oscillation on the HT-7 tokamak using the Fourier-Bessel harmonic reconstruction method, and using the singular value decomposition to analyse the data from soft x-ray cameras. By these two arrays, 2D image reconstruction of soft x-ray emissivity can be obtained without assumption of plasma rigid rotation. Tomographic reconstruction of the m=1 mode structure is obtained during the precursor oscillation of the sawtooth crash. The crescent-shaped mode structure appearing on the contour map of the soft x-ray emissivity is consistent with the quasiinterchange mode. The characteristics of the m=1/n=1 mode structure observed in the soft x-ray tomography are as follows： the magnetic surface is made up of the crescent-shaped “hot core” and the circular “cold bubble”. The structure of the magnetic surface rotates in the direction of the electron diamagnetic drift and the rotation frequency is the oscillation frequency of soft x-ray signals.

The Soft X-ray Imager(SXI)on the SMILE Mission

The Soft X-ray Imager(SXI)is part of the scientific payload of the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.SMILE is a joint science mission between the European Space Agency(ESA)and the Chinese Academy of Sciences(CAS)and is due for launch in 2025.SXI is a compact X-ray telescope with a wide field-of-view(FOV)capable of encompassing large portions of Earth’s magnetosphere from the vantage point of the SMILE orbit.SXI is sensitive to the soft X-rays produced by the Solar Wind Charge eXchange(SWCX)process produced when heavy ions of solar wind origin interact with neutral particles in Earth’s exosphere.SWCX provides a mechanism for boundary detection within the magnetosphere,such as the position of Earth’s magnetopause,because the solar wind heavy ions have a very low density in regions of closed magnetic field lines.The sensitivity of the SXI is such that it can potentially track movements of the magnetopause on timescales of a few minutes and the orbit of SMILE will enable such movements to be tracked for segments lasting many hours.SXI is led by the University of Leicester in the United Kingdom(UK)with collaborating organisations on hardware,software and science support within the UK,Europe,China and the United States.

Finding the magnetopause location using soft X-ray observations and a statistical inverse method

Variability in the location and shape of the dayside magnetopause is attributed to magnetic reconnection,a fundamental process that enables the transfer of mass,energy,and momentum from the solar wind into the magnetosphere.The spatial and temporal properties of the magnetopause,under varying solar and magnetospheric conditions,remain largely unknown because empirical studies using in-situ observations are challenging to interpret.Global wide field-of-view(FOV)imaging is the only means to simultaneously observe the spatial distribution of the plasma properties over the vast dayside magnetospheric region and,subsequently,quantify the energy transport from the interplanetary medium into the terrestrial magnetosphere.Two upcoming missions,ESA/CAS SMILE and NASA’s LEXI will provide wide-field imagery of the dayside magnetosheath in soft X-rays,an emission generated by charge exchange interactions between high charge-state heavy ions of solar wind origin and exospheric neutral atoms.High-cadence two-dimensional observations of the magnetosheath will allow the estimation of dynamic properties of its inner boundary,the magnetopause,and enable studies of its response to changes in the solar wind dynamic pressure and interplanetary magnetic field orientation.This work introduces a statistically-based estimation approach based on inverse theory to estimate the spatial distribution of magnetosheath soft X-ray emissivities and,with this,identify the location of the magnetopause over the Sun−Earth line.To do so,we simulate the magnetosheath structure using the MHD-based OpenGGCM model and generate synthetic soft X-ray images using LEXI’s orbit and attitude information.Our results show that 3-D estimations using the described statistically-based technique are robust against Poisson-distributed shot noise inherent to soft X-ray images.Also,our proposed methodology shows that the accuracy of both three-dimensional(3-D)estimation and the magnetopause standoff distance calculation highly depends on the observational point.

Two methods for separating the magnetospheric solar wind charge exchange soft X-ray emission from the diffuse X-ray background

Solar wind charge exchange(SWCX)is the process of solar wind high-valence ions exchanging charges with neutral components and generating soft X-rays.Recently,detecting the SWCX emission from the magnetosphere is proposed as a new technique to study the magnetosphere using panoramic soft X-ray imaging.To better prepare for the data analysis of upcoming magnetospheric soft X-ray imaging missions,this paper compares the magnetospheric SWCX emission obtained by two methods in an XMM-Newton observation,during which the solar wind changed dramatically.The two methods differ in the data used to fit the diffuse X-ray background(DXB)parameters in spectral analysis.The method adding data from the ROSAT All-Sky Survey(RASS)is called the RASS method.The method using the quiet observation data is called the Quiet method,where quiet observations usually refer to observations made by the same satellite with the same target but under weaker solar wind conditions.Results show that the spectral compositions of magnetospheric SWCX emission obtained by the two methods are very similar,and the changes in intensity over time are highly consistent,although the intensity obtained by the RASS method is about 2.68±0.56 keV cm^(-2)s^(-1)sr^(-1)higher than that obtained by the Quiet method.Since the DXB intensity obtained by the RASS method is about 2.84±0.74 keV cm^(-2)s^(-1)sr^(-1)lower than that obtained by the Quiet method,and the linear correlation coefficient between the difference of SWCX and DXB obtained by the two methods in diffe rent energy band is close to-1,the diffe rences in magnetospheric SWCX can be fully attributed to the diffe rences in the fitted DXB.The difference between the two methods is most significant when the energy is less than 0.7 keV,which is also the main energy band of SWCX emission.In addition,the difference between the two methods is not related to the SWCX intensity and,to some extent,to solar wind conditions,because SWCX intensity typically va ries with the solar wind.In summary,both methods are robust and reliable,and should be considered based on the best available options.

Hybrid-Vlasov simulation of soft X-ray emissions at the Earth’s dayside magnetospheric boundaries

Solar wind charge exchange produces emissions in the soft X-ray energy range which can enable the study of near-Earth space regions such as the magnetopause,the magnetosheath and the polar cusps by remote sensing techniques.The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)and Lunar Environment heliospheric X-ray Imager(LEXI)missions aim to obtain soft Xray images of near-Earth space thanks to their Soft X-ray Imager(SXI)instruments.While earlier modeling works have already simulated soft X-ray images as might be obtained by SMILE SXI during its mission,the numerical models used so far are all based on the magnetohydrodynamics description of the space plasma.To investigate the possible signatures of ion-kinetic-scale processes in soft Xray images,we use for the first time a global hybrid-Vlasov simulation of the geospace from the Vlasiator model.The simulation is driven by fast and tenuous solar wind conditions and purely southward interplanetary magnetic field.We first produce global X-ray images of the dayside near-Earth space by placing a virtual imaging satellite at two different locations,providing meridional and equatorial views.We then analyze regional features present in the images and show that they correspond to signatures in soft X-ray emissions of mirrormode wave structures in the magnetosheath and flux transfer events(FTEs)at the magnetopause.Our results suggest that,although the time scales associated with the motion of those transient phenomena will likely be significantly smaller than the integration time of the SMILE and LEXI imagers,mirror-mode structures and FTEs can cumulatively produce detectable signatures in the soft X-ray images.For instance,a local increase by 30%in the proton density at the dayside magnetopause resulting from the transit of multiple FTEs leads to a 12%enhancement in the line-of-sight-and time-integrated soft X-ray emissivity originating from this region.Likewise,a proton density increase by 14%in the magnetosheath associated with mirror-mode structures can result in an enhancement in the soft X-ray signal by 4%.These are likely conservative estimates,given that the solar wind conditions used in the Vlasiator run can be expected to generate weaker soft X-ray emissions than the more common denser solar wind.These results will contribute to the preparatory work for the SMILE and LEXI missions by providing the community with quantitative estimates of the effects of small-scale,transient phenomena occurring on the dayside.

Estimating the subsolar magnetopause position from soft X-ray images using a low-pass image filter

The Lunar Environment heliospheric X-ray Imager(LEXI)and Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)missions will image the Earth’s dayside magneto pause and cusps in soft X-rays after their respective launches in the near future,to specify glo bal magnetic reconnection modes for varying solar wind conditions.To suppo rt the success of these scientific missions,it is critical to develop techniques that extract the magnetopause locations from the observed soft X-ray images.In this research,we introduce a new geometric equation that calculates the subsolar magnetopause position(RS)from a satellite position,the look direction of the instrument,and the angle at which the X-ray emission is maximized.Two assumptions are used in this method:(1)The look direction where soft X-ray emissions are maximized lies tangent to the magnetopause,and(2)the magnetopause surface near the subsolar point is almost spherical and thus RSis nea rly equal to the radius of the magneto pause curvature.We create synthetic soft X-ray images by using the Open Geospace General Circulation Model(OpenGGCM)global magnetohydrodynamic model,the galactic background,the instrument point spread function,and Poisson noise.We then apply the fast Fourier transform and Gaussian low-pass filte rs to the synthetic images to re move noise and obtain accurate look angles for the soft X-ray pea ks.From the filte red images,we calculate RS and its accuracy for different LEXI locations,look directions,and solar wind densities by using the OpenGGCM subsolar magnetopause location as ground truth.Our method estimates RS with an accuracy of<0.3 RE when the solar wind density exceeds>10 cm-3.The accuracy improves for greater solar wind densities and during southward interplanetary magnetic fields.The method ca ptures the magnetopause motion during southwa rd interplaneta ry magnetic field turnings.Consequently,the technique will enable quantitative analysis of the magnetopause motion and help reveal the dayside reconnection modes for dynamic solar wind conditions.This technique will suppo rt the LEXI and SMILE missions in achieving their scientific o bjectives.

SMILE soft X-ray Imager flight model CCD370 pre-flight device characterisation

Throughout the SMILE mission the satellite will be bombarded by radiation which gradually damages the focal plane devices and degrades their performance.In order to understand the changes of the CCD370s within the soft X-ray Imager,an initial characterisation of the devices has been carried out to give a baseline performance level.Three CCDs have been characterised,the two flight devices and the flight spa re.This has been carried out at the Open University in a bespo ke cleanroom measure ment facility.The results show that there is a cluster of bright pixels in the flight spa re which increases in size with tempe rature.However at the nominal ope rating tempe rature(-120℃) it is within the procure ment specifications.Overall,the devices meet the specifications when ope rating at -120℃ in 6 × 6 binned frame transfer science mode.The se rial charge transfer inefficiency degrades with temperature in full frame mode.However any charge losses are recovered when binning/frame transfer is implemented.

Simulation of the SMILE Soft X-ray Imager response to a southward interplanetary magnetic field turning

The Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)Soft X-ray Imager(SXI)will shine a spotlight on magnetopause dynamics during magnetic reconnection.We simulate an event with a southward interplanetary magnetic field turning and produce SXI count maps with a 5-minute integration time.By making assumptions about the magnetopause shape,we find the magnetopause standoff distance from the count maps and compare it with the one obtained directly from the magnetohydrodynamic(MHD)simulation.The root mean square deviations between the reconstructed and MHD standoff distances do not exceed 0.2 RE(Earth radius)and the maximal difference equals 0.24 RE during the 25-minute interval around the southward turning.

Simultaneous fluorescence and Compton scattering computed tomography based on linear polarization X-ray

Purpose To propose a method for simultaneous fluorescence and Compton scattering computed tomography by using linearly polarized X-rays.Methods Monte Carlo simulations were adopted to demonstrate the feasibility of the proposed method.In the simulations,the phantom is a polytetrafluoroethylene cylinder inside which are cylindrical columns containing aluminum,water,and gold(Au)-loaded water solutions with Au concentrations ranging between 0.5 and 4.0 wt%,and a parallel-hole collimator imaging geometry was adopted.The light source was modeled based on a Thomson scattering X-ray source.The phantom images for both imaging modalities were reconstructed using a maximumlikelihood expectation maximization algorithm.Results Both the X-ray fluorescence computed tomography(XFCT)and Compton scattering computed tomography(CSCT)images of the phantom were accurately reconstructed.A similar attenuation contrast problem for the different cylindrical columns in the phantom can be resolved in the XFCT and CSCT images.The interplay between XFCT and CSCT was analyzed,and the contrast-to-noise ratio(CNR)of the reconstruction was improved by correcting for the mutual influence between the two imaging modalities.Compared with K-edge subtraction imaging,XFCT exhibits a CNR advantage for the phantom.Conclusion Simultaneous XFCT and CSCT can be realized by using linearly polarized X-rays.The synergy between the two imaging modalities would have an important application in cancer radiation therapy.

A machine learning-based strategy for predicting the mechanical strength of coral reef limestone using X-ray computed tomography

Different sedimentary zones in coral reefs lead to significant anisotropy in the pore structure of coral reef limestone(CRL),making it difficult to study mechanical behaviors.With X-ray computed tomography(CT),112 CRL samples were utilized for training the support vector machine(SVM)-,random forest(RF)-,and back propagation neural network(BPNN)-based models,respectively.Simultaneously,the machine learning model was embedded into genetic algorithm(GA)for parameter optimization to effectively predict uniaxial compressive strength(UCS)of CRL.Results indicate that the BPNN model with five hidden layers presents the best training effect in the data set of CRL.The SVM-based model shows a tendency to overfitting in the training set and poor generalization ability in the testing set.The RF-based model is suitable for training CRL samples with large data.Analysis of Pearson correlation coefficient matrix and the percentage increment method of performance metrics shows that the dry density,pore structure,and porosity of CRL are strongly correlated to UCS.However,the P-wave velocity is almost uncorrelated to the UCS,which is significantly distinct from the law for homogenous geomaterials.In addition,the pore tensor proposed in this paper can effectively reflect the pore structure of coral framework limestone(CFL)and coral boulder limestone(CBL),realizing the quantitative characterization of the heterogeneity and anisotropy of pore.The pore tensor provides a feasible idea to establish the relationship between pore structure and mechanical behavior of CRL.

Researches on operating region of Tokamak device with soft X-ray tomography

The structures of three operating regions in HT-6B Tokamak have been studied by soft X-ray tomo-graphic system with high sensibility and high time-space resolution. One of the requisites for forming sawtooth discharge is the effective heating action in the central region. In the sawtooth region there are five evolutional phases and five types of magnetic surface structures correspondingly; that is, the concentric, the eccentric, the double-core, the 'MHD-type' and the 'ultra-MHD type' magnetic surface structures. In the MHD oscillation region, there is a stable 'MHD-type' magnetic surface structure. It consists of a crescent 'hot core' and a circular 'cold bubble' and rotates in the diamagnetic direction of electrons. In the resonant region, the resonant helical field improves the heating status and suppresses the MHD disturbances; therefore the single 'MHD-type' magnetic surface changes into a sawtooth type one

Soft x-ray tomographic reconstruction of Heliotron J plasma for the study of magnetohydrodynamic equilibrium and stability

The preparation for an experimental soft x-ray tomography study on the Heliotron J (H-J) machine is carried out,with the objectives of evaluating the capability of the current soft x-ray tomographic system in terms of the identification of different mode structures and their poloidal rotation,and the axis shift with different plasma and machine parameters,and fixing the physics goals for the experimental study.These preparations were carried out via a simulated soft x-ray data set arising from different plasma conditions,such as magnetic islands,low beta and high beta.Soft x-ray tomography (SXT) is performed by the discrete pixel method including singular value decomposition and Phillips-Tikhonov regularization,to obtain clear and smooth images.The H-J soft x-ray tomography results from simulated soft x-rays for the equilibrium H-J plasma sensed the magnetic axis shift clearly and an estimate was also achieved.Successful reconstruction for mode structure m =1/n =1 was obtained along with the realization of the poloidal rotation of the structure.The reconstruction for the m =2/n =1 mode was not very clear for the current soft x-ray diagnostic design.Effective mode identification was not possible due to the lack of measurements.The SXT from the current soft x-ray diagnostic on H-J,the magnetic axis shift can be estimated and the m =1/n =1 mode can be studied.Study of higher poloidal modes is difficult with the current design.

Identification of m = 2 competent mode of complex magneto-hydro-dynamics activities during internal soft disruption based on singular value decomposition and tomography of soft-X-ray emission on the HT-7 tokamak

In this paper,the singular value decomposition（SVD） method as a filter is applied before the tomographic inversion of soft-X-ray emission.Series of ＇filtered＇ signals including specific chronos and topos are obtained.（Here,chronos and topos are the decomposed spatial vectors and the decomposed temporal vectors,respectively）.Given specific magnetic flux function with coupling m = 1 and m = 2 modes,the line-integrated soft-X-ray signals at all chords have been obtained.Then m = 1 and m = 2 modes have been identified by tomography of simulated ＇filtered＇ signals extracted by the SVD method.Finaly,using the experimental line-integrated soft-X-ray signals,m = 2 competent mode of complex magnetohydrodynamics（MHD） activities during internal soft disruption is observed.This result demonstrates that m = 2 mode plays an important role in internal disruption（Here,m is the poloidal mode number）.

Identification of m = 2 competent mode of complex magneto-hydro-dynamics activities during internal soft disruption based on singular value decomposition and tomography of soft-X-ray emission on the HT-7 tokamak

In this paper,the singular value decomposition(SVD) method as a filter is applied before the tomographic inversion of soft-X-ray emission.Series of ’filtered’ signals including specific chronos and topos are obtained.(Here,chronos and topos are the decomposed spatial vectors and the decomposed temporal vectors,respectively).Given specific magnetic flux function with coupling m = 1 and m = 2 modes,the line-integrated soft-X-ray signals at all chords have been obtained.Then m = 1 and m = 2 modes have been identified by tomography of simulated ’filtered’ signals extracted by the SVD method.Finaly,using the experimental line-integrated soft-X-ray signals,m = 2 competent mode of complex magnetohydrodynamics(MHD) activities during internal soft disruption is observed.This result demonstrates that m = 2 mode plays an important role in internal disruption(Here,m is the poloidal mode number).

Use of X-ray computed tomography to study structures and particle contacts of granite residual soil

A small problem about soil particle regularization and contacts but essential to geotechnical engineering was studied.The soils sourced from Guangzhou and Xiamen were sieved into five different particle scale ranges(d<0.075 mm,0.075 mm≤d<0.1 mm,0.1 mm≤d<0.2 mm,0.2 mm≤d<0.5 mm and 0.5 mm≤d<1.0 mm)to study the structures and particle contacts of granite residual soil.The X-ray micro computed tomography method was used to reconstruct the microstructure of granite residual soil.The particle was identified and regularized using principal component analysis(PCA).The particle contacts and geometrical characteristics in 3D space were analyzed and summarized using statistical analyses.The results demonstrate that the main types of contact among the particles are face-face,face-angle,face-edge,edge-edge,edge-angle and angle-angle contacts for particle sizes less than 0.2 mm.When the particle sizes are greater than 0.2 mm,the contacts are effectively summarized as face-face,face-angle,face-edge,edge-edge,edge-angle,angle-angle,sphere-sphere,sphere-face,sphere-edge and sphere-angle contacts.The differences in porosity among the original sample,reconstructed sample and regularized sample are closely related to the water-swelling and water-disintegrable characteristics of granite residual soil.

Use of high-resolution X-ray computed tomography and 3D image analysis to quantify mineral dissemination and pore space in oxide copper ore particles

Mineral dissemination and pore space distribution in ore particles are important features that influence heap leaching performance. To quantify the mineral dissemination and pore space distribution of an ore particle, a cylindrical copper oxide ore sample (I center dot 4.6 mm x 5.6 mm) was scanned using high-resolution X-ray computed tomography (HRXCT), a nondestructive imaging technology, at a spatial resolution of 4.85 mu m. Combined with three-dimensional (3D) image analysis techniques, the main mineral phases and pore space were segmented and the volume fraction of each phase was calculated. In addition, the mass fraction of each mineral phase was estimated and the result was validated with that obtained using traditional techniques. Furthermore, the pore phase features, including the pore size distribution, pore surface area, pore fractal dimension, pore centerline, and the pore connectivity, were investigated quantitatively. The pore space analysis results indicate that the pore size distribution closely fits a log-normal distribution and that the pore space morphology is complicated, with a large surface area and low connectivity. This study demonstrates that the combination of HRXCT and 3D image analysis is an effective tool for acquiring 3D mineralogical and pore structural data.

Effects of dry-wet cycles on three-dimensional pore structure and permeability characteristics of granite residual soil using X-ray micro computed tomography

Due to seasonal climate alterations,the microstructure and permeability of granite residual soil are easily affected by multiple dry-wet cycles.The X-ray micro computed tomography(micro-CT)acted as a nondestructive tool for characterizing the microstructure of soil samples exposed to a range of damage levels induced by dry-wet cycles.Subsequently,the variations of pore distribution and permeability due to drywet cycling effects were revealed based on three-dimensional(3D)pore distribution analysis and seepage simulations.According to the results,granite residual soils could be separated into four different components,namely,pores,clay,quartz,and hematite,from micro-CT images.The reconstructed 3D pore models dynamically demonstrated the expanding and connecting patterns of pore structures during drywet cycles.The values of porosity and connectivity are positively correlated with the number of dry-wet cycles,which were expressed by exponential and linear functions,respectively.The pore volume probability distribution curves of granite residual soil coincide with the χ^(2)distribution curve,which verifies the effectiveness of the assumption of χ^(2)distribution probability.The pore volume distribution curves suggest that the pores in soils were divided into four types based on their volumes,i.e.micropores,mesopores,macropores,and cracks.From a quantitative and visual perspective,considerable small pores are gradually transformed into cracks with a large volume and a high connectivity.Under the action of dry-wet cycles,the number of seepage flow streamlines which contribute to water permeation in seepage simulation increases distinctly,as well as the permeability and hydraulic conductivity.The calculated hydraulic conductivity is comparable with measured ones with an acceptable error margin in general,verifying the accuracy of seepage simulations based on micro-CT results.

Microstructure Evolution of Sandstone Cemented by Microbe Cement Using X-ray Computed Tomography

The bio-sandstone, which was cemented by microbe cement, was firstly prepared, and then the microstructure evolution process was studied by X-ray computed tomography （X-CT） technique. The experimental results indicate that the microstructure of bio-sandstone becomes dense with the development of age. The evolution of inner structure at different positions is different due to the different contents of microbial induced precipitation calcite. Besides, the increase rate of microbial induced precipitation calcite gradually decreases because of the reduction of microbe absorption content with the decreasing pore size in bio-sandstone.