Simulation of plunging wave impact on a vertical wall

1996 ◽  
Vol 327 ◽  
pp. 221-254 ◽  
Author(s):  
Sheguang Zhang ◽  
Dick K. P. Yue ◽  
Katsuji Tanizawa

We present a numerical study of the impact of a two-dimensional plunging wave on a rigid vertical wall in the context of potential flow. The plunging wave impinging the wall is generated using a mixed-Eulerian-Lagrangian (MEL) boundary-integral scheme. The initial stage of the impact is characterized by an oblique impact of a liquid wedge on the wall and is solved using a similarity solution. Following the initial impact, the MEL simulation is continued to capture the transient impact process. The effect of an air cushion trapped between the plunger and the wall is considered. In addition to details such as temporal evolutions and surface profiles, the main interests are the maximum impact pressure on the wall and its rise time. To arrive at appropriate scaling laws for these, simulations are performed and correlations are explored for a broad range of local plunging wave kinematic and geometric parameters. To assess the present results, direct comparisons are made with the experiment of Chan & Melville (1988). Reasonable quantitative agreement is obtained and likely sources for discrepancies are identified and discussed.

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.


2020 ◽  
Vol 61 (9) ◽  
Author(s):  
T. Mai ◽  
C. Mai ◽  
A. Raby ◽  
D. M. Greaves

Abstract Local and global loadings, which may cause the local damage and/or global failure and collapse of offshore structures and ships, are experimentally investigated in this study. The research question is how the elasticity of the structural section affects loading during severe environmental conditions. Two different experiments were undertaken in this study to try to answer this question: (i) vertical slamming impacts of a square flat plate, which represents a plate section of the bottom or bow of a ship structure, onto water surface with zero degree deadrise angle; (ii) wave impacts on a truncated vertical wall in water, where the wall represents a plate section of a hull. The plate and wall are constructed such that they can be either rigid or elastic by virtue of a specially designed spring system. The experiments were carried out in the University of Plymouth’s COAST Laboratory. For the cases considered here, elasticity of the impact plate and/or wall has an effect on the slamming and wave impact loads. Here the slamming impact loads (both pressure and force) were considerably reduced for the elastic plate compared to the rigid one, though only at high impact velocities. The total impact force on the elastic wall was found to reduce for the high aeration, flip-through and slightly breaking wave impacts. However, the impact pressure decreased on the elastic wall only under flip-through wave impact. Due to the elasticity of the plates, the impulse of the first positive phase of pressure and force decreases significantly for the vertical slamming impact tests. This significant effect of hydroelasticity is also found for the total force impulse on the vertical wall under wave impacts. Graphic abstract Hydroelasticity effects on water-structure impacts: a impact pressures on dropped plates; b impact forces on dropped plates; c, d, e, f wave impact pressures on the vertical walls; g wave impact forces on the vertical walls; h wave force impulses on the vertical walls: elastic wall 1 vs. rigid wall (filled markers); elastic wall 2 vs. rigid wall (empty markers)


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.


2015 ◽  
Vol 12 (11) ◽  
pp. 2045-2060 ◽  
Author(s):  
Buyun Su ◽  
Zhiwei Zhou ◽  
Jianjun Zhang ◽  
Zhihua Wang ◽  
Xuefeng Shu ◽  
...  

2017 ◽  
Author(s):  
Xizeng Zhao ◽  
Yong Chen ◽  
Zhenhua Huang ◽  
Yangyang Gao

Abstract. There is a general lack of the understanding of tsunami wave interacting with complex geographies, especially the process of inundation. Numerical simulations are performed to understand the effects of several factors on tsunami wave impact and run-up in the presence of submarine gentle slopes and coastal cliffs, using an in-house code, named a Constrained Interpolation Profile (CIP)-based model in Zhejiang University (CIP-ZJU). The model employs a high-order finite difference method, the CIP method as the flow solver, utilizes a VOF-type method, the Tangent of hyperbola for interface capturing/Slope weighting (THINC/SW) scheme to capture the free surface, and treats the solid boundary by an immersed boundary method. A series of incident waves are arranged to interact with varying coastal geographies. Numerical results are compared with experimental data and good agreement is obtained. The influences of submarine gentle slope, coastal cliff and incident wave height are discussed. It is found that the rule of tsunami amplification factor varying with incident wave is affected by angle of cliff slope, and there is a critical angle about 45°. The run-up on a toe-erosion cliff is smaller than that on a normal cliff. The run-up is also related to the length of submarine gentle slope with a critical about 2.292 m in the present study. The impact pressure on the cliff is extremely large and concentrated, and the backflow effect is nonnegligible. Results of our work are in high precision and helpful in inversing tsunami source and forecasting disaster.


2016 ◽  
Vol 792 ◽  
pp. 850-868 ◽  
Author(s):  
Wilco Bouwhuis ◽  
Xin Huang ◽  
Chon U Chan ◽  
Philipp E. Frommhold ◽  
Claus-Dieter Ohl ◽  
...  

