No Arabic abstract
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.
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.
Monitoring the dynamics processes in combustors is crucial for safe and efficient operations. However, in practice, only limited data can be obtained due to limitations in the measurable quantities, visualization window, and temporal resolution. This work proposes an approach based on neural differential equations to approximate the unknown quantities from available sparse measurements. The approach tackles the challenges of nonlinearity and the curse of dimensionality in inverse modeling by representing the dynamic signal using neural network models. In addition, we augment physical models for combustion with neural differential equations to enable learning from sparse measurements. We demonstrated the inverse modeling approach in a model combustor system by simulating the oscillation of an industrial combustor with a perfectly stirred reactor. Given the sparse measurements of the temperature inside the combustor, upstream fluctuations in compositions and/or flow rates can be inferred. Various types of fluctuations in the upstream, as well as the responses in the combustor, were synthesized to train and validate the algorithm. The results demonstrated that the approach can efficiently and accurately infer the dynamics of the unknown inlet boundary conditions, even without assuming the types of fluctuations. Those demonstrations shall open a lot of opportunities in utilizing neural differential equations for fault diagnostics and model-based dynamic control of industrial power systems.
In studies on instabilities of flowfield in rotating detonation, one of the most common concerns is the instability at the slip line originating from the conjunction of the detonation wave and oblique shock. Using Euler equations associated with 7-species-and-8-reaction finite-rate chemical reaction model of hydrogen/air mixtures, further studies are performed to simulate the 2-D rotating detonation, and the flow mechanism of instability at the slip line is investigated in depth. The results show that the distinct wake profile exists at the slip line, which is different from the typical mixing layer. Analysis indicates that the generation of wake is caused by the transition shock between the detonation wave and oblique shock. Because of the wake profile, the vorticity distribution therein appears in a double-layer layout, and different evolution exist in different vorticity layers. Based on the velocity profile across the slip line, the analysis by the linear stability theory is made, and two unstable modes which have different shape profiles and phase velocities are found. Discrete Fourier transformation is utilized to analyze the numerical results, and similar shape profiles are obtained. A general coincidence in velocity of vortex movement is also attained between the theoretical predictions and simulations. Investigations show that the wake instability is responsible for the unstable mechanism, and corresponding unstable structures differs from the canonical ones in typical mixing layers.
Conflict between formation of a cyclonic vortex and isotropization in forced homogeneous rotating turbulence is numerically investigated. It is well known that a large rotation rate of the system induces columnar vortices to result in quasi-two-dimensional (Q2D) flow, while a small rotation rate allows turbulence to be three-dimensional (3D). It is found that the transition from the Q2D turbulent flow to the 3D turbulent flow and the reverse transition occur at different values of the rotation rates. At the intermediate rotation rates, bistability of these two statistically steady states is observed. Such hysteretic behavior is also observed for the variation of the amplitude of an external force.
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.