The effect of numerical methods on the simulation of mid-ocean ridge hydrothermal models
Tóm tắt
This work considers the effect of the numerical method on the simulation of a 2D model of hydrothermal systems located in the high-permeability axial plane of mid-ocean ridges. The behavior of hot plumes, formed in a porous medium between volcanic lava and the ocean floor, is very irregular due to convective instabilities. Therefore, we discuss and compare two different numerical methods for solving the mathematical model of this system. In concrete, we consider two ways to treat the temperature equation of the model: a semi-Lagrangian formulation of the advective terms in combination with a Galerkin finite element method for the parabolic part of the equations and a stabilized finite element scheme. Both methods are very robust and accurate. However, due to physical instabilities in the system at high Rayleigh number, the effect of the numerical method is significant with regard to the temperature distribution at a certain time instant. The good news is that relevant statistical quantities remain relatively stable and coincide for the two numerical schemes. The agreement is larger in the case of a mathematical model with constant water properties. In the case of a model with nonlinear dependence of the water properties on the temperature and pressure, the agreement in the statistics is clearly less pronounced. Hence, the presented work accentuates the need for a strengthened validation of the compatibility between numerical scheme (accuracy/resolution) and complex (realistic/nonlinear) models.
Tài liệu tham khảo
Cherkaoui A.S.M., Wilcock W.S.D.: Characteristics of high Rayleigh number two-dimensional convection in an open-top porous layer heated from below. J. Fluid Mech. 394, 241–260 (1999)
Jupp T., Schultz A.: A thermodynamic explanation for black smoker temperatures. Nature 403(6772), 880–883 (2000)
German, C.R., Lin, J.: The thermal structure of the oceanic crust, ridge-spreading and hydrothermal circulation: how well do we understand their inter-connections? In: German, C.R., Lin, J., Parson, L.M. (eds.) Mid-ocean Ridges: Hydrothermal Interactions between the Lithosphere and Oceans, Geophys. Monogr., vol. 148, pp. 1–18 (2005)
Coumou D., Driesner T., Geiger S., Heinrich C.A., Matthäi S.: The dynamics of mid-ocean ridge hydrothermal systems: splitting plumes rand fluctuating vent temperatures. Earth Planet. Sci. Lett. 245, 218–231 (2006)
Driesner, T., Heinrich, C.A.: Accurate P–T–X–V–H correlations for the system NaCl–H2O from 0 to 800°C, 0 to 500MPa, and 0 to 1 XNaCl, Acta Mineralog. Petrograph. Abstr. Ser. 2, pp. 55–56 (2003)
Driesner, T., Heinrich, C.A.: The system NaCl–H2O. I. Correlation formulae for phase relations in temperature-pressure-composition space from 0 to 1000°C, 0 to 5000 bar, and 0 to 1000°C, 0 to 5000 bar, and 0 to 1 XNaCl, Geochim. Cosmochim. Acta (2005)
Driesner, T.: The system NaCl–H2O... II. Molar volume, enthalpy, and isobaric heat capacity from 0 to 1000°C, 0 to 5000 bar, and 0 to 1000°C, 0 to 5000 bar, and 0 to 1 XNaCl, Geochim. Cosmochim. Acta (2005)
Haar, L., Gallagher, J.S., Kell, G.S.: NBS/NRC Steam Tables. Hemisphere Publishing Corporation, Washington, DC (1984). (www.iapws.org)
Bauer, O.: PROST 4.1, PROperties of Water and STeam. Technical University of Hamburg-Harburg, Germany. http://www.enersim.de/download/PROST4.html (2002)
Wagner W., Pruss A.: The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use. J. Phys. Chem. Ref. Data 31(2), 387–535 (2002)
Méheust Y., Schmittbuhl J.: Flow enhancement of a rough fracture. Geophys. Res. Lett. 27(18), 2989–2992 (2000)
Méheust, Y., Schmittbuhl, J.: Geometrical heterogeneities and permeability anisotropy of rough fractures. J. Geophys. Res. 106(B2), 2089–2102
Trangenstein J.A., Bell J.B.: Mathematical structure of compositional reservoir simulation. SIAM J. Sci. Stat. Comput. 10, 817–845 (1989)
Hughes T.J.R., Brooks A.N.: Streamline upwind Petrov Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equation. Comput. Methods Appl. Mech. Eng. 32, 199–259 (1982)
Johnson C., Saranen J.: Streamline diffusion methods for the incompressible Euler and Navier-Stokes equations. Math. Comp. 47, 1–18 (1986)
Allievi A., Bermejo R.: Finite element modified method of characteristics for the Navier-Stokes equations. Int. J. Numer. Methods Fluids 32, 439–464 (2000)
Schilders, W.H.A., ter Maten, E.J.W. (eds.): Handbook of Numerical Analysis. vol. XIII. Special volume: Numerical Methods in Electromagnetics. Handbook of Numerical Analysis 13. Elsevier/North Holland, Amsterdam (2005)
Roos, H.-G., Stynes, M., Tobiska, L.: Robust Numerical Methods for Singularly Perturbed Differential Equations. Springer Series in Comp. Math., (2008)
Douglas J., Russell T.F.: Numerical methods for convection dominated diffusion problems based on combining the method of the characteristics with finite element of finite difference procedures. SIAM J. Numer. Anal. 19, 871–885 (1982)
Pironneau O.: On the transport diffusion algorithm and its applications to the Navier-Stokes equations. Numer. Math. 38, 309–332 (1982)
Benque, J.P., Labadie, G., Ronat, J.: A new finite element method for the Navier-Stokes equations coupled with a temperature equation. In: Kawai, T. (ed.) Proceedings of 4th International Symposium on Finite Element Methods in Flow Problems, pp. 295–301. North-Holland, Amsterdam (1982)
Robert A.: A stable numerical integration scheme for the primitive meteorological equations. Atmos-Ocean. 19, 35–46 (1981)
Bermejo, R., Carpio, J.: A semi-lagrangian-Galerkin projection scheme for convection equations. IMA J. Numer. Anal. doi:10.1093/imanum/drn044
Aziz K., Settari A.: Petroleum Reservoir Simulation. Applied Science Publishers, Barking (1979)
Geiger S., Driesner T., Matthai S.K., Heinrich C.A.: Multiphase thermohaline convection in the Earth’s crust I. A New Finite Element-Finite Volume solution technique combined with a new equation of state for NaCl−H2O. Transport Porous Med 63, 399–434 (2006)
Haymon, R.M.: The response of ridge-crest hydrothermal systems to segmented episodic magma supply. In: MacLeod, C.J., Taylor, P.A., Walker, C.L. (eds.) Tectonic; Magmatic; Hydrothermal; and Biological Segmentation of Mid-Ocean Ridges, pp. 157–168 (1996)
Kelley S.K., Baross J.A., Delaney J.R.: Volcanoes, fluids, and life at mid-ocean ridge speading centers. Annu. Rev. Earth Planet. Sci. 30, 385–491 (2002)
Schoofs S., Hansen U.: Depletion of a brine layer at the base of the ridge-crest hydrothermal system. Earth Planet Sci. Lett. 180(3-4), 341–353 (2000)
Otero J., Lubomira A.D., Johnston H., Worthing R.A., Kurganov A., Petrova G., Doering C.R.: High-Rayleigh-number convection in a fluid-saturated porous layer. J. Fluid Mech. 500, 263–281 (2004)
Tan C.T., Homsy G.M.: Stability of miscible displacements in porous media: rectilinear flow. Phys. Fluids 31(6), 3549–3556 (1988)