Improvement of Euler＇s Method Using Particle Swarm Optimization

Many problems in applied mathematics lead to ordinary differential equation. In this paper, a considerable refinement and improvement of the Euler＇s method obtained using PSO （particle swarm optimization） was presented. PSO is a technique based on the cooperation between particles. The exchange of information between these particles allows to resolve difficult problems. This approach is carefully handled and tested with an illustrated example.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

The forecasting efficiency under different selected regions by Pattern Informatics Method and seismic potential estimation in the North-South Seismic Zone

In 2022,four earthquakes with M_(S)≥6.0 including the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes occurred in the North-South Seismic Zone(NSSZ),which demonstrated high and strong seismicity.Pattern Informatics(PI)method,as an effective long and medium term earthquake forecasting method,has been applied to the strong earthquake forecasting in Chinese mainland and results have shown the positive performance.The earthquake catalog with magnitude above M_(S)3.0 since 1970 provided by China Earthquake Networks Center was employed in this study and the Receiver Operating Characteristic(ROC)method was applied to test the forecasting efficiency of the PI method in each selected region related to the North-South Seismic Zone systematically.Based on this,we selected the area with the best ROC testing result and analyzed the evolution process of the PI hotspot map reflecting the small seismic activity pattern prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes.A“forward”forecast for the area was carried out to assess seismic risk.The study shows the following.1)PI forecasting has higher forecasting efficiency in the selected study region where the difference of seismicity in any place of the region is smaller.2)In areas with smaller differences of seismicity,the activity pattern of small earthquakes prior to the Menyuan M_(S)6.9 and Luding M_(S)6.8 earthquakes can be obtained by analyzing the spatio-temporal evolution process of the PI hotspot map.3)The hotspot evolution in and around the southern Tazang fault in the study area is similar to that prior to the strong earthquakes,which suggests the possible seismic hazard in the future.This study could provide some ideas to the seismic hazard assessment in other regions with high seismicity,such as Japan,Californi,Turkey,and Indonesia.

Research on the Application of Montessori Education Method in Cognitive Training of Patients with Alzheimer’s Disease

Objective:To study the application of the Montessori education method in cognitive training in patients with Alzheimer’s disease(AD).Methods:40 cases of senile dementia patients who were admitted to our hospital from January 2022 to January 2023 were selected and randomly divided into an intervention group and a control group according to the single and double number table method,with 20 cases in each group.The intervention group used the Montessori education method,the principle of which was to implement individualized health interventions based on the individual conditions of the patients,for a period of 6 months;the control group was given conventional treatment and nursing of the disease.The Mini-Mental State Examination(MMSE)was used to compare the effects of the two groups of patients before and after health intervention and conduct statistical analysis.Results:The score of the intervention group was higher than that of the control group,and there was a statistical difference between the two(P<0.05).Conclusion:Implementing the Montessori education method for diagnosed Alzheimer’s patients can effectively improve their cognitive function and delay the progress of further dementia.

ADASYN与类别逆比例加权法在阿尔茨海默病不平衡数据中的应用

目的利用自适应合成抽样(adaptive synthetic sampling,ADASYN)与类别逆比例加权法处理类别不平衡数据,结合分类器构建模型对阿尔茨海默病(alzheimer′s disease,AD)患者疾病进程进行分类预测。方法数据源自阿尔茨海默病神经影像学计划(Alzheimer′s disease neuroimaging initiative,ADNI),经随机森林填补缺失值,弹性网络筛选特征子集后,利用ADASYN与类别逆比例加权法处理类别不平衡数据。分别结合随机森林(random forest,RF)、支持向量机(support vector machine,SVM)构建四种模型:ADASYN-RF、ADASYN-SVM、加权随机森林(weighted random forest,WRF)、加权支持向量机(weighted support vector machine,WSVM),与RF、SVM比较分类性能。模型评价指标为宏观平均精确率(macro-average of precision,macro-P)、宏观平均召回率(macro-average of recall,macro-R)、宏观平均F1值(macro-average of F1-score,macro-F1)、准确率(accuracy,ACC)、Kappa值和AUC(area under the ROC curve)。结果ADASYN-RF的分类性能最优(Kappa值为0.938,AUC为0.980),ADASYN-SVM次之。利用ADASYN-RF预测得到的重要分类特征分别为CDRSB、LDELTOTAL、MMSE,在临床上均可得到证实。结论ADASYN与类别逆比例加权法都能辅助提升分类器性能,但ADASYN算法更优。

