Numerical Modelling of Diesel Spray Using the Eulerian Multiphase Approach

This research investigates high pressure diesel fuel injection into the combustion chamber by performing computational simulations using the Euler-Eulerian multiphase approach. Six diesel-like conditions were simulated for which the liquid fuel jet was injected into a pressurised inert environment (100 % N2) through a 205 µm nozzle hole. The analysis was focused on the liquid jet and vapour penetration, describing spatial and temporal spray evolution. For this purpose, an Eulerian multiphase model was implemented, variations of the sub-model coefficients were performed, and their impact on the spray formation was investigated. The final set of sub-model coefficients was applied to all operating points. Several simulations of high pressure diesel injections (50, 80, and 120 MPa) combined with different chamber pressures (5.4 and 7.2 MPa) were carried out and results were compared to the experimental data. The predicted results share a similar spray cloud shape for all conditions with the different vapour and liquid penetration length. The liquid penetration is shortened with the increase in chamber pressure, whilst the vapour penetration is more pronounced by elevating the injection pressure. Finally, the results showed good agreement when compared to the measured data, and yielded the correct trends for both the liquid and vapour penetrations under different operating conditions

Numerical modelling of landslide-tsunami propagation in a wide range of idealised water body geometries

© 2019 Elsevier B.V. Large landslide-tsunamis are caused by mass movements such as landslides or rock falls impacting into a water body. Research of these phenomena is essentially based on the two idealised water body geometries (i) wave flume (2D, laterally confined wave propagation) and (ii) wave basin (3D, unconfined wave propagation). The wave height in 2D and 3D differs by over one order of magnitude in the far field. Further, the wave characteristics in intermediate geometries are currently not well understood. This article focuses on numerical landslide-tsunami propagation in the far field to quantify the effect of the water body geometry. The hydrodynamic numerical model SWASH, based on the non-hydrostatic non-linear shallow water equations, was used to simulate approximate linear, Stokes, cnoidal and solitary waves in 6 different idealised water body geometries. This includes 2D, 3D as well as intermediate geometries consisting of “channels” with diverging side walls. The wavefront length was found to be an excellent parameter to correlate the wave decay along the slide axis in all these geometries in agreement with Green's law and with diffraction theory in 3D. Semi-theoretical equations to predict the wave magnitude of the idealised waves at any desired point of the water bodies are also presented. Further, simulations of experimental landslide-tsunami time series were performed in 2D to quantify the effect of frequency dispersion. This process may be negligible for solitary- and cnoidal-like waves for initial landslide-tsunami hazard assessment but becomes more important for Stokes-like waves in deeper water. The findings herein significantly improve the reliability of initial landslide-tsunami hazard assessment in water body geometries between 2D and 3D, as demonstrated with the 2014 landslide-tsunami event in Lake Askja

Numerical Modelling of Optical Trapping

Optical trapping is a widely used technique, with many important applications in biology and metrology. Complete modelling of trapping requires calculation of optical forces, primarily a scattering problem, and non-optical forces. The T-matrix method is used to calculate forces acting on spheroidal and cylindrical particles.Comment: 4 pages, 4 figure

Numerical Modelling of Grate Combustion

