高雷诺数流动的基本控制方程体系和扩散跑物化Navier-Stokes方程组的意义和用途

在计算机发达的时代, 高雷诺($Re$)数绕流计算中有无必要使用简化NS方程组, 本文讨论这个问题. 主要内容如下: (1)高$Re$数绕流包含3种基本流动: 所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动), 3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊, NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes (NS)和扩散抛物化(DP) NS方程组的数学性质彼此不同, 因此利用Euler-DPNS-NS方程组体系分析计算高$Re$数绕流流动就是一个合乎逻辑的选择, 该法与利用单一NS方程组的常用方法可以彼此检验和补充. (2)流体之间以及流体与外界的动量、能量和质量交换, 流态从层流到湍流的演化主要发生在干扰剪切流动中, 干扰剪切流及其最简单控制方程------DPNS方程组具有基础意义; DPNS方程组笔者在1967年已提出. (3)诸简化NS方程组: DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系, 它们的历史贡献和今后的用途; 它们的数学性质均为扩散抛物型, 但它们包含的黏性项彼此有所不同; 从流体力学角度来看, 它们中只有DPNS方程组能够准确描述干扰剪切流动. 提出把诸简化NS方程组统一为DPNS方程组的建议. (4)干扰剪切流------DPNS方程组与无干扰剪切流------边界层方程组之间的关系以及进一步研究干扰剪切流的意义.      

层流射流流动的简化Navier-Stokes方程(SNSE)解

本文利用十一种简化 Navier-Stokes 方程(SNSE) 求解已知Navier-Stokes(NS)方程准确解的层射流流动,表明:多数SNSE~([1-6])的解与NS方程的准确解不一致;少数SNSE~([7,8])的解与NS方程的准确解一致,文中在射流的喉部和拐点位置,给出几种SNSE解与准确解的相对偏差,并把粘性及惯性诸项加以定量比较,强调指出:按照边界层理论量级分析为Re~(1/2)和Re~1量级的惯性项以及Re~(-1/2)量级的粘性项具有重要影响;据此从力学角度论证了简化 NS 方程时,保留全部惯性项和合理取舍粘性项的必要性。     

叶轮机械流道中粘性可压缩流动的数值计算——(Ⅱ)旋转叶列中轴对称流动的计算

本文把第一部分中所提出的计算叶轮机械流道中粘性可压缩定常层流问题的数值方法应用到旋转坐标系中。首先导得旋转叶列中粘性可压缩流动所满足的基本方程组,然后通过在任意非正交曲线方向上的焓梯度方程、能量方程和熵方程之间的迭代计算,得到叶轮机械旋转叶列流道中粘性轴对称流动的数值解。     

粘流-无粘干扰流动(IF)理论

对不可压缩层流二维干扰流动,本文提出一个干扰流动(IF)理论。IF理论要点为:1)干扰流动沿主流的法向被分为三层即粘性层、干扰层和无粘层,引进了法向动量交换为主导过程的干扰层概念。2)利用力学守恒律、三层匹配关系及文中引进的干扰模型,把三层的空间尺度及惯性-粘性诸力的数置级表示为单参数m的函数,m<1/2·3)导出描述各层流动的控制方程、导出描述全城流动的控制方程为简化Navie-Stokes(SNS)方程。IF理论适用于不存在分离的附着干扰流动以及存在分离的大范围干扰流动,经典边界层(CBL)理论和流动分离局部区域Triple-Deck(TD)理论分别是本文理论在参数m=O和1/4时的两个特例,本文理论容易推广到可压缩、三维及湍流流动。     

NND格式在理想磁流体方程组中的应用

针对守恒型磁流体力学方程组（MHD）和流体力学方程组（HD）通量项不同特点,提出了一种能够采用无振荡、无自由参数（NND）格式离散MHD方程组的通量分裂方法,并首先在一维模型方程中验证了方法的可行性,进一步全三维离散了MHD方程组,在轴对称盔形磁场位形太阳风流动的数值试验中,选取46个太阳半径（Rs)的计算域,其中能够反映行星际空间物理参数在径向有大到8－9个量级变化的特点。计算结果表明针对气动力学跨音速流动的NND格式可以推广到磁体力学方程组中,并有很好的稳定性。     

对流-扩散相互作用结构的不变性

本文提出并证明了不可压缩剪切层流中对流-扩散相互作用结构不变性诸定理:即二难剪切层流与其线性化及非线性扰动存在同一的对流-扩散相互作用结构,且物理尺度(指时间、空间和速度尺度)相同。给出十个推论,例如:对流-扩散相互作用可在剪切层流及其扰动场内“激发“快时间尺度和小空间尺度结构,线性化稳定性原理的约定对剪切流体系统成立等。应用题例导出计及时间-空间尺度效应和非平行流效应的广义Orr-Sommerfeld(GOS)方程,证实它有两个粘性解:阻尼层解和干扰层解;经典OS方程及其两个粘性解:边界层解和Heisenberg临界层解,Triple-deck稳定性理论基本方程及其两个粘性解,均是本文GOS方程及其两个粘性解的特例。     

