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.

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 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.

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.

Adomian Decomposition Method for Solving Time Fractional Burgers Equation Using Maple

In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.

An Extended Riccati Equation Method to Find New Solitary Wave Solutions of the Burgers-Fisher Equation

In this paper, our objective is to explore novel solitary wave solutions of the Burgers-Fisher equation, which characterizes the interplay between diffusion and reaction phenomena. Understanding this equation is crucial for addressing challenges in fluid, chemical kinetics and population dynamics. We tackle this task by employing the Riccati equation and employing various function transformations to solve the Burgers-Fisher equation. By adopting different coefficients in the Riccati equation, we obtain a wide range of exact solutions, many of which have not been previously documented. These abundant solitary wave solutions serve as valuable tools for comprehending the Burgers-Fisher equation and contribute to expanding our knowledge in this field.

基于Burger’s蠕变模型的采矿车行驶海底稀软底质下陷研究

在陆地矿产资源日渐枯竭的今天,深海矿产资源已成为全球各个国家争相开采与利用的焦点,深海采矿车是实现深海矿产资源开采的重要装备。海底稀软底质是一种承载力与抗剪强度极低的特殊底质,在采矿作业中,深海稀软底质的物理力学特性直接影响采矿车行走的稳定性。文章选取Burger’s接触模型作为深海稀软底质的本构模型,对某海域海底稀软原状土开展室内三轴试验,通过PFC3D颗粒流数值模拟实验对比实际三轴试验,对稀软底质的Burger’s蠕变模型进行参数标定,同时依据标定结果改变相应参数,针对5种不同底质条件的工况,建立海底采矿车的数字仿真模型,模拟各工况下采矿车在不同行驶速度时的下陷深度。结果显示,下陷深度会随行驶速度呈非线性变化,在一定范围内随着行驶速度的增大而减少并逐渐趋于稳定。同时结果还表明,该区域海底稀软底质具有更高的黏粒含量(38.1%~48.4%)、含水率(88.13%~137.79%)和压缩性(压缩系数:1.86~3.73 MPa^(-1),压缩模量:1.26~2.13 MPa),具有更低的密度(1.3~1.5 g/cm^(3))和强度特性(贯入阻力:0.19~1.32 N,黏聚力:3.7~6.9 kPa,内摩擦角:2.4°~3.9°),即承载力较低,蠕变性能较强。本研究在宏观上做了一般的探讨,为类似参数的稀软底质下海底采矿车的运行安全控制提供了较好借鉴与依据。

一种风向监督双流神经网络--以一维Burgers方程求解为例

针对一维Burgers方程下单一建模方式难以充分考虑不同阶段风向对系数的影响比重,无法有效获得各节点间的关联信息的问题,本文提出了一种风向监督双流神经网络分别预测上下风向的有限差分系数。同时设计了一种风向判断模块,实现了对预测得到有限差分系数的权重融合。通过风向监督双流神经网络,并结合先验知识对学得的系数分配一定的权重,以突出上下风向对预测结果的不同影响,可以有效实现对不同风向上的点分别进行预测,使得空间结构特征信息挖掘更加充分,从而提高差分系数预测的精度。在比传统数值求解方法网格分辨率粗4~8倍的同时,提高了谷歌团队工作的精度,以此提高了计算的速度。

定常Burgers方程的周期解

本文首先对Burgers方程的周期解关于粘性系数ε作幂级数展开,得到逼近函数之间的递推关系。然后,通过对逼近函数作Fourier级数展开,给出了逼近函数Fourier展开的通项形式,并证明了每个逼近函数的Fourier级数是一致且绝对收敛的。

Auto-Bcklund transformation and exact solutions of Wick-type stochastic Burgers equation

Burgers equation in random environment is studied. In order to give the exact solutions of random Burgers equation, we only consider the Wick-type stochastic Burgers equation which is the perturbation of the Burgers equation with variable coefficients by white noise W（t）=Bt, where Bt is a Brown motion. The auto-Baecklund transformation and stochastic soliton solutions of the Wick-type stochastic Burgers equation are shown by the homogeneous balance and Hermite transform. The generalization of the Wick-type stochastic Burgers equation is also studied.

一类非线性Burgers型问题的预测校正紧差分方法

