Development of accurate, robust liquid equations of state for multi-phase CFD simulations with a modified AUSM+-up scheme

Eight equations of state (EOS) have been evaluated for the simulation of compressible liquid water properties, based on empirical correlations, the principle of corresponding states and thermodynamic relations. The IAPWS-IF97 EOS for water was employed as the reference case. These EOSs were coupled to a modified AUSM+-up convective flux solver to determine flow profiles for three test cases of differing flow conditions. The impact of the non-viscous interaction term discretisation scheme, interfacial pressure method and selection of low-Mach number diffusion were also compared. It was shown that a consistent discretisation scheme using the AUSM+-up solver for both the convective flux and the non-viscous interfacial term demonstrated both robustness and accuracy whilst facilitating a computationally cheaper solution than discretisation of the interfacial term independently by a central scheme. The simple empirical correlations gave excellent results in comparison to the reference IAPWS-IF97 EOS and were recommended for developmental work involving water as a cheaper and more accurate EOS than the more commonly used stiffened-gas model. The correlations based on the principles of corresponding-states and the modified Peng-Robinson cubic EOS also demonstrated a high degree of accuracy, which is promising for future work with generic fluids. Further work will encompass extension of the solver to multiple dimensions and to account for other source terms such as surface tension, along with the incorporation of phase changes. © 2013.

Direct numerical simulation of turbulent flows over superhydrophobic surfaces with gas pockets using linearized boundary conditions

Superhydrophobic surfaces are shown to be effective for surface drag reduction under laminar regime by both experiments and simulations (see for example, Ou and Rothstein, Phys. Fluids 17:103606, 2005). However, such drag reduction for fully developed turbulent flow maintaining the Cassie-Baxter state remains an open problem due to high shear rates and flow unsteadiness of turbulent boundary layer. Our work aims to develop an understanding of mechanisms leading to interface breaking and loss of gas pockets due to interactions with turbulent boundary layers. We take advantage of direct numerical simulation of turbulence with slip and no-slip patterned boundary conditions mimicking the superhydrophobic surface. In addition, we capture the dynamics of gas-water interface, by deriving a proper linearized boundary condition taking into account the surface tension of the interface and kinematic matching of interface deformation and normal velocity conditions on the wall. We will show results from our simulations predicting the dynamical behavior of gas pocket interfaces over a wide range of dimensionless surface tensions.

Drop speeds from drop-on-demand ink-jet print heads

Measured drop speeds from a range of industrial drop-on-demand (DoD) ink-jet print head designs scale with the predictions of very simple physical models and results of numerical simulations. The main drop/jet speeds at a specified stand-off depend on fluid properties, nozzle exit diameter, and print head drive amplitude for fixed waveform timescales. Drop speeds from the Xaar, Spectra Dimatix, and MicroFab DoD print heads tested with (i) Newtonian, (ii) weakly elastic, and (iii) highly shear-thinning fluids all show a characteristic linear rise with drive voltage (setting) above an apparent threshold drive voltage. Jetting, simple modeling approaches, and numerical simulations of Newtonian fluids over the typical DoD printing range of surface tensions and viscosities were studied to determine how this threshold drive value and the slope of the characteristic linear rise depend on these fluid properties and nozzle exit area. The final speed is inversely proportional to the nozzle exit area, as expected from volume conservation. These results should assist specialist users in the development and optimization of DoD applications and print head design. For a given density, the drive threshold is determined primarily by viscosity, and the constant of proportionality k linking speed with drive above a drive threshold becomes independent of viscosity and surface tension for more viscous DoD fluid jetting. © 2013 Society for Imaging Science and Technology.

Effects of the biodiesel fuel physical properties on the swirl stabilised spray combustion characteristics

An increasin g interest in biofuel applications in modern engines requires a better understanding of biodiesel combustion behaviour. Many numerical studies have been carried out on unsteady combustion of biodiesel in situations similar to diesel engines, but very few studies have been done on the steady combustion of biodiesel in situations similar to a gas turbine combustor environment. The study of biodiesel spray combustion in gas turbine applications is of special interest due to the possible use of biodiesel in the power generation and aviation industries. In modelling spray combustion, an accurate representation of the physical properties of the fuel is a first important step, since spray formation is largely influenced by fuel properties such as viscosity, density, surface tension and vapour pressure. In the present work, a calculated biodiesel properties database based on the measured composition of Fatty Acid Methyl Esters (FAME) has been implemented in a multi-dimensional Computational Fluid Dynamics (CFD) spray simulation code. Simulations of non-reacting and reacting atmospheric-pressure sprays of both diesel and biodiesel have been carried out using a spray burner configuration for which experimental data is available. A pre-defined droplet size probability density function (pdf) has been implemented together with droplet dynamics based on phase Doppler anemometry (PDA) measurements in the near-nozzle region. The gas phase boundary condition for the reacting spray cases is similar to that of the experiment which employs a plain air-blast atomiser and a straight-vane axial swirler for flame stabilisation. A reaction mechanism for heptane has been used to represent the chemistry for both diesel and biodiesel. Simulated flame heights, spray characteristics and gas phase velocities have been found to compare well with the experimental results. In the reacting spray cases, biodiesel shows a smaller mean droplet size compared to that of diesel at a constant fuel mass flow rate. A lack of sensitivity towards different fuel properties has been observed based on the non-reacting spray simulations, which indicates a need for improved models of secondary breakup. By comparing the results of the non-reacting and reacting spray simulations, an improvement in the complexity of the physical modelling is achieved which is necessary in the understanding of the complex physical processes involved in spray combustion simulation. Copyright © 2012 SAE International.

Spontaneous formation of highly ordered nanostructures: thermal instability and mode selection in surface-capped polymer films.

We demonstrate a controllable formation process of wave-like patterns in thermally unstable surface-capped polymer films on a rigid substrate. Self-ordered wave-like structures over a large area can be created by applying a small lateral tension to the film, whereupon it becomes unstable. A clear mode selection process which includes creation, decay and interference between coexisting waves at different annealing conditions has been observed, which makes it possible to restrain the patterns which are formed finally. Our results provide a clear and new evidence of spinodal behaviour in such a film due to thermal instability. Furthermore, we show that the well-controlled patterns generated in such a process can be used to fabricate nanostructures for various applications.

Important aspects when modelling the interaction between surface structures and tunnelling in sand

Tunnelling in urban areas continues to increase and has highlighted the need for a better understanding of the impact of tunnel excavations on existing buildings. This paper considers the influence of surface structures on ground displacements caused by tunnelling in sand through finite element modelling and centrifuge testing. First, the importance of modelling assumptions is evaluated by comparing centrifuge modelling results to finite element modelling results for various soil constitutive models: both a Young's modulus that linearly increases with depth and a power law relation between the soil stiffness and stresses are considered. Second, the most effective soil constitutive model was used to perform a sensitivity study on the effect of different factors governing the structural response. In particular, the effect of the building stiffness and weight on the modification of soil displacements is investigated by introducing a simple surface structure. The use of a no-tension interface between the building and the soil was found to be essential to investigate the effect of weight on gap formation between the soil and the structure, as observed during the experimental tests. Results show the importance of considering the relation between the building weight and the relative stiffness between the building and the soil when assessing the structural response. © 2014 Korean Geotechnical Society.