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.

3D seismic forward modeling from the multiphysical inversion at the Ketzin CO_(2) storage site

From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.

UAV-based transient electromagnetic 3D forward modeling and inversion and analysis of exploration capability

Unmanned aerial vehicle transient electromagnetic(UAV-TEM)is a novel airborne exploration method that offers advantages such as low cost,simple operation,high exploration efficiency and suitability for near-surface exploration in complex terrain areas.To improve the accuracy of data interpretation in this method,the authors conducted a systematic three-dimensional(3D)forward modeling and inversion of the UAV-TEM.This study utilized the finite element method based on unstructured tetrahedral elements and employed the second-order backward Euler method for time discretization.This allowed for accurate 3D modeling and accounted for the effects of complex terrain.Based on these,the influence characteristics of flight altitudes and the sizes,burial depths,and resistivities of anomalies are compared and analyzed to explore the UAV-TEM systems’exploration capability.Lastly,four typical geoelectrical models of landslides are designed,and the inversion method based on the Gauss-Newton optimization method is used to image the landslide models and analyze the imaging effect of the UAV-TEM method on landslide geohazards.Numerical results showed that UAV-TEM could have better exploration resolution and fine imaging of nearsurface structures,providing important technical support for monitoring,early warning,and preventing landslides and other geological hazards.

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.

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

The application of a gravimetric forward modelling of the lithospheric structure for an estimate of the average density of the upper asthenosphere

The average density of 3300 kg m-3 is often attributed for the asthenosphere. In this study, we inspect this value by estimating the average value of the(upper) asthenosphere based on applying the gravimetric forward modelling of major known lithospheric density structures. The LITHO1.0 global seismic model of the lithospheric density structure is used for this purpose, while considering that the lithosphere-asthenosphere boundary(LAB) is rheological, conventionally taken at the 1300C isotherm,above which the mantle behaves in a rigid fashion and below which it behaves in a ductile fashion.According to our result, the average density of the upper asthenosphere is roughly 3400 kg m-3. This density value closely agrees with the corresponding average value 3371 kg m-3 computed based on an empirical density model provided in the Preliminary Reference Earth Model(PREM), while using the LITHO1.0 LAB depth data. We also demonstrate that the sub-lithospheric mantle gravity map exhibits mainly a thermal signature. The most prominent features in this gravity map are mid-oceanic spearing ridges marked by gravity lows, while oceanic subductions in the West Pacific are characterized by the most pronounced gravity highs.

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.

Estimation of Displacement and Extension due to Reverse Drag of Normal Faults: Forward Method

In the case of reverse drag of normal faulting, the displacement and horizontal extension are determined based on the established equations for the three mechanisms: rigid body, vertical shear and inclined shear. There are three sub-cases of basal detachment for the rigid body model: horizontal detachment, antithetic detachment and synthetic detachment. For the rigid body model, the established equations indicate that the total displacement on the synthetic base (Dt2) is the largest, that on the horizontal base (Dt1) is moderate, and that on the antithetic base (Dt3) is the smallest. On the other hand, the value of (Dt1) is larger than the displacement for the vertical shear (Dt4). The value of (Dt1) is larger than or less than the displacement for the inclined shear (Dt5) depending on the original fault dip δ0, bedding angle θ, and the angle of shear direction β. For all original parameters, the value of Dt5 is less than the value of Dt4. Also, by comparing three rotation mechanisms, we find that the inclined shear produces largest extension, the rigid body model with horizontal detachment produces the smallest extension, and the vertical shear model produces moderate extension.

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.

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.

Variable-coordinate forward modeling of irregular surface based on dual-variable grid

The mapping method is a forward-modeling method that transforms the irregular surface to horizontal by mapping the rectangular grid as curved; moreover, the wave field calculations move from the physical domain to the calculation domain. The mapping method deals with the irregular surface and the low-velocity layer underneath it using a fine grid. For the deeper high-velocity layers, the use of a fine grid causes local oversampling. In addition, when the irregular surface is transformed to horizontal, the flattened interface below the surface is transformed to curved, which produces inaccurate modeling results because of the presence of ladder-like burrs in the simulated seismic wave. Thus, we propose the mapping method based on the dual-variable finite-difference staggered grid. The proposed method uses different size grid spacings in different regions and locally variable time steps to match the size variability of grid spacings. Numerical examples suggest that the proposed method requires less memory storage capacity and improves the computational efficiency compared with forward modeling methods based on the conventional grid.

