基于Bernstein多项式的二阶线性奇异边值问题的数值解

提出了一种用Bernstein多项式来构造一类线性奇异两点边值问题的数值解方法.该方法不需要事先对方程进行非奇异化,且若方程的精确解为多项式时,利用这种方法可得方程的精确解.本文包含一些数值实例,并且与三次样条法的数值计算结果进行了比较,从而说明我们提出方法的可靠性和有效性.     

三次样条插值及其在第一类积分方程求解中的应用

从三次样条插值的定义出发,通过研究第一类积分方程中未知函数的三次样条函数逼近,给出了第一类积分方程的三次样条插值离散化.利用该离散化形式,将第一类积分方程转化成线性方程组的形式.由于第一类积分方程的求解通常是不适定的,进而引起线性方程组的病态性.最后,为克服线性方程组的病态性,通过引入未知函数的多重光滑化约束,得到第一...     

多变量样条元法分析弹性地基板的弯曲,振动与稳定问题

本文应用双三次乘积型二元B样条函数来构造弹性地基板的位移、弯矩和扭矩等多种场函数，由混合变分原理导出多变量样条无法方程．文中，对弹性地基板的弯曲、振动与稳定问题作了分析与计算．由于，本文方法设定了各自独立的场函数，因此，所算得的场未知量如位移、弯矩和扭矩值的精度均比较高。     

可调形三次三角Cardinal插值样条曲线

在三次Cardinal插值样条曲线的基础上,引入了三角函数多项式,得到一组带调形参数的三次三角Cardinal样条基函数,以此构造一种可调形的三次三角Cardinal插值样条曲线.该插值样条可以精确表示直线、圆弧、椭圆以及自由曲线,改变调形参数可以调控插值曲线的形状.该插值样条避免了使用有理形式,其表达式较为简洁,计算量也相对较少,从而为多种线段的构造与处理提供了一种通用与简便的方法.     

差分样条插值与差分样面插值及其应用于板弯曲的计算

给出一种简单的有理分式插值——差分样条插值。其变量个数类似于三次样条,光滑性与逼近阶都比三次样条降低一次,即类似于二次样条。推广到二维情况能应用于任意区域,因此对板弯曲的计算能做到解算变量少,效率较高。     

The Cubic B-Spline Method for a Class of Caputo-Fabrizio Fractional Differential Equations北大核心CSCD

基于分数阶微积分基本定理和三次B样条理论,构造了求解线性Caputo-Fabrizio型分数阶微分方程数值解的三次B样条方法,利用分数阶微积分基本定理将初值问题转化为关于解函数的表达式,再使用三次B样条函数逼近表达式中积分项的被积函数,进而计算了一类Caputo-Fabrizio型分数阶微分方程的数值解.给出了所构造的三次B样条方法的误差估计、收敛性和稳定性的理论证明.数值实验表明,该文数值方法在求解一类Caputo-Fabrizio型分数阶微分方程数值解时具有一定的可行性和有效性,且计算精度和计算效率优于现有的两种数值方法.     

非负严格对角占优三对角矩阵逆元素的估计

估计非负三对角矩阵的逆矩阵元素的值,在样条插值的凸性研究方面十分重要.文献[1]给出三对角阵 A=diag(α_i,ρ_i,β_i)在 x_i+β_i=1条件下的估计式.但在样条插值的连续性方程中有许多三对角阵只是非负严格对角占优的,如象在归范样条曲线和局部张力样条曲线插值的情形.因此本文给出非负严格对角占优三对角阵逆元素新的一类估计式.     

单调三次样条的一个充分条件

在理论和实践中,人们对保形拟合问题特别有兴趣、[1]与[2]讨论了二次样条的保形问题,[3]则给出有关三次样条保凸问题的一系列充分条件。 给定点列{(x_i),y_i)}(i=0,1,…,N)(也称型值).通过此点列的三次样条函数在 x=x_j处的一阶导数m_j满足方程组:     

正算参数三次均匀B样条曲线的品质分析

<正> 由[1]可得出平面上正算参数三次均匀 B 样条曲线出现奇点和拐点的充要条件.据已知顶点 V 直接算出曲线值 C 的参数三次均匀 B 样条曲线方程     

