Search

In this study, a species-clustered ordinary differential equations (ODE) solver for chemical kinetics with large detailed mechanisms based on operator-splitting is presented. The ODE system is split into clusters of species by using graph partition methods which has been intensively studied in areas of model reduction, parameterization and coarse-graining, etc. , such as diffusion maps based on the concept of Markov random walk. Definition of the weight (similarity) matrix is application-driven and according to chemical kinetics. Each cluster of species is then integrated by VODE, an implicit solver which is intractable and costly for large systems of many species and reactions. Expected speedup in computational efficiency is observed by numerical experiments on three zero-dimensional (0D) auto-ignition problems, considering the detailed hydrocarbon/air combustion mechanisms in varying scales, from 53 species with 325 reactions of methane to 2115 species with 8157 reactions of n-hexadecane., Comment: 28 pages and 14 figures

In this paper we develop a Lagrangian Inertial Centroidal Voronoi Particle (LICVP) method to extend the original CVP method \cite{fu2017physics} to dynamic load balancing in particle-based simulations. Two new concepts are proposed to address the additional problems encountered in repartitioning the system. First, a background velocity is introduced to transport Voronoi particle according to the local fluid field, which facilitates data reuse and lower data redistribution cost during rebalancing. Second, in order to handle problems with skew-aligned computational load and large void space, we develop an inertial-based partitioning strategy, where the inertial matrix is utilized to characterize the load distribution, and to confine the motion of Voronoi particles dynamically adapting to the physical simulation. Intensive numerical tests in fluid dynamics simulations reveal that the underlying LICVP method improves the incremental property remarkably without sacrifices on other objectives, i.e. the inter-processor communication is optimized simultaneously, and the repartitioning procedure is highly efficient., Comment: 31 pages, 12 figures

In this paper, a new fractional step method is proposed for simulating stiff and nonstiff chemically reacting flows. In stiff cases, a well-known spurious numerical phenomenon, i.e. the incorrect propagation speed of discontinuities, may be produced by general fractional step methods due to the under-resolved discretization in both space and time. The previous random projection method has been successfully applied for stiff detonation capturing in under-resolved conditions. Not to randomly project the intermediate state into two presumed equilibrium states (completely burnt or unburnt) as in the random projection method, the present study is to randomly choose the time-dependent advance or stop of a reaction process. Each one-way reaction has been decoupled from the multi-reaction kinetics using operator splitting and the local smeared temperature due to numerical dissipation of shock-capturing schemes is compared with a random one within two limited temperatures corresponding to the advance and its inverse states, respectively, to control the random reaction. The random activation or deactivation in the reaction step is thus promising to correct the deterministic accumulative error of the propagation of discontinuities. Extensive numerical experiments, including model problems and realistic reacting flows in one and two dimensions, demonstrate this expectation as well as the effectiveness and robustness of the method. Meanwhile, for nonstiff problems when spatial and temporal resolutions are fine, the proposed random method recovers the results as general fractional step methods, owing to the increasing possibility of activation with diminishing randomness by adding a shift term., Comment: 48 pages, 14 figures

An incremental-stencil WENO reconstruction method, which uses low-order candidate stencils with incrementally increasing width, is proposed for finite-volume simulation of compressible two-phase flow with the quasi-conservative interface model. While recovering the original 5th-order WENO reconstruction in smooth region of the solution, due to the presence of 2-point candidate stencils, the present method is able to handle closely located discontinuities, which is a typical scenario of shock-interface interaction. Furthermore, a MOOD-type positivity preserving approach is applied to ensure physical meaningful reconstruction. Compared with the hybrid method which switches between with the 5th-order WENO and 2nd-order reconstructions, the present method is free of problem-dependent tunable parameters. A number of numerical examples show that the present method achieves small numerical dissipation and good robustness for simulating two-phase flow problems with strong shock-interface interaction and large density ratio., Comment: 46 pages, 19 figures