成型充填过程的ALE有限元模拟

在ALE框架中提出了一个用于成型充填过程有限元数值模拟的模型。应用ALE参考构形及ALE参考粒子速度描写充填过程中的熔体质量运动。摒弃了Hele-Shaw近似假定,因而所提出的模型能用于非薄壁型腔中高分子材料充填过程的数值模拟。应用基于时域分步算法的Taylor-Galerkin方法,对控制成型充填过程的守恒方程建立了弱形式。对移动自由面附近的充填材料区构造了网格生成算法与网格重划分方案。给出了在几种不同形状的典型腔体中充填过程的数值模拟结果,表明了所提出的ALE有限元模型模拟充填过程的有效性。     

一类格心型ALE有限体积格式方法

现在国内外流行的ALE有限体积格式基本上都基于交错网榕进行格式的离散.该类格武在进行重映时,速度、密度和能量需要分别进行重映计算,效率较低,而且由于速度在网格角点.而密度、能量在网格中心,重映时会出现动能和内能不协调现泉.本文在巳有格心型Lagrange有限体积格式研究的基础上,结合Abgrall R.等关于榕心型格式下的网格角点速度的计算方法,利用最小二乘法进行高阶插值多项式重构,构造了一类新的格心型的高精度Lagrangian有限体积格式,并结合有效的高精度ENO守恒重映方法,获得了一类格心型的高精度ALE有限体积格式.数值试验的结果说明本文的格式是有效的,高精度的,收敛的,并且避免了物理量的不协调现象.     

自由活塞压缩管ALE方法数值模拟

当前国际上实现高焓气体流动的实验手段之一是自由活塞驱动类脉冲设备,包括自由活塞激波风洞和自由活塞膨胀管.采用自由活塞压缩管作为激波风洞和膨胀管的驱动段时,其驱动能力在很大程度上决定了该类设备的性能.本文采用计算流体力学中任意拉格朗日——欧拉方法(arbitrary Lagrangian Eulerian)数值模拟了压缩管内部的自由活塞运动和气体流动特征.采用移动网格技术来适应活塞运动边界,耦合求解网格运动和气体流动过程,并通过双时间步长方法进行流体运动的时间积分.为了满足几何守恒律(geometric conservation law),对移动网格的法向矢量和表面面积计算进行了修正.不同时刻的活塞位置试验测量结果及欧拉方法预测结果,以及基于简单波理论获得的运动活塞底部气体压力、活塞速度与活塞位置都与当前的ALE方法十分一致.该工作为下一步数值模拟自由活塞激波风洞和自由活塞膨胀管中包括压缩管、激波管和喷管等不同部位的耦合流动提供了基础.      

基于ALE描述的压缩流动三维粘弹性数值模拟

塑料熔体压缩流动中在主流动平面及厚度方向均有速度变化,结合了剪切和拉伸两种流动特征。为了准确描述和模拟塑料熔体的压缩流动,本文基于粘弹性及ALE原理建立了熔体三维流动理论模型,构造了有限元求解的变分方程。为了避免整体求解计算量大、稳定性差的缺陷,提出了两重迭代解耦合算法分别求解耦合的连续方程、动量方程、本构方程、能量方程,开发了模拟程序。开展了等厚度板及变厚度板的注压成型实验及相应的数值模拟,结果表明：压缩过程中出现压力变化小于4.57%的平台现象;温度呈指数规律下降;塑料入口、流动末端第一法向应力差比平均值分别高1.73MPa、0.87MPa,变厚度区域第一法向应力差比平均值高1.16MPa。本文提出的理论模型和数值算法能够较好地表征压缩过程中熔体的压力、温度变化、应力演化。     

基于NMM的EFG方法及其裂纹扩展模拟

数值流形方法（numerucal manifold method,NMM）通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题. 通过用移动最小二乘插值（moving least squares interpolation,MLS）中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法（element free Galerkin,EFG）. 该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题. 建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹（kinked crack）的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性. 大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.     

断裂力学问题的杂交边界点方法

提出了一种求解断裂力学的新的边界类型无网格方法-杂交边界点法.以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性,求解时仅仅需要边界上离散点的信息.该文将杂交边界点方法应用到弹性断裂问题中,将移动最小二乘方法中的基函数扩充,能更好的模拟裂纹尖端应力场的奇异性,推导了求解断裂力学的杂交边界点法方程,与传统的元网格方法相比,文中方法具有后处理简单,计算精度高的优点.数值算例表明了该方法的稳定性和有效性.     

基于ALE方法的飞机水面降落过程

