A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.

STUDY OF CRYSTAL GROWTH AND SOLUTE PRECIPITATION THROUGH FRONT TRACKING METHOD

Crystal growth and solute precipitation is a Stefan problem. It is a free boundary problem for a parabolic partial differential equation with a time-dependent phase interface. The velocity of the moving interface between solute and crystal is a local function. The dendritic structure of the crystal interface, which develops dynamically, requires high resolution of the interface geometry. These facts make the Lagrangian front tracking method well suited for the problem. In this paper, we introduce an upgraded version of the front tracking code and its associated algorithms for the numerical study of crystal formation. We compare our results with the smoothed particle hydrodynamics method （SPH） in terms of the crystal fractal dimension with its dependence on the Damkohler number and density ratio.

Multi-relaxation-time lattice Boltzmann front tracking method for two-phase flow with surface tension

In this paper, an improved incompressible multi-relaxation-time lattice Boltzmann-front tracking approach is proposed to simulate two-phase flow with a sharp interface, where the surface tension is implemented. The lattice Boltzmann method is used to simulate the incompressible flow with a stationary Eulerian grid, an additional moving Lagrangian grid is adopted to track explicitly the motion of the interface, and an indicator function is introduced to update the fluid properties accurately. The interface is represented by using a four-order Lagrange polynomial through fitting a set of discrete marker points, and then the surface tension is directly computed by using the normal vector and curvature of the interface. Two benchmark problems, including Laplace＇s law for a stationary bubble and the dispersion relation of the capillary wave between two fluids are conducted for validation. Excellent agreement is obtained between the numerical simulations and the theoretical results in the two cases.

CONSERVATIVE FRONT TRACKING IN HIGHER SPACE DIMENSIONS

September 12,2001 we propose a fully conservative front tracking algorithm in two space dimension. This algorithm first uses the point shifted algorithm on two adjacent time levels and then constructs space time hexahedra as computational units. We develop and prove a successful geometric construction under certain interface requirement. The algorithm has a first order local truncation error for cells near the tracked discontinuity, which is an improvement by one order of accuracy over most finite difference schemes, which have O (1) local truncation errors near discontinuities.

MARCHING CUBES BASED FRONT TRACKING METHOD AND ITS APPLICATION TO SOME INTERFACE INSTABILITY PROBLEMS

A front tracking method based on a marching cubes isosurface extractor, which is related filter generating isosurfaces from a structured point set, is provided to achieve sharp resolution for the simulation of non-diffusive interfacial flow. Compared with the traditional topology processing procedure, the current front tracking method is easier to be implemented and presents high performance in terms of computational resources. The numerical tests for 2-D highly-shearing flows and 3-D bubbles merging process are conducted to numerically examine the performance of the current methodology for tracking interfaces between two immiscible fluids The Rayleigh-Taylor （RT） and Richtmyer-Meshkov （RM） instability problems are successfully investigated with the present marching cubes based front tracking method.

Runge-Kutta Discontinuous Galerkin Method with Front Tracking Method for Solving the Compressible Two-Medium Flow on Unstructured Meshes

In this paper,we extend using the Runge-Kutta discontinuous Galerkin method together with the front tracking method to simulate the compressible twomedium flow on unstructured meshes.A Riemann problem is constructed in the normal direction in the material interfacial region,with the goal of obtaining a compact,robust and efficient procedure to track the explicit sharp interface precisely.Extensive numerical tests including the gas-gas and gas-liquid flows are provided to show the proposed methodologies possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow and the interfacial vicinities of the two-medium flow in many occasions.

A Front TrackingMethod for the Simulation of Compressible Multimedium Flows

A front trackingmethod combinedwith the real ghost fluidmethod(RGFM)is proposed for simulations of fluid interfaces in two-dimensional compressible flows.In this paper the Riemann problem is constructed along the normal direction of interface and the corresponding Riemann solutions are used to track fluid interfaces.The interface boundary conditions are defined by the RGFM,and the fluid interfaces are explicitly tracked by several connected marker points.The Riemann solutions are also used directly to update the flow states on both sides of the interface in the RGFM.In order to validate the accuracy and capacity of the new method,extensive numerical tests including the bubble advection,the Sod tube,the shock-bubble interaction,the Richtmyer-Meshkov instability and the gas-water interface,are simulated by using the Euler equations.The computational results are also compared with earlier computational studies and it shows good agreements including the compressible gas-water system with large density differences.

Numerical Method of Fabric Dynamics Using Front Tracking and Spring Model

We use front tracking data structures and functions to model the dynamic evolution of fabric surface.We represent the fabric surface by a triangulated mesh with preset equilibrium side length.The stretching and wrinkling of the surface are modeled by the mass-spring system.The external driving force is added to the fabric motion through the"Impulse method"which computes the velocity of the point mass by superposition of momentum.The mass-spring system is a nonlinear ODE system.Added by the numerical and computational analysis,we show that the spring system has an upper bound of the eigen frequency.We analyzed the system by considering two spring models and we proved in one case that all eigenvalues are imaginary and there exists an upper bound for the eigen-frequency.This upper bound plays an important role in determining the numerical stability and accuracy of the ODE system.Based on this analysis,we analyzed the numerical accuracy and stability of the nonlinear spring mass system for fabric surface and its tangential and normal motion.We used the fourth order Runge-Kutta method to solve the ODE system and showed that the time step is linearly dependent on the mesh size for the system.

