Fast and Accurate Predictor-Corrector Methods Using Feedback-Accelerated Picard Iteration for Strongly Nonlinear Problems

Although predictor-corrector methods have been extensively applied,they might not meet the requirements of practical applications and engineering tasks,particularly when high accuracy and efficiency are necessary.A novel class of correctors based on feedback-accelerated Picard iteration(FAPI)is proposed to further enhance computational performance.With optimal feedback terms that do not require inversion of matrices,significantly faster convergence speed and higher numerical accuracy are achieved by these correctors compared with their counterparts;however,the computational complexities are comparably low.These advantages enable nonlinear engineering problems to be solved quickly and accurately,even with rough initial guesses from elementary predictors.The proposed method offers flexibility,enabling the use of the generated correctors for either bulk processing of collocation nodes in a domain or successive corrections of a single node in a finite difference approach.In our method,the functional formulas of FAPI are discretized into numerical forms using the collocation approach.These collocated iteration formulas can directly solve nonlinear problems,but they may require significant computational resources because of the manipulation of high-dimensionalmatrices.To address this,the collocated iteration formulas are further converted into finite difference forms,enabling the design of lightweight predictor-corrector algorithms for real-time computation.The generality of the proposed method is illustrated by deriving new correctors for three commonly employed finite-difference approaches:the modified Euler approach,the Adams-Bashforth-Moulton approach,and the implicit Runge-Kutta approach.Subsequently,the updated approaches are tested in solving strongly nonlinear problems,including the Matthieu equation,the Duffing equation,and the low-earth-orbit tracking problem.The numerical findings confirm the computational accuracy and efficiency of the derived predictor-corrector algorithms.

Calculation of Mass Concrete Temperature Containing Cooling Water Pipe Based on Substructure and Iteration Algorithm

Mathematical physics equations are often utilized to describe physical phenomena in various fields of science and engineering.One such equation is the Fourier equation,which is a commonly used and effective method for evaluating the effectiveness of temperature control measures for mass concrete.One important measure for temperature control in mass concrete is the use of cooling water pipes.However,the mismatch of grids between large-scale concrete models and small-scale cooling pipe models can result in a significant waste of calculation time when using the finite element method.Moreover,the temperature of the water in the cooling pipe needs to be iteratively calculated during the thermal transfer process.The substructure method can effectively solve this problem,and it has been validated by scholars.The Abaqus/Python secondary development technology provides engineers with enough flexibility to combine the substructure method with an iteration algorithm,which enables the creation of a parametric modeling calculation for cooling water pipes.This paper proposes such a method,which involves iterating the water pipe boundary and establishing the water pipe unit substructure to numerically simulate the concrete temperature field that contains a cooling water pipe.To verify the feasibility and accuracy of the proposed method,two classic numerical examples were analyzed.The results showed that this method has good applicability in cooling pipe calculations.When the value of the iteration parameterαis 0.4,the boundary temperature of the cooling water pipes can meet the accuracy requirements after 4∼5 iterations,effectively improving the computational efficiency.Overall,this approach provides a useful tool for engineers to analyze the temperature control measures accurately and efficiently for mass concrete,such as cooling water pipes,using Abaqus/Python secondary development.

Value Iteration-Based Cooperative Adaptive Optimal Control for Multi-Player Differential Games With Incomplete Information

This paper presents a novel cooperative value iteration(VI)-based adaptive dynamic programming method for multi-player differential game models with a convergence proof.The players are divided into two groups in the learning process and adapt their policies sequentially.Our method removes the dependence of admissible initial policies,which is one of the main drawbacks of the PI-based frameworks.Furthermore,this algorithm enables the players to adapt their control policies without full knowledge of others’ system parameters or control laws.The efficacy of our method is illustrated by three examples.

Another SSOR Iteration Method

Kellogg gave a version of the Peaceman-Radford method. In this paper, we introduce a SSOR iteration method which uses Kellogg’s method. The new algorithm has some advantages over the traditional SSOR algorithm. A Cyclic Reduction algorithm is introduced via a decoupling in Kellogg’s method.

Iteration and evaluation of shale oil development technology for continental rift lake basins

By benchmarking with the iteration of drilling technology,fracturing technology and well placement mode for shale oil and gas development in the United States and considering the geological characteristics and development difficulties of shale oil in the Jiyang continental rift lake basin,East China,the development technology system suitable for the geological characteristics of shale oil in continental rift lake basins has been primarily formed through innovation and iteration of the development,drilling and fracturing technologies.The technology system supports the rapid growth of shale oil production and reduces the development investment cost.By comparing it with the shale oil development technology in the United States,the prospect of the shale oil development technology iteration in continental rift lake basins is proposed.It is suggested to continuously strengthen the overall three-dimensional development,improve the precision level of engineering technology,upgrade the engineering technical indicator system,accelerate the intelligent optimization of engineering equipment,explore the application of complex structure wells,form a whole-process integrated quality management system from design to implementation,and constantly innovate the concept and technology of shale oil development,so as to promote the realization of extensive,beneficial and high-quality development of shale oil in continental rift lake basins.

