Galerkin local maximum entropy method

We propose a new particle‐based method for simulating incompressible Navier‐Stokes flows. It is based on a reinterpretation of the optimal transportation meshfree method within the context of Galerkin discretization. This enables us to introduce the incompressibility constraint into the formulation. Furthermore, we present convergence test and illustrate the usability of the method along several test problems.     

Applications of maximum entropy method in capillary viscometry

Maximum entropy method (MEM) is presented as a technique for processing data obtained from capillary viscometers. The performance of MEM is assessed by comparing the viscosity versus shear rate curves generated by MEM against that obtained by the standard method based on the Weissenberg-Rabinowitsch-Mooney equation. In all the cases in vestigated, MEM proved to be a reliable technique in coping with the experimental noise in the capillary data.     

一类改进最大熵方法的可调参数分析

在结构可靠性分析中,引入含可调参数的转换函数能对传统的最大熵方法进行改进,获得更高的失效概率预测精度。但是,此可调参数的最佳取值很难确定。针对这一问题,引入概率守恒方程,从功能函数转换前后所得概率密度函数出发,建立其最大熵值的变化关系,给出转换前后最大熵值之差的理论形式。通过对三种典型单调非线性转换函数开展算例研究,发现功能函数转换前后的最大熵值之差与转换函数的最佳可调参数值有关。改变可调参数值驱使最大熵值之差变化的同时,改进最大熵方法能遍历到更好的失效概率估计值。     

薄壁圆筒壳初始几何缺陷不确定性量化的极大熵方法

针对薄壁圆筒壳结构轴压屈曲载荷的缺陷敏感性以及真实几何缺陷的不确定性,提出一种基于实测缺陷数据和极大熵原理的初始缺陷建模与屈曲载荷预测方法。首先,将初始几何缺陷视为二维随机场,并利用实测缺陷数据和Karhunen-Loève展开法将初始缺陷的随机场建模转化为随机向量的建模;其次,利用极大熵方法确定随机向量的概率分布;最后,基于所构建的初始缺陷随机模型,利用MCMC抽样方法和确定性屈曲分析方法,进行随机屈曲分析并给出基于可靠度的屈曲载荷折减因子。数值算例表明,与直接假设随机场相关结构的方法相比,本文方法的结果是对薄壁圆筒壳屈曲载荷的一个更无偏估计。     

Pressure‐stabilized maximum‐entropy methods for incompressible Stokes

We present a parameter‐free stable maximum‐entropy method for incompressible Stokes flow. Derived from a least‐biased optimization inspired by information theory, the meshfree maximum‐entropy method appears as an interesting alternative to classical approximation schemes like the finite element method. Especially compared with other meshfree methods, e.g. the moving least‐squares method, it allows for a straightforward imposition of boundary conditions. However, no Eulerian approach has yet been presented for real incompressible flow, encountering the convective and pressure instabilities. In this paper, we exclusively address the pressure instabilities caused by the mixed velocity‐pressure formulation of incompressible Stokes flow. In a preparatory discussion, existing stable and stabilized methods are investigated and evaluated. This is used to develop different approaches towards a stable maximum‐entropy formulation. We show results for two analytical tests, including a presentation of the convergence behavior. As a typical benchmark problem, results are also shown for the leaky lid‐driven cavity. The already presented information‐flux method for convection‐dominated problems in mind, we see this as the last step towards a maximum‐entropy method capable of simulating full incompressible flow problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Direct numerical simulation of modulation of isotropic turbulence by poly‐dispersed particles

