Hamilton体系下环扇形域的Stokes流动问题

基于极坐标下Stokes流的基本方程,将径向坐标模拟为时间坐标,推导了Hamilton体系下Stokes流动问题的对偶方程,采用本征向量展开法对环扇形域Stokes流动问题进行了分析,并给出了相应的实际算例,其结果说明了本文方法的有效性。     

平行平板间Couette流起动过程的运算微积解

 求解平行平板间Couette流的起动过程通常使用分离变量法，得到的速度分布在平板间距 趋于无穷时难以逼近Stokes第一问题的解，而两者在物理上却是完全一致的. 本文利用运 算微积法和Heaviside算子进行求解，不但解决了这一问题，而且具有运算简单的 特点.     

固支矩形Mindlin板弯曲问题的辛求解体系

本文利用二类变量广义变分原理推出了Mindlin板弯曲问题的Hamilton体系,利用辛几何方法对全状态向量进行分离变量,得到相应的横向本征问题,在求出其本征值后,按本征函数展开法导出了原问题的辛本征通解。给出了一个承受集中载荷的四边固支矩形薄板的算例,按本文求解体系得到的解与经典解吻合较好。本文直接从Mindlin板弯曲问题出发,在其Hamilton体系内使用辛几何方法给出了的一套新的求解体系,突破了传统解法的局限性,具有一般性及较高的理论推广价值。     

四边固支矩形弹性薄板的精确解析解

利用辛几何方法本文推导出了四边固支矩形弹性薄板弯曲问题的精确解析解.由于在求解过程中不需要事先人为的选取挠度函数,而是从弹性薄板的基本方程出发,首先将矩形薄板弯曲问题表示成Hamilton正则方程,然后利用分离变量和本征函数展开的方法求出可以完全满足四边固支边界条件的精确解析解.本文中所采用的方法突破了传统的半逆法的限制,使得问题的求解更加合理化.文中还给出了计算实例来证明推导结果的正确性.     

矩形梁受分布荷载的辛解答

在弹性力学平面直角坐标辛体系中,采用分离变量法,放弃齐次边界条件,得到了矩形梁侧边受幂函数形式分布荷载问题的辛解答,给出了这类问题在辛体系中的一般解法,分别对矩形梁受法向和切向分布荷载的问题进行了求解,显示了此方法的有效性.辛解法采用对偶的二类变量进行求解,可同时给出位移和应力;由于辛解法能较好地处理各种边界条件,因此不仅能求解静定问题,也能直接求解静不定问题.     

横向力作用下悬臂梁固定端应力分布研究

在弹性力学Hamilton体系中,利用解析法,考虑圣维南原理所覆盖的解,对横向力作用下悬臂梁固定端应力分布问题进行研究,并对计算结果进行分析。研究结果表明,辛体系解析法采用对偶的二类变量求解,能很好地处理各种复杂边界条件,并且对此类问题的分析具有优越性,计算精度较高。该方法对其他边界问题的研究也具有指导意义。     

矩形悬臂薄板的解析解

首先把弹性薄板弯曲问题的控制方程表示成为Hamilton正则方程,然后利用辛几何方法对全状态相变量进行分离变量,求出其本征值后,再按本征函数展开的方法求出矩形悬臂薄板的解析解。由于在求解过程中不需要事先人为地选取挠度函数,而是从薄板弯曲的基本方程出发,直接利用数学的方法求出可以满足其边界条件的这类问题的解析解,使得问题的求解更加理论化和合理化。文中的最后还给出了计算实例来验证本文所采用的方法以及所推导出的公式的正确性。     

广义Stokes方程的完全边界积分表示式及其在求解N-S方程中的应用

为了分解N-S方程组各变量相互偶合,本文采用Peaceman-Rachford算子分裂法,将时间相依的N-S方程组分解成不存在上述偶合特性的线性和非线性的子问题。线性子问题具有广义Stokes方程类型。本文采用多重互易法,即采用多阶拉普拉斯算子基本解逐步变换,将其解表示成完全边界积分形式,从而使问题的计算维数降低一维。广义Stokes方程的算例以及二维圆柱在剪切流中的Stokes绕流解,都表明多重互易算法具有高效特点,而且后者与文[3]解析解吻合得非常好。     

变厚度圆柱壳的轴对称自由振动

本文借助于状态空间法研究变厚度圆柱壳的轴对称自由振动问题.引进状态变量,建立状态方程,用状态空间法求解具有任意边界条件和厚度变化形式的圆柱壳的固有频率和振型。     

用分离变量法求解Stokes第一问题

探讨如何用较为通用的分离变量法求解Stokes第一问题。     

Weak imposition of boundary conditions for the Navier–Stokes equations by a penalty method

We prove convergence of the finite element method for the Navier–Stokes equations in which the no‐slip condition and no‐penetration condition on the flow boundary are imposed via a penalty method. This approach has been previously studied for the Stokes problem by Liakos (Weak imposition of boundary conditions in the Stokes problem. Ph.D. Thesis, University of Pittsburgh, 1999). Since, in most realistic applications, inertial effects dominate, it is crucial to extend the validity of the method to the nonlinear Navier–Stokes case. This report includes the analysis of this extension, as well as numerical results validating their analytical counterparts. Specifically, we show that optimal order of convergence can be achieved if the computational boundary follows the real flow boundary exactly. Copyright © 2008 John Wiley & Sons, Ltd.     

An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Higher order finite element methods and multigrid solvers in a benchmark problem for the 3D Navier–Stokes equations