Variational Iteration Method for Solving Time Fractional Burgers Equation Using Maple

The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.

Policy Iteration for Optimal Control of Discrete-Time Time-Varying Nonlinear Systems

Aimed at infinite horizon optimal control problems of discrete time-varying nonlinear systems,in this paper,a new iterative adaptive dynamic programming algorithm,which is the discrete-time time-varying policy iteration(DTTV)algorithm,is developed.The iterative control law is designed to update the iterative value function which approximates the index function of optimal performance.The admissibility of the iterative control law is analyzed.The results show that the iterative value function is non-increasingly convergent to the Bellman-equation optimal solution.To implement the algorithm,neural networks are employed and a new implementation structure is established,which avoids solving the generalized Bellman equation in each iteration.Finally,the optimal control laws for torsional pendulum and inverted pendulum systems are obtained by using the DTTV policy iteration algorithm,where the mass and pendulum bar length are permitted to be time-varying parameters.The effectiveness of the developed method is illustrated by numerical results and comparisons.

全钨ITER朗缪尔探针的高热负荷测试及损伤行为

本文研究了全钨ITER朗缪尔探针在高热负荷条件下的耐热性能、绝缘性能和损伤行为。利用核工业西南物理研究院的60 kW电子束测试平台,对探针在10 MW/m^(2)稳态热负荷、20 MW/m^(2)瞬态热负荷及15 MW/m^(2)高参数稳态热负荷下进行模拟测试。实验结果表明,全钨朗缪尔探针能够承受10 MW/m^(2)热负荷100分钟、20 MW/m^(2)热负荷300次循环以及15 MW/m^(2)热负荷3000次循环。在15 MW/m^(2)条件下,探针表面在前1000次循环中无明显损伤,但在3000次循环后出现裂纹,且部分探针表面发生再结晶现象。研究表明,现有全钨朗缪尔探针的设计能够满足当前ITER装置的服役条件,但面对未来DEMO堆更加苛刻的条件可能还需要进一步优化。

Iteration dependent interval based open‐closed‐loop iterative learning control for time varying systems with vector relative degree

For linear time varying(LTV)multiple input multiple output(MIMO)systems with vector relative degree,an open‐closed‐loop iterative learning control(ILC)strategy is developed in this article,where the time interval of operation is iteration dependent.To compensate the missing tracking signal caused by iteration dependent interval,the feedback control is introduced in ILC design.As the tracking signal of many continuous iterations is lost in a certain interval,the feedback control part can employ the tracking signal of current iteration for compensation.Under the assumption that the initial state vibrates around the desired initial state uniformly in mathematical expectation sense,the expectation of ILC tracking error can converge to zero as the number of iteration tends to infinity.Under the circumstance that the initial state varies around the desired initial state with a bound,as the number of iteration tends to infinity,the expectation of ILC tracking error can be driven to a bounded range,whose upper bound is proportional to the fluctuation.It is revealed that the convergence condition is dependent on the feed-forward control gains,while the feedback control can accelerate convergence speed by selecting appropriate feedback control gains.As a special case,the controlled system with integrated high relative degree is also addressed by proposing a simplified iteration dependent interval based open‐closed‐loop ILC method.Finally,the effectiveness of the developed iteration dependent interval based open‐closed‐loop ILC is illustrated by a simulation example with two cases on initial state.

A Novel Method to Enhance the Inversion Speed and Precision of the NMR T_(2) Spectrum by the TSVD Based Linearized Bregman Iteration

The low-field nuclear magnetic resonance(NMR)technique has been used to probe the pore size distribution and the fluid composition in geophysical prospecting and related fields.However,the speed and accuracy of the existing numerical inversion methods are still challenging due to the ill-posed nature of the first kind Fredholm integral equation and the contamination of the noises.This paper proposes a novel inversion algorithmto accelerate the convergence and enhance the precision using empirical truncated singular value decompositions(TSVD)and the linearized Bregman iteration.The L1 penalty term is applied to construct the objective function,and then the linearized Bregman iteration is utilized to obtain fast convergence.To reduce the complexity of the computation,empirical TSVD is proposed to compress the kernel matrix and determine the appropriate truncated position.This novel inversion method is validated using numerical simulations.The results indicate that the proposed novel method is significantly efficient and can achieve quick and effective data solutions with low signal-to-noise ratios.

