Control of the dynamics of a boiling vapour bubble using pressure-modulated high intensity focused ultrasound without the shock scattering effect: A first proof-of-concept study

Boiling histotripsy is a promising High-Intensity Focused Ultrasound (HIFU) technique that can be used to induce mechanical tissue fractionation at the HIFU focus via cavitation. Two different types of cavitation produced during boiling histotripsy exposure can contribute towards mechanical tissue destruction: (1) a boiling vapour bubble at the HIFU focus and (2) cavitation clouds in between the boiling bubble and the HIFU source. Control of the extent and degree of mechanical damage produced by boiling histotripsy is necessary when treating a solid tumour adjacent to normal tissue or major blood vessels. This is, however, difficult to achieve with boiling histotripsy due to the stochastic formation of the shock scattering-induced inertial cavitation clouds. In the present study, a new histotripsy method termed pressure-modulated shockwave histotripsy is proposed as an alternative to or in addition to boiling histotripsy without inducing the shock scattering effect. The proposed concept is (a) to generate a boiling vapour bubble via localised shockwave heating and (b) subsequently control its extent and lifetime through manipulating peak pressure magnitudes and a HIFU pulse length. To demonstrate the feasibility of the proposed method, bubble dynamics induced at the HIFU focus in an optically transparent liver tissue phantom were investigated using a high speed camera and a passive cavitation detection systems under a single 10, 50 or 100 ms-long 2, 3.5 or 5 MHz pressure-modulated HIFU pulse with varying peak positive and negative pressure amplitudes from 5 to 89 MPa and −3.7 to −14.6 MPa at the focus. Furthermore, a numerical simulation of 2D nonlinear wave propagation with the presence of a boiling bubble at the focus of a HIFU field was conducted by numerically solving the generalised Westervelt equation. The high speed camera experimental results showed that, with the proposed pressure-modulated shockwave histotripsy, boiling bubbles generated by shockwave heating merged together, forming a larger bubble (of the order of a few hundred micron) at the HIFU focus. This coalesced boiling bubble then persisted and maintained within the HIFU focal zone until the end of the exposure (10, 50, or 100 ms). Furthermore, and most importantly, no violent cavitation clouds which typically appear in boiling histotripsy occurred during the proposed histotripsy excitation (i.e. no shock scattering effect). This was likely because that the peak negative pressure magnitude of the backscattered acoustic field by the boiling bubble was below the cavitation cloud intrinsic threshold. The size of the coalesced boiling bubble gradually increased with the peak pressure magnitudes. In addition, with the proposed method, an oval shaped lesion with a length of 0.6 mm and a width of 0.1 mm appeared at the HIFU focus in the tissue phantom, whereas a larger lesion in the form of a tadpole (length: 2.7 mm, width: 0.3 mm) was produced by boiling histotripsy. Taken together, these results suggest that the proposed pressure-modulated shockwave histotripsy could potentially be used to induce a more spatially localised tissue destruction with a desired degree of mechanical damage through controlling the size and lifetime of a boiling bubble without the shock scattering effect.     

Moving shocks through metallic grids: their interaction and potential for blast wave mitigation

Numerical simulations and laboratory measurements have been used to illuminate the interaction of a moving shock wave impacting on metallic grids at various shock strengths and grid solidities. The experimental work was carried out in a large scale shock tube facility while computational work simulated the flow field with a time-dependent inviscid and a time-dependent viscous model. The pressure drop measured across the grids is a result of two phenomena which are associated with the impact of the shock on the metallic grids. First are the reflection and refraction of the incoming shock on the grid itself. This appears to be the main inviscid mechanism associated with the reduction of the strength of the transmitted shock. Second, viscous phenomena are present during the reflection and refraction of the wave as well as during the passage of the induced flow of the air through the grid. The experimental data of pressure drop across the grid obtained in the present investigation are compared with those obtained from computations. The numerical results slightly overpredict the experimental data of relative pressure drop which increases substantially with grid solidity at fixed flow Mach numbers. The processes of shock reflection and refraction are continuous and they can be extended in duration by using thicker grids that will result in lower compression rates of the structural loading and increase the viscous losses associated with these phenomena which will further attenuate the impacting shock. Preliminary theoretical analysis suggests that the use of a graded porosity/solidity material will result in higher pressure drop than a constant porosity/solidity material and thus provide effective blast mitigation.      

