On the elliptic mesh generation in domains containing multiple inclusions and undergoing large deformations

We present an improved method to generate a sequence of structured meshes even when the physical domain contains deforming inclusions. This method belongs to the class of Arbitrary Lagrangian-Eulerian (ALE) methods for solving moving boundary problems. Its tools are either (a) separate mappings of the domain boundaries and enforcing the node distribution on lines emanating from singular points or (b) domain decomposition and separate mappings of each subdomain using suitable coordinate systems. The latter is shown to be more versatile and general. In both cases a set of elliptic equations is used to generate the grid extending in this way the method advanced by Dimakopoulos and Tsamopoulos [Y. Dimakopoulos, J.A. Tsamopoulos, A quasi-elliptic transformation for moving boundary problems with large anisotropic deformations, J. Comput. Phys. 192 (2003) 494-522]. We shall present examples where this earlier method and all other mesh generating methods which are based on a conformal mapping or solving a quasi-elliptic set of PDEs fail to produce an acceptable mesh and accurate solutions in such geometries. Furthermore, in contrast to other methods, appropriate boundary conditions and constraints such as, orthogonality of specific mesh lines and prespecified node distributions on them, can be easily implemented along a specific part of the domain or its boundary. Hence, no attractive terms at specific corners or singular points are needed. To increase the mesh resolution around the moving interfaces while keeping low the memory requirements and the computational time, a local mesh refinement technique has been incorporated as well. The method is demonstrated in two challenging examples where no remeshing is required in spite of the large domain deformations. In the first one, the transient growth of two bubbles embedded in a viscoelastic filament undergoing stretching in the axial direction is examined, while in the second one the linear and non-linear dynamics of two bubbles in a viscous medium are determined in an acoustic field. The large elasticity of the filament in the first case or the large inertia in the second case coupled with the externally induced large deformations of the liquid domain requires the accurate calculation which is achieved by the method we propose herein. The governing equations are solved using the finite element/Galerkin method with appropriate modifications to solve the hyperbolic constitutive equation of a viscoelastic fluid. These are coupled with an implicit Euler method for time integration or with Arnoldi's algorithm for normal mode analysis.

[1]  A. Dvinsky.  Adaptive grid generation from harmonic maps on Reimannian manifolds , 1991 .

[2]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[3]  M. Fortin,et al.  A new approach for the FEM simulation of viscoelastic flows , 1989 .

[4]  Timothy J. Baker,et al.  Mesh generation: Art or science? , 2005 .

[5]  M. S Carvalho,et al.  Regular ArticleThree-Dimensional Stability Analysis of Free Surface Flows: Application to Forward Deformable Roll Coating☆ , 1999 .

[6]  J. Tsamopoulos,et al.  On the gas-penetration in straight tubes completely filled with a viscoelastic fluid , 2004 .

[7]  C. Macosko,et al.  Bubble growth and collapse in viscoelastic liquids analyzed , 1984 .

[8]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[9]  V. Mavrantzas,et al.  Numerical simulation of bubble growth in Newtonian and viscoelastic filaments undergoing stretching , 2004 .

[10]  Transient displacement of a Newtonian fluid by air in straight or suddenly constricted tubes , 2003 .

[11]  Robert C. Armstrong,et al.  On the numerical stability of mixed finite-element methods for viscoelastic flows governed by differential constitutive equations , 1993 .

[12]  R. Codina.  Comparison of some finite element methods for solving the diffusion-convection-reaction equation , 1998 .

[13]  Bernardo Cockburn,et al.  Devising discontinuous Galerkin methods for non-linear hyperbolic conversation laws , 2001 .

[14]  John Tsamopoulos,et al.  Bjerknes forces between two bubbles. Part 2. Response to an oscillatory pressure field , 1993, Journal of Fluid Mechanics.

[15]  J. Tsamopoulos,et al.  Transient displacement of a viscoplastic material by air in straight and suddenly constricted tubes , 2003 .

[16]  Y. Kallinderis,et al.  Numerical study of incompressible flow about fixed cylinder pairs , 2003 .

[17]  Roland Keunings,et al.  Transient deformation of an inviscid inclusion in a viscoelastic extensional flow , 1988 .

[18]  J. Tsamopoulos,et al.  Steady bubble rise and deformation in Newtonian and viscoplastic fluids and conditions for bubble entrapment , 2008, Journal of Fluid Mechanics.

[19]  Yannis Kallinderis,et al.  Finite element method for incompressible viscous flows with adaptive hybrid grids , 1994 .

[20]  J. Tsamopoulos,et al.  A quasi-elliptic transformation for moving boundary problems with large anisotropic deformations , 2003 .

[21]  V. Mavrantzas,et al.  Numerical simulation of multiple bubbles growing in a Newtonian liquid filament undergoing stretching , 2006 .

