ID#: 211
Abstract Title: Numerical Simulation of Multiphase Reactive Boundary Layer Flow
Session Title: Multi-phase Combustion I
Session Date: 7/31/01
Session Start Time: 3:40 PM
Contributing Author: Lu, S.
Organization: University of Science and Technology of China
Country: People's Republic of China
Authors: Shou-Xiang Lu
Short Abstract: A comprehensive model of multiphase reactive boundary layer flow behind a shock wave is presented. Because of the multiphase reactions, the governing equations of the boundary layer contain source terms. As a result, the Leeds transformation which is valid to solve the problems of ordinary boundary layer can not be employed to solve the equations of the present boundary layer. Based on the view of the mass conservation, the gas can be divided into two parts. One part is the intrinsic part which is from the initial gas phase and another part is the extrinsic part which is generated during the process of multiphase reaction. And the concepts of extrinsic gas density and intrinsic gas density are defined to deal with the present boundary layer problem. The extrinsic gas density is the gas density increment due to the multiphase reaction in the boundary layer and the intrinsic gas density equals to the gas density of the boundary layer flow presuming without any chemical reaction. Therefore, the mass equation can be resolved to two mass equations, the extrinsic and intrinsic mass equations. According to the intrinsic mass equation, the flow function can be introduced as easily as the problem of ordinary boundary layer so that the Leeds-Dorodnitzyn transformation is valid for the present model. The sets of the equations adds one extrinsic equation after above transformed but the numerical solve for them become possible. In the present work, the numerical method for the laminar and turbulent boundary layer presented by Cebeci and Smith is employed. The Box scheme employed is of two order accuracy. To make a coupling of gas phase and solid phase equations, the method of Particle-Source-in Cell (PSIC) is modified because the PSIC method is based on the SIMPLER scheme which is constructed on the cells but BOX scheme on the junctions. The results by the present methods reveal the basic structure of the multiphase reactive boundary layer and give a good agreement with some of the experimental results. It is show that the present model is proper for the simulation of the phenomenon mentioned above.