Calculation and prediction of divertor detachment via impurity seeding by using one-dimensional model

Achieving the detachment of divertor can help to alleviate excessive heat load and sputtering problems on the target plates,thereby extending the lifetime of divertor components for fusion devices.In order to provide a fast but relatively reliable prediction of plasma parameters along the flux tube for future device design,a one-dimensional(1D)modeling code for the operating point of impurity seeded detached divertor is developed based on Python language,which is a fluid model based on previous work(Plasma Phys.Control.Fusion 58045013(2016)).The experimental observation of the onset of divertor detachment by neon(Ne)and argon(Ar)seeding in EAST is well reproduced by using the 1D modeling code.The comparison between the 1D modeling and two-dimensional(2D)simulation by the SOLPS-ITER code for CFETR detachment operation with Ne and Ar seeding also shows that they are in good agreement.We also predict the radiative power loss and corresponding impurity concentration requirement for achieving divertor detachment via different impurity seeding under high heating power conditions in EAST and CFETR phase II by using the 1D model.Based on the predictions,the optimized parameter space for divertor detachment operation on EAST and CFETR is also determined.Such a simple but reliable 1D model can provide a reasonable parameter input for a detailed and accurate analysis by 2D or three-dimensional(3D)modeling tools through rapid parameter scanning.

A study on the numerical prediction method for the vertical thermal structure in the Bohai Sea and the Huanghai Sea-I.One-dimensional numerical prediction model

In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.

Reconstructing the upper ocean thermal profiles using one-dimensional numerical model

The observation data for 5 d at a station in the South China Sea is presented. After brief anMysis of the wind speed, air temperature from the ship-borne meteorological instruments and temperature and salinity profiles from the CTD （conductivity, temperature, depth recorder） data, the authors find that the CTD casts are too sparse for us to understand the diurnal evolution of the thermal structure in the upper ocean. A one-dimensional （1D） numericM code based on Mellor-Yamada turbulence closure model is used to reconstruct the upper ocean thermal structure, utilizing the atmospheric forcing data from ship-borne weather station. The simulation results show good agreement with the observational data; the significance of breaking waves is also briefly discussed. The evolution of turbulence kinetic energy （TKE） and the contribution from shear production and buoy- ancy production are discussed respectively. Finally, several possible factors which might influence the numerical results are briefly analyzed.

A Dugdale-Barenblatt model for elliptical orifice problem with asymmetric cracks in one-dimensional orthorhombic quasicrystals

By means of Muskhelishvili’s method and the technique of generalized conformal mapping,the physical plane problems are transformed into regular mathematical problems in quasicrystals(QCs).The analytical solution to an elliptical orifice problem with asymmetric cracks in one-dimensional(1D)orthorhombic QCs is obtained.By using the Dugdale-Barenblatt model,the plastic simulation at the crack tip of the elliptical orifice with asymmetric cracks in 1D orthorhombic QCs is performed.Finally,the size of the atomic cohesive force zone is determined precisely,and the size of the atomic cohesive force zone around the crack tip of an elliptical orifice with a single crack or two symmetric cracks is obtained.

Research on the Application of the Radiative Transfer Model Based on Deep Neural Network in One-dimensional Variational Algorithm

