In this paper we develop a conservative sharp-interface method dedicated to simulating multiple compressible fluids. Numerical treatments for a cut cell shared by more than two materials are proposed. First, we simplify the interface interaction inside such a cell with a reduced model to avoid explicit interface reconstruction and complex flux calculation. Second, conservation is strictly preserved by an efficient conservation correction procedure for the cut cell. To improve the robustness, a multi-material scale separation model is developed to consistently remove non-resolved interface scales. In addition, the multi-resolution method and local time-stepping scheme are incorporated into the proposed multi-material method to speed up the high-resolution simulations.Various numerical test cases, including the multi-material shock tube problem, inertial confinement fusion implosion, triple-point shock interaction and shock interaction with multi-material bubbles, show that the method is suitable for a wide range of complex compressible multi-material flows.
The paper presents phenomenology of interaction and penetration of liquid-liquid material interfaces initiated by shock-driven collapse of single and multiple microbubbles situated near the material interface. Previous experimental studies have established such a generic setting as relevant for the investigation of sonoporation, i.e., the perforation of live cells by microbubble collapses. We consider a planar or spherical, single-or dual-layer, material interface between a gelatin material and water. A single or several ideal-gas microbubbles are positioned near the interface. Bubble collapse is initiated by a shock wave with a pressure profile specific to laser generation and is flat when hitting the gas-water interface. The interfacial acoustic impedance match singles out the collapseinduced re-entrant jet as main event. High-resolution sharp-interface numerical methods are employed to ensure that wave dynamics, hydrodynamics, and interface transporting are accurately resolved. Bubble configurations are varied between single and double and between attached and with standoff distance. Parameters varied are shock-wave peak pressure and viscosity ratio between single and double layers of gelatin and the surrounding water. For inertia-dominated cases, two regimes are observed, the first characterized by linear growth of the penetration depth and the second by a t 2/3 scaling. The latter range is affected by viscosity which reduces penetration speed. The results show that process parameters, in particular shock overpressure, control not only penetration depth but also the size of the interface perforation, indicating means to steer processes in biomedical applications.
In this paper, a new fractional step method is proposed for simulating stiff and nonstiff chemically reacting flows. In stiff cases, a well-known spurious numerical phenomenon, i.e. the incorrect propagation speed of discontinuities, may be produced by general fractional step methods due to the under-resolved discretization in both space and time. The previous random projection method has been successfully applied for stiff detonation capturing in under-resolved conditions. Not to randomly project the intermediate state into two presumed equilibrium states (completely burnt or unburnt) as in the random projection method, the present study is to randomly choose the time-dependent advance or stop of a reaction process. Each one-way reaction has been decoupled from the multi-reaction kinetics using operator splitting and the local smeared temperature due to numerical dissipation of shock-capturing schemes is compared with a random one within two limited temperatures corresponding to the advance and its inverse states, respectively, to control the random reaction. The random activation or deactivation in the reaction step is thus promising to correct the deterministic accumulative error of the propagation of discontinuities. Extensive numerical experiments, including model problems and realistic reacting flows in one and two dimensions, demonstrate this expectation as well as the effectiveness and robustness of the method. Meanwhile, for nonstiff problems when spatial and temporal resolutions are fine, the proposed random method recovers the results as general fractional step methods, owing to the increasing possibility of activation with diminishing randomness by adding a shift term.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.