ترغب بنشر مسار تعليمي؟ اضغط هنا

Effects of wave interaction on ignition and deflagration-to-detonation transition in ethylene/air mixtures behind a reflected shock

146   0   0.0 ( 0 )
 نشر من قبل Huangwei Zhang
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Dynamics of ethylene autoignition and Deflagration-to-Detonation Transition (DDT) in a one-dimensional shock tube are numerically investigated using a skeletal chemistry including 10 species and 10 reactions. Different combustion modes are investigated through considering various premixed gas equivalence ratios (0.2 to 2.0) and incident shock wave Mach numbers (1.8 to 3.2). Four ignition and DDT modes are observed from the studied cases, i.e., no ignition, deflagration combustion, detonation after reflected shock and deflagration behind the incident shock. For detonation development behind the reflected shock, three autoignition hot spots are formed. The first one occurs at the wall surface after the re-compression of the reflected shock and contact surface, which further develops to a reaction shock because of the explosion in the explosion regime. The other two are off the wall, respectively caused by the reflected shock rarefaction wave interaction and reaction induction in the compressed mixture. The last hot spot develops to a reaction wave and couples with the reflected shock after a DDT process, which eventually leads to detonation combustion. For deflagration development behind the reflected shock, the wave interactions, wall surface autoignition hot spot as well as its induction of reaction shock are qualitatively similar to the mode of detonation after incident shock reflection, before the reflected shock rarefaction wave collision point. However, only one hot spot is induced after the collision, which also develops to a reaction wave but cannot catch up with the reflected shock. For deflagration behind the incident shock, deflagration combustion is induced by the incident shock compression whereas detonation occurs after the shock reflection.

قيم البحث

اقرأ أيضاً

