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 MOVING WATER EQUILIBRIA PRESERVING NONSTAGGERED CENTRAL SCHEME ACHIEVED VIA FLUX GLOBALIZATION FOR THE RIPA MODEL

In this paper,we propose a second-order moving-water equilibria preserving nonstaggered central scheme to solve the Ripa model via flux globalization.To maintain the moving-water steady states,we use the discrete source terms proposed by Britton et al.(J Sci Comput,2020,82(2):Art 30)by incorporating the expression of the source terms as a whole into the flux gradient,which directly avoids the discrete complexity of the source terms in order to maintain the well-balanced properties of the scheme.In addition,since the nonstaggered central scheme requires re-projecting the updated values of the nonstaggered cells onto the staggered cells,we modify the calculation of the global variables by constructing ghost cells and alternating the values of the global variables with the water depths obtained from the solution through the nonlinear relationship between the global flux and the water depth.In order to maintain the second-order accuracy of the scheme on the time scale,we incorporate a new Runge-Kutta type time discretization in the evolution of the numerical solution for the nonstaggered cells.Meanwhile,we introduce the"draining"time step technique to ensure that the water depth is positive and that it satisfies mass conservation.Numerical experiments verify that the scheme is well-balanced,positivity-preserving and robust.

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.

Modeling of surface flux in Tongyu using the Simple Biosphere Model 2 (SiB2)

The modeling ability of a stand-alone version of the Simple Biosphere Model 2（SiB2） was tested mainly through diagnosing the simulated latent heat（LE）,sensible heat（H）,CO2 flux,and air temperature at the Tongyu field observation station（44°25＇N,122°52＇E,184 m elevation） of Coordinated Enhanced Observing Period（CEOP）,where the land cover is cropland and grassland.In the whole year of 2003,the canopy height and the leaf area index was variable.During non-growth period,the surface would become bare,while during the growth period,the canopy height could reach 2.0 m high over cropland and 0.8 m high over grassland,respectively,and max leaf area index could reach 4.2 and 2.4,respectively.The model was initialized with measurement and driven by half-hourly atmospheric observations.The simulation values for 2003 were compared against measurements.Results show that the model is of a good ability of simulating the hourly latent heat（LE）,sensible heat（H）,CO2 flux and temperature during the growth period.Moreover,the daily LE,H and CO2 flux simulated by SiB2 could reflect their yearly change reasonably.However,the model may overestimate the H generally.

A discrete model for prediction of radon flux from fractured rocks

Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.

Inter-comparisons of thermodynamic sea-ice modeling results using various parameterizations of radiative flux

Radiative fluxes are of primary importance in the energy and mass balance of the sea-ice cover. Various parameterizations of the radiative fluxes are studied in a thermodynamic sea-ice model. Model outputs of the surface radiative and heat fluxes and mass balance are compared with observations. The contribution of short-wave radiation is limited to a long part of winter. Therefore, simple schemes are often sufficient. Errors in estimations of the short-wave radiation are due mainly to cloud effects and occasionally to multi-reflection between surface and ice crystals in the air. The long-wave radiation plays an important role in the ice surface heat and mass balance during most part of a winter. The effect of clouds on the accuracy of the simple radiative schemes is critical, which needs further attention. In general, the accuracy of an ice model depends on that of the radiative fluxes.

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.

Evaluation of CMIP5 Climate Models in Simulating 1979–2005 Oceanic Latent Heat Flux over the Pacific

The climatological mean state, seasonal variation and long-term upward trend of 1979-2005 latent heat flux （LHF） in historical runs of 14 coupled general circulation models from CMIP5 （Coupled Model Intercomparison Project Phase 5） are evaluated against OAFlux （Objectively Analyzed air-sea Fluxes） data. Inter-model diversity of these models in simulating the annual mean climatological LHF is discussed. Results show that the models can capture the climatological LHF fairly well, but the amplitudes are generally overestimated. Model-simulated seasonal variations of LHF match well with observations with overestimated amplitudes. The possible origins of these biases are wind speed biases in the CMIP5 models. Inter-model diversity analysis shows that the overall stronger or weaker LHF over the tropical and subtropical Pacific region, and the meridional variability of LHF, are the two most notable diversities of the CMIP5 models. Regression analysis indicates that the inter-model diversity may come from the diversity of simulated SST and near-surface atmospheric specific humidity. Comparing the observed long-term upward trend, the trends of LHF and wind speed are largely underestimated, while trends of SST and air specific humidity are grossly overestimated, which may be the origins of the model biases in reproducing the trend of LHF.

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.

Calibration and validation of SiBcrop Model for simulating LAI and surface heat fluxes of winter wheat in the North China Plain

