Numerical study of a planar shock interacting with a cylindrical water column embedded with an air cavity

2017 ◽  
Vol 825 ◽  
pp. 825-852 ◽  
Author(s):  
Gaoming Xiang ◽  
Bing Wang

This paper performs a numerical study on the interaction of a planar shock wave with a water column embedded with/without a cavity of different sizes at high Weber numbers. The conservative-type Euler and non-conservative scalar two-equations representing the transportation of two-phase properties consist of the diffusion interface capture models. The numerical fluxes are computed by the Godunov-type Harten-Lax–van Leer contact Riemann solver coupled with an incremental fifth-order weighted essentially non-oscillatory (WENO) scheme. A third-order total variation diminishing (TVD) Runge–Kutta scheme is used to advance the solution in time. The morphology and dynamical characteristics are analysed qualitatively and quantitatively to demonstrate the breakup mechanism of the water column and formation of transverse jets under different incident shock intensities and embedded-cavity sizes. The jet tip velocities are extracted by analysing the interface evolution. The liquid column is prone to aerodynamic breakup with the formation of micro-mist at later stages instead of liquid evaporation because of the weakly heating effects of the surrounding air. It is numerically confirmed that the liquid-phase pressure will drop below the saturated vapour pressure, and the low pressure can be sustained for a certain time because of the focusing of the expansion wave, which accounts for the cavitation inside the liquid water column. The geometrical parameters of the deformed water column are identified, showing that the centreline width decreases but the transverse height increases nonlinearly with time. The deformation rates are nonlinearly correlated under different Mach numbers. The first transverse jet is found for a water column with an embedded cavity, whereas the water hammer shock and second jet do not occur under the impact of low intensity incident shock waves. The $x$-velocity component recorded at the rear stagnation point can remain unchanged for a comparable time after a declined evolution, which indicates that the downstream wall of the shocked water ring somehow moves uniformly. It can be explained that the acceleration of the downstream wall is balanced by the trailing shedding vortex, and this effect is more evident under higher Mach numbers. The increased enstrophy, mainly generated at the interface, demonstrates the competition of the baroclinic effects of the shock wave impact over dilatation.

1996 ◽  
Vol 313 ◽  
pp. 105-130 ◽  
Author(s):  
O. Igra ◽  
J. Falcovitz ◽  
H. Reichenbach ◽  
W. Heilig

The interaction of a planar shock wave with a square cavity is studied experimentally and numerically. It is shown that such a complex, time-dependent, process can be modelled in a relatively simple manner. The proposed physical model is the Euler equations which are solved numerically, using the second-order-accurate high-resolution GRP scheme, resulting in very good agreement with experimentally obtained findings. Specifically, the wave pattern is numerically simulated throughout the entire interaction process. Excellent agreement is found between the experimentally obtained shadowgraphs and numerical simulations of the various flow discontinuities inside and around the cavity at all times. As could be expected, it is confirmed that the highest pressure acts on the cavity wall which experiences a head-on collision with the incident shock wave while the lowest pressures are encountered on the wall along which the incident shock wave diffracts. The proposed physical model and the numerical simulation used in the present work can be employed in solving shock wave interactions with other complex boundaries.


2021 ◽  
Vol 11 (11) ◽  
pp. 4934
Author(s):  
Viola Rossano ◽  
Giuliano De Stefano

Computational fluid dynamics was employed to predict the early stages of the aerodynamic breakup of a cylindrical water column, due to the impact of a traveling plane shock wave. The unsteady Reynolds-averaged Navier–Stokes approach was used to simulate the mean turbulent flow in a virtual shock tube device. The compressible flow governing equations were solved by means of a finite volume-based numerical method, where the volume of fluid technique was employed to track the air–water interface on the fixed numerical mesh. The present computational modeling approach for industrial gas dynamics applications was verified by making a comparison with reference experimental and numerical results for the same flow configuration. The engineering analysis of the shock–column interaction was performed in the shear-stripping regime, where an acceptably accurate prediction of the interface deformation was achieved. Both column flattening and sheet shearing at the column equator were correctly reproduced, along with the water body drift.


Author(s):  
S. M. French