Unidimensional SPH simulations of reactive shock tubes in an astrophysical perspective

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method widely used for the modelling of a large variety of astrophysical fluid flows in more than one dimension. Simulations of thermonuclear explosions in stars require, besides the hydrodynamic equations, a realistic equation of state, an energy source term, and a set of nuclear kinetic equations to follow the composition changes of the gas during the explosion. The implementation of a realistic stellar equation of state, and the coupling of hydrodynamics and nuclear burning are investigated in the framework of the simple shock tube geometry. We present and discuss the results of a series of SPH simulations of a detonation in the presence of (1) a single exothermic nuclear reaction, and (2) a restricted network of nuclear reactions. Our results are compared to those of identical simulations performed by other authors using a different hydrodynamic method.     

Shock tunnel flow visualization using planar laser-induced fluorescence imaging of NO and OH

Temporal sequences of planar laser-induced fluorescence (PLIF) images of several high-speed, transient flowfields created in a reflection-type shock tunnel facility were acquired. In each case, the test gas contained either nitric oxide or the hydroxyl radical, the fluorescent species. The processes of shock reflection from an endwall with a converging nozzle and of underexpanded free jet formation were examined. A comparison was also made between PLIF imaging and shadow photography. The investigation demonstrated some of the capabilities of PLIF imaging diagnostics in complex, transient, hypersonic flowfields, including those with combustion.Nomenclature A spontaneous emission rate - A _las cross sectional area of laser sheet - B laser absorption rate - C _opt constant dependent on optical arrangement, collection efficiency, etc. - D nozzle throat diameter - E _p laser pulse energy - f _J Boltzmann fraction of absorbing state - g spectral convolution of laser and absorption lineshapes - k Boltzmann constant - M _s incident shock Mach number - N noise, root-mean-square signal fluctuation - P static pressure - P ₁ initial pressure of test gas in shock tube - P _a free jet ambient pressure - P _s stagnation pressure - Q electronic quenching rate of excited state - S PLIF signal - t time between shock reflection and image acquisition - T static temperature - T _s stagnation temperature - _a mole fraction of absorbing species     

Limit states of elastic–viscoplastic plate deformations caused by repeated shock wave-loadings. Part 1: Experimental observation

In the present work, the evolution of the inelastic centre deflections of shock wave-loaded circular metal plates due to repeated loadings is studied experimentally and numerically. These displacements are compared to those of quasi-statically deformed plates loaded by a pressure equal to the peak pressure of the impulsive loading. Thereby three types of permanent centre deflections are observed: (1) The quasi-statically obtained deflection is exceeded by the middle point displacement of a dynamically loaded structure already after the first impulse and tends towards a limit state after repeated shock wave-loadings. (2) The centre deflection of the impulsively loaded plate exceeds also the quasi-statically caused deflection and does not increase after repeated impulsive loadings any more. (3) The permanent middle point displacement of a dynamically loaded plate is smaller than the deflection of a quasi-statically loaded one and tends towards the middle point displacement of the quasi-static counterpart after repeated shock wave loadings. This phenomenon is known in the literature as ‘Pseudo-shakedown’.

In Part 1 of this study the experimental observation is described, followed by a theoretical study in Part 2.     

Experimental investigation on tunnel sonic boom

Upon the entrance of a high-speed train into a relatively long train tunnel, compression waves are generated in front of the train. These compression waves subsequently coalesce into a weak shock wave so that a unpleasant sonic boom is emitted from the tunnel exit. In order to investigate the generation of the weak shock wave in train tunnels and the emission of the resulting sonic boom from the train tunnel exit and to search for methods for the reduction of these sonic booms, a 1300 scaled train tunnel simulator was constructed and simulation experiments were carried out using this facility.In the train tunnel simulator, an 18 mm dia. and 200 mm long plastic piston moves along a 40 mm dia. and 25 m long test section with speed ranging from 60 to 100 m/s. The tunnel simulator was tilted 8° to the floor so that the attenuation of the piston speed was not more than 10 % of its entrance speed. Pressure measurements along the tunnel simulator and holographic interferometric optical flow visualization of weak shock waves in the tunnel simulator clearly showed that compression waves, with propagation, coalesced into a weak shock wave. Although, for reduction of the sonic boom in prototype train tunnels, the installation of a hood at the entrance of the tunnels was known to be useful for their suppression, this effect was confirmed in the present experiment and found to be effective particularly for low piston speeds. The installation of a partially perforated wall at the exit of the tunnel simulator was found to smear pressure gradients at the shock. This effect is significant for higher piston speeds. Throughout the series of train tunnel simulator experiments, the combination of both the entrance hood and the perforated wall significantly reduces shock overpressures for piston speeds ofu _p ranging from 60 to 100 m/s. These experimental findings were then applied to a real train tunnel and good agreement was obtained between the tunnel simulator result and the real tunnel measurements.     

