scholarly journals Detailed Simulation of Complex Hydraulic Problems with Macroscopic and Mesoscopic Mathematical Methods

2013 ◽  
Vol 2013 ◽  
pp. 1-14 ◽  
Author(s):  
Chiara Biscarini ◽  
Silvia Di Francesco ◽  
Fernando Nardi ◽  
Piergiorgio Manciola

The numerical simulation of fast-moving fronts originating from dam or levee breaches is a challenging task for small scale engineering projects. In this work, the use of fully three-dimensional Navier-Stokes (NS) equations and lattice Boltzmann method (LBM) is proposed for testing the validity of, respectively, macroscopic and mesoscopic mathematical models. Macroscopic simulations are performed employing an open-source computational fluid dynamics (CFD) code that solves the NS combined with the volume of fluid (VOF) multiphase method to represent free-surface flows. The mesoscopic model is a front-tracking experimental variant of the LBM. In the proposed LBM the air-gas interface is represented as a surface with zero thickness that handles the passage of the density field from the light to the dense phase and vice versa. A single set of LBM equations represents the liquid phase, while the free surface is characterized by an additional variable, the liquid volume fraction. Case studies show advantages and disadvantages of the proposed LBM and NS with specific regard to the computational efficiency and accuracy in dealing with the simulation of flows through complex geometries. In particular, the validation of the model application is developed by simulating the flow propagating through a synthetic urban setting and comparing results with analytical and experimental laboratory measurements.

2019 ◽  
Vol 21 (2) ◽  
pp. 248-262 ◽  
Author(s):  
David Markt ◽  
Ashish Pathak ◽  
Mehdi Raessi ◽  
Seong-Young Lee ◽  
Roberto Torelli

This article uniquely characterizes the secondary droplets formed during the impingement of a train of ethanol drops, using three-dimensional direct numerical simulations performed under conditions studied experimentally by Yarin and Weiss. Our numerical results have been previously validated against experimental data demonstrating the ability to accurately capture the splashing dynamics. In this work, the predictive ability of the model is leveraged to gain further insight into secondary droplet formation. We present a robust post-processing algorithm, which scrutinizes the liquid volume fraction field in the volume-of-fluid method and quantifies the number, volume and velocity of secondary droplets. The high-resolution computational simulations enable secondary droplet characterization within close proximity of the impingement point at small length and time scales, which is extremely challenging to achieve experimentally. By studying the temporal evolution of secondary droplet formation, direct connections are made between liquid structures seen in the simulation and the instantaneous distribution of secondary droplets, leading to detailed insight into the instability-driven breakup process of lamellae. Time-averaged secondary droplet characteristics are also studied to describe the global distribution of secondary droplets. Such analysis is vital to understanding fuel drop impingement in direct injection engines, facilitating the development of highly accurate spray–wall interaction models for use in Lagrangian solvers.


2020 ◽  
Vol 10 (4) ◽  
pp. 1369 ◽  
Author(s):  
Yogeshwar Nath Mishra ◽  
Timo Tscharntke ◽  
Elias Kristensson ◽  
Edouard Berrocal

Structured laser illumination planar imaging (SLIPI)-based techniques have been employed during the past decade for addressing multiple light scattering issues in spray imaging. In this article, SLIPI droplet sizing based on the intensity ratio of laser-induced fluorescence (LIF) over Mie scattering (SLIPI-LIF/Mie) and SLIPI-Scan for extinction-coefficient (µe) mapping are applied simultaneously. In addition, phase Doppler anemometry (PDA) and numerical calculations based on the Lorenz–Mie theory are also employed in order to extract the droplets Sauter mean diameter (SMD), the droplets number density (N), and the liquid volume fraction (LVF) in a steady asymmetric hollow cone water spray. The SLIPI-LIF/Mie ratio is converted to droplets SMD by means of a calibration procedure based on PDA measurements. The droplet SMD for the investigated spray varies from 20 µm to 60 µm, the N values range from 5 to 60 droplets per mm3, and the LVF varies between 0.05 × 10−4 and 5.5 × 10−4 within the probed region of the spray. To generate a series of two-dimensional images at different planes, the spray scanning procedure is operated in a “bread slicing” manner by moving the spray perpendicularly to the light sheet axis. From the resulting series of images, the procedure described here shows the possibility of obtaining three-dimensional reconstructions of each scalar quantity, allowing a more complete characterization of droplet clouds forming the spray region.


2010 ◽  
Vol 97-101 ◽  
pp. 2736-2743
Author(s):  
Shi Xiong Ren ◽  
Sha Sha Dang ◽  
Tao Lu ◽  
Kui Sheng Wang

Three-dimensional models of heat transfer have been established and numerically solved using a commercial software package, Fluent, in order to obtain distributions of temperature, velocity, pressure, and liquid volume fraction of the polymer. The influences of the boundary conditions on the phase change of the polymer and the temperature distribution in the die have been evaluated. The results show that the temperature of the region close to the pelletizing surface is relatively low due to the cooling effect of the cool water, while the temperature deeper inside the die is higher, with a lower temperature gradient, as a result of the heating effect of the hot thermal oil and the polymer. A solidification phase change of the polymer occurs near the polymer outlet due to heat loss from the polymer to the water, while deeper inside the hole the polymer remains fluid without solidification, due to heating by the thermal oil. Numerical simulation provides a reliable method to optimize the design of the die, the choice of metallic material for the die, and the operating conditions of the polymer pelletizing under water.


