scholarly journals Subgrid Liquid Flux and interface modelling for LES of Atomization

Author(s):  
Aqeel Ahmed ◽  
Javier Anez ◽  
Stefano Puggelli ◽  
Julien Reveillon ◽  
Jorge César Brändle de Motta ◽  
...  

Traditional Discrete Particle Methods (DPM) such as the Euler-Lagrange approaches for modelling atomization, even if widely used in technical literature, are not suitable in the near injector region. Indeed, the first step of atomization process is to separate the continuous liquid phase in a set of individual liquid parcels, the so-called primary break-up. Describing two-phase flow by DPM is to define a carrier phase and a discrete phase, hence they cannot be used for primary breakup. On the other hand, full scale simulations (direct simulation of the dynamic DNS, and interface capturing method ICM) are powerful numerical tools to study atomization, however, computational costs limit their application to academic cases for understanding and complementing partial experimental data. In an industrial environment, models that are computationally cheap and still accurate enough are required to meet new challenges of fuel consumption and pollutant reduction. Application of DNS-ICM methods without fairly enough resolution to solve all length scales are currently used for industrial purpose. Nevertheless, effects of unresolved scales are generally cast aside. The Euler-Lagrange Spray Atomization model family (namely, ELSA, also call, Σ−𝑌 or Ω−𝑌) developed by Vallet and Borghi pioneering work [1], and [2], at the contrary aims to model those unresolved terms. This approach is actually complementary to DNS-ICM method since the importance of the unresolved term depends directly on mesh resolution. For full interface resolution the unclosed terms are negligible, except in the far-field spray when the unresolved terms become dominant. Depending on the complexity of the flow and the available computational resources, a Large Eddy Simulation (LES) formalism could be employed as modelling approach. This work focus on the two main terms that drive these different modelling approaches namely the subgrid turbulent liquid flux and the resolved interface. Thanks to the open source library OpenFoam® this work is an attempt to review and to release an adapted modelling strategy depending on the available mesh resolution. For validation, these solvers are tested against realistic experimental data to see the overall effect of each model proposal.DOI: http://dx.doi.org/10.4995/ILASS2017.2017.4694

Author(s):  
J. C. Jepsen ◽  
J. L. Ralph

The object of this study was to obtain data on the radial variation of gas- and liquid-phase mass flux profiles in two-phase upflow in vertical pipelines. Experimental data were obtained on the radial gas-liquid flux, impact pressure, and linear liquid-phase velocity profiles for superficial gas- and liquid-phase velocities ranging from 20 to 125 ft/s and from 5 to 15 ft/s, respectively. Studies were made on 1-in, 4-in, and 8·4-in vertical pipelines and in a 1-in i.d. by 8·4-in o.d. vertical annular flow pipeline. Gas-liquid systems studied were air-water, air-aqueous glycerol, and air-tetrabromoethane. In Part 2 time-averaged radial liquid hold-up and linear-phase velocities were estimated from radial mass flux and impact pressure data. Estimates were also made on the magnitude, frequency, and velocity of the flow disturbances. Empirical correlations for mean liquid hold-up and pressure drop were developed from experimental data.


Author(s):  
C.-L. Ng ◽  
K. A. Sallam

The deformation of laminar liquid jets in gaseous crossflow before the onset of primary breakup is studied motivated by its application to fuel injection in jet afterburners and agricultural sprays, among others. Three crossflow Weber numbers that represent three different liquid jet breakup regimes; column, bag, and shear breakup regimes, were studied at large liquid/gas density ratios and small Ohnesorge numbers. In each case the liquid jet was simulated from the jet exit and ended before the location where the experimental data indicated the onset of breakup. The results show that in column and bag breakup, the reduced pressures along the sides of the jet cause the liquid to move to the sides of the jet and enhance the jet deformation. In shear breakup, the flattened upwind surface pushes the liquid towards the two sides of the jet and causing the gaseous crossflow to separate near the edges of the liquid jet thus preventing further deformation before the onset of breakup. It was also found out that in shear breakup regime, the liquid phase velocity inside the liquid jet was large enough to cause onset of ligament formation along the jet side, which was not the case in the column and bag breakup regimes. In bag breakup, downwind surface waves were observed to grow along the sides of the liquid jet triggered a complimentary experimental study that confirmed the existence of those waves for the first time.


