Boiling flow simulations on adaptive octree grids

Boiling flow simulations are conducted on adaptive octree grids. A phase change model consistent with the mixture formulation, in conjunction with the Volume-of-Fluid (VOF) model, is used to track the liquid–vapor interface. Test cases including Rayleigh Taylor instability and bubble growth in a uniform superheat are conducted to validate the phase change model on adaptive grids. The validated model is then used to conduct film boiling simulations on both two-dimensional and three-dimensional adaptive grids. The average wall Nusselt number agrees well with the widely accepted correlations of Berenson (1961) and Klimenko (1981) and Klimenko and Shelepen (1982) for film boiling on a horizontal surface. For the test cases presented, the efficiency of the adaptive technique, as measured by the adaptive mesh refinement (AMR) efficiency, is mostly in the range of 50–80%. Although this efficiency is a function of the nature and dimensionality of the problem, this range of efficiency is comparable to those obtained in the simulations of primary jet atomization conducted by Fuster et al. (2009). This work opens the prospect of conducting more realistic (three-dimensional) multi-modal boiling flow simulations, and problems of similar complexity, in an efficient manner.     

入水后爆炸问题的数值模拟技术

应用各种数值计算方法及计算处理技术,编制程序实现了对弹体入水后爆炸问题全过程的数值模拟.其中弹体与水,爆轰产物气体与弹壳之间的相互作用通过流固耦合技术来描述;水面与空气、爆轰产物与水、空气之间的相互作用采用VOF方法(Volume-of-Fluid)来描述;采用了刚体-柔体转换、单元失效删除等计算技术以更高效、更好地模...     

A two-phase solver for complex fluids: Studies of the Weissenberg effect

In this work a new two-phase solver is presented and described, with a particular interest in the solution of highly elastic flows of viscoelastic fluids. The proposed code is based on a combination of classical Volume-of-Fluid and Continuum Surface Force methods, along with a generic kernel-conformation tensor transformation to represent the rheological characteristics of the (multi)-fluid phases. Benchmark test problems are solved in order to assess the numerical accuracy of distinct levels of physical complexities, such as the interface representation, the influence of advection schemes, the influence of surface tension and the role of fluid rheology. In order to demonstrate the new features and capabilities of the solver in simulating of complex fluids in transient regime, we have performed a set of simulations for the problem of a rotating rod inserted into a container with a viscoelastic fluid, known as the Weissenberg or Rod-Climbing effect. Firstly, our results are compared with numerical and experimental data from the literature for low angular velocities. Secondly, we have presented results obtained for high angular velocities (high elasticity) using the Oldroyd-B model which displayed very elevated climbing heights. Furthermore, above a critical value for the angular velocity, it was observed the onset of elastic instabilities driven by the combination of elastic stresses, interfacial curvature and secondary flows, that to the authors best knowledge, were not yet reported in the literature.     

Modelling far-field entrainment in compressible flows

A procedure that generates steady-states with accurate far-field entrainment is described here. It can be applied to any existing code originally designed for unsteady simulations. Such steady-states are obtained with physical-time damping, introduced through a dual time stepping methodology to effectively eliminate all transient fluctuations without affecting spatial resolution. Far-field entrainment is guaranteed by using properly selected local boundary conditions. Hence, reference profiles with accurate far-field entrainment are generated from the very same code that will eventually use them in unsteady simulations. This new procedure is tested on an unsteady code designed for compressible planar mixing-layer simulations at arbitrary Mach numbers. It is numerically stable for a wide range of Mach numbers and velocity ratios.     

Numerical simulation of ligament-growth on a spinning wheel

Liquid ligaments can grow from perturbations in liquid film spread on a spinning wheel due to the centrifugal force acting on the film. Typically, the growth is strongly influenced by the surface tension on the evolving liquid–air interface. This phenomenon, frequently exploited in industry for the production of fibers, was investigated numerically and the Volume-of-Fluid (VOF) methodology was used to model the interface. The freely available Gerris code with adaptive mesh refinement (AMR) was used to achieve the fine resolution of the computational grid required at the evolving liquid–air interface. The results were compared with the experimental data captured by a high-speed camera. The influence of the process operating variables on the ligament growth is also presented.     

