Analytical and Numerical Methods for Solving Partial Differential Equations and Integral Equations Arising in Physical Models
View this Special IssueResearch Article  Open Access
Yong Cao, Yuchuan Chu, Xiaoming He, Mingzhen Wei, "Decoupling the Stationary NavierStokesDarcy System with the BeaversJosephSaffman Interface Condition", Abstract and Applied Analysis, vol. 2013, Article ID 136483, 10 pages, 2013. https://doi.org/10.1155/2013/136483
Decoupling the Stationary NavierStokesDarcy System with the BeaversJosephSaffman Interface Condition
Abstract
This paper proposes a domain decomposition method for the coupled stationary NavierStokes and Darcy equations with the BeaversJosephSaffman interface condition in order to improve the efficiency of the finite element method. The physical interface conditions are directly utilized to construct the boundary conditions on the interface and then decouple the NavierStokes and Darcy equations. Newton iteration will be used to deal with the nonlinear systems. Numerical results are presented to illustrate the features of the proposed method.
1. Introduction
The StokesDarcy model has been extensively studied in the recent years due to its wide range of applications in many natural world problems and industrial settings, such as the subsurface flow in karst aquifers, oil flow in vuggy porous media, industrial filtrations, and the interaction between surface and subsurface flows [1–8]. Since the problem domain naturally consists of two different physical subdomains, several different numerical methods have been developed to decouple the Stokes and Darcy equations [6, 9–26]. For other works on the numerical methods and analysis of the StokesDarcy model, we refer the readers to [27–45].
Recently the more physically valid NavierStokesDarcy model has attracted scientists’ attention, and several coupled finite element methods have been studied for it [46–51]. On the other hand, the advantages of the domain decomposition methods (DDMs) in parallel computation and natural preconditioning have motivated the development of different DDMs for solving the StokeDarcy model [6, 10–18, 21, 22]. In this paper, we will develop a domain decomposition method for the NavierStokesDarcy model based on Robin boundary conditions constructed from the interface conditions. This physicsbased DDM is different from the traditional ones in the sense that they focus on decomposing different physical domains by directly utilizing the given physical interface conditions.
The rest of paper is organized as follows. In Section 2, we introduce the NavierStokesDarcy model with the BeaversJosephSaffman interface condition. In Section 3, we recall the coupled weak formulation and the corresponding coupled finite element method for the NavierStokesDarcy model. In Section 4, a parallel domain decomposition method and its finite element discretization are proposed to decouple the NavierStokesDarcy system by using the Robintype boundary conditions constructed from the physical interface conditions. Finally, in Section 5, we present a numerical example to illustrate the features of the proposed method.
2. Stationary NavierStokesDarcy Model
In this section we introduce the following coupled NavierStokesDarcy model on a bounded domain , (); see Figure 1. In the porous media region , the flow is governed by the Darcy system Here, is the fluid discharge rate in the porous media, is the hydraulic conductivity tensor, is a sink/source term, and is the hydraulic head defined as , where denotes the dynamic pressure, the height, the density, and the gravitational acceleration. We will consider the following secondorder formulation, which eliminates in the Darcy system:
In the fluid region , the fluid flow is assumed to be governed by the NavierStokes equations: where is the fluid velocity, is the kinematic pressure, is the external body force, is the kinematic viscosity of the fluid, is the stress tensor, and is the deformation tensor.
Let denote the interface between the fluid and porous media regions. On the interface , we impose the following three interface conditions: where and denote the unit outer normal to the fluid and the porous media regions at the interface , respectively, denote mutually orthogonal unit tangential vectors to the interface , and . The third condition (7) is referred to as the BeaversJosephSaffman (BJS) interface condition [52–55].
In this paper, for simplification, we assume that the hydraulic head and the fluid velocity satisfy the homogeneous Dirichlet boundary condition except on , that is, on the boundary and on the boundary .
3. Coupled Weak Formulation and Finite Element Method
In this section we will recall the coupled weak formulation and the corresponding coupled finite element method for the NavierStokesDarcy model with BeaversJosephSaffman condition. Let denote the inner product on the domain ( or ) and let denote the inner product on the interface or the duality pairing between and [5]. Define the spaces the bilinear forms and the projection onto the tangent space on :
With these notations, the weak formulation of the coupled NavierStokesDarcy model with BJS interface condition is given as follows [46–51]: find such that
Assume that we have in hand regular subdivisions of and into finite elements with mesh size . Then one can define finite element spaces , and . We assume that and satisfy the infsup condition [56, 57] where is a constant independent of . This condition is needed in order to ensure that the spatial discretizations of the NavierStokes equations used here are stable. See [56, 57] for more details of finite element spaces , , and that satisfy (12). One example is the TaylorHood element pair that we use in the numerical experiments; for that pair, consists of continuous piecewise quadratic polynomials and consists of continuous piecewise linear polynomials.
Then a coupled finite element method with Newton iteration for the coupled NavierStokesDarcy model is given as follows [46]: find in the following procedure.(1)The initial value is chosen.(2)For , solve (3)Set , , and .
4. PhysicsBased Domain Decomposition Method
The coupled finite element method may end up with a huge algebraic system, which combines all parts from the NavierStokes equations, Darcy equation, and interface conditions together into one sparse matrix. Hence it is often impractical to directly apply this method to largescale real world applications. In order to develop a more efficient numerical method in this section, we will directly utilize the three physical interface conditions to construct a physicsbased parallel domain decomposition method to decouple the NavierStokes and Darcy equations.
Let us first consider the following Robin condition for the Darcy system: for a given constant and a given function defined on , Then, the corresponding weak formulation for the Darcy part is given by the following: for , find such that
Second, we can propose the following two Robintype conditions for the NavierStokes equations: for a given constant and given functions and defined on ,
Then, the corresponding weak formulation for the NavierStokes equation is given by the following: for , find such that
Our next step is to show that, for appropriate choices of , , , and , the solutions of the coupled system (11) are equivalent to the solutions of the decoupled equations (15) and (17), and hence we may solve the latter system instead of the former.
Lemma 1. Let be the solution of the coupled NavierStokesDarcy system (11) and let be the solution of the decoupled NavierStokes and Darcy equations (15) and (17) with Robin boundary conditions at the interface. Then, if and only if , , , , and satisfy the following compatibility conditions:
Proof. Adding (15) and (17) together, we obtain the following: given , find such that
For the necessity of the lemma, we pick and such that in (11) and (20); then by subtracting (20) from (11), we get
which implies (19). The necessity of (18) can be derived in a similar fashion.
As for the sufficiency of the lemma, by substituting the compatibility conditions (18)(19), we easily see that solves the coupled NavierStokesDarcy system (11), which completes the proof.
Now we use the decoupled weak formulation constructed above to propose a physicsbased parallel domain decomposition method with Newton iteration as follows. (1) Initial values and are guessed. They may be taken to be zero.(2) For , independently solve the Darcy and NavierStokes equations constructed above. More precisely, is computed from and and are computed from the following Newton iteration.(i) Initial value is chosen for the Newton iteration. For instance, it may be taken to be and for .(ii) For , solve (iii) Set and .(3) and are updated in the following manner:
Then the corresponding domain decomposition finite element method is proposed as follows.(1) Initial values and are guessed. They may be taken to be zero. (2) For , independently solve the Darcy and NavierStokes equations with the Robin boundary conditions on the interface, which are constructed previously. More precisely, is computed from and and are computed from the following Newton iteration.(i) Initial value is chosen for the Newton iteration. For instance, it may be taken to be and for .(ii) For , solve (iii) Set and .(3) and are updated in the following manner:
5. Numerical Example
Example 1. Consider the model problem (2)–(6) with the BJS interface condition (7) on with and . Choose , , , , and , where is the identity matrix and . The boundary condition data functions and the source terms are chosen such that the exact solution is given by We divide and into rectangles of height and width , where is a positive integer, and then subdivide each rectangle into two triangles by drawing a diagonal. The TaylorHood element pair is used for the NavierStokes equations, and the quadratic finite element is used for the secondorder formulation of the Darcy equation.
For the coupled finite element method of the steady NavierStokesDarcy model with BJS interface condition, Table 1 provides errors for different choices of . Using linear regression, the errors in Table 1 satisfy These rates of convergence are consistent with the approximation capability of the TaylorHood element and quadratic element, which is third order with respect to norm of and , second order with respect to norm of and , and second order with respect to norms of . In particular, the thirdorder convergence rate of observed above, which is better than the approximation capability of the linear element, is mainly due to the special analytic solution .
For the parallel DDM with , , , and , Figures 2 and 3 show the errors of hydraulic head, velocity, pressure, and . We can see that the parallel domain decomposition method is convergent for . Moreover, Figures 4 and 5 show that a smaller leads to faster convergence.
(a)
(b)
(a)
(b)
(a)
(b)
(a)
(b)
Then Tables 2 and 3 list some errors in velocity, hydraulic head, pressure, and for the parallel domain decomposition method with and . The data in these two tables indicate the geometric convergence rate since all the error ratios are less than .


