Reduced augmented Lagrangian bi-conjugate gradient method for impact-contact problems

To avoid the numerical oscillation of the penalty method and non-compatibility with explicit operators of conventional Lagrange multiplier methods used in transient contact problems to enforce surface contact conditions, a new approach to enforcing surface contact constraints for the transient nonlinear finite element problems, referred to as “the reduced augmented Lagrangian bi-conjugate gradient method (ALCG)”, is developed in this paper. Based on the nonlinear constrained optimization theory and is compatible with the explicit time integration scheme, this approach can also be used in implicit scheme naturally. The new surface contact constraint method presented has significant advantages over the widely adopted penalty function methods and the conventional Lagrangian multiplier methods. The surface contact constraints are satisfied more accurately for each step by the algorithm, so the oscillation of numerical solution for the explicit scheme is depressed. Through the development of new iteration strategy for solving nonlinear equations, ALCG method improves the computational efficiency greatly. Project supported by State Education Commission Doctoral Foundation and Natural Science Foundation of Liaoning Province.     

UNIFIED COMPUTATION OF FLOW WITH COMPRESSIBLE AND INCOMPRESSIBLE FLUID BASED ON ROE＇S SCHEME

A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe＇s type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe＇s scheme.     

一种基于Associated Hermite正交函数求解对流扩散方程的算法

提出了一种基于AH(Associated Hermite)正交基函数求解对流扩散方程的无条件稳定算法。该算法将方程的时间项通过Hermite多项式作为正交基函数进行展开,利用Galerkin方法消除时间变量项,从而导出有限维AH域隐式差分方程,突破了传统显式差分格式稳定性条件的限制,最后通过对AH域展开系数的求解得到该对流扩散方程的数值解。在数值算例中,将该算法与传统显示差分法和交替方向隐式差分法进行对比分析,数值计算结果表明,算法无条件稳定且其计算精度与时间步长无关,对于具有精细结构的对流换热问题,该算法具有明显的效率优势,且保持了较高的精度。     

Galerkin meshless methods based on partition of unity quadrature

Introduction Meshlessmethodsarenewmethodsofnumericalcomputationwhichhavebeendeveloped rapidlyinrecentyears.Inthesemethods,onlynodesareneeded,meshinformationistotally unnecessary.Thiscanavoidorpartlyavoidthedifficultyofmeshgeneration.Duetohigh accuracyandstability,Galerkinmeshlessmethodsareappliedbroadly,butitisunavoidable tocomputetheintegrationoverthewholephysicaldomaininGalerkinweakform,whichisa greatchallengeforGalerkinmeshlessmethodsbecauseoftheabsenceofmesh.TocarryouttheintegrationinGal…     

Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder atRe=200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number. The project supported by the State Key Fundamental Research Project of “Large Scale Scientific Computation Research” (G199903281)     

精细时程积分法及其数值衍生格式应用评估

旋翼气动弹性耦合动力学方程本质上是一组刚性比较大的非线性偏微分方程。在有限元结构离散后，可改写为非齐次微分方程组，其中非齐次项是桨叶运动量（位移与速度）和气动载荷的函数。针对这类方程，本文尝试引入精细积分法及其衍生格式，借助数值方法计算Duhamel积分项。从积分精度与数值稳定性方面比较研究具有代表性的精细库塔法和高精度直接积分法。结合隐式积分算法，评估精细积分法应用于旋翼动力学方程的可行性。算例表明，精细积分法对矩形直桨叶动力学方程具有足够的求解精度。     

An iterative stabilized fractional step algorithm for numerical solution of incompressible N–S equations

Stabilized fractional step algorithm has been widely employed for numerical solution of incompressible Navier–Stokes equations. However, smaller time step sizes are required to use for existing explicit and semi‐implicit versions of the algorithm due to their fully or partially explicit nature particularly for highly viscous flow problems. The purpose of this paper is to present two modified versions of the fractional step algorithm using characteristic based split and Taylor–Galerkin like based split. The proposed modified versions of the algorithm are based on introducing an iterative procedure into the algorithm and allow much larger time step sizes than those required to the preceding ones. A numerical study of stability at acceptable convergence rate and accuracy as well as capability in circumventing the restriction imposed by the LBB condition for the proposed iterative versions of the algorithm is carried out with the plane Poisseuille flow problem under different Reynolds numbers ranging from low to high viscosities. Numerical experiments in the plane Poisseuille flow and the lid‐driven cavity flow problems demonstrate the improved performance of the proposed versions of the algorithm, which are further applied to numerical simulation of the polymer injection moulding process. Copyright © 2005 John Wiley & Sons, Ltd.     

