In the present study, a numerical model based on the lattice Boltzmann method (LBM) is proposed to simulate multiphase mass transfer, referred to as the CST-LB model. This model introduced continuum species transfer (CST) formulation by an additional collision term to model the mass transfer across the multiphase interface. The boundary condition treatment of this model is also discussed. In order to verify the applicability, the CST-LB model is combined with the pseudopotential multiphase model to simulate a series of benchmark cases, including concentration jump near the interface, gas dissolution in a closed system, species transport during drainage in a capillary tube, and multiphase species transport in the porous media. This CST-LB model can also be coupled with other multiphase LBMs since the model depends on the phase fraction field, which is not explicitly limited to specified multiphase models.

As advancements in aerospace technology continue, reentry conditions pose increasingly rigorous requirements for thermal protection materials. Among these, fiber-reinforced phenolic resin composites have drawn substantial interest for their robust thermal insulation capabilities and enhanced ablation resistance, mechanical strength, and long-term reliability. This paper provides a comprehensive review of recent developments in fiber-reinforced phenolic resin composites, examining factors such as resin density, nanopore size within the matrix, resin cross-linking density, fiber–resin interfacial bonding, fiber length, fiber crystallization degree, and fiber weave structures that collectively influence composite performance. The typical applications of these composites in ultrahigh-speed aircraft are also discussed. Furthermore, the paper offers recommendations for future advancements focusing on optimizing mechanical, ablative, and insulation properties to meet the multifunctional demands of thermal protection materials.

Methane hydrate is a promising source of alternative energy. An in-depth understanding of the hydrate dissociation mechanism is crucial for the efficient extraction. In the present work, a comprehensive set of pore-scale numerical studies of hydrate dissociation mechanisms is presented. Pore-scale lattice Boltzmann (LB) models are proposed to simulate the multiphysics process during methane hydrate dissociation. The numerical simulations employ the actual hydrate sediment pore structure obtained by the micro-CT imaging. Experimental results of xenon hydrate dissociation are compared with the numerical simulations, indicating that the observed hydrate pore habits evolution is accurately captured by the proposed LB models. Furthermore, simulations of methane hydrate dissociation under different sediment water saturations, fluid flow rates and thermal conditions are conducted. Heat and mass transfer limitations both have significant effects on the methane hydrate dissociation rate. The bubble movement can further influence the dissociation process. Dissociation patterns can be divided into three categories, uniform, non-uniform and wormholing. The fluid flow impacts hydrate dissociation rates differently in three-dimensional real structures compared to two-dimensional idealized ones, influenced by variations in hydrate pore habits and flow properties. Finally, upscaling investigations are conducted to provide the permeability and kinetic models for the representative elementary volume (REV)-scale production forecast. Due to the difference in the hydrate pore habits and dissociation mechanisms, the three-dimensional upscaling results contrast with prior findings from two-dimensional studies. The present work provides a paradigm for pore-scale numerical simulation studies on the hydrate dissociation, which can offer theoretical guidance on efficient hydrate extraction.

Non-isothermal reactive transport in complicated porous media is diverse in nature and industrial applications. There are challenges in the modelling of multiple physicochemical processes in multiscale pore structures with various length scales ranging from nanometres to micrometres. This study focuses on coke combustion during in situ crude oil combustion techniques. A micro-continuum model was developed to perform an image-based simulation of coke combustion through a multiscale porous medium. The simulation coupled weakly compressible gas flow, species transport, conjugate heat transfer, heterogeneous coke oxidation kinetics and structural evolution. The unresolved nanoporous coke region was treated as a continuum, for which the random pore model, permeability model and species diffusivity model were integrated as sub-grid models to account for the sub-resolution reactive surface area, Darcy flow and Knudsen diffusion, respectively. A Pe–Da diagram was provided to present five characteristic combustion regimes covering the ignition temperature and air flux in realistic field operations and laboratory measurements. The present model proved to achieve more accurate predictions of the feasible ignition temperature than previous models. Compared with the air flux of $\phi \sim O\textrm{(1) s}{\textrm{m}^\textrm{3}}(\textrm{air})\;{({\textrm{m}^\textrm{2}}\ \textrm{h})^{ - 1}}$ in the field, the increasing air flux in the laboratory transformed the combustion regime from diffusion-limited to convection-limited, which led to an overpredicted burning temperature. Reactive fingering combustion was analysed to understand the potential risks in some experimental measurements. The findings provide a better understanding of coke combustion and can help engineers design sustainable combustion methods. The developed image-based model allows other types of multiscale and nonlinear reactive transport to be simulated.