Abstract

In this paper, we describe an interface reconstruction method in two dimension. This method is an extension of DPIR [1], which reconstructs continuous interfaces and preserves partial volumes using dynamic programming. First we extend the method to curved interfaces. Then, we present tools to improve its robustness in order to apply it to unstructured grid. Finally, we describe an extension to three materials.

Highlights

  • Interface reconstruction (IR) methods are encountered in numerical simulation of multi-material or multifluid flows

  • This method consists in assuming that the interface in each mixed cell is made of a segment joining two of its edges. The normal to this segment is colinear to the gradient of the volume fraction ∇C and its position is obtained by assuming an exact conservation of partial volumes

  • If the volume fraction between Mi and Mi+1 is too small with respect to the correction, the second penalization p becomes big and those points are not chosen by the minimization

Read more

Summary

Introduction

Interface reconstruction (IR) methods are encountered in numerical simulation of multi-material or multifluid flows. This method consists in assuming that the interface in each mixed cell (that is such 0 < C < 1) is made of a segment joining two of its edges The normal to this segment is colinear to the gradient of the volume fraction ∇C and its position is obtained by assuming an exact conservation of partial volumes. In [1], a new reconstruction method which ensures continuity of the interface and preserves volume fractions has been introduced This new interface reconstruction method, called DPIR (Dynamic Programming Interface Reconstruction), is introduced and will be used as a starting point for the presented work. The addition of a control point in each cell in order to recover the correct volume fractions. The method is tested on two test cases with triple point configurations on cartesian meshes, giving encouraging results for future applications to unstructured meshes

Interface reconstruction with DPIR
First step: minimization of J with dynamic programming
Second step: local correction of volume fractions
DPIR extension to curved interfaces reconstruction
Robustness improvement of DPIR
A new search direction for the control point
A new penalty term
Numerical results
DPIR extension to triple point reconstruction
Classification of triple cells
The new algorithm for interface reconstruction
First step: one-against-all approach
Second step: mixed cells correction
Complexity
A example of DPIR reconstruction on a triple point configuration
Perspectives on the filament issue
Conclusion and perspectives
Full Text
Published version (Free)

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