配置点谱方法-人工压缩法(SCM-ACM)求解同心圆筒内流体流动

发展了配置点谱方法SCM(Spectral collocation method)和人工压缩法ACM(Artificial compressibility method)相结合的SCM-ACM数值方法,计算了柱坐标系下稳态不可压缩流动N-S方程组。选取典型的同心圆筒间旋转流动Taylor-Couette流作为测试对象,首先,采用人工压缩法获得人工压缩格式的非稳态可压缩流动控制方程;再将控制方程中的空间偏微分项用配置点谱方法进行离散,得到矩阵形式的代数方程;编写了SCM-ACM求解不可压缩流动问题的程序;最后,通过与公开发表的Taylor-Couette流的计算结果对比,验证了求解程序的有效性。结果证明,本文发展的SCM-ACM数值方法能够用于求解圆筒内不可压缩流体流动问题,该方法既保留了谱方法指数收敛的特性,也具有ACM形式简单和易于实施的特点。本文发展的SCM-ACM数值方法为求解柱坐标下不可压缩流体流动问题提供了一种新的选择。     

锥形端部弹体在岩石(混凝土)介质层中侵彻实用计算方法

采用刚塑性不可压缩的介质模型 ,用极限分析理论的上限方法 ,通过建立动力学许可速度场得到介质对弹体侵彻的静阻力分量 ,通过破碎介质动量守恒条件得到弹体侵彻的动阻力分量 ,在模型中还考虑了介质的尺度关系和弹体的弹头形状两种影响因素。侵彻过程中的速度、加速度、阻力和经历时间通过弹体的运动方程增量计算得到。通过与几种常用经验公式比较说明了本文方法的实用和可靠性。     

湍流中的层次结构和标度律

湍流是由各种不同尺度、不同幅度和不同相干度的层次结构组成的．在从积分尺度到耗散尺度连续分布的每个尺度上,都存在具有最大振幅和最高相干度的结构──最强间歇结构．幅度较小的结构按照一定的层次对称关系与最强间歇结构相联系．这就是连接多尺度和多幅度脉动结构的层次结构模型（She和Leveque,Phys．Rev．Lett,1994;72,336）．本文对湍流层次结构理论的内容、基本观念及其新发展进行了一些综述和讨论,希望对国内的湍流基础研究能起一定的推动作用．      

叶轮机械流道中粘性可压缩流动的数值计算——(Ⅰ)静止叶列中轴对称流动的计算

本文提出一种计算叶轮机械流道中粘性可压缩定常层流问题的数值方法,通过在任意非正交曲线方向上的焓梯度方程、能量方程和熵方程之间的迭代计算,可以得到整个流道中粘性可压缩定常流动问题的数值解。本文首先在静止坐标系中进行分析和讨论,并描述了叶轮机械的静止叶列中轴对称流动的计算方法。计算实例表明,本方法的特点是简单明了,计算速度快,可以广泛地应用于工程设计之中。     

Hierarchial structure theory for basic equations of fluid mechanics (BEFM) and the simplified Navier-Stokes equations (SNSE)

A hierarchial structure for the basic equations of fluid mechanics (BEFM) is found through the analysis of scales of length and time that proves a measure of the rate of change of the quantities describing the motion of the fluid as well as an estimation of the order of magnitude of various terms included in BEFM. The hierarchial structure theory shows that if (1) the characteristic Reynolds numbersRe is larger than unity and (2) the length scale in one coordinate direction is larger than that in other coordinate directions. BEFM can be classified into some levels according to the estimation of the order of magnitude of various terms included in BEFM. The hierarchial structure of BEFM has two branches: one is from BLE- to BEFM inner hierarchy, the other is from EE- to BEFM outer hierarchy, where BLE and EE are abbreviations of the boundary-layer equations and of Euler equations, respectively. The relationship between the two branches of the hierarchial structure, the characteristics, subcharacteristics and mathematical properties of the hierarchial equations are studied. A comparison between the present hierarchial equations and the Simplified Navier-Stokes equations (SNSE) appeared in literatures is also made. BLE-, EE-and Inner-outer-matched (IOM) equations hierarchies are the most important and useful three levels for solving viscous flow-fields approximately.     

SIMILARITY SOLUTIONS OF BOUNDARY LAYER EQUATIONS FOR A SPECIAL NON-NEWTONIAN FLUID IN A SPECIAL COORDINATE SYSTME