2014 ◽  
Vol 17 (1) ◽  
pp. 1-6 ◽  
Author(s):  
Silvia Di Francesco ◽  
Chiara Biscarini ◽  
Piergiorgio Manciola

Three-dimensional (3D) hydraulic modelling of rapidly varying surface flows is a challenging task for practical engineering applications. One example is represented by the fast-moving fronts originating from dam breaches that proceed downstream through artificial channels. In this work, a fully 3D lattice Boltzmann method (LBM) is tested. The numerical model is a front-tracking variant of the LBM, being the free surface tracked through the liquid volume fraction. Model performances are evaluated simulating the effect of dam-break flows on synthetic settings schematically represented by an artificial domain and comparing results with analytical data and experimental laboratory measurements. Obtained results are promising for the use of LBM for practical engineering applications.


2005 ◽  
Vol 40 ◽  
pp. 179-184 ◽  
Author(s):  
C. Bock ◽  
H. Eicken

AbstractThe microstructural evolution of brine inclusions in granular and columnar sea ice has been investigated through magnetic resonance imaging (MRI) for temperatures between –28 and –3˚C. Thin-section and salinity measurements were completed on core samples obtained from winter sea ice near Barrow, Alaska, USA. Subsamples of granular (2–5cm depth in core) and columnar sea ice (20–23 cm depth) were investigated with morphological spin-echo and diffusion-weighted imaging in a Bruker 4.7T MRI system operating at field gradients of 200 mTm–1 at temperatures of approximately –28, –15, –6 and –3˚C. Average linear pore dimensions range from 0.2 to 1 mm and increase with bulk liquid volume fraction as temperatures rise from –15 to –3˚C. Granular ice pores are significantly larger than columnar ice pores and exhibit a higher degree of connectivity. No evidence is found of strongly non-linear increases in pore connectivity based on the MRI data. This might be explained by shortcomings in resolution, sensitivity and lack of truly three-dimensional data, differences between laboratory and field conditions or the absence of a percolation transition. Pore connectivity increases between –6 and –3˚C. Pore-number densities average at 1.4±1.2mm–2. The pore-number density distribution as a function of cross-sectional area conforms with power-law and lognormal distributions previously identified, although significant variations occur as a function of ice type and temperature. At low temperatures (< –26˚C), pore sizes were estimated from 1H self-diffusivity measurements, with self-diffusivity lower by up to an order of magnitude than in the free liquid. Analysis of diffusional length scales suggests characteristic pore dimensions of <1 μm at < –26˚C.


Author(s):  
Luis San Andrés ◽  
Xueliang Lu

Wet gas compression systems and multiphase pumps are enabling technologies for the deep sea oil and gas industry. This extreme environment determines both machine types have to handle mixtures with a gas in liquid volume fraction (GVF) varying over a wide range (0 to 1). The gas (or liquid) content affects the system pumping (or compression) efficiency and reliability, and places a penalty in leakage and rotordynamic performance in secondary flow components, namely seals. In 2015, tests were conducted with a short length smooth surface annular seal (L/D = 0.36, radial clearance = 0.127 mm) operating with an oil in air mixture whose liquid volume fraction (LVF) varied to 4%. The test results with a stationary journal show the dramatic effect of a few droplets of liquid on the production of large damping coefficients. This paper presents further measurements and predictions of leakage, drag power, and rotordynamic force coefficients conducted with the same test seal and a rotating journal. The seal is supplied with a mixture (air in ISO VG 10 oil), varying from a pure liquid to an inlet GVF = 0.9 (mostly gas), a typical range in multiphase pumps. For operation with a supply pressure (Ps) up to 3.5 bar (a), discharge pressure (Pa) = 1 bar (a), and various shaft speed (Ω) to 3.5 krpm (ΩR = 23.3 m/s), the flow is laminar with either a pure oil or a mixture. As the inlet GVF increases to 0.9 the mass flow rate and drag power decrease monotonically by 25% and 85% when compared to the pure liquid case, respectively. For operation with Ps = 2.5 bar (a) and Ω to 3.5 krpm, dynamic load tests with frequency 0 < ω < 110 Hz are conducted to procure rotordynamic force coefficients. A direct stiffness (K), an added mass (M) and a viscous damping coefficient (C) represent well the seal lubricated with a pure oil. For tests with a mixture (GVFmax = 0.9), the seal dynamic complex stiffness Re(H) increases with whirl frequency (ω); that is, Re(H) differs from (K-ω2M). Both the seal cross coupled stiffnesses (KXY and −KYX) and direct damping coefficients (CXX and CYY) decrease by approximately 75% as the inlet GVF increases to 0.9. The finding reveals that the frequency at which the effective damping coefficient (CXXeff = CXX-KXY/ω) changes from negative to positive (i.e., a crossover frequency) drops from 50% of the rotor speed (ω = 1/2 Ω) for a seal with pure oil to a lesser magnitude for operation with a mixture. Predictions for leakage and drag power based on a homogeneous bulk flow model match well the test data for operation with inlet GVF up to 0.9. Predicted force coefficients correlate well with the test data for mixtures with GVF up to 0.6. For a mixture with a larger GVF, the model under predicts the direct damping coefficients by as much as 40%. The tests also reveal the appearance of a self-excited seal motion with a low frequency; its amplitude and broad band frequency (centered at around ∼12 Hz) persist and increase as the gas content in the mixture increase. The test results show that an accurate quantification of wet seals dynamic force response is necessary for the design of robust subsea flow assurance systems.


Sign in / Sign up

Export Citation Format

Share Document