Multicomponent seismic forward modeling of gas hydrates beneath the seafloor

We investigated the effect of microscopic distribution modes of hydrates in porous sediments, and the saturation of hydrates and free gas on the elastic properties of saturated sediments. We simulated the propagation of seismic waves in gas hydrate-bearing sediments beneath the seafloor, and obtained the common receiver gathers of compressional waves（P-waves） and shear waves（S-waves）. The numerical results suggest that the interface between sediments containing gas hydrates and free gas produces a large-amplitude bottomsimulating reflector. The analysis of multicomponent common receiver data suggests that ocean-bottom seismometers receive the converted waves of upgoing P- and S-waves, which increases the complexity of the wavefield record.

基于Forward-Modeling方法的黑河流域水储量变化特征研究

黑河流域陆地水储量变化对流域下游等周边区域水资源的合理利用以及经济和社会发展等有着重要的意义.本文利用2003年1月至2013年12月的GRACE RL05数据反演了黑河流域陆地水储量长时间序列的变化,并针对重力场模型和数据处理中产生的信号泄漏问题,采用Forward-Modeling方法进行了改正并恢复泄漏信号;将GRACE获得的泄漏信号恢复前后的黑河流域水储量变化结果与全球水文模型GLDAS和CPC进行比较分析,结果表明泄漏信号改正后的结果与水文模型结果的时间序列相关性均有明显提高,从其空间分布结果可以看出Forward-Modeling方法有效地恢复初始信号、增强被湮没的信号,泄漏信号误差减小;通过分析黑河流域水储量变化的长时间序列结果,发现其具有明显的阶段性变化特征,即2003—2006年呈明显下降趋势,约为-0.86cm·a-1,在2007—2010年趋于平衡状态,而2011—2013年则呈现缓慢上升趋势约为0.14cm·a-1;联合GRACE数据和GLDAS数据反演了黑河流域地下水储量变化,并与全球降雨数据GPCC进行了比较分析,两者相关性可达到0.88以上.

Forward modeling of marine DC resistivity method for a layered anisotropic earth

Since the ocean bottom is a sedimentary environment wherein stratification is well developed, the use of an anisotropic model is best for studying its geology. Beginning with Maxwell＇s equations for an anisotropic model, we introduce scalar potentials based on the divergence-free characteristic of the electric and magnetic （EM） fields. We then continue the EM fields down into the deep earth and upward into the seawater and couple them at the ocean bottom to the transmitting source. By studying both the DC apparent resistivity curves and their polar plots, we can resolve the anisotropy of the ocean bottom. Forward modeling of a high-resistivity thin layer in an anisotropic half-space demonstrates that the marine DC resistivity method in shallow water is very sensitive to the resistive reservoir but is not influenced by airwaves. As such, it is very suitable for oil and gas exploration in shallowwater areas but, to date, most modeling algorithms for studying marine DC resistivity are based on isotropic models. In this paper, we investigate one-dimensional anisotropic forward modeling for marine DC resistivity method, prove the algorithm to have high accuracy, and thus provide a theoretical basis for 2D and 3D forward modeling.

AVO forwarding modeling in two-phase media: multiconstrained matrix mineral modulus inversion

AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.

Behavioral modeling of RF power amplifiers with time-delay feed-forward neural networks

A novel behavioral model using three-layer time-delay feed-forward neural networks （TDFFNN）is adopted to model radio frequency （RF）power amplifiers exhibiting memory nonlinearities. In order to extract the parameters, the back- propagation algorithm is applied to train the proposed neural networks. The proposed model is verified by the typical odd- order-only memory polynomial model in simulation, and the performance is compared with different numbers of taped delay lines（TDLs） and perceptrons of the hidden layer. For validating the TDFFNN model by experiments, a digital test bench is set up to collect input and output data of power amplifiers at a 60 × 10^6 sample/s sampling rate. The 3.75 MHz 16-QAM signal generated in the vector signal generator（VSG） is chosen as the input signal, when measuring the dynamic AM/AM and AM/PM characteristics of power amplifiers. By comparisons and analyses, the presented model provides a good performance in convergence, accuracy and efficiency, which is approved by simulation results and experimental results in the time domain and frequency domain.