Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMoment method as a numerical solver: Challenge from shock structure problems
393
0
0.0
(
0
)
Authors
Zhenning Cai
Ask ChatGPT about the research
No Arabic abstract
We survey a number of moment hierarchies and test their performances in computing one-dimensional shock structures. It is found that for high Mach numbers, the moment hierarchies are either computationally expensive or hard to converge, making these methods questionable for the simulation of highly non-equilibrium flows. By examining the convergence issue of Grads moment methods, we propose a new moment hierarchy to bridge the hydrodynamic models and the kinetic equation, allowing nonlinear moment methods to be used as a numerical tool to discretize the velocity space for high-speed flows. For the case of one-dimensional velocity, the method is formulated for odd number of moments, and it can be extended seamlessly to the three-dimensional case. Numerical tests show that the method is capable of predicting shock structures with high Mach numbers accurately, and the results converge to the solution of the Boltzmann equation as the number of moments increases. Some applications beyond the shock structure problem are also considered, indicating that the proposed method is suitable for the computation of transitional flows.
rate research
Read More
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.
Turbulence structure resulting from multi-fluid or multi-species, variable-density isotropic turbulence interaction with a Mach 2 shock is studied using turbulence-resolving shock-capturing simulations and Eulerian (grid) and Lagrangian (particle) methods. The complex roles density play in the modification of turbulence by the shock wave are identified. Statistical analyses of the velocity gradient tensor (VGT) show that the density variations significantly change the turbulence structure and flow topology. Specifically, a stronger symmetrization of the joint probability density function (PDF) of second and third invariants of the anisotropic velocity gradient tensor, PDF$(Q^ast, R^ast)$, as well as the PDF of the vortex stretching contribution to the enstrophy equation, are observed in the multi-species case. Furthermore, subsequent to the interaction with the shock, turbulent statistics also acquire a differential distribution in regions having different densities. This results in a nearly symmetrical PDF$(Q^ast, R^ast)$ in heavy fluid regions, while the light fluid regions retain the characteristic tear-drop shape. To understand this behavior and the return to standard turbulence structure as the flow evolves away from the shock, Lagrangian dynamics of the VGT and its invariants are studied by considering particle residence times and conditional particle variables in different flow regions. The pressure Hessian contributions to the VGT invariants transport equations are shown to be not only affected by the shock wave, but also by the density in the multi-fluid case, making them critically important to the flow dynamics and turbulence structure.
The performance of interFoam (a widely used solver within OpenFOAM package) in simulating the propagation of water waves has been reported to be sensitive to the temporal and spatial resolution. To facilitate more accurate simulations, a numerical wave tank is built based on a high-order accurate Navier-Stokes model, which employs the VPM (volume-average/point-value multi-moment) scheme as the fluid solver and the THINC/QQ method (THINC method with quadratic surface representation and Gaussian quadrature) for the free-surface capturing. Simulations of regular waves in an intermediate water depth are conducted and the results are assessed via comparing with the analytical solutions. The performance of the present model and interFoam solver in simulating the wave propagation is systematically compared in this work. The results clearly demonstrate that compared with interFoam solver, the present model significantly improves the dissipation properties of the propagating wave, where the waveforms as well as the velocity distribution can be substantially maintained while the waves propagating over long distances even with large time steps and coarse grids. It is also shown that the present model requires much less computation time to reach a given error level in comparison with interFoam solver.
Lattice Boltzmann (LB) models used for the computation of fluid flows represented by the Navier-Stokes (NS) equations on standard lattices can lead to non-Galilean invariant (GI) viscous stress involving cubic velocity errors. This arises from the dependence of their third order diagonal moments on the first order moments for standard lattices, and strategies have recently been introduced to restore GI without such errors using a modified collision operator involving either corrections to the relaxation times or to the moment equilibria. Convergence acceleration in the simulation of steady flows can be achieved by solving the preconditioned NS equations, which contain a preconditioning parameter that alleviates the numerical stiffness. In the present study, we present a GI formulation of the preconditioned cascaded central moment LB method used to solve the preconditioned NS equations, which is free of cubic velocity errors on a standard lattice. A Chapman-Enskog analysis reveals the structure of the spurious non-GI defect terms and it is demonstrated that the anisotropy of the resulting viscous stress is dependent on the preconditioning parameter, in addition to the fluid velocity. It is shown that partial correction to eliminate the cubic velocity defects is achieved by scaling the cubic velocity terms in the off-diagonal third-order moment equilibria with the square of the preconditioning parameter. Furthermore, we develop additional corrections based on the extended moment equilibria involving gradient terms with coefficients dependent locally on the fluid velocity and the preconditioning parameter. Several conclusions are drawn from the analysis of the structure of the non-GI errors and the associated corrections, with particular emphasis on their dependence on the preconditioning parameter. Improvements in accuracy and convergence acceleration are demonstrated.
In this numerical study, an original approach to simulate non-isothermal viscoelastic fluid flows at high Weissenberg numbers is presented. Stable computations over a wide range of Weissenberg numbers are assured by using the root conformation approach in a finite volume framework on general unstructured meshes. The numerical stabilization framework is extended to consider thermo-rheological properties in Oldroyd-B type viscoelastic fluids. The temperature dependence of the viscoelastic fluid is modeled with the time-temperature superposition principle. Both Arrhenius and WLF shift factors can be chosen, depending on the flow characteristics. The internal energy balance takes into account both energy and entropy elasticity. Partitioning is achieved by a constant split factor. An analytical solution of the balance equations in planar channel flow is derived to verify the results of the main field variables and to estimate the numerical error. The more complex entry flow of a polyisobutylene-based polymer solution in an axisymmetric 4:1 contraction is studied and compared to experimental data from the literature. We demonstrate the stability of the method in the experimentally relevant range of high Weissenberg numbers. The results at different imposed wall temperatures, as well as Weissenberg numbers, are found to be in good agreement with experimental data. Furthermore, the division between energy and entropy elasticity is investigated in detail with regard to the experimental setup.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
