Search

In Lagrangian stochastic collision models, a fictitious particle is generated to act as a collision partner, with a velocity correlated to the velocity of the real colliding particle. However, most often, the fluid velocity seen by this fictitious particles is not accounted for in the generation of the fictitious particle velocity, leading to a de-correlation between the fictitious particle velocity and the local fluid velocity, which, after collision, leads to an unrealistic de-correlation of the real particle velocity and the fluid velocity as seen by the particle. This de-correlation, in turn, causes a spurious decrease of the particle kinetic energy, even though the collisions are assumed perfectly elastic. In this paper, we propose a new model in which the generated fictitious particle velocity is correctly correlated to both the real particle velocity and the local fluid velocity at the particle, hence preventing the spurious loss of the total particle kinetic energy. The model is suitable for small inertial particles. Two algorithms for integrating the collision frequency are also compared to each other. The models are validated using large eddy simulation (LES) of mono-dispersed particle-laden stationary homogeneous isotropic turbulence. Simulations are conducted with spherical particles with different turbulent Stokes number, S t t = [ 0.75 - 5.8 ] , and volume fractions, α p = [ 0.014 - 0.044 ] , and are compared to the results of the LES using a deterministic discrete particle simulation model. [ABSTRACT FROM AUTHOR]

In this work, the reconstruction of sub-grid scales in large eddy simulation (LES) of turbulent flows in stratocumulus clouds is addressed. The approach is based on the fractality assumption of turbulent velocity field. The fractal model reconstructs sub-grid velocity field from known filtered values on LES grid, by means of fractal interpolation, proposed by Scotti and Meneveau (Physica D 127, 198–232 1999). The characteristics of the reconstructed signal depend on the stretching parameter d, which is related to the fractal dimension of the signal. In many previous studies, the stretching parameter values were assumed to be constant in space and time. To improve the fractal interpolation approach, we account for the stretching parameter variability. The local stretching parameter is calculated from direct numerical simulation (DNS) data with an algorithm proposed by Mazel and Hayes (IEEE Trans. Signal Process 40(7), 1724–1734, 1992), and its probability density function (PDF) is determined. It is found that the PDFs of d have a universal form when the velocity field is filtered to wave-numbers within the inertial range. In order to investigate Reynolds number (Re) dependence, we compare the inertial-range PDFs of d in DNS and large eddy simulation (LES) of stratocumulus cloud-top and experimental airborne data from Physics of Stratocumulus Top (POST) research campaign. Next, fractal reconstruction of the subgrid velocity is performed and energy spectra and statistics of velocity increments are compared with DNS data. It is assumed that the stretching parameter d is a random variable with the prescribed PDF. We show that the agreement with the DNS is in such case better and the error in mass conservation is smaller compared to the use of constant values of d. The motivation of this study is to reproduce effect of sub-grid scales on a motion of Lagrangian particles (e.g. droplets) in clouds. [ABSTRACT FROM AUTHOR]

Airworthiness standards require a fire resistance demonstration for aircraft or helicopter engines to obtain a type certificate. This demonstration relies on tests performed with prototype engine parts in the late stages of the development. In hardest tests, a kerosene standardized flame with imposed burnt gas temperature and heat flux is placed next to the engine casing during a given time. The aim of this work is to provide a better characterization of a kerosene/air certification burner in order to reach a better understanding of the thermal environment during fire tests. To this purpose, Large-Eddy Simulation (LES) of the certification burner is carried out. Spray combustion, forced convection on walls and conduction in the solid parts of the burner are coupled to achieve a detailed description of heat transfer. In a first place, physical aspects involved inside the burner in an adiabatic case are described. Then, differences that exist with a conjugate convective and conductive heat transfer case are analyzed. To a larger extent, the aim is to have a better characterization of the flow impinging the casing and to progress on fire test modeling so as to minimize the risks of test failure. [ABSTRACT FROM AUTHOR]

This paper presents highly resolved large eddy simulations (LES) of an internal combustion engine (ICE) using an immersed boundary method (IBM), which can describe moving and stationary boundaries in a simple and efficient manner. In this novel approach, the motion of the valves and the piston is modeled by Lagrangian particles, whilst the stationary parts of the engine are described by a computationally efficient IBM. The proposed mesh-free technique of boundary representation is simple for parallelization and suitable for high performance computing (HPC). To demonstrate the method, LES results are presented for the flow and the combustion in an internal combustion engine. The Favre-filtered Navier-Stokes equations are solved for a compressible flow employing a finite volume method on Cartesian grids. Non-reflecting boundary conditions are applied at the intake and the exhaust ports. Combustion is described using a flame surface density (FSD) model with an algebraic reaction rate closure. A simplified engine with a fixed axisymmetric valve (see Appendix A) is employed to show the correctness of the method while avoiding the uncertainties which may be induced by the complex engine geometry. Three test-cases using a real engine geometry are investigated on different grids to evaluate the impact of the cell size and the filter width. The simulation results are compared against the experimental data. A good overall agreement was found between the measurements and the simulation data. The presented method has particular advantages in the efficient generation of the grid, high resolution and low numerical dissipation throughout the domain and an excellent suitability for massively parallel simulations. [ABSTRACT FROM AUTHOR]