Prismatic versus tetrahedral elements in three-dimensional finite element analyses of subsurface systems

A mesh of prismatic or tetrahedral elements automatically generated from an initial triangular grid is used to integrate 3-D flow equation in space. Many numerical comparisons between these two models have been performed. The results show that integration with tetrahedrons is as accurate as integration with prisms but much more efficient. The CPU time of solution with prismatic elements is about three times greater than that required employing tetrahedral elements.     

广义二阶流体管内轴向流动

在流体的本构关系中引入分数阶导数运算，对于介于粘性与弹性之间的流体的描述更具有合理性。本文将这种关系引入二阶流体，研究其管内轴向流动。我们先求出了１／２阶导数的解析解，用以验证Ｌａｐｌａｃｅ数值反演的ＣＲＵＭＰ方法的有效性。然后用ＣＲＵＭＰ法分析二阶流体管内轴向流动的特征。分析表明粘弹性特征越明显的流体，其速度与应力对分数导数的阶数越具有敏感性。     

自由表面与粘性尾迹的相互作用

解析地研究了无限深不可压粘性流体中运动物体产生层流尾迹与自由表面波的相互作用．以定常的Oseen方程模拟受扰流动,对于小振幅自由表面波则采用线性化的运动学和动力学边界条件．在数学描述上,运动物体以Oseen极子模拟,受扰流场分解成表述粘性尾迹的无界奇异Oseen流和描述自由面效应的有界正则Oseen流之和．通过积分变换法,得到自由表面波的精确解．借助Lighthill的两步格式,导出了自由面波高带有附加校正项的渐近解．所得对称解显示了波动的振幅因粘性和潜深的存在而呈指数衰减．     

Discrete groups generated by reflections in the faces of symplicial prisms in Lobachevskian spaces

The paper provides the Coxeter schemes of all the discrete groups generated by reflections in the faces of simplicial prisms the combinatory type of which is the product of a simplex and a segment in Lobachevskian spaces.     

A composite Level Set and Extended-Domain-Eigenfunction Method for simulating 2D Stokes flow involving a free surface

In this paper, the Extended-Domain-Eigenfunction-Method (EDEM) is combined with the Level Set Method in a composite numerical scheme for simulating a moving boundary problem. The liquid velocity is obtained by formulating the problem in terms of the EDEM methodology and solved using a least square approach. The propagation of the free surface is effected by a narrow band Level Set Method. The two methods both pass information to each other via a bridging process, which allows the position of the interface to be updated. The numerical scheme is applied to a series of problems involving a gas bubble submerged in a viscous liquid moving subject to both an externally generated flow and the influence of surface tension.     

四段式有限元法分析热耦合的流体-固体相互作用问题

给出了一种流(体)-热-结构综合的分析方法,固体中的热传导耦合了粘性流体中的热对流,因而在固体中产生热应力.应用四段式有限元法和流线逆风Petrov-Galerkin法分析热粘性流动,应用Galerkin法分析固体中的热传导和热应力.应用二阶半隐式Crank-Nicolson格式对时间积分,提高了非线性方程线性化后的计算效率.为了简化所有有限元公式,采用3节点的三角形单元,对所有的变量:流体的速度分量、压力、温度和固体的位移,使用同阶次的插值函数.这样做的主要优点是,使流体-固体介面处的热传导连接成一体.数个测试问题的结果表明,这种有限元法是有效的,且能加深对流(体)-热-结构相互作用现象的理解.     

Operator splitting for the numerical solution of free surface flow at low capillary numbers

Free surface flows are pervasive in engineering and biomedical applications. In many interesting cases—particularly when small length scales are involved—surface forces (capillarity) dominate the flow dynamics. In these cases, computing the flow together with the shape of the surfaces, requires specialized solution techniques. This article investigates the capabilities of an operator splitting/finite elements method at handling accurately incompressible viscous flow with free surfaces at low capillary numbers. The test case of flow in the downstream section of a slot coater is used for three reasons: (1) it is an established benchmark; (2) it represents an idealized, yet industrially relevant flow; (3) high-fidelity results obtained with monolithic algorithms are available in literature. The flow and free surface shape attained with the new operator splitting scheme agree very satisfactorily with the results obtained with monolithic solvers. Because of its inherent computational simplicity, the new operator splitting scheme is attractive for large-scale simulations, three-dimensional flows, and flows of complex fluids.     

A vertex‐centered and positivity‐preserving scheme for anisotropic diffusion equations on general polyhedral meshes

We propose a new nonlinear positivity‐preserving finite volume scheme for anisotropic diffusion problems on general polyhedral meshes with possibly nonplanar faces. The scheme is a vertex‐centered one where the edge‐centered, face‐centered, and cell‐centered unknowns are treated as auxiliary ones that can be computed by simple second‐order and positivity‐preserving interpolation algorithms. Different from most existing positivity‐preserving schemes, the presented scheme is based on a special nonlinear two‐point flux approximation that has a fixed stencil and does not require the convex decomposition of the co‐normal. More interesting is that the flux discretization is actually performed on a fixed tetrahedral subcell of the primary cell, which makes the scheme very easy to be implemented on polyhedral meshes with star‐shaped cells. Moreover, it is suitable for polyhedral meshes with nonplanar faces, and it does not suffer the so‐called numerical heat‐barrier issue. The truncation error is analyzed rigorously, while the Picard method and its Anderson acceleration are used for the solution of the resulting nonlinear system. Numerical experiments are also provided to demonstrate the second‐order accuracy and well positivity of the numerical solution for heterogeneous and anisotropic diffusion problems on severely distorted grids.     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