The accurate representation of surface characteristic is an important process to simulate surface energy and water flux in land-atmosphere boundary layer.Coupling crop growth model in land surface model is an important method to accurately express the surface characteristics and biophysical processes in farmland.However,the previous work mainly focused on crops in single cropping system,less work was done in multiple cropping systems.This article described how to modify the sub-model in the SiBcrop to realize the accuracy simulation of leaf area index(LAI),latent heat flux(LHF)and sensible heat flux(SHF)of winter wheat growing in double cropping system in the North China Plain(NCP).The seeding date of winter wheat was firstly reset according to the actual growing environment in the NCP.The phenophases,LAI and heat fluxes in 2004–2006 at Yucheng Station,Shandong Province,China were used to calibrate the model.The validations of LHF and SHF were based on the measurements at Yucheng Station in 2007–2010 and at Guantao Station,Hebei Province,China in 2009–2010.The results showed the significant accuracy of the calibrated model in simulating these variables,with which the R2,root mean square error(RMSE)and index of agreement(IOA)between simulated and observed variables were obviously improved than the original code.The sensitivities of the above variables to seeding date were also displayed to further explain the simulation error of the SiBcrop Model.Overall,the research results indicated the modified SiBcrop Model can be applied to simulate the growth and flux process of winter wheat growing in double cropping system in the NCP.

Multi-layer Quasi Three-dimensional Equivalent Model of Axial-Flux Permanent Magnet Synchronous Machine

Axial-flux permanent magnet synchronous machine(AFPMSM)enjoys the merits of high torque density and high efficiency,which make it one good candidate in the direct-drive application.The AFPMSM is usually analyzed based on the three-dimensional finite element method(3D FEM)due to its three-dimensional magnetic field distribution.However,the 3D FEM suffers large amount of calculation,time-consuming and is not suitable for the optimization of AFPMSM.Addressing this issue,a multi-layer quasi three-dimensional equivalent model of the AFPMSM is investigated in this paper,which could take the end leakage into consideration.Firstly,the multi-layer quasi three-dimensional equivalent model of the AFPMSM with single stator and single rotor is derived in details,including the equivalent processes and conversions of structure dimensions,motion conditions and electromagnetic parameters.Then,to consider the influence of end leakage on the performance,a correction factor is introduced in the multi-layer quasi three-dimensional equivalent model.Finally,the proposed multi-layer quasi three-dimensional equivalent model is verified by the 3D FEM based on an AFPMSM under different structure parameters.It demonstrates that the errors of flux linkage and average torque obtained by the multi-layer quasi three-dimensional equivalent model and 3D FEM are only around 2%although the structure parameters of the AFPMSM are varied.Besides,the computation time of one case based on the multi-layer quasi three-dimensional equivalent model is only 6 min,which is much less than that of the 3D FEM,1.8 h,under the same conditions.Thus,the proposed multi-layer quasi three-dimensional equivalent model could be used to optimize the AFPMSM and much time could be saved by this method compared with the 3D FEM.

Simulation of internal nitrogen release from bottom sediments in an urban lake using a nitrogen flux model

Nutrient release from sediment is considered a significant source for overlying water. Given that nutrient release mechanisms in sediment are complex and difficult to simulate, traditional approaches commonly use assigned parameter values to simulate these processes. In this study, a nitrogen flux model was developed and coupled with the water quality model of an urban lake. After parameter sensitivity analyses and model calibration and validation, this model was used to simulate nitrogen exchange at the sediment–water interface in eight scenarios. The results showed that sediment acted as a buffer in the sediment–water system. It could store or release nitrogen at any time, regulate the distribution of nitrogen between sediment and the water column, and provide algae with nitrogen. The most effective way to reduce nitrogen levels in urban lakes within a short time is to reduce external nitrogen loadings. However, sediment release might continue to contribute to the water column until a new balance is achieved. Therefore, effective measures for reducing sediment nitrogen should be developed as supplementary measures. Furthermore, model parameter sensitivity should be individually examined for different research subjects.

A data assimilation-based forecast model of outer radiation belt electron fluxes

Because radiation belt electrons can pose a potential threat to the safety of satellites orbiting in space,it is of great importance to develop a reliable model that can predict the highly dynamic variations in outer radiation belt electron fluxes.In the present study,we develop a forecast model of radiation belt electron fluxes based on the data assimilation method,in terms of Van Allen Probe measurements combined with three-dimensional radiation belt numerical simulations.Our forecast model can cover the entire outer radiation belt with a high temporal resolution(1 hour)and a spatial resolution of 0.25 L over a wide range of both electron energy(0.1-5.0 MeV)and pitch angle(5°-90°).On the basis of this model,we forecast hourly electron fluxes for the next 1,2,and 3 days during an intense geomagnetic storm and evaluate the corresponding prediction performance.Our model can reasonably predict the stormtime evolution of radiation belt electrons with high prediction efficiency(up to~0.8-1).The best prediction performance is found for~0.3-3 MeV electrons at L=~3.25-4.5,which extends to higher L and lower energies with increasing pitch angle.Our results demonstrate that the forecast model developed can be a powerful tool to predict the spatiotemporal changes in outer radiation belt electron fluxes,and the model has both scientific significance and practical implications.