A high-order nodal discontinuous Galerkin method to solve preconditioned multiphase Euler/Navier-Stokes equations for inviscid/viscous cavitating flows

In this study, a high-order accurate numerical method is applied and examined for the simulation of the inviscid/viscous cavitating flows by solving the preconditioned multiphase Euler/Navier-Stokes equations on triangle elements. The formulation used here is based on the homogeneous equilibrium model considering the continuity and momentum equations together with the transport equation for the vapor phase with applying appropriate mass transfer terms for calculating the evaporation/condensation of the liquid/vapor phase. The spatial derivative terms in the resulting system of equations are discretized by the nodal discontinuous Galerkin method (NDGM) and an implicit dual-time stepping method is used for the time integration. An artificial viscosity approach is implemented and assessed for capturing the steep discontinuities in the interface between the two phases. The accuracy and robustness of the proposed method in solving the preconditioned multiphase Euler/Navier-Stokes equations are examined by the simulation of different two-dimensional and axisymmetric cavitating flows. A sensitivity study is also performed to examine the effects of different numerical parameters on the accuracy and performance of the solution of the NDGM. Indications are that the solution methodology proposed and applied here is based on the NDGM with the implicit dual-time stepping method and the artificial viscosity approach is accurate and robust for the simulation of the inviscid and viscous cavitating flows.     

Local discontinuous galerkin method for radial porous flow with dispersion and adsorption

Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu, according to the form of the generalized convection-diffusion equations which model the radial porous flow with dispersion and adsorption, a local discontinuous Galerkin method for radial porous flow with dispersion and adsorption was developed, a high order accurary new scheme for radial porous flow is obtained. The presented method was applied to the numerical tests of two cases of radial porous, i. e. , the convection-dispersion flow and the convection-dispersion-adsorption flow, the corresponding parts of the numerical results are in good agreement with the published solutions, so the presented method is reliable. Reckoning of the computational cost also shows that the method is practicable.     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

广义对流扩散方程的隐式特征线混合元方法

提出了一个广义对流扩散方程的混合有限元方法，方程的基本变量及其空间梯度和流量在单 元内均作为独立变量分别插值. 基于胡海昌-Washizu三变量广义变分原理结合特征线法给 出了控制方程的单元弱形式. 混合元方法采用基于一点积分方案并结合可以滤掉虚假的 数值震荡的隐式特征线法. 数值结果证明了所提出的方法可以提供和四点积分同样的数 值计算结果，并能够提高计算效率.     

Nonlinear vibrations of viscoelastic cylindrical shells taking into account shear deformation and rotatory inertia

The vibration problem of a viscoelastic cylindrical shell is studied in a geometrically nonlinear formulation using the refined Timoshenko theory. The problem is solved by the Bubnov–Galerkin procedure combined with a numerical method based on quadrature formulas. The choice of relaxation kernels is substantiated for solving dynamic problems of viscoelastic systems. The numerical convergence of the Bubnov–Galerkin procedure is examined. The effect of viscoelastic properties of the material on the response of the cylindrical shell is discussed. The results obtained by various theories are compared.     

率相关本构方程积分新算法

提出一种积分率相关本构方程的隐式积分新算法,引入0～1范围内的缩放因子λ对本构方程进行间接求解,可以避免直接求解等效塑性应变或等效塑性应变率时,由于其数值过大或过小而造成的收敛困难或收敛失败,实现对率相关本构方程的快速准确求解.以B-P统一本构方程及双曲正弦本构方程为例,验证了算法的可行性.结果表明,新算法对于准静态变形条件下的无硬化本构方程也可以得出准确的解.     

An adaptive algorithm of precise integration for transient analysis

This paper presents an improved precise integration algorithm for transient analysis of heat transfer and some other problems. The original precise integration method is improved by means of the inverse accuracy analysis so that the parameterN, which has been taken as a constant and an independent parameter without consideration of the problems in the original method, can be generated automatically by the algorithm itself. Thus, the improved algorithm is adaptive and the accuracy of the algorithm is not dependent on the length of the time step in the integration process. It is shown that the numerical results obtained by the method proposed are more accurate than those obtained by the conventional time integration methods such as the difference method and others. Four examples are given to demonstrate the validity, accuracy and efficiency of the new method. Project supported by the National Natural Science Foundation of China (No. 19872016, 19872017), the National Key Basic Research Special Foundation (G1999032805) and the Foundation for University Key Teachers by the Ministry of Education of China.     

High-order time integration applied to metal powder plasticity

