FracSym: Automated symbolic computation of Lie symmetries of fractional differential equations

In this paper, we present an algorithm for the systematic calculation of Lie point symmetries for fractional order differential equations (FDEs) using the method as described by Buckwar & Luchko (1998) and Gazizov, Kasatkin & Lukashchuk (2007, 2009, 2011). The method has been generalised here to allow for the determination of symmetries for FDEs with n$n$ independent variables and for systems of partial FDEs. The algorithm has been implemented in the new MAPLE package FracSym (Jefferson and Carminati 2013) which uses routines from the MAPLE symmetry packages DESOLVII (Vu, Jefferson and Carminati, 2012) and ASP (Jefferson and Carminati, 2013). We introduce FracSym by investigating the symmetries of a number of FDEs; specific forms of any arbitrary functions, which may extend the symmetry algebras, are also determined. For each of the FDEs discussed, selected invariant solutions are then presented.     

A conservative discontinuous Galerkin scheme for the 2D incompressible Navier–Stokes equations

In this paper we consider a conservative discretization of the two-dimensional incompressible Navier–Stokes equations. We propose an extension of Arakawa’s classical finite difference scheme for fluid flow in the vorticity–stream function formulation to a high order discontinuous Galerkin approximation. In addition, we show numerical simulations that demonstrate the accuracy of the scheme and verify the conservation properties, which are essential for long time integration. Furthermore, we discuss the massively parallel implementation on graphic processing units.     

Jet bundle formulation of infinite-dimensional port-Hamiltonian systems using differential operators

We consider infinite-dimensional port-Hamiltonian systems described on jet bundles. Based on a power balance relation we introduce the port-Hamiltonian system representation using differential operators regarding the structural mapping, the dissipation mapping and the input mapping. In contrast to the well-known representation on the basis of the underlying Stokes–Dirac structure our approach is not necessarily based on using energy-variables which leads to a different port-Hamiltonian representation of the analyzed partial differential equations. The presented constructions will be specialized to mechanical systems to which class also the presented examples belong.     

Modelling the dynamics of the students’ academic performance in the German region of the North Rhine-Westphalia: an epidemiological approach with uncertainty

Student academic underachievement is a concern of paramount importance in Europe, where around 15% of the students in the last high school courses do not achieve the minimum knowledge academic requirement. In this paper, we propose a model based on a system of differential equations to study the dynamics of the students’ academic performance in the German region of the North Rhine-Westphalia. This approach is supported by the idea that both good and bad study habits, are a mixture of personal decisions and influence of classmates. This model allows us to forecast the student academic performance by means of confidence intervals over the next few years.     

PHASE-SPACE DIFFUSION EQUATIONS FOR SINGLE BROWNIAN PARTICLE MOTION NEAR SURFACES

In numerous physical processes involving the motion of micron and submicron sized particles near surfaces, such as the filtration of hydrosols and aerosols, the particle motion is the net result of the combined effects of fluid convection, external forces, particle inertia, Brownian particle motion, and particle-surface fluid dynamic interactions. The most general method of describing particle motion under the combined action of these effects is through the so-called Fokker-Planck equation. In the absence of particle-surface fluid dynamic interactions, the Fokker-Planck equation is well-known, and it has been applied in a general way to problems involving the adsorption or deposition of Brownian particles onto surfaces through a solution technique known as the Brownian dynamics simulation method.

In this study, the Fokker-Planck equation for Brownian particle motion near surfaces is generalized to include particle-surface fluid dynamic interactions. The Fokker-Planck equation is shown to follow from the Liouville equation for the Brownian particle and n-fluid molecules present in the system, thus, establishing a firm theoretical foundation for the Fokker-Planck equation and the various other phase-space diffusion equations that follow from it.

Based on diagonalization of the Fokker-Planck equation, its short-time behavior is also derived here which enables a generalization of the Brownian dynamics method for the study of particle motion near surfaces including fluid dynamic interactions. Additionally, a perturbation solution of the Fokker-Planck equation under the conditions of small, but finite particle Stokes number is also derived. These solutions are shown to agree with previously given representations of the Smoluchowski or convective-diffusion equation for Brownian particle motion near surfaces, as well as with inertial corrections to the Smoluchowski equation available in the literature. This latter equation is also generalized here to include particle-surface fluid dynamic interactions.     

改进型Boussinesq方程高精度紧致差分显格式

采用一种高精度的紧致差分显格式对改进型Boussinesq方程进行数值求解;采用具有TVD性质的三阶Runge-Kutta方法进行预报,用三次样条函数进行校正,时间精度可达到四阶;在空间离散上采用六阶精度的三点紧致显格式进行计算;运用以上数值格式对Beji和Nadaoka改进型Boussinesq方程进行了求解,求解证明:高精度的数值结果和已知的试验结果吻合良好.作为验证算例,同时对波浪在台阶上的传播进行了模拟,从效果对比上可以看出,所得结果明显比Kittitanasuan的计算结果更靠近试验值.     

电路模拟中非线性方程的周期波形响应

用范数估计方法对非线性高阶微分方程的周期边值问题进行了讨论,通过对非线性二阶微分方程周期边值问题的详细讨论,给出了系统函数对某些变量偏导数的某种范数小于1时,非线性二阶微分方程的波形松弛算法产生的迭代序列收敛到该方程的周期解.用类似的方法给出了非线性高阶微分方程的波形松弛算法产生的迭代序列收敛到该方程周期解的充分性条件.     

利用Liapanov直接方法判定差分方程稳定性

研究差分方程稳定性的强有力的方法是Liapanov直接方法,这一方法的核心思想是针对所研究的方程,构造出其相应的Liapanov函数,利用它来判定方程的稳定性．     

高维力学系统的能量方程

在低维系统中,能量方程因其明晰的物理意义而得到广泛的应用,而研究高维力学系统的能量方程也同样具有理论价值和实际意义,如位移对时间的导数是速度,速度对时间的导数是加速度,这些具有明显的物理意义,而加速度的导数其物理意义就不是很清晰了,但具有理论上的意义.本文应用广义经典力学中关于广义拉格朗日函数、广义动量和广义哈密顿函数等概念,推导了高维系统的能量方程,文中举了实例具体说明新方程的应用,为力学数学系统能量方程的推广提供了一种途经.     

基于物理的波浪实时模拟方法

为解决广域范围内波浪状态的实时计算和可视化问题,结合无粘滞流体力学的物理模型,提出了一种以三角形为控制单元的有限体积简化算法。该算法的优点在于：以不规则边界区域的非规则三角网格为基础,通过简化通量向量分裂方法获得三角形控制单元的边界数值通量,能快速逼近二维浅水方程的解进而模拟非规则边界浅水的实时流动。实验结果显示,所提方法能在符合现实世界物理规律的前提下较好地实现大规模波浪的实时可视化模拟。