Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver
96
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock-bubble interaction, and gas bubble cavitation. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock-bubble, shock-droplet, and shock-water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas-liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock-bubble-vessel-wall and acoustic-bubble-net interactions are used to demonstrate the full capabilities of MFC.
rate research
Read More
A high-performance gas kinetic solver using multi-level parallelization is developed to enable pore-scale simulations of rarefied flows in porous media. The Boltzmann model equation is solved by the discrete velocity method with an iterative scheme. The multi-level MPI/OpenMP parallelization is implemented with the aim to efficiently utilise the computational resources to allow direct simulation of rarefied gas flows in porous media based on digital rock images for the first time. The multi-level parallel approach is analyzed in details confirming its better performance than the commonly-used MPI processing alone for an iterative scheme. With high communication efficiency and appropriate load balancing among CPU processes, parallel efficiency of 94% is achieved for 1536 cores in the 2D simulations, and 81% for 12288 cores in the 3D simulations. While decomposition in the spatial space does not affect the simulation results, one additional benefit of this approach is that the number of subdomains can be kept minimal to avoid deterioration of the convergence rate of the iteration process. This multi-level parallel approach can be readily extended to solve other Boltzmann model equations.
In this paper, we present an open-source multi-resolution and multi-physics library: SPHinXsys (pronunciation: sfinksis) which is an acronym for underline{S}moothed underline{P}article underline{H}ydrodynamics (SPH) for underline{in}dustrial compleunderline{X} underline{sys}tems. As an open-source library, SPHinXsys is developed and released under the terms of Apache License (2.0). Along with the source code, a complete documentation is also distributed to make the compilation and execution easy. SPHinXsys aims at modeling coupled multi-physics industrial dynamic systems including fluids, solids, multi-body dynamics and beyond, in a multi-resolution unified SPH framework. As an SPH solver, SPHinXsys has many advantages namely, (1) the generic design provides a C++ API showing a very good flexibility when building domain-specific applications, (2) numerous industrial or scientific applications can be coupled within the same framework and (3) with the open-source philosophy, the community of users can collaborate and improve the library. SPHinXsys presently (v0.2.0) includes validations and applications in the fields of fluid dynamics, solid dynamics, thermal and mass diffusion, reaction-diffusion, electromechanics and fluid-structure interactions (FSI).
Current multi-component, multiphase pseudo-potential lattice Boltzmann models have thermodynamic inconsistencies that prevent them to correctly predict the thermodynamic phase behavior of partially miscible multi-component mixtures, such as hydrocarbon mixtures. This paper identifies these inconsistencies and attempts to design a thermodynamically consistent multi-component, multiphase pseudo-potential lattice Boltzmann model that allows mass transfer across the phase interfaces and is capable to predict the phase behavior of typically partially miscible hydrocarbon mixtures. The designed model defines the total interaction force for the entire phase and split the force into individual components. Through a properly derived force split factor associated with the volatility of each component, the model can achieve precise thermodynamic consistency in multi-component hydrocarbon mixtures, which is described by the iso-fugacity rule.
The general synthetic iteration scheme (GSIS) is extended to find the steady-state solution of nonlinear gas kinetic equation, removing the long-standing problems of slow convergence and requirement of ultra-fine grids in near-continuum flows. The key ingredients of GSIS are that the gas kinetic equation and macroscopic synthetic equations are tightly coupled, and the constitutive relations in macroscopic synthetic equations explicitly contain Newtons law of shear stress and Fouriers law of heat conduction. The higher-order constitutive relations describing rarefaction effects are calculated from the velocity distribution function, however, their constructions are simpler than our previous work (Su et al. Journal of Computational Physics 407 (2020) 109245) for linearized gas kinetic equations. On the other hand, solutions of macroscopic synthetic equations are used to inform the evolution of gas kinetic equation at the next iteration step. A rigorous linear Fourier stability analysis in periodic system shows that the error decay rate of GSIS can be smaller than 0.5, which means that the deviation to steady-state solution can be reduced by 3 orders of magnitude in 10 iterations. Other important advantages of the GSIS are (i) it does not rely on the specific form of Boltzmann collision operator and (ii) it can be solved by sophisticated techniques in computational fluid dynamics, making it amenable to large scale engineering applications. In this paper, the efficiency and accuracy of GSIS is demonstrated by a number of canonical test cases in rarefied gas dynamics.
We propose a multi-resolution strategy that is compatible with the lattice Greens function (LGF) technique for solving viscous, incompressible flows on unbounded domains. The LGF method exploits the regularity of a finite-volume scheme on a formally unbounded Cartesian mesh to yield robust and computationally efficient solutions. The original method is spatially adaptive, but challenging to integrate with embedded mesh refinement as the underlying LGF is only defined for a fixed resolution. We present an ansatz for adaptive mesh refinement, where the solutions to the pressure Poisson equation are approximated using the LGF technique on a composite mesh constructed from a series of infinite lattices of differing resolution. To solve the incompressible Navier-Stokes equations, this is further combined with an integrating factor for the viscous terms and an appropriate Runge-Kutta scheme for the resulting differential-algebraic equations. The parallelized algorithm is verified through with numerical simulations of vortex rings, and the collision of vortex rings at high Reynolds number is simulated to demonstrate the reduction in computational cells achievable with both spatial and refinement adaptivity.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
