Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy

Mohammed Ahmadi

    Research output: Chapter in Book/Report/Conference proceedingConference contribution

    Abstract

    Numerical simulation of thermal recovery processes like steam injection often involves localized phenomena such as saturation and temperature fronts due to hyperbolic features of governing conservation laws, Treating more efficiently convective terms could help to diminish spurious oscillation and/or numerical dispersion and better tracking of discontinuity shocks. But in regions near the shock numerical dispersion can only be removed by the use of very fine uniform grids with many grid blocks. To avoid expensive solution of such a finely girded domain, we develop a moving mesh approach combined with higher order up-winding schemes. Numerical solver here have been employed is Finite volume method. A MMPDE(moving mesh PDE) is solved associated with physical PDE's of steam injection process in order to relocates the mesh nodes to concentrate them in regions of sharp discontinuity and Equi-distribute a measure of error-estimate (monitor function) over the meshes. Solution will advance more rapidly on course meshes and fluxes at the coarse-fine grid interfaces are refined to guarantee mass conservation. Since the region surrounding the sharp discontinuity and requiring high resolution consists of only a small fraction of the entire domain, prescribed locally time stepping results in a great saving in computational time. Specific features of moving mesh methods like monitor-function smoothing, control of mesh widths and readjustment of solutions further to mesh movement are addressed. Numerical experiments are carried out to demonstrate the efficiency and robustness of the proposed method in 1-D and 2-D. However numerical results for moving coordinates are compared with those obtained from simulation on non-adapted mesh framework. Preferences of higher-order solvers over lower-order ones in terms of shock capturing is being investigated. Although we have limited our modeling to steam flooding process, but simulation demonstrates main features of our approach, applicable to other EOR processes such as VAPEX, SAGD, and In Situ Combustion Process. Copyright 2008, SPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium.

    Original languageEnglish
    Title of host publicationSociety of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces"
    Pages37-57
    Number of pages21
    Volume1
    Publication statusPublished - 2008
    EventSPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium - Calgary, Alberta, Canada
    Duration: 20 Oct 200823 Oct 2008

    Conference

    ConferenceSPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium
    CountryCanada
    CityCalgary, Alberta
    Period20/10/0823/10/08

    Fingerprint

    Thermal oil recovery
    Finite volume method
    Steam
    Conservation
    In situ combustion
    Crude oil
    Fluxes
    Recovery
    Computer simulation
    Experiments
    Temperature
    Hot Temperature

    Cite this

    Ahmadi, M. (2008). Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy. In Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces" (Vol. 1, pp. 37-57)
    Ahmadi, Mohammed. / Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy. Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces". Vol. 1 2008. pp. 37-57
    @inproceedings{d1a399273e364eee855e747f3e1703f3,
    title = "Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy",
    abstract = "Numerical simulation of thermal recovery processes like steam injection often involves localized phenomena such as saturation and temperature fronts due to hyperbolic features of governing conservation laws, Treating more efficiently convective terms could help to diminish spurious oscillation and/or numerical dispersion and better tracking of discontinuity shocks. But in regions near the shock numerical dispersion can only be removed by the use of very fine uniform grids with many grid blocks. To avoid expensive solution of such a finely girded domain, we develop a moving mesh approach combined with higher order up-winding schemes. Numerical solver here have been employed is Finite volume method. A MMPDE(moving mesh PDE) is solved associated with physical PDE's of steam injection process in order to relocates the mesh nodes to concentrate them in regions of sharp discontinuity and Equi-distribute a measure of error-estimate (monitor function) over the meshes. Solution will advance more rapidly on course meshes and fluxes at the coarse-fine grid interfaces are refined to guarantee mass conservation. Since the region surrounding the sharp discontinuity and requiring high resolution consists of only a small fraction of the entire domain, prescribed locally time stepping results in a great saving in computational time. Specific features of moving mesh methods like monitor-function smoothing, control of mesh widths and readjustment of solutions further to mesh movement are addressed. Numerical experiments are carried out to demonstrate the efficiency and robustness of the proposed method in 1-D and 2-D. However numerical results for moving coordinates are compared with those obtained from simulation on non-adapted mesh framework. Preferences of higher-order solvers over lower-order ones in terms of shock capturing is being investigated. Although we have limited our modeling to steam flooding process, but simulation demonstrates main features of our approach, applicable to other EOR processes such as VAPEX, SAGD, and In Situ Combustion Process. Copyright 2008, SPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium.",
    author = "Mohammed Ahmadi",
    year = "2008",
    language = "English",
    isbn = "9781605606767",
    volume = "1",
    pages = "37--57",
    booktitle = "Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - {"}Heavy Oil: Integrating the Pieces{"}",

    }

    Ahmadi, M 2008, Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy. in Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces". vol. 1, pp. 37-57, SPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium, Calgary, Alberta, Canada, 20/10/08.

    Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy. / Ahmadi, Mohammed.

    Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces". Vol. 1 2008. p. 37-57.

    Research output: Chapter in Book/Report/Conference proceedingConference contribution

    TY - GEN

    T1 - Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy

    AU - Ahmadi, Mohammed

    PY - 2008

    Y1 - 2008

    N2 - Numerical simulation of thermal recovery processes like steam injection often involves localized phenomena such as saturation and temperature fronts due to hyperbolic features of governing conservation laws, Treating more efficiently convective terms could help to diminish spurious oscillation and/or numerical dispersion and better tracking of discontinuity shocks. But in regions near the shock numerical dispersion can only be removed by the use of very fine uniform grids with many grid blocks. To avoid expensive solution of such a finely girded domain, we develop a moving mesh approach combined with higher order up-winding schemes. Numerical solver here have been employed is Finite volume method. A MMPDE(moving mesh PDE) is solved associated with physical PDE's of steam injection process in order to relocates the mesh nodes to concentrate them in regions of sharp discontinuity and Equi-distribute a measure of error-estimate (monitor function) over the meshes. Solution will advance more rapidly on course meshes and fluxes at the coarse-fine grid interfaces are refined to guarantee mass conservation. Since the region surrounding the sharp discontinuity and requiring high resolution consists of only a small fraction of the entire domain, prescribed locally time stepping results in a great saving in computational time. Specific features of moving mesh methods like monitor-function smoothing, control of mesh widths and readjustment of solutions further to mesh movement are addressed. Numerical experiments are carried out to demonstrate the efficiency and robustness of the proposed method in 1-D and 2-D. However numerical results for moving coordinates are compared with those obtained from simulation on non-adapted mesh framework. Preferences of higher-order solvers over lower-order ones in terms of shock capturing is being investigated. Although we have limited our modeling to steam flooding process, but simulation demonstrates main features of our approach, applicable to other EOR processes such as VAPEX, SAGD, and In Situ Combustion Process. Copyright 2008, SPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium.

    AB - Numerical simulation of thermal recovery processes like steam injection often involves localized phenomena such as saturation and temperature fronts due to hyperbolic features of governing conservation laws, Treating more efficiently convective terms could help to diminish spurious oscillation and/or numerical dispersion and better tracking of discontinuity shocks. But in regions near the shock numerical dispersion can only be removed by the use of very fine uniform grids with many grid blocks. To avoid expensive solution of such a finely girded domain, we develop a moving mesh approach combined with higher order up-winding schemes. Numerical solver here have been employed is Finite volume method. A MMPDE(moving mesh PDE) is solved associated with physical PDE's of steam injection process in order to relocates the mesh nodes to concentrate them in regions of sharp discontinuity and Equi-distribute a measure of error-estimate (monitor function) over the meshes. Solution will advance more rapidly on course meshes and fluxes at the coarse-fine grid interfaces are refined to guarantee mass conservation. Since the region surrounding the sharp discontinuity and requiring high resolution consists of only a small fraction of the entire domain, prescribed locally time stepping results in a great saving in computational time. Specific features of moving mesh methods like monitor-function smoothing, control of mesh widths and readjustment of solutions further to mesh movement are addressed. Numerical experiments are carried out to demonstrate the efficiency and robustness of the proposed method in 1-D and 2-D. However numerical results for moving coordinates are compared with those obtained from simulation on non-adapted mesh framework. Preferences of higher-order solvers over lower-order ones in terms of shock capturing is being investigated. Although we have limited our modeling to steam flooding process, but simulation demonstrates main features of our approach, applicable to other EOR processes such as VAPEX, SAGD, and In Situ Combustion Process. Copyright 2008, SPE/PS/CHOA International Thermal Operations and Heavy Oil Symposium.

    UR - http://www.scopus.com/inward/record.url?scp=63849177648&partnerID=8YFLogxK

    M3 - Conference contribution

    SN - 9781605606767

    VL - 1

    SP - 37

    EP - 57

    BT - Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces"

    ER -

    Ahmadi M. Development of a higher-order finite volume method for simulation of thermal oil recovery process using moving mesh strategy. In Society of Petroleum Engineers - International Thermal Operations and Heavy Oil Symposium, ITOHOS 2008 - "Heavy Oil: Integrating the Pieces". Vol. 1. 2008. p. 37-57