Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImproving the staggered grid Lagrangian hydrodynamics for modeling multi-material flows
77
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this work, we make two improvements on the staggered grid hydrodynamics (SGH) Lagrangian scheme for modeling 2-dimensional compressible multi-material flows on triangular mesh. The first improvement is the construction of a dynamic local remeshing scheme for preventing mesh distortion. The remeshing scheme is similar to many published algorithms except that it introduces some special operations for treating grids around multi-material interfaces. This makes the simulation of extremely deforming and topology-variable multi-material processes possible, such as the complete process of a heavy fluid dipping into a light fluid. The second improvement is the construction of an Euler-like flow on each edge of the mesh to count for the edge-bending effect, so as to mitigate the checkerboard oscillation that commonly exists in Lagrangian simulations, especially the triangular mesh based simulations. Several typical hydrodynamic problems are simulated by the improved staggered grid Lagrangian hydrodynamic method to test its performance.
rate research
Read More
We present the second-order multidimensional Staggered Grid Hydrodynamics Residual Distribution (SGH RD) scheme for Lagrangian hydrodynamics. The SGH RD scheme is based on the staggered finite element discretizations as in [Dobrev et al., SISC, 2012]. However, the advantage of the residual formulation over classical FEM approaches consists in the natural mass matrix diagonalization which allows one to avoid the solution of the linear system with the global sparse mass matrix while retaining the desired order of accuracy. This is achieved by using Bernstein polynomials as finite element shape functions and coupling the space discretization with the deferred correction type timestepping method. Moreover, it can be shown that for the Lagrangian formulation written in non-conservative form, our residual distribution scheme ensures the exact conservation of the mass, momentum and total energy. In this paper we also discuss construction of numerical viscosity approximations for the SGH RD scheme allowing to reduce the dissipation of the numerical solution. Thanks to the generic formulation of the staggered grid residual distribution scheme, it can be directly applied to both single- and multimaterial and multiphase models. Finally, we demonstrate computational results obtained with the proposed residual distribution scheme for several challenging test problems.
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.
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.
Non-Newtonian fluid flows, especially in three dimensions (3D), arise in numerous settings of interest to physics. Prior studies using the lattice Boltzmann method (LBM) of such flows have so far been limited to mainly to two dimensions and used less robust collision models. In this paper, we develop a new 3D cascaded LBM based on central moments and multiple relaxation times on a three-dimensional, nineteen velocity (D3Q19) lattice for simulation of generalized Newtonian (power law) fluid flows. The relaxation times of the second order moments are varied locally based on the local shear rate and parameterized by the consistency coefficient and the power law index of the nonlinear constitutive relation of the power law fluid. Numerical validation study of the 3D cascaded LBM for various benchmark problems, including the complex 3D non-Newtonian flow in a cubic cavity at different Reynolds numbers and power law index magnitudes encompassing shear thinning and shear thickening fluids, are presented. Furthermore, numerical stability comparisons of the proposed advanced LBM scheme against the LBM based on other collision models, such as the SRT model and MRT model based on raw moments, are made. Numerical results demonstrate the accuracy, second order grid convergence and significant improvements in stability of the 3D cascaded LBM for simulation of 3D non-Newtonian flows of power law fluids.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
