非牛顿磁性纳米流体旋转流动

研究非牛顿磁性纳米磁体定常旋转流动。考虑到磁化强度的特殊贡献,应用修正的非牛顿磁性流体运动方程,磁场作用是纯扩散,不受流体运动速度场的影响。给出了多种流动的解析解。应用计算机符号运算技术,软件Maple16,可以得岀非牛顿上随体Maxwell磁性纳米流体和Oldroyd B磁性纳米流体非定常流动方程的各级近似解的常微分方程,获得解析和数值解。在研究中对纳米流体的物质常数作纳米修正。     

利用Tanh方法产生的(2+1)和(3+1)维非线性可积方程的精确解

主要利用Tanh函数方法,对两个高维五阶非线性可积方程进行了讨论,通过行波约化,分别将(2+1)和(3+1)维非线性可积方程转化为常微分方程.结合Riccati方程的性质,分别得到关于若干参变量的代数系统,借助于Mathematica软件符号运算功能,最终得到了上述两个高维方程的精确解.     

Oldroyd B流体依时性管内流动的变分解析方法

在本文中，研究上随体Ｏｌｄｒｏｙｄ Ｂ流体在水平管内依时性流动，该问题可归结为无量纲速度分量三阶偏微分方程的初边值问题，采用改进的Ｋａｎｔｏｒｏｖｉｃｈ方法，将该方程化为各级近似的二阶常微分方程组的初值问题，通过Ｌａｐｌａｃｅ变换，求得其二阶常微分方程的解析解。在本文中，提出了变分解析的新概念，获得了二级近似变分解析解，其中包括常压力力梯度和周期性压力梯度两种情形，应用计算机符呈处理和Ｌａｐｌａｃ     

随机分数阶微分方程初值问题基于模拟方程法的数值求解

基于模拟方程法,提出了一种求解随机分数阶微分方程初值问题的数值方法.考虑含两个分数阶导数项的微分方程,引入两个线性的、非耦合的随机模拟方程,利用它们解构原方程,借助Laplace变换及逆变换,得到方程解的积分表达式,同时建立起两个模拟方程之间的联系,结合初始状态,得到求解随机微分方程初值问题的数值迭代算法.作为特例,对于含两个分数阶导数项线性常微分方程的初值问题,给出了基于模拟方程法的数值解法的显式结果.该方法是稳定的,它的误差仅存在于积分近似时的截断误差和计算软件的舍入误差.应用实例说明了数值方法在确定和随机情形的有效性和准确性.     

常微分方程近似解的LS-SVM改进求法

提出一种基于最小二乘支持向量机(LS-SVMs)的求解常微分方程近似解的改进方法.该方法首先通过离散计算域,将常微分方程转换为有约束条件的目标优化问题,然后利用径向基(RBF)核函数可导的性质,用带有导数形式的LS-SVM模型将此优化问题转化为LS-SVM回归问题,进而进行求解.最终得到的闭式近似解具有精度高、连续可微、结构简单且形式固定的特点.该方法适用于任意阶非刚性和奇异的线性常微分方程初值问题和边值问题,以及一阶非线性常微分方程问题.仿真结果验证了该方法具有良好的有效性.     

输入带有时滞的线性系统的镇定

考虑带有输入时滞的线性系统的镇定问题.通过把时滞写成一阶传播方程,带有输入时滞的镇定问题转化为常微分方程和一阶双曲方程组成的串联系统的镇定问题.与现有Backstepping方法不同,文章给出了新的变换,其核函数是一阶倒向向量值常微分方程,这使得控制的设计更加简单.文章给出了新的状态反馈控制器,并证明了闭环系统解的适定性和指数稳定性.数值模拟说明,给出的方法是非常有效的.     

三维热传导方程的高精度有限差分方法

提出了数值求解三维热传导方程的一个四阶精度的有限差分格式,首先对三个空间方向上的二阶导数项,采用四次样条函数来近似,从而得到半离散的常微分方程.然后利用常微分方程的解析解表达式,时间矩阵利用Padé近似,得到时间和空间均为四阶精度的差分格式.最后利用方法计算了两个数值算例,并与文献中结果进行了对比,从而验证了高精度格式的性能.     

不可压流体的边界层问题

研究三维有界区域在边界上有流动的不可压流体的边界层问题,导出了Navier-Stokes方程区域内部的近似方程(Euler方程和线性化的Euler方程)和边界附近近似的方程(零阶边界层方程与一阶边界层方程),证明了这种近似的合理性.     

分数阶常微分方程迭代方法的解

通过采用分数阶积分与导数的复合,把分数阶常微分方程转化为积分方程.构造出迭代格式,证明它的收敛性,进一步给出近似解的误差估计.并给出数值例子.     

特殊坐标系中特殊非牛顿流边界层方程的相似解

给出了在一个特殊坐标系中三阶流体的二维定常运动方程组．该坐标系中由无粘流体的势流确定，即以环绕任意物体的非粘性流动的流线为Ф-坐标，速度势线为ψ-坐标，构成正交曲线坐标系．结果表明，边界层方程与浸没在流体中的物体的形状无关．第一次近似假定第二梯度项与粘性项和第三梯度项相比，可以忽略不计．第二梯度项的存在，将防碍第三梯度流相似解的比例变换的导出．利用李群方法计算了边界层方程的无穷小生成元．将边界层方程组变换为常微分方程组．利用Runge-Kutta法结合打靶技术求解了该非线性微分方程组的数值解．     

冠状动脉狭窄情况下的非牛顿血液流动和大分子传质

