有限长圆柱轴承非稳态油膜力的级数解法

本文用分离变量法求解雷诺方程，在π油膜的假设下，求得雷诺方程应满足的特征值与用傅立叶级数表达的特征函数，进而求得有限长轴的非特急油膜力解析表达式。为分析轴承转子系统的非线性动力特性提供了帮助。     

有限长滑动轴承非线性油膜力的近似解法

本文中提出了一种求解有限长径向滑动轴承非线性油膜力的近似解析方法.在滑动轴承-转子系统非线性动力行为分析中,油膜力计算模型通常采用"π"油膜假设,但是,实际工况中油膜的存在区域并非是"π"区域,运行时油膜中出现气穴,破裂成条纹状(即具有Reynolds边界条件).本文中的近似解析方法采用Reynolds边界条件,基于变分原理,运用分离变量法求解油膜的压力分布,其中油膜压力的周向分离函数通过无限长轴承的油膜压力分布获得,油膜的破裂终止位置角通过连续条件确定,轴向分离函数运用变分原理并结合周向函数求得.计算结果表明:本文中提出的方法和有限元方法的结果吻合得很好.在此基础上,分析了一些轴承参数对油膜压力分布的影响.     

用瑞利-李兹法求解瞬时非稳态滑动轴承油膜力的新算法

提出了一个利用瑞利一李兹方法求解Reynolds边界条件下非稳态滑动轴承油膜力的近似算法,充分利用油膜力分布的特性,用双曲余弦函数来表示油膜力的轴向压力分布,而用多项式函数插值法来求解油膜力的周向压力分布,并同时计算出油膜力的破裂边界。算例表明本算法达到了相当高的精度,可用于转子系统的非线性数值分析,能大大降低数值求解瞬态油膜力的计算时间。     

基于有限差分-谱方法的分数阶Burgers流体的非稳态驻点流动

研究了分数阶Burgers流体通过拉伸平板的非稳态驻点流动问题。将分数阶导数引入Burgers流体模型可以更好地模拟流动过程,但也增加了模型的复杂性和求解难度。首次运用有限差分-谱方法求解分数阶Burgers流体模型,离散格式构造简单有效。采用谱方法对控制方程中的空间项进行离散,利用有限差分方法分别结合L-1和L-2算法离散控制方程中的时间项,给出了两种离散格式,并且通过构造数值算例证明了离散格式的收敛性。结果表明,在靠近平板处,速度随着分数阶导数的增加而减小,而无穷远处的流体速度呈现出相反的趋势,体现了分数阶导数的记忆特性。此外,雷诺数越小,流体的粘度越大,导致流体速度越大。由于松弛时间参数的松弛特性,靠近平板处松弛时间参数对速度分布有抑制作用,远离平板处松弛时间促进流体流动。     

非稳态沙波纹流场的数值模拟

风成沙波纹是沙质地表在风力作用下形成的最基本的地貌类型.本文研究了沙波纹在形成过程中,风场对波纹的形态变化的响应.采用目前比较流行的有限体积法,通过求解非稳态情况下的Navier-Stokes方程来计算风场;地表形态采用阶梯网格进行逼近,为了对地表形态变化进行有效的控制,根据相关参数采用正弦函数作为沙波纹的基本形状;最后根据计算结果详细讨论了沙波纹从形成到稳定的发展过程中地表流场的行为变化特征.     

非稳态不可压缩自由表面粘性流体的三维数值计算

扩展了相对体积算法，计算了变速旋转的敞口圆筒内水的真实非稳态流动。采用了原始变量、交错非等分网格和显式迭代。计算结果与实验现象相吻合。当圆筒长时间等速旋转，其内流体与筒体一起作刚性旋转时，计算自由表面形状与流体力学理论公式的预言吻合得很好。     

钢筋混凝土板壳材料非线性分析的有限条法

本文发展了变厚度任意配筋的钢筋混凝土板壳问题中考虑材料非线性分析的有限条法。文中建立了沿厚度方向作分层离散的三维蜕化二次厚壳条，而在计算上，沿壳条跨度方向采用Ｓｉｍｐｓｏｎ积分技术可以方便地处理该方向的任意配筋并给出了满足精度要求的选取Ｓｉｍｐｓｏｎ积分点数的表达式。文中提出的“拟刚度矩阵”法大大简化了非线性有限条法的求解过程。给出两个数值算例表明了理论方法的有效性及其精度。     

基于不同时间步长时域非结构有限体积法模拟声-弹性耦合问题