A train of high-speed microdrops impacting on a liquid pool can create a very deep and narrow cavity, reaching depths more than 1000 times the size of the individual drops. The impact of such a droplet train is studied numerically using boundary integral simulations. In these simulations, we solve the potential flow in the pool and in the impacting drops, taking into account the influence of liquid inertia, gravity and surface tension. We show that for microdrops the cavity shape and maximum depth primarily depend on the balance of inertia and surface tension and discuss how these are influenced by the spacing between the drops in the train. Finally, we derive simple scaling laws for the cavity depth and width.


Bioimpacts ◽  
2017 ◽  
Vol 8 (4) ◽  
pp. 271-279
Author(s):  
Mohammad T. Shervani-Tabar ◽  
Babak Farzaneh ◽  
Reza Ahrabi ◽  
Seyed E. Razavi

Introduction: Employing of gaseous plugs inside a vein for preventing of blood flow to the damaged or cancerous tissues has been known as a gas embolism in the medicine. In this research, a numerical investigation has been carried out on the delivery of the liquid drug DDFP, encapsulated in the microlipid-coated spheres (MLCSs), to target the human vein for construction of the gaseous plug inside the veins. Methods: The encapsulated liquid drug DDFP were delivered to the vein by injection of an emulsion. Releasing of the liquid drug DDFP results in an explosive growth of a gaseous plug inside the vein. The targeted vein was served as a rigid cylinder with a compliant coating. The boundary integral equation method has been employed for the numerical simulation of the hydrodynamic behavior of the gaseous plug inside the vein. Results: Numerical results showed that in the case of a rigid cylinder vein with an internal compliant coating, the maximum volume of the gaseous plug was smaller than the case of just a rigid cylinder vein. Furthermore, its elapsed time from the instant of bubble generation to the instant when the bubble reaches its maximum volume was shorter. Numerical results also showed that the compliant coating on the internal wall of the rigid cylindrical vein had a tendency of reducing the impact of the explosive growth of the gaseous plug. Conclusion: This numerical research showed that the compliant coating on the internal wall of the rigid cylindrical vein had the tendency of reducing the impact of the impulsive growth of the gaseous plug. Therefore, in the case of having severed arteriosclerosis, treatment of the cancerous or damaged tissues by use of the gaseous embolism must be done very carefully in order to prevent the hazardous effects of the gaseous plug’s impulsive growth.


2021 ◽  
pp. 174425912098418
Author(s):  
Toivo Säwén ◽  
Martina Stockhaus ◽  
Carl-Eric Hagentoft ◽  
Nora Schjøth Bunkholt ◽  
Paula Wahlgren

Timber roof constructions are commonly ventilated through an air cavity beneath the roof sheathing in order to remove heat and moisture from the construction. The driving forces for this ventilation are wind pressure and thermal buoyancy. The wind driven ventilation has been studied extensively, while models for predicting buoyant flow are less developed. In the present study, a novel analytical model is presented to predict the air flow caused by thermal buoyancy in a ventilated roof construction. The model provides means to calculate the cavity Rayleigh number for the roof construction, which is then correlated with the air flow rate. The model predictions are compared to the results of an experimental and a numerical study examining the effect of different cavity designs and inclinations on the air flow rate in a ventilated roof subjected to varying heat loads. Over 80 different test set-ups, the analytical model was found to replicate both experimental and numerical results within an acceptable margin. The effect of an increased total roof height, air cavity height and solar heat load for a given construction is an increased air flow rate through the air cavity. On average, the analytical model predicts a 3% higher air flow rate than found in the numerical study, and a 20% lower air flow rate than found in the experimental study, for comparable test set-ups. The model provided can be used to predict the air flow rate in cavities of varying design, and to quantify the impact of suggested roof design changes. The result can be used as a basis for estimating the moisture safety of a roof construction.


Nanomaterials ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 1425
Author(s):  
Tarek Bouzennada ◽  
Farid Mechighel ◽  
Kaouther Ghachem ◽  
Lioua Kolsi

A 2D-symmetric numerical study of a new design of Nano-Enhanced Phase change material (NEPCM)-filled enclosure is presented in this paper. The enclosure is equipped with an inner tube allowing the circulation of the heat transfer fluid (HTF); n-Octadecane is chosen as phase change material (PCM). Comsol-Multiphysics commercial code was used to solve the governing equations. This study has been performed to examine the heat distribution and melting rate under the influence of the inner-tube position and the concentration of the nanoparticles dispersed in the PCM. The inner tube was located at three different vertical positions and the nanoparticle concentration was varied from 0 to 0.06. The results revealed that both heat transfer/melting rates are improved when the inner tube is located at the bottom region of the enclosure and by increasing the concentration of the nanoparticles. The addition of the nanoparticles enhances the heat transfer due to the considerable increase in conductivity. On the other hand, by placing the tube in the bottom area of the enclosure, the liquid PCM gets a wider space, allowing the intensification of the natural convection.


Sign in / Sign up

Export Citation Format

Share Document