scholarly journals Numerical Study of Violent Impact Flow Using a CIP-Based Model

2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Qiao-ling Ji ◽  
Xi-zeng Zhao ◽  
Sheng Dong

A two-phase flow model is developed to study violent impact flow problem. The model governed by the Navier-Stokes equations with free surface boundary conditions is solved by a Constrained Interpolation Profile (CIP)-based high-order finite difference method on a fixed Cartesian grid system. The free surface is immersed in the computation domain and expressed by a one-fluid density function. An accurate Volume of Fluid (VOF)-type scheme, the Tangent of Hyperbola for Interface Capturing (THINC), is combined for the free surface treatment. Results of another two free surface capturing methods, the original VOF and CIP, are also presented for comparison. The validity and utility of the numerical model are demonstrated by applying it to two dam-break problems: a small-scale two-dimensional (2D) and three-dimensional (3D) full scale simulations and a large-scale 2D simulation. Main attention is paid to the water elevations and impact pressure, and the numerical results show relatively good agreement with available experimental measurements. It is shown that the present numerical model can give a satisfactory prediction for violent impact flow.

Energies ◽  
2019 ◽  
Vol 12 (22) ◽  
pp. 4244 ◽  
Author(s):  
Shuo Mi ◽  
Zongliu Huang ◽  
Xin Jin ◽  
Mahdi Tabatabaei Malazi ◽  
Mingming Liu

The highly viscous liquid (glycerin) sloshing is investigated numerically in this study. The full-scale membrane-type tank is considered. The numerical investigation is performed by applying a two-phase numerical model based on the spatially averaged Navier-Stokes equations. Firstly, the numerical model is validated against the available numerical model and a self-conducted experiment then is applied to systematically investigate the full-scale sloshing. In this study, two filling levels (50% and 70% of the tank height) are considered. The fluid kinematic viscosity is fixed at a value being 6.0 × 10−5 m2/s with comparative value to that of the crude oil. A wide range of forcing periods varying from 8.0 s to 12.0 s are used to identify the response process of pressures as well as free surface displacements. The pressures are analyzed along with breaking free surface snapshots and corresponding pressure distributions. The slamming effects are also demonstrated. Finally, the frequency response is further identified by the fast Fourier transformation technology.


Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 350
Author(s):  
Georgios A. Leftheriotis ◽  
Iason A. Chalmoukis ◽  
Guillermo Oyarzun ◽  
Athanassios A. Dimas

An advanced numerical model is presented for the simulation of wave-induced free-surface flow, utilizing an efficient hybrid parallel implementation. The model is based on the solution of the Navier–Stokes equations using large-eddy simulation of large-scale coastal free-surface flows. The three-dimensional immersed boundary method was used for the enforcement of the no-slip boundary condition on the bed surface. The water-air interface was tracked using the level-set method. The numerical model was effectively validated against laboratory measurements involving wave propagation over a flatbed with an elliptical shoal, whose presence induces combined wave refraction and diffraction phenomena. The parallel implementation of the model enabled the efficient simulation of depth-resolved, wave-induced, three-dimensional, free-surface flow; the model parallel efficiency and strong scaling are quantitatively demonstrated.


Energies ◽  
2021 ◽  
Vol 14 (3) ◽  
pp. 639
Author(s):  
Dennis Powalla ◽  
Stefan Hoerner ◽  
Olivier Cleynen ◽  
Nadine Müller ◽  
Jürgen Stamm ◽  
...  

The objective of the present paper is to develop a validated numerical model of a water vortex power plant that serves as a digital twin for further studies such as assessments of the ethohydraulic characteristics or the performance of such devices. The reference for the validation process is a large-scale hydraulic installation equipped with a full-scale water vortex power plant prototype installed in Dresden (Germany), where flow field measurements were carried out using three-dimensional Acoustic Doppler Velocimetry. The numerical model was implemented within the software package Star-CCM+. The unsteady, two-phase flow was solved with the Reynolds-Averaged Navier–Stokes equations in a Eulerian Multiphase approach, deploying a Volume of Fluid method to describe the free-surface flow. Water level and flow velocities were systematically compared in key areas of the device, demonstrating that the simulation is in good agreement with experimental observations. Relative differences are limited to at most 4% regarding water height in the system, and even the much more challenging velocity fields are reproduced with typical relative errors of roughly 10%. This validates the ability of the model to model the challenging flow conditions found in a water vortex power plant, enabling subsequent studies of the characteristics of this power plant concerning fish migration.