A direct numerical simulation technique based on two‐way coupling is presented to study a particle‐laden, decaying isotropic turbulent flow. Physical characteristics of turbulence modulation because of the mono‐dispersed (i.e., particles with single Stokes number) and poly‐dispersed particles (i.e., particles with more than one Stokes number) were investigated. A scale dependent effective viscosity that summarizes the aspects of the interaction between the velocity field and particles is defined in the study. Particles of Stokes number (St) 3.2,6.4 and 12.8 were used in performing the simulations. Poly‐dispersed particles were acquired by mixing particles of two different Stokes numbers at a time. As a whole, decay of turbulence because of the poly‐dispersed particles is observed to be larger than that of the decay of turbulence because of the mono‐dispersed particles. Simulations of poly‐dispersed particle indicate nonlinear characteristics in the modification of the temporal evolution of turbulence energy and dissipation. The scale dependent effective viscosity, which correlates with the energy spectrum plot, indicates that the decay of turbulence is mostly observed at the intermediate scales of turbulence. The effective viscosity for the simulations of the poly‐dispersed particles was calculated to be higher than that of the simulations of the mono‐dispersed particles. Copyright © 2011 John Wiley & Sons, Ltd.     

A filter‐based,mass‐conserving lattice Boltzmann method for immiscible multiphase flows

Some issues of He–Chen–Zhang lattice Boltzmann equation (LBE) method (referred as HCZ model) (J. Comput. Physics 1999; 152 :642–663) for immiscible multiphase flows with large density ratio are assessed in this paper. An extended HCZ model with a filter technique and mass correction procedure is proposed based on HCZ's LBE multiphase model. The original HCZ model is capable of maintaining a thin interface but is prone to generating unphysical oscillations in surface tension and index function at moderate values of density ratio. With a filtering technique, the monotonic variation of the index function across the interface is maintained with larger density ratio. Kim's surface tension formulation for diffuse–interface method (J. Comput. Physics 2005; 204 :784–804) is then used to remove unphysical oscillation in the surface tension. Furthermore, as the density ratio increases, the effect of velocity divergence term neglected in the original HCZ model causes significant unphysical mass sources near the interface. By keeping the velocity divergence term, the unphysical mass sources near the interface can be removed with large density ratio. The long‐time accumulation of the modeling and/or numerical errors in the HCZ model also results in the error of mass conservation of each dispersed phase. A mass correction procedure is devised to improve the performance of the method in this regard. For flows over a stationary and a rising bubble, and capillary waves with density ratio up to 100, the present approach yields solutions with interface thickness of about five to six lattices and no long‐time diffusion, significantly advancing the performance of the LBE method for multiphase flow simulations. Copyright © 2010 John Wiley & Sons, Ltd.     

A curvature‐free multiphase flow solver via surface stress‐based formulation

In this note, we show the link between the classical continuous surface stress and continuous surface force approaches together with special finite element method techniques toward a fully implicit level set method. Based on a modified surface stress formulation, neither normals nor curvature has to be explicitly calculated. The method is space‐dimension independent. Prototypical numerical tests of benchmarking character for a rising 2D bubble are provided for validating the accuracy of this new approach. We show additionally that the explicit redistancing can be avoided using a nonlinear PDE so that a fully implicit and even monolithic formulation of the corresponding multiphase problem gets feasible.     

A CIP/multi‐moment finite volume method for shallow water equations with source terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume‐integrated average (VIA) for each mesh cell, the surface‐integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi‐Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux‐based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non‐oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright © 2007 John Wiley & Sons, Ltd.     

Chebyshev super spectral viscosity solution of a two‐dimensional fluidized‐bed model

The numerical solution of a model describing a two‐dimensional fluidized bed by a Chebyshev super spectral viscosity (SSV) method is considered. The model is in the form of a hyperbolic system of conservation laws with a source term, coupled with an elliptic equation for determining a stream function. The coupled elliptic equation is solved by a finite‐difference method. The mixed SSV/finite‐difference method produces physically shaped bubbles, on a very coarse grid. Fine scale details, which were not present in previous finite‐difference solutions, are present in the solution. Copyright © 2003 John Wiley & Sons, Ltd.     

Quantitative comparison of Taylor flow simulations based on sharp‐interface and diffuse‐interface models

