CFD-DEM simulation of propagation of sound waves in fluid particles fluidised medium

In this work, speed of sound in 2 phase mixture has been explored using CFD-DEM (Computational Fluid Dynamcis - Discrete Element Modelling). In this method volume averaged Navier Stokes, continuity and energy equations are solved for fluid. Particles are simulated as individual entities; their behaviour is captured by Newton's laws of motion and classical contact mechanics. Particle-fluid interaction is captured using drag laws given in literature.The speed of sound in a medium depends on physical properties. It has been found experimentally that speed of sound drops significantly in 2 phase mixture of fluidised particles because of its increased density relative to gas while maintaining its compressibility. Due to the high rate of heat transfer within 2 phase medium as given in Roy et al. (1990), it has been assumed that the fluidised gas-particle medium is isothermal.The similar phenomenon has been tried to be captured using CFD-DEM numerical simulation. The disturbance is introduced and fundamental frequency in the medium is noted to measure the speed of sound for e.g. organ pipe. It has been found that speed of sound is in agreement with the relationship given in Roy et al. (1990). Their assumption that the system is isothermal also appears to be valid.

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.

Model based control for closed loop testing of 1-D engine simulation models

1-D engine simulation models are widely used for the analysis and verification of air-path design concepts and prediction of the resulting engine transient response. The latter often requires closed loop control over the model to ensure operation within physical limits and tracking of reference signals. For this purpose, a particular implementation of Model Predictive Control (MPC) based on a corresponding Mean Value Engine Model (MVEM) is reported here. The MVEM is linearised on-line at each operating point to allow for the formulation of quadratic programming (QP) problems, which are solved as the part of the proposed MPC algorithm. The MPC output is used to control a 1-D engine model. The closed loop performance of such a system is benchmarked against the solution of a related optimal control problem (OCP). As an example this study is focused on the transient response of a light-duty car Diesel engine. For the cases examined the proposed controller implementation gives a more systematic procedure than other ad-hoc approaches that require considerable tuning effort. © 2012 IFAC.

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.

A flame surface density approach to large-eddy simulation of premixed turbulent combustion

The flame surface density approach to the modeling of premixed turbulent combustion is well established in the context of Reynolds-averaged simulations. For the future, it is necessary to consider large-eddy simulation (LES), which is likely to offer major advantages in terms of physical accuracy, particularly for unsteady combustion problems. LES relies on spatial filtering for the removal of unresolved phenomena whose characteristic length scales are smaller than the computational grid scale. Thus, there is a need for soundly based physical modeling at the subgrid scales. The aim of this paper is to explore the usefulness of the flame surface density concept as a basis for LES modeling of premixed turbulent combustion. A transport equation for the filtered flame surface density is presented, and models are proposed for unclosed terms. Comparison with Reynolds-averaged modeling is shown to reveal some interesting similarities and differences. These were exploited together with known physics and statistical results from experiment and from direct numerical stimulation in order to gain insight and refine the modeling. The model has been implemented in a combustion LES code together with standard models for scalar and momentum transport. Computational results were obtained for a simple three-dimensional flame propagation test problem, and the relative importance of contributing terms in the modeled equation for flame surface density was assessed. Straining and curvature are shown to have a major influence at both the resolved and subgrid levels.

Modelling of settlement induced building damage

This thesis focuses on the modelling of settlement induced damage to masonry buildings. In densely populated areas, the need for new space is nowadays producing a rapid increment of underground excavations. Due to the construction of new metro lines, tunnelling activity in urban areas is growing. One of the consequences is a greater attention to the risk of damage on existing structures. Thus, the assessment of potential damage of surface buildings has become an essential stage in the excavation projects in urban areas (Chapter 1). The current damage risk assessment procedure is based on strong simplifications, which not always lead to conservative results. Object of this thesis is the development of an improved damage classification system, which takes into account the parameters influencing the structural response to settlement, like the non-linear behaviour of masonry and the soil-structure interaction. The methodology used in this research is based on experimental and numerical modelling. The design and execution of an experimental benchmark test representative of the problem allows to identify the principal factors and mechanisms involved. The numerical simulations enable to generalize the results to a broader range of physical scenarios. The methodological choice is based on a critical review of the currently available procedures for the assessment of settlement-induced building damage (Chapter 2). A new experimental test on a 1/10th masonry façade with a rubber base interface is specifically designed to investigate the effect of soil-structure interaction on the tunnelling-induced damage (Chapter 3). The experimental results are used to validate a 2D semi-coupled finite element model for the simulation of the structural response (Chapter 4). The numerical approach, which includes a continuum cracking model for the masonry and a non-linear interface to simulate the soil-structure interaction, is then used to perform a sensitivity study on the effect of openings, material properties, initial damage, initial conditions, normal and shear behaviour of the base interface and applied settlement profile (Chapter 5). The results assess quantitatively the major role played by the normal stiffness of the soil-structure interaction and by the material parameters defining the quasi-brittle masonry behaviour. The limitation of the 2D modelling approach in simulating the progressive 3D displacement field induced by the excavation and the consequent torsional response of the building are overcome by the development of a 3D coupled model of building, foundation, soil and tunnel (Chapter 6). Following the same method applied to the 2D semi-coupled approach, the 3D model is validated through comparison with the monitoring data of a literature case study. The model is then used to carry out a series of parametric analyses on geometrical factors: the aspect ratio of horizontal building dimensions with respect to the tunnel axis direction, the presence of adjacent structures and the position and alignment of the building with respect to the excavation (Chapter 7). The results show the governing effect of the 3D building response, proving the relevance of 3D modelling. Finally, the results from the 2D and 3D parametric analyses are used to set the framework of an overall damage model which correlates the analysed structural features with the risk for the building of being damaged by a certain settlement (Chapter 8). This research therefore provides an increased experimental and numerical understanding of the building response to excavation-induced settlements, and sets the basis for an operational tool for the risk assessment of structural damage (Chapter 9).

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.