为了研究飞机水面降落过程中的动态响应规律,采用arbitrary Lagrangian-Eulerian(ALE)方法开展了飞机水面降落的动力学分析。研究了飞机水面降落过程的速度变化规律,结果表明,降落速度和飞行速度在飞机入水的初始阶段变化较快,随后变化幅度趋于平缓。分析了不同降落速度、飞行速度和降落仰角下的机身结构响应,得到了飞机水面降落时结构响应随时间的变化规律。机身结构应力在入水的初始阶段达到最大值,随后迅速下降,最后保持稳定。飞机结构的最大变形也出现在入水的初始阶段,随后迅速回复到初始状态。对比了降落速度、飞行速度和降落仰角对飞机结构响应的影响程度,结果表明降落速度对结构响应的影响程度最大,降落仰角次之,飞行速度的影响最小。 更多还原     

无网格方法中本质边界条件的处理

基于新的离散模定义,采用移动最小二乘近似方法,给出了场变量的近似表达形式;从约束变分原理出发,给出了无网格方法中新的本质边界条件的处理方案--罚函数法;对权函数、罚函数的选择对计算精度的影响进行了研究.数值计算结果表明该方法不仅简单合理,而且具有较高的精度和稳定性.     

基于拟合求解数字图像相关亚像素位移算法

数字图像相关方法是一种非接触式的变形场测量方法。虽然目前该算法能取得较高的计算精度,但是通常存在计算量较大、计算时间较长的问题。为了提高算法的计算效率,本文提出了一种基于移动最小二乘算法拟合整像素位移求解亚像素位移的数字图像相关方法。因为采用移动最小二乘算法有效利用了局部特征,所以可以采用相对简单的多项式拟合求解复杂的周期型余弦变形,同时计算结果具有较高的精度。计算结果的平均误差为0.0238pixel,标准误差为0.0791pixel,与常用的空间相关亚像素位移计算结果精度相当。实验结果验证了本文算法的正确性和有效性。通过与等计算精度的NR算法对比分析,结果表明本文算法在亚像素位移计算阶段具有更高的计算效率。计算效率至少提升了58.1倍,最多提升了437.8倍。     

弹性力学问题的局部 Petrov-Galerkin方法

提出了弹性力学平面问题的局部Petrov-Galerkin方法,这是一种真正的无网格方法。这种方法采和移动最小二乘近似函数作为试函数,并且采用移动最小二乘近似函数的权函数作为加权残值法加权函数;同时这种方法只包含中心在所考虑点处的规则局部区域上以及局部边界上的积分,所得系统矩阵是一个带状稀疏矩阵,该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。还计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模及其相对误差。所得计算结果证明：该方法是一种具有收敛快、精度高、简便有效的通用方法;在工程中具有广阔的应用前景。     

Prediction of extrudate swell in polymer melt extrusion using an Arbitrary Lagrangian Eulerian (ALE) based finite element method

Accurate prediction of extrudate (die) swell in polymer melt extrusion is important as this helps in appropriate die design for profile extrusion applications. Extrudate swell prediction has shown significant difficulties due to two key reasons. The first is the appropriate representation of the constitutive behavior of the polymer melt. The second is regarding the simulation of the free surface, which requires special techniques in the traditionally used Eulerian framework. In this paper we propose a method for simulation of extrudate swell using an Arbitrary Lagrangian Eulerian (ALE) technique based finite element formulation. The ALE technique provides advantages of both Lagrangian and Eulerian frameworks by allowing the computational mesh to move in an arbitrary manner, independent of the material motion. In the present method, a fractional-step ALE technique is employed in which the Lagrangian phase of material motion and convection arising out of mesh motion are decoupled. In the first step, the relevant flow and constitutive equations are solved in Lagrangian framework. The simpler representation of polymer constitutive equations in a Lagrangian framework avoids the difficulties associated with convective terms thereby resulting in a robust numerical formulation besides allowing for natural evolution of the free surface with the flow. In the second step, mesh is moved in ALE mode and the associated convection of the variables due to relative motion of the mesh is performed using a Godunov type scheme. While the mesh is fixed in space in the die region, the nodal points of the mesh on the extrudate free surface are allowed to move normal to flow direction with special rules to facilitate the simulation of swell. A differential exponential Phan Thien Tanner (PTT) model is used to represent the constitutive behavior of the melt. Using this method we simulate extrudate swell in planar and axisymmetric extrusion with abrupt contraction ahead of the die exit. This geometry allows the extrudate to have significant memory for shorter die lengths and acts as a good test for swell predictions. We demonstrate that our predictions of extrudate swell match well with reported experimental and numerical simulations.     

Modelling of two‐dimensional laminar flow using finite element method

In this paper, a Galerkin weighted residual finite element numerical solution method, with velocity material time derivative discretisation, is applied to solve for a classical fluid mechanics system of partial differential equations modelling two‐dimensional stationary incompressible Newtonian fluid flow. Classical examples of driven cavity laminar flow and laminar flow past a cylinder are presented. Numerical results are compared with data found in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

Dynamic simulation of free surfaces in capillaries with the finite element method

