GEL方法在一维双界面问题中的应用

用ghost-fluid Euler-Lagrange(GEL)耦合数值方法,对二维GEL计算程序(其中Euler流场计算采用以SCB格式编制的二阶计算程序,Lagrange域计算采用DEFEL二维动力有限元程序)进行扩展,将程序应用于多界面问题.算例为空气和水的一维双界面黎曼问题,获得了密度、速度和压力结果,并通过i...     

二维流体力学有限差分和有限元程序的耦合计算

根据压力和法向速度连续准则,将Euler方法为基础的MFPPM（Piecewise-Parabolic Method）程序和Lagrange方法为基础的DEFEL（2-D Finite ElementsCode,二维流体弹塑性动力有限元）程序进行耦合,发展了基于levelset的GEL（Ghost—fluid Euler-Lagrange）方法。该方法在处理大变形流场与小变形结构以及复杂流动与多物体相互作用等问题具有优越性。通过二维算例的计算结果与文献比较,检验了GEL方法和耦合程序的正确性,并对水下爆炸形成的流场对多物体作用过程进行了数值模拟。     

二维磁流体力学ALE数值模拟

介绍采用ALE方法进行的电磁内爆二维磁流体（MHD）力学模拟。二维MHD力学模拟的ALE方法分成拉格朗日、网格重分、对流输运三步。拉格朗日步采用文献[4]的时间分裂法,分成辐射冷却、热扩散、磁扩散、Lagrangian流动四步求解,并对Lagrangian流动步的方法进行了改进,消除了其拉氏计算的质量损失等。编制了二维计算程序并通过验证,获得的二维电磁内爆不稳定性发展与文献[2,5,11]吻合。     

多体系统Euler-Lagrange方程的最小二乘法与违约修正

针对受完整约束的多体系统动力学Euler-Lagrange方程,在其传统的直接增广算法和零空间方法基础上提出了当系统存在冗余约束情形下的最小二乘法.同时,对应于最小二乘法提出了改进的约束违约修正方法.本文还针对Euler-Lagrange方程的计算过程给出了相应Jacobi矩阵的QR分解和零空间连续正交基的算法.最后,以平行五连杆机构给出了数值结果并与部分现有方法进行比较.     

二维流体力学计算中的人为热流方法

采用Ｗ．Ｆ．Ｎｏｈ提出的人为粘性和人为热流方法，构造了计算人为热流的二维差分格式。应用Ｌａｇｒａｎｇｅ型流体力学计算程序，计算了一种二维轴对称高速碰撞模型。计算结果表明，人为热流方法可以很好地消除人为粘性在碰撞界面处的误差。     

采用NND方法计算三维喷管气流场

本文运用NND显式差分格式，计算了三维喷管气流场。气流场计算的基本方程为一般贴体坐标系下三维守恒型的欧拉方程。采用了时间分裂法和Steger-Warming矢通量分裂技术。在喷管内沿周向的每个由轴线和壁面构成的子午面上根据泊松方程生成贴体网格。本文运用三维程序计算了轴对称JPL喷管，同时与实验结果和前人采用轴对称二维程序所计算的结果做了对比。最后，本文还计算了三维矢量喷管，计算结果与现有的实验结果一致。通过轴对称JPL喷管和三维矢量喷管的计算考核，表明建立的算法和编写的计算程序是正确的。文中提出了采用子午面形式的贴体网格时奇性轴的处理方法。计算结果表明在喷管壁面处，马赫数与压强的计算结果与实验值吻合较好，而在喷管轴线处，只有当网格较密时，才能得出与实验结果接近的计算结果。     

基于Updated Lagrangian法的三维粘弹性大变形问题的有限元分析

本文基于Updated Lagrangian增量迭加方法,采用以现时Kirchhoff应力增量和现时Green应变增量表示的虚功方程和Kirchhoff应力张量-Green应变张量的积分本构关系,导出粘弹性大变形的虚功方程。依此采用二十节点三维立方等参数元编制相应的计算程序。三个算例结果与以往一维、二维的计算结果完全符合。     

移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

液体三维晃动特征问题的有限元数值计算方法

本文采用有限元方法数值求解任意刚性容器内液体三维晃动的固有频率和模态。通过建立液体晃动特征问题的泛函极值原理，编制了四面体等参单元有限元程序，计算了平放圆柱腔内三维液体晃动的特征频率，并将矩形容器、球腔、带“十”字隔板球形容器内的液体三维晃动计算结果与解析解、实验结果和二维有限元数值解进行了比较，程序的正确性得到了验证。     