A pressure‐driven flow of elongated bullet‐shaped bubbles in a narrow channel is known as Taylor flow or bubble‐train flow. This process is of relevance in various applications of chemical engineering. In this paper, we describe a typical simplified experimental setting, with surface tension, density and viscosity as prescribed input parameters. We compare a sharp‐interface model based on a moving grid aligned with the bubble boundary (ALE coordinates) and a diffuse‐interface model where the bubble shape is implicitly given by a phase‐field function. Four independent implementations based on the two modeling approaches are introduced and described briefly. Besides the simulation of the bubble shapes, we compare some resulting quantities such as pressure difference and film widths within the implementations and to existing analytical and experimental results. The simulations were conducted in 2D and 3D (rotationally symmetric). Good accordance of the results indicate the applicability and the usability of all approaches. Differences between the models and their implementations are visible but in no contradiction to theoretical results. Copyright © 2013 John Wiley & Sons, Ltd.     

Improving the convergence behaviour of a fixed‐point‐iteration solver for multiphase flow in porous media

A new method to admit large Courant numbers in the numerical simulation of multiphase flow is presented. The governing equations are discretized in time using an adaptive θ‐method. However, the use of implicit discretizations does not guarantee convergence of the nonlinear solver for large Courant numbers. In this work, a double‐fixed point iteration method with backtracking is presented, which improves both convergence and convergence rate. Moreover, acceleration techniques are presented to yield a more robust nonlinear solver with increased effective convergence rate. The new method reduces the computational effort by strengthening the coupling between saturation and velocity, obtaining an efficient backtracking parameter, using a modified version of Anderson's acceleration and adding vanishing artificial diffusion. © 2016 The Authors. International Journal for Numerical Methods in Fluids Published by John Wiley & Sons Ltd.     

An extended mixture model for the simultaneous treatment of small‐scale and large‐scale interfaces

The aim of this work is to present a new model based on the volume of fluid method and the algebraic slip mixture model in order to solve multiphase gas–fluid flows with different interface scales and the transition among them. The interface scale is characterized by a measure of the grid, which acts as a geometrical filter and is related with the accuracy in the solution; in this sense, the presented coupled model allows to reduce the grid requirements for a given accuracy. With this objective in mind, a generalization of the algebraic slip mixture model is proposed to solve problems involving small‐scale and large‐scale interfaces in an unified framework taking special care in preserving the conservativeness of the fluxes. This model is implemented using the OpenFOAM^® libraries to generate a tool capable of solving large problems on high‐performance computing facilities. Several examples are solved as a validation for the presented model, including new quantitative measurements to assess the advantages of the method. Copyright © 2014 John Wiley & Sons, Ltd.     

A robust multi‐grid pressure‐based algorithm for multi‐fluid flow at all speeds

This paper reports on the implementation and testing, within a full non‐linear multi‐grid environment, of a new pressure‐based algorithm for the prediction of multi‐fluid flow at all speeds. The algorithm is part of the mass conservation‐based algorithms (MCBA) group in which the pressure correction equation is derived from overall mass conservation. The performance of the new method is assessed by solving a series of two‐dimensional two‐fluid flow test problems varying from turbulent low Mach number to supersonic flows, and from very low to high fluid density ratios. Solutions are generated for several grid sizes using the single grid (SG), the prolongation grid (PG), and the full non‐linear multi‐grid (FMG) methods. The main outcomes of this study are: (i) a clear demonstration of the ability of the FMG method to tackle the added non‐linearity of multi‐fluid flows, which is manifested through the performance jump observed when using the non‐linear multi‐grid approach as compared to the SG and PG methods; (ii) the extension of the FMG method to predict turbulent multi‐fluid flows at all speeds. The convergence history plots and CPU‐times presented indicate that the FMG method is far more efficient than the PG method and accelerates the convergence rate over the SG method, for the problems solved and the grids used, by a factor reaching a value as high as 15. Copyright © 2003 John Wiley & Sons, Ltd.     

A numerical method of K-S entropy calculation for a strange attractor

Based directly on the original definition of K-S entropy, a new algorithm for calculating K-S entropy from chaotic time series is developed by using some techniques of coding and code operation. The project supported by National Natural Science Foundation of China     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