As a typical physical retrieval algorithm for retrieving atmospheric parameters,one-dimensional variational(1 DVAR)algorithm is widely used in various climate and meteorological communities and enjoys an important position in the field of microwave remote sensing.Among algorithm parameters affecting the performance of the 1 DVAR algorithm,the accuracy of the microwave radiative transfer model for calculating the simulated brightness temperature is the fundamental constraint on the retrieval accuracies of the 1 DVAR algorithm for retrieving atmospheric parameters.In this study,a deep neural network(DNN)is used to describe the nonlinear relationship between atmospheric parameters and satellite-based microwave radiometer observations,and a DNN-based radiative transfer model is developed and applied to the 1 DVAR algorithm to carry out retrieval experiments of the atmospheric temperature and humidity profiles.The retrieval results of the temperature and humidity profiles from the Microwave Humidity and Temperature Sounder(MWHTS)onboard the Feng-Yun-3(FY-3)satellite show that the DNN-based radiative transfer model can obtain higher accuracy for simulating MWHTS observations than that of the operational radiative transfer model RTTOV,and also enables the 1 DVAR algorithm to obtain higher retrieval accuracies of the temperature and humidity profiles.In this study,the DNN-based radiative transfer model applied to the 1 DVAR algorithm can fundamentally improve the retrieval accuracies of atmospheric parameters,which may provide important reference for various applied studies in atmospheric sciences.

Simulation of the Influence of Ion-Produced NO_X and HO_X Radicals on the Antarctic Ozone Depletion with a One-Dimensional Model

A one -dimensional time-dependent photochemical model is used to simulate the influence of ion-produced NOx and HOx radicals on the Antarctic ozone depletion in polar night and polar spring at a latitude of 73 degrees south.Vertical transport and nitrogen-oxygen (NOx). hydrogen-oxygen (HOx) production by ionic reactions have been introduced into the model.NOx and HOx produced by precipitating ions are transported into the lower stratosphere by vertical motion and have some effects in the development of the Antarctic ozone depletion.From winter through spring the calculated ozone column decreases to 269.4 DU. However, this value is significantly higher than the total ozone observed at several Antarctic ozone stations.

Entanglement in the Quantum Phase Transition of the Half-Integer Spin One-Dimensional Heisenberg Model

We use the Bethe’s ansatz method to study the entanglement of spinons in the quantum phase transition of half integer spin one-dimensional magnetic chains known as quantum wires. We calculate the entanglement in the limit of the number of particles . We obtain an abrupt change in the entanglement next the quantum phase transition point of the anisotropy parameter ?from the gapped phase ?to gapless phase .

Atmosphere-Ocean Coupling Schemes in a One-Dimensional Climate Model

In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.

A SOLVABLE ONE-DIMENSIONAL ALLOCATION DECISION MODEL OF RENEWABLE RESOURCES

In the paper, the determinate atlecation decision model and the probabilistic allocation decision model of a kind of renewable resource are separatly studied by means of dynamic programming, and the optimal allocation policy is given under some special conditions.

Process Modeling and Dynamic Simulation for EAST Helium Refrigerator

In this paper,the process modeling and dynamic simulation for the EAST helium refrigerator has been completed.The cryogenic process model is described and the main components are customized in detail.The process model is controlled by the PLC simulator,and the realtime communication between the process model and the controllers is achieved by a customized interface.Validation of the process model has been confirmed based on EAST experimental data during the cool down process of 300-80 K.Simulation results indicate that this process simulator is able to reproduce dynamic behaviors of the EAST helium refrigerator very well for the operation of long pulsed plasma discharge.The cryogenic process simulator based on control architecture is available for operation optimization and control design of EAST cryogenic systems to cope with the long pulsed heat loads in the future.

Evolution of Bond-Order-Wave Phase in One-Dimensional Mott Insulators

In this paper, by using the level spectroscopy method and bosonization theory, we discuss the evolution of the bond-order-wave （BOW） phase in a one-dimensional half-filled extended Hubbard model wlth the on-site Coulomb repulsion U as well as the inter-site Coulomb repulsion V and antiferromagnetic exchange J. After clarifying the generic phase diagrams in three limiting cases with one of the parameters being fixed at zero individually, we find that the BOW phase in the U-V phase diagram is initially enlarged as J increases from zero but is eventually suppressed as J increases further in the strong-coupling regime. A three-dimensional phase diagram is suggested where the BOW phase exists in an extended region separated from the spin-density-wave and charge-density-wave phases.

One-dimensional horizontal infiltration experiment for determining permeability coefficient of loamy sand