以Kirchhoff应力张量-Green应变张量表示本构关系的粘弹性大变形平面问题的有限元法

本文基于Updated Lagrangian增量叠加方法,采用以现时Kirchhoff应力增量和现时Green应变增量表示的虚功方程和Kirchhoff应力-Green应变的积分本构关系,导出粘弹性大变形的变分方程.依此采用八节点二维等参数元编制相应的平面问题的计算程序,算例的计算结果与以往的工作完全符合.     

On the use of characteristic‐based split meshfree method for solving flow problems

This study presents characteristic‐based split (CBS) algorithm in the meshfree context. This algorithm is the extension of general CBS method which was initially introduced in finite element framework. In this work, the general equations of flow have been represented in the meshfree context. A new finite element and MFree code is developed for solving flow problems. This computational code is capable of solving both time‐dependent and steady‐state flow problems. Numerical simulation of some known benchmark flow problems has been studied. Computational results of MFree method have been compared to those of finite element method. The results obtained have been verified by known numerical, analytical and experimental data in the literature. A number of shape functions are used for field variable interpolation. The performance of each interpolation method is discussed. It is concluded that the MFree method is more accurate than FEM if the same numbers of nodes are used for each solver. Meshfree CBS algorithm is completely stable even at high Reynolds numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

Towards numerical simulations of supersonic liquid jets using ghost fluid method

A computational tool based on the ghost fluid method (GFM) is developed to study supersonic liquid jets involving strong shocks and contact discontinuities with high density ratios. The solver utilizes constrained reinitialization method and is capable of switching between the exact and approximate Riemann solvers to increase the robustness. The numerical methodology is validated through several benchmark test problems; these include one-dimensional multiphase shock tube problem, shock–bubble interaction, air cavity collapse in water, and underwater-explosion. A comparison between our results and numerical and experimental observations indicate that the developed solver performs well investigating these problems. The code is then used to simulate the emergence of a supersonic liquid jet into a quiescent gaseous medium, which is the very first time to be studied by a ghost fluid method. The results of simulations are in good agreement with the experimental investigations. Also some of the famous flow characteristics, like the propagation of pressure-waves from the liquid jet interface and dependence of the Mach cone structure on the inlet Mach number, are reproduced numerically. The numerical simulations conducted here suggest that the ghost fluid method is an affordable and reliable scheme to study complicated interfacial evolutions in complex multiphase systems such as supersonic liquid jets.     

A computationally efficient 3D finite‐volume scheme for violent liquid–gas sloshing

We describe a semi‐implicit volume‐of‐fluid free‐surface‐modelling methodology for flow problems involving violent free‐surface motion. For efficient computation, a hybrid‐unstructured edge‐based vertex‐centred finite volume discretisation is employed, while the solution methodology is entirely matrix free. Pressures are solved using a matrix‐free preconditioned generalised minimum residual algorithm and explicit time‐stepping is employed for the momentum and interface‐tracking equations. The high resolution artificial compressive (HiRAC) volume‐of‐fluid method is used for accurate capturing of the free surface in violent flow regimes while allowing natural applicability to hybrid‐unstructured meshes. The code is parallelised for solution on distributed‐memory architectures and evaluated against 2D and 3D benchmark problems. Good parallel scaling is demonstrated, with almost linear speed‐up down to 6000 cells per core. Finally, the code is applied to an industrial‐type problem involving resonant excitation of a fuel tank, and a comparison with experimental results is made in this violent sloshing regime. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of three solvers for viscoelastic fluid flow problems

The computational efficiency of three numerical schemes has been examined for the solution of a linearized system of equations resulting from the finite element discretization of a viscoelastic fluid flow problem. The first scheme is a modified frontal solver, which solves the linear system of equations directly. The other two, one based on a biconjugate gradient stabilized (BiCGStab) method and another based on a generalized minimal residual (GMRES) method, are iterative schemes. The stick-slip problem and the four-to-one contraction problem were analyzed and the viscoelastic fluid was assumed to obey the Oldroyd-B model. The two iterative schemes are superior to the direct scheme in terms of CPU time consumed and the BiCGStab scheme is even faster than the GMRES scheme. The range of convergence for both iterative schemes is compatible with that of the direct scheme.     