不同核函数对光滑粒子流体动力学计算结果的影响分析

分别基于β-样条、新四次函数及五次样条三类光滑核函数,利用光滑粒子流体动力学方法,计算分析典型方腔剪切流和泊肃叶流动.通过对运动粒子空间分布、流动流速特性等分析,研究了不同光滑核函数对计算结果的影响.结果认为:基于β-样条及五次样条核函数的光滑粒子流体动力学计算结果标准误差分别为1.02×10~(-4)和1.54×10~(-4),能恰当的反应流体的流动状态,但基于五次样条核函数的计算结果易导致粒子分布的不均匀性,破坏计算的稳定性.而新四次核函数的计算结果精度及运算效率较低,与流体的流动特性存在较大差异.研究结果可为提高光滑粒子流体动力学的计算精度提供重要参考.     

二锥形分离杯之间血液流动稳定性的窄间隙理论

本文在获得血液分离器锥形分离杯内(其中一杯静止,另一杯以ω等角速旋转)血液流动的边界摄动解的基础上[1],采用窄间隙稳定性理论,证明了带轴向流的二锥形分离杯(其中一杯静止,另一杯以ω等角速旋转)之间旋转密度分层血液流动的稳定性.     

Exact solutions for the flow of an Oldroyd-B fluid due to an infinite flat plate

The unsteady flow of an Oldroyd-B fluid due to an infinite flat plate, subject to a translation motion of linear time-dependent velocity in its plane, is studied by means of the Laplace transform. The velocity field and the associated tangential stress corresponding to the flow induced by the constantly accelerating plate as well as those produced by the impulsive motion of the plate are obtained as special cases. The solutions that have been determined, in all accordance with the solutions established using the Fourier transform, reduce to those for a Newtonian fluid as a limiting case. The similar solutions for a Maxwell fluid are also obtained.     

Analytical solution for the space fractional diffusion equation by two-step Adomian Decomposition Method

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are increasingly used in modeling practical superdiffusive problems in fluid flow, finance and other areas of application. This paper presents the analytical solutions of the space fractional diffusion equations by two-step Adomian Decomposition Method (TSADM). By using initial conditions, the explicit solutions of the equations have been presented in the closed form and then their solutions have been represented graphically. Two examples, the first one is one-dimensional and the second one is two-dimensional fractional diffusion equation, are presented to show the application of the present technique. The solutions obtained by the standard decomposition method have been numerically evaluated and presented in the form of tables and then compared with those obtained by TSADM. The present TSADM performs extremely well in terms of efficiency and simplicity.     

Exact analytical solutions for the Poiseuille and Couette-Poiseuille flow of third grade fluid between parallel plates

Exact analytical solutions for the velocity profiles and flow rates have been obtained in explicit forms for the Poiseuille and Couette-Poiseuille flow of a third grade fluid between two parallel plates. These exact solutions match well with their numerical counter parts and are better than the recently developed approximate analytical solutions. Besides, effects of various parameters on the velocity profile and flow rate have been studied.     

Analytic solution for MHD Transient rotating flow of a second grade fluid in a porous space

This article looks into the unsteady rotating magnetohydrodynamic (MHD) flow of an incompressible second grade fluid in a porous half space. The flow is induced by a suddenly moved plate in its own plane. Both the fluid and plate rotate in unison with the same angular velocity. Analytic solution of the governing flow problem is obtained by using Fourier sine transform. Based on the modified Darcy's law, expression for velocity is obtained. The influence of pertinent parameters on the flow is delineated and appropriate conclusions are drawn. Several existing solutions of Newtonian fluid have been also deduced as limiting cases.     

Transition of MHD boundary layer flow past a stretching sheet

In this paper a study is carried out to understand the transition effect of boundary layer flow: (1) due to a suddenly imposed magnetic field over a viscous flow past a stretching sheet and (2) due to sudden withdrawal of magnetic field over a viscous flow past a stretching sheet under a magnetic field. In both the cases the sheet stretches linearly along the direction of the fluid flow. Governing equations have been non-dimensionalised and the non-dimensionalised equations have been solved using the implicit finite difference method of Crank–Nicholson type. Comparison between the steady state exact solutions and the steady state computed solutions has been carried out. Graphical representation of the dimensionless horizontal velocity, vertical velocity and local skin friction profiles of the steady state and unsteady state has been presented. Computation has been carried out for various values of the magnetic parameter M. The obtained results has been interpreted and discussed.     