A knowledge of soil permeability is essential to evaluate hydrologic characteristics of soil, such as water storage and water movement, and soil permeability coefficient is an important parameter that reflects soil permeability. In order to confirm the acceptability of the one-dimensional horizontal infiltration method（one-D method） for simultaneously determining both the saturated and unsaturated permeability coefficients of loamy sand, we first measured the cumulative infiltration and the wetting front distance under various infiltration heads through a series of one-dimensional horizontal infiltration experiments, and then analyzed the relationships of the cumulative horizontal infiltration with the wetting front distance and the square root of infiltration time. We finally compared the permeability results from Gardner model based on the one-D method with the results from other two commonly-used methods（i.e., constant head method and van Genuchten model） to evaluate the acceptability and applicability of the one-D method. The results showed that there was a robust linear relationship between the cumulative horizontal infiltration and the wetting front distance, suggesting that it is more appropriate to take the soil moisture content after infiltration in the entire wetted zone as the average soil moisture content than as the saturated soil moisture content. The results also showed that there was a robust linear relationship between the cumulative horizontal infiltration and the square root of infiltration time, suggesting that the Philip infiltration formula can better reflect the characteristics of cumulative horizontal infiltration under different infiltration heads. The following two facts indicate that it is feasible to use the one-D method for simultaneously determining the saturated and unsaturated permeability coefficients of loamy sand. First, the saturated permeability coefficient（prescribed in the Gardner model） of loamy sand obtained from the one-D method well agreed with the value obtained from the constant head method. Second, the relationship of unsaturated permeability coefficient with soil water suction for loamy sand calculated using Gardner model based on the one-D method was nearly identical with the same relationship calculated using van Genuchten model.

Modeling tungsten response under helium plasma irradiation:a review

Tungsten,a leading candidate for plasma-facing materials(PFM) in future fusion devices,will be exposed to high-flux low-energy helium plasma under the anticipated fusion operation conditions.In the past two decades,experiments have revealed that exposure to helium plasma strongly modifies the surface morphology and hence the sputtering,thermal and other properties of tungsten,posing a serious danger to the performance and lifetime of tungsten and the steadystate operation of plasma.In this article,we provide a review of modeling and simulation efforts on the long-term evolution of helium bubbles,surface morphology,and property changes of tungsten exposed to low-energy helium plasma.The current gap and outstanding challenges to establish a predictive modeling capability for dynamic evolution of PFM are discussed.

Anomalous energy diffusion and heat conduction in one-dimensional system

We propose a new concept, the centre of energy, to study energy diffusion and heat conduction in a one-dimensional hard-point model. For the diatom model, we find an anomalous energy diffusion as （x2） - tβ with β = 1.33, which is independent of initial condition and mass rate. The present model can be viewed as the model composed by independent quasi-particles, the centre of energy. In this way, heat current can be calculated. Based on the theory of dynamic billiard, the divergent exponent of heat conductivity is estimated to be α = 0.33, which is confirmed by a simple numerical calculation.

Theoretical Research on Scattering Resonance States of Reaction I＋HI（v=O）→IH（v＇=0）＋I： Partial Potential Energy Surface and One-dimensional Quantum Reactive Scattering Calculation

Based on the vibrational potential curves coupled with the minimum energy reaction path, the partial potential energy surface of the reaction I＋HI→IH＋I was constructed at the QCISD（T）//MP4SDQ level with pseudo potential method. And the formation mechanism of the scattering resonance states of this reaction was well interpreted with the partial potential energy surface. The scattering resonance states of this reaction should belong to Feshbach resonance because of the coupling of the vibrational mode and the translational mode. With the one-dimensional square potential well model, the resonance width and lifetime of the I＋HI（v=0）→IH（v＇=0）＋I state-to-state reaction were calculated, which preferably explained the high-resolved threshold photodetachment spectroscopy of the IHI- anion performed by Neumark et al..