Modeling Three-Dimensional Groundwater Flows by the Body-Fitted Coordinate (BFC) Method: II. Free and Moving Boundary Problems

This paper presents a methodology and solution procedure of the time-dependent body-fitted coordinate (BFC) method for the analysis of transient, three-dimensional groundwater flow problems characterized by free and moving boundaries. The technique consists of numerical grid generation, time-dependent body-fitted coordinate transformation, and application of the finite difference method (FDM) to the transformed partial differential equations. Based on the time-dependent BFC method, a three-dimensional finite-difference computer code, BFC3DGW, was developed and used to solve two unconfined flow problems. The code was verified by comparing numerical results with analytical solutions for a steady-state seepage problem. In order to demonstrate capability of the method in dealing with flow problems with irregular and moving boundary surfaces, an unconfined well-flow problem was solved by the developed code. Difficulties associated with the free and moving irregular boundary have been successfully overcome by employing this method.     

Nonlinear convection of a viscoelastic fluid with inclined temperature gradient

The energy method is used to analyze the viscoelastic fluid convection problem in a thin horizontal layer, subjected to an applied inclined temperature gradient. The boundaries are considered to be rigid and perfectly conducting. Both linear and nonlinear stability analyses are carried out. The eigenvalue problem is solved by the Chebyshev Tau-QZ method and comparisons are reported between the results of the linear theory and energy stability theory.Received: 12 March 2004, Accepted: 19 April 2004, Published online: 17 September 2004PACS: 47.20 Ky, 47-27 Te, 83.60 Wc Correspondence to: P.N. Kaloni     

基于OpenMP的三维显式物质点法并行化研究

基于OpenMP技术开发了三维显式物质点并行程序MPM3DMP。为了避免节点更新阶段的数据竞争,采用区域分解法将背景网格分解为均匀的子域,每个线程负责一个子域的节点变量更新,然后将更新后的节点变量装配到整体。在质点更新阶段采用了循环分解方法进行并行。针对Taylor杆碰撞的三种计算模型,在双Intel Woodcrest 4核CPU服务器下进行了测试:粗模型在4核下加速比为3.82,在8核下为6.23,中模型在4核下加速比为3.79,在8核下加速比为6.23;细模型在4核下加速比为3.75,8核下加速比为6.26。因此,本文的并行程序具有较好的并行效率和可扩展性。     

Critical evaluation of CFD codes for interfacial simulation of bubble‐train flow in a narrow channel

Computational fluid dynamics (CFD) codes that are able to describe in detail the dynamic evolution of the deformable interface in gas–liquid or liquid–liquid flows may be a valuable tool to explore the potential of multi‐fluid flow in narrow channels for process intensification. In the present paper, a computational exercise for co‐current bubble‐train flow in a square vertical mini‐channel is performed to investigate the performance of well‐known CFD codes for this type of flows. The computations are based on the volume‐of‐fluid method (VOF) where the transport equation for the liquid volumetric fraction is solved either by the methods involving a geometrical reconstruction of the interface or by the methods that use higher‐order difference schemes instead. The codes contributing to the present code‐to‐code comparison are an in‐house code and the commercial CFD packages CFX, FLUENT and STAR‐CD. Results are presented for two basic cases. In the first one, the flow is driven by buoyancy only, while in the second case the flow is additionally forced by an external pressure gradient. The results of the code‐to‐code comparison show that only the VOF method with interface reconstruction leads to physically sound and consistent results, whereas the use of difference schemes for the volume fraction equation shows some deficiencies. Copyright © 2007 John Wiley & Sons, Ltd.     

A combined level set/ghost cell immersed boundary representation for floating body simulations

A six degrees of freedom (6DOF) algorithm is implemented in the open‐source CFD code REEF3D. The model solves the incompressible Navier–Stokes equations. Complex free surface dynamics are modeled with the level set method based on a two‐phase flow approach. The convection terms of the velocities and the level set method are treated with a high‐order weighted essentially non‐oscillatory discretization scheme. Together with the level set method for the free surface capturing, this algorithm can model the movement of rigid floating bodies and their interaction with the fluid. The 6DOF algorithm is implemented on a fixed grid. The solid‐fluid interface is represented with a combination of the level set method and ghost cell immersed boundary method. As a result, re‐meshing or overset grids are not necessary. The capability, accuracy, and numerical stability of the new algorithm is shown through benchmark applications for the fluid‐body interaction problem. Copyright © 2016 John Wiley & Sons, Ltd.