Three-dimensional simulation of unstable gravity-driven infiltration of water into a porous medium

Journal of Computational Physics

Volumn 238, Issue , 2013, Pages 217-239

  Gomez, Hector ^a   Cueto Felgueroso, Luis ^b   Juanes, Ruben ^b  

^a UNIVERSITY OF A CORUÑA   (Spain)

^b MASSACHUSETTS INSTITUTE OF TECHNOLOGY   (United States)

Author keywords

Gravity fingering; Infiltration; Isogeometric analysis; Phase field

Indexed keywords

COMPUTATION THEORY; FINITE ELEMENT METHOD; FLOW OF WATER; INFILTRATION; MATHEMATICAL OPERATORS; NUMERICAL METHODS; PERTURBATION TECHNIQUES; SOLUTE TRANSPORT;

DISCRETIZATIONS; GRAVITY FINGERING; HIGH-ORDER; HIGHER-ORDER; ISOGEOMETRIC ANALYSIS; PHASE FIELDS; POROUS MEDIUM; RICHARDS EQUATION; THREE DIMENSIONS; UNSATURATED FLOWS;

POROUS MATERIALS;

EID: 84873297512     PISSN: 00219991     EISSN: 10902716     Source Type: Journal    
DOI: 10.1016/j.jcp.2012.12.018     Document Type: Article

References (77)