Unsteady MHD two-phase Couette flow of fluid–particle suspension

The problem of two-phase unsteady MHD Couette flow between two parallel infinite plates has been studied taking the viscosity effect of the two phases into consideration. Unified closed form expressions are obtained for the velocities and the skin frictions for both cases of the applied magnetic field being fixed to either the fluid or the moving plate. The novelty of this study is that we have obtained the solution of the unsteady flow using the Laplace transform technique, D’Alemberts method and the Riemann-sum approximation method. The solution obtained is validated by assenting comparisons with the closed form solutions obtained for the steady states which have been derived separately and also by the implicit finite difference method. Graphical result for the velocity of both phases based on the semi-analytical solutions are presented and discussed. A parametric study of some of the physical parameters involved in the problem is conducted. The skin friction for both the fluid and the particle phases decreases with time on both plates until a steady state is reached, it is also observed to decrease with increase in the particle viscosity on the moving plate while an opposite behaviour has been noticed on the stationary plate.     

A Study of Multiple Solutions for the Navier-Stokes Equations by a Finite Element Method

In this paper, a finite element method is proposed to investigate multiple solutions of the Navier-Stokes equations for an unsteady, laminar, incompressible flow in a porous expanding channel. Dual or triple solutions for the fixed values of the wall suction Reynolds number $R$ and the expansion ratio $α$ are obtained numerically. The computed multiple solutions for the symmetric flow are validated by comparing them with approximate analytic solutions obtained by the similarity transformation and homotopy analysis method. Unlike previous works, our method deals with the Navier-Stokes equations directly and thus has no similarity and other restrictions as in previous works. Finally we use the method to study multiple solutions for three cases of the asymmetric flow (which has not been studied before using the similarity-type techniques).     

Exact solution of the Navier-Stokes equations for the oscillating flow in a duct of a cross-section of right-angled isosceles triangle

The Navier-Stokes equations have been solved in order to obtain an analytical solution of the fully developed laminar flow in a duct having a cross section of a right-angled, isosceles triangle. We obtained a solution for the case of oscillating pressure gradient flow. The pulsating flow is obtained by the superposition of the steady and oscillating pressure gradient solutions.     

Numerical solution of the free‐surface viscous flow on a horizontal rotating elliptical cylinder

The numerical solution of the free‐surface fluid flow on a rotating elliptical cylinder is presented. Up to the present, research has concentrated on the circular cylinder for which steady solutions are the main interest. However, for noncircular cylinders, such as the ellipse, steady solutions are no longer possible, but there will be periodic solutions in which the solution is repeated after one full revolution of the cylinder. It is this new aspect that makes the investigation of noncircular cylinders novel. Here we consider both the time‐dependent and periodic solutions for zero Reynolds number fluid flow. The numerical solution is expedited by first mapping the fluid film domain onto a rectangle such that the position of the free‐surface is determined as part of the solution. For the time‐dependent case a simple time‐marching method of lines approach is adopted. For the periodic solution the discretised nonlinear equations have to be solved simultaneously over a time period. The resulting large system of equations is solved using Newton's method in which the form of the Jacobian enables a straightforward decomposition to be implemented, which makes matrix inversion manageable. In the periodic case all derivatives have been approximated pseudospectrally with the time derivative approximated by a differentiation matrix which has been specially derived so that the weight of fluid is algebraically conserved. Of interest is the solution for which the weight of fluid is at its maximum possible value, and this has been obtained by increasing the weight until a consistency break‐down occurs. Time‐dependent solutions do not produce the periodic solution after a long time‐scale but have protuberances which are constantly appearing and disappearing. Periodic solutions exhibit spectral accuracy solutions and maximum supportable weight solutions have been obtained for ranges of eccentricity and angular velocity. The maximum weights are less than and approximately proportional to those obtained for the circular case. The shapes of maximum weight solutions is distinctly different from sub‐maximum weight solutions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008