A Novel Method for Linear Systems of Fractional Ordinary Differential Equations with Applications to Time-Fractional PDEs

This paper presents an efficient numerical technique for solving multi-term linear systems of fractional ordinary differential equations(FODEs)which have been widely used in modeling various phenomena in engineering and science.An approximate solution of the system is sought in the formof the finite series over the Müntz polynomials.By using the collocation procedure in the time interval,one gets the linear algebraic system for the coefficient of the expansion which can be easily solved numerically by a standard procedure.This technique also serves as the basis for solving the time-fractional partial differential equations(PDEs).The modified radial basis functions are used for spatial approximation of the solution.The collocation in the solution domain transforms the equation into a system of fractional ordinary differential equations similar to the one mentioned above.Several examples have verified the performance of the proposed novel technique with high accuracy and efficiency.

基于高维PDE投影恢复的低剂量CT重建方法

目的提出一种基于高维偏微分方程(PDE)投影恢复的低剂量CT重建方法。方法先将原始的投影数据映射到高维空间中,构造投影数据的高维表示,通过移动高维空间中的点来对高维表示进行更新,再使用偏微分方程对投影数据进行滤波,最后将恢复后的数据使用FBP算法重建出最终CT图像。结果在Shepp-Logan体模实验中,与FBP,PWLS-QM和TGV-WLS方法相比,新方法在相对均方根误差指标上分别降低了68.87%、50.15%和27.36%,结构相似性上分别提高了23.50%,8.83%和1.62%,特征相似性上分别提高了17.30%、2.71%和2.82%。在腹部临床数据实验中,与FBP,PWLS-QM和TGV-WLS方法相比,新方法在相对均方根误差中分别降低了42.09%、31.04%和21.93%,结构相似性上分别提高了18.33%、13.45%和4.63%,特征相似性上分别提高了3.13%、1.46%和1.10%。结论本研究提出的新方法在有效去除低剂量CT图像中的条形伪影和噪声的同时,可以保持图像的空间分辨率。

基于自适应神经网络的PDEs求解研究

针对当前基于神经网络的PDEs求解方法效率和精度均不够理想的缺陷,研究提出一种基于改进BP神经网络(BP neural network,BPNN)的PDEs求解模型。首先,参照自适应网格法来改进神经网络结构,构建自适应神经网络,改进模型的输出精度;其次,提出一种引入Levy飞行机制和鲸鱼优化算法(Whale optimization algorithm,WOA)的改进海鸥优化算法(Improved Seagull Optimization Algorithm,ISOA)来优化BPNN,寻找BPNN的最佳参数,提高模型的性能;基于上述内容,构建基于ISOA-BPNNPDEs智能求解模型。结果显示,该模型的F1值为95.74%,准确率达到97.96%,输出误差为0.021,优于当前最先进的两种PDEs求解模型。上述内容表明,研究构建的基于ISOA-BPNNPDEs智能求解模型能够高效、准确地实现PDEs求解,为PDEs求解研究提供了新的路径。

Some Modified Equations of the Sine-Hilbert Type

Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.

基于Matlab PDE工具箱的土体冻结温度场模拟研究

分析并运用Matlab PDE tool模拟模型槽土体冻结的温度场并验证模拟值的准确性。首先,通过合理的边界假设,使三维的模型槽土体冻结热传导问题转化为二维空间;然后,运用Matlab PDE tool进行模拟并对比测温点的模拟温度和实测温度。结果表明:土体冻结中温度场导热可以通过Matlab PDE tool模拟,且操作简单易于掌握;测温点模拟值的准确性与测温点距边界距离和冻结时间有关。

Identification of partial differential equations from noisy data with integrated knowledge discovery and embedding using evolutionary neural networks

Identification of underlying partial differential equations(PDEs)for complex systems remains a formidable challenge.In the present study,a robust PDE identification method is proposed,demonstrating the ability to extract accurate governing equations under noisy conditions without prior knowledge.Specifically,the proposed method combines gene expression programming,one type of evolutionary algorithm capable of generating unseen terms based solely on basic operators and functional terms,with symbolic regression neural networks.These networks are designed to represent explicit functional expressions and optimize them with data gradients.In particular,the specifically designed neural networks can be easily transformed to physical constraints for the training data,embedding the discovered PDEs to further optimize the metadata used for iterative PDE identification.The proposed method has been tested in four canonical PDE cases,validating its effectiveness without preliminary information and confirming its suitability for practical applications across various noise levels.

Data-Driven Ai-and Bi-Soliton of the Cylindrical Korteweg-de Vries Equation via Prior-Information Physics-Informed Neural Networks