介绍一种改进的时域非结构有限体积法(FVM),并将其应用于声-弹性耦合问题.在流体与固体介质中分别求解声波动方程与弹性波方程,根据交界面上的力平衡与质点振速连续条件考虑二者的相互作用.同时考虑双线性四边形单元的线性变化项及常数项,并结合常应变三角形单元处理混合网格问题.分别对三角形单元和四边形单元进行色散分析,给出声波动方程的稳定性条件.在不同介质中采用不同时间步长,提高计算效率.求解弹性波问题、声-弹性耦合问题,结果表明,改进后的方法求解声-弹性耦合问题是有效和准确的,并且具有良好的数值稳定性.     

静水压力下环加肋圆柱壳总体与局部失稳分析的复合有限条法

本文提出了一种计算环加肋圆柱壳稳定性的新方法─—复合有限条法。该方法能够在一个复合有限条元中计入若干个横向加肋的影响，井能考虑肋骨的偏心，是一种子结构方法。算例表明，复合有限条法具有很高的求解精度和效率，是一种分析加肋柱壳结构的较理想方法。该方法不仅能够计算环加肋圆柱壳的总体失稳临界载荷，而且可计算其局部失稳临界载荷。     

移动荷载下路面板-层状多相弹性地基的稳态响应

考虑路面板和地基的相互作用,将路面板作为三维弹性体,建立路面板-层状多相弹性地基动力响应的三维多相层状弹性半空间体模型,分析了移动荷载作用下路面板-多相弹性地基系统的动力响应.将移动单元法引入到两相饱和弹性介质的半解析方法中,构造了随荷载按照相同速度运动的移动层单元,基于移动坐标下两相饱和弹性介质的动力控制方程和边界条件,应用加权残数法建立了两相弹性介质移动层单元动力方程,该方程可退化为单相弹性介质移动层单元动力方程,基于此,建立了移动荷载下多相弹性地基与路面板系统的三维动态响应的统一的半解析方程.并以两相饱和弹性半空间地基-单相弹性地基-路面板所形成的路面结构为例,数值分析了荷载速度、饱和层渗透系数、弹性层厚度等参数对路面板位移和土体孔压响应的影响.研究结果表明移动单元法是研究移动荷载下路面结构动态响应的一种有效的方法.     

Dynamic response of plates due to moving vehicles using finite strip method

ＩｎｔｒｏｄｕｃｔｉｏｎＤｙｎａｍｉｃｒｅｓｐｏｎｓｅｏｆｂｅａｍ＿ｌｉｋｅｓｔｒｕｃｔｕｒｅｓｔｏｍｏｖｉｎｇｖｅｈｉｃｌｅｓｈａｓｂｅｅｎｅｘｔｅｎｓｉｖｅｌｙｓｔｕｄｉｅｄ .Ｈｏｗｅｖｅｒ,ｔｈｅｓｔｕｄｙｏｎｄｙｎａｍｉｃｒｅｓｐｏｎｓｅｏｆｐｌａｔｅｓｔｏｍｏｖｉｎｇｖｅｈｉｃｌｅｓｈａｓｓｏｆａｒｒｅｃｅｉｖｅｄｂｕｔｓｃａｎｔａｔｔｅｎｔｉｏｎ .Ｗｕｅｔａｌ.[1],ＴａｈｅｒｉａｎｄＴｉｎｇ[2 ],ＹｅｎｅｒａｎｄＣｈｏｍｐｏｏｍｉｎｇ[3],ＨｕｍａｒａｎｄＫａｓｈｉｆ[4 ]ｅｔａｌ.ｓｅｐａ…     

Upwind finite difference method for miscible oil and water displacement problem with moving boundary values

The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.     

Fixed domain methods for free and moving boundary flows in porous media

A survey is presented concerning fixed domain methods used to solve mathematical models of free and moving boundary flow problems in porous media. These include the following: variational inequality or quasi-variational inequality formulations; general inequality formulations which have been set and solved in fixed domains; and the residual flow procedure. Finally, some parallel computing methods and mesh adaptation methods are discussed to demonstrate how these fixed domain formulations can be solved with current technology.The fixed domain methods that are referenced herein can be classified into two groups: the variational inequality method and the extended pressure head method. Baiocchi was the first to apply the variational inequality method to free boundary problems of flows through porous media. This method in general also uses an extension of the pressure head but adds an application of an integral transformation (a Baiocchi transformation) to the problem. The method possesses a beautiful mathematical structure for its theory and yields simple numerical solution algorithms. However, application of the method is difficult if not impossible in some cases depending upon the regularity of the seepage domain.The extended pressure head method is based on the concept that the pressure is extended smoothly across the free or moving boundary into the unsaturated region from the flow domain. The extension of the pressure head to the entire porous medium yields an extended coefficient of permeability of the medium which is equal to the saturated coefficient in the seepage region and is equal to zero or some small value (for computational purposes) in the unsaturated region.     