Multiscale simulations and experiments on water jet atomization

Numerical simulation of primary atomization at high Reynolds number is still a challenging problem. In this work a multiscale approach for the numerical simulation of liquid jet primary atomization is applied, using an Eulerian-Lagrangian coupling. In this approach, an Eulerian volume of fluid (VOF) method, where the Reynolds stresses are closed by a Reynolds stress model is applied to model the global spreading of the liquid jet. The formation of the micro-scale droplets, which are usually smaller than the grid spacing in the computational domain, is modelled by an energy-based sub-grid model. Where the disruptive forces (turbulence and surface pressure) of turbulent eddies near the surface of the jet overcome the capillary forces, droplets are released with the local properties of the corresponding eddies. The dynamics of the generated droplets are modelled using Lagrangian particle tracking (LPT). A numerical coupling between the Eulerian and Lagrangian frames is then established via source terms in conservation equations. As a follow-up study to our investigation in Saeedipour et al. (2016a), the present paper aims at modelling drop formation from liquid jets at high Reynolds numbers in the atomization regime and validating the simulation results against in-house experiments. For this purpose, phase-Doppler anemometry (PDA) was used to measure the droplet size and velocity distributions in sprays produced by water jet breakup at different Reynolds numbers in the atomization regime. The spray properties, such as droplet size spectra, local and global Sauter-mean drop sizes and velocity distributions obtained from the simulations are compared with experiment at various locations with very good agreement.     

Coaxial atomization of a round liquid jet in a high speed gas stream: A phenomenological study

Coaxial injectors have proven to be advantageous for the injection, atomization and mixing of propellants in cryogenic H₂/O₂ rocket engines. Thereby, a round liquid oxygen jet is atomized by a fast, coaxial gaseous hydrogen jet. This article summarizes phenomenological studies of coaxial spray generation under a broad variation of influencing parameters including injector design, inflow, and fluid conditions. The experimental investigations, performed using spark light photography and high speed cinematography in a shadow graph setup as main diagnostic means, illuminate the most important processes leading to atomization. These are identified as turbulence in the liquid jet, surface instability, surface wave growth and droplet detachment. Numerical simulations including free surface flow phenomena are a further diagnostic tool to elucidate some atomization particulars. The results of the study are of general importance in the field of liquid atomization.     

Parametric study on the fuel film breakup of a cold start PFI engine

In order to provide more insight on improving the cold start fuel atomization for reducing unburned hydrocarbon emissions, the liquid fuel film breakup phenomenon in the intake valve/port region was investigated in depth for port-fuel-injected engines. Experiments were conducted using high-speed high-resolution imaging techniques to visualize the liquid film atomization and airflow patterns in an axisymmetric steady flow apparatus. The impact of valve/port seat geometry, surface roughness, and fuel properties on airflow separation and fuel film breakup were determined through a parametric study. CFD simulations were also performed with FLUENT to help understand the airflow behavior inside the intake port and valve gap region and its potential impact on fuel film atomization.     

Numerical simulation of primary break-up and atomization: DNS and modelling study