The sensitizing effect of cavities in the form of microbubbles on the shock initiation of a homogeneous liquid explosive is studied computationally. While the presence of voids in an explosive has long been known to induce so-called hot spots that gr eatly accelerate the global reaction rate, the ability to computationally resolve the details of the interaction of the shock front with heterogeneities existing on the scale of the detonation reaction zone has only recently become feasible. In this study, the influence of the spatial distribution of air-filled cavities has been examined, enabled by the use of graphic processing unit (GPU) accelerated computations that can resolve shock initiation and detonation propagation through an explosive while fully resolving features at the mesoscale. Different spatial distributions of cavities are examined in two-dimensional simulations, including regular arrays of cavities, slightly perturbed arrays, random arrays (with varying minimum spacing being imposed on the cavities), and randomly distributed clusters of cavities. The presence of the cavities is able to reduce the time required to initiate detonation---for a given input shock strength---by greater than 50%, in agreement with previous experimental results. Randomly distributing the cavities results in a 15-20% decrease in detonation initiation time in comparison to a regular array of cavities. Clustering the cavities---as would occur in the case of agglomeration---results in an additional 10% decrease in detonation initiation time in comparison to random arrays. The effect of clustering is shown not to be a result of the clusters forming an effectively larger cavity, but rather due to interactions between clusters upon shock loading occurring on the microscale.
The mechanisms governing the low-frequency unsteadiness in the shock wave/turbulent boundary layer interaction at Mach 2 are considered. The investigation is conducted based on the numerical database issued from large-eddy simulations covering approx imately 300 cycles of the low-frequency shock fluctuations. The evaluation of the spectrum in the interaction zone indicates that the broadband low-frequency unsteadiness is predominantly two-dimensional, and can be isolated via spanwise averaging. Empirically derived transfer functions are computed using the averaged flow field, and indicate the occurrence of a feedback mechanism between downstream flow regions and shock fluctuations. The transfer functions are also used as an estimation tool to predict the shock motion accurately; for the largest streamwise separation between input and output signals, correlations above 0.6 are observed between predicted and LES data. Computation of spectral proper orthogonal decomposition (SPOD) modes confirms the existence of upstream traveling waves in the leading spectral mode. Finally, the spectral modes obtained using selected flow regions downstream of the shock enable the reconstruction of a significant portion of the energy in the interaction zone. The current results shed further light on the physical mechanisms driving the shock motion, pointing towards a causal behavior between downstream areas and the characteristic unsteady fluctuations at the approximate shock position.
Two-dimensional, meso-resolved numerical simulations are performed to investigate the complete shock-to-detonation transition (SDT) process in a mixture of liquid nitromethane (NM) and air-filled, circular cavities. The shock-induced initiation behav iors resulting from the cases with neat NM, NM with an array of regularly spaced cavities, and NM with randomly distributed cavities are examined. For the case with randomly distributed cavities, hundreds of cavities are explicitly resolved in the simulations using a diffuse-interface approach to treat two immiscible fluids and GPU-enabled parallel computing. Without invoking any empirically calibrated, phenomenological models, the reaction rate in the simulations is governed by Arrhenius kinetics. For the cases with neat NM, the resulting SDT process features a superdetonation that evolves from a thermal explosion after a delay following the passage of the incident shock wave and eventually catches up with the leading shock front. For the cases wherein mesoscale heterogeneities are explicitly considered, a gradual SDT process is captured. These two distinct initiation behaviors for neat NM and heterogeneous NM mixtures agree with experimental findings. Via examining the global reaction rate of the mixture, a unique time scale characterizing the SDT process, i.e., the overtake time, is measured for each simulation. For an input shock pressure less than approximately $9.4~mathrm{GPa}$, the overtake time resulting from a heterogeneous mixture is shorter than that for neat NM. This sensitizing effect is more pronounced for lower input shock pressures. A random distribution of cavities is found to be more effective in enhancing the SDT process than a regular array of cavities. Statistical analysis on the meso-resolved simulation data provides more insights into the mechanism of energy release underlying the SDT process.
The delayed detonation model describes the observational properties of the majority of type Ia supernovae very well. Using numerical data from a three-dimensional deflagration model for type Ia supernovae, the intermittency of the turbulent velocity field and its implications on the probability of a deflagration-to-detonation (DDT) transition are investigated. From structure functions of the turbulent velocity fluctuations, we determine intermittency parameters based on the log-normal and the log-Poisson models. On the other hand, the analysis of the turbulent velocity fluctuations in the vicinity of the flame front by Roepke suggests a much higher probability of large velocity fluctuations on the grid scale in comparison to the log-normal intermittency model. Following Pan et al., we computed probability density functions for a DDT for the different distributions. Assuming that a DDT can occur in the stirred flame regime, as proposed by Woosley et al., the log-normal model would imply a delayed detonation between 0.7 and 0.8 seconds after the beginning of the deflagration phase for the multi-spot ignition scenario used in the simulation. However, the probability drops to virtually zero if a DDT is further constrained by the requirement that the turbulent velocity fluctuations reach about 500 km/s. Under this condition, delayed detonations are only possible if the distribution of the velocity fluctuations is not log-normal. From our calculations follows that the distribution obtained by Roepke allow for multiple DDTs around 0.8 seconds after ignition at a transition density close to 1x10^7 g/cm^3.
It has long been recognized that the neutrinos detected from the next core-collapse supernova in the Galaxy have the potential to reveal important information about the dynamics of the explosion and the nucleosynthesis conditions as well as allowing us to probe the properties of the neutrino itself. The neutrinos emitted from thermonuclear - type Ia - supernovae also possess the same potential, although these supernovae are dimmer neutrino sources. For the first time, we calculate the time, energy, line of sight, and neutrino-flavor-dependent features of the neutrino signal expected from a three-dimensional delayed-detonation explosion simulation, where a deflagration-to-detonation transition triggers the complete disruption of a near-Chandrasekhar mass carbon-oxygen white dwarf. We also calculate the neutrino flavor evolution along eight lines of sight through the simulation as a function of time and energy using an exact three-flavor transformation code. We identify a characteristic spectral peak at $sim 10$ MeV as a signature of electron captures on copper. This peak is a potentially distinguishing feature of explosion models since it reflects the nucleosynthesis conditions early in the explosion. We simulate the event rates in the Super-K, Hyper-K, JUNO, and DUNE neutrino detectors with the SNOwGLoBES event rate calculation software and also compute the IceCube signal. Hyper-K will be able to detect neutrinos from our model out to a distance of $sim 10$ kpc. At 1 kpc, JUNO, Super-K, and DUNE would register a few events while IceCube and Hyper-K would register several tens of events.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا