No abstract

a b s t r a c tFlows in which shock waves and turbulence are present and interact dynamically occur in a wide range of applications, including inertial confinement fusion, supernovae explosion, and scramjet propulsion. Accurate simulations of such problems are challenging because of the contradictory requirements of numerical methods used to simulate turbulence, which must minimize any numerical dissipation that would otherwise overwhelm the small scales, and shock-capturing schemes, which introduce numerical dissipation to stabilize the solution. The objective of the present work is to evaluate the performance of several numerical methods capable of simultaneously handling turbulence and shock waves. A comprehensive range of high-resolution methods (WENO, hybrid WENO/central difference, artificial diffusivity, adaptive characteristic-based filter, and shock fitting) and suite of test cases (Taylor-Green vortex, Shu-Osher problem, shock-vorticity/entropy wave interaction, Noh problem, compressible isotropic turbulence) relevant to problems with shocks and turbulence are considered. The results indicate that the WENO methods provide sharp shock profiles, but overwhelm the physical dissipation. The hybrid method is minimally dissipative and leads to sharp shocks and well-resolved broadband turbulence, but relies on an appropriate shock sensor. Artificial diffusivity methods in which the artificial bulk viscosity is based on the magnitude of the strain-rate tensor resolve vortical structures well but damp dilatational modes in compressible turbulence; dilatation-based artificial bulk viscosity methods significantly improve this behavior. For well-defined shocks, the shock fitting approach yields good results.

Histotripsy is an ultrasound ablation method that depends on the initiation of a cavitation bubble cloud to fractionate soft tissue. Previous work has demonstrated a cavitation cloud can be formed by a single pulse with one high amplitude negative cycle, when the negative pressure amplitude directly exceeds a pressure threshold intrinsic to the medium. We hypothesize that the intrinsic threshold in water-based tissues is determined by the properties of the water inside the tissue and changes in tissue stiffness or ultrasound frequency will have a minimal impact on the histotripsy intrinsic threshold. To test this hypothesis, the histotripsy intrinsic threshold was investigated both experimentally and theoretically. The probability of cavitation was measured by subjecting tissue phantoms with adjustable mechanical properties and ex vivo tissues to a histotripsy pulse of 1–2 cycles produced by 345 kHz, 500 kHz, 1.5 MHz, and 3 MHz histotripsy transducers. Cavitation was detected and characterized by passive cavitation detection and high-speed photography, from which the probability of cavitation was measured vs. pressure amplitude. The results demonstrated that the intrinsic threshold (the negative pressure at which probability=0.5) is independent of stiffness for Young’s moduli (E) < 1 MPa with only a small increase (~2–3 MPa) in the intrinsic threshold for tendon (E=380 MPa). Additionally, results for all samples showed only a small increase of ~2–3 MPa when the frequency was increased from 345 kHz to 3 MHz. The intrinsic threshold was measured to be between 24.7–30.6 MPa for all samples and frequencies tested in this study. Overall, the results of this study indicate that the intrinsic threshold to initiate a histotripsy bubble cloud is not significantly impacted by tissue stiffness or ultrasound frequency in hundreds of kHz to MHz range.

Atmospheric pressure plasma jets (APPJs) are being investigated in the context plasma medicine and biotechnology applications, and surface functionalization. The composition of the surface being treated ranges from plastics, liquids, and biological tissue, to metals. The dielectric constant of these materials ranges from as low as 1.5 for plastics to near 80 for liquids, and essentially infinite for metals. The electrical properties of the surface are not independent variables as the permittivity of the material being treated has an effect on the dynamics of the incident APPJ. In this paper, results are discussed from a computational investigation of the interaction of an APPJ incident onto materials of varying permittivity, and their impact on the discharge dynamics of the plasma jet. The computer model used in this investigation solves Poisson's equation, transport equations for charged and neutral species, the electron energy equation, and the Navier-Stokes equations for the neutral gas flow. The APPJ is sustained in He/O 2 ¼ 99.8/0.2 flowing into humid air, and is directed onto dielectric surfaces in contact with ground with dielectric constants ranging from 2 to 80, and a grounded metal surface. Low values of relative permittivity encourage propagation of the electric field into the treated material and formation and propagation of a surface ionization wave. High values of relative permittivity promote the restrike of the ionization wave and the formation of a conduction channel between the plasma discharge and the treated surface. The distribution of space charge surrounding the APPJ is discussed. V C 2015 AIP Publishing LLC.

A high-order accurate shock- and interface-capturing scheme is used to simulate the collapse of a gas bubble in water. In order to better understand the damage caused by collapsing bubbles, the dynamics of the shock-induced and Rayleigh collapse of a bubble near a planar rigid surface and in a free field are analysed. Collapse times, bubble displacements, interfacial velocities and surface pressures are quantified as a function of the pressure ratio driving the collapse and of the initial bubble stand-off distance from the wall; these quantities are compared to the available theory and experiments and show good agreement with the data for both the bubble dynamics and the propagation of the shock emitted upon the collapse. Non-spherical collapse involves the formation of a re-entrant jet directed towards the wall or in the direction of propagation of the incoming shock. In shock-induced collapse, very high jet velocities can be achieved, and the finite time for shock propagation through the bubble may be non-negligible compared to the collapse time for the pressure ratios of interest. Several types of shock waves are generated during the collapse, including precursor and water-hammer shocks that arise from the re-entrant jet formation and its impact upon the distal side of the bubble, respectively. The water-hammer shock can generate very high pressures on the wall, far exceeding those from the incident shock. The potential damage to the neighbouring surface is quantified by measuring the wall pressure. The range of stand-off distances and the surface area for which amplification of the incident shock due to bubble collapse occurs is determined.