Journal article 1259 views 315 downloads
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model
Raheel Ahmed,
Michael G. Edwards,
Sadok Lamine,
Bastiaan Huisman,
Mayur Pal
Journal of Computational Physics, Volume: 303, Pages: 470 - 497
Swansea University Author: Michael G. Edwards
-
PDF | Accepted Manuscript
Download (6.63MB)
DOI (Published version): 10.1016/j.jcp.2015.10.001
Abstract
A novel cell-centred control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulation is presented for discrete fracture–matrix simulations on unstructured grids in three-dimensions (3D). The grid is aligned with fractures and barriers which are then modelled as lower-d...
Published in: | Journal of Computational Physics |
---|---|
ISSN: | 0021-9991 |
Published: |
2015
|
Online Access: |
Check full text
|
URI: | https://cronfa.swan.ac.uk/Record/cronfa24695 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
first_indexed |
2015-11-24T02:01:17Z |
---|---|
last_indexed |
2021-01-07T03:39:38Z |
id |
cronfa24695 |
recordtype |
SURis |
fullrecord |
<?xml version="1.0"?><rfc1807><datestamp>2021-01-06T12:17:40.0548488</datestamp><bib-version>v2</bib-version><id>24695</id><entry>2015-11-23</entry><title>Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model</title><swanseaauthors><author><sid>8903caf3d43fca03602a72ed31d17c59</sid><firstname>Michael G.</firstname><surname>Edwards</surname><name>Michael G. Edwards</name><active>true</active><ethesisStudent>false</ethesisStudent></author></swanseaauthors><date>2015-11-23</date><deptcode>FGSEN</deptcode><abstract>A novel cell-centred control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulation is presented for discrete fracture–matrix simulations on unstructured grids in three-dimensions (3D). The grid is aligned with fractures and barriers which are then modelled as lower-dimensional surface interfaces located between the matrix cells in the physical domain. The three-dimensional pressure equation is solved in the matrix domain coupled with a two-dimensional (2D) surface pressure equation solved over fracture networks via a novel surface CVD-MPFA formulation. The CVD-MPFA formulation naturally handles fractures with anisotropic permeabilities on unstructured grids. Matrix–fracture fluxes are expressed in terms of matrix and fracture pressures and define the transfer function, which is added to the lower-dimensional flow equation and couples the three-dimensional and surface systems. An additional transmission condition is used between matrix cells adjacent to low permeable fractures to couple the velocity and pressure jump across the fractures. Convergence and accuracy of the lower-dimensional fracture model is assessed for highly anisotropic fractures having a range of apertures and permeability tensors. A transport equation for tracer flow is coupled via the Darcy flux for single and intersecting fractures. The lower-dimensional approximation for intersecting fractures avoids the more restrictive CFL condition corresponding to the equi-dimensional approximation with explicit time discretisation. Lower-dimensional fracture model results are compared with equi-dimensional model results. Fractures and barriers are efficiently modelled by lower-dimensional interfaces which yield comparable results to those of the equi-dimensional model. Pressure continuity is built into the model across highly conductive fractures, leading to reduced local degrees of freedom in the CVD-MPFA approximation. The formulation is applied to geologically complex fracture networks in three-dimensions. The effects of the fracture permeability, aperture and grid resolution are also assessed with respect to convergence and computational cost.</abstract><type>Journal Article</type><journal>Journal of Computational Physics</journal><volume>303</volume><journalNumber/><paginationStart>470</paginationStart><paginationEnd>497</paginationEnd><publisher/><placeOfPublication/><isbnPrint/><isbnElectronic/><issnPrint>0021-9991</issnPrint><issnElectronic/><keywords/><publishedDay>15</publishedDay><publishedMonth>12</publishedMonth><publishedYear>2015</publishedYear><publishedDate>2015-12-15</publishedDate><doi>10.1016/j.jcp.2015.10.001</doi><url/><notes/><college>COLLEGE NANME</college><department>Science and Engineering - Faculty</department><CollegeCode>COLLEGE CODE</CollegeCode><DepartmentCode>FGSEN</DepartmentCode><institution>Swansea University</institution><apcterm/><lastEdited>2021-01-06T12:17:40.0548488</lastEdited><Created>2015-11-23T12:57:20.7382778</Created><path><level id="1">Faculty of Science and Engineering</level><level id="2">School of Engineering and Applied Sciences - Uncategorised</level></path><authors><author><firstname>Raheel</firstname><surname>Ahmed</surname><order>1</order></author><author><firstname>Michael G.</firstname><surname>Edwards</surname><order>2</order></author><author><firstname>Sadok</firstname><surname>Lamine</surname><order>3</order></author><author><firstname>Bastiaan</firstname><surname>Huisman</surname><order>4</order></author><author><firstname>Mayur</firstname><surname>Pal</surname><order>5</order></author></authors><documents><document><filename>0024695-11022016172015.pdf</filename><originalFilename>AhmedThreeDimensionalControlVolume2015Postprint.pdf</originalFilename><uploaded>2016-02-11T17:20:15.0700000</uploaded><type>Output</type><contentLength>6881291</contentLength><contentType>application/pdf</contentType><version>Accepted Manuscript</version><cronfaStatus>true</cronfaStatus><embargoDate>2016-10-08T00:00:00.0000000</embargoDate><copyrightCorrect>true</copyrightCorrect></document></documents><OutputDurs/></rfc1807> |
spelling |
2021-01-06T12:17:40.0548488 v2 24695 2015-11-23 Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model 8903caf3d43fca03602a72ed31d17c59 Michael G. Edwards Michael G. Edwards true false 2015-11-23 FGSEN A novel cell-centred control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulation is presented for discrete fracture–matrix simulations on unstructured grids in three-dimensions (3D). The grid is aligned with fractures and barriers which are then modelled as lower-dimensional surface interfaces located between the matrix cells in the physical domain. The three-dimensional pressure equation is solved in the matrix domain coupled with a two-dimensional (2D) surface pressure equation solved over fracture networks via a novel surface CVD-MPFA formulation. The CVD-MPFA formulation naturally handles fractures with anisotropic permeabilities on unstructured grids. Matrix–fracture fluxes are expressed in terms of matrix and fracture pressures and define the transfer function, which is added to the lower-dimensional flow equation and couples the three-dimensional and surface systems. An additional transmission condition is used between matrix cells adjacent to low permeable fractures to couple the velocity and pressure jump across the fractures. Convergence and accuracy of the lower-dimensional fracture model is assessed for highly anisotropic fractures having a range of apertures and permeability tensors. A transport equation for tracer flow is coupled via the Darcy flux for single and intersecting fractures. The lower-dimensional approximation for intersecting fractures avoids the more restrictive CFL condition corresponding to the equi-dimensional approximation with explicit time discretisation. Lower-dimensional fracture model results are compared with equi-dimensional model results. Fractures and barriers are efficiently modelled by lower-dimensional interfaces which yield comparable results to those of the equi-dimensional model. Pressure continuity is built into the model across highly conductive fractures, leading to reduced local degrees of freedom in the CVD-MPFA approximation. The formulation is applied to geologically complex fracture networks in three-dimensions. The effects of the fracture permeability, aperture and grid resolution are also assessed with respect to convergence and computational cost. Journal Article Journal of Computational Physics 303 470 497 0021-9991 15 12 2015 2015-12-15 10.1016/j.jcp.2015.10.001 COLLEGE NANME Science and Engineering - Faculty COLLEGE CODE FGSEN Swansea University 2021-01-06T12:17:40.0548488 2015-11-23T12:57:20.7382778 Faculty of Science and Engineering School of Engineering and Applied Sciences - Uncategorised Raheel Ahmed 1 Michael G. Edwards 2 Sadok Lamine 3 Bastiaan Huisman 4 Mayur Pal 5 0024695-11022016172015.pdf AhmedThreeDimensionalControlVolume2015Postprint.pdf 2016-02-11T17:20:15.0700000 Output 6881291 application/pdf Accepted Manuscript true 2016-10-08T00:00:00.0000000 true |
title |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
spellingShingle |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model Michael G. Edwards |
title_short |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
title_full |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
title_fullStr |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
title_full_unstemmed |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
title_sort |
Three-dimensional control-volume distributed multi-point flux approximation coupled with a lower-dimensional surface fracture model |
author_id_str_mv |
8903caf3d43fca03602a72ed31d17c59 |
author_id_fullname_str_mv |
8903caf3d43fca03602a72ed31d17c59_***_Michael G. Edwards |
author |
Michael G. Edwards |
author2 |
Raheel Ahmed Michael G. Edwards Sadok Lamine Bastiaan Huisman Mayur Pal |
format |
Journal article |
container_title |
Journal of Computational Physics |
container_volume |
303 |
container_start_page |
470 |
publishDate |
2015 |
institution |
Swansea University |
issn |
0021-9991 |
doi_str_mv |
10.1016/j.jcp.2015.10.001 |
college_str |
Faculty of Science and Engineering |
hierarchytype |
|
hierarchy_top_id |
facultyofscienceandengineering |
hierarchy_top_title |
Faculty of Science and Engineering |
hierarchy_parent_id |
facultyofscienceandengineering |
hierarchy_parent_title |
Faculty of Science and Engineering |
department_str |
School of Engineering and Applied Sciences - Uncategorised{{{_:::_}}}Faculty of Science and Engineering{{{_:::_}}}School of Engineering and Applied Sciences - Uncategorised |
document_store_str |
1 |
active_str |
0 |
description |
A novel cell-centred control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulation is presented for discrete fracture–matrix simulations on unstructured grids in three-dimensions (3D). The grid is aligned with fractures and barriers which are then modelled as lower-dimensional surface interfaces located between the matrix cells in the physical domain. The three-dimensional pressure equation is solved in the matrix domain coupled with a two-dimensional (2D) surface pressure equation solved over fracture networks via a novel surface CVD-MPFA formulation. The CVD-MPFA formulation naturally handles fractures with anisotropic permeabilities on unstructured grids. Matrix–fracture fluxes are expressed in terms of matrix and fracture pressures and define the transfer function, which is added to the lower-dimensional flow equation and couples the three-dimensional and surface systems. An additional transmission condition is used between matrix cells adjacent to low permeable fractures to couple the velocity and pressure jump across the fractures. Convergence and accuracy of the lower-dimensional fracture model is assessed for highly anisotropic fractures having a range of apertures and permeability tensors. A transport equation for tracer flow is coupled via the Darcy flux for single and intersecting fractures. The lower-dimensional approximation for intersecting fractures avoids the more restrictive CFL condition corresponding to the equi-dimensional approximation with explicit time discretisation. Lower-dimensional fracture model results are compared with equi-dimensional model results. Fractures and barriers are efficiently modelled by lower-dimensional interfaces which yield comparable results to those of the equi-dimensional model. Pressure continuity is built into the model across highly conductive fractures, leading to reduced local degrees of freedom in the CVD-MPFA approximation. The formulation is applied to geologically complex fracture networks in three-dimensions. The effects of the fracture permeability, aperture and grid resolution are also assessed with respect to convergence and computational cost. |
published_date |
2015-12-15T03:29:21Z |
_version_ |
1763751138676441088 |
score |
11.037056 |