Abstract

A general, higher-order, conservative and bounded interpolation for the dynamic and adaptive meshing of control-volume fields dual to continuous and discontinuous finite element representations is presented. Existing techniques such as node-wise interpolation are not conservative and do not readily generalise to discontinuous fields, whilst conservative methods such as Grandy interpolation are often too diffusive. The new method uses control-volume Galerkin projection to interpolate between control-volume fields. Bounded solutions are ensured by using a post-interpolation diffusive correction. Example applications of the method to interface capturing during advection and also to the modelling of multiphase porous media flow are presented to demonstrate the generality and robustness of the approach.

Highlights

  • Dynamic, adaptive meshing is often used during numerical simulations of transient fluid flows to improve accuracy [1,2]

  • This is achieved by first mapping the control-volume field into a finite element representation on the donor mesh by Galerkin projection, interpolating onto the target mesh by constructing a finite element supermesh from the intersection of the donor and target meshes and projecting back to a control-volume representation on the target mesh

  • A series of test cases have demonstrated that the method is superior to both node-wise and Grandy interpolation for control-volume fields with a discontinuous dual finite element representation

Read more

Summary

Introduction

Adaptive meshing is often used during numerical simulations of transient fluid flows to improve accuracy [1,2]. The mesh may change every time-step depending on the error metrics used to control refinement and coarsening This inevitably means that data must be mapped from one mesh to another using an interpolation algorithm [3]. Interfaces between fluids as well as material property boundaries in both inertial and porous media flow [4,5,6,7,8] In the latter example, discontinuous representations of control-volume fields (that is to say when the dual finite element representation is discontinuous) may further help to capture abrupt changes in material properties (e.g. rock permeability) [9,10].

Node-wise interpolation
Grandy interpolation
Control-volume Galerkin interpolation
Three projections
Step two – projection from old to new mesh by building a supermesh
Implementation of bounded solution corrections
Interpolation during mesh refinement
Higher-order convergence
Advection of a square scalar field
Advection of a circular scalar field
Viscous instabilities in two phase immiscible displacement in a porous medium
Findings
Conclusions
Full Text
Paper version not known

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.