This paper presents a numerical study of the 3D flow around a cylinder which was defined as a benchmark problem for the steady state Navier–Stokes equations within the DFG high‐priority research program flow simulation with high‐performance computers by Schafer and Turek (Vol. 52, Vieweg: Braunschweig, 1996). The first part of the study is a comparison of several finite element discretizations with respect to the accuracy of the computed benchmark parameters. It turns out that boundary fitted higher order finite element methods are in general most accurate. Our numerical study improves the hitherto existing reference values for the benchmark parameters considerably. The second part of the study deals with efficient and robust solvers for the discrete saddle point problems. All considered solvers are based on coupled multigrid methods. The flexible GMRES method with a multiple discretization multigrid method proves to be the best solver. Copyright © 2002 John Wiley & Sons, Ltd.     

谱元法和高阶时间分裂法求解方腔顶盖驱动流

详细推导了谱元方法的具体计算公式和时间分裂法的具体计算过程 ;对一般的时间分裂法进行了改进 ,即对非线性步分别用 3阶 Adams-Bashforth方法和 4阶显式 Runge-Kutta法 ,粘性步采用 3阶隐式 Adams-Moulton形式 ,提高了时间方向的离散精度 ,同时还改进了压力边界条件 ,采用 3阶的压力边界条件 ;利用改进的时间分裂方法分解不可压缩 Navier-Stokes方程 ,并结合谱元法计算了移动顶盖方腔驱动流 ,提高了方法可以计算的 Re数 ,缩短了达到收敛的时间 ,并将结果与基准解进行比较 ;分析了移动顶盖方腔驱动流中 Re数对流场分布的影响。     

Spectral/hp penalty least‐squares finite element formulation for unsteady incompressible flows

In this paper, we present spectral/hp penalty least‐squares finite element formulation for the numerical solution of unsteady incompressible Navier–Stokes equations. Pressure is eliminated from Navier–Stokes equations using penalty method, and finite element model is developed in terms of velocity, vorticity and dilatation. High‐order element expansions are used to construct discrete form. Unlike other penalty finite element formulations, equal‐order Gauss integration is used for both viscous and penalty terms of the coefficient matrix. For time integration, space–time decoupled schemes are implemented. Second‐order accuracy of the time integration scheme is established using the method of manufactured solution. Numerical results are presented for impulsively started lid‐driven cavity flow at Reynolds number of 5000 and transient flow over a backward‐facing step. The effect of penalty parameter on the accuracy is investigated thoroughly in this paper and results are presented for a range of penalty parameter. Present formulation produces very accurate results for even very low penalty parameters (10–50). Copyright © 2008 John Wiley & Sons, Ltd.     

Finite element implementation of boundary conditions for the pressure Poisson equation of incompressible flow

In this paper we address the problem of the implementation of boundary conditions for the derived pressure Poisson equation of incompressible flow. It is shown that the direct Galerkin finite element formulation of the pressure Poisson equation automatically satisfies the inhomogeneous Neumann boundary conditions, thus avoiding the difficulty in specifying boundary conditions for pressure. This ensures that only physically meaningful pressure boundary conditions consistent with the Navier-Stokes equations are imposed. Since second derivatives appear in this formulation, the conforming finite element method requires C¹ continuity. However, for many problems of practical interest (i.e. high Reynolds numbers) the second derivatives need not be included, thus allowing the use of more conventional C⁰ elements. Numerical results using this approach for a wall-driven contained flow within a square cavity verify the validity of the approach. Although the results were obtained for a two-dimensional problem using the p-version of the finite element method, the approach presented here is general and remains valid for the conventional h-version as well as three-dimensional problems.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

共轭梯度法求解非线性多宗量稳态传热反问题

应用共轭梯度法求解非线性多宗量稳态热传导反问题。采用八节点的等参单元在空间上进行离散,建立了便于敏度分析的非线性正演和反演的有限元模型,可直接求导进行敏度分析。给出了相关的数值验证,对测量误差及测点数目的影响作了初步探讨,结果表明,采用的算法能够对非线性稳态热传导中导热系数和边界条件联合反问题进行有效的求解,并具有较高精度。     

Boundary element analysis of driven cavity flow for low and moderate Reynolds numbers

A boundary element method for steady two‐dimensional low‐to‐moderate‐Reynolds number flows of incompressible fluids, using primitive variables, is presented. The velocity gradients in the Navier–Stokes equations are evaluated using the alternatives of upwind and central finite difference approximations, and derivatives of finite element shape functions. A direct iterative scheme is used to cope with the non‐linear character of the integral equations. In order to achieve convergence, an underrelaxation technique is employed at relatively high Reynolds numbers. Driven cavity flow in a square domain is considered to validate the proposed method by comparison with other published data. Copyright © 2001 John Wiley & Sons, Ltd.     

Potential flow around obstacles using the scaled boundary finite‐element method

The scaled boundary finite‐element method is a novel semi‐analytical technique, combining the advantages of the finite element and the boundary element methods with unique properties of its own. The method works by weakening the governing differential equations in one co‐ordinate direction through the introduction of shape functions, then solving the weakened equations analytically in the other (radial) co‐ordinate direction. These co‐ordinate directions are defined by the geometry of the domain and a scaling centre. The method can be employed for both bounded and unbounded domains. This paper applies the method to problems of potential flow around streamlined and bluff obstacles in an infinite domain. The method is derived using a weighted residual approach and extended to include the necessary velocity boundary conditions at infinity. The ability of the method to model unbounded problems is demonstrated, together with its ability to model singular points in the near field in the case of bluff obstacles. Flow fields around circular and square cylinders are computed, graphically illustrating the accuracy of the technique, and two further practical examples are also presented. Comparisons are made with boundary element and finite difference solutions. Copyright © 2003 John Wiley & Sons, Ltd.