Two dimensional equations of steady motion for third order fluids are expressed in a special coordinate system generated by the potential flow corresponding to an inviscid fluid. For the inviscid flow around an arbitrary object, the streamlines are the phicoordinates and velocity potential lines are psi-coordinates which form an orthogonal curvilinear set of coordinates. The outcome, boundary layer equations, is then shown to be independent of the body shape immersed into the flow. As a first approximation, assumption that second grade terms are negligible compared to viscous and third grade terms. Second grade terms spoil scaling transformation which is only transformation leading to similarity solutions for third grade fluid. By ~sing Lie group methods, infinitesimal generators of boundary layer equations are calculated. The equations are transformed into an ordinary differential system. Numerical solutions of outcoming nonlinear differential equations are found by using combination of a Runge-Kutta algorithm and shooting technique.     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

具有精确色散性的非线性波浪数学模型

以完全非线性的自由表面边界条件为基础,以波面升高\eta和自由表面速度势\phi _\eta为待求变量,建立了新的波浪方程．方程在色散性上是完全精确的,非线性近似至三阶．与缓坡方程相比较,两者都具有精确的色散性,但该方程属于非线性模型,可模拟波浪的非线性效应,且适用于不规则波．方程的特点是属于微分-积分方程,对如何处理方程中积分项进行了讨论,并数值模拟了不同周期的线性波和二阶Stokes波,也模拟了波群的非线性演化,以对模型进行验证．      

Magnetohydrodynamic (MHD) flow of a second grade fluid in a channel with porous wall

An analysis has been carried out to study the effect of magnetic field on an electrically conducting fluid of second grade in a parallel channel. The coolant fluid is injected into the porous channel through one side of the channel wall into the other heated impermeable wall. The combined effect of inertia, viscous, viscoelastic and magnetic forces are studied. The basic equations governing the flow and heat transfer are reduced to a set of ordinary differential equations by using appropriate transformations for velocity and temperature. Numerical solutions of these equations are obtained with the help of Runge-Kutta fourth order method in association with quasi-linear shooting technique. Numerical results for velocity field, temperature field, skin friction and Nusselt number are presented in terms of elastic parameter, Hartmann number, Prandtl number and Reynolds number. Special case of our results is in good agreement with earlier published work.     

论庞加莱-契达耶夫方程

研究表明：庞加莱－契达耶夫正则方程是非正则变量下相当普遍的哈密顿方程．这表明，多余坐标下的广义拉格朗日方程和广义哈密顿方程（其阶数低于带有不定乘子的方程），以及准坐标下的欧拉－拉格朗日方程，都是庞加莱－契达耶夫方程的特殊情况；从而，可将其理论推广到上述系统．而且还研讨了庞加莱－契达耶夫方程在非完整系动力学中的应用问题．     

扩散抛物化Navier-Stokes方程数值解法评述

20世纪60年代末期在边界层理论基础上发展起来的各种简化Navier-Stokes (N-S)方程(统称为扩散抛物化N-S方程)及其算法, 较为彻底地解决了无黏流及黏流的相互干扰问题, 并为高雷诺数大型复杂黏性流场的数值模拟开辟了新的途径. 本文将系统地评述这一领域的主要成果, 包括各种简化N-S模型的优缺点; 数学奇性及正则化方法; 代表性的数值解法以及最近几年的新进展.      

Analyzing the 2D fracture problems via the enriched boundary element-free method

In this paper, the enriched boundary element-free method for two-dimensional fracture problems is presented. An improved moving least-squares (IMLS) approximation, in which the orthogonal function system with a weight function is used as the basis function, is used to obtain the shape functions. The IMLS approximation has greater computational efficiency and precision than the existing moving least-squares (MLS) approximation, and does not lead to an ill-conditioned system of equations. Combining the boundary integral equation (BIE) method and the IMLS approximation, a boundary element-free method (BEFM), for two-dimensional fracture problems is obtained. For two-dimensional fracture problems, the enriched basis function is used at the tip of the crack, and then the enriched BEFM is presented. In comparison with other existing meshless boundary integral equation methods, the BEFM is a direct numerical method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be implemented easily, which leads to a greater computational precision. When the enriched BEFM is used, the singularity of the stresses at the tip of the crack can be shown better than that in the BEFM. For the purposes of demonstration, some selected numerical examples are solved using the enriched BEFM.     

Computer-aided heuristic analysis of fluid models

A new kind of symbolic program to aid the heuristic simplification of fluid models is presented. The program, AOM, employs order of magnitude analysis and method of dominant balance to generate simplified models. It has two novel features: (1) it uses heuristic techniques to decide what equations to solve and what algebra to do, and (2) it explains its deduction steps. The basic operation of AOM consists of five steps: (1) assign order of magnitude estimates to terms in the equations, (2) find maximal terms of each equation, i.e., terms that are not dominated by any other terms in the same equation, (3) consider all possible n-term dominant balance assumptions, (4) propagate the effects of the balance assumptions, and (5) remove partial models based on inconsistent balance assumptions. AOM also exploits constraints among equations and submodels to simplify complicated fluid models such as the triple-deck equations. Three annotated examples are presented to explain the operations of AOM. The implications for the development of computer-aided analysis programs for fluid dynamics and education are discussed.This research was funded in part by NSF NYI Award ECS-9357773 and CCR-9109567.