分析力学初值问题的一种变分原理形式

明确了分析力学初值问题的控制方程，按照广义力和广义位移之间的对应关系，将 各控制方程卷乘上相应的虚量，代数相加，进而在 原空间中建立了分析力学初值问题的一种变分原理形式，即建立了分析力学初值问题的卷积 型变分原理和卷积型广义变分原理. 推导了分析力学初值问题卷积型变分原理和卷积型广义 变分原理的驻值条件. 在建立分析力学初值问题的一种变分原理形式的同时, 将变积方法推广为卷变积方法.     

The Opposing Effect of Viscous Dissipation Allows for a Parallel Free Convection Boundary-Layer Flow Along a Cold Vertical Flat Plate

External free convection boundary-layer flows are usually treated by neglecting the effect of viscous dissipation. This assumption always results in a non-parallel flow, besides a strong parallel component also a weak transversal component of the (steady) velocity field occurs. The present paper shows, however, that the weak opposing effect of the buoyancy forces due to heat release by viscous dissipation, can give rise along a cold vertical plate adjacent to a fluid saturated porous medium to a strictly parallel steady free convection flow. This boundary-layer flow is described by an algebraically decaying exact analytical solution of the basic balance equations.     

用离散速度方法计算浅水长波方程

用离散速度法计算浅水波方程，将空气动力学方程和浅水波方程作了比较，用Nadiga提出的近平衡流动方法模拟浅水波方程的连续和间断解。计算了一维的溃坝波问题和Thacker提出的连续解问题，结果与精确解作了比较，并且计算了水流跃过障碍物的问题。     

Numerical solutions of incompressible Euler and Navier-Stokes equations by efficient discrete singular convolution method

An efficient discrete singular convolution (DSC) method is introduced to the numerical solutions of incompressible Euler and Navier-Stokes equations with periodic boundary conditions. Two numerical tests of two-dimensional Navier-Stokes equations with periodic boundary conditions and Euler equations for doubly periodic shear layer flows are carried out by using the DSC method for spatial derivatives and fourth-order Runge-Kutta method for time advancement, respectively. The computational results show that the DSC method is efficient and robust for solving the problems of incompressible flows, and has the potential of being extended to numerically solve much broader problems in fluid dynamics. The project supported by the National Natural Science Foundation of China (No.19902010).     

Comparison of Two Mathematical Models for 3D Groundwater Flow: Block-Centered Heads and Edge-Based Stream Functions

Traditionally, groundwater flow models, as well as oil reservoir models, are based on the block-centered finite difference method. Well-known models based on this approach are MODFLOW (groundwater) and ECLIPSE (oil and gas). Such models are well proven and robust; their underlying principles are well understood by hydrologists and petroleum reservoir engineers. Nevertheless, the desire to improve the block-centered finite difference paradigm has always been alive, for instance, to be able to apply deformed grid blocks, or to model anisotropy that is not aligned along the coordinate axes. This article introduces the edge-based stream function as a potential alternative to the paradigmatic model, not only to mitigate the above mentioned limitations, but especially for its promise to inverse modeling. Computer programs have been developed for the discrete analog equations of the stream function method and the conventional method. The two methods are tested by using synthetic forward modeling problems of uniform and radial flow. The theoretical formulation and the numerical results show that the two methods are algebraically equivalent and yield the same flux output. However, for rectangular grid blocks and anisotropy aligned along the coordinate axes, the block-centered method is shown to be computationally more efficient than the edge-based stream function method. The major advantage of the stream function method is that it is linear in the resistivities, proving it an ideal candidate for direct inverse modeling. Moreover, any arbitrary specification of stream functions yields a solution that satisfies the mass balance.     

Analytical development of disturbed matrix eigenvalue problem applied to mixed convection stability analysis in Darcy media