Water ◽  
2019 ◽  
Vol 11 (12) ◽  
pp. 2511
Author(s):  
Jintao Liu ◽  
Di Xu ◽  
Shaohui Zhang ◽  
Meijian Bai

This paper investigates the physical processes involved in the water filling and air expelling process of a pipe with multiple air valves under water slow filling condition, and develops a fully coupledwater–air two-phase stratified numerical model for simulating the process. In this model, the Saint-Venant equations and the Vertical Average Navier–Stokes equations (VANS) are respectively applied to describe the water and air in pipe, and the air valve model is introduced into the VANS equations of air as the source term. The finite-volume method and implicit dual time-stepping method (IDTS) with two-order accuracy are simultaneously used to solve this numerical model to realize the full coupling between water and air movement. Then, the model is validated by using the experimental data of the pressure evolution in pipe and the air velocity evolution of air valves, which respectively characterize the water filling and air expelling process. The results show that the model performs well in capturing the physical processes, and a reasonable agreement is obtained between numerical and experimental results. This agreement demonstrates that the proposed model in this paper offers a practical method for simulating water filling and air expelling process in a pipe with multiple air valves under water slow filling condition.


1999 ◽  
Vol 396 ◽  
pp. 37-71 ◽  
Author(s):  
LEONID BREVDO ◽  
PATRICE LAURE ◽  
FREDERIC DIAS ◽  
THOMAS J. BRIDGES

The film flow down an inclined plane has several features that make it an interesting prototype for studying transition in a shear flow: the basic parallel state is an exact explicit solution of the Navier–Stokes equations; the experimentally observed transition of this flow shows many properties in common with boundary-layer transition; and it has a free surface, leading to more than one class of modes. In this paper, unstable wavepackets – associated with the full Navier–Stokes equations with viscous free-surface boundary conditions – are analysed by using the formalism of absolute and convective instabilities based on the exact Briggs collision criterion for multiple k-roots of D(k, ω) = 0; where k is a wavenumber, ω is a frequency and D(k, ω) is the dispersion relation function.The main results of this paper are threefold. First, we work with the full Navier–Stokes equations with viscous free-surface boundary conditions, rather than a model partial differential equation, and, guided by experiments, explore a large region of the parameter space to see if absolute instability – as predicted by some model equations – is possible. Secondly, our numerical results find only convective instability, in complete agreement with experiments. Thirdly, we find a curious saddle-point bifurcation which affects dramatically the interpretation of the convective instability. This is the first finding of this type of bifurcation in a fluids problem and it may have implications for the analysis of wavepackets in other flows, in particular for three-dimensional instabilities. The numerical results of the wavepacket analysis compare well with the available experimental data, confirming the importance of convective instability for this problem.The numerical results on the position of a dominant saddle point obtained by using the exact collision criterion are also compared to the results based on a steepest-descent method coupled with a continuation procedure for tracking convective instability that until now was considered as reliable. While for two-dimensional instabilities a numerical implementation of the collision criterion is readily available, the only existing numerical procedure for studying three-dimensional wavepackets is based on the tracking technique. For the present flow, the comparison shows a failure of the tracking treatment to recover a subinterval of the interval of unstable ray velocities V whose length constitutes 29% of the length of the entire unstable interval of V. The failure occurs due to a bifurcation of the saddle point, where V is a bifurcation parameter. We argue that this bifurcation of unstable ray velocities should be observable in experiments because of the abrupt increase by a factor of about 5.3 of the wavelength across the wavepacket associated with the appearance of the bifurcating branch. Further implications for experiments including the effect on spatial amplification rate are also discussed.