By the modifying loss function MSE and training area of physics-informed neural networks(PINNs),we propose a neural networks model,namely prior-information PINNs(PIPINNs).We demonstrate the advantages of PIPINNs by simulating Ai-and Bi-soliton solutions of the cylindrical Korteweg-de Vries(cKdV)equation.

PDE Standardization Analysis and Solution of TypicalMechanics Problems

A numerical approach is an effective means of solving boundary value problems(BVPs).This study focuses on physical problems with general partial differential equations(PDEs).It investigates the solution approach through the standard forms of the PDE module in COMSOL.Two typical mechanics problems are exemplified:The deflection of a thin plate,which can be addressed with the dedicated finite element module,and the stress of a pure bending beamthat cannot be tackled.The procedure for the two problems regarding the three standard forms required by the PDE module is detailed.The results were in good agreement with the literature,indicating that the PDE module provides a promising means to solve complex PDEs,especially for those a dedicated finite element module has yet to be developed.

Matrix Riccati Equations in Optimal Control

In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.

THE STABILITY OF BOUSSINESQ EQUATIONS WITH PARTIAL DISSIPATION AROUND THE HYDROSTATIC BALANCE

This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).

The Maxwell-Heaviside Equations Explained by the Theory of Informatons

In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.

THE SMOOTHING EFFECT IN SHARP GEVREY SPACE FOR THE SPATIALLY HOMOGENEOUS NON-CUTOFF BOLTZMANN EQUATIONS WITH A HARDPOTENTIAL

In this article, we study the smoothing effect of the Cauchy problem for the spatially homogeneous non-cutoff Boltzmann equation for hard potentials. It has long been suspected that the non-cutoff Boltzmann equation enjoys similar regularity properties as to whose of the fractional heat equation. We prove that any solution with mild regularity will become smooth in Gevrey class at positive time, with a sharp Gevrey index, depending on the angular singularity. Our proof relies on the elementary L^(2) weighted estimates.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

Equivalence between the internal observability and equation

This paper is concerned with a third order in time linear Moore-Gibson-Thompson equation which describes the acoustic velocity potential in ultrasound wave program.Influenced by the work of Kaltenbacher,Lasiecka and Marchand(Control Cybernet.2011,40:971-988),we establish an observability inequality of the conservative problem,and then discuss the equivalence between the exponential stabilization of a dissipative system and the internal observational inequality of the corresponding conservative system.

Multi-soliton solutions, breather-like and bound-state solitons for complex modified Korteweg–de Vries equation in optical fibers

Under investigation in this paper is a complex modified Korteweg–de Vries(KdV) equation, which describes the propagation of short pulses in optical fibers. Bilinear forms and multi-soliton solutions are obtained through the Hirota method and symbolic computation. Breather-like and bound-state solitons are constructed in which the signs of the imaginary parts of the complex wave numbers and the initial separations of the two parallel solitons are important factors for the interaction patterns. The periodic structures and position-induced phase shift of some solutions are introduced.

ON THE STABILITY OF PERIODIC SOLUTIONS OF PIECEWISE SMOOTH PERIODIC DIFFERENTIAL EQUATIONS

In this paper,we address the stability of periodic solutions of piecewise smooth periodic differential equations.By studying the Poincarémap,we give a sufficient condition to judge the stability of a periodic solution.We also present examples of some applications.

Exact and heuristic formulae to compute the geodetic height from the ellipse equation

The conversion of the cartesian coordinates of a point to its geodetic equivalent coordinates in reference to the geodetic ellipsoid is one of the main challenges in geodesy.The ellipse equation in the meridian plane significantly influences the value of the geodetic coordinates.This research analyzes this influence and how it can contribute to their solutions.The study investigates the mathematical relation between them and presents an exact formula relating to the geodetic height and the ellipse equation.In addition,a heuristic formula for the relation between the geodetic height and the ellipse equation is proposed,which is independent of the geodetic latitude and has a relative accuracy better than 99.9 %.The calculation is stable,and the cost is low.

A deep learning method based on prior knowledge with dual training for solving FPK equation

The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.

Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

On Two Types of Stability of Solutions to a Pair of Damped Coupled Nonlinear Evolution Equations

The stability of a set of spatially constant plane wave solutions to a pair of damped coupled nonlinear Schrödinger evolution equations is considered. The equations could model physical phenomena arising in fluid dynamics, fibre optics or electron plasmas. The main result is that any small perturbation to the solution remains small for all time. Here small is interpreted as being both in the supremum sense and the square integrable sense.