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.

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.

考虑工序相关性的动态Job shop调度问题启发式算法

提出一类考虑工序相关性的、工件批量到达的动态Job shop调度问题,在对工序相关性进行了定义和数学描述的基础上,进一步建立了动态Job shop调度问题的优化模型。设计了一种组合式调度规则RAN(FCFS,ODD),并提出了基于规则的启发式算法以及该类动态Job shop调度问题的算例生成方法。为验证算法和比较评估调度规则的性能,对算例采用文献提出的7种调度规则和RAN(FCFS,ODD)进行了仿真调度,对调度结果的分析表明了算法的有效性和RAN(FCFS,ODD)调度规则求解所提出的动态Job Shop调度问题的优越性能。

基于遗传算法求解Job Shop调度优化的新方法

针对Job Shop调度问题,提出了一种遗传算法编码新方法和矩阵解码方法。该方法根据问题的特点,采用一种按工序进行总体排序染色体编码方案,并采用矩阵解码,解码时体现了编码与调度方案一一对应,并且该编码方案有多种交叉操作算子可用,不需要专门设计算子。算例计算结果表明,基于该编码方案的遗传算法是有效的,能适用解决Job Shop调度问题,通过比较,用该编码方案的遗传算法优化Job Shop调度操作简单并且收敛速度快。

基于遗传算法求解Job Shop调度的编码新方法

针对Job Shop调度问题,提出了一种新的遗传算法编码新方法。该方法根据问题的特点,采用一种按工序用不同编号进行的染色体编码方案,每一个编号包含工件工序号、对应的机器号、加工时间等所有信息,此编码与调度方案一一对应,并且该编码方案有多种交叉操作算子可用,不需要专门设计算子。算例计算结果表明,基于该编码方案的遗传算法是有效的,能适用解决Job Shop调度问题。通过比较,用该编码方案的遗传算法优化Job Shop调度操作简单并且收敛速度快。

调整时间与工序顺序相关的Job Shop调度优化方法研究

以多目标多约束Job Shop调度问题为背景,考虑工序调整、作业相似性等特点,利用成组技术和人工智能树,提出了与工序排序有关的调整时间描述方法及其调度规则,给出了此问题的混合整数规划模型,并构造了基于生物免疫机理的调度算法以提高求解此类问题的效果。介绍了保证抗体的多样性和有效性的基于自然数的抗体编码方法和基于抗体浓度的抗体激增和抑制策略。用毛纺后整理车间的一个调度实例和La31 案例证明了免疫算法用于求解复杂生产调度问题的速度和效果。

基于重排S-method的多分量辐射源信号分析方法

为实现多分量辐射源信号的有效检测和识别,提出了基于重排S-method(RSM)的多分量辐射源信号分析方法.针对Wigner-Ville分布交叉项严重和计算复杂度高的问题,采用RSM分析密集、交叠的辐射源信号,给出了其算法的硬件实现和计算复杂度分析.该方法能处理线性和非线性调频信号,算法简单,易于实现.复杂体制多分量辐射源信号分析的实验结果表明,该方法能有效分析多分量辐射源信号,时频分辨率高,不受交叉项干扰,具有比Wigner-Ville分布和S-method更强的噪声抑制能力.

基于S-method的多分量雷达辐射源信号检测

针对传统时频方法在处理多分量雷达辐射源信号时存在交叉项,不能检测各分量信号时域参数,难以适应低信噪比环境的问题,提出一种基于S-method(SM)的多分量雷达辐射源信号检测新方法。该方法首先计算信号的SM时频分布,然后在时频面的基础上检测各信号分量的瞬时频率和脉冲起止时间。实验结果表明,该方法能处理线性及非线性调频信号、时频分辨率高且不受交叉项干扰,时域检测精度大于98.60%,频域检测精度大于99.48%,信噪比降低时仍然保持强检测能力。

Comparative Analysis of Improving Effects of 12 Maize Populations with Improved S_1 Selection Method