This work deals with numerical simulations of atomization with high Weber and Reynolds values. A special attention has been devoted to the modelling of primary break-up. Due to progress of direct numerical simulations (DNS) of two phase flows it is now possible to simulate the primary break-up of a Diesel spray [Menard, T., Tanguy, S., Berlemont, A., 2007. Coupling level set/VOF/ghost fluid methods: validation and application to 3D simulation of the primary break-up of a liquid jet. Int. J. Multiphase Flow 33 (5), 510–524]. The present formulation of the so-called ELSA (Eulerian–Lagrangian Spray Atomization model) [Vallet, A., Borghi, R., 1999. Modélisation Eulerienne de L’atomisation d’un Jet Liquide. C. R. Acad. Sci. Paris Sér. II b 327, 1015–1020] for atomization is presented and evaluated in the dense zone of the spray by comparison to a DNS based on a coupled level set/VOF/ghost fluid method. Once constants and parameters of the model are fixed thanks to comparisons with DNS, the model is tested with experimental data. The liquid and vapour penetrations show a good agreement when they are compared to experiments of Diesel atomization. In particular the influence of the gas temperature is well recovered. For different temperatures, similarly to the experiments, vapour penetrations are unchanged, but the corresponding equivalent ratio fields are strongly modified. Finally, the combustion model ECFM-3Z [Colin, O., Benkenida, A., 2004. The 3-zones extended coherent flame model (ecfm-3z) for computing premixed/diffusion combustion. Oil Gas Sci. Technol. 59 (6) 593–609] is joined to the ELSA model and the effect of gas temperature changes on a Diesel spray flame is reproduced.     

Computational analysis of binary collisions of shear-thinning droplets

Scale-reduced models of transport processes and reactive mixing in sprays require improved closure laws, taking into account the characteristic features of elementary spray processes. The present paper investigates binary droplet collisions as such an elementary process. In the case of shear-thinning liquids considered here, this requires a profound understanding of the influence of the non-Newtonian fluid rheology on the flow inside the colliding drops and the collision complex dynamics. We employ direct numerical simulations based on the Volume-of-Fluid method to study these collisions. The results give a quantitative prediction of the resulting droplet collision diameter as well as a qualitative prediction of the complete time evolution. During collisions, extremely thin fluid lamellae appear inside the expanding complex. These have to be accounted for in a physically sound simulation and we apply a stabilization of the lamella to keep it from rupturing. The simulations show that in all considered cases an effective constant viscosity can be found a posteriori which leads to the same collision dynamics. But this effective viscosity is neither the mean nor the minimum viscosity.     

The effect of viscoelasticity on stress fields within polyethylene melt flow for a cross-slot and contraction–expansion slit geometry

The sensitivity of the principal stress difference (PSD) profiles to material viscoelasticity is demonstrated for two flow geometries using three different polyethylenes. Studies were performed using both experimental optical techniques and computational simulations, in the latter case to evaluate the ability to model these complex flows. The materials were characterised using linear and extensional rheology which was fitted to a multimode POM-POM model implemented in the Lagrangian–Eulerian code flowSolve. A contraction–expansion (CE) slit geometry was used to create a mixed, but primarily simple shear flow, whilst a cross-slot geometry provided a region of high extensional shear and high strain. In both flows, the PSD developed from an initial Newtonian profile to increasing levels of asymmetry between the inlet and the outlet flow. More specific phenomena, such as downstream stress fangs in the CE slit and the formation of centreline cusps and “W”-shaped cusps in the cross-slot, were also observed. The simulations of PSD development within the CE slit geometry quantitatively captured the experimental results. In the case of the cross-slot geometry, the qualitative features of the PSD development were well captured, although the results were quantitatively less accurate.     

Parallel adaptive high-order CFD simulations characterising SOFIA cavity acoustics

ABSTRACT

This paper presents large-scale parallel computational fluid dynamics simulations for the Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA is an airborne, 2.5-m infrared telescope mounted in an open cavity in the aft fuselage of a Boeing 747SP. These simulations focus on how the unsteady flow field inside and over the cavity interferes with the optical path and mounting structure of the telescope. A temporally fourth-order accurate Runge–Kutta, and a spatially fifth-order accurate WENO-5Z scheme were used to perform implicit large eddy simulations. An immersed boundary method provides automated gridding for complex geometries and natural coupling to a block-structured Cartesian adaptive mesh refinement framework. Strong scaling studies using NASA's Pleiades supercomputer with up to 32 k CPU cores and 4 billion computational cells show excellent scaling. Dynamic load balancing based on execution time on individual adaptive mesh refinement (AMR) blocks addresses irregular numerical cost associated with blocks containing boundaries. Limits to scaling beyond 32 k cores are identified, and targeted code optimisations are discussed.     