This work consists in evaluating algebraically and numerically the influence of a disturbance on the spectral values of a diagonalizable matrix. Thus, two approaches will be possible; to use the theorem of disturbances of a matrix depending on a parameter, due to Lidskii and primarily based on the structure of Jordan of the no disturbed matrix. The second approach consists in factorizing the matrix system, and then carrying out a numerical calculation of the roots of the disturbances matrix characteristic polynomial. This problem can be a standard model in the equations of the continuous media mechanics. During this work, we chose to use the second approach and in order to illustrate the application, we choose the Rayleigh–Bénard problem in Darcy media, disturbed by a filtering through flow. The matrix form of the problem is calculated starting from a linear stability analysis by a finite elements method. We show that it is possible to break up the general phenomenon into other elementary ones described respectively by a disturbed matrix and a disturbance. A good agreement between the two methods was seen. To cite this article: H.B. Hamed, R. Bennacer, C. R. Mecanique 336 (2008).     

The Edge-Based Face Element Method for 3D-Stream Function and Flux Calculations in Porous Media Flow

We present a velocity-oriented discrete analog of the partial differential equations governing porous media flow: the edge-based face element method. Conventional finite element techniques calculate pressures in the nodes of the grid. However, such methods do not satisfy the requirement of flux continuity at the faces. In contrast, the edge-based method calculates vector potentials along the edges, leading to continuity of fluxes. The method is algebraically equivalent with the popular block-centered finite difference method and with the mixed-hybrid finite element method, but is algorithmically different and has the same robustness as the more conventional node-based velocity-oriented method. The numerical examples are computed analytically and may, therefore, be considered as an 'heuristic proof' of the theory and its practical applicability for reservoir engineering and geohydrology.     

Thermodynamic properties of equilibrium air by an exact two-equation model for Euler and Navier–Stokes CFD algorithms

This paper details an exact two-equation procedure to generate pressure, temperature and mass and mole fractions as well as their thermodynamic and Jacobian partial derivatives for five-species neutral equilibrium air. Applicable for arbitrary forms of equilibrium constants and especially designed for explicit and implicit CFD algorithms, the procedure algebraically reduces to two equations the six-equation thermodynamic system comprising the equations for internal energy, law of mass action and conservation of species mass and ratio of oxygen and nitrogen nuclei. This exact algebraic reduction explicitly expresses four mass fractions in terms of nitric oxide mass fraction and temperature, which are then determined through a rapidly converging numerical solution of the internal energy and nitric oxide mass action equations. The procedure then exactly determines the partial derivatives of pressure, temperature and mass fractions analytically. The mathematical formulation also introduces a convenient system non-dimensionalization that makes the procedure uniformly applicable to flows ranging from shock tube flows with zero initial velocity to aerothermodynamic flows with supersonic/hypersonic freestream Mach numbers. Over a wide range of density and internal energy the predicted distributions of mole fractions for the model five species agree with independent published results, while pressure and temperature as well as their partial derivatives remain continuous, smooth and physically meaningful. © 1997 John Wiley & Sons, Ltd.     

求解变分型积分方程的一种新型数值方法——有限变分法

将作者提出的多虚拟裂纹扩展法(MVCE法)拓展为求解变分型积分方程问题的一种新型数值方法——有限变分法(FVM)。它的基本思想是,给定有限个(N个)局部变分模式,将所求解的未知量用适当的方法离散化,针对这N个局部变分模式列出N个方程,求解N个未知系数,从而求得未知量。单一未知变量FVM的最终方程组的系数矩阵通常是一个对称的窄带矩阵,对角元是大数,有很好的数值计算性能。用FVM求解了三维I型裂纹前缘的应力强度因子(SIF)分布。利用基于FVM的通用权函数法计算程序,可以高精度和高效率地求解表面力、体积力和温度载荷共同作用情况下三维裂纹前缘SIF的分布及其时间历程。FVM可以被推广到更广泛的领域,是一个求解变分型积分方程问题的普遍适用的新型数值方法。     

数值模拟离子强度敏感水凝胶的多场特性