Mathematical models and expert system for grate-kiln process of iron ore oxide pellet production.Part Ⅱ:Rotary kiln process control

Rotary kiln process for iron ore oxide pellet production is hard to detect and control.Construction of one-dimensional model of temperature field in rotary kiln was described.And the results lay a solid foundation for online control.Establishment of kiln process control expert system was presented,with maximum temperature of pellet and gas temperature at the feed end as control cores,and interval estimate as control strategy.Software was developed and put into application in a pellet plant.The results show that control guidance of this system is accurate and effective.After production application for nearly one year,the compressive strength and first grade rate of pellet are increased by 86 N and 2.54%,respectively,while FeO content is 0.05% lowered.This system can reveal detailed information of real time kiln process,and provide a powerful tool for online control of pellet production.

The Accumulation of He on a W Surface During keV-He Irradiation:Cluster Dynamics Modeling

The accumulation of He on a W surface during keV-He ion irradiation has been simulated using cluster dynamics modeling. This is based mainly on rate theory and improved by involving different types of objects, adopting up-to-date parameters and complex reaction processes, as well as considering the diffusion process along with depth. These new features make the simulated results compare very well with the experimental ones. The accumulation and diffusion processes are analyzed, and the depth and size dependence of the He concentrations contributed by different types of He clusters is also discussed. The exploration of the trapping and diffusion effects of the He atoms is helpful in understanding the evolution of the damages in the near-surface of plasma-facing materials under He ion irradiation.

Performance Optimization of Torque Converters Based on Modified 1D Flow Model

A methodology for performance optimization of torque converters is put forward based on the one-dimensional (1D) flow model. It is found that the inaccuracy of 1D flow model for predicting hydraulic performance at the low speed ratio is mainly caused by the separation phenomenon at the stator cascade which is induced by large flow impinging at the pressure side of the stator blades. A semi-empirical separation model is presented and incorporated to the original 1D flow model. It is illustrated that the improved model is able to predict the circumferential velocity components accurately, which can be applied to performance optimization. Then, the Pareto front is obtained by using the genetic algorithm (GA) in order to inspect the coupled relationship among stalling impeller torque capacity, stalling torque ratio and efficiency. The efficiency is maximized on the premise that a target stalling impeller torque capacity and torque ratio are achieved. Finally, the optimized result is verified by the computational fluid dynamics(CFD) simulation, which indicates that the maximal efficiency is increased by 0.96%.

Numerical simulation of a direct current glow discharge in atmospheric pressure helium

Characteristics of a direct current （DC） discharge in atmospheric pressure helium are numerically investigated based on a one-dimensional fluid model. The results indicate that the discharge does not reach its steady state till it takes a period of time. Moreover, the required time increases and the current density of the steady state decreases with increasing the gap width. Through analyzing the spatial distributions of the electron density, the ion density and the electric field at different discharge moments, it is found that the DC discharge starts with a Townsend regime, then transits to a glow regime. In addition, the discharge operates in a normal glow mode or an abnormal glow one under different parameters, such as the gap width, the ballast resistors, and the secondary electron emission coefficients, judged by its voltage-current characteristics.

Numerical Simulation for One Dimensional Steady Quasineutral Hybrid Model of Stationary Plasma Thruster

Based on the analysis of the physical mechanism of the Stationary Plasma Thruster （SPT）, an integral equation describing the ion density of the steady SPT and the ion velocity distribution function at an arbitrary axial position of the steady SPT channel are derived. The integral equation is equivalent to the Vlasov equation, but the former is simpler than the latter. A one dimensional steady quasineutral hybrid model is established. In this model, ions are described by the above integral equation, and neutrals and electrons are described by hydrodynamic equations. The transferred equivalency to the differential equation and the integral equation, together with other equations, are solved by an ordinary differential equation （ODE） solver in the Matlab. The numerical simulation results show that under various circumstances, the ion average velocity would be different and needs to be deduced separately.