A Numerical Model for Trickle Bed Reactors（滴流床反应器的数学模型）.pdf

Journal of Computational Physics 165, 311–333 (2000) doi:10.1006/jcph.2000.6604, available online at on A Numerical Model for Trickle Bed Reactors1 Richard M. Propp, Phillip Colella, William Y. Crutchfield, and Marcus S. Day Lawrence Berkeley National Laboratory, Berkeley, California 94720 E-mail: RMPropp@, PColella@, WYCrutchfield@, MSDay@ Received March 18, 1999; revised January 25, 2000 Trickle bed reactors are governed by equations of flow in porous media such as Darcy’s law and the conservation of mass. Our numerical method for solving these equations is based on a total-velocity splitting, sequential formulation which leads to an implicit pressure equation and a semi-implicit mass conservation equation. We use high-resolution finite-difference methods to discretize these equations. Our solution scheme extends previous work in modeling porous media flows in two ways. First, we incorporate physical effects due to capillary pressure, a nonlinear inlet boundary condition, spatial porosity variations, and inertial effects on phase mobilities. In particular, capillary forces introduce a parabolic component into the recast evolution equation, and the inertial effects give rise to hyperbolic nonconvexity. Second, we introduce a modification of the slope-limiting algorithm to prevent our numerical method from producing spurious shocks. We present a numerical algorithm for accommodating these difficulties, show the algorithm is second-order accurate, and demonstrate its performance on a number of simplified problems relevant to trickle bed reactor modeling. c 2000 Academic Press Key Words: trickle bed reactor; conservation laws; porous media flows;