针对冠状动脉狭窄的情况,采用数值模拟方法求解了牛顿流体与非牛顿流体(幂次律流体和Casson流体)的定常与脉动的流场。在此基础上,求解了LDL(低密度脂肪蛋白)和Albumin(血清白蛋白)的浓度场。根据计算结果,详细讨论了壁面剪应力、非牛顿流效应、分子大小等因素对大分子传质的影响;并对牛顿流体与非牛顿流体、定常流动与脉动流动的大分子浓度场进行了比较,这些结果对于了解动脉硬化成因与流动特性和大分子传质的联系提供了较为丰富的信息。     

精确解析法的子结构算法

文[1]提出精确解析法,用以求解任意变系数常微分方程,并利用初参数算法给出一个解的解析表达式.但利用初参数算法,对某一类问题,如长柱壳弯曲和振动等,它们的解将难以在计算机上得到.本文通过非均匀轴对称长圆柱壳弯曲问题,给出精确解析法的子结构算法,它能够计算初参数算法在计算机上不能解决的问题.问题最后和初参数算法一样能归结为求解一个低阶代数方程组.文末给出算例,表明本文算法的正确性,并和初参数算法作了比较.     

Steady flow of incompressible fluid between two co-rotating disks

This article provides an analytical solution of the Navier–Stokes equations for the case of the steady flow of an incompressible fluid between two uniformly co-rotating disks. The solution is derived from the asymptotical evolution of unknown components of velocity and pressure in a radial direction – in contrast to the Briter–Pohlhausen analytical solution, which is supported by simplified Navier–Stokes equations. The obtained infinite system of ordinary differential equations forms recurrent relations from which unknown functions can be calculated successively. The first and second approximations of solution are solved analytically and the third and fourth approximations of solutions are solved numerically. The numerical example demonstrates agreements with results obtained by other authors using different methods.     

Computational treatment of free convection effects on perfectly conducting viscoelastic fluid

We introduced a magnetohydrodynamic model of boundary-layer equations for a perfectly conducting viscoelastic fluid. This model is applied to study the effects of free convection currents with one relaxation time on the flow of a perfectly conducting viscoelastic fluid through a porous medium, which is bounded by a vertical plane surface. The state space approach is adopted for the solution of one-dimensional problems. The resulting formulation together with the Laplace transform technique is applied to a thermal shock problem and a problem for the flow between two parallel fixed plates, both without heat sources. Also a problem for the semi-infinite space in the presence of heat sources is considered. A discussion of the effects of cooling and heating on a perfectly conducting viscoelastic fluid is given. Numerical results are illustrated graphically for each problem considered.     

精细积分区段混合能矩阵的一种计算方法

根据结构力学与最优控制的模拟理论中阐述的各混合能矩阵的力学意义,介绍了一种利用微分方程组的状态转移矩阵计算区段混合能矩阵的方法,其计算结果与泰勒级数展开法是一致的。     

The impact of finite precision arithmetic and sensitivity on the numerical solution of partial differential equations

In this paper, we address some fundamental issues concerning “time marching” numerical schemes for computing steady state solutions of boundary value problems for nonlinear partial differential equations. Simple examples are used to illustrate that even theoretically convergent schemes can produce numerical steady state solutions that do not correspond to steady state solutions of the boundary value problem. This phenomenon must be considered in any computational study of nonunique solutions to partial differential equations that govern physical systems such as fluid flows. In particular, numerical calculations have been used to “suggest” that certain Euler equations do not have a unique solution. For Burgers' equation on a finite spatial interval with Neumann boundary conditions the only steady state solutions are constant (in space) functions. Moreover, according to recent theoretical results, for any initial condition the corresponding solution to Burgers' equation must converge to a constant as t → ∞. However, we present a convergent finite difference scheme that produces false nonconstant numerical steady state “solutions.” These erroneous solutions arise out of the necessary finite floating point arithmetic inherent in every digital computer. We suggest the resulting numerical steady state solution may be viewed as a solution to a “nearby” boundary value problem with high sensitivity to changes in the boundary conditions. Finally, we close with some comments on the relevance of this paper to some recent “numerical based proofs” of the existence of nonunique solutions to Euler equations and to aerodynamic design.     

A non-standard symmetry-preserving method to compute bounded solutions of a generalized Newell-Whitehead-Segel equation

In this work, we propose a finite-difference scheme to approximate the solutions of a generalization of the classical, one-dimensional, Newell-Whitehead-Segel equation from fluid mechanics, which is an equation for which the existence of bounded solutions is a well-known fact. The numerical method preserves the skew-symmetry of the problem of interest, and it is a non-standard technique which consistently approximates the solutions of the equation under investigation, with a consistency of the first order in time and of the second order in space. We prove that, under relatively flexible conditions on the computational parameters of the method, our technique yields bounded numerical approximations for every set of bounded initial estimates. Some simulations are provided in order to verify the validity of our analytical results. In turn, the validity of the computational constraints under which the method guarantees the preservation of the boundedness of the approximations, is successfully tested by means of computational experiments in some particular instances.     

On the identification of coefficients of elliptic problems by asymptotic regularization

The problem of recovering coefficients of elliptic problems from measured data is considered. An algorithm is developed to identify the unknown coefficients without a minimization technique. The method is based on the construction of certain time-dependent problems which contain the original equation as asymptotic steady state. A Liapunovtype a-priori estimate is fundamental to prove that the solution of the time-dependent regularized equations approach a solution of the original problem as t →∞. A related behavior is proved for the solution of corresponding finite-dimensional Galerkin approximations. A stability result is proved for the Galerkin approximations.