发展了多物理模型来研究溶液中离子强度敏感水凝胶电-力-化学多场耦合的特性. 模 型的主要控制方程包括：计算水凝胶内外离子浓度分布的Nernst-Planck化学场方程；描述 膨胀变形的力学场方程和描述电场的泊松方程. 无网格有限云团法和牛顿迭代法用来数值离 散和求解控制方程. 通过对比多场耦合的响应, 包括胶的膨胀率和胶内外离子浓度和电势的分布，探讨了影响胶体变形的主要因素. 对数值模拟结果和实验结果进行了对照.     

Compatibility between finite volumes and finite elements using solutions of shallow water equations for substance transport

This paper formulates a finite volume analogue of a finite element schematization of three‐dimensional shallow water equations. The resulting finite volume schematization, when applied to the continuity equation, exactly reproduces the set of matrix equations that is obtained by the application of the corresponding finite element schematization to the continuity equation. The procedure allows the consistent and mass conserving coupling of the finite element Telemac model for three‐dimensional flow with the finite volume Delft3D‐WAQ model for water quality. The work has been carried out as part of a joint development by LNHE and WL∣Delft Hydraulics to explore the mutual interaction of their software. Copyright © 2006 John Wiley & Sons, Ltd.     

非压缩性流体的湍流理论

 本文是作者50年的湍流研究工作的一个总结. 整个研究工作分成 4个阶段. 第一阶段是在上一世纪的40年代，提出研究湍流必须研究脉动 速度场，给出了湍流模式理论的基础. 第二阶段是上一世纪的50年代到 60年代，提出了先求解后平均的旋涡结构的湍流统计理论. 第三阶段是 上一世纪的70年代，他和黄永念提出了准相似性假设，完成了均匀各向 同性湍流统计理论的工作. 第四阶段是上一世纪的80年代到90年代初， 提出了广义准相似性假设，给出了处理一般具有剪应力的湍流流动的奇 阶截断法和逐级逼近法.     

On a new algorithm of constructing solitary wave solutions for systems of nonlinear evolution equations in mathematical physics

ＩｎｔｒｏｄｕｃｔｉｏｎＳｏｌｖｉｎｇｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓｉｓａｌｗａｙｓａｖｅｒｙｉｎｔｅｒｅｓｔｉｎｇｓｕｂｊｅｃｔｆｏｒｍａｔｈｅｍａｔｉｃｉａｎａｎｄｐｈｙｓｉｃｉａｎ ,ｉｎｐａｒｔｉｃｕｌａｒ,ｓｏｌｉｔａｒｙｗａｖｅｓｏｌｕｔｉｏｎｓｆｏｒｎｏｎｌｉｎｅａｒｅｑｕａｔｉｏｎｓａｒｅｏｆｂｏｔｈｔｈｅｏｒｅｔｉｃａｌａｎｄｐｒａｃｔｉｃａｌｉｍｐｏｒｔａｎｃｅ.Ｒｅｃｅｎｔｌｙ ,Ｙａｎ[1]ｏｂｔａｉｎｅｄａｔｒａｎｓｆｏｒｍａｔｉｏｎｄｉｒｅｃｔｌｙｆｒｏｍｔｈｅｆａｍｏ…     

MHD flow of upper-convected Maxwell fluid over porous stretching sheet using successive Taylor series linearization method

This paper investigates the magnetohydrodynamic(MHD) boundary layer flow of an incompressible upper-convected Maxwell(UCM) fluid over a porous stretching surface.Similarity transformations are used to reduce the governing partial differential equations into a kind of nonlinear ordinary differential equations.The nonlinear problem is solved by using the successive Taylor series linearization method(STSLM).The computations for velocity components are carried out for the emerging parameters.The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem.     

Flow and heat transfer over a generalized stretching/shrinking wall problem—Exact solutions of the Navier-Stokes equations