2020 ◽  
Author(s):  
Nicholas S. Tavouktsoglou ◽  
Aggelos Dimakopoulos ◽  
Jeremy Spearman ◽  
Richard J. S. Whitehouse

Abstract Submerged water jet causing soil excavation is a typical water-soil interaction process that occurs widely in many engineering disciplines. In hydraulic engineering for instance, a typical example would be scour downstream of headcuts, culverts, or dam spillways. In port and waterway engineering, erosion of the channel bed or quay wall by the propellers of passing ships are also typical water jet/soil interaction problems. In ocean engineering, trenching by impinging high-velocity water jets has been used as an efficient method for cable and pipeline burial. At present, physical modelling and simple prediction equations have been the main practical engineering tool for evaluating scour in these situations. However, with the increasing computational power of modern computers and the development of new Computational Fluid Dynamics (CFD) solvers, scour prediction in such engineering problems has become possible. In the present work three-dimensional (3D) numerical modelling has been applied to reproduce the capability of a pair of water jets to backfill an excavated trench. The simulations are carried out using a state-of-the-art three-dimensional Eulerian two-phase scour model based on the open source CFD software OpenFOAM. The fluid phase is resolved by solving modified Navier-Stokes equations, which take into consideration the influence of the solid phase, i.e., the soil particles. This paper first presents a validation of the numerical model against vertical jet erosion tests from the literature and conducted at HR Wallingford. The results of the model show good agreement with the experimental tests, with the numerical model predicting the scour hole depth and extent with good accuracy. The paper then presents a validation of the model’s ability to reproduce deposition which is evaluated through a comparison with settling velocity data and empirical formulations found in literature, again with the model showing good agreement. Finally, the model is applied to a prototype cable burial problem using a commercially available controlled flow jet excavator. The study found that the use of water jets can be effective (subject to confirmation of the time-scale required for real operations) for performing backfill operations but that the effectiveness is closely related to the type of sediment and selection of an appropriate jet discharge. As a result, in order for the water jet method to be effective for backfill, there is a requirement for a good description of the variation in sediment type along the trench and a requirement for the jet discharge to be varied as different sediment types are encountered.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-14
Author(s):  
Tamara de Riese ◽  
Paul D. Bons ◽  
Enrique Gomez-Rivas ◽  
Till Sachau

Crustal-scale fluid flow can be regarded as a bimodal transport mechanism. At low hydraulic head gradients, fluid flow through rock porosity is slow and can be described as diffusional. Structures such as hydraulic breccias and hydrothermal veins both form when fluid velocities and pressures are high, which can be achieved by localized fluid transport in space and time, via hydrofractures. Hydrofracture propagation and simultaneous fluid flow can be regarded as a “ballistic” transport mechanism, which is activated when transport by diffusion alone is insufficient to release the local fluid overpressure. The activation of a ballistic system locally reduces the driving force, through allowing the escape of fluid. We use a numerical model to investigate the properties of the two transport modes in general and the transition between them in particular. We developed a numerical model in order to study patterns that result from bimodal transport. When hydrofractures are activated due to low permeability relative to fluid flux, many hydrofractures form that do not extend through the whole system. These abundant hydrofractures follow a power-law size distribution. A Hurst factor of ~0.9 indicates that the system self-organizes. The abundant small-scale hydrofractures organize the formation of large-scale hydrofractures that ascend through the whole system and drain fluids in large bursts. As the relative contribution of porous flow increases, escaping fluid bursts become less frequent, but more regular in time and larger in volume. We propose that metamorphic rocks with abundant veins, such as in the Kodiak accretionary prism (Alaska) and Otago schists (New Zealand), represent regions with abundant hydrofractures near the fluid source, while hydrothermal breccias are formed by the large fluid bursts that can ascend the crust to shallower levels.


Author(s):  
D. C. Hong ◽  
T. B. Ha ◽  
K. H. Song

