A Cartesian grid method for two‐phase gel dynamics on an irregular domain

We develop a numerical method for simulating models of two‐phase gel dynamics in an irregular domain using a regular Cartesian grid. The models consist of transport equations for the volume fractions of the two phases, polymer network and solvent; coupled momentum equations for the two phases; and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box relaxation scheme is used as a preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Ghost points are used to enforce no‐slip boundary conditions for the velocity field of each phase, and no‐flux boundary conditions for the volume fractions. The behavior of the new method, including its rate of convergence, is explored through numerical experiments for a problem in which strong phase separation develops from an initially (almost) homogeneous phase distribution. We also use the method to explore situations, motivated by biology, which show that imposed boundary velocities can cause substantial redistribution of network and solvent. Copyright © 2010 John Wiley & Sons, Ltd.     

基于格子Boltzmann方法的粘弹性流体圆柱绕流分析

将光滑界面法引入到格子Boltzmann方法中分析粘弹性流体绕流问题,分别采用单松弛模型和对流扩散模型求解运动方程和Oldroyd-B本构方程,针对圆形和椭圆内部边界条件,给出连续界面插值函数,在此基础上,运用光滑界面法将内部边界转换为作用力项施加到演化方程中。首先分析圆柱绕流问题,给出不同材料参数情况下的流场分布和阻力系数计算结果,比较发现与宏观数值模拟结果相吻合。将模型拓展到绕椭圆流动中,分析椭圆形状和材料参数对粘弹性流体绕柱流的影响,发现随着椭圆长轴与短轴比值的增加和维森伯格数的增加,阻力系数逐渐下降,并且长短轴比对迭代收敛有较大影响。     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

Numerical prediction of extrudate swell of branched polymer melts

This paper is concerned with the numerical prediction of the extrudate swell behaviour of branched polymer melts in a planar configuration. The multi-mode extended pom-pom (XPP) model is used to describe the polymer dynamics. A second-order operator-integration-factor splitting scheme is used for the temporal discretisation of the problem, whilst a spectral element scheme is used in space. The free surface is evolved in a Lagrangian manner using the third-order conditionally stable Adams–Bashforth method. A thorough mesh convergence study is performed with respect to the temporal and spatial discretisation parameters. The influence of the nature of the discrete relaxation spectrum on the swelling ratio and as an indicator of polydispersity is investigated. The predictions of numerical simulations are also compared with a selection of experimental results from the literature. The parameters in the XPP model are determined from rheological data. Good agreement is obtained for branched low-density polyethylenes. The ability to model a melt with a high molecular weight tail using a discrete relaxation spectrum for which the largest relaxation time is isolated from the others is also investigated.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

A novel implicit numerical treatment for non‐linear turbulence models using high and low Reynolds number formulations

Non‐linear turbulence models can be seen as an improvement of the classical eddy‐viscosity concept due to their better capacity to simulate characteristics of important flows. However, application of non‐linear models demand robustness of the numerical method applied, requiring a stable discretization scheme for convergence of all variables involved. Usually, non‐linear terms are handled in an explicit manner leading to possible numerical instabilities. Thus, the present work shows the steps taken to adapt a general non‐linear constitutive equation using a new semi‐implicit numerical treatment for the non‐linear diffusion terms. The objective is to increase the degree of implicitness of the solution algorithm to enhance convergence characteristics. Flow over a backward‐facing step was computed using the control volume method applied to a boundary‐fitted coordinate system. The SIMPLE algorithm was used to relax the algebraic equations. Classical wall function and a low Reynolds number model were employed to describe the flow near the wall. The results showed that for certain combination of relaxation parameters, the semi‐implicit treatment proposed here was the sole successful treatment in order to achieve solution convergence. Also, application of the implicit method described here shows that the stability of the solution either increases (high Reynolds with non‐orthogonal mesh) or preserves the same (low Reynolds number applications). Additional advantages of the procedure proposed here lie in the possibility of testing different non‐linear expressions if one considers the enhanced robustness and stability obtained for the entire numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

考虑侧面摩擦的侧限压缩模型的数值解法

海上箱筒型基础结构下沉到位后,舱内土体的受力变形可以看作是侧限压缩模型问题,上部荷载引起的侧壁摩擦力与土体应力之间的关系较为复杂。本文通过简化模型,列出了侧限压缩模型的平衡微分方程,分析了方程解析计算时存在的难点,结合摩擦应力边界条件的特点,构造函数迭代法对该方程进行数值求解,展示了计算的收敛过程,得到接触侧面上水平应力和土体表面位移的非线性分布规律,并将其结果与ABAQUS数值模拟结果进行了对比。结果表明,函数迭代法收敛速度较快,计算过程稳定,对其他存在变量耦合边界条件的微分方程数值求解有一定的参考价值。     

Analytical and numerical studies of the boundary slip in the immersed boundary‐thermal lattice Boltzmann method