Two damaged final reheat tubes from a 30 year old supercritical unit were submitted to the laboratory for evaluation following the discovering of a failure of one of the tubes after deslagging operations; a third, dented tube was left in service. The 304H stainless steel tubes were installed in 1990 when the reheater was replaced. The bulk microstructure of both tubes shows evidence of sensitization, which is not unusual given this application (reheater). The failed tube appears to be an intergranular separation that started either subsurface or at the ID, propagating to the OD surface. The sensitization of the steel apparently made the material susceptible to corrosion as well as significantly reduced the impact strength of the material to 10–15% of its estimated original level (verified by Charpy impact test). Examination of the dented tube (#101A) showed a subsurface plane of damage some 30 mils from the ID surface, running parallel to the surface. The damage consisted of intergranular separation, caused by the impact loading event, and referred to in the literature as an “attached spalling failure”. Spalling failures occur when the shock wave is reflected from the back surface (the ID surface of the tube), interacting with the incident shock wave as a stress wave. When the magnitude of this tensile stress exceeds the inherent strength of the material, failure occurs. The overall area of the attached spalling failure is being investigated; the concern there is if it is exceptionally large, it may provide a thermal barrier to heat transfer from the OD to the ID and result in a local overheating failure. Within the metallographic sample, however, the damage area was quite small and therefore did not appear to be an immediate issue. The long-term suitability of tube 105A, which remains in service with a dent induced by the same deslagging process that damaged tubes 101A and 103A, is doubtful and should be addressed during the Fall 2006 boiler overhaul. For the shortterm, the assumption was made that cracking due to the deslagging impact would be oriented similar to non-failed tube and extension of these fissures to failure between Spring 2006 and the Fall outage is not expected.


Author(s):  
H. T. C. Pedro ◽  
K.-W. Leung ◽  
M. H. Kobayashi ◽  
H. R. Riggs

This work concerns the numerical investigation of the impact of a wave on a square column. The wave is generated by a dam break in a wave tank. Two turbulence models were used: Large Eddy Simulations (LES) and Unsteady Reynolds Averaged Navier-Stokes (URANS). The numerical simulations were carried out using a finite volume approximation and the SIMPLE algorithm for the solution of the governing equations. Turbulence was modeled with the standard Smagorinsky-Lilly subgrid-model for the LES and the standard κ-ε model for the URANS. The results are validated against experimental data for the wave impact on a square column facing the flow. The results, especially for LES, show very good agreement between the predictions and experimental results. The overall accuracy of the LES, as expected, is superior to the URANS. However, if computational resources are limited, URANS can still provide satisfactory results for structural design.


Author(s):  
Matthieu Ancellin ◽  
Laurent Brosset ◽  
Jean-Michel Ghidaglia

Understanding the physics of sloshing wave impacts is necessary for the improvement of sloshing assessment methodology based on sloshing model tests, for LNG membrane tanks on floating structures. The phase change between natural gas and liquefied natural gas is one of the physical phenomena involved during a LNG wave impact but is not taken into account during sloshing model tests. In this paper, some recent numerical and analytical works on the influence of phase change are summarized and discussed. For the impact of an ideally shaped wave, phase change influences two different steps of the impact in different ways: during the gas escape phase, phase change leads to a higher impact velocity; for entrapped gas pockets, phase change causes a reduction of the pressure in the gas pocket. However, this influence is quantitatively small. The generalization to more realistic wave shapes (including e.g. liquid aeration) should be the focus of future works.


2017 ◽  
Vol 813 ◽  
pp. 70-84 ◽  
Author(s):  
V. Soni ◽  
A. Hadjadj ◽  
A. Chaudhuri ◽  
G. Ben-Dor

Numerical simulations were conducted to understand the different wave configurations associated with the shock-wave reflections over double-concave cylindrical surfaces. The reflectors were generated computationally by changing different geometrical parameters, such as the radii of curvature and the initial wedge angles. The incident-shock-wave Mach number was varied such as to cover subsonic, transonic and supersonic regimes of the flows induced by the incident shock. The study revealed a number of interesting wave features starting from the early stage of the shock interaction and transition to transitioned regular reflection (TRR) over the first concave surface, followed by complex shock reflections over the second one. Two new shock bifurcations have been found over the second wedge reflector, depending on the velocity of the additional wave that appears during the TRR over the first wedge reflector. Unlike the first reflector, the transition from a single-triple-point wave configuration (STP) to a double-triple-point wave configuration (DTP) and back occurred several times on the second reflector, indicating that the flow was capable of retaining the memory of the past events over the entire process.


2004 ◽  
Vol 126 (3) ◽  
pp. 399-409 ◽  
Author(s):  
A. Britan ◽  
A. V. Karpov ◽  
E. I. Vasilev ◽  
O. Igra ◽  
G. Ben-Dor ◽  
...  

The flow developed behind shock wave transmitted through a screen or a perforated plat is initially highly unsteady and nonuniform. It contains multiple shock reflections and interactions with vortices shed from the open spaces of the barrier. The present paper studies experimentally and theoretically/numerically the flow and wave pattern resulted from the interaction of an incident shock wave with a few different types of barriers, all having the same porosity but different geometries. It is shown that in all investigated cases the flow downstream of the barrier can be divided into two different zones. Due immediately behind the barrier, where the flow is highly unsteady and nonuniform in the other, placed further downstream from the barrier, the flow approaches a steady and uniform state. It is also shown that most of the attenuation experienced by the transmitted shock wave occurs in the zone where the flow is highly unsteady. When solving the flow developed behind the shock wave transmitted through the barrier while ignoring energy losses (i.e., assuming the fluid to be a perfect fluid and therefore employing the Euler equation instead of the Navier-Stokes equation) leads to non-physical results in the unsteady flow zone.


Sign in / Sign up

Export Citation Format

Share Document