The Oscillating Water Column (OWC) wave energy converter has been shown to have high potential, thus rendering extensive development in recent years. In order to further accelerate its development, ...highly accurate yet computationally efficient tools are necessary particularly when studying the interaction of multiple OWC devices. This paper proposes a new framework for fixed OWC devices with an orifice, that uses the input from a high fidelity non-linear numerical model to improve the accuracy of a low fidelity linear numerical model keeping computational costs low. This is done by accounting for the non-linearities in the pressure-flow of an orifice in the input to the linear numerical model. Experimental data is used to validate the framework, thus providing an accurate and computationally efficient linear numerical model, that can be used for the preliminary analysis of fixed OWC devices.
•A new framework combining low-fidelity and high-fidelity numerical methods for fixed OWC systems.•Accounting for non-linearities in pressure-flow in an orifice in a linear BEM model.•Validation with experimental data.•Use of validated open-source codes.
•Cavitating flow around a 3-D Hydrofoil with Wavy Leading Edges (WLE) is investigated using the LES approach in OpenFOAM.•Comparison of cavity features of wavy leading edges (WLE) hydrofoil with ...straight leading edge (SLE) one.•Force coefficients, re-entrant jet and separation are assessed during the cloud cavity evolution around WLE and SLE hydrofoils.•The mechanisms of the laminar separation bubble (LSB) and low-pressure zone behind the WLE hydrofoil are illustrated.•Role of vorticity stretching and vorticity dilatation with the presence of cavity for the WLE and SLE hydrofoils are investigated.
The present study seeks to conduct numerical investigations of the cavitating flow characteristics around a sinusoidal wavy leading edge (WLE) 3-D hydrofoil underlying a NACA 634–021 profile with an aspect ratio of 4.3. Cavitational and non-cavitational characteristics of hydrofoils are numerically examined at a chord-based Reynolds number of 7.2 × 105. The sinusoidal leading edge geometries include two WLE amplitudes of 5% and 25% and two WLE wavelengths of 25% and 50% of the mean chord length. We examined the cavitating flow around the hydrofoils in different cavitation numbers, namely σ = 0.8 and σ = 1.2. The flow over the protuberances of the WLE hydrofoil is considered at varying chord lengths and a constant angle of attack α = 6°, where significant spanwise variations in all flow properties, in contrast to the straight leading edge (SLE) hydrofoil, were observed. Large eddy simulation (LES) and Kunz mass transfer models are employed to simulate the dynamic and unsteady behavior of the cavitating flow. Besides, the compressive volume of fluid (VOF) method is used to track the cavity interface. Simulation is performed under the two-phase flow solver —interPhaseChangeFoam— of the OpenFOAM package. Compared to the SLE hydrofoil, we provided an exhaustive report of the time-averaged and instantaneous fluid dynamic characteristics of the cavitating flow around the sinusoidal leading edge hydrofoil, i.e., pressure, velocity, and vorticity fields, as well as lift and drag coefficients, and turbulent kinetic energy are reported. Furthermore, detailed analyses of the instantaneous cavity leading edge and flow separation treatment, vortical structure of the flow, vorticity stretching and dilatation, details of the spanwise flow, the formation of a low-pressure zone behind the WLE hydrofoil, streamwise velocity fluctuation, and evolution of the cavity dynamics through a complete cycle are reported. Results show that early development of the laminar separation bubble (LSBs) on the suction side of WLE hydrofoil prevents significant flow separation. Furthermore, the WLE cases exhibit a significantly reduced level of unsteady fluctuations in aerodynamic forces at the frequency of periodic vortex shedding.
This paper numerically investigates the effect of concentration gradient on the re-initiation process of H2/O2 detonation propagating from a vertical channel to a horizontal channel. Two different ...types of concentration gradient distribution are considered in this study: from fuel-rich near the upper wall to fuel-lean near the lower wall (gradient > 0), from fuel-rich near the lower wall to fuel-lean near the upper wall (gradient < 0). The results show that the detonation can be re-initiated successfully for the two different types of concentration gradient, but the formation modes of the transverse waves in the re-initiation process are different. When the gradient > 0, the transverse waves are mainly formed after the 2nd reflection. When the gradient < 0, the transverse waves are mainly formed after the 1st reflection. However, whether the concentration gradient > 0 or < 0, the transverse waves are mainly formed on the fuel-lean side. The earlier transverse waves are mainly originated from the transverse waves of the transverse detonation. The detonation speed under fuel-rich condition is much higher as predicted by 1D simulation in the current study, and a local reactant distribution from fuel-lean to fuel-rich is formed by the leading shock in the horizontal direction. Thus, the transverse waves of the transverse detonation on the fuel-lean side can overtake the Mach detonation. Furthermore, the number of the transverse waves decreases with steeper concentration gradient, leading to the cell structure more irregular and the detonation transition from multi-head mode to single-head mode.
•Unsteady sheet/cloud cavitating flows are studied using the compressible cavitation solver developed on OpenFOAM.•Cavitation fluid compressibility effects on transient flow structures and dynamics ...are illustrated.•Effects of density variance on velocity divergence are proposed.•The physics for the compressibility effects on cavitation vortex dynamics are investigated based on VTE budget analysis.
The objective of the paper is to investigate the physics involved in the compressible cavitating flows, with emphasis on the compressibility effects. 3-D numerical simulations were conducted on the open source software platform OpenFOAM, using both the native incompressible cavitation solver interPhaseChangeFOAM and implemented compressible cavitation solver, where the cavitation model and turbulence model are kept the same and differences of the two approaches mainly root in the density variances of pure liquid and pure vapor. Results are presented for the transient sheet/cloud cavitating flows around a Clark-Y hydrofoil fixed at attack of angle α=8° at inlet velocity U=10m/s and cavitation number σ=0.8, where both ensemble averaged statistics and transient characteristics are analyzed. Good agreement can be obtained using both the incompressible and compressible approaches when compared with the experiment data. While it is found that compared with the incompressible approach, the compressible approach can predict the unsteady cavitation evolution and cavity shedding frequency better. With the compressibility effects considered, the time averaged void fraction distribution decreases, and the cavity size (i.e. cavitation area) becomes smaller. The re-entrant flow thickness normalized by local cavity thickness predicted by the compressible approach is larger than that by the incompressible approach, indicating that the compressible approach can predict the re-entrant jet dynamics well. The velocity divergence analysis show that compared with that in incompressible approach, where velocity divergence mainly comes from the mass transfer between phases, in compressible approach, the velocity divergence originates from both the cavitation two-phase fluid compressibility and mass transfer, and the fluid density variance dominates in compressible results. Following, the budget analysis of vorticity transport equation (VTE) show that the vortex stretching term dominates the cavitation vortex dynamics. Compressibility effects will significantly increase the dilatation term and decrease baroclinic term by decreasing the misalignment between density gradient and pressure gradient. Finally, the temperature and density variance in different cavitation structures are presented.
In this work, the performance of a buoyancy-modified turbulence model is shown for simulating wave breaking in a numerical wave flume. Reynolds-Averaged Navier-Stokes (RANS) modelling is performed by ...applying both a k-ω and a k-ω SST turbulence model using the Computational Fluid Dynamics (CFD) toolbox OpenFOAM. In previous work of the authors (Devolder et al., 2017), the observed significant decrease in wave height over the length of the numerical wave flume based on RANS turbulence modelling for the case of propagating waves has been avoided by developing a buoyancy-modified k-ω SST model in which (i) the density is explicitly included in the turbulence transport equations and (ii) a buoyancy term is added to the turbulent kinetic energy (TKE) equation. In this paper, two buoyancy-modified turbulence models are applied for the case of wave breaking simulations: k-ω and k-ω SST. Numerical results of wave breaking under regular waves are validated with experimental data measured in a wave flume by Ting and Kirby (1994). The numerical results show a good agreement with the experimental measurements for the surface elevations, undertow profiles of the horizontal velocity and turbulent kinetic energy profiles. Moreover, the underlying motivations for the concept of a buoyancy-modified turbulence model are demonstrated by the numerical results and confirmed by the experimental observations. Firstly, the buoyancy term forces the solution of the flow field near the free water surface to a laminar solution in case of wave propagation. Secondly in the surf zone where waves break, the buoyancy term goes to zero and a fully turbulent solution of the flow field is calculated. Finally and most importantly, the buoyancy-modified turbulence models significantly reduce the common overestimation of TKE in the flow field.
•Simulations of wave breaking in a numerical wave flume using OpenFOAM®.•Performance of a buoyancy-modified k-ω and k-ω SST turbulence model are evaluated.•The buoyancy-modified turbulence models are validated with experimental results.•Accurate predictions of the turbulent flow field at the wave breaking zone.
Display omitted
•Three node CHyQMOM is successfully implemented in OpenFOAM-7.•The algorithm can capture particle trajectory crossing without delta shocks.•For wall bounded problems results are not ...accurate.•Attention should be paid to prevent over-accumulation.
The modeling of dilute gas-particle flows is challenging due to particle-trajectory crossing (PTC). Lagrangian particle tracking can be used, but requires a large number of parcels resulting in high computational costs. A less costly method is the Eulerian number density function (NDF) approach, based on the Boltzmann equation, often solved in terms of lower-order moments of the NDF. In this context the conditional hyperbolic quadrature method of moments (CHyQMOM) was developed and is here implemented for the first time in the OpenFOAM-7, together with high-order advection schemes and a new operator splitting procedure. The resulting solver is used to simulate different test cases: phase segregation problems, collision-less and weakly-collisional PTC flows, asymmetric and symmetric Taylor-Green vortex flow and a dilute gas-particle riser. Results, validated against analytical solutions and predictions obtained with Lagrangian particle tracking, show that the implemented CHyQMOM can be used to handle highly non-equilibrium flows.
Although OpenFOAM is a widely-used open source computational fluid dynamics (CFD) tool, it is limited to numerical simulations of multi-dimensional reacting/nonreacting flows at relatively-low ...pressures. This is not only because real-fluid models that can evaluate thermophysical properties at high pressures are not available in the thermophysicalModels library of OpenFOAM, but also because the existing mixing model cannot handle various mixing rules of real-fluid models. In the present study, we develop a novel algorithm applicable for a mixture model incorporating various mixing rules in OpenFOAM. Based on the new algorithm, we update the thermophysicalModels library of OpenFOAM 6.0 by implementing a set of real-fluid models such as the Soave-Redlich-Kwong/Peng-Robinson equation of state, Chung's model for dynamic viscosity and thermal conductivity, mixture averaged model for mass diffusivity using Takahashi's correction for binary diffusion coefficients at high pressure. The new library is validated against experimental data and is further assessed for compressible reacting flows by performing two-dimensional numerical simulations of axisymmetric laminar non-premixed counterflow flames and one-dimensional numerical simulations of premixed CH4/air flames at high pressures. The developed library can be used for any reacting flow solvers in OpenFOAM 6.0 that adopt a set of implemented real-fluid models.
Program title: Real-fluid thermophysicalModels
CPC Library link to program files:https://doi.org/10.17632/n8zb2wpjp6.1
Developer's repository link:https://github.com/danhnam11/realFluidThermophysicalModels-6
Licensing provisions: GPLv3
Programming language: C++
Nature of problem: The mixing rules required for evaluating the real-fluid based thermophysical properties of mixtures are relatively complicated and different from each other, while the existing thermophysicalModels library offers only a mixing rule based on the mass fraction weighted average of each species, for which overloading operator functions (“*” and “+=”) are utilized to update parameters for a mixture. However, this approach is not applicable for complicated mixing rules in the real-fluid models.
Solution method: To build a mixture class containing different mixing rules in different models such as the equation of states and thermodynamics/transport properties in OpenFOAM, we adopt void-type functions (parameter updating functions) to update parameters of a mixture instead of the original overloading operator functions in OpenFOAM. The parameters in the mixing rules that depend on pressure and temperature are handled by decomposition such that they can be updated in the mixture class.
This paper presents a new methodology allowing the discretisation of phase specific transport equations within the Volume-of-Fluid interface capturing method framework. The method uses a sharp ...interface algorithm to compute the transport of species concentration. The interface reconstruction and advection is provided by the geometric advection scheme isoAdvector1 implemented in the OpenFOAM® library. When discretising the transport equation for species, lack of consistency with the free surface advection scheme can lead to numerical errors, causing conservation or boundedness issues. This work addresses the issue of consistency in convective transport of species and is divided in two parts. First, a new interpolation procedure is used to compute face values from cell-centered values. Then, the diffusive operator of the transport equation is corrected. Finally, a set of test cases are presented to validate the transport equation's consistency with the free surface advection. Species transport across the interface is not part of the scope of this article, however, this methodology can further be used to study mass transfer at the gas-liquid interface using additional mass source terms that are not discussed here.
•A new methodology to discretize a transport equation for species concentration has been described.•This methodology is consistent with the geometric Volume-of-Fluid method isoAdvector which is available in OpenFOAM.•The developments made in this work avoid artificial mass transfer between the phases.•Diffusive fluxes were also corrected to avoid non-physical fluxes using the species concentration inside the phase.
The Lead-bismuth cooled Fast Reactors (LFR) utilizes wrapped hexagonal fuel assemblies arranged in a honeycomb pattern. The inter-wrapper flow (IWF) between adjacent assemblies has a significant ...influence on the core thermal-hydraulics. In this paper, the LFR core assembly analysis model, inter-wrapper flow and heat transfer model, and multi-scale coupling analysis model were developed and were implemented into the previously self-developed subchannel-level three-dimensional thermal hydraulic analysis code CorTAF-LBE for whole core of LFR, which is based on the OpenFOAM computational fluid dynamics simulation platform. The inter-wrapper flow and heat transfer model was validated against the KALLA-IWF experiment. Whole core steady state operating condition and low flow rate condition simulations were performed with reference to the 19-rod bundle geometry from KALLA laboratory and typical LFR core assembly arrangement. The whole core steady state distribution laws of key thermal-hydraulic parameters were obtained, and the influence of low flow rate operation on safety parameters such as cladding and pellet temperatures was analyzed. Under low flow rate operating conditions, the IWF heat transfer was accounted for 4.87% of the core thermal power, which was 315.6 kW higher than that under normal operating conditions, playing an important role in core heat removal and safe operation. This work provides important insights into LFR core design and the assessment of the role of IWF under low flow rate conditions.
•The Multi-scale coupling analysis model is established and implanted into the CorTAF-LBE code.•The steady-state thermal-hydraulic parameter distribution at whole core scales is analyzed.•The effect of inter-wrapper flow on core heat transfer at low flow rate condition is evaluated.
Liquid Metal Batteries (LMBs) are a promising concept for cheap electrical energy storage at grid level. These are built as a stable density stratification of three liquid layers, with two liquid ...metals separated by a molten salt. In order to ensure a safe and efficient operation, the understanding of transport phenomena in LMBs is essential. With this motivation we study thermal convection induced by internal heat generation. We consider the electrochemical nature of the cell in order to define the heat balance and the operating parameters. Moreover we develop a simple 1D heat conduction model as well as a fully 3D thermo-fluid dynamics model. The latter is implemented in the CFD library OpenFOAM, extending the volume of fluid solver, and validated against a pseudo-spectral code. Both models are used to study a rectangular 10×10 cm² Li||Bi LMB cell at three different states of charge.