We analytically and numerically investigate the boundary slip, including the velocity slip and the temperature jump, in immersed boundary‐thermal lattice Boltzmann methods (IB‐TLBMs) with the two‐relaxation‐time collision operator. We derive the theoretical equation for the relaxation parameters considering the effect of the advection velocity on the temperature jump of the IB‐TLBMs. The analytical and numerical solutions demonstrate that the proposed iterative correction methods without the computational cost of the sparse matrix solver reduce the boundary slip and boundary‐value deviation as effectively as the implicit correction method for any relaxation time. Because the commonly used multi‐direct forcing method does not consider the contributions of the body force to the momentum flux, it cannot completely eliminate the boundary slip because of the numerical instability for a long relaxation time. Both types of proposed iterative correction methods are more numerically stable than the implicit correction method. In simulations of flow past a circular cylinder and of natural convection, the present iterative correction methods yield adequate results without the errors of the velocity slip, the temperature jump, and the boundary‐value deviation for any relaxation time parameters and for any number of Lagrangian points per length. The combination of the present methods and the two‐relaxation‐time collision operator is suitable for simulating fluid flow with thermal convection in the multiblock method in which the relaxation time increases in inverse proportion to the grid size.     

基于Kriging模型和改进MCMC算法的随机有限元模型修正

针对待修正参数维数较高时,标准马尔可夫链蒙特卡罗MCMC (Markov Chain Monte Carlo)算法不易收敛、拒绝率高的问题,提出了基于Kriging模型和在MCMC中融合花朵授粉算法的修正方法.首先,以待修正参数作为输入,以应变模态作为输出,建立Kriging模型,通过蝙蝠算法确定Kriging模型的相关系数;然后,采用最大熵的贝叶斯方法估计参数的后验概率密度函数,将花朵授粉算法融入MH (M etropolis-Hasting)抽样算法,提高局部寻优和全局寻优能力;最后,通过三自由度弹簧-质量系统和三维桁架结构的数值算例验证所提模型修正方法,修正后参数相对误差均低于0.86％.结果 表明,所提方法修正后较高维参数的马尔可夫链能够快速收敛且样本接受率也有所提高,该方法也对随机噪声具有一定的鲁棒性.     

二维弹性新型快速多极虚边界元的最小二乘法

在虚边界元最小二乘法的方程求解中采用新型的快速多极展开和广义极小残值法,提出了一种二维弹性新型快速多极虚边界元最小二乘法的求解思想。基于二维弹性问题原有的快速多极虚边界元最小二乘法的展开格式,通过引入对角化的概念,以更新展开传递格式;相对于原有快速多极算法,该方法可进一步提高计算效率且仍能保证具有较高的计算精度。数值算例说明了该方法的可行性、计算效率、计算精度均较高。     

Simulation of incompressible multiphase flows with complex geometry using etching multiblock method

The incompressible two-phase flows are simulated using combination of an etching multiblock method and a diffuse interface (DI) model, particularly in the complex domain that can be decomposed into multiple rectangular subdomains. The etching multiblock method allows natural communications between the connected subdomains and the efficient parallel computation. The DI model can consider two-phase flows with a large density ratio, and simulate the flows with the moving contact line (MCL) when a geometric formulation of the MCL model is included. Therefore, combination of the etching method and the DI model has potential to deal with a variety of two-phase flows in industrial applications. The performance is examined through a series of numerical experiments. The convergence of the etching method is firstly tested by simulating single-phase flows past a square cylinder, and the method for the multiphase flow simulation is validated by investing drops dripping from a pore. The numerical results are compared with either those from other researchers or experimental data. Good agreement is achieved. The method is also used to investigate the impact of a droplet on a grooved substrate and droplet generation in flow focusing devices.     

A hybrid phase field multiple relaxation time lattice Boltzmann method for the incompressible multiphase flow with large density contrast

A hybrid phase field multiple relaxation time lattice Boltzmann method (LBM) is presented in this paper for simulation of multiphase flows with large density contrast. In the present method, the flow field is solved by a lattice Boltzmann equation. Concurrently, the interface of two fluids is captured by solving the macroscopic Cahn‐Hilliard equation using the upwind scheme. To be specific, for simulation of the flow field, an lattice Boltzmann equation (LBE) model developed in Shao et al. (Physical Review E, 89 (2014), 033309) for consideration of density contrast in the momentum equation is used. Moreover, in the present work, the multiple relaxation time collision operator is applied to this LBE to enable simulation of problems with large viscosity contrast or high Reynolds number. For the interface capturing, instead of solving another set of LBE as in many phase field LBMs, the macroscopic Cahn‐Hilliard equation is directly solved by using a weighted essentially non‐oscillatory scheme. In this way, the present hybrid phase field LBM shares full advantages of the phase field LBM while enhancing numerical stability. The ability of the present method to simulate multiphase flow problems with large density contrast is demonstrated by several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Large eddy simulation of high-Reynolds-number atmospheric boundary layer flow with improved near-wall correction

