Commun. Comput. Phys., 8 (2010), pp. 797-822.


Central Schemes and Second Order Boundary Conditions for 1D Interface and Piston Problems in Lagrangian Coordinates

Riccardo Fazio 1*, Giovanni Russo 2

1 Department of Mathematics, University of Messina, Salita Sperone 31, 98166 Messina, Italy.
2 Department of Mathematics and Computer Science, University of Catania, Viale Andrea Doria 6, 95125 Catania, Italy.

Received 31 January 2009; Accepted (in revised version) 22 January 2010
Available online 17 May 2010
doi:10.4208/cicp.310109.220110a

Abstract

We study high-resolution central schemes in Lagrangian coordinates for the one-dimensional system of conservation laws describing the evolution of two gases in slab geometry separated by an interface. By using Lagrangian coordinates, the interface is transformed to a fixed coordinate in the computational domain and, as a consequence, the movement of the interface is obtained as a byproduct of the numerical solution. The main contribution is the derivation of a special equation of state to be imposed at the interface in order to avoid non-physical oscillations. Suitable boundary conditions at the piston that guarantee second order convergence are described. We compare the solution of the piston problem to other results available in the literature and to a reference solution obtained within the adiabatic approximation. A shock-interface interaction problem is also treated. The results on these tests are in good agreement with those obtained by other methods.

AMS subject classifications: 65M06, 65M99, 76T05

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Lagrangian central schemes, Euler equations, interface conditions, boundary conditions.

*Corresponding author.
Email: rfazio@dipmat.unime.it (R. Fazio), russo@dmi.unict.it (G. Russo)
 

The Global Science Journal