Three-dimensional unsteady flow simulations: Alternative strategies for a volume-averaged calculation

This work builds on a SIMPLE-type code to produce two numerical codes of greatly improved speed and accuracy for solution of the Navier–Stokes equations. Both implicit and explicit codes employ an improved QUICK (quadratic upstream interpolation for convective kinematics) scheme to finite difference convective terms for non-uniform grids. The PRIME (update pressure implicit, momentum explicit) algorithm is used as the computational procedure for the implicit code. Use of both the ICCG (incomplete Cholesky decomposition, conjugate gradient) method and the MG (multigrid) technique to enhance solution execution speed is illustrated. While the implicit code is first-order in time, the explicit is second-order accurate. Two- and three-dimensional forced convection and sidewall-heated natural convection flows in a cavity are chosen as test cases. Predictions with the new schemes show substantial computational savings and very good agreement when compared to previous simulations and experimental data.     

Numerical simulation of the disintegration of forced liquid jets using volume-of-fluid method

The present numerical study investigates the effect of finite sinusoidal velocity modulations imposed on an otherwise unperturbed cylindrical liquid jet issuing into stagnant gas using Volume-of-Fluid (VOF) methodology. Variation of the simulation parameters, comprising of the mean liquid jet velocity, modulation amplitude and frequency grouped together using a set of non-dimensional parameters, leads to the formation of a wide gamut of reproducible liquid structures such as surface waves, upstream/downstream directed bells and chains of droplets similar to those observed in experiments. The computations efficiently capture the diverse flow structures generated by the evolving modulated liquid jet inclusive of several nonlinear dynamics such as growth of surface waves, ligament interaction with shear vortices and its subsequent thinning process. The simulations identify the deterministic behaviour of modulated liquid jets to predict liquid disintegration modes under given set of non-dimensional parameters.     

Development of numerical techniques for near‐field aeroacoustic computations

This work is concerned with the development of a numerical scheme capable of producing accurate simulations of sound propagation in the presence of a mean flow field. The method is based on the concept of variable decomposition, which leads to two separate sets of equations. These equations are the linearised Euler equations and the Reynolds‐averaged Navier–Stokes equations. This paper concentrates on the development of numerical schemes for the linearised Euler equations that leads to a computational aeroacoustics (CAA) code. The resulting CAA code is a non‐diffusive, time‐ and space‐staggered finite volume code for the acoustic perturbation, and it is validated against analytic results for pure 1D sound propagation and 2D benchmark problems involving sound scattering from a cylindrical obstacle. Predictions are also given for the case of prescribed source sound propagation in a laminar boundary layer as an illustration of the effects of mean convection. Copyright © 1999 John Wiley & Sons, Ltd.     

A comparison study of discrete-ordinates and flux-limited diffusion methods for modeling radiation transport

The Center for Radiative Shock Hydrodynamics (CRASH) is investigating methods of improving the predictive capability of numerical simulations for radiative shock waves that are produced in Omega laser experiments. The laser is used to shock, ionize, and accelerate a beryllium foil into a xenon-filled shock tube. These shock waves, when driven above a threshold velocity of about 60 km/s, become strongly radiative and convert much of the incident energy flux into radiation.Radiative shocks have properties that are significantly different from purely hydrodynamic shocks and, in modeling this phenomenon numerically, it is important to compute radiative effects accurately. In this article, we examine approaches to modeling radiation transport by comparing two methods: (i) a computationally efficient, multigroup, flux-limited-diffusion approximation, currently in use in the CRASH radiation-hydrodynamics code, with (ii) a more accurate discrete-ordinates treatment that is offered by the radiation-transport code PDT. We present a selection of results from a growing suite of code-to-code comparison tests, showing both results for idealized problems and for those that are representative of conditions found in the CRASH experiment.     