Fluids ◽  
2021 ◽  
Vol 6 (2) ◽  
pp. 80
Author(s):  
Yuria Okagaki ◽  
Taisuke Yonomoto ◽  
Masahiro Ishigaki ◽  
Yoshiyasu Hirose

Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in the convection term of the volume fraction equation. Thus, in this study, we focused on an interface compression (IC) method for such a VOF approach; this technique prevents numerical diffusion issues and maintains boundedness and conservation with negative diffusion. First, on a sufficiently high mesh resolution and without the IC method, the validation process was considered by comparing the amplitude growth of the interfacial wave between a two-dimensional gas sheet and a quiescent liquid using the linear theory. The disturbance growth rates were consistent with the linear theory, and the validation process was considered appropriate. Then, this validation process confirmed the effects of the IC method on numerical diffusion, and we derived the optimum value of the IC coefficient, which is the parameter that controls the numerical diffusion.


Author(s):  
Stephan Uhkoetter ◽  
Stefan aus der Wiesche ◽  
Michael Kursch ◽  
Christian Beck

The traditional method for hydrodynamic journal bearing analysis usually applies the lubrication theory based on the Reynolds equation and suitable empirical modifications to cover turbulence, heat transfer, and cavitation. In cases of complex bearing geometries for steam and heavy-duty gas turbines this approach has its obvious restrictions in regard to detail flow recirculation, mixing, mass balance, and filling level phenomena. These limitations could be circumvented by applying a computational fluid dynamics (CFD) approach resting closer to the fundamental physical laws. The present contribution reports about the state of the art of such a fully three-dimensional multiphase-flow CFD approach including cavitation and air entrainment for high-speed turbo-machinery journal bearings. It has been developed and validated using experimental data. Due to the high ambient shear rates in bearings, the multiphase-flow model for journal bearings requires substantial modifications in comparison to common two-phase flow simulations. Based on experimental data, it is found, that particular cavitation phenomena are essential for the understanding of steam and heavy-duty type gas turbine journal bearings.


2011 ◽  
Vol 133 (5) ◽  
Author(s):  
F. P. P. Tan ◽  
N. B. Wood ◽  
G. Tabor ◽  
X. Y. Xu

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y+ less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y+ less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.


Processes ◽  
2021 ◽  
Vol 9 (3) ◽  
pp. 413
Author(s):  
Sandra Lopez-Zamora ◽  
Jeonghoon Kong ◽  
Salvador Escobedo ◽  
Hugo de Lasa

The prediction of phase equilibria for hydrocarbon/water blends in separators, is a subject of considerable importance for chemical processes. Despite its relevance, there are still pending questions. Among them, is the prediction of the correct number of phases. While a stability analysis using the Gibbs Free Energy of mixing and the NRTL model, provide a good understanding with calculation issues, when using HYSYS V9 and Aspen Plus V9 software, this shows that significant phase equilibrium uncertainties still exist. To clarify these matters, n-octane and water blends, are good surrogates of naphtha/water mixtures. Runs were developed in a CREC vapor–liquid (VL_ Cell operated with octane–water mixtures under dynamic conditions and used to establish the two-phase (liquid–vapor) and three phase (liquid–liquid–vapor) domains. Results obtained demonstrate that the two phase region (full solubility in the liquid phase) of n-octane in water at 100 °C is in the 10-4 mol fraction range, and it is larger than the 10-5 mol fraction predicted by Aspen Plus and the 10-7 mol fraction reported in the technical literature. Furthermore, and to provide an effective and accurate method for predicting the number of phases, a machine learning (ML) technique was implemented and successfully demonstrated, in the present study.


Sign in / Sign up

Export Citation Format

Share Document