An arbitrary Lagrangian-Eulerian method for simulating bubble growth in polymer foaming

We present a sharp-interface algorithm for simulating the diffusion-driven bubble growth in polymer foaming. A moving mesh of unstructured triangular elements tracks the expanding and deforming bubble surface. In the interior of the liquid, the mesh velocity is determined by solving a Laplace equation to ensure spatially smooth mesh movement. When mesh distortion becomes severe, remeshing and interpolation are performed. The governing equations are solved using a Galerkin finite-element formalism, with fully implicit time marching that requires iteration among the bubble and mesh deformation, gas diffusion and the flow and stress fields. Besides numerical stability, the implicit scheme also guarantees a smooth interfacial curvature as numerical disturbances on the interface are automatically relaxed through the iterations. The polymer melt is modeled as a viscoelastic Oldroyd-B fluid. First, we compute three benchmark problems to validate various aspects of the algorithm. Then we use a periodic hexagonal cell to simulate bubble growth in an isothermal two-dimensional foam, fed by a gaseous blowing agent initially dissolved in the melt to supersaturation. Results show two distinct stages: a rapid initial expansion followed by slow drainage of the liquid film between bubbles driven by capillarity. The effect of viscoelastic rheology is to enhance the speed of bubble growth in the first stage, and hinder film drainage in the second. Finally, we use axisymmetric simulations to investigate the thinning film between a bubble and a free surface. Melt viscoelasticity is shown to initially enhance film thinning but later resist it. An important insight from the simulations is that polymer strain-hardening, namely the steep increase of elongational viscosity with strain, helps stabilize the foam structure by suppressing bubble-bubble coalescence and bubble burst at the foam surface. This confirms prior observations in foam extrusion experiments.

[1]  R. Aseeva,et al.  Handbook of polymeric foams and foam technology. , 2004 .

[2]  Hongwei Zheng,et al.  A lattice Boltzmann model for multiphase flows with large density ratio , 2006, J. Comput. Phys..

[3]  Oliver G. Harlen,et al.  Bubble growth in a two-dimensional viscoelastic foam , 2006 .

[4]  L. G. Leal,et al.  Numerical solution of free-boundary problems in fluid mechanics. Part 2. Buoyancy-driven motion of a gas bubble through a quiescent liquid , 1984, Journal of Fluid Mechanics.

[5]  Mark Sussman,et al.  A computational study of the effect of initial bubble conditions on the motion of a gas bubble rising in viscous liquids , 2005 .

[6]  James J. Feng,et al.  Diffuse-interface simulations of drop coalescence and retraction in viscoelastic fluids , 2005 .

[7]  Wei Shyy,et al.  A fixed-grid, sharp-interface method for bubble dynamics and phase change , 2001 .

[8]  James J. Feng,et al.  Prediction of bubble growth and size distribution in polymer foaming based on a new heterogeneous nucleation model , 2003 .

[9]  P. George,et al.  Mesh Generation: Application to Finite Elements , 2007 .

[10]  Xavier Frank,et al.  Complex flow around a bubble rising in a non-Newtonian fluid. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  L. G. Leal,et al.  Numerical solution of free-boundary problems in fluid mechanics. Part 1. The finite-difference technique , 1984, Journal of Fluid Mechanics.

[12]  R. G. Fenton,et al.  Computer Simulation of Bubble-Growth Phenomena in Foaming , 2006 .

[13]  Eli Ruckenstein,et al.  Decay of standing foams: drainage, coalescence and collapse , 1997 .

[14]  I. Sevostianov,et al.  Elastic and electric properties of closed-cell aluminum foams Cross-property connection , 2006 .

[15]  R. Gendron,et al.  Rheological Behavior of Mixtures of Polystyrene with HCFC 142b and HFC 134a , 1999 .

[16]  Howard H. Hu Direct simulation of flows of solid-liquid mixtures , 1996 .

[17]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[18]  Bharat K. Soni,et al.  Mesh Generation , 2020, Handbook of Computational Geometry.

[19]  A. Caboussat A Numerical Method for the Simulation of Free Surface Flows with Surface Tension , 2006 .

[20]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[21]  Durian,et al.  Shear-Induced "Melting" of an Aqueous Foam. , 1999, Journal of colloid and interface science.

[22]  Barry Bernstein,et al.  Analysis of diffusion-induced bubble growth in viscoelastic liquids , 1998 .

[23]  Stéphane Popinet,et al.  Bubble collapse near a solid boundary: a numerical study of the influence of viscosity , 2002, Journal of Fluid Mechanics.

[24]  D. Reinelt,et al.  Structure of random monodisperse foam. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[25]  David P. Schmidt,et al.  DIRECT INTERFACE TRACKING OF DROPLET DEFORMATION , 2002 .

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

[27]  Natarajan S. Ramesh,et al.  A Non-Isothermal Model to Study the Influence of Blowing Agent Concentration on Polymer Viscosity and Gas Diffusivity in Thermoplastic Foam Extrusion , 1999 .

[28]  C. Han,et al.  Studies on structural foam processing. IV. Bubble growth during mold filling , 1981 .

[29]  Scott D. Phillips,et al.  Computational and experimental analysis of dynamics of drop formation , 1999 .

[30]  Hua-Xin Peng,et al.  Microstructure of ceramic foams , 2000 .

[31]  Robert C. Armstrong,et al.  Dynamics of polymeric liquids: Fluid mechanics , 1987 .

[32]  James J. Feng,et al.  Plasticization effects on bubble growth during polymer foaming , 2006 .

[33]  J. G. Hnat,et al.  Spherical cap bubbles and skirt formation , 1976 .

[34]  Jan Kurt Walter Sandler,et al.  Principal correlation of PVC melt elongational properties with foam cell morphology , 2000 .

[35]  Chunfeng Zhou,et al.  Phase-field simulations of interfacial dynamics in viscoelastic fluids using finite elements with adaptive meshing , 2006, J. Comput. Phys..

[36]  J. Rappaz,et al.  Numerical simulation of free surface incompressible liquid flows surrounded by compressible gas , 2005 .

[37]  M. Filipiak Mesh Generation , 2007 .

[38]  Ng Niels Deen,et al.  Numerical simulation of gas bubbles behaviour using a three-dimensional volume of fluid method , 2005 .

[39]  Malcolm R. Mackley,et al.  Experimental observations and modelling relating to foaming and bubble growth from pentane loaded polystyrene melts , 2004 .

[40]  Howard H. Hu,et al.  Direct numerical simulations of fluid-solid systems using the arbitrary Langrangian-Eulerian technique , 2001 .

[41]  Eric Saltel,et al.  Emc un logiciel d'edition de maillages et de contours bidimensionnels , 1990 .

[42]  M. Krasowska,et al.  Influence of surface active substances on bubble motion and collision with various interfaces. , 2005, Advances in colloid and interface science.

[43]  Peter Smereka,et al.  Axisymmetric free boundary problems , 1997, Journal of Fluid Mechanics.

[44]  W. Shyy,et al.  Computation of Solid-Liquid Phase Fronts in the Sharp Interface Limit on Fixed Grids , 1999 .

[45]  Thomas E. Beechem,et al.  Bubble growth mechanism in carbon foams , 2005 .

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

[47]  Oliver G. Harlen,et al.  Competition and interaction of polydisperse bubbles in polymer foams , 2006 .

[48]  Daniel Read,et al.  Bubble dynamics in viscoelastic fluids with application to reacting and non-reacting polymer foams , 2003 .