The present article treats two objectives. In the first investigation attention is focused on the application of time-adaptive finite elements formulated on the basis of a high-order time integration procedure on a constitutive model for compressible finite strain viscoplasticity for metal powder. In this connection, it has to be emphasized that the integration procedure is not only applied to the evolution equations on Gauss-point level but on the total system of differential–algebraic equations resulting from the application of the vertical line method on the quasi-static finite element equations. The specific application emerges from the field of metal powder compaction. Particular studies are carried out using stiffly accurate, diagonally implicit Runge–Kutta methods in combination with the Multilevel-Newton algorithm for solving the DAE-system. In this respect, the effort vs. accuracy behavior is investigated which is also related to order reduction known in elastoplasticity. The second topic treats the local stress algorithm for taking into account the yield function based finite strain viscoplasticity model, where the classical Newton–Raphson method fails. This is the reason why most constitutive models of powder materials are implemented into explicit finite element codes. Thus, the proposed investigations compare different methods in view of a stable and efficient integration process in implicit finite element formulations.     

注塑成型三维流动的数值模拟

建立了非等温、粘性、不可压缩、非牛顿流体流动的控制方程。为了避免同时求解耦合的压力场、速度场,本文通过修改Galerkin方法的变分方程,导出了关于压力场的拟Poisson方程,用迭代法独立地求解连续性方程、动量方程,并进行速度一粘度迭代求出最终的压力场、速度场。由于直接使用Galerkin方法求解能量方程容易引起温度场的振荡,本文采用隐式格式及“上风”法离散能量方程,用超松驰迭代法求解温度场的代数方程组。比较了模拟结果与等温管道流动的解析解及法兰的实际注射结果,算例表明本文方法可以预测注射成型流动过程中的一些重要特征。与传统Galerkin方法相比,本文方法可以减少内存,提高数值方法的稳定性。     

Direct numerical simulation of disturbances generated by periodic suction-blowing in a hypersonic boundary layer

A numerical algorithm and code are developed and applied to direct numerical simulation (DNS) of unsteady two-dimensional flow fields relevant to stability of the hypersonic boundary layer. An implicit second-order finite-volume technique is used for solving the compressible Navier–Stokes equations. Numerical simulation of disturbances generated by a periodic suction-blowing on a flat plate is performed at free-stream Mach number 6. For small forcing amplitudes, the second-mode growth rates predicted by DNS agree well with the growth rates resulted from the linear stability theory (LST) including nonparallel effects. This shows that numerical method allows for simulation of unstable processes despite its dissipative features. Calculations at large forcing amplitudes illustrate nonlinear dynamics of the disturbance flow field. DNS predicts a nonlinear saturation of fundamental harmonic and rapid growth of higher harmonics. These results are consistent with the experimental data of Stetson and Kimmel obtained on a sharp cone at the free-stream Mach number 8.     

基于Fourier级数的时变周期系数Riccati微分方程精细积分

结合Fourier级数展开方法,本文提出了基于精细积分的时变周期系数Riccati微分方程求解高效算法.首先,利用Fourier级数展开方法将周期系统表示成三角级数形式,在一个积分步内使用精细积分方法得到对应Hamilton系统状态转移矩阵的表达式.然后,通过Riccati变换的方法,得到含有状态转移矩阵的时变周期系数Riccati微分方程解的递推格式.本文方法充分利用了方程本身的周期性特点,文中的数值算例表明算法具有计算效率高、结果可靠等优势.     

Lagrangian finite element analysis applied to viscous free surface fluid flow

A new Lagrangian finite element formulation is presented for time-dependent incompressible free surface fluid flow problems described by the Navier-Stokes equations. The partial differential equations describing the continuum motion of the fluid are discretized using a Galerkin procedure in conjunction with the finite element approximation. Triangular finite elements are used to represent the dependent variables of the problem. An effective time integration procedure is introduced and provides a viable computational method for solving problems with equality of representation of the pressure and velocity fields. Its success has been attributed to the strict enforcement of the continuity constraint at every stage of the iterative process. The capabilities of the analysis procedure and the computer programs are demonstrated through the solution of several problems in viscous free surface fluid flow. Comparisons of results are presented with previous theoretical, numerical and experimental results.     

二维非局部线弹性平面问题的辛分析

将二维非局部线弹性理论引入到Hamilton体系下,基于变分原理推导得出了二维线弹性理论的对偶方程和相应的边界条件.在分析验证对偶方程的准确性的基础上,该套方法被应用于二维弹性平面波问题的求解.将精细积分与扩展的W-W算法相结合在Hamilton体系下建立了求解平面Rayleigh波的数值算法.从推导到计算的保辛性确保了辛体系非局部理论与算法的准确性.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别,并进一步指出了该套算法的适用性和优势所在.