针对一类非线性Burgers型方程,提出一种预测-校正紧差分方法。首先,对时间一阶导数采用一阶Euler格式,时间积分项运用一阶卷积求积公式进行离散,并以MacCormack方法的两步预测-校正方法处理非线性项;然后采用四阶紧差分离散空间的一阶和二阶导数,构造了Euler预测-校正紧差分全离散格式。最后通过案例验证了所提出算法的有效性。

BBM-Burgers方程的非协调有限元方法的超收敛分析

研究Benjamin-Bona-Mahony-Burgers(BBM-Burgers)方程的非协调EQ_(1)^(rot)元的线性化BDF格式下的超收敛性质。通过使用数学归纳法来处理非线性项,并利用该单元已有的高精度结果及插值后处理技术,得到了在对空间剖分尺度和时间步长无网格比约束的前提下,关于离散H^(1)-模意义下具有O(h^(2)+τ^(2))阶的超逼近和超收敛结果。最后,通过给出数值算例验证了理论分析的正确性。

基于PINN的Burgers方程求解模型

传统的数值求解方法面临维数灾难和效率与精度平衡问题,而基于数据驱动的神经网络求解方法又存在训练量冗余和不可解释性问题。针对此问题,物理信息神经网络(Physical Information Neural Networks,PINNs)关注了训练数据中隐含的物理先验知识,融合了神经网络拟合复杂变量的能力,赋予了传统神经网络所缺乏的物理可解释性。应用该算法模型,提出了一种基于PINN的Burgers方程求解模型,该算法模型在训练中施加物理信息约束,因此能用少量的训练样本学习预测到分布在时空域上的偏微分方程模型。实验结果表明,在1+1维Burgers方程算例下,所提方法相比于经典的机器学习算法能有效捕抓到方程的变化并进行精确模拟,相比于有限差分法,可以大幅度缩短模拟时间。通过对不同的网络参数进行比较实验,所提方法在10%的噪声破坏下能产生合理的识别准确度,网络逼近方程的待定系数误差在0.001以内。

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.

Tiêu Equation Experimentation of Understanding by Energy Transfer Quantum Mechanics

Background: The Tiêu equation has a ground roots approach to the process of Quantum Biology and goes deeper through the incorporation of Quantum Mechanics. The process can be measured in plant, animal, and human usage through a variety of experimental or testing forms. Animal studies were conducted for which, in the first day of the study all the animals consistently gained dramatic weight, even as a toxic substance was introduced as described in the introduction of the paper to harm animal subjects which induced weight loss through toxicity. Tests can be made by incorporating blood report results. Human patients were also observed to show improvement to their health as administration of the substance was introduced to the biological mechanism and plants were initially exposed to the substance to observe results. This is consistent with the Tiêu equation which provides that wave function is created as the introduction of the substance to the biological mechanism which supports Quantum Mechanics. The Tiêu equation demonstrates that Quantum Mechanics moves a particle by temperature producing energy thru the blood-brain barrier for example. Methods: The methods for the Tiêu equation incorporate animal studies to include the substance administered through laboratory standards using Good Laboratory Practices under Title 40 C.F.R. § 158. Human patients were treated with the substance by medical professionals who are experts in their field and have knowledge to the response of patients. Plant applications were acquired for observation and guidance of ongoing experiments of animals’ representative for the biologics mechanism. Results: The animal studies along with patient blood testing results have been an impressive line that has followed the Tiêu equation to consistently show improvement in the introduction of the innovation to biologic mechanisms. The mechanism responds to the substance by producing energy to the mechanism with efficient effect. For plant observations, plant organisms responded, and were seen as showing improvement thru visual observation.

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.

Verification of the Landau Equation and Hardy’s Inequality

We prove the L estimate for the isotropic version of the homogeneous landau problem, which was explored by M. Gualdani and N. Guillen. As shown in a region of the smooth potentials range under values of the interaction exponent (2), a weighted Poincaré inequality is a natural consequence of the traditional weighted Hardy inequality, which in turn implies that the norms of solutions propagate in the L1 space. Now, the L estimate is based on the work of De Giorgi, Nash, and Moser, as well as a few weighted Sobolev inequalities.