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.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

High-Order Solitons and Hybrid Behavior of (3 + 1)-Dimensional Potential Yu-Toda-Sasa-Fukuyama Equation with Variable Coefficients

In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton solution. By using the long wave limit method, the N-order rational solution can be obtained from N-order soliton solution. Then, through the paired complexification of parameters, the lump solution is obtained from N-order rational solution. Meanwhile, we obtained a hybrid solution between 1-lump solution and N-soliton (N=1,2) by using the long wave limit method and parameter complex. Furthermore, four different sets of three-dimensional graphs of solitons, lump solutions and hybrid solutions are drawn by selecting four different sets of coefficient functions which include one set of constant coefficient function and three sets of variable coefficient functions.

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.

Solution of Laguerre’s Differential Equations via Modified Adomian Decomposition Method

This paper presents a technique for obtaining an exact solution for the well-known Laguerre’s differential equations that arise in the modeling of several phenomena in quantum mechanics and engineering. We utilize an efficient procedure based on the modified Adomian decomposition method to obtain closed-form solutions of the Laguerre’s and the associated Laguerre’s differential equations. The proposed technique makes sense as the attitudes of the acquired solutions towards the neighboring singular points are correctly taken care of.

Adomian Modification Methods for the Solution of Chebyshev’s Differential Equations

The current study examines the important class of Chebyshev’s differential equations via the application of the efficient Adomian Decomposition Method (ADM) and its modifications. We have proved the effectiveness of the employed methods by acquiring exact analytical solutions for the governing equations in most cases;while minimal noisy error terms have been observed in a particular method modification. Above all, the presented approaches have rightly affirmed the exactitude of the available literature. More to the point, the application of this methodology could be extended to examine various forms of high-order differential equations, as approximate exact solutions are rapidly attained with less computation stress.

Pythagorician Divisors and Applications to Some Diophantine Equations

We consider the Pythagoras equation X2 +Y2 = Z2, and for any solution of the type (a,b = 2sb1 ≠0,c) ∈ N*3, s ≥ 2, b1odd, (a,b,c) ≡ (±1,0,1)(mod 4), c > a , c > b, and gcd(a,b,c) = 1, we then prove the Pythagorician divisors Theorem, which results in the following: , where (d,d′′) (resp. (e,en)) are unique particular divisors of a and b, such that a = dd′′ (resp. b = ee′′ ), these divisors are called: Pythagorician divisors from a, (resp. from b). Let’s put λ ∈{0,1}, defined by: and S = s -λ (s -1). Then such that . Moreover the map is a bijection. We apply this new tool to obtain a new classification of the primitive, positive and non-trivial solutions of the Pythagoras equations: a2 + b2 = c2 via the Pythagorician parameters (d,e,S ). We obtain for (d, e) fixed, the equivalence class of any Pythagorician solution (a,b,c), checking , namely: . We also update the solutions of some Diophantine equations of degree 2, already known, but very important for the resolution of other equations. With this tool of Pythagorean divisors, we have obtained (in another paper) new recurrent methods to solve Fermat’s equation: a4 + b4 = c4, other than usual infinite descent method;and to solve congruent numbers problem. We believe that this tool can bring new arguments, for Diophantine resolution, of the general equations of Fermat: a2p + b2p = c2p and ap + bp = cp. MSC2020-Mathematical Sciences Classification System: 11A05-11A51-11D25-11D41-11D72.

高压管汇材料疲劳性能测试及P-S-N模型曲线的拟合

高压管汇作为压裂设备中的主要易损件之一,其失效危害较大。它的失效原因主要是疲劳、冲蚀、腐蚀或者材料缺陷引起的刺漏和爆裂,其中尤以疲劳失效最不可预估。目前,对于高压管汇材料的疲劳性能研究不够深入,为解决高压管汇材料疲劳寿命的准确描述问题,以某国产高压管汇材料为例,进行了一系列疲劳试验,并基于试验数据,采用多种分布模型和不同S-N模型进行拟合分析,得出综合评价拟合能力最强的P-S-N模型。结果表明,该材料在中长疲劳寿命区,Weibull三参数模型在7级应力水平下综合评价能力最好;在存活率分别为50%、90%、99%、99.9%时,指数S-N模型的拟合系数均大于0.98,拟合能力最好。得出的P-S-N模型曲线可以为高压管汇的疲劳寿命以及安全设计提供依据。

A Family of Global Attractors for the Generalized Kirchhoff-Beam Equations