ITER气体注入系统注入管线和气体流量控制阀特性研究

为验证国际热核聚变实验堆计划(ITER)的气体注入系统(GIS)的设计是否能够满足装置运行的要求,以及为优化气体注入系统的气体加料控制模式提供直接而可靠的数据支撑。该研究在基于ITER GIS最终设计阶段提供的设计方案的基础上,在模拟实验平台上进行实验研究,该研究涉及燃料气体H2、D2以及杂质气体Ne、Ar、N2及其混合气体的注入特性,包括其供气压力、注入管线结构对响应时间的影响。实验结论证明,气体注入系统的设计可以满足ITER装置的运行要求,且实验数据为ITER装置的气体控制模式,尤其是实现杂质气体快速、精确的控制,从而保障ITER装置安全运行,提供有力的技术支持。

A novel three-step implicit iteration process for threeΦ-hemicontractive mapping in the intermediate sense

In this paper,we introduce a three-step composite implicit iteration process for approximating the common fixed point of three uniformly continuous and asymptotically generalizedΦ-hemicontractive mappings in the intermediate sense.We prove that our proposed iteration process converges to the common fixed point of three finite family of asymptotically generalizedΦ-hemicontractive mappings in the intermediate sense.Our results extends,improves and complements several known results in literature.

Data-Driven Active Disturbance Rejection Control of Plant-Protection Unmanned Ground Vehicle Prototype: A Fuzzy Indirect Iterative Learning Approach

Dear Editor,This letter proposes a fuzzy indirect iterative learning(FIIL)active disturbance rejection control(ADRC)scheme to address the impact of uncertain factors of plant-protection unmanned ground vehicle(UGV),in which ADRC is a data-driven model-free control algorithm that only relies on the input and output data of the system.Based on the established nonlinear time-varying dynamic model including dynamic load(medicine box),the FIIL technology is adopted to turn the bandwidth and control channel gain online,in which the fuzzy logic system is used to update the gain parameters of iterative learning in real time.Simulation and experiment show the FIIL-ADRC scheme has better control performance.

西物院ITER项目强矩阵管理模式探索与成效

ITER计划是目前全球规模最大、影响最深远的国际大科学工程,由欧盟、中国、美国、日本、俄罗斯、韩国、印度等七方共35个国家参与,也是我国参与的最大的国际合作项目,我国承担任务在所有七方中完成最好,多次受到ITER总干事的表扬。文章主要针对西物院ITER项目强矩阵管理模式探索进行了总结,希望为未来参与大科学工程项目提供借鉴。

DQ变换和MUSIC算法在ITER磁体电源信号间谐波检测中的应用

随着国际热核聚变实验堆(International Thermonuclear Experimental Reactor,ITER)计划的逐步开展,保证ITER磁体电源系统的稳定运行显得尤为重要。文章采用将DQ变换和多信号分类(multiple signal classification,MUSIC)算法相结合的方法进行间谐波频率检测,信号的幅度和相位由最小二乘法来估计。DQ变换可以消除大幅度ITER基波分量,MUSIC算法可以通过矩阵特征分解检测出短数据条件下的谐波和间谐波,适用短时平稳的间谐波检测,两者相结合可以有效检测出大幅度基波附近存在小幅度间谐波。仿真实验表明,计算经DQ变换后检测出的ITER信号谐波频率时,取中间信号计算真实频谱较为正确,两侧信号则有较大的误差。

Gyrokinetic simulations of the kinetic electron effects on the electrostatic instabilities on the ITER baseline scenario

The linear and nonlinear simulations are carried out using the gyrokinetic code NLT for the electrostatic instabilities in the core region of a deuterium plasma based on the International Thermonuclear Experimental Reactor(ITER)baseline scenario.The kinetic electron effects on the linear frequency and nonlinear transport are studied by adopting the adiabatic electron model and the fully drift-kinetic electron model in the NLT code,respectively.The linear simulations focus on the dependence of linear frequency on the plasma parameters,such as the ion and electron temperature gradientsκ_(Ti,e)≡R=L_(Ti,e),the density gradientκ_(n)≡R/L_(n)and the ion-electron temperature ratioτ=T_(e)=T_(i).Here,is the major radius,and T_(e)and T_(i)denote the electron and ion temperatures,respectively.L_(A)=-(δ_(r)lnA)^(-1)is the gradient scale length,with denoting the density,the ion and electron temperatures,respectively.In the kinetic electron model,the ion temperature gradient(ITG)instability and the trapped electron mode(TEM)dominate in the small and large k_(θ)region,respectively,wherek_(θ)is the poloidal wavenumber.The TEMdominant region becomes wider by increasing(decreasing)κ_(T_(e))(κ_(T_(i)))or by decreasingκ_(n).For the nominal parameters of the ITER baseline scenario,the maximum growth rate of dominant ITG instability in the kinetic electron model is about three times larger than that in the adiabatic electron model.The normalized linear frequency depends on the value ofτ,rather than the value of T_(e)or T_(i),in both the adiabatic and kinetic electron models.The nonlinear simulation results show that the ion heat diffusivity in the kinetic electron model is quite a lot larger than that in the adiabatic electron model,the radial structure is finer and the time oscillation is more rapid.In addition,the magnitude of the fluctuated potential at the saturated stage peaks in the ITGdominated region,and contributions from the TEM(dominating in the higher k_(θ)region)to the nonlinear transport can be neglected.In the adiabatic electron model,the zonal radial electric field is found to be mainly driven by the turbulent energy flux,and the contribution of turbulent poloidal Reynolds stress is quite small due to the toroidal shielding effect.However,in the kinetic electron model,the turbulent energy flux is not strong enough to drive the zonal radial electric field in the nonlinear saturated stage.The kinetic electron effects on the mechanism of the turbulence-driven zonal radial electric field should be further investigated.

