Abstract:Jacobian clustering is proposed to reduce the computational cost associated with matrix operations en-countered in the Newton iteration in fully coupled, fully implicit schemes for unsteady reactive flow simulations with detailed chemistry. The iterative solver is based on the Lower-Upper Symmetric Gauss -Seidel (LUSGS) algorithm and sparse matrix technique. The evaluation and sparse Lower-Upper (LU) fac-torization of the diagonal block of the system Jacobian are performed within clusters rather than indi-vidual cells for Computational Fluid Dynamics (CFD) simulations. Cells that are close in the state space are clustered to provide the averaged states for calculating the Jacobians of chemical source terms and transport fluxes. The cells retrieve the factorized sparse matrices from the belonging clusters to perform the necessary iterations. For the purpose of clustering, the spatial dependency of transport Jacobian in the diagonal block is eliminated. To further reduce the computational cost, the sparsity of chemical Ja-cobian is augmented by removing the insignificant matrix elements. The method is tested in various one-dimensional hydrocarbon flames with both the second order Crank-Nicolson scheme and a third or-der implicit Runge-Kutta scheme. Various chemical mechanisms with 9 to 111 species are used to test the performance of the iterative solver. Fast convergence of Newton iteration is achieved, and the formal order of accuracy is demonstrated with Jacobian clustering. The overall costs of evaluation and factoriza-tion of the block diagonal Jacobian are negligible compared to the cost of calculating transport fluxes and chemical source terms. The averaged costs of Jacobian evaluation, LU factorization and Newton iteration, all increase only linearly with the number of chemical species. The fully coupled, fully implicit Crank-Nicolson scheme with Jacobian clustering shows 4 to 42 times speedup in computational time compared to the decoupled implicit scheme with Strang operator splitting. Jacobian clustering is promising to in-crease the computational efficiency of high order fully coupled, fully implicit schemes for unsteady reac-tive flow simulations with detailed chemistry.(c) 2022 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

A Fully Explicit Integrator for Modeling Astrophysical Reactive Flows

The Challenges of Modeling Astrophysical Reacting Flows

Numerical Simulation of Hypersonic Reaction Flows with Nonlinear Coupled Constitutive Relations

Efficient time stepping for reactive turbulent simulations with stiff chemistry

IMPLICIT-EXPLICIT RUNGE-KUTTA METHOD FOR COMBUSTION SIMULATION

A Spectral Radius Scaling Semi-Implicit Iterative Time Stepping Method for Reactive Flow Simulations with Detailed Chemistry

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

Krylov-based Adaptive-Rank Implicit Time Integrators for Stiff Problems with Application to Nonlinear Fokker-Planck Kinetic Models

An implicit adaptive unified gas-kinetic scheme for steady-state solutions of non-equilibrium flows

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker-Planck-Boltzmann Model

Application of the Explicitly Iterative Scheme to Simulating Subsonic Reacting Gas Flows

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker–Planck–Boltzmann Model

An Asymptotically-Correct Implicit-Explicit Time Integration Scheme for Finite Volume Radiation-Hydrodynamics

Lagrangian diffusive reactor for detailed thermochemical computations of plasma flows

Variable time-stepping exponential integrators for chemical reactors with analytical Jacobians

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: gray diffusion

A fully coupled, fully implicit simulation method for unsteady flames using Jacobian approximation and clustering

An implicit kinetic inviscid flux for predicting continuum flows in all speed regimes

A component-splitting implicit time integration for multicomponent reacting flows simulations

An asymptotically correct implicit–explicit time integration scheme for finite volume radiation-hydrodynamics

Mass-Flux-Based Implicit Multigrid Method for Modeling Multidimensional Combustion