In this paper, we investigate the steady momentum and heat transfer of a viscous fluid flow over a stretching/shrinking sheet. Exact solutions are presented for the Navier-Stokes equations. The new solutions provide a more general formulation including the linearly stretching and shrinking wall problems as well as the asymptotic suction velocity profiles over a moving plate. Interesting non-linear phenomena are observed in the current results including both exponentially decaying solution and algebraically decaying solution, multiple solutions with infinite number of solutions for the flow field, and velocity overshoot. The energy equation ignoring viscous dissipation is solved exactly and the effects of the mass transfer parameter, the Prandtl number, and the wall stretching/shrinking strength on the temperature profiles and wall heat flux are also presented and discussed. The exact solution of this general flow configuration is a rare case for the Navier-Stokes equation.     

Analytical solutions for some nonlinear evolution equations

IntroductionItiswell_knownthatmanyimportantdynamicsprocessescanbedescribedbyspecificnonlinearpartialdifferentialequations .Whenanonlinearpartialdifferentialequationisusedtodescribeaphysicalparameterthatshowssomekindsofpropagationoraggregationproperties,oneofthemostimportantphysicalmotivationsistosolvethepartialdifferentialequationwithacertaintypeoftravellingwavesolution .Inthepastseveraldecades,therehavebeenmanyattemptsinthisfieldbothbymathematiciansandphysicists[1]- [16 ],however,duetothecomp…     

局部与整体水压法在水下饱和边坡稳定分析中的应用

针对以考虑土条侧向条间力来提高瑞典条分法的计算精度方法, 作者基于物理学思想, 提出了局部与整体水压法以进行水下饱和边坡稳定分析, 并采用该法对瑞典条分法中一对相互矛盾的"真假李逵" 式进行了辨析. 指出李鬼式精度低的原因在于直接提取土条顶面和底面水压力沿着竖向的分量来建立平衡, 而忽略了其沿着水平方向分力的影响, 并证明了局部与整体水压法在提高有水边坡计算精度和便捷度时具有的普适性.     

分析力学中Santilli方法和Engels第一方法的意义和局限性

讨论变分法逆问题理论中的两种构造拉格朗日函数的基本方法:Santilli方法和Engels第一方法.(1)指出Santilli方法的理论意义在于直接用构造法证明自伴随微分方程能够从变分原理导出,即表示为欧勒-拉格朗日方程形式.(2)提出利用Santilli方法构造的结果,不是唯一的拉格朗日函数,而是一规范等效的拉格朗日函数族,为此修正了该方法.(3)指出在实际应用中Santilli方法的局限性,特别是对某些力学系统,可能因对参变量的定积分发散,而不能有效构造拉格朗日函数.(4)分析Engels第一方法的意义和优越性,同时指出这种方法存在与Santilli方法相似的局限性.(5)以两个力学系统为例,说明上述讨论的结论.     

复合材料中的渐近均匀化方法

本文将非均质弹性体的渐近均匀化方法应用于复合材料的宏观与细观分析之中。该方法基于平均化的思想，将复合材料视作由周期性的细观结构所构成，其场变量依赖于宏观和细观两个尺度的坐标变量而变化。通过建立位移和应力的渐近表达式，推导出关于周期性基元的细观平衡方程和细观本构关系，并与有限元数值方法相结合，得到材料的宏观等效性能和细观应力分布。对典型算例的分析，反映出该方法的有效性及准确性。     

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered finite volume/finite element method for solution of the two dimensional (2D) incompressible Navier-Stokes equations on unstructured grids.An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling.The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by joining the centroid of cells sharing the common vertex.For the temporal integration of the momentum equations,an implicit second-order scheme is utilized to enhance the computational stability and eliminate the time step limit due to the diffusion term.The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite element method (FEM).The momentum interpolation is used to damp out the spurious pressure wiggles.The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both velocity and pressure.The classic test cases,the lid-driven cavity flow,the skew cavity flow and the backward-facing step flow,show that numerical results are in good agreement with the published benchmark solutions.