In this paper, we discuss the existence and uniqueness of global solutions, the existence of the family of global attractors and its dimension estimation for generalized Beam-Kirchhoff equation under initial conditions and boundary conditions, using the previous research results for reference. Firstly, the existence of bounded absorption set is proved by using a prior estimation, then the existence and uniqueness of the global solution of the problem is proved by using the classical Galerkin’s method. Finally, Housdorff dimension and fractal dimension of the family of global attractors are estimated by linear variational method and generalized Sobolev-Lieb-Thirring inequality.

Long Time Behavior of a Class of Generalized Beam-Kirchhoff Equations

In this paper, we study the long time behavior of a class of generalized Beam-Kirchhoff equation , and prove the existence and uniqueness of the global solution of this class of equation by Galerkin method by making some assumptions about the nonlinear function term . The existence of the family of global attractor and its Hausdorff dimension and Fractal dimension estimation are proved.

Complex Maxwell’s Equations

Maxwell’s equations in electromagnetism can be categorized into three dis-tinct groups based on the electromagnetic source when employing quaterni-ons. Each group represents a self-contained system in which Maxwell’s equations are applied and validated concurrently, in contrast to the previous approach that did not account for this. It has been noted that the formulation of these Maxwell equations ultimately results in the formulation of Max-well’s equations utilizing the scalar function.

基于断裂力学的腐蚀钢丝疲劳S-N曲线研究

高强冷拔钢丝因其优异的力学性能被广泛应用于缆索承重结构,服役过程中钢丝腐蚀与断裂威胁着结构安全。钢丝腐蚀后表面的蚀坑因应力集中效应易产生疲劳裂纹,因此腐蚀钢丝剩余疲劳寿命主要取决于裂纹的扩展阶段,且裂纹扩展受到腐蚀程度(初始损伤)的显著影响。精确表征腐蚀程度并准确评估服役期钢丝腐蚀后的剩余寿命,仍是亟待解决的维修难题。首先采用蚀坑最大深度表征钢丝腐蚀程度,以应力强度因子幅作为裂纹扩展驱动力指标,基于断裂力学的Paris公式建立受拉钢丝裂纹扩展理论模型,进而推导腐蚀钢丝的疲劳S-N曲线。大量既有的腐蚀钢丝疲劳试验结果与理论预测值的比较验证了所提理论模型及疲劳S-N曲线的合理性和有效性。含多个蚀坑的钢丝裂纹扩展有限元模拟表明,等效单裂纹适用于多腐蚀坑钢丝疲劳寿命的简化分析。参数分析表明,蚀坑深度和腐蚀环境是腐蚀钢丝疲劳强度的关键影响因素,而应力比和疲劳加载频率影响较小,所提出的S-N曲线在钢丝常规服役条件下适用于评估其剩余寿命,可得到合理的预测结果。基于断裂力学所提的理论方法和疲劳S-N曲线可为腐蚀钢丝的疲劳强度和寿命评估提供参考。

温度对UNS N08825合金在高含H_(2)S环境中耐腐蚀性能的影响

在不同温度(80,90,110,132℃)、高压、高含H2S和高含氯离子环境中对UNS N08825合金进行了72 h浸泡和电化学腐蚀试验,研究了温度对合金耐腐蚀性能的影响,分析了其影响机理。结果表明:在80℃下试验合金几乎不发生腐蚀,90,110,132℃下合金表面出现黑色腐蚀产物,并且110,132℃下的腐蚀产物增多且呈疏松多孔特征;随着温度升高,合金表面的点蚀坑数量增多且深度增加,最大点蚀速率和均匀腐蚀速率均增大;随着温度升高,试验合金的自腐蚀电位、电荷转移电阻和钝化膜电阻减小,自腐蚀电流密度增大。当温度低于90℃时,合金表面钝化膜均匀致密,点蚀敏感性低,具有较好的耐腐蚀性能;当温度不低于90℃时,元素硫的水解加剧,合金表面钝化膜发生破坏,点蚀更易发生,耐腐蚀性能变差。

采用两步炭化法和熔盐模板法制备N、S共掺杂煤基硬炭及共储钠性能