The use of Hestenes' method of multipliers to resolve dual gaps in engineering system optimization

For nonconvex problems, the saddle point equivalence of the Lagrangian approach need not hold. The nonexistence of a saddle point causes the generation of a dual gap at the solution point, and the Lagrangian approach then fails to give the solution to the original problem. Unfortunately, dual gaps are a fairly common phenomenon for engineering system design problems. Methods which are available to resolve the dual gaps destroy the separability of separable systems. The present work employs the method of multipliers by Hestenes to resolve the dual gaps of engineering system design problems; it then develops an algorithmic procedure which preserves the separability characteristics of the system. The theoretical foundations of the proposed algorithm are developed, and examples are provided to clarify the approach taken.     

DTM-BF方法和可渗透收缩壁面上带滑移速度的非稳态磁流体力学流动

研究了运动的粘性导电流体中可渗透收缩壁面上非稳态磁流体边界层流动,利用解析和数值方法对问题进行了研究,并考虑了壁面速度滑移的影响．提出了一个新的解析方法（DTM-BF）,并将其应用于求解带有无穷远边界条件的非线性控制方程的近似解析解．对所有的解析结果和数值结果进行了对比,结果显示两者非常吻合,从而证明了DTM-BF方法的有效性．另外,对不同的参数,得到了控制方程双解和单解的存在范围．最后,分别讨论了滑移参数、非稳态参数、磁场参数、抽吸/喷注参数和速度比例参数对壁面摩擦、唯一解速度分布和双解速度分布的影响．     

椭圆型方程四面体线元的超逼近与外推

重新讨论了三角线元的积分恒等式,使之适用于三维区域的拟一致四面体元,借此证明了椭圆型方程有限元解梯度有超逼近现象,函数值Richardson外推可以提高精度.     

Third-order accurate finite-volume method on a triangular grid

A third-order accurate finite-volume method on triangular grids is proposed for the numerical solution of conservation law systems. The method is described in detail as applied to the advection equation. The accuracy of its numerical solutions obtained on grids of various degrees of detail is compared. The viscous flow over a plate and the unsteady flow around a cylinder are solved. For comparison purposes, the latter problem is also solved using third- and fifth-order accurate compact approximations.     

Computational fluid dynamics applied to the aerodynamic design of a land‐based supersonic vehicle

The BLOODHOUND SSC project was publicly announced in October 2008, with a primary engineering objective of designing, constructing and running a vehicle capable of achieving a speed of 1000 mph on land. The aerodynamic design of this vehicle is to be accomplished using computational simulation only and this paper describes the development and application of the approach adopted. The computational model employs a cell vertex finite volume algorithm for the solution of compressible viscous flow problems on unstructured hybrid meshes. A one equation turbulence model is adopted and the solution of the steady flow equations is obtained by explicit relaxation. For the combination of high Mach number, complex geometry and complex boundary conditions, involving rotating surfaces and a rolling ground, a consistent HLLC numerical flux function is adopted to ensure a stable procedure. To illustrate the impact of the approach upon the final configuration, a number of simulations undertaken to aid the aerodynamic design are described. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

黏性不可压流体的自适应网格技术和基本特性方程分离算法的联合分析

组合基本特性方程分离算法和自适应网格技术,分析二维黏性不可压流体.该方法使用3节点三角单元,对速度分量和压力等变量分析,使用等阶次的插值函数.组合解法的主要优点在于,在自适应网格技术中,对解梯度变化大的区域,通过耦合误差估计生成小的单元,利于提高解的精度,在其它区域生成大单元,可以节省时间.最后,通过对一个黏性流体圆柱体绕流问题的瞬态和稳态特性分析,给出了组合解法性能的评价.     

HIGH ORDER FINITE DIFFERENCE/SPECTRAL METHODS TO A WATER WAVE MODEL WITH NONLOCAL VISCOSITY

In this paper, efficient numerical schemes are proposed for solving the water wave model with nonlocal viscous term that describe the propagation of surface water wave. By using the Caputo fractional derivative definition to approximate the nonlocal fractional operator, finite difference method in time and spectral method in space are constructed for the considered model. The proposed method employs known 5/2 order scheme for fractional derivative and a mixed linearization for the nonlinear term. The analysis shows that the proposed numerical scheme is unconditionally stable and error estimates are provided to predict that the second order backward differentiation plus 5/2 order scheme converges with order 2 in time, and spectral accuracy in space. Several numerical results are provided to verify the efficiency and accuracy of our theoretical claims. Finally, the decay rate of solutions is investigated.     

三维矩形域上泊松方程四面体线元的超逼近与外推

改进三角元的积分恒等式,使之适用于拟一致四面体元,借此证明了泊松方程四面体线元梯度有超逼近现象,函数值Richardson外推可以提高精度.     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

Explicit Hybrid Numerical Method for the Allen-Cahn Type Equations on Curved Surfaces

We present a simple and fast explicit hybrid numerical scheme for the motion by mean curvature on curved surfaces in three-dimensional (3D) space. We numerically solve the Allen-Cahn (AC) and conservative Allen-Cahn (CAC) equations on a triangular surface mesh. We use the operator splitting method and an explicit hybrid numerical method. For the AC equation, we solve the diffusion term using a discrete Laplace-Beltrami operator on the triangular surface mesh and solve the reaction term using the closed-form solution, which is obtained using the separation of variables. Next, for the CAC equation, we additionally solve the time-space dependent Lagrange multiplier using an explicit scheme. Our numerical scheme is computationally fast and efficient because we use an explicit hybrid numerical scheme. We perform various numerical experiments to demonstrate the robustness and efficiency of the proposed scheme.