A Comparison of the Ghost Cell Technique for Front Tracking Method

The treatment of moving material interfaces and their vicinity is very important for compressible multifluids. In this paper, we propose one type of ghost fluid method based on Riemann solutions for front tracking method. The accuracy of the interface boundary condition is discussed for the gas-gas Riemann problem. It is shown that the solution of the ghost fluid method approximates the exact solution to second-order accuracy in the sense of comparing to the exact solution of a Riemann problem at the material interface. Numerical examples suggest that the present scheme is able to handle multifluids problems with large density differences and has the property of reduced conservation error.


Introduction
The dynamics of interfaces separating different fluids in compressible flows is of interest in several scientific fields as diverse as astrophysics and geophysics.It is also of significant importance in many engineering applications.A relatively dominant difficulty for simulating compressible multi-medium flow is the treatment of moving material interfaces and their vicinity.In general, there are two basic issues that need to be taken into account.One is to capture the interface location and topological changes accurately.Researchers can take all kinds of effective measures such as volume of fluid method [1] or level set technique [2] or front tracking technique [3] to deal with it.The other is to faithfully simulate the interface state including physical nonlinear wave interaction occurring at the interface.
The ghost fluid method (GFM) based techniques [4][5][6][7][8] provide us simple and flexible ways for handling multimedium flows with immiscible material interfaces.A closely related ghost cell level set method [4] was proposed by Fedkiw.But the GFM is problem-related and not suitable for some cases like high speed jet impacting [5].To take into consideration the influence of both, wave interaction and material properties on the interfacial evolution led to the development of a modified ghost fluid method (MGFM) [5,7].Wang [8] also presented a real ghost fluid method (RGFM) by predicting the flow states for the real fluid nodes just next to the interface and the ghost fluid nodes using the Riemann problem solver.
In this paper, a relatively class of GFM based on Riemann solutions is developed.One solves a system of a two-shock approximation to the Riemann problem for prediction of the condition at the interface.Then one uses the predicted pressure and velocity as those for the ghost fluid and the real fluid nodes just next to the interface.The isobaric entropy is employed at the interface to fix the real fluid density at node just next to the interface to suppress the possible "overheating" and also to assign density for the ghost fluid.The RGFM is one of this type of GFM.We combine the GFM with a front tracking method (GFM-FT).The front tracking method [12] is used here.The accuracy of the interface boundary condition is discussed for the gas-gas Riemann problem.
The conservation errors of this class of GFM and MGFM are analyzed for FT method.Numerical tests show that the GFM-FT has the property of reduced conservation error comparing to MGFM-FT.

