饱和多孔弹性杆流固耦合动力响应的保结构算法

研究了不可压饱和多孔弹性杆的流固耦合动力响应问题.基于多孔介质理论,根据多孔介质流固混合物动量方程、孔隙流体动量方程及体积分数方程,建立流固耦合不可压饱和多孔弹性杆的轴向振动方程;引入正则变量,构造饱和多孔弹性杆轴向振动方程的广义多辛保结构形式、广义多辛守恒律及广义多辛局部动量误差;采用中点Box离散方法得到轴向振动方程的广义多辛离散格式、广义多辛守恒律数值误差及局部动量数值误差;数值模拟不可压饱和多孔弹性杆的轴向振动过程及流相渗流速度分布,考察了流固两相耦合系数对轴向振动过程及广义多辛守恒律误差和局部动量误差的影响.结果表明,已构造的广义多辛保结构算法具有很高的精确性和长时间的数值稳定性.     

Landau-Ginzburg-Higgs方程的多辛Runge-Kutta方法

非线性波动方程作为一类重要的数学物理方程吸引着众多的研究者,基于Hamilton空间体系的多辛理论研究了Landau-Ginzburg-Higgs方程的多辛算法,讨论了利用Runge-Kutta方法构造离散多辛格式的途径,并构造了一种典型的半隐式的多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

膜自由振动的多辛方法

基于Hamilton空间体系的多辛理论研究了膜自由振动问题,讨论了构造复合离散多辛格式的方法,并构造了一种典型的9×3点半隐式的多辛复合离散格式,该格式满足多辛守恒律、能量守恒律和动量守恒律．数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性．     

广义Boussinesq方程的多辛方法

广义Boussinesq方程作为一类重要的非线性方程有着许多有趣的性质,基于Hamilton空间体系的多辛理论研究了广义Boussinesq方程的数值解法,构造了一种等价于多辛Box格式的新隐式多辛格式,该格式满足多辛守恒律、局部能量守恒律和局部动量守恒律．对广义Boussinesq方程孤子解的数值模拟结果表明,该多辛离散格式具有较好的长时间数值稳定性．     

饱和多孔弹性Timoshenko悬臂梁的动、静力弯曲

在经典单相Timoshenko梁变形和孔隙流体仅沿饱和多孔弹性梁轴向运动的假定下,基于不可压饱和多孔介质的三维Gurtin型变分原理,首先建立了饱和多孔弹性Timoshenko悬臂梁动力响应的一维数学模型．在若干特殊情形下,该模型可分别退化为饱和多孔弹性梁的Euler-Bernoulli模型、Rayleigh模型和Shear模型等．其次,利用Laplace变换,分析了固定端不可渗透、自由端可渗透的饱和多孔弹性Timoshenko悬臂梁在自由端阶梯载荷作用下的动静力响应,给出了梁自由端处挠度随时间的响应曲线,考察了固相与流相相互作用系数、梁长细比等参数对悬臂梁动静力行为的影响．结果表明:饱和多孔弹性梁的拟静态挠度具有与粘弹性梁挠度类似的蠕变特征．在动力响应中,随着梁长细比的增大,自由端挠度的振动周期和幅值增大,且趋于稳态值的时间增长,而随着两相相互作用系数的增大,梁挠度振动衰减加快,并最终趋于经典单相弹性Timoshenko梁的静态挠度．     

一类强色散DGH方程的多辛普雷斯曼格式

DGH方程作为一类重要的非线性水波方程有着许多广泛的应用前景.基于Hamilton系统的多辛理论研究了一类强色散DGH方程的数值解法,利用多辛普雷斯曼方法构造了一种典型的半隐式的多辛格式.分析了该格式的局部能量和动量守恒律误差,并给出了数值算例.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.     

偏心冲击荷载作用下薄圆板动力学响应的保结构分析

关注动力学系统的局部几何性质,采用多辛分析方法研究了偏心冲击荷载作用下薄圆板振动特性.在探索偏心冲击荷载作用下薄圆板振动问题动力学控制方程的对称性和守恒律的对应关系基础上,对动力学控制方程在多辛体系下重新描述,并采用显式中点差分离散方法构造其多辛格式,通过对存在不同相对偏心距冲击荷载作用下的薄圆板振动过程的数值模拟,研究了相对偏心距对薄圆板振动特性的影响,同时,数值模拟结果也充分体现了多辛算法的良好保结构性能.该研究结果不仅为由于荷载作用位置误差带来的动力学响应偏差估计提供了依据,而且为偏心冲击动力学问题的研究提供了新的途径.     

二维阻尼非线性sine-Gordon方程的共形多辛Fourier拟谱格式

为二维阻尼非线性sine-Gordon方程构造了一个新的共形多辛Fourier拟谱格式.基于原系统的共形多辛哈密尔顿形式,首先在时间和空间方向上分别用辛中点和Fourier拟谱方法进行离散,得到一个全离散格式.随后证明了构造的格式保持离散的共形多辛守恒律.最后数值实验验证了格式的有效性.     