Finally, for the preconditioning feature of the domain decomposition method, Table 4 shows the number of iterations is independent of the grid size . Here, we set , , , and . Let , , and denote the finite element solutions of , , and at the th step of the domain decomposition algorithm. The criterion used to stop the iteration, that is, to determine the value , is , where the tolerance .
6. Conclusions
In this paper, a parallel physicsbased domain decomposition method is proposed for the stationary NavierStokesDarcy model with the BJS interface condition. This method is based on the Robin boundary conditions constructed from the three physical interface conditions. Moreover, it is convergent with geometric convergence rates if the relaxation parameter is selected properly. The number of iteration steps is independent of the grid size due to the natural preconditioning advantage of the domain decomposition methods.
Acknowledgments
This work is partially supported by DOE Grant DEFE0009843, National Natural Science Foundation of China (11175052).
References
 T. Arbogast and D. S. Brunson, “A computational method for approximating a DarcyStokes system governing a vuggy porous medium,” Computational Geosciences, vol. 11, no. 3, pp. 207–218, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Arbogast and M. S. M. Gomez, “A discretization and multigrid solver for a DarcyStokes system of three dimensional vuggy porous media,” Computational Geosciences, vol. 13, no. 3, pp. 331–348, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 T. Arbogast and H. L. Lehr, “Homogenization of a DarcyStokes system modeling vuggy porous media,” Computational Geosciences, vol. 10, no. 3, pp. 291–302, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Cao, M. Gunzburger, X. Hu, F. Hua, X. Wang, and W. Zhao, “Finite element approximations for StokesDarcy flow with BeaversJoseph interface conditions,” SIAM Journal on Numerical Analysis, vol. 47, no. 6, pp. 4239–4256, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Cao, M. Gunzburger, F. Hua, and X. Wang, “Coupled StokesDarcy model with BeaversJoseph interface boundary condition,” Communications in Mathematical Sciences, vol. 8, no. 1, pp. 1–25, 2010. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Discacciati, Domain decomposition methods for the coupling of surface and groundwater flows [Ph.D. thesis], École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 2004.
 V. J. Ervin, E. W. Jenkins, and S. Sun, “Coupled generalized nonlinear Stokes flow with flow through a porous medium,” SIAM Journal on Numerical Analysis, vol. 47, no. 2, pp. 929–952, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 M. Moraiti, “On the quasistatic approximation in the StokesDarcy model of groundwatersurface water flows,” Journal of Mathematical Analysis and Applications, vol. 394, no. 2, pp. 796–808, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 I. Babuška and G. N. Gatica, “A residualbased a posteriori error estimator for the StokesDarcy coupled problem,” SIAM Journal on Numerical Analysis, vol. 48, no. 2, pp. 498–523, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Cao, M. Gunzburger, X. He, and X. Wang, “RobinRobin domain decomposition methods for the steadystate StokesDarcy system with the BeaversJoseph interface condition,” Numerische Mathematik, vol. 117, no. 4, pp. 601–629, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Cao, M. Gunzburger, X.M. He, and X. Wang, “Parallel, noniterative, multiphysics domain decomposition methods for timedependent StokesDarcy systems,” Mathematics of Computation. In press. View at: Google Scholar
 W. Chen, M. Gunzburger, F. Hua, and X. Wang, “A parallel RobinRobin domain decomposition method for the StokesDarcy system,” SIAM Journal on Numerical Analysis, vol. 49, no. 3, pp. 1064–1084, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M. Discacciati, “Iterative methods for Stokes/Darcy coupling,” in Domain Decomposition Methods in Science and Engineering, vol. 40 of Lecture Notes in Computational Science and Engineering, pp. 563–570, Springer, Berlin, Germany, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 M. Discacciati, E. Miglio, and A. Quarteroni, “Mathematical and numerical models for coupling surface and groundwater flows,” Applied Numerical Mathematics, vol. 43, no. 12, pp. 57–74, 2002, 19th Dundee Biennial Conference on Numerical Analysis (2001). View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Discacciati and A. Quarteroni, “Analysis of a domain decomposition method for the coupling of Stokes and Darcy equations,” in Numerical Mathematics and Advanced Applications, pp. 3–20, Springer, Milan, Italy, 2003. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Discacciati and A. Quarteroni, “Convergence analysis of a subdomain iterative method for the finite element approximation of the coupling of Stokes and Darcy equations,” Computing and Visualization in Science, vol. 6, no. 23, pp. 93–103, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 M. Discacciati, A. Quarteroni, and A. Valli, “RobinRobin domain decomposition methods for the StokesDarcy coupling,” SIAM Journal on Numerical Analysis, vol. 45, no. 3, pp. 1246–1268, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Feng, X. He, Z. Wang, and X. Zhang, “Noniterative domain decomposition methods for a nonstationary StokesDarcy model with BeaversJoseph interface condition,” Applied Mathematics and Computation, vol. 219, no. 2, pp. 453–463, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 G. N. Gatica, S. Meddahi, and R. Oyarzúa, “A conforming mixed finiteelement method for the coupling of fluid flow with porous media flow,” IMA Journal of Numerical Analysis, vol. 29, no. 1, pp. 86–108, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. N. Gatica, R. Oyarzúa, and F.J. Sayas, “A residualbased a posteriori error estimator for a fullymixed formulation of the StokesDarcy coupled problem,” Computer Methods in Applied Mechanics and Engineering, vol. 200, no. 2122, pp. 1877–1891, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. H. W. Hoppe, P. Porta, and Y. Vassilevski, “Computational issues related to iterative coupling of subsurface and channel flows,” Calcolo, vol. 44, no. 1, pp. 1–20, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Jiang, “A parallel domain decomposition method for coupling of surface and groundwater flows,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 9–12, pp. 947–957, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Layton, H. Tran, and X. Xiong, “Long time stability of four methods for splitting the evolutionary StokesDarcy problem into Stokes and Darcy subproblems,” Journal of Computational and Applied Mathematics, vol. 236, no. 13, pp. 3198–3217, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. J. Layton, F. Schieweck, and I. Yotov, “Coupling fluid flow with porous media flow,” SIAM Journal on Numerical Analysis, vol. 40, no. 6, pp. 2195–2218, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 L. Shan, H. Zheng, and W. J. Layton, “A decoupling method with different subdomain time steps for the nonstationary StokesDarcy model,” Numerical Methods for Partial Differential Equations, vol. 29, no. 2, pp. 549–583, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. Mu and X. Zhu, “Decoupled schemes for a nonstationary mixed StokesDarcy model,” Mathematics of Computation, vol. 79, no. 270, pp. 707–731, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 C. Bernardi, T. C. Rebollo, F. Hecht, and Z. Mghazli, “Mortar finite element discretization of a model coupling Darcy and Stokes equations,” Mathematical Modelling and Numerical Analysis, vol. 42, no. 3, pp. 375–410, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Boubendir and S. Tlupova, “StokesDarcy boundary integral solutions using preconditioners,” Journal of Computational Physics, vol. 228, no. 23, pp. 8627–8641, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 W. Chen, P. Chen, M. Gunzburger, and N. Yan, “Superconvergence analysis of FEMs for the StokesDarcy system,” Mathematical Methods in the Applied Sciences, vol. 33, no. 13, pp. 1605–1617, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 C. D'Angelo and P. Zunino, “Robust numerical approximation of coupled Stokes' and Darcy's flows applied to vascular hemodynamics and biochemical transport,” Mathematical Modelling and Numerical Analysis, vol. 45, no. 3, pp. 447–476, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 M.F. Feng, R.S. Qi, R. Zhu, and B.T. Ju, “Stabilized CrouzeixRaviart element for the coupled Stokes and Darcy problem,” Applied Mathematics and Mechanics, vol. 31, no. 3, pp. 393–404, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 J. Galvis and M. Sarkis, “Balancing domain decomposition methods for mortar coupling StokesDarcy systems,” in Domain Decomposition Methods in Science and Engineering XVI, vol. 55 of Lecture Notes in Computational Science and Engineering, pp. 373–380, Springer, Berlin, Germany, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 J. Galvis and M. Sarkis, “Nonmatching mortar discretization analysis for the coupling StokesDarcy equations,” Electronic Transactions on Numerical Analysis, vol. 26, pp. 350–384, 2007. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J. Galvis and M. Sarkis, “FETI and BDD preconditioners for StokesMortarDarcy systems,” Communications in Applied Mathematics and Computational Science, vol. 5, pp. 1–30, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. N. Gatica, R. Oyarzúa, and F.J. Sayas, “Convergence of a family of Galerkin discretizations for the StokesDarcy coupled problem,” Numerical Methods for Partial Differential Equations. An International Journal, vol. 27, no. 3, pp. 721–748, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 G. Kanschat and B. Rivière, “A strongly conservative finite element method for the coupling of Stokes and Darcy flow,” Journal of Computational Physics, vol. 229, no. 17, pp. 5933–5943, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 T. Karper, K.A. Mardal, and R. Winther, “Unified finite element discretizations of coupled DarcyStokes flow,” Numerical Methods for Partial Differential Equations, vol. 25, no. 2, pp. 311–326, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Khabthani, L. Elasmi, and F. Feuillebois, “Perturbation solution of the coupled StokesDarcy problem,” Discrete and Continuous Dynamical Systems B, vol. 15, no. 4, pp. 971–990, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Mu and J. Xu, “A twogrid method of a mixed StokesDarcy model for coupling fluid flow with porous media flow,” SIAM Journal on Numerical Analysis, vol. 45, no. 5, pp. 1801–1813, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Münzenmaier and G. Starke, “Firstorder system least squares for coupled StokesDarcy flow,” SIAM Journal on Numerical Analysis, vol. 49, no. 1, pp. 387–404, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 B. Rivière, “Analysis of a discontinuous finite element method for the coupled Stokes and Darcy problems,” Journal of Scientific Computing, vol. 2223, pp. 479–500, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Rivière and I. Yotov, “Locally conservative coupling of Stokes and Darcy flows,” SIAM Journal on Numerical Analysis, vol. 42, no. 5, pp. 1959–1977, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. Rui and R. Zhang, “A unified stabilized mixed finite element method for coupling Stokes and Darcy flows,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 33–36, pp. 2692–2699, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Tlupova and R. Cortez, “Boundary integral solutions of coupled Stokes and Darcy flows,” Journal of Computational Physics, vol. 228, no. 1, pp. 158–179, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. M. Urquiza, D. N'Dri, A. Garon, and M. C. Delfour, “Coupling Stokes and Darcy equations,” Applied Numerical Mathematics, vol. 58, no. 5, pp. 525–538, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. Badea, M. Discacciati, and A. Quarteroni, “Numerical analysis of the NavierStokes/Darcy coupling,” Numerische Mathematik, vol. 115, no. 2, pp. 195–227, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 M. Cai, M. Mu, and J. Xu, “Numerical solution to a mixed NavierStokes/Darcy model by the twogrid approach,” SIAM Journal on Numerical Analysis, vol. 47, no. 5, pp. 3325–3338, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Çeşmelioğlu and B. Rivière, “Analysis of timedependent NavierStokes flow coupled with Darcy flow,” Journal of Numerical Mathematics, vol. 16, no. 4, pp. 249–280, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Çeşmelioğlu and B. Rivière, “Primal discontinuous Galerkin methods for timedependent coupled surface and subsurface flow,” Journal of Scientific Computing, vol. 40, no. 1–3, pp. 115–140, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Chidyagwai and B. Rivière, “On the solution of the coupled NavierStokes and Darcy equations,” Computer Methods in Applied Mechanics and Engineering, vol. 198, no. 4748, pp. 3806–3820, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 V. Girault and B. Rivière, “DG approximation of coupled NavierStokes and Darcy equations by BeaverJosephSaffman interface condition,” SIAM Journal on Numerical Analysis, vol. 47, no. 3, pp. 2052–2089, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 G. Beavers and D. Joseph, “Boundary conditions at a naturally permeable wall,” Journal of Fluid Mechanics, vol. 30, pp. 197–207, 1967. View at: Google Scholar
 W. Jäger and A. Mikelić, “On the interface boundary condition of Beavers, Joseph, and Saffman,” SIAM Journal on Applied Mathematics, vol. 60, no. 4, pp. 1111–1127, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 I. P. Jones, “Low Reynolds number flow past a porous spherical shell,” Proceedings of the Cambridge Philosophical Society, vol. 73, pp. 231–238, 1973. View at: Google Scholar
 P. Saffman, “On the boundary condition at the interface of a porous medium,” Studies in Applied Mathematics, vol. 1, pp. 77–84, 1971. View at: Google Scholar
 V. Girault and P.A. Raviart, Finite Element Methods for NavierStokes Equations. Theory and Algorithms, vol. 5 of Springer Series in Computational Mathematics, Springer, Berlin, Germany, 1986. View at: Publisher Site  MathSciNet
 M. D. Gunzburger, Finite Element Methods for Viscous Incompressible Flows. A Guide to Theory, Practice, and Algorithms, Computer Science and Scientific Computing, Academic Press, Boston, Mass, USA, 1989. View at: MathSciNet
Copyright
Copyright © 2013 Yong Cao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.