Fault Estimation for a Class of Markov Jump Piecewise-Affine Systems: Current Feedback Based Iterative Learning Approach

In this paper, the issues of stochastic stability analysis and fault estimation are investigated for a class of continuoustime Markov jump piecewise-affine(PWA) systems against actuator and sensor faults. Firstly, a novel mode-dependent PWA iterative learning observer with current feedback is designed to estimate the system states and faults, simultaneously, which contains both the previous iteration information and the current feedback mechanism. The auxiliary feedback channel optimizes the response speed of the observer, therefore the estimation error would converge to zero rapidly. Then, sufficient conditions for stochastic stability with guaranteed performance are demonstrated for the estimation error system, and the equivalence relations between the system information and the estimated information can be established via iterative accumulating representation.Finally, two illustrative examples containing a class of tunnel diode circuit systems are presented to fully demonstrate the effectiveness and superiority of the proposed iterative learning observer with current feedback.

On traceable iterated line graph and hamiltonian path index

Xiong and Liu[21]gave a characterization of the graphs G for which the n-iterated line graph L^(n)(G)is hamiltonian,for n≥2.In this paper,we study the existence of a hamiltonian path in L^(n)(G),and give a characterization of G for which L^(n)(G)has a hamiltonian path.As applications,we use this characterization to give several upper bounds on the hamiltonian path index of a graph.

Transient Stability Analysis of Converter-based Islanded Microgrids Based on Iterative Equal Area Criterion Considering Reactive Power Loop Dynamics and Varying Damping

With the rapid increase in the installed capacity of renewable energy in modern power systems,the stable operation of power systems with considerable power electronic equipment requires further investigation.In converter-based islanded mi-crogrid(CIM)systems equipped with grid-following(GFL)and grid-forming(GFM)voltage-source converters(VSCs),it is chal-lenging to maintain stability due to the mutual coupling effects between different VSCs and the loss of voltage and frequency sup-port from the power system.In previous studies,quantitative transient stability analysis was primarily used to assess the active power loop of GFM-VSCs.However,frequency and voltage dy-namics are found to be strongly coupled,which strongly affects the estimation result of stability boundary.In addition,the vary-ing damping terms have not been fully captured.To bridge these gaps,this paper investigates the transient stability of CIM consid-ering reactive power loop dynamics and varying damping.First,an accuracy-enhanced nonlinear model of the CIM is derived based on the effects of reactive power loop and post-disturbance frequency jump phenomena.Considering these effects will elimi-nates the risk of misjudgment.The reactive power loop dynamics make the model coefficients be no longer constant and thus vary with the power angle.To evaluate quantitatively the effects of re-active power loop and varying damping on the transient stability of CIM,an iterative criterion based on the equal area criterion theory is proposed.In addition,the effects of parameters on the stable boundary of power system are analyzed,and the dynamic interaction mechanisms are revealed.Simulation and experiment results verify the merits of the proposed method.

Iterative physical optics method based on efficient occlusion judgment with bounding volume hierarchy technology

This paper builds a binary tree for the target based on the bounding volume hierarchy technology,thereby achieving strict acceleration of the shadow judgment process and reducing the computational complexity from the original O(N^(3))to O(N^(2)logN).Numerical results show that the proposed method is more efficient than the traditional method.It is verified in multiple examples that the proposed method can complete the convergence of the current.Moreover,the proposed method avoids the error of judging the lit-shadow relationship based on the normal vector,which is beneficial to current iteration and convergence.Compared with the brute force method,the current method can improve the simulation efficiency by 2 orders of magnitude.The proposed method is more suitable for scattering problems in electrically large cavities and complex scenarios.