二维非线性Schr?dinger方程的两类局部守恒算法

本文针对二维非线性Schr?dinger方程,提出两类局部守恒算法.不需要考虑边界条件,即可保持任意时空区域上相应的局部能量守恒律和局部动量守恒律.在合适的边界条件下,它们能自然地保持电荷、全局能量或全局动量守恒律.本文同时对算法进行了守恒分析和误差分析.在数值实验部分,本文构造了类似的多辛Preissman算法进行比较,数值结果验证了其长时间计算的优势.     

Sine-Gordon方程的多辛Leap-frog格式

非线性发展方程由于具有多种形式的解析解而吸引着众多的研究者,借助多辛保结构理论研究了Sine-Gordon方程的多辛算法．利用Hamilton变分原理,构造出了sine-Gordon方程的多辛格式;采用显辛离散方法得到了Leap-frog多辛离散格式,该格式满足多辛守恒律;数值结果表明leap-frog多辛离散格式能够精确地模拟sine-Gordon方程的孤子解和周期解,模拟结果证实了该离散格式具有良好的数值稳定性．     

GSDBM方程的多辛几何和多辛Preissman格式研究（英文）

The multi-symplectic geometry for the GSDBM equation is presented in this paper. The multi-symplectic formulations for the GSDBM equation are presented and the local conservation laws are shown to correspond to certain well-known Hamiltonian functionals. The multi-symplectic discretization of each formulation is exemplified by the multisymplectic Preissmann scheme. The numerical experiments are given, and the results verify the efficiency of the Preissmann scheme.     

Multi-symplectic Geometry and Preissmann Scheme for GSDBM Equation

The multi-symplectic geometry for the GSDBM equation is presented in this paper. The multi-symplectic formulations for the GSDBM equation are presented and the local conservation laws are shown to correspond to certain well-known Hamiltonian func-tionals. The multi-symplectic discretization of each formulation is exemplified by the multi-symplectic Preissmann scheme. The numerical experiments are given, and the results verify the efficiency of the Preissmann scheme.     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

The multi-symplectic Fourier pseudospectral method for solving two-dimensional Hamiltonian PDEs

In this paper, the multi-symplectic Fourier pseudospectral (MSFP) method is generalized to solve two-dimensional Hamiltonian PDEs with periodic boundary conditions. Using the Fourier pseudospectral method in the space of the two-dimensional Hamiltonian PDE (2D-HPDE), the semi-discrete system obtained is proved to have semi-discrete multi-symplectic conservation laws and a global symplecticity conservation law. Then, the implicit midpoint rule is employed for time integration to obtain the MSFP method for the 2D-HPDE. The fully discrete multi-symplectic conservation laws are also obtained. In addition, the proposed method is applied to solve the Zakharov–Kuznetsov (ZK) equation and the Kadomtsev–Petviashvili (KP) equation. Numerical experiments on soliton solutions of the ZK equation and the KP equation show the high accuracy and effectiveness of the proposed method.     

The multi-symplectic Fourier pseudospectral method for solving two-dimensional Hamiltonian PDEs

In this paper, the multi-symplectic Fourier pseudospectral (MSFP) method is generalized to solve two-dimensional Hamiltonian PDEs with periodic boundary conditions. Using the Fourier pseudospectral method in the space of the two-dimensional Hamiltonian PDE (2D-HPDE), the semi-discrete system obtained is proved to have semi-discrete multi-symplectic conservation laws and a global symplecticity conservation law. Then, the implicit midpoint rule is employed for time integration to obtain the MSFP method for the 2D-HPDE. The fully discrete multi-symplectic conservation laws are also obtained. In addition, the proposed method is applied to solve the Zakharov-Kuznetsov (ZK) equation and the Kadomtsev-Petviashvili (KP) equation. Numerical experiments on soliton solutions of the ZK equation and the KP equation show the high accuracy and effectiveness of the proposed method.     

Explicit Multi-symplectic Method for a High Order Wave Equation of KdV Type

In this paper, we consider multi-symplectic Fourier pseudospectral method for a high order integrable equation of KdV type, which describes many important physical phenomena. The multi-symplectic structure are constructed for the equation, and the conservation laws of the continuous equation are presented. The multisymplectic discretization of each formulation is exemplified by the multi-symplectic Fourier pseudospectral scheme. The numerical experiments are given, and the results verify the efficiency of the Fourier pseudospectral method.     

On the well-posedness of incompressible flow in porous media with supercritical diffusion

In this article we study the heat transfer equation with a supercritical diffusion term of an incompressible fluid in porous media governed by Darcy's law. We obtain the global well-posedness for small initial data belonging to critical Besov spaces and the local well-posedness for arbitrary initial data. We further show the pointwise blowup criterion.