Equations
The Euler equations governing one-dimensional compressible flows are written as where Here  is the density, is the velocity, is the pressure, and is the total energy per unit volume.The total energy is the sum of internal energy and kinetic energy, where is the internal energy per unit mass.
e For closure of Equation ( 1), the equation of state (EOS) is required.The EOS for compressible gases and water can be expressed in the following consistent form as where  and are treated as constants.For the air, B 1.4 . For the water medium (Tait's equation),

Tracking Fluid Interfaces
Suppose that the interface lies between node i and node i + 1.To define the Riemann problem at the interface, the two initial constant states of the Riemann problem are simply given as . Once the Riemann problem is defined, the two shock Riemann problem solver [5]

Ghost Fluid
The initial conditions for the governing Equation (1) are given as Regardless of the numerical scheme used, an algorithm based on the GFM essentially consists of solving two separate Riemann problems in the two respective single media with one-sided ghost fluid.One is in medium 1 (on the left side of the interface) with the initial conditions of and it solves from the grid point 1 on the left-end to the ghost point, The other is in medium 2 (on the right side of the interface) with the initial conditions of and it solves from the ghost point to the end point on the right.Here, "*" indicates the fluid states have been replaced by the ghost fluid states.
With the method of characteristics, we have the relations where ,    are the densities and sound speeds of the ghost fluid.
The ghost fluid states can be chosen by solving Equation (7) and (8).The simplest case is that the ghost fluid pressure be defined as the interface pressure, i.e.
Hence the integral in Equation ( 7) and ( 8) become zero and the ghost fluid velocity is Furthermore, one can find that any ghost fluid density can satisfy Equation ( 5) or (6).
For Equation ( 5) or ( 6), Hu [6] had considered two algorithms.The focus is on defining the ghost fluid states while the pressure and velocity in the real fluid side are taken for granted, except for the correction made to the density at the real fluid nodes next to the interface to overcome the overheating.Indeed, the real fluid states next to the interface instead of the ghost fluid states should be predicted by solving a Riemann problem [8].This will result in the complete redefinition of the real fluid next to the interface.In doing so, one will find that the nonphysical reflection for the shock impedance matching problem can be greatly and further suppressed.
For medium 1, we take I p and I u as the pressure and velocity at nodes , , For medium 2, we also take 2, i  I p and I u as the pressure and velocity at nodes 1 i  , , , 1 i i   .The ghost fluid density is based on isentropic fixing.Here, we shall consider two cases: Case 1: For medium 1, the predicted isobaric entropy at the interface is employed to fix the real fluid density at point i and the ghost fluid density at point , .A similar procedure is used for computation in medium 2. This is the RGFM [8] which can be seen from Figure 1.
For medium 1, one employs isobaric entropy 1 1 1 to fix the real and ghost fluid density.
The rest of the procedures are the same as those for Case 1.We take this method as SGFM which can be seen from Figure 2.

The Accuracy of SGFM and RGFM
We only consider the gas-gas Riemann problem.The accuracy of MGFM had been discussed [13].It may be noted that the accuracy discussed should be interpreted as how accurate the boundary conditions are implicitly imposed at the material interface and how accurate the interface states are approximated by the GFM technique; it should not be explained or analyzed as the accuracy of final numerical solution and/or the numerical scheme used because the detailed numerical scheme is never involved in the discussion.An approximate Riemann problem solver (ARPS) based on a doubled shock structure is used for the Riemann problem.For the SGFM  applied to wing error the Riemann problem (4), the follo estimations are held [13] According to SGFM, one has Using Equation (11) For the GFM Riemann problem (6), in a similar way, on onservation Errors rvation for each me-e can show that error estimates (C) and (D) in Theorem 1 are true.
The states of the ghost cell for RGFM are in accord with those for SGFM except the density.So Theorem 1 is also true for RGFM.It is shown that the solution of the SGFM and RGFM approximate the exact solution to at least second-order accuracy in the sense of comparing to the exact solution of a Riemann problem at the material interface.

The C
Both overall conservation and conse dium will be evaluated using Equation (1) over the com- x x .We denote RHSL(n), RHSR(n) and RHST(n conservation errors for the medium on the left, right side of the interface and over the whole computational domain, respectively, as follows: ) as the he overall average conservation error can be taken as


We shall use two cases to study the cons of the RGFM and SGFM.ervation errors Case 1: a strong shock impacting on an air-air interface (i.e., Problem 1 in Section 7) From Table 1, it is clear that very large errors incur for the MUSCL-based MGFM-FT in the first 60 steps for this specific case, while these are suppressed very well by the MUSCL-based SGFM-FT and RGFM-FT.The conservative errors of the RGFM-FT and SGFM-FT are tend to zero faster than the MGFM-FT.
Case 2: This is a problem of a water column separated by the air and oscillating in the middle of a closed tube [8].The problem is sketched in Figure 3.The initial status of the air and water are defined as


With the fluids moving to the right, the air on the left of the water column expands while the air on the right of the water column is compressed.The pressure decreases at the former and increases at the latter; this leads to the deceleration and eventually stagnation of the fluids in the tube.Then the fluids move in the reverse direction due to the pressure gradient, and so on with the water column oscillating.Pressure waves are evident in the tube due to the impact of the air on the end boundaries and interacttion with the water column.The conservation errors of the air mass M(t) are given as where   a m t is the total mass of the air in the tube at time t.It can indeed be observed that the conservation errors associated with the air mass appear to be oscillating similar to that observed in [8]. Figure 4(a) depicts that MGFM has an oscillatory M(t) with a much larger trend of net conversion of the air mass into water when

GFM-FT, SGFM-F
Figure 5 shows the time evolution of the pr essure coef-

of t ectio
ficients at the wall boundaries calculated using MGFM, RGFM and SGFM. Figure 5(c) shows that there is a increase of the pressure amplitude with time for MGFM, while the pressure amplitude with RGFM and SGFM is changeless on the whole as show in Figures 5(a

Applications
For all the one-dimensional problems the computational [0,1] with 201 nodes uniformly dis-


This problem is solved at time .The shock limiter is used [11].Figures 6(a)-(c) show that the densi d by t For this case, a very strong shock is physically reflected back into the air.Figures 7(a)-(c) show that the den  This test has solution consisting of a rarefaction wave, two contact discontinuities and three shock waves.The domain is taken as tributed unless stated otherwise.
Problem 1: Gaseous shock in impedance-matching medium [8].In this case a strong shock on the left side of the interface impacts on a gas-gas interface; The shock strength is 100 and the initial position of the shock is same as the interface, which is located at x 0 = 0.2.The states on the left and right sides of the interface are defined respectively as MUSCL-based SGFM-FT method is used.Figures 8(a)- (c) describe the surprisingly good solution profiles at time 0.00028.

Conclusion
In this paper, we developed a kind of GFM based on Riemann solutions.The R M is one of this type of GFM.A simplified front tracking algorithm had been used to track contact discontinuity.It is noted that different isentropic entropy can produce little effect on the conservation error and this kind of GFM has the property of reduced conservation error.A rigorous analysis is car- ried out on the accuracy of the GFM when applied to the gas-gas Riemann problem.It is shown that at the material interface this type of GFM solution approximates the exact solution to at least second-order accuracy in the sense of comparing to the exact solution of a Riem n an problem.When the shock tor is used, discernible non-physical hump and trough which can not be avoided in shock impedance matching (-like) problems [7,8] do not appear for our proposed method.The ghost cell nique can be used for simulating compressible multimedium flow.
the pressure and veloc node the interface are also redefined.Following conclusions are held for the SGFM when applied to the Riemann problem (4).Theorem 1.The foll ity at the real fluid next to owing error estimates are held for th (A) e respective GFM Riemann problems (5) and (6) using the SGFM are the exact interfacial rem 2.1[13] to the GFM Riem velocity and pressu GFM Riemann problem (5) and (6) using the SGFM.Proof: Applying Theo ann problem (5), one gets 18) Here, I x denotes the interface position; A F and B F are es at A flux x and B x , respectively.IL F and IR F are fluxes at the respective left and right side f the rface.

Figure 3 .
Figure 3. Initial status of the oscillating water column problem.

Figure 4 .
Figure 4. Comparison of c rvation errors amo the MGFM, the RGFM and the SGFM for Case 2 in Section 6. onse ng (a) the MGFM-FT and the SGFM-FT; (b) the RGFM-FT; (c) the SGFM-FT.

Figure 5 .
Figure 5.Time evolution he pressure coefficients at boundaries for Case 2 in S n 6. Solid line: at the left boundary; dashed line: at the right boundary.

Figure 8 .
Figure 8. Problem 3: The "o" denotes the numerical results.The solid line denotes the exact solution.