Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid particle volume-of-fluid method for curvature estimation in multiphase flows
185
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a particle method for estimating the curvature of interfaces in volume-of-fluid simulations of multiphase flows. The method is well suited for under-resolved interfaces, and it is shown to be more accurate than the parabolic fitting that is employed in such cases. The curvature is computed from the equilibrium positions of particles constrained to circular arcs and attracted to the interface. The proposed particle method is combined with the method of height functions at higher resolutions, and it is shown to outperform the current combinations of height functions and parabolic fitting. The algorithm is conceptually simple and straightforward to implement on new and existing software frameworks for multiphase flow simulations thus enhancing their capabilities in challenging flow problems. We evaluate the proposed hybrid method on a number of two- and three-dimensional benchmark flow problems and illustrate its capabilities on simulations of flows involving bubble coalescence and turbulent multiphase flows.
rate research
Read More
In this paper, we develop a simplified hybrid weighted essentially non-oscillatory (WENO) method combined with the modified ghost fluid method (MGFM) [28] to simulate the compressible two-medium flow problems. The MGFM can turn the two-medium flow problems into two single-medium cases by defining the ghost fluids status in terms of the predicted the interface status, which makes the material interface invisible. For the single medium flow case, we adapt between the linear upwind scheme and the WENO scheme automatically by identifying the regions of the extreme points for the reconstruction polynomial as same as the hybrid WENO scheme [50]. Instead of calculating their exact locations, we only need to know the regions of the extreme points based on the zero point existence theorem, which is simpler for implementation and saves computation time. Meanwhile, it still keeps the robustness and has high efficiency. Extensive numerical results for both one and two dimensional two-medium flow problems are performed to demonstrate the good performances of the proposed method.
A computational fluid dynamics (CFD) simulation framework for predicting complex flows is developed on the Tensor Processing Unit (TPU) platform. The TPU architecture is featured with accelerated performance of dense matrix multiplication, large high bandwidth memory, and a fast inter-chip interconnect, which makes it attractive for high-performance scientific computing. The CFD framework solves the variable-density Navier-Stokes equation using a Low-Mach approximation, and the governing equations are discretized by a finite difference method on a collocated structured mesh. It uses the graph-based TensorFlow as the programming paradigm. The accuracy and performance of this framework is studied both numerically and analytically, specifically focusing on effects of TPU-native single precision floating point arithmetic on solution accuracy. The algorithm and implementation are validated with canonical 2D and 3D Taylor Green vortex simulations. To demonstrate the capability for simulating turbulent flows, simulations are conducted for two configurations, namely the decaying homogeneous isotropic turbulence and a turbulent planar jet. Both simulations show good statistical agreement with reference solutions. The performance analysis shows a linear weak scaling and a super-linear strong scaling up to a full TPU v3 pod with 2048 cores.
We present a novel computational modeling framework to numerically investigate fluid-structure interaction in viscous fluids using the phase field embedding method. Each rigid body or elastic structure immersed in the incompressible viscous fluid matrix, grossly referred to as the particle in this paper, is identified by a volume preserving phase field. The motion of the particle is driven by the fluid velocity in the matrix for passive particles or combined with its self-propelling velocity for active particles. The excluded volume effect between a pair of particles or between a particle and the boundary is modeled by a repulsive potential force. The drag exerted to the fluid by a particle is assumed proportional to its velocity. When the particle is rigid, its state is described by a zero velocity gradient tensor within the nonzero phase field that defines its profile and a constraining stress exists therein. While the particle is elastic, a linear constitutive equation for the elastic stress is provided within the particle domain. A hybrid, thermodynamically consistent hydrodynamic model valid in the entire computational domain is then derived for the fluid-particle ensemble using the generalized Onsager principle accounting for both rigid and elastic particles. Structure-preserving numerical algorithms are subsequently developed for the thermodynamically consistent model. Numerical tests in 2D and 3D space are carried out to verify the rate of convergence and numerical examples are given to demonstrate the usefulness of the computational framework for simulating fluid-structure interactions for passive as well as self-propelling active particles in a viscous fluid matrix.
In this paper, a high order quasi-conservative discontinuous Galerkin (DG) method using the non-oscillatory kinetic flux is proposed for the 5-equation model of compressible multi-component flows with Mie-Gruneisen equation of state. The method mainly consists of three steps: firstly, the DG method with the non-oscillatory kinetic flux is used to solve the conservative equations of the model; secondly, inspired by Abgralls idea, we derive a DG scheme for the volume fraction equation which can avoid the unphysical oscillations near the material interfaces; finally, a multi-resolution WENO limiter and a maximum-principle-satisfying limiter are employed to ensure oscillation-free near the discontinuities, and preserve the physical bounds for the volume fraction, respectively. Numerical tests show that the method can achieve high order for smooth solutions and keep non-oscillatory at discontinuities. Moreover, the velocity and pressure are oscillation-free at the interface and the volume fraction can stay in the interval [0,1].
In this paper, authors focus effort on improving the conventional discrete velocity method (DVM) into a multiscale scheme in finite volume framework for gas flow in all flow regimes. Unlike the typical multiscale kinetic methods unified gas-kinetic scheme (UGKS) and discrete unified gas-kinetic scheme (DUGKS), which concentrate on the evolution of the distribution function at the cell interface, in the present scheme the flux for macroscopic variables is split into the equilibrium part and the nonequilibrium part, and the nonequilibrium flux is calculated by integrating the discrete distribution function at the cell center, which overcomes the excess numerical dissipation of the conventional DVM in the continuum flow regime. Afterwards, the macroscopic variables are finally updated by simply integrating the discrete distribution function at the cell center, or by a blend of the increments based on the macroscopic and the microscopic systems, and the multiscale property is achieved. Several test cases, involving unsteady, steady, high speed, low speed gas flows in all flow regimes, have been performed, demonstrating the good performance of the multiscale DVM from free molecule to continuum Navier-Stokes solutions and the multiscale property of the scheme is proved.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
