Detailed numerical simulations of pore competition in idealized micro-spall using the VOF method

被引:6
作者
Malan, L. C. [1 ,2 ,3 ]
Ling, Y. [1 ,2 ,4 ]
Scardovelli, R. [5 ]
Llor, A. [6 ]
Zaleski, S. [1 ,2 ]
机构
[1] UPMC Univ Paris 06, Sorbonne Univ, UMR 7190, Inst Jean Rond Alembert, F-75005 Paris, France
[2] CNRS, Inst Jean Rond Alembert, UMR 7190, F-75005 Paris, France
[3] Univ Cape Town, Dept Mech Engn, InCFD, ZA-7701 Cape Town, South Africa
[4] Baylor Univ, Dept Mech Engn, Waco, TX 76798 USA
[5] Univ Bologna, DIN Lab Montecuccolino, I-40136 Bologna, Italy
[6] CEA, DAM, DIF, F-91297 Arpqjon, France
关键词
Micro-spall; Pore competition; Volume-of-Fluid; Bubbles; Cavitation; Free-surface; DYNAMIC FRAGMENTATION; MELTED METALS; LEVEL SET; VOLUME; FLUID; INTERFACES; BOUNDARY; COLLAPSE; FLOW;
D O I
10.1016/j.compfluid.2019.05.011
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Under shock loading, metals have been found to melt and with reflection of the shock wave from the material free surface, cavities nucleate and grow. This process is referred to as microspalling and has been the topic of several experimental investigations. Measurements during the cavity growth phase are not possible at present and we present here a Detailed Numerical Simulation of an idealized problem where we assume an inviscid, incompressible liquid subject to a constant expansion rate with cavities at a vanishing vapour pressure. To allow for a time-varying gas volume a free-surface interface condition has been implemented in an existing incompressible multiphase Navier-Stokes solver, PARIS, using a Volume-Of-Fluid method. The gas flow remains unsolved and is instead assumed to have a fixed pressure which is applied to the liquid through a Dirichlet boundary condition on the liquid-gas interface. Gas bubbles are tracked individually, allowing the gas pressure to be prescribed using a suitable equation of state. Simulations with hundreds of bubbles have been performed in a fixed domain under a constant rate of expansion. A bubble competition is observed: larger bubbles tend to expand more rapidly at the demise of smaller ones. The time scale of this competition is shown to depend on the Weber number. (C) 2019 Elsevier Ltd. All rights reserved.
引用
收藏
页码:60 / 72
页数:13
相关论文
共 42 条
[1]  
Andriot P, 1983, INFLUENCE MELTING SH, P277
[2]  
[Anonymous], 2013, SC 13
[3]   Interface reconstruction with least-squares fit and split advection in three-dimensional Cartesian geometry [J].
Aulisa, E. ;
Manservisi, S. ;
Scardovelli, R. ;
Zaleski, S. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2007, 225 (02) :2301-2319
[4]   VOFI - A library to initialize the volume fraction scalar field [J].
Bna, S. ;
Manservisi, S. ;
Scardovelli, R. ;
Yecko, P. ;
Zaleski, S. .
COMPUTER PHYSICS COMMUNICATIONS, 2016, 200 :291-299
[5]   Numerical integration of implicit functions for the initialization of the VOF function [J].
Bna, S. ;
Manservisi, S. ;
Scardovelli, R. ;
Yecko, P. ;
Zaleski, S. .
COMPUTERS & FLUIDS, 2015, 113 :42-52
[6]   On the properties and limitations of the height function method in two-dimensional Cartesian geometry [J].
Bornia, G. ;
Cervone, A. ;
Manservisi, S. ;
Scardovelli, R. ;
Zaleski, S. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 230 (04) :851-862
[7]   A level set method for vapor bubble dynamics [J].
Can, Edip ;
Prosperetti, Andrea .
JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 231 (04) :1533-1552
[8]  
Chan R. K.-C., 1970, Journal of Computational Physics, V6, P68, DOI 10.1016/0021-9991(70)90005-7
[9]  
Childs H., 2012, HIGH PERFORMANCE VIS, P357
[10]   ON CONVERGENCE OF DISCRETE APPROXIMATIONS TO NAVIER-STOKES EQUATIONS [J].
CHORIN, AJ .
MATHEMATICS OF COMPUTATION, 1969, 23 (106) :341-&