1. Anderson D.M., McFadden G.B., Wheeler A.A. Diffuse-interface methods in fluid mechanics. Annu. Rev. Fluid Mech. 1998, 30:139-165.
2. Akkerman I., Bazilevs Y., Calo V.M., Hughes T.J.R., Hulshoff S. The role of continuity in residual-based variational multiscale modeling of turbulence. Comput. Mech. 2007, 41:371-378.
3. Barret J.W., Blowey J.F., Garcke H. Finite element approximation of the Cahn-Hilliard equation with degenerate mobility. SIAM J. Numer. Anal. 1999, 37:286-318.
4. Bazilevs Y., Calo V.M., Cottrell J.A., Hughes T.J.R., Reali A., Scovazzi G. Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows. Comput. Methods Appl. Mech. Eng. 2007, 197:173-201.
5. Bazilevs Y., Calo V.M., Cottrell J.A., Evans J.A., Hughes T.J.R., Lipton S., Scott M.A., Sederberg T.W. Isogeometric analysis using T-splines. Comput. Methods Appl. Mech. Eng. 2010, 199:229-263.
6. Bazilevs Y., Calo V.M., Hughes T.J.R., Zhang Y. Isogeometric fluid-structure interaction: theory, algorithms and computations. Computational Mechanics 2008, 43:3-37.
7. Bear J. Dynamics of Fluids in Porous Media 1972, Dover Publications.
8. Bearer E.L., Lowengrub J.S., Frieboes H.B., Chuang Y.L., Jin F., Wise S.M., Ferrari M., Agus D.B., Cristini V. Multiparameter computational modeling of tumor invasion. Cancer Res. 2009, 69:4493-4501.
9. Blowey J.F., Elliott C.M. The Cahn-Hilliard gradient theory for phase separation with non-smooth free energy. Part II: numerical analysis. Eur. J. Appl. Math. 1992, 3:147-149.
10. Borden M.J., Verhoosel C.V., Scott M.A., Hughes T.J.R., Landis C.M. A phase-field description of dynamic brittle fracture. ICES Rep. 2011.
11. Brenner S.C., Scott L.R. The Mathematical Theory of Finite Element Methods 1994, Springer-Verlag, New York.
12. Brooks R.H., Corey A.T. Properties of porous media affecting fluid flow. J. Irrig. Drain. Div. Am. Soc. Civ. Eng. 1966, 92:61-88.
13. Buffa A., Sangalli G., Vázquez R. Isogeometric analysis in electromagnetics: B-splines approximation. Comput. Methods Appl. Mech. Eng. 2010, 199:1143-1152.
14. Cahn J.W. On spinodal decomposition. Acta Met. 1961, 9:795-801.
15. Cahn J.W., Hilliard J.E. Free energy of a non-uniform system. I. Interfacial free energy. J. Chem. Phys. 1958, 28:258-267.
16. Cahn J.W., Hilliard J.E. Free energy of a non-uniform system. III. Nucleation in a two-component incompressible fluid. J. Chem. Phys. 1959, 31:688-699.
17. Chen L.Q. Phase-field models for microstructural evolution. Ann. Rev. Mater. Res. 2002, 32:113-140.
18. Chung J., Hulbert G.M. A time integration algorithm for structural dynamics with improved numerical dissipation: the generalized-α method. J. Appl. Mech. 1993, 60:371-375.
19. Cottrell J.A., Hughes T.J.R., Bazilevs Y. Isogeometric Analysis: Toward Integration of CAD and FEA 2009, Wiley.
20. Cottrell J.A., Reali A., Bazilevs Y., Hughes T.J.R. Isogeometric analysis of structural vibrations. Comput. Methods Appl. Mech. Eng. 2006, 195:5257-5296.
21. Cristini V., Li X., Lowengrub J.S., Wise S.M. Nonlinear simulations of solid tumor growth using a mixture model: invasion and branching. J. Math. Biol. 2009, 58:723-763.
22. Cueto-Felgueroso L., Juanes R. Nonlocal interface dynamics and pattern formation in gravity-driven unsaturated flow through porous media. Phys. Rev. Lett. 2008, 101:244504.
23. Cueto-Felgueroso L., Juanes R. A phase-field model of unsaturated flow. Water Resour. Res. 2009, 45:W10409.
24. Cueto-Felgueroso L., Juanes R. Stability analysis of a phase-field model of gravity-driven unsaturated flow through porous media. Phys. Rev. E 2009, 79:036301.
25. Cueto-Felgueroso L., Peraire J. A time-adaptive finite volume method for the Cahn-Hilliard and Kuramoto-Sivashinsky equations. J. Comput. Phys. 2009, 227:9985-10017.
26. Cueto-Felgueroso L., Juanes R. Adaptive rational spectral methods for the linear stability analysis of nonlinear fourth-order problems. J. Comput. Phys. 2009, 228:6536-6552.
27. DiCarlo D.A. Experimental measurements of saturation overshoot on infiltration. Water Resour. Res. 2004, 40:W04215.
28. Elliott C.M., French D.A., Milner F.A. A 2nd-order splitting method for the Cahn-Hilliard equation. Numer. Math. 1989, 54:575-590.
29. Elguedj T., Bazilevs Y., Calo V.M., Hughes T.J.R. B̄ and F̄ projection methods for nearly incompressible linear and non-linear elasticity and plasticity using higher-order NURBS elements. Comput. Methods Appl. Mech. Eng. 2008, 197:2732-2762.
30. Eliassi M., Glass R.J. On the continuum scale modeling of gravity-driven fingers in unsaturated porous media: the inadequacy of the Richards equation with standard monotonic constitutive relations and hysteretic equations of state. Water Resour. Res. 2001, 37:2019-2035.
31. Evans J.A., Bazilevs Y., BabuŠka I., Hughes T.J.R. n-widths, sup infs, and optimality ratios for the k-version of the isogeometric finite element method. Comput. Methods Appl. Mech. Eng. 2009, 198:1726-1741.
32. Feng X., Prohl A. Analysis of a fully discrete finite element method for phase field model and approximation of its sharp interface limits. Math. Comput. 2003, 73:541-547.
33. Furihata D. A stable and conservative finite difference scheme for the Cahn-Hilliard equation. Numer. Math. 2001, 87:675-699.
34. Fürst T., Vodák R., Ši{dotless}́r M., Bi{dotless}́l M. On the incompatibility of Richards' equation and finger-like infiltration in unsaturated homogeneous porous media. Water Resour. Res. 2009, 45:W03408.
35. Geiger S.L., Durnford D.S. Infiltration in homogeneous sands and a mechanistic model of unstable flow. Soil Sci. Soc. Am. J. 2000, 64:460-469.
36. Gelhar L.W., Axness C.L. Three-dimensional stochastic analysis of macrodispersion in aquifers. Water Resour. Res. 1983, 19:161-180.
37. Glass R.J., Cann S., King J., Baily N., Parlange J.-Y., Steenhuis T.S. Wetting front instability in unsaturated porous media: a three-dimensional study in initially dry sand. Transp. Porous Media 1990, 5:247-268.
38. Glass R.J., Parlange J.-Y., Steenhuis T.S. Immiscible displacement in porous media: stability analysis of three-dimensional axisymmetric disturbances with application to gravity-driven wetting front instability. Water Resour. Res. 1991, 27:1947-1956.
39. Glass R.J., Parlange J.-Y., Steenhuis T.S. Wetting front instability: 2. Experimental determination of relationships between system parameters and two-dimensional unstable flow field behaviour in initially dry porous media. Water Resources Research 1989, 25:1195-1207.
40. Gomez H., Calo V.M., Bazilevs Y., Hughes T.J.R. Isogeometric analysis of the Cahn-Hilliard phase-field model. Comput. Methods Appl. Mech. Eng. 2008, 197:4333-4352.
41. Gomez H., Hughes T.J.R. Provably unconditionally stable, second-order time-accurate, mixed variational methods for phase-field models. J. Comput. Phys. 2011, 230:5310-5327.
42. Gomez H., Hughes T.J.R., Nogueira X., Calo V.M. Isogeometric analysis of the isothermal Navier-Stokes-Korteweg equations. Comput. Methods Appl. Mech. Eng. 2010, 199:1828-1840.
43. Gomez H., París J. Numerical simulation of asymptotic states of the damped Kuramoto-Sivashinsky equation. Phys. Rev. E 2011, 83:046702.
44. He Y., Liu Y., Tang T. On large time-stepping methods for the Cahn-Hilliard equation. Appl. Numer. Math. 2007, 57:616-628.
45. E Hill D., Parlange J.-Y. Wetting front instability in layered soils. Soil Sci. Soc. Am. J. 1972, 36:697-702.
46. Hughes T.J.R., Reali A., Sangalli G. Duality and unified analysis of discrete approximations in structural dynamics and wave propagation: comparison of p-method finite elements with k-method NURBS. Comput. Methods Appl. Mech. Eng. 2008, 197:4104-4124.
47. Hughes T.J.R., Cottrell J.A., Bazilevs Y. Isogeometric analysis: CAD finite elements NURBS exact geometry and mesh refinement. Comput. Methods Appl. Mech. Eng. 2005, 194:4135-4195.
48. Huppert H.E. Flow and instability of a viscous current down a slope. Nature 1982, 300:427-429.
49. Jansen K.E., Whiting C.H., Hulbert G.M. A generalized-α method for integrating the filtered Navier-Stokes equations with a stabilized finite element method. Comput. Methods Appl. Mech. Eng. 1999, 190:305-319.
50. Kim J., Lowengrub J. Phase field modeling and simulation of three-phase flows. Interfaces Free Bound. 2005, 7:435-466.
51. Kim J., Kang K., Lowengrub J. Conservative multigrid methods for Cahn-Hilliard fluids. J. Comput. Phys. 2004, 193:511-543.
52. Kuramoto Y., Tsuzuki T. Persistent propagation of concentration waves in dissipative media far from thermal equilibrium. Prog. Theor. Phys. 1976, 55:356-369.
53. Lee H.-G., Lowengrub J.S., Goodman J. Modeling pinchoff and reconnection in a Hele-Shaw cell. II. Analysis and simulation in the nonlinear regime. Phys. Fluids 2002, 14:514-545.
54. Leverett M.C. Capillary behavior of porous solids. Trans. AIME 1941, 142:152-169.
55. Lipton S., Evans J.A., Bazilevs Y., Elguedj T., Hughes T.J.R. Robustness of isogeometric structural discretizations under severe mesh distortion. Comput. Methods Appl. Mech. Eng. 2010, 199:357-373.
56. Liu C., Shen J. A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method. Physica D 2003, 179:211-228.
57. Lowengrub J., Frieboes H.B., Jin F., Chuang Y.L., Li X., Macklin P., Wise S.M., Cristini V. Nonlinear modeling of cancer: bridging the gap between cells and tomours. Nonlinearity 2010, 23:1-91.
58. Lowengrub J., Truskinovsky L. Quasi-incompressible Cahn-Hilliard fluids and topological transitions. Proc. Roy. Soc. London Ser. A 1998, 454:2617-2654.
59. L Nieber J., Dautov R.Z., Egorov A.G., Sheshukov A.Y. Dynamic capillary pressure mechanism for gravity-driven flows; review and extension to very dry conditions. Transp. Porous Media 2005, 58:147-152.
60. Oden J.T., Hawkins A., Prudhomme S. General diffuse-interface theories and an approach to predictive tumor growth modeling. Math. Models Methods Appl. Sci. 2010, 20:477-517.
61. Piegl L., Tiller W. The NURBS book 1996, Springer.
62. Richards L.A. Capillary conduction of liquids through porous mediums. Physics 1931, 1:318-333.
63. Pons A.J., Karma A. Helical crack-front instability in mixed-mode fracture. Nature 2010, 464:85-89.
64. Rodriguez-Iturbe I., D'Odorico P., Porporato A., Ridolfi L. On the spatial and temporal links between vegetation, climate and soil moisture. Water Resources Research 1999, 35:3709-3722.
65. Rogers D.F. An Introduction to NURBS: With Historical Perspective 2000, Morgan Kaufman.
66. Saad Y., Schultz M.H. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 1986, 7:856-869.
67. Schillinger D., Ruess M., Zander N., Bazilevs Y., Düster A., Rank E. Small and large deformation analysis with the p- and the B-spline versions of the Finite Cell Method. Computational Mechanics 2012, 10.1007/s00466-012-0684-z.
68. Sivashinsky G.I. Nonlinear analysis of hydrodynamic instability in laminar flames. I. Derivation of basic equations. Acta Astron. 1977, 4:1177-1206.
69. C1 macroelements in adaptive finite element methods. Int. J. Numer. Methods Eng. 2007, 70:1076-1095.
70. Sun Z.Z. A second order accurate linearized difference scheme for the two-dimensional Cahn-Hilliard equation. Math. Comput. 1995, 64:1463-1471.
71. Teigen K.E., Song P., Lowengrub J., Voigt A. A diffuse-interface method for two-phase flows with soluble surfactants. J. Comput. Phys. 2011, 230:375-393.
72. Trefethen L.N., Embree M. Spectra and Pseudospectra: The Behavior of Nonnormal Matrices and Operators 2005, Princeton University Press.
73. Tullis B.P., Wright S.J. Wetting front instabilities: a three-dimensional experimental investigation. Transp. Porous Media 2007, 70:335-353.
74. Verhoosel C.V., Scott M.A., Hughes T.J.R., de Borst R. An isogeometric analysis approach to gradient damage models. Int. J. Numer. Methods Eng. 2011, 86:115-134.
75. Wise S.M., Lowengrub J.S., Cristini V. An adaptive multigrid algorithm for simulating solid tumor growth using mixture models. Math. Comput. Model. 2011, 53:1-20.
76. Ye X., Cheng X. The Fourier spectral method for the Cahn-Hilliard equation. Appl. Math. Comput. 2005, 171:345-357.
77. Zhu J., Chen L.-Q., Shen J., Tikare V. Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. Phys. Rev. E 1999, 60:3564-3572.