The reflected pressure field in the interaction of weak shock waves with a compressible foam

This paper deals with the waves that are reflected from slabs of porous compressible foam attached to a rigid wall when impacted by a weak shock wave. The interest is in establishing possible attenuation of the pressure field after a shock or blast wave has struck the surface. Foam densities from 14 to 38 kg/m³ were tested over a range of shock wave Mach numbers less than 1.4. It is shown that the initial reflected shock wave strength is accurately predicted by the pseudo-gas model of Gelfand et al. (1983), with a pressure ratio of approximately 80% of the value for reflection off a rigid wall. Evidence is presented of gas entering the foam during the early stages of the process. A second wave emerges from the foam at a later stage and is separated from the first by a region of constant velocity and pressure. This second wave is not a shock wave but a compression front of significant thickness, which emerges from the foam earlier than predicted by the pseudo-gas model. Analysis of the origin of this wave points to much more complex flows within the foam than previously assumed, particularly in an apparent decrease in average wave front speed as the foam is compressed. It is shown that the pressure ratio across both these waves taken together is slightly higher than that for reflection off a rigid wall. In some cases this compression wave train is followed by a weak expansion wave.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Over forty years of continuous research at UTIAS on nonstationary flows and shock waves

Analytical and experimental research on non-stationary shock waves, rarefaction waves and contact surfaces has been conducted continuously at UTIAS since its inception in 1948. Some unique facilities were used to study the properties of planar, cylindrical and spherical shock waves and their interactions. Investigations were also performed on shock-wave structure and boundary layers in ionizing argon, water-vapour condensation in rarefaction waves, magnetogasdynamic flows, and the regions of regular and various types of Mach reflections of oblique shock waves. Explosively-driven implosions have been employed as drivers for projectile launchers and shock tubes, and as a means of producing industrial-type diamonds from graphite, and fusion plasmas in deuterium. The effects of sonic-boom on humans, animals and structures have also formed an important part of the investigations. More recently, interest has focussed on shock waves in dusty gases, the viscous and vibrational structure of weak spherical blast waves in air, and oblique shock-wave reflections. In all of these studies instrumentation and computational methods have played a very important role. A brief survey of this work is given herein and in more detail in the relevant references.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

GPU accelerated CESE method for 1D shock tube problems

In the present study, a GPU accelerated 1D space–time CESE method is developed and applied to shock tube problems with and without condensation. We have demonstrated how to implement the CESE algorithm to solve 1D shock tube problems using an older generation GPU (the NVIDIA 9800 GT) with relatively limited memory. To optimize the code performance, we used Shared Memory and solved the inter-Block boundary problem in two ways, namely the branch scheme and the overlapping scheme. The implementations of these schemes are discussed in detail and their performances are compared for the Sod shock tube problems. For the Sod problem without condensation, the speedup over an Intel CPU E7300 is 23 for the branch scheme and 41 for the overlapping scheme, respectively. While for problems with condensation, both schemes achieve higher acceleration ratios, 53 and 71, respectively. The higher speedup of the condensation case can be ascribed to the source term calculation which has a local dependence on the mesh point and the SOURCE kernel has a higher acceleration ratio.     

Adaptive hierarchic transformations for dynamically p-enriched slope-limiting over discontinuous Galerkin systems of generalized equations

We study a family of generalized slope limiters in two dimensions for Runge–Kutta discontinuous Galerkin (RKDG) solutions of advection-diffusion systems. We analyze the numerical behavior of these limiters applied to a pair of model problems, comparing the error of the approximate solutions, and discuss each limiter’s advantages and disadvantages. We then introduce a series of coupled p-enrichment schemes that may be used as standalone dynamic p-enrichment strategies, or may be augmented via any in the family of variable-in-p slope limiters presented.