It is highly attractive to develop an efficient and flexible large eddy simulation(LES) technique for high-Reynolds-number atmospheric boundary layer(ABL) simulation using the low-order numerical scheme on a relatively coarse grid, that could reproduce the logarithmic profile of the mean velocity and some key features of large-scale coherent structures in the outer layer. In this study, an improved near-wall correction scheme for the vertical gradient of the resolved streamwise velocity in the s...     

Finite element framework for describing dynamic wetting phenomena

The finite element simulation of dynamic wetting phenomena, requiring the computation of flow in a domain confined by intersecting a liquid–fluid free surface and a liquid–solid interface, with the three‐phase contact line moving across the solid, is considered. For this class of flows, different finite element method (FEM) implementations have been used in the literature, and in some cases, these produced apparently contradictory results. In the present paper, a robust framework for the FEM simulation of dynamic wetting flows is developed, which, by consistently adhering to the FEM methodology, leaves no room for ad hoc ‘optional’ variations in the numerical handling of these flows. The developed approach makes it possible to conduct a convergence study, assess the spatial resolution required to achieve a preset accuracy and provide the corresponding benchmark calculations. This analysis allows one to identify numerical artefacts, which had previously been interpreted as physical effects, and demonstrates that suppressing numerical errors using a ‘strong’ implementation of a boundary condition creates bigger and less detectable errors elsewhere in the computational domain. We provide practical recommendations on the spatial resolution required by a numerical scheme for a given set of non‐dimensional similarity parameters and give a user‐friendly step‐by‐step guide specifying the entire implementation, which allows the reader to easily reproduce all presented results including the benchmark calculations. It is also shown how the developed framework accommodates generalizations of the mathematical model accounting for additional physical effects, such as gradients in surface tensions. Copyright © 2011 John Wiley & Sons, Ltd.     

Unsupervised neural network model optimized with evolutionary computations for solving variants of nonlinear MHD Jeffery-Hamel problem

A heuristic technique is developed for a nonlinear magnetohydrodynamics(MHD) Jeffery-Hamel problem with the help of the feed-forward artificial neural network(ANN) optimized with the genetic algorithm(GA) and the sequential quadratic programming(SQP) method. The two-dimensional(2D) MHD Jeffery-Hamel problem is transformed into a higher order boundary value problem(BVP) of ordinary differential equations(ODEs). The mathematical model of the transformed BVP is formulated with the ANN in an unsupervised manner. The training of the weights of the ANN is carried out with the evolutionary calculation based on the GA hybridized with the SQP method for the rapid local convergence. The proposed scheme is evaluated on the variants of the Jeffery-Hamel flow by varying the Reynold number, the Hartmann number, and the angles of the walls. A large number of simulations are performed with an extensive analysis to validate the accuracy, convergence, and effectiveness of the scheme. The comparison of the standard numerical solution and the analytic solution establishes the correctness of the proposed designed methodologies.     

快速多极虚边界元法对含圆孔薄板有效弹性模量的模拟分析

针对虚边界元法,引入快速多极展开和广义极小残值法(GMRES)的思想,以形成快速多极虚边界元法的求解思想,并将此方法用于含圆孔薄板有效弹性模量的模拟分析.由于本文方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,从而使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.本文工作的研究目的在于:提高虚边界元法在普通台式机上的运算能力和拓宽虚边界元法对大规模复杂问题的求解(或数值模拟).文中给出了均布圆孔的正方形薄板和之字形分布圆孔薄板二个算例,以验证该方法的可行性,计算精度和计算效率.     

湍流横掠固定颗粒的全尺度直接数值模拟

气固两相流模拟中,当固相尺度接近或大于Kolmogorov尺度时,普通的点源模型将不再适用,固体相的体积效应和表面效应将对流体相产生显著的影响。通过采用直接数值模拟方法,结合内嵌边界方法对湍流中不同湍流强度流体横掠大于Kolmogorov尺度的固相颗粒进行了全尺度模拟,讨论分析了在两种湍流度下方形颗粒对湍流的调制影响以及颗粒的受力情况。     

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.     

使用新型物态方程的超高速碰撞物质点法模拟

为更准确地对超高速碰撞进行数值模拟、获得与实验结果相似度更高的碎片云形态,利用分子动力学方法求解材料的冷能、冷压,并结合Grover定标律方程,建立了一种表达形式简洁、可处理相变影响的新型物态方程,并代入自编柱坐标物质点法计算程序,使用新型物态方程计算所得的碎片云与使用Mie-Grüneisen、Tillotson等传统物态方程的计算结果相比,在尺寸、形态方面均能够与实验结果更好地吻合,证明了新型物态方程的有效性。