Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

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

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

流向强迫振荡圆柱绕流的涡脱落模态分析

通过求解采用ALE方法描述的运动坐标系Navier-Stokes方程组，分析均匀来流下雷诺 数为150的静止和流向振荡的圆柱绕流. 主要研究了强迫振荡频率和较大振幅比 (A/D=0.3-1.2)对圆柱升力、阻力变化特性以及涡脱落模态的影响. 研究表 明，流向振荡圆柱绕流存在多种涡脱落模态，如对称S以及反对称A-I, A-III, A-IV等多种形式；比较研究结果，拓展了各模态下对应的锁定区域，并将其分为5个 子区；A-I模态中圆柱受力较以前所知更复杂；通过分析计算结果，发现最大加速度 比Af_{c}^{2}/Df_{s0}^{2}可能是涡脱落模态(尤其是对称S模态)最有效的控制参数.     

充液系统液体-多体耦合动力响应分析

提出了充液系统的液体-多体耦合力学模型，基于ALE有限元法和多体系统动力学理论，发展了液体-多体耦合动力响应分析的一种有效方法. 对于液体子系统，将其运动分解为随同贮箱的大位移运动和相对贮箱的大幅晃动，引入贮箱固连参考系中的任意拉格朗日-欧拉(ALE)运动学描述，建立了贮箱固连非惯性参考系中液体的ALE有限元方程，对液体有限元方程的缩聚大大减少了液体子系统的计算规模. 为了计及液体阻尼的影响，引入了液体修正的Rayleigh阻尼，避免了伪阻尼力的出现. 对于多体子系统，应用多体系统动力学理论建立动力学方程. 在此基础上详细导出了液体-多体耦合动力学方程，并采用预估-多重校正算法(PMA)和时间步长控制算法进行迭代求解，既保证了迭代收敛，又提高了计算效率. 所给算例成功求解了液体运输车辆系统的液体-多体耦合动力响应，深入分析了有关参数对系统动力响应的影响，获得了一些结论.     

大型渡槽动力设计方法研究

根据大型渡槽结构的特点,在缺乏适用的抗震分析模型前提下,提出支撑-渡槽-水伪三维流固耦合动力分析方法。模型基于结构动力分析理论和流固耦合分析理论,整体结构分析采用Newmark方法进行动力时程分析,分析中考虑渡槽结构整体受力性能,包括下部支撑结构、地基等多种因素的影响作用;渡槽-水耦合体子结构分析采用任意拉格朗日-欧拉(Arbitrary Lagrangian Eulerian,ALE)有限元方法进行动力分析,考虑渡槽-水的耦合相互作用,模拟渡槽中水体的大幅晃动作用,相应得出渡槽中动力压力值及其分布。以某排架支撑矩形渡槽为例,采用本文提出的伪三维模型进行地震时程分析,计算结果与三维模型计算结果吻合良好。本文提出的伪三维动力分析方法计算效率高,满足渡槽结构抗震设计要求,便于工程应用,是大型渡槽结构动力分析研究和抗震设计的实用静、动力分析方法。     

Provably stable and time‐accurate extensions of Runge–Kutta schemes for CFD computations on moving grids

Extending fixed‐grid time integration schemes for unsteady CFD applications to moving grids, while formally preserving their numerical stability and time accuracy properties, is a nontrivial task. A general computational framework for constructing stability‐preserving ALE extensions of Eulerian multistep time integration schemes can be found in the literature. A complementary framework for designing accuracy‐preserving ALE extensions of such schemes is also available. However, the application of neither of these two computational frameworks to a multistage method such as a Runge–Kutta (RK) scheme is straightforward. Yet, the RK methods are an important family of explicit and implicit schemes for the approximation of solutions of ordinary differential equations in general and a popular one in CFD applications. This paper presents a methodology for filling this gap. It also applies it to the design of ALE extensions of fixed‐grid explicit and implicit second‐order time‐accurate RK (RK2) methods. To this end, it presents the discrete geometric conservation law associated with ALE RK2 schemes and a method for enforcing it. It also proves, in the context of the nonlinear scalar conservation law, that satisfying this discrete geometric conservation law is a necessary and sufficient condition for a proposed ALE extension of an RK2 scheme to preserve on moving grids the nonlinear stability properties of its fixed‐grid counterpart. All theoretical findings reported in this paper are illustrated with the ALE solution of inviscid and viscous unsteady, nonlinear flow problems associated with vibrations of the AGARD Wing 445.6. Copyright © 2011 John Wiley & Sons, Ltd.     

An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface

In this article, an ALE finite element method to simulate the partial melting of a workpiece of metal is presented. The model includes the heat transport in both the solid and liquid part, fluid flow in the liquid phase by the Navier–Stokes equations, tracking of the melt interface solid/liquid by the Stefan condition, treatment of the capillary boundary accounting for surface tension effects and a radiative boundary condition. We show that an accurate treatment of the moving boundaries is crucial to resolve their respective influences on the flow field and thus on the overall energy transport correctly. This is achieved by a mesh‐moving method, which explicitly tracks the phase boundary and makes it possible to use a sharp interface model without singularities in the boundary conditions at the triple junction. A numerical example describing the welding of a thin‐steel wire end by a laser, where all aforementioned effects have to be taken into account, proves the effectiveness of the approach.Copyright © 2012 John Wiley & Sons, Ltd.     

A thermodynamic and dynamic subgrid closure model for two‐material cells

A general and robust subgrid closure model for two‐material cells is proposed. The conservative quantities of the entire cell are apportioned between two materials, and then, pressure and velocity are fully or partially equilibrated by modeling subgrid wave interactions. An unconditionally stable and entropy‐satisfying solution of the processes has been successfully found. The solution is valid for arbitrary level of relaxation. The model is numerically designed with care for general materials and is computationally efficient without recourse to subgrid iterations or subcycling in time. The model is implemented and tested in the Lagrange‐remap framework. Two interesting results are observed in 1D tests. First, on the basis of the closure model without any pressure and velocity relaxation, a material interface can be resolved without creating numerical oscillations and/or large nonphysical jumps in the problem of the modified Sod shock tube. Second, the overheating problem seen near the wall surface can be solved by the present entropy‐satisfying closure model. The generality, robustness, and efficiency of the model make it useful in principle in algorithms, such as ALE methods, volume of fluid methods, and even some mixture models, for compressible two‐phase flow computations. Copyright © 2013 John Wiley & Sons, Ltd.     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.