Předkládaná práce je zaměřena na numerické modelování spalování tuhých paliv na roštu metodami výpočtové dynamiky tekutin (CFD). Jelikož výsledky CFD simulací roštového spalování závisí na kvalitě vstupních dat, která zahrnují i údaje o teplotě, hmotnostním toku a chemickém složení spalin vystupujících z lože, pozornost je věnována především procesům, probíhajícím v loži během spalování na roštu. Velká část práce je věnována vývoji spolehlivého modelu spalování v sypaných ložích, jelikož může napomoci zkvalitnit výsledky simulací i rozšířit znalosti principů spalování tuhých paliv v sypaných ložích. V rámci práce byl vyvinut jednorozměrný nestacionární model spalování v experimentálním reaktoru a implementován do počítačového programu GRATECAL 1.3 včetně grafického uživatelského rozhraní. Zvláštní důraz byl kladen na konzervativnost modelu. Proto byla vyvinuta metoda pro kontrolu hmotnostní a energetické bilance systému a následně aplikována v řadě studií, v rámci nichž byly odhaleny některé chyby týkající se definic zdrojových členů, které byly převzaty z literatury a opraveny. Pomocí modelu byla provedena analýza šíření čela sušení a reakce hoření koksu po výšce lože pšeničné slámy. Na základě výsledků těchto analýz bylo doporučeno zahrnout i modelování změny porozity částic paliva, aby šířka reakční zóny byla predikována korektně v případě, že je uvažována změna porozity celého lože. Rovněž vyvinutá bilanční metoda byla použita k analýze vlivu kritérií konvergence na hmotnostní a energetickou nerovnováhu simulovaného systému. Bylo zjištěno, že škálovaná rezidua rovnic všech veličin by měla poklesnout aspoň na hodnotu $10^{-6}$, aby bylo dosaženo nízké hmotnostní a energetické nerovnováhy a tudíž uspokojivě přesných výsledků ze simulací v loži. Druhá část práce je věnována vývoji a implementaci knihovny uživatelem definovaných funkcí pro komerční CFD nástroj ANSYS FLUENT, které slouží k propojení modelu lože s modelem komory reálné spalovací jednotky, aby byla umožněna dynamická změna okrajových podmínek na vstupu do komory v závislosti na výstupech ze simulací v loži. Vytvořené rozhraní pro propojení těchto dvou modelů je dostatečně obecné pro aplikaci na širokou škálu modelů roštových kotlů. Popsané výsledky přispívají k lepšímu porozumění numerickému modelování spalování na roštu, a to zejména ve fázi sestavování numerického modelu a nastavení parametrů řešiče pro kontrolu konvergence.The present work is focused on numerical modelling of grate combustion of solid fuels by means of computational fluid dynamics (CFD) methods. Since CFD results from simulations of grate combustion depend on the quality of input data including information on temperature, mass flux and chemical composition of flue gas leaving a fuel bed, the attention is turned to modelling of processes, that take place within the fuel bed on a grate. A great part of the work is devoted to development of a reliable numerical model of packed-bed combustion as it may help improve both results from simulations and knowledge of principles of solid fuel combustion in fixed or moving beds. A one-dimensional transient numerical model of combustion in an experimental reactor is developed and implemented into a computer program called GRATECAL 1.3 with a grapical user interface. A special emphasis is put on the conservativeness property of the model. Therefore, a method for control of mass and energy balance over the system is developed and applied to a series of case studies, which have revealed certain errors in definitions of mass source terms, so that data adopted from literature have been reconciled. The model is used for analysis of propagation of drying and char combustion reaction fronts in a bed of wheat straw particles. It is suggested to include modelling of particle internal porosity change in order to obtain correct reaction zone thickness, if porosity of the bed is allowed to change during combustion. The balance-based method is also used to analyse effects of convergence criteria on mass and energy imbalance of the modelled system. It is found that all the scaled residuals must drop to as low as $10^{-6}$ or lower in order to obtain sufficiently accurate results from in-bed simulations in terms of mass and energy conservation within the packed bed. The second part of the work is devoted to development of a library of user-defined functions for the commercial CFD software ANSYS FLUENT for coupling the bed model with a freeboard model of a real combustion unit in order to specify the boundary conditions indirectly using results from in-bed simulations. The created interface is general enough to be used for a wide range of models of grate furnaces. The presented results contribute to better understanding of numerical modelling of grate combustion, especially in the setup of a numerical model and parameters of solver for the control of the convergence.

Numerical modelling of metal melt refining process in ladle with rotating impeller and breakwaters

The paper describes research and development of aluminium melt refining technology in a ladle with rotating impeller and breakwaters using numerical modelling of a finite volume/element method. The theoretical aspects of refining technology are outlined. The design of the numerical model is described and discussed. The differences between real process conditions and numerical model limitations are mentioned. Based on the hypothesis and the results of numerical modelling, the most appropriate setting of the numerical model is recommended. Also, the possibilities of monitoring of degassing are explained. The results of numerical modelling allow to improve the refining technology of metal melts and to control the final quality under different boundary conditions, such as rotating speed, shape and position of rotating impeller, breakwaters and intensity of inert gas blowing through the impeller.Web of Science64266465

Numerical modelling of Bose-Einstein correlations

We propose extension of the algorithm for numerical modelling of Bose-Einstein correlations (BEC), which was presented some time ago in the literature. It is formulated on quantum statistical level for a single event and uses the fact that identical particles subjected to Bose statistics do bunch themselves, in a maximal possible way, in the same cells in phase-space. The bunching effect is in our case obtained in novel way allowing for broad applications and fast numerical calculations. First comparison with $e^+e^-$ annihilations data performed by using simple cascade hadronization model is very encouraging.Comment: LaTeX file and 5 eps file with figures, 9 pages altogethe

Proposition of numerical modelling of BEC

We propose extension of the numerical method to model effect of Bose-Einstein correlations (BEC) observed in hadronization processes which allows for calculations not only correlation functions $C_2(Q_{inv})$ (one-dimensional) but also corresponding to them $C_2(Q_{x,y,z})$ (i.e., three-dimensional). The method is based on the bunching of identical bosonic particles in elementary emitting cells (EEC) in phase space in manner leading to proper Bose-Einstein form of distribution of energy (this was enough to calculate $C_2(Q_{inv})$). To obtain also $C_2(Q_{x,y,z})$ one has to add to it also symmetrization of the multiparticle wave function to properly correlate space-time locations of produced particles with their energy-momentum characteristics.Comment: 7 pages, 2 figures, poster presented at QM2005, to be published in Nukleonika (2006); acknowledgements adde

Numerical modelling of functionally graded coatings

Peer reviewedPostprin