Cantera  2.1.2
TortuosityMaxwell.h
1 /**
2  * @file TortuosityBase.h
3  * Virtual base class to compute the increase in diffusive path length associated with
4  * tortuous path diffusion through, for example, porous media.
5  */
6 
7 /*
8  * Copyright (2005) Sandia Corporation. Under the terms of
9  * Contract DE-AC04-94AL85000 with Sandia Corporation, the
10  * U.S. Government retains certain rights in this software.
11  */
12 
13 #ifndef CT_TORTUOSITYBRUGGEMAN_H
14 #define CT_TORTUOSITYBRUGGEMAN_H
15 
16 #include "TortuosityBase.h"
17 
18 namespace Cantera
19 {
20 
21 //! Maxwell model for tortuosity
22 /*!
23  *
24  * This class implements transport coefficient corrections
25  * appropriate for porous media with a dispersed phase.
26  * This model goes back to Maxwell. The formula for the
27  * conductivity is expressed in terms of the volume fraction
28  * of the continuous phase, \f$ \phi \f$, and the relative
29  * conductivities of the dispersed and continuous phases,
30  * \f$ r = \kappa_d / \kappa_0 \f$. For dilute particle
31  * suspensions the effective conductivity is
32  *
33  * \f[
34  * \kappa / \kappa_0 = 1 + 3 ( 1 - \phi ) ( r - 1 ) / ( r + 2 )
35  * + O(\phi^2)
36  * \f]
37  *
38  * The class is derived from the TortuosityBase class.
39  *
40  */
42 {
43 
44 public:
45  //! Default constructor uses Maxwelln exponent of 1.5
46  /*!
47  * @param setPower Exponent in the Maxwell factor. The default is 1.5
48  */
49  TortuosityMaxwell(double relativeConductivites = 0.0);
50 
51  //! Copy Constructor
52  /*!
53  * @param right Object to be copied
54  */
56 
57  //! Assignment operator
58  /*!
59  * @param right Object to be copied
60  */
62 
63  //! Duplication operator
64  /*!
65  * @return Returns a pointer to a duplicate of the current object given a
66  * base class pointer
67  */
69 
70  //! The tortuosity factor models the effective increase in the
71  //! diffusive transport length.
72  /*!
73  * This method returns \f$ 1/\tau^2 \f$ in the description of the flux
74  *
75  * \f$ C_T D_i \nabla X_i / \tau^2 \f$.
76  *
77  *
78  */
79  virtual doublereal tortuosityFactor(doublereal porosity);
80 
81  //! The McMillan number is the ratio of the flux-like
82  //! variable to the value it would have without porous flow.
83  /**
84  * The McMillan number combines the effect of tortuosity
85  * and volume fraction of the transported phase. The net flux
86  * observed is then the product of the McMillan number and the
87  * non-porous transport rate. For a conductivity in a non-porous
88  * media, \f$ \kappa_0 \f$, the conductivity in the porous media
89  * would be \f$ \kappa = (\rm McMillan) \kappa_0 \f$.
90  */
91  virtual doublereal McMillanFactor(doublereal porosity);
92 
93 
94 protected:
95 
96  //! Relative conductivities of the dispersed and continuous phases,
97  /*!
98  *
99  * \f[
100  * \code{relativeConductivites_} = \kappa_d / \kappa_0
101  * \f]
102  */
104 
105 };
106 
107 
108 
109 }
110 
111 #endif
112 
Virtual base class to compute the increase in diffusive path length associated with tortuous path dif...
doublereal relativeConductivities_
Relative conductivities of the dispersed and continuous phases,.
Base case to handle tortuosity corrections for diffusive transport in porous media.
virtual TortuosityBase * duplMyselfAsTortuosityBase() const
Duplication operator.
TortuosityMaxwell(double relativeConductivites=0.0)
Default constructor uses Maxwelln exponent of 1.5.
TortuosityMaxwell & operator=(const TortuosityMaxwell &right)
Assignment operator.
virtual doublereal tortuosityFactor(doublereal porosity)
The tortuosity factor models the effective increase in the diffusive transport length.
virtual doublereal McMillanFactor(doublereal porosity)
The McMillan number is the ratio of the flux-like variable to the value it would have without porous ...
Maxwell model for tortuosity.