A non-iterative transformation method for Newton's free boundary problem

In book II of Newton's Principia Mathematica of 1687 several applicative problems are introduced and solved. There, we can find the formulation of the first calculus of variations problem that leads to the first free boundary problem of history. The general calculus of variations problem is concerned with the optimal shape design for the motion of projectiles subject to air resistance. Here, for Newton's optimal nose cone free boundary problem, we define a non-iterative initial value method which is referred in the literature as a transformation method. To define this method we apply invariance properties of Newton's free boundary problem under a scaling group of point transformations. Finally, we compare our non-iterative numerical results with those available in the literature and obtained via an iterative shooting method. We emphasize that our non-iterative method is faster than shooting or collocation methods and does not need any preliminary computation to test the target function as the iterative method or even provide any initial iterate. Moreover, applying Buckingham Pi-Theorem we get the functional relation between the unknown free boundary and the nose cone radius and height.     

Characteristic finite difference method and application for moving boundary value problem of coupled system

The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l^2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.     

SOLVING THE FREE BOUNDARY PROBLEM IN CONTINUOUS CASTING BY USING BOUNDARY ELEMENT METHOD

ＳＯＬＶＩＮＧＴＨＥＦＲＥＥＢＯＵＮＤＡＲＹＰＲＯＢＬＥＭＩＮＣＯＮＴＩＮＵＯＵＳＣＡＳＴＩＮＧＢＹＵＳＩＮＧＢＯＵＮＤＡＲＹＥＬＥＭＥＮＴＭＥＴＨＯＤＬｉＹａｏｙｏｎｇ（李耀勇）；ＺｈａｎｇＺｈｉｌｉ（张自立）（ＲｅｃｅｉｖｅｄＪｕｎｅ，１８，１９...     

Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.     

Modified characteristic finite difference fractional step method for moving boundary value problem of nonlinear percolation system

A fractional step scheme with modified characteristic finite differences running in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of difference operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in l ² norm is displayed to complete the convergence analysis of the numerical algorithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.     

基于移动相似中心的比例边界有限元方法

在传统的比例边界有限元中,相似中心是固定的,难以用其求解关于偏心域的场问题。本文引入移动相似中心的概念,建立新的比例边界坐标变换,并利用加权余量法将控制方程半弱化为关于径向坐标的二阶常微分方程,引入对偶变量,将其降为系数矩阵为Hamilton矩阵的一阶常微分方程。对Hamilton矩阵进行Schur分解,得到微分方程的通解,代入边值条件可得关于积分常数的代数方程。此方法将比例边界有限元扩展到偏心域的边值问题,同时在径向是半解析的,解的精度高;仅需要离散求解域的一个边界,数据量小;在计算中仅需要对Hamilton矩阵进行Schur分解以及求解关于积分常数的代数方程,运算量少。将偏心环形域静电场边值问题的算例与解析解或其他数值方法计算结果的比较,表明此方法具有精度高、数据量小及运算量小的优点。     

A diffuse-interface immersed boundary method for simulation of compressible viscous flows with stationary and moving boundaries

In this paper, a diffuse-interface immersed boundary method (IBM) is proposed for simulation of compressible viscous flows with stationary and moving boundaries. In the method, the solution of flow field and the implementation of boundary conditions are decoupled into two steps by applying the fractional step technique, ie, the predictor step and the corrector step. Firstly, in the predictor step, the intermediate flow field is resolved by a recently developed gas kinetic flux solver (GKFS) without consideration of the solid boundary. The GKFS is a finite volume approach that solves the Navier-Stokes equations for the flow variables at cell centers. In GKFS, the inviscid and viscous fluxes are evaluated as a single entity by reconstructing the local solution of continuous Boltzmann equation. Secondly, in the corrector step, the intermediate flow field is corrected by the present diffuse-interface IBM. During this process, the velocity field is firstly corrected by the implicit boundary condition–enforced IBM so that the no-slip boundary condition can be accurately satisfied. After that, the density correction is made by an iterative approach with the help of the continuity equation. Finally, the correction of the temperature field is made in the same way as that of the velocity field. Good agreements between the present simulations and the reference data in literature demonstrate the reliability of the proposed method.