The mathematical formulation of the dynamics of free liquid surfaces including the effects of surface tension is governed by a non-linear system of elliptic differential equations. The major difficulty of getting unique closed solutions only in trivial cases is overcome by numerical methods. This paper considers transient simulations of liquid–gas menisci in vertical capillary tubes and gaps in the presence of gravity. Therefore the CFD code FIDAP 7.52 based on the Galerkin finite element method (FEM) is used. Calculations using the free surface model are presented for a variety of contact angles and cross-sections with experimental and theoretical verification. The liquid column oscillations are compared for numerical accuracy with a mechanical mathematical model, and the sensitivity with respect to the node density is investigated. The efficiency of the numerical treatment of geometric non-trivial problems is demonstrated by a prismatic capillary. Present restrictions limiting efficient transient simulations with irregularly shaped calculational domains are stated. © 1998 John Wiley & Sons, Ltd.     

A finite element method for two‐dimensional water‐wave problems

In this paper, the authors treat the free‐surface waves generated by a moving disturbance with a constant speed in water of finite and constant depth. Specifically, the case when the disturbance is moving with the critical speed is investigated. The water is assumed inviscid and its motion irrotational. The surface tension is neglected. It is well‐known that the linear theory breaks down when a disturbance is moving with the critical speed. As a remedy to overcome the invalid linear theory, approximate non‐linear theories have been applied with success in the past, i.e. Boussinesq and Korteweg de Vries equations, for example. In the present paper, the authors describe a finite element method applied to the non‐linear water‐wave problems in two dimensions. The present numerical method solves the exact non‐linear formulation in the scope of potential theory without any additional assumptions on the magnitude of the disturbances. The present numerical results are compared with those obtained by other approximate non‐linear theories. Also presented are the discussions on the validity of the existing approximate theories applied to two types of the disturbances, i.e. the bottom bump and the pressure patch on the free‐surface at the critical speed. Copyright © 1999 John Wiley & Sons, Ltd.     

Efficient simulation of free surface flows with discrete least-squares meshless method using a priori error estimator

In this article, a priori error estimate is employed to improve the efficiency of simulating free surface flows with discrete least-squares meshless (DLSM) method. DLSM is a fully least-squares approach in which both function approximation and the discretisation of the governing differential equations is carried out using a least-squares concept. The meshless shape functions are derived using the moving least-squares (MLS) method of function approximation. The discretised equations are obtained via a discrete least-squares method in which the sum of the squared residuals are minimised with respect to unknown nodal parameters. The governing equations of mass and momentum conservation are solved in a Lagrangian form using a pressure projection method. The proposed simulation strategy is composed of error estimation and a node moving refinement method. Since in free surface problems, the position of the free surface is of primary interest, a priori error estimate is used which automatically associates higher error to the nodes near the free surface. The node moving refinement method is used to construct a nodal configuration with dense nodal arrangement near the free surface. Four test problems namely dam break, evolution of a water bubble, solitary wave propagation and wave run-up on slope are investigated to test the ability and efficiency of the proposed efficient simulation method.     

近水面水下爆炸二维Level-set数值模拟

水下爆炸是一个多介质、动边界、瞬态非线性过程,捕捉介质之间交界面的变化是数值模拟的一个难点.本文采用Level-set方法来描述水下爆炸气体产物和流体交界面以及自由表面变化,应用TVD (Total Variation Diminishing) 计算技术求解Level-set流场,首先建立了近水面水下爆炸的二维数值方法,然后,将该方法推广到流场内有刚性结构的情况,考虑了冲击波和刚性结构物的相互作用.为了模拟无限边界流场问题,引进了Thompson建立的无反射边界条件.     

Thermal solution of cylindrical composite systems using meshless method

This paper presents the thermal solution of cylindrical composite systems using meshless element free Galerkin (EFG) method. The EFG method utilizes the moving least square approximants, which are constructed by using a weight function, a basis function and a set of non-constant coefficients to approximate the unknown function of temperature. Dirichlet (essential) boundary conditions have been enforced using Lagrange multiplier and penalty methods. Existing rational weight function has been modified and used in the present analysis. MATLAB codes have been developed to obtain the numerical solution. The EFG results have been obtained using cubicspline, quarticspline, Gaussian, quadratic, hyperbolic, exponential, rational and cosine weight functions for a model problem. The results obtained using different EFG weight functions are also compared with those obtained by finite element method. The effect of scaling and penalty parameters has also been studied in detail.     

A numerical study of projection algorithms in the finite element simulation of three‐dimensional viscous incompressible flow

In this work a comparative study of two versions of the projection algorithm used either for time integration or as an iterative method to solve the three‐dimensional incompressible Navier–Stokes equations is presented. It is also shown that these projection algorithms combined with the finite element method are particularly suited for the treatment of outflow boundary conditions in the context of external flows. This assertion is illustrated by means of some numerical examples in which five types of boundary conditions are compared. The scheme is applied to simulate the flow past a cylinder clamped on two fixed parallel solid walls. Comparison with experimental data available for this problem shows good agreement of the velocity and pressure fields, both computed with continuous piecewise linear elements. Copyright © 1999 John Wiley & Sons, Ltd.