Atomization and spray characteristics of bioethanol and bioethanol blended gasoline fuel injected through a direct injection gasoline injector

The focus of this study was to investigate the spray characteristics and atomization performance of gasoline fuel (G100), bioethanol fuel (E100), and bioethanol blended gasoline fuel (E85) in a direct injection gasoline injector in a gasoline engine. The overall spray and atomization characteristics such as an axial spray tip penetration, spray width, and overall SMD were measured experimentally and predicted by using KIVA-3V code.The development process and the appearance timing of the vortices in the test fuels were very similar. In addition, the numerical results accurately described the experimentally observed spray development pattern and shape, the beginning position of the vortex, and the spray breakup on the spray surface. Moreover, the increased injection pressure induced the occurrence of a clear circular shape in the downstream spray and a uniform mixture between the injected spray droplets and ambient air. The axial spray tip penetrations of the test fuels were similar, while the spray width and spray cone angle of E100 were slightly larger than the other fuels. In terms of atomization performance, the E100 fuel among the tested fuels had the largest droplet size because E100 has a high kinematic viscosity and surface tension.     

Properties of acoustic resonance in double-actuator ultra-sonic gas nozzle: numerical study

The ultra-sonic gas atomization(USGA) nozzle is an important apparatus in the metal liquid air-blast atomization process.It can generate oscillating supersonic gas efflux,which is proved to be effective to enforce the atomization and produce narrow-band particle distributions.A double-actuator ultra-sonic gas nozzle is proposed in the present paper by joining up two active signals at the ends of the resonance tubes.Numerical simulations are adopted to study the effects of the flow development on the acoustic resonant properties inside the Hartmann resonance cavity with/without actuators.Comparisons show that the strength and the onset process of oscillation are enhanced remarkably with the actuators.The multiple oscillating amplitude peaks are found on the response curves,and two kinds of typical behaviors,i.e.,the Hartmann mode and the global mode,are discussed for the corresponding frequencies.The results for two driving actuators are also investigated.When the amplitudes,the frequencies,or the phase difference of the input signals of the actuators are changed,the oscillating amplitudes of gas efflux can be altered effectively.     

High-order overset grid method for detecting particle impaction on a cylinder in a cross flow

An overset grid method was developed to investigate the interaction between a particle-laden flow and a circular cylinder. The method is implemented in the Pencil Code, a high-order finite-difference code for compressible flow simulation. High-order summation-by-parts operators were used at the cylinder boundary, and both bi-linear Lagrangian and bi-quadratic spline interpolation were used to communicate between the Cartesian background grid and the body-conformal cylindrical grid. The performance of the overset grid method was assessed to benchmark cases of steady and unsteady flows past a cylinder. Results show high-order accuracy and good agreement to the literature. Particle-laden flow simulations were performed, with inertial point particles impacting on a cylinder. The simulations reproduced results from the literature at a significantly reduced cost. Further, an investigation into blockage effects on particle impaction revealing that the previously published DNS data is less accurate than assumed for particles with very small Stokes numbers.     

A computational study of shock speeds in high-performance shock tubes

This paper describes U2DE, a finite-volume code that numerically solves the Euler equations. The code was used to perform multi-dimensional simulations of the gradual opening of a primary diaphragm in a shock tube. From the simulations, the speed of the developing shock wave was recorded and compared with other estimates. The ability of U2DE to compute shock speed was confirmed by comparing numerical results with the analytic solution for an ideal shock tube. For high initial pressure ratios across the diaphragm, previous experiments have shown that the measured shock speed can exceed the shock speed predicted by one-dimensional models. The shock speeds computed with the present multi-dimensional simulation were higher than those estimated by previous one-dimensional models and, thus, were closer to the experimental measurements. This indicates that multi-dimensional flow effects were partly responsible for the relatively high shock speeds measured in the experiments. Received 15 November 1996 / Accepted 3 February 1997