No Arabic abstract
A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem, and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model, the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition. The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.
A quasi-one-dimensional analytic model is proposed for the internal fluid of rotating detonation combustors (RDCs). This model uses the shock-tube model that constrains the flow to have only a longitudinal component, while admitting the propagation of the detonation wave in the azimuthal direction. The proposed model is able to compute the thruster performance and two-dimensional distributions of gas properties. The calculation process of the model excludes iterative calculation or space discretization. The case calculations of the hydrogen-air RDC and the ethylene-oxygen RDC are conducted, and the results calculated by the analytic model are compared with those simulated by computational fluid dynamics (CFD). Good agreement has been observed between the results obtained with the proposed model and CFD, in terms of both of the qualitative and quantitative comparisons. The proposed model is simple and fast, and also maintains the fundamental characteristics of RDCs.
Historical experimental testing of high-altitude nuclear explosions (HANEs) are known to cause severe and detrimental effects to radio frequency signals and communications infrastructure. In order to study and predict the impact of HANEs, tractable computational approaches are required to model the complex physical processes involved in the detonation wave physics. Modern reduced-order models (ROMs) can enable long-time and many-parameter simulations with minimal computational cost. However, translational and scale invariances inherent to this type of wave propagation problem are known to limit traditional ROM approaches. Specifically, dimensionality reduction methods are typically ineffective in producing low-rank models when invariances are present in the data. In this work, an unsupervised machine learning method is used to discover coordinate systems that make such invariances amenable to traditional dimensionality reduction methods. The method, which has previously been demonstrated on one-dimensional translations, is extended to higher dimensions and additional invariances. A surrogate HANE system, i.e. a HANE-ROM, with one detonation wave is captured well at extremely low-rank. Two detonation-waves are also considered with various amounts of interaction between the waves, with improvements to low-rank models for multiple wave quantities with limited interaction.
The present study addresses the reaction zone structure and burning mechanism of unstable detonations. Experiments investigated mainly two-dimensional methane-oxygen cellular detonations in a thin channel geometry. The sufficiently high temporal resolution permitted to determine the PDF of the shock distribution, a power-law with an exponent of -3, and the burning rate of unreacted pockets from their edges - through surface turbulent flames with a speed approximately 3-7 times larger than the laminar one at the local conditions. Numerical simulations were performed using a novel Large Eddy Simulation method where the reactions due to both auto-ignition and turbulent transport and treated exactly at the sub-grid scale in a reaction-diffusion formulation. The model is an extension of Kerstein & Menons Linear Eddy Model for Large Eddy Simulation to treat flows with shock waves and rapid gasdynamic transients. The two-dimensional simulations recovered well the amplification of the laminar flame speed owing to the turbulence generated mainly by the shear layers originating from the triple points and subsequent Richtmyer-Meshkov instability associated with the internal pressure waves. The simulations clarified how the level of turbulence generated controlled the burning rate of the pockets, the hydrodynamic thickness of the wave, the cellular structure and its distribution. Three-dimensional simulations were found in general good agreement with the two-dimensional ones, in that the sub-grid scale model captured the ensuing turbulent burning once the scales associated with the cellular dynamics, where turbulent kinetic energy is injected, are well resolved.
We begin with the theoretical study of spectral energy cascade due to the propagation of high amplitude sound in the absence of thermal sources. To this end, a first-principles-based system of governing equations, correct up to second order in perturbation variables is derived. The exact energy corollary of such second-order system of equations is then formulated and used to elucidate the spectral energy dynamics of nonlinear acoustic waves. We then extend this analysis to thermoacoustically unstable waves -- i.e. amplified as a result of thermoacoustic instability. We drive such instability up until the generation of shock waves. We further study the nonlinear wave propagation in geometrically complex case of waves induced by the spark plasma between the electrodes. This case adds the geometrical complexity of a curved, three-dimensional shock, yielding vorticity production due to baroclinic torque. Finally, detonation waves are simulated by using a low-order approach, in a periodic setup subjected to high pressure inlet and exhaust of combustible gaseous mixture. An order adaptive fully compressible and unstructured Navier Stokes solver is currently under development to enable higher fidelity studies of both the spark plasma and detonation wave problem in the future.
Discrete Boltzmann model (DBM) is a type of coarse-grained mesoscale kinetic model derived from the Boltzmann equation. Physically, it is roughly equivalent to a hydrodynamic model supplemented by a coarse-grained model for the relevant thermodynamic non-equilibrium (TNE) behaviours. The Navier-Stokes (NS) model is a traditional macroscopic hydrodynamic model based on continuity hypothesis and conservation laws. In this study, the two models are compared from two aspects, physical capability and computational cost, by simulating two kinds of flow problems including the thermal Couette flow and a Mach 3 step problem. In the cases where the TNE effects are weak, both the two models give accurate results for the hydrodynamic behaviour. Besides, DBM can provide more detailed non-equilibrium information, while the NS is more efficient if concern only the density, momentum, energy and their derived quantities. It is concluded that, if the TNE effects are strong or are to be investigated, the NS is insufficient while DBM is a good choice. While in the cases where the TNE effects are weak and only the macro flow fields are to be studied, the NS is more preferable.