[22]  Joe F. Thompson,et al.  Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies , 1974 .

[23]  J. Tsamopoulos,et al.  Transient squeeze flow of viscoplastic materials , 2006 .

[24]  V. Bjerknes.  Fields of Force , 2010 .

[25]  O. Schenk,et al.  ON FAST FACTORIZATION PIVOTING METHODS FOR SPARSE SYMMETRI C INDEFINITE SYSTEMS , 2006 .

[26]  Otilio Rojas,et al.  Generation of curvilinear coordinates on multiply connected regions with boundary singularities , 2007, J. Comput. Phys..

[27]  L. E. Scriven,et al.  The oscillations of a fluid droplet immersed in another fluid , 1968, Journal of Fluid Mechanics.

[28]  Nikolaos A. Malamataris,et al.  A new outflow boundary condition , 1992 .

[29]  J. Tsamopoulos,et al.  Steady extrusion of viscoelastic materials from an annular die , 2008 .

[30]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[31]  R. Sani,et al.  On the time-dependent solution of the incompressible Navier-Stokes equations in two and three dimensions , 1980 .

[32]  R. Bird.  Dynamics of Polymeric Liquids , 1977 .

[33]  R. Cairncross,et al.  Effect of substrate flexibility on dynamic wetting: a finite element model , 2003 .

[34]  J. F. Thompson,et al.  A general three-dimensional elliptic grid generation system on a composite block structure , 1987 .

[35]  Emil-Alexandru Brujan.  A first-order model for bubble dynamics in a compressible viscoelastic liquid , 1999 .

[36]  J. Tsamopoulos,et al.  Transient displacement of Newtonian and viscoplastic liquids by air in complex tubes , 2007 .

[37]  Yannis Dimakopoulos,et al.  A direct comparison between volume and surface tracking methods with a boundary-fitted coordinate transformation and third-order upwinding , 2007, J. Comput. Phys..

[38]  U. Parlitz,et al.  Spatio–temporal dynamics of acoustic cavitation bubble clouds , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[39]  L. E. Scriven,et al.  Discretization of free surface flows and other moving boundary problems , 1992 .

[40]  Jeremiah Brackbill,et al.  Applications and generalizations of variational methods for generating adaptive meshes , 1982 .

[41]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2004, Future Gener. Comput. Syst..

[42]  R. Tanner,et al.  A new constitutive equation derived from network theory , 1977 .

[43]  D. Birchall,et al.  Computational Fluid Dynamics , 2020, Radial Flow Turbocompressors.

[44]  L. E. Scriven,et al.  Three-Dimensional Stability Analysis of Free Surface Flows , 1999 .

[45]  Yannis Kallinderis,et al.  A priori mesh quality estimation via direct relation between truncation error and mesh distortion , 2009, J. Comput. Phys..

[46]  Joe F. Thompson,et al.  Boundary-fitted coordinate systems for numerical solution of partial differential equations—A review , 1982 .

[47]  Robert E. Smith,et al.  Algebraic grid generation , 1982 .

[48]  Chongyoup Kim.  Collapse of spherical bubbles in Maxwell fluids , 1994 .

[49]  Joe F. Thompson,et al.  Numerical grid generation , 1985 .

[50]  I. Babuska,et al.  Finite Element Analysis , 2021 .

[51]  Rekha Ranjana Rao,et al.  A Newton-Raphson Pseudo-Solid Domain Mapping Technique for Free and Moving Boundary Problems , 1996 .

[52]  D. Peric,et al.  On finite element modelling of surface tension , 2006 .

[53]  Bharat K. Soni,et al.  Handbook of Grid Generation , 1998 .

[54]  D. Perić,et al.  On finite element modelling of surface tension: Variational formulation and applications – Part II: Dynamic problems , 2006 .

[55]  Robert A. Brown,et al.  BOUNDARY-CONFORMING MAPPING APPLIED TO COMPUTATIONS OF HIGHLY DEFORMED SOLIDIFICATION INTERFACES , 1992 .

[56]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[57]  Robert A. Brown,et al.  Solution of free-boundary problems using finite-element/Newton methods and locally refined grids: Application to analysis of solidification microstructure , 1993 .

[58]  A. Prosperetti,et al.  Bubble Dynamics and Cavitation , 1977 .

[59]  Andrew P. Kuprat,et al.  Boundary Orthogonality in Elliptic Grid Generation , 1998 .

[60]  K. Morgan,et al.  Recent advances in numerical methods in fluids , 1980 .

[61]  John Tsamopoulos,et al.  Bjerknes forces between two bubbles. Part 1. Response to a step change in pressure , 1993, Journal of Fluid Mechanics.

[62]  N. Phan-Thien.  A Nonlinear Network Viscoelastic Model , 1978 .