ln this study, 12 maize populations were improved with the improved S1 selection method, aiming to increase the improvement efficiency of maize popula-tions. The results showed that number of rows per spike, number of grains per row and 100-grain weight were the three important component traits of maize yield. The highest genetic increment was found in Mengqun 2, fol owed by Mengqun 4, Mengqun 1 and other 7 maize populations. Negative genetic increment was shown in Mengqun 3 and 3 introduced foreign maize populations. Some changes were shown in spikes, plant traits and genetic diversity of maize populations. Based on our results, we concluded that Mengqun 2, Mengqun 4, Zhongzong 7, Mengqun C and Mengqun A could be directly used for the line breeding by selfing for their higher genetic increment of yield and better improvement effects of other agricultural traits. Compared with these 5 populations, the improvement potential of other maize populations was limited for their lower genetic increment.

一种求解Job-Shop调度问题的快速禁忌搜索算法

介绍了一种解决Job-Shop调度最短完工时间的有效的快速禁忌搜索算法,该算法首先利用插入算法构造尽可能好的初始解,然后使用禁忌搜索算法改进当前解,用基准实例进行仿真,实验结果表明该算法是可行的和有效的.

Stoichiometry and Stability Constant Values for Copper (II) Chelates with Ethylene Diamine in Deep Eutectic Solvents (DES) (Ethaline) Solutions

In this study we used the deep eutectic solvents (ionic liquids) to investigate the reaction between copper (II) with ethylene diamine (en). Two of the existing methods for analyzing spectrophotometric measurements have been applied for establishing, the stoichiometry and whenever possible, the stability constants of the chelates formed. The method of continuous variations was necessary to determine first whether, the metal ion and the ligand ethylene diamine form one or more than one chelate, when more than one chelate formed, the results obtained depend on the wavelength and for meaningful conclusions the wavelengths were carefully selected. The empirical formulae of the chelates were further substantiated by the molar ratio method. The effect of time and temperature on the formation and stability of these chelates in solution is also studied. The stability constants, K1 and K2 for the copper (II) chelates were calculated, though reliable, and are comparable to literature values.

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算法更优。

基于S-R和分解定理的二维几何非线性问题的虚单元法求解

应变-旋转(Strain-Rotation,S-R)和分解定理为分析几何非线性问题提供了合理可靠的理论基础,但用有限元求解时会遇到大变形发生后的网格畸变问题。近年提出的虚单元法(Virtual element method,VEM)适用于一般的多边形网格,因此,该文尝试使用一阶虚单元求解基于S-R和分解定理的二维几何非线性问题,以克服网格畸变的影响。基于重新定义的多项式位移空间基函数,推演获得一阶虚单元分析线弹性力学问题时允许位移空间向多项式位移空间的投影表达式;按照虚单元法双线性格式的计算规则,分析处理基于更新拖带坐标法和势能率原理的增量变分方程;进而建立离散系统方程及其矩阵表达形式,并编制MATLAB求解程序;采用常规多边形网格和畸变网格,应用该文算法分析均布荷载下的悬臂梁和均匀内压下的厚壁圆筒变形。结果与已有文献和ANSYS软件的对比表明:该文算法在两种网格中均可有效执行且具备足够数值精度。总体该文算法为基于S-R和分解定理的二维几何非线性问题求解提供了一种鲁棒方法。

基于S-Method分布的微多普勒特征分析

微多普勒特征是雷达目标所具有的独特特征之一,对目标的分类、识别具有特殊的意义。研究高精度时频分析方法在分析目标微动特性中的作用,可以为后续目标识别提供很好的支撑。S-Method分布作为一种新型的时频分析方法,它基于短时傅里叶变换来实现,减少了分析过程中的运算量,同时能较好地解决交叉项问题。首先对弹道导弹弹头的微动模型进行建模,推导得到微动模型的理论微多普勒频率,然后采用S-Method分布对回波信号进行时频分析仿真实验,获得弹头目标章动的高精度的时间-微多普勒频率图。通过比较其在分析过程中的时频分辨率、交叉项,具体阐述S-Method分布在时频分析中的优势。因此可以将S-Method分布应用于雷达目标微多普勒分析中,分析实时变化的微多普勒频率特征。

Q-factor estimation in CMP gather and the continuous spectral ratio slope method

The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.

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

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

广义概率密度演化方程的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)条件的限制,因此该方法具有良好的数值稳定性和计算效率.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.