基于蒙特卡洛法的Euler-Bernoulli梁基频和振型求解方法

将Rayleigh法和蒙特卡洛法相结合,在Euler-Bernoulli梁理论假设下求解了均匀梁、变截面梁和附带集中质量的变截面梁自由振动问题.对原本连续的梁结构模型进行离散化处理,利用蒙特卡洛法给出梁结构的假设振型.将假设得到的梁结构振型函数代入Rayleigh法,多次计算过程中,将历次基频所得值与计算所得最小值进行比较,根据其相对误差判断是否满足收敛条件,进而求得基频及对应的振型.结果表明,不同计算模型中基频最大误差不超过10%,能够满足工程需求,且精度和时间的控制参数调整灵活,使用者可根据自身需要自行调节.该方法理论简明,适用范围广泛,能够快速准确地求解诸多类型的梁结构基频和振型.

广义概率密度演化方程的Chebyshev拟谱法

概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation,GDEE).对于GDEE的求解通常采用有限差分法,然而,由于GDEE是初始条件间断的变系数一阶双曲偏微分方程,通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象.文章从全局逼近的角度出发,基于Chebyshev拟谱法为GDEE构造了全局插值格式,解决了数值色散、数值耗散以及网格敏感性问题.考虑GDEE的系数在每个时间步长均为常数,推导了GDEE在每一个时间步长内时域上的序列矩阵指数解.由于序列矩阵指数解形式上是解析的,从而很好地克服了数值稳定性问题.两个数值算例表明,通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合,且数值耗散和数值色散现象几乎可以忽略.此外,拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制,因此该方法具有良好的数值稳定性和计算效率.

基于S-CAD方法的宅基地“三权分置”改革试点政策评估——以新疆试点市县为例

采用S-CAD政策评估方法对新疆宅基地“三权分置”政策试点地区伊宁市和奇台县的政策执行效果进行评价,研究发现伊宁市和奇台县两试点政策整体逻辑一致性较好、政策要素之间具备充分性与必要性,利益相关者支持程度较高,有效实现了宅基地“三权分置”改革的政策目标,但仍存在政策手段与政策目标充要性不足,预期政策效果与实际政策效果有一定偏差等问题。针对此,提出宅基地“三权分置”改革政策应适应城乡融合发展趋势,加强农户利益保障机制建设;坚持因地制宜,创新农村宅基地“三权分置”改革思路;兼顾多方利益与公平,提升综合价值效应的建议。

电液3-UPS/S并联稳定平台参数振动特性分析

针对电液3-UPS/S并联稳定平台驱动液压缸的压力脉动所产生的参数振动,建立了稳定平台的参数振动方程并利用多尺度法求解了主共振响应与组合共振响应的一次近似解;分析了主共振与组合共振响应特性以及振动幅值在初始工作空间内的变化规律,最后采用四阶龙格库塔法与模态试验对参数振动模型进行验证。结果表明:数值解与理论解之间的最大误差为4.20%,固有频率理论值与试验值之间的最大误差为4.66%,可验证参数振动模型的正确性。

2024年1月23日乌什M_(S)7.1地震序列频谱偏移特征分析

采用P波初动方法计算了乌什M_(S)7.1地震序列中M_(S)≥4.0地震的震源机制解,震源机制解在时间上存在前期一致性较好,后期一致性紊乱的特点。空间上,乌什县内震源机制解类型多为逆冲型,阿合奇县内震源机制解类型多为走滑型。利用新疆地震台网数字波形资料和快速傅里叶变换方法,计算乌什M_(S)7.1地震序列中M_(S)≥5.0地震序列的频率谱,计算结果显示,乌什县内发生的5级余震前期有频谱偏移特征,最后一个5级地震频带变宽,频谱偏移特征消失;阿合奇县内发生的5级余震均出现了频谱偏移特征。

基于D-S证据理论改进AHP-熵权的流域洪涝灾害评估研究

