11.07.2015 Views

Boreskov Institute of Catalysis of the Siberian Branch of Russian ...

Boreskov Institute of Catalysis of the Siberian Branch of Russian ...

Boreskov Institute of Catalysis of the Siberian Branch of Russian ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

OP-III-9mass dispersion and to analyze transversal symmetry breaking <strong>of</strong> moving and stationaryfronts using <strong>the</strong> relation between <strong>the</strong> velocity (V f ) and <strong>the</strong> local curvature (K) <strong>of</strong> <strong>the</strong> front.AnalysisAs a first stage we derived for a planar front two algebraic equations with two unknowns(∆T m ,V f ), which and in <strong>the</strong> limiting case <strong>of</strong> Pe C →∞ can be reduced to <strong>the</strong> approximaterelations known for a 1-D «ideal» front in <strong>the</strong> case <strong>of</strong> negligible mass dispersion 5 . Numericalsimulations <strong>of</strong> a 1-D PBR show a good agreement between <strong>the</strong> simulated and approximatedvalues for moderate and large Pe C .At <strong>the</strong> next stage we considered a curvilinear front propagation in a local polar coordinatesystem and derived approximate relations for <strong>the</strong> propagation velocity (V C f ) and <strong>the</strong> maximaltemperature rise for <strong>the</strong> case <strong>of</strong> a finite curvature. Applying <strong>the</strong> condition dV C f /dK| K=0 =0 weobtained a necessary condition for symmetry breaking in <strong>the</strong> following form:C(α,∆T ad ,Da,Pe T ,Pe C ,∆T C m ,V C f ,γ)[1- α] 0and <strong>the</strong> planar 1-D front is stable. Analysis <strong>of</strong> simplified Eq. (2) with α >1-δ (δ>0, small)shows that C>0 at least for 1< α 1 coincides with <strong>the</strong> previously obtained criterion (1), while refinedcondition (2) can define <strong>the</strong> upper α boundary <strong>of</strong> transversal patterns to emerge.SimulationsThe derived criterion was verified by direct numericalsimulations <strong>of</strong> <strong>the</strong> 2-D cylindrical shell and a full 3-D PBRmodels showing various types <strong>of</strong> moving transversal patternsin upstream propagating fronts (Fig. 1).In conclusion we would like to point that transversalthree-dimensional patterns were predicted and simulated inPBR's with a first order activated kinetics within <strong>the</strong> feasibledomain <strong>of</strong> parameters for <strong>the</strong> first time.Fig. 1.Typical quasi-«frozen»spiral structure <strong>of</strong> a frontpropagating in a 3-D PBRshowing <strong>the</strong> temperature patternin a cross-section normal to <strong>the</strong>flow (behind <strong>the</strong> front).References1. V. Balakotaiah, E.L. Christaforatou, & D.H. West, Chem. Eng. Sci., 54, 1725–1734 (1999).2. G. Viswanathan, A. Bindal, J. Khinast, & D. Luss, D, AIChE J., 51, 3028-3038 (2005).3. V. Balakotaiah, N. Gupta, & D.H. West, Chem. Eng. Sci., 57, 435–448 (2002).4. O. Nekhamkina, & M. Sheintuch, Chem. Eng. Sci, submitted.5. O.V. Kiselev, & Yu.Sh. Matros, Combustion, Explos. And Shock Waves, 16, 152 (1980).108

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!