硬炭因资源丰富、结构稳定及安全性高等优势,已成为钠离子电池常用阳极材料。其中,煤基衍生硬炭受到了广泛的关注。本工作以长焰煤为碳源,硫脲为氮硫源,NaCl为模板,通过两步炭化工艺和杂原子掺杂相结合的方法合成了N和S共掺杂的煤基硬炭(NSPC1200)。两步炭化过程在调节碳微晶结构和扩大层间距方面发挥了重要的作用。N和S的共掺杂调节了炭材料的电子结构,赋予其更多的活性位点;此外,引入NaCl作为模板有助于孔结构的构建,有利于电极和电解质之间的接触,从而实现Na+和电子的有效传输。在协同作用下,样品NSPC1200表现出优异的储钠能力,在20 mA g^(−1)电流密度下呈现314.2 mAh g^(−1)的可逆容量。即使在100 mA g^(−1)下循环200次,仍保持224.4 mAh g^(−1)的比容量。这项工作成功实现了策略性调整煤基炭材料微观结构的目标,最终获得了具有优异的电化学性能的硬炭阳极。

N,S-CQDs/Mn_(0.5)Cd_(0.5)S的制备及光催化产氢性能

采用水热法制备了氮、硫共掺杂碳量子点(N,S-CQDs),并用超声辅助法合成了N,S-CQDs修饰的Mn_(0.5)Cd_(0.5)S(MCS)复合光催化材料.研究表明,在模拟太阳光照射下,质量分数为1%的N,S-CQDs/MCS的光催化产氢速率显著提高,为34159.25μmol/(g·h),是MCS的1.63倍,是质量分数为1%的CQDs/MCS的1.14倍.此复合物光催化产氢性能的提高主要归因于N,S-CQDs优异的电子迁移能力以及MCS与N,S-CQDs之间的紧密接触.

Tuning interface mechanism of FeCo alloy embedded N,S-codoped carbon substrate for rechargeable Zn-air battery

The interface mechanism between catalyst and carbon substrate has been the focus of research.In this paper,the FeCo alloy embedded N,S co-doped carbon substrate bifunctional catalyst(FeCo/S-NC)is obtained by a simple one-step pyrolysis strategy.The experimental results and density functional theory(DFT)calculation show that the formation of FeCo alloy is conducive to promoting electron transfer,and the introduction of S atom can enhance the interaction between FeCo alloy and carbon substrate,thus inhibiting the migration and agglomeration of particles on the surface of carbon material.The FeCo/SNC catalysts show outstanding performance for oxygen reduction reaction(ORR)and oxygen evolution reaction(OER).FeCo/S-NC shows a high half-wave potential(E_(1/2)=0.8823 V)for ORR and a low overpotential at 10 mA cm^(-2)(E_(j=10)=299 mV)for OER.In addition,compared with Pt/C+RuO_(2) assembled Zn-air battery(ZAB),the FeCo/S-NC assembled ZAB exhibits a larger power density(198.8 mW cm^(-2)),a higher specific capacity(786.1 mA h g_(zn)~(-1)),and ultra-stable cycle performance.These results confirm that the optimized composition and the interfacial interaction between catalyst and carbon substrate synergistically enhance the electrochemical performance.

基于N掺杂Ti_(3)C_(2)MXene量子点的荧光探针用于Hg2+和S2-的传感检测

基于N掺杂Ti_(3)C_(2) MXene量子点(N-Ti_(3)C_(2) MQDs)荧光探针和配位相互作用,构建了一种检测Hg^(2+)和S^(2-)的“开-关-开”型荧光传感新方法.研究发现,制备的N-Ti_(3)C_(2) MQDs发射蓝色荧光(λem=440 nm),荧光量子产率为15.7%.Hg^(2+)与N-Ti_(3)C_(2) MQDs表面的—NH2,—COOH,—OH等官能团产生选择性配位作用,导致N-Ti_(3)C_(2) MQDs体系荧光猝灭.当加入S^(2-)后,由于S^(2-)与Hg^(2+)之间强的结合力,形成HgS沉淀,从而使N-Ti_(3)C_(2) MQDs体系荧光恢复.基于该原理,构建了一种“开-关-开”型荧光传感方法,实现了对Hg^(2+)和S^(2-)的定量检测.N-Ti_(3)C_(2) MQDs探针的荧光强度与Hg^(2+)浓度在0.02~200μmol/L范围内呈良好线性关系,检出限为10 nmol/L(S/N=3);与S^(2-)浓度在0.07~150μmol/L范围内呈良好线性关系,检出限为30 nmol/L(S/N=3).该方法具有成本低、操作简单、灵敏度高和选择性好等特点,并可用于水样中Hg^(2+)和S^(2-)的检测.