考虑致灾因子危险性、孕灾环境敏感性以及承灾体易损性,选取指标构建小清河流域洪涝灾害风险评估指标体系,提出一种基于D-S证据理论的改进AHP-熵权法计算指标权重,求取洪涝灾害风险指数,运用自然断点分级法确定洪涝灾害风险等级,分析小清河流域洪涝灾害风险空间分布情况。结果表明:小清河流域洪涝灾害风险总体上表现出南低北高的趋势,其中高风险区和较高风险区分别占流域面积的8.7%和14.3%,主要分布在小清河干流以及主要支流两岸。所得评估结果同“利奇马”台风发生期间实际洪灾风险分布情况一致,对比证明基于D-S证据理论的改进AHP-熵权法优于AHP和熵权法,可为小清河流域防洪减灾决策提供依据。

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

Global Existence of Smooth Solutions for the One-Dimensional Full Euler System for a Dusty Gas

We study the existence of global-in-time classical solutions for the one-dimensional nonisentropic compressible Euler system for a dusty gas with large initial data.Using the characteristic decomposition method proposed by Li et al.(Commun Math Phys 267:1–12,2006),we derive a group of characteristic decompositions for the system.Using these characteristic decompositions,we find a sufficient condition on the initial data to ensure the existence of global-in-time classical solutions.

On the Approximation of Fractal-Fractional Differential Equations Using Numerical Inverse Laplace Transform Methods

Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.

The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate

In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.

Geophysical Study: Estimation of Deposit Depth Using Gravimetric Data and Euler Method (Jalalabad Iron Mine, Kerman Province of IRAN)

Mineral exploration is done by different methods. Geophysical and geochemical studies are two powerful tools in this field. In integrated studies, the results of each study are used to determine the location of the drilling boreholes. The purpose of this study is to use field geophysics to calculate the depth of mineral reserve. The study area is located 38 km from Zarand city called Jalalabad iron mine. In this study, gravimetric data were measured and mineral depth was calculated using the Euler method. 1314 readings have been performed in this area. The rocks of the region include volcanic and sedimentary. The source of the mineralization in the area is hydrothermal processes. After gravity measuring in the region, the data were corrected, then various methods such as anomalous map remaining in levels one and two, upward expansion, first and second-degree vertical derivatives, analytical method, and analytical signal were drawn, and finally, the depth of the deposit was estimated by Euler method. As a result, the depth of the mineral deposit was calculated to be between 20 and 30 meters on average.

A Fourth-Order Unstructured NURBS-Enhanced Finite Volume WENO Scheme for Steady Euler Equations in Curved Geometries

In Li and Ren(Int.J.Numer.Methods Fluids 70:742–763,2012),a high-order k-exact WENO finite volume scheme based on secondary reconstructions was proposed to solve the two-dimensional time-dependent Euler equations in a polygonal domain,in which the high-order numerical accuracy and the oscillations-free property can be achieved.In this paper,the method is extended to solve steady state problems imposed in a curved physical domain.The numerical framework consists of a Newton type finite volume method to linearize the nonlinear governing equations,and a geometrical multigrid method to solve the derived linear system.To achieve high-order non-oscillatory numerical solutions,the classical k-exact reconstruction with k=3 and the efficient secondary reconstructions are used to perform the WENO reconstruction for the conservative variables.The non-uniform rational B-splines(NURBS)curve is used to provide an exact or a high-order representation of the curved wall boundary.Furthermore,an enlarged reconstruction patch is constructed for every element of mesh to significantly improve the convergence to steady state.A variety of numerical examples are presented to show the effectiveness and robustness of the proposed method.

基于S-CAD方法的四川省单元式乡村国土空间规划政策评价与优化研究

国土空间规划改革以来,各地开展了探索与实践,2021年四川在全省推行以片区为单元编制乡村国土空间规划,是全国第一个在全省范围开展单元式乡村国土空间规划的省域。通过运用S-CAD方法和政策工具分析法对其政策体系进行分析,并结合规划编制成果对政策进行综合评价,发现其存在以下不足:一是政策体系不完善,主要出台了规划编制相关政策,而在规划实施和监督方面缺乏专门政策;二是政策价值观和政策预期结果超出了政策目标和政策手段所能够实现的范围,政策元素之间具有不一致性;三是政策工具结构不平衡,在规划区域衔接、对上衔接等方面非常薄弱。因此,需要进一步完善政策体系、调适政策元素关系、优化政策工具结构。通过这一研究,以期为深化乡村国土空间规划探索提供借鉴。