AN EMBEDDED BOUNDARY METHOD FOR ELLIPTIC AND PARABOLIC PROBLEMS WITH INTERFACES AND APPLICATION TO MULTI-MATERIAL SYSTEMS WITH PHASE TRANSITIONS

The embedded boundary method for solving elliptic and parabolic problems in geometrically complex domains using Cartesian meshes by Johansen and Colella （1998, J. Comput. Phys. 147, 60） has been extended for elliptic and parabolic problems with interior boundaries or interfaces of discontinuities of material properties or solutions. Second order accuracy is achieved in space and time for both stationary and moving interface problems. The method is conservative for elliptic and parabolic problems with fixed interfaces. Based on this method, a front tracking algorithm for the Stefan problem has been developed. The accuracy of the method is measured through comparison with exact solution to a two-dimensional Stefan problem. The algorithm has been used for the study of melting and solidification problems.

Thermocapillary migration and interaction of two nondeformable drops

A numerical study on the interaction of two spherical drops in the thermocapillary migration is presented in the microgravity environment. Finite-difference methods are adopted. The interfaces of the drops are captured by the front-tracking technique. It is found that the arrangement of the drops directly influences their migration and interaction, and the motion of one drop is mainly determined by the disturbed temperature field because of the existence of the other drop.

ON THE CAUCHY PROBLEM FOR AW-RASCLE SYSTEM WITH LINEAR DAMPING

The existence of global BV solutions for the Aw-Rascle system with linear damping is considered.In order to get approximate solutions we consider the system in Lagrangian coordinates,then by using the wave front tracking method coupling with and suitable splitting algorithm and the ideas of[1]we get a sequence of approximate solutions.Finally we show the convergence of this approximate sequence to the weak entropic solution.

An Implementation ofMAC Grid-Based IIM-Stokes Solver for Incompressible Two-Phase Flows

Abstract.In this paper,a novel implementation of immersed interface method combined with Stokes solver on a MAC staggered grid for solving the steady two-fluid Stokes equations with interfaces.The velocity components along the interface are introduced as two augmented variables and the resulting augmented equation is then solved by the GMRES method.The augmented variables and/or the forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity,and are interpolated using cubic splines and are then applied to the fluid through the jump conditions.The Stokes equations are discretized on a staggered Cartesian grid via a second order finite difference method and solved by the conjugate gradient Uzawa-typemethod.The numerical results show that the overall scheme is second order accurate.The major advantages of the present IIM-Stokes solver are the efficiency and flexibility in terms of types of fluid flow and different boundary conditions.The proposed method avoids solution of the pressure Poisson equation,and comparisons are made to show the advantages of time savings by the present method.The generalized two-phase Stokes solver with correction terms has also been applied to incompressible two-phase Navier-Stokes flow.

Global Stability of Multi-wave Configurations for the Compressible Non-isentropic Euler System

This paper is contributed to the structural stability of multi-wave configurations to Cauchy problem for the compressible non-isentropic Euler system with adiabatic exponentγ∈(1,3].Given some small BV perturbations of the initial state,the author employs a modified wave front tracking method,constructs a new Glimm functional,and proves its monotone decreasing based on the possible local wave interaction estimates,then establishes the global stability of the multi-wave configurations,consisting of a strong 1-shock wave,a strong 2-contact discontinuity,and a strong 3-shock wave,without restrictions on their strengths.

One-Sided Exact Boundary Null Controllability of Entropy Solutions to a Class of Hyperbolic Systems of Conservation Laws with Characteristics with Constant Multiplicity

This paper proves the local exact one-sided boundary null controllability of entropy solutions to a class of hyperbolic systems of conservation laws with characteristics with constant multiplicity. This generalizes the results in [Li, T. and Yu, L., One-sided exact boundary null controllability of entropy solutions to a class of hyperbolic systems of conservation laws, To appear in Journal de Mathematiques Pures et Appliquees, 2016.] for a class of strictly hyperbolic systems of conservation laws.

Geometric Evolution Laws for Thin Crystalline Films: Modeling and Numerics

Geometrical evolution laws are widely used in continuum modeling of surface and interface motion in materials science.In this article,we first give a brief review of various kinds of geometrical evolution laws and their variational derivations,with an emphasis on strong anisotropy.We then survey some of the finite element based numerical methods for simulating the motion of interfaces focusing on the field of thin film growth.We discuss the finite element method applied to front-tracking,phase-field and level-set methods.We describe various applications of these geometrical evolution laws to materials science problems,and in particular,the growth and morphologies of thin crystalline films.

An Immersed Interface Method for the Simulation of Inextensible Interfaces in Viscous Fluids

In this paper,an immersed interface method is presented to simulate the dynamics of inextensible interfaces in an incompressible flow.The tension is introduced as an augmented variable to satisfy the constraint of interface inextensibility,and the resulting augmented system is solved by the GMRES method.In this work,the arclength of the interface is locally and globally conserved as the enclosed region undergoes deformation.The forces at the interface are calculated from the configuration of the interface and the computed augmented variable,and then applied to the fluid through the related jump conditions.The governing equations are discretized on a MAC grid via a second-order finite difference scheme which incorporates jump contributions and solved by the conjugate gradient Uzawa-type method.The proposed method is applied to several examples including the deformation of a liquid capsule with inextensible interfaces in a shear flow.Numerical results reveal that both the area enclosed by interface and arclength of interface are conserved well simultaneously.These provide further evidence on the capability of the present method to simulate incompressible flows involving inextensible interfaces.