The added resistance of a ship was calculated using Maruo’s formula [1] involving the three-dimensional Kochin function obtained using the source and normal doublet distribution over the wetted surface of the ship. The density of the doublet distribution was obtained as the solution of the three-dimensional frequency-domain forward-speed Green integral equation containing the exact line integral along the waterline. Numerical results of the Wigley ship models II and III in head seas, obtained by making use of the inner-collocation 9-node second-order boundary element method have been compared with the experimental results reported by Journée [2]. The forward-speed hydrodynamic coefficients of the Wigley models have shown no irregular-frequencylike behavior. The steady disturbance potential due to the constant forward speed of the ship has also been calculated using the Green integral equation associated with the steady forward-speed free-surface Green function since the so-called mj-terms [3] appearing in the body boundary conditions contain the first and second derivatives of the steady potential over the wetted surface of the ship. However, the free-surface boundary condition was kept linear in the present study. The added resistances of the Wigley II and III models in head seas obtained using Maruo’s formula showing acceptable comparison with experimental results, have been presented. The added resistances in following seas obtained using Maruo’s formula have also been presented.


2017 ◽  
Vol 11 (01) ◽  
pp. 1740006 ◽  
Author(s):  
Changbo Jiang ◽  
Xiaojian Liu ◽  
Yu Yao ◽  
Bin Deng ◽  
Jie Chen

Seawall is a most commonly used structure in coastal areas to protect the landscape and coastal facilities. The studies of interactions between the tsunami-like solitary waves and the seawalls are relatively rare in the literature. In this study, a three-dimensional numerical model based on OpenFOAM® was developed to investigate the tsunami-like solitary waves propagating over a rectangular seawall. The Navier–Stokes equations for two-phase incompressible flow, combining with methods of [Formula: see text] for turbulence closure and Volume of Fluid (VOF) for tracking the free surface, were solved. Laboratory experiments were performed to measure some of the hydrodynamic feature associated with solitary waves. The model was then validated by the laboratory data, and good agreements were found for free surface, velocity and dynamic pressure around the seawall. Finally, a series of numerical experiments were conducted to analyze the evolution of both wave and flow fields, the overtopping discharge as well as wave pressure (force) around the seawall, special attention is given to the effects of seawall crest width. Our findings will help to improve the understanding in the occurrences of tsunami-induced damages in the vicinity of seawall such as wave impact and local scouring.


1996 ◽  
Vol 306 ◽  
pp. 293-323 ◽  
Author(s):  
Vadim Borue ◽  
Steven A. Orszag

High-resolution numerical simulations (with up to 2563 modes) are performed for three-dimensional flow driven by the large-scale constant force fy = F cos(x) in a periodic box of size L = 2π (Kolmogorov flow). High Reynolds number is attained by solving the Navier-Stokes equations with hyperviscosity (-1)h+1Δh (h = 8). It is shown that the mean velocity profile of Kolmogorov flow is nearly independent of Reynolds number and has the ‘laminar’ form vy = V cos(x) with a nearly constant eddy viscosity. Nevertheless, the flow is highly turbulent and intermittent even at large scales. The turbulent intensities, energy dissipation rate and various terms in the energy balance equation have the simple coordinate dependence a + b cos(2x) (with a, b constants). This makes Kolmogorov flow a good model to explore the applicability of turbulence transport approximations in open time-dependent flows. It turns out that the standard expression for effective (eddy) viscosity used in K-[Escr ] transport models overpredicts the effective viscosity in regions of high shear rate and should be modified to account for the non-equilibrium character of the flow. Also at large scales the flow is anisotropic but for large Reynolds number the flow is isotropic at small scales. The important problem of local isotropy is systematically studied by measuring longitudinal and transverse components of the energy spectra and crosscorrelation spectra of velocities and velocity-pressure-gradient spectra. Cross-spectra which should vanish in the case of isotropic turbulence decay only algebraically but somewhat faster than corresponding isotropic correlations. It is verified that the pressure plays a crucial role in making the flow locally isotropic. It is demonstrated that anisotropic large-scale flow may be considered locally isotropic at scales which are approximately ten times smaller than the scale of the flow.


Sign in / Sign up

Export Citation Format

Share Document