Convergence acceleration by self‐adjusted time stepsize using Bi‐CGSTAB method for turbulent separated flow computation

Poor convergence behavior is usually encountered when numerical computations on turbulent separated flow are performed. A design of self‐adjusted stepsize concept both in time span and spatial coordinate systems to achieve faster convergence is demonstrated in this study. The determination of the time stepsize based on the concept of minimization of residuals using the Bi‐CGSTAB algorithm is proposed. The numerical results show that the time stepsize adjusted by the proposed method indeed improves the convergence rate for turbulent separated flow computations using advanced turbulence models in low‐Reynolds number forms. Copyright © 2002 John Wiley & Sons, Ltd.     

Pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist

We propose a pressure‐based unified solver for gas‐liquid two‐phase flows where compressible and incompressible flows coexist. Unlike the original thermo–Cubic Interpolated Propagation Combined Unified Procedure (CIP‐CUP) method proposed by Himeno et al (Transactions of the Japan Society of Mechanical Engineers, Series B, 2003), we split the advection term of the governing equations into a conservation part and into the rest. The splitting of advection term has two advantages. One is the high degree of freedom in choosing discretization schemes such as central‐difference schemes, upwind schemes, and Total Variation Diminishing (TVD) schemes. The other is the ease of implementation on unstructured grids. The advantages enable the analyses of various flows such as turbulent and supersonic ones in actual complicated boundaries. Therefore, the solver is useful for practical analyses. The solver was validated on the following test cases: subsonic single‐phase flows, incompressible single‐phase turbulent flows, and incompressible gas‐liquid two‐phase flows. With unstructured grids, we obtained the equivalent results as the ones with structured grids. After the validations, subsonic jet impinging on a water pool was calculated and compared with experimental results. It was confirmed that the calculated results were consistent with the experimental ones.     

An evaluation of a 3D free‐energy‐based lattice Boltzmann model for multiphase flows with large density ratio

In this paper, the 3D Navier–Stokes (N–S) equation and Cahn–Hilliard (C–H) equations were solved using a free‐energy‐based lattice Boltzmann (LB) model. In this model, a LB equation with a D3Q19 velocity model is used to recover continuity and N–S equations while another LB equation with D3Q7 velocity model for solving C–H equation (Int. J. Numer. Meth. Fluids, 2008; 56 :1653–1671) is applied to solve the 3D C–H equation. To avoid the excessive use of computational resources, a moving reference frame is adopted to allow long‐time simulation of a bubble rising. How to handle the inlet/outlet and moving‐wall boundary conditions are suggested. These boundary conditions are simple and easy for implementation. This model's performance on two‐phase flows was investigated and the mass conservation of this model was evaluated. The model is validated by its application to simulate the 3D air bubble rising in viscous liquid (density ratio is 1000). Good agreement was obtained between the present numerical results and experimental results when Re is small. However, for high‐Re cases, the mass conservation seems not so good as the low‐Re case. Copyright © 2009 John Wiley & Sons, Ltd.     

An analytical interface reconstruction algorithm in the PLIC‐VOF method for 2D polygonal unstructured meshes

Piecewise linear interface calculation (PLIC) schemes have been extensively employed in the volume‐of‐fluid (VOF) method for interface capturing in numerical simulations of multiphase flows. Polygonal unstructured meshes are often adopted because of their geometric flexibility and superiority in gradient calculation. An analytical interface reconstruction algorithm in the PLIC‐VOF method for arbitrary convex polygonal cells has been proposed in this study. The line interface at a given orientation within a polygonal cell is located by an analytical technique. It has been tested successfully for four different geometric shapes that are common in polygonal meshes. The computational efficiency of the present algorithm has been compared with several published schemes in the literature. The proposed algorithm has been shown to yield higher accuracy with reduction in computational complexity. A numerical simulation of a dam‐breaking problem has been performed using the proposed analytical PLIC technique on polygonal meshes. The results are in good agreement with experimental data available in the literature, which serves as a demonstration of its performance in a real multiphase flow.