A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

三维复杂流动的间断有限元方法模拟

本文基于三维可压缩Euler方程,采用基于Runge-Kutta时间离散的间断有限元方法(RKDG方法),对三维前台阶、三维Riemann问题和球Riemann等问题进行了模拟。结果表明,本文的RKDG方法能够在很少的网格内清晰地捕捉到三维复杂流场中的激波和接触间断;同时,将球Riemann问题中z=0.4平面压强沿到对称轴距离的分布与文献中的近似精确解相比,吻合较好,这也验证了本文的RKDG方法不仅能够进行三维复杂流场的定性描述,也能够应用于三维复杂流场的定量计算。     

Invariant discrete flows

In this paper, we investigate the evolution of joint invariants under invariant geometric flows using the theory of equivariant moving frames and the induced invariant discrete variational complex. For certain arc length preserving planar curve flows invariant under the special Euclidean group , the special linear group , and the semidirect group , we find that the induced evolution of the discrete curvature satisfies the differential‐difference mKdV, KdV, and Burgers' equations, respectively. These three equations are completely integrable, and we show that a recursion operator can be constructed by precomposing the characteristic operator of the curvature by a certain invariant difference operator. Finally, we derive the constraint for the integrability of the discrete curvature evolution to lift to the evolution of the discrete curve itself.     

Time‐periodic solution to the compressible nematic liquid crystal flows in periodic domain

In this paper, we consider the time‐periodic solution to a simplified version of Ericksen‐Leslie equations modeling the compressible hydrodynamic flow of nematic liquid crystals with a time‐periodic external force in a periodic domain in . By using an approach of parabolic regularization and combining with the topology degree theory, we establish the existence of the time‐periodic solution to the model under some smallness and symmetry assumptions on the external force. Then, we give the uniqueness of the periodic solution of this model.     

A large deviations principle for stochastic flows of viscous fluids

We study the well-posedness of a stochastic differential equation on the two dimensional torus ${\mathbb{T}}^2$, driven by an infinite dimensional Wiener process with drift in the Sobolev space $L^2(0,T;H^1\left({\mathbb{T}}^2\right))$. The solution corresponds to a stochastic Lagrangian flow in the sense of DiPerna Lions. By taking into account that the motion of a viscous incompressible fluid on the torus can be described through a suitable stochastic differential equation of the previous type, we study the inviscid limit. By establishing a large deviations principle, we show that, as the viscosity goes to zero, the Lagrangian stochastic Navier–Stokes flow approaches the Euler deterministic Lagrangian flow with an exponential rate function.     

The asymptotic limits of Riemann solutions for the isentropic drift-flux model of compressible two-phase flows

The formation of vacuum state and delta shock wave are observed and studied in the limits of Riemann solutions for the one-dimensional isentropic drift-flux model of compressible two-phase flows by letting the pressure in the mixture momentum equation tend to zero. It is shown that the Riemann solution containing two rarefaction waves and one contact discontinuity turns out to be the solution containing two contact discontinuities with the vacuum state between them in the limiting situation. By comparison, it is also proved rigorously in the sense of distributions that the Riemann solution containing two shock waves and one contact discontinuity converges to a delta shock wave solution under this vanishing pressure limit.