PArallel, Robust, Interface Simulator (PARIS)

Paris (PArallel, Robust, Interface Simulator) is a finite volume code for simulations of immiscible multifluid or multiphase flows. It is based on the ?one-fluid?formulation of the Navier?Stokes equations where different fluids are treated as one material with variable properties, and surface tension is added as a singular interface force. The fluid equations are solved on a regular structured staggered grid using an explicit projection method with a first-order or second-order time integration scheme. The interface separating the different fluids is tracked by a Front-Tracking (FT) method, where the interface is represented by connected marker points, or by a Volume-of-Fluid (VOF) method, where the marker function is advected directly on the fixed grid. Paris is written in Fortran95/2002 and parallelized using MPI and domain decomposition. It is based on several earlier FT or VOF codes such as Ftc3D, Surfer or Gerris. These codes and similar ones, as well as Paris, have been used to simulate a wide range of multifluid and multiphase flows.Program summaryProgram Title: PArallel Robust Interface Simulator ? ParisCPC Library link to program files: https://doi.org/10.17632/5cb2yrfx7r.1Licensing provisions: GPLv3. Programming language: Fortran95/2002. Parallelized using MPI and domain decomposition. Nature of problem: Paris is a free code, or software, for computational fluid dynamics (CFD) of multiphase flows (or computational multiphase fluid dynamics (CMFD)), specialized in the numerical simulation of interfacial fluid flows, involving droplets, bubbles and waves, as described in the book by Tryggvason, Scardovelli and Zaleski [1]. It solves the Euler or Navier?Stokes equations in the one-fluid formulation of two-phase flow, including a surface tension force. It computes complex flows such as fast atomizing jets or droplets, expanding cavitation bubble clusters and multiphase flow through porous media.Solution method: The code mostly implements the methods described in the book by Tryggvason, Scardovelli and Zaleski [1]. Time stepping is performed using a first-order or a second-order in time predictor?corrector method with an explicit projection step for the pressure. Spatial discretization is by finite volumes on a regular cuboid grid. Interface tracking is performed with the Front-Tracking (FT) method or the Volume-of-Fluid (VOF) method. In the VOF version Paris uses either the LagrangianExplicit (LE) advection method or the exactly mass-conserving method of Weymouth and Yue [2]. The normal computation is performed using the Mixed-Youngs-Centered (MYC) scheme. A mass? momentum advection method has been also implemented that is consistent with the VOF advection [3]. Curvature is computed with the Height Function (HF) method. This is combined with the balanced Continuous Surface Force (CSF) method to compute surface tension forces.