High Temperature Region Research Articles

Numerical modeling and parametric analysis of temperature distribution in soil based on water content variation under radio frequency heating for in-situ thermal remediation.

Radio frequency heating (RFH) is recognized as an efficient and economical method for volatilizing organic pollutants from contaminated soils. Although some numerical simulations have been conducted to predict temperature changes during RFH, models that account for effects of water content (WC) variation on dielectric properties and temperature evolution in soil were rarely reported. To address this, a three-dimensional numerical model integrating electromagnetic-thermal conversion, heat transfer, and mass transfer was developed with coupling dynamic changes in WC and corresponding impacts on soil properties. The effects of soil types and RFH engineering parameters on temperature evolution and energy utilization were evaluated. The model's temperature prediction was validated using experimental data, indicating that temperature decreases outward from hot electrode, with high-temperature regions forming at the top and bottom regions near electrode. Water evaporation in soil causes a temperature plateau with >76.25 % of energy utilized for heating and evaporating water during 60-d RFH. Moreover, higher soil WC could increase electromagnetic-thermal conversion capability, counteracting the negative impact of increased specific heat capacity on temperature rise. Thermophysical properties of soil are crucial in RFH: sand, with lower specific heat capacity and higher thermal conductivity, results in faster heating rates. While, clay, with higher density and specific heat capacity, results in lower heating rates but greater energy utilization efficiency. Low-power RHF and wider well spacing improve the uniformity of temperature distribution but require more time to reach target temperatures. Moreover, input power influences the magnitude and timing of optimal energy utilization efficiency, while well spacing primarily affects timing. The models developed in this study provide essential guidance for RFH-based soil remediation projects.

Reinforcing the Efficiency of Plastic Upgrading through Full-Spectrum Photothermal Effect Integration of Heat Isolator.

Photoreforming of polyethylene terephthalate (PET) to H2 is practically attractive strategy for upgrading waste plastics. The major challenge is to utilize the infrared energy in the solar spectrum to improve the efficiency for photoreforming of PET to H2. Herein, through the ingenious integration of tungsten phosphide nanoparticles and tungsten single atoms (WP/W SAs) with carbon nitride (g-C3N4), the constructed hybrid inherits both the desirable properties and structural merits of the respective building blocks. Specifically, the photothermal effect of WP/W SAs couples with the "heat isolator" role of g-C3N4 due to its low thermal conductivity, thereby forming localized high-temperature regions, reducing the activation energy and improving the kinetics in the photoreforming of PET to H2. Additionally, the green pretreatment of PET using alkali-free hydrothermal strategy is reported, achieving direct separation of the ethylene glycol and terephthalic acid. This work not only provides an alkali-free hydrothermal pretreatment for PET, but also integrates the photothermal effect with the thermal insulation and opens a new avenue for harnessing solar energy into to convert plastics into H2.

Design and performance analysis of a coupled burner for the Solid Oxide Fuel Cell system

The Solid Oxide Fuel Cell (SOFC) system is characterized by high energy conversion efficiency and low emissions. The energy supply and recovery of the SOFC system are facilitated through the ignition burner and the off-gas burner, respectively. This study proposes a coupled design of the ignition burner and the off-gas burner to reduce the burner volume, thereby enhancing the power density of the SOFC system. The ignition burner employs radially opposed jet combustion with fully premixed natural gas, while the off-gas burner uses swirl diffusion combustion with the anode off gas. When the premixed gas and lean anode off gas are unable to sustain the flame, catalytic combustion is initiated. Simulation and experiment methods are utilized to analyze the performance of the coupled burner in this paper. The results indicate that uniform mixing is achieved with a high swirl number (s1 = 2.6) for the premixed gas, and the maximum excess air coefficient for stable flame combustion at 400 °C is 2.47. As the temperature of the cathode off gas increases, the flameout boundary for the premixed gas gradually increases; conversely, with increasing burner power, the flameout boundary gradually decreases. The intersection cooling effect of mixed air and circumferential high-temperature flue gas is effective, and the flame does not exceed the tolerance temperature of the catalytic carrier. The diffusion combustion of anode off gas with a swirl number (s2 = 0.5) exhibits great gas mixing uniformity without a high-temperature core region. When transitioning from flame to catalytic combustion with premixed natural gas and anode off gas, there is a temperature lag of approximately 60 s. The isotropic cylindrical cordierite catalytic carrier with a 400 mesh and a coating of 60 g/ft3 Pt0.1Pd0.5 achieves stable catalytic combustion in an environment with preheated air temperatures ranging from 400 to 600 °C and H2 concentrations of 0.25–2.91 %.

Enhancing artificial permafrost table predictions using integrated climate and ground temperature data: A case study from the Qinghai-Xizang highway

The complexities of permafrost changes, driven by climate warming and engineering activities, coupled with challenges in data acquisition, make it crucial and challenging to accurately predict the artificial permafrost table, particularly for subgrades in high-temperature unstable permafrost regions. To address this, this study developed a hybrid machine learning model (RF-LSTM-XGBoost) for permafrost table prediction. By analyzing climate change and ground temperature data from various positions and depths along the subgrade in the Tuotuo River section of the Qinghai-Xizang Highway, the Spearman correlation coefficient method was initially used to determine the important influencing factors. Random Forest (RF), Long Short-Term Memory Neural Network (LSTM), and Extreme Gradient Boosting (XGBoost) models were used to predict the artificial permafrost table, and grid search and cross-validation methods were employed to optimize the hyperparameters of each model. A linear weighted combination method based on the minimum cumulative absolute error was utilized to merge the models, and its performance was compared with the individual RF, LSTM, and XGBoost models. Subsequently, the feature importance of the variables in the machine learning model was analyzed. The results indicated a strong correlation between artificial permafrost table changes and factors such as daily average atmospheric temperature, subgrade surface ground temperature, and subgrade surface ground heat flux during the freezing-thawing cycle. The combined model highlighted daily atmospheric temperature as the most influential predictor, followed by ground heat flux, with the surface ground temperature being less impactful. The combined model demonstrated improved predictive accuracy, with MSE, MAPE, RMSE, MAE, and R2 values of 0.003, 0.052, 0.0085, 0.029, and 0.989, respectively, surpassing those of individual models. This model offers a rapid, accurate, and reliable approach for permafrost table prediction, advancing subgrade stability research in challenging permafrost environments.

CFD Investigation on Air-Staged Combustion Characteristics of Deep Peak Shaving Wall-Fired Burner Boilers

ABSTRACT In alignment with China’s ambitious goals of reaching peak carbon emissions and achieving carbon neutrality, stringent load peak regulations are being enforced for coal-fired power generation. This underlines the significance of deploying deep load peaking strategies in coal-fired units. The current research focus is on wall-fired boilers operating under ultra-low loads of 600 MWe. It delves into the flow dynamics, combustion characteristics, component concentration field, and NOx emission characteristics under varying air coefficients in the primary combustion zones (α). The findings suggest the formation of a consistently high-temperature region in the primary combustion zone, where temperatures surpass 1500 K. The middle burner exhibits superior combustion stability compared to the bottom burner. As α diminishes, several changes occur. The flow field intensity in the primary combustion zone lessens, the ignition distance contracts, and the O2 concentration in this area drops. This engenders a reducing atmosphere characterized by low O2 and high CO levels, which hampers the formation of NOx. At α = 0.95, the NO concentration peaks, with the primary generation mechanism being temperature-driven. At this juncture, the NO concentration at the furnace flue gas outlet is 287.36 mg/m3.

Optimized emamectin benzoate trunk injection: addressing temperature limitations for pine wilt disease control.

Pine wilt disease (PWD), caused by the pinewood nematode Bursaphelenchus xylophilus, poses a significant threat to global forestry, resulting in extensive economic and ecological damage. Traditional trunk injection agents (TIAs) are limited by environmental factors, such as high winter temperatures and vigorous resin flow, particularly in southern China. This study aimed to develop an ordinary temperature trunk injection agent (OTTIA) for year-round application and to evaluate its performance under high-temperature conditions in comparison to commercial TIAs. Extensive laboratory and field tests identified N,N-Dimethylformamide (DMF) and benzyl acetate as optimal solvents and Tween-40 as an effective emulsifier for the OTTIA formulation. The new agent was fully absorbed by Pinus massoniana within 3 h at temperatures exceeding 30 °C, outperforming existing commercial agents. Treated trees maintained comparable emamectin benzoate (EB) levels but achieved lower LC90 lethal concentration values (27.65 mg L-1), indicating greater efficacy. The OTTIA provided protection against PWD for ≤360 days postinjection. The development of OTTIA marks a significant advancement in PWD management, offering an effective solution for high-temperature regions and potentially transforming year-round treatment strategies. The results highlight the critical role of optimizing solvent and emulsifier combinations to enhance the efficacy and application of trunk injection agents, contributing to sustainable forest management and improved protection against PWD. © 2024 Society of Chemical Industry.

Temperature dependences of the magnetic anisotropy constants of single-crystal inclusions MnSb in an InSb matrix

The contributions of the first order K1 and second order K2 magnetic anisotropy constants to the effective constant are separated. Their competition determines the type of magnetic anisotropy «easy plane». Extrapolation of the dependences K1(T) and K2(T) to the region of high temperatures made it possible to predict the temperature TSR = 570 K, which corresponds to the spin-reorientation transition, at which the easy-axis magnetic anisotropy is formed.

Damage assessment of the CRTS III ballastless slab track in high temperature regions

Temperature-induced interfacial damage between the prefabricated slab and filling layers is a common ailment faced by slab-type ballastless tracks systems. In this study, the damage behaviour of the CRTS III ballastless track in high temperature regions is explored. To this end, a 3-dimensional finite element model of the CRTS III ballastless track is built in Abaqus, which takes into account the non-linear temperature distribution along the track depth, the non-linear cohesive relationship between the concrete to self-compacting concrete interface, and the constitutive stress-strain relation of concrete. Upon subjecting the model to various temperature and wheel loading conditions, the following main conclusions were drawn: (1) Upward warpage of the prefabricated slab undergoes an average percentage increase of 37.62 % with each 10 °C/m increase in temperature gradient, signalling a greater up-warping effect in higher temperature regions; (2) the stresses along the interface’s first and second shear directions have the greatest influence on damage initiation of the interface; (3) under the effect of equally increasing temperature gradient loads, interface nodes where damage initiation had taken place showed a steady increase with an average difference of 5.21 %; (4) applied wheel loads significantly increase the damage initiation area of the track, but the effect of varying concentrated forces typically experienced during the track’s service life does not affect the change in damage initiation area as severely.

Influence of temperature gradients on helium diffusion and clustering in tungsten: A molecular dynamics study

Tungsten and its alloys, used as plasma-facing materials in fusion reactors, endure high-flux, low-energy helium ion irradiation and temperature gradients induced by thermal load and shocks. This study utilizes molecular dynamics (MD) simulations to explore the diffusion and clustering behavior of helium (He) in tungsten (W) under temperature gradients. The migration behavior of isolated helium atoms, small helium clusters, and helium self-interstitial atom (He-SIA) clusters within tungsten is analyzed. It is revealed that both helium and He-SIA clusters tend to migrate towards higher temperature regions, exhibiting negative thermophoresis. The nucleation of helium clusters, the formation of Frenkel pairs, and the interactions between self-interstitial clusters and helium clusters are also investigated. The results indicate that helium concentration significantly impacts He nucleation behavior. Directional diffusion may lead to areas of high concentration over time, potentially promoting the formation of He clusters and bubbles. These insights enhance the understanding of tungsten's performance and durability as a plasma-facing component in fusion reactors.

Effect of transition metals (Fe, Ru and Os) on the elastic, thermodynamic, electronic, magnetic and thermoelectric properties of three Nd-filled arsenide skutterudites based on first-principles computations

Elastic, electronic, magnetic, thermodynamic, and thermoelectric properties of three neodymium-filled skutterudites NdT4As12 (T = Fe, Ru and Os) have been investigated using density functional theory. The bulk moduli increase with the atomic numbers of the transition metals. All materials are stiff, while the Os compound is the stiffest. A considerable ionic contribution is seen in their inter-atomic bonding. The Fe compound possesses the highest Debye temperature, while the Os compound has the highest melting temperature. GGA+U formalism predicts the semimetallic nature of all the compounds, but their total DOS differs near the Fermi level in both spin channels. The contributions of the Nd-f and T-d electrons differ in the three compounds and two methods. GGA+U predicts the highest magnetic moment in the Fe compound, followed by the Ru and Os compounds. The thermoelectric figure of merit shows no systematic variation with z of the transition metals and in the high temperature region it increases with temperature pretty fast. The highest value of ZT is 0.32 in the Ru compound in the spin-down channel while it is 0.24 in the spin-up channel in the Fe compound, both at 1000 K.

Continuous Casting Slab Mold: Key Role of Nozzle Immersion Depth.

Based on a physical water model with a scaling factor of 0.5 and a coupled flow-heat transfer-solidification numerical model, this study investigates the influence of the submerged entry nozzle (SEN) depth on the mold surface behavior, slag entrapment, internal flow field, temperature distribution, and initial solidification behavior in slab casting. The results indicate that when the SEN depth is too shallow (80 mm), the slag layer on the narrow face is thin, leading to slag entrapment. Within a certain range of SEN depths (less than 170 mm), increasing the SEN depth reduces the impact on the mold walls, shortening the "plateau period" of stagnated growth on the narrow face shell. This allows the upper recirculation flow to develop more fully, resulting in an increase in the surface flow velocity and an expansion in the high-temperature region near the meniscus, which promotes uniform slag melting but also heightens the risk of slag entrainment due to shear stress at the liquid surface (with 110 mm being the most stable condition). As the SEN depth continues to increase, the surface flow velocity gradually decreases, and the maximum fluctuation in the liquid surface diminishes, while the full development of the upper recirculation zone leads to a higher and more uniform meniscus temperature. This suggests that in practical production, it is advisable to avoid this critical SEN depth. Instead, the immersion depth should be controlled at a slightly shallower position (around 110 mm) or a deeper position (around 190 mm).

Thermodynamic behavior of high-power inductively coupled plasma quartz tube wall

High-power inductively coupled plasmas are commonly used in planetary entry simulations and are increasingly being used in electric propulsion applications. However, during the operation of the system, the walls of the quartz tube will crack and melt. Its thermodynamic behavior is key to ensuring the safe and reliable operation of the system, which is directly related to the distribution of thermal energy within the discharge volume. In this paper, the temperature and stress distribution of the quartz tube wall of an inductively coupled plasma generator at 27 kW–85 kW are described. A numerical simulation model was established to depict the interaction between the plasma and the quartz tube wall. In the field of experimental research, the temperature of the outer wall of the quartz tube was obtained by using a thermal imager, and a non-uniform B-spline difference method was proposed to fit the outer wall temperature of the quartz tube to eliminate the influence of the induction coil. It is found that the numerical simulation and experimental results show that the temperature is stable region, temperature rise area, temperature drop zone, and the high temperature region of the quartz tube wall is located in the coil area, and the high stress area is also located in this region. On this basis, the outer wall temperature and thermal stress of quartz tubes under different heat fluxes are studied. When the heat flux exceeds 18.6 kW/m2, the stresses in the coil area and downstream of the coil exceed the limit stress. Mechanical failures may occur in areas where the ultimate stresses are exceeded, and these results can provide theoretical data for the optimal design of high-power inductively coupled plasma generators.

Interparticle interaction effect on superparamagnetic properties of ferrimagnetic particles

Uncoated and oleic acid coated ferrimagnetic particles are synthesized. Both samples contain similar magnetite particles. Structural characterization of these particles is performed using x-ray diffraction, Fourier transform infrared spectra, dynamic light scattering, transmission electron microscopy and thermogravimetric analysis. Average crystallite size of magnetite is 7 nm. Thickness of oleic acid layer on each magnetite particle is 1 nm. The oleic acid coated sample contains 23% oleic acid. Temperature dependence of the zero field cooled and field cooled susceptibility of the uncoated magnetite particles peaks at 120 K and both the curves bifurcate at 195 K. The peak temperature and the bifurcation temperature of the system decrease due to the oleic acid coating. The inverse susceptibility does not vary linearly with temperature in high temperature region. Magnetization versus applied magnetic field data at 250 K are analyzed in three different ways to study the interparticle interaction and the particle size distribution effects on magnetization. We find that the both factors affect the magnetization process of ferrimagnetic nanoparticles.

Improved thermoelectric properties of n-type ZrNiCu0.05Sn by doping Y2O3

Half-Heusler (HH)-based alloys are promising candidates for high-temperature thermoelectric (TE) energy conversion materials due to their excellent adaptability in high-temperature regions (900 - 1100 K). However, their TE performance is greatly limited owing to the inherently high lattice thermal conductivity (κl). In this study, we introduce a synergistic strategy to reduce thermal conductivity while modifying electrical transport properties through yttrium oxide (Y2O3) dispersion. The binary NiSn compound precipitations with sub-micron dimension together with Y2O3 nanoparticles synergistically intensify the phonon scattering and enhance the seebeck coefficient (S) effectively via energy filtering effect (EFE) at hetero-interfaces, leading to a remarkable reduction in κl and a slight improvement in the power factor (PF). Consequently, upon 1.5 wt% Y2O3 introduction, the sample incorporating Y2O3 exhibits a ∼13% reduction in total thermal conductivity (κt) as compared to the HH matrix, reaching a peak ZT of 0.86 at 938 K, thereby successfully optimizing its TE performance. The presented strategy provides a feasible approach to simultaneously modulate the thermal and electrical transport behaviors to improve the TE performance of HH-based materials.

Heat transfer performance of regenerative cooling in the presence of supersonic combustion based on two-way decoupling method

As a promising active cooling method for aero engines, regenerative cooling with endothermic hydrocarbon fuel plays an important role in thermal protection systems. However, simply assuming the thermal boundary conditions for the numerical simulation of regenerative cooling as either Dirichlet or Neumann boundary conditions is often unsuitable, as the combustion process exhibits complex three-dimensional (3D) characteristics. We propose a novel multi-step iterative calculation method in this work to simulate the coupled flow and heat transfer between supersonic combustion in the combustor and regenerative cooling outside. The method captures the main features of the supersonic combustion including different shockwave structures, and generates a realistic 3D distribution of heat flux as the boundary condition for the regenerative cooling study. With improved simulation, the result shows that regenerative cooling significantly reduces the peak temperature of the combustor by 58.5% with a mass flow rate of 1 g/s. A local high-temperature region is observed on the coupled wall inside the combustor near the injector, where the flame is primarily concentrated. Increasing the mass flow rate by 100%, from 1 g/s to 2 g/s, leads to a reduction in the maximum wall temperature by up to 15.7%.

Structural modulation in colored polymer bands of methacrylic acid-based frontal polymerization

The periodic colored polymer bands have been investigated in frontal polymerization (FP) reaction containing methacrylic acid (MAA), hydroquinone, Benzoyl peroxide (BP), and N, N Dimethyl aniline (DMA) system. The MAA acts as a monomer, which was diluted adequately with methanol (MeOH), ethanol (EtOH), and butanol (BuOH) separately to observe the patterning characteristics in different reaction schemes. The structural modulations that include the periodicity of colored polymer bands, morphological changes, and spacings between two adjacent layers were studied. The colored variations of polymer bands and their dependency on the concentration of BP, hydroquinone and alcoholic compositions have been studied. The highly consistent pink-white colored polymer bands, loosely-packed band structures, and perfectly ordered band structures resulted in MeOH, EtOH, and BuOH systems, respectively, as discussed. The evolution of tiny bubbles and convection-type instability has been observed in different conditions of the FP reactions that significantly affect the planar movement of the polymer fronts. The hot spots, which usually represent the high-temperature regions of a typical frontal surface, have also been demonstrated in all reacting systems, which may cause spin mode propagation of the fronts and change the surface geometry, as described. The materials characterization was carried out using UV-visible spectrophotometer, NMR, FTIR, and FESEM techniques, providing information about the polymer phases involved in band structures, composition, and surface properties. The analytical data and results obtained during the study further emphasized that two different coloured polymers, namely 1, 4-dihydroxy anthraquinone methacrylic acid and poly-methacrylic acid dihydroxy, anthraquinone produced simultaneously and crystallized periodically, results in the development of a coloured polymer band structures. The possible reactions and associated chemical mechanisms concerning the observations, the nature of the monomer, and solvent characteristics have been proposed, which show a close agreement with the analytical data and results obtained during the study

Numerical Simulation of Spectral Radiation for Hypersonic Vehicles

The spectral radiation of hypersonic vehicles and surrounding flow fields is crucial for optical target detection. A novel approach combines Low-Discrepancy Sequences with the Reverse Monte Carlo Method to simulate the spectral radiation of hypersonic missiles. The effects of chemical non-equilibrium reactions and high-emissivity coating failures on the spectral radiation of a typical biconical hypersonic missile were investigated. Significant differences were found among chemical equilibrium, non-equilibrium, non-reactive, and catalytic wall models. High-emissivity coating failures occur mainly in the high-temperature regions of the nose cone and fins. The non-equilibrium model shows a transition peak distribution of NO in the head shock layer within the ultraviolet gamma band, with integrated ultraviolet radiation (200–300 nm) three times higher than the equilibrium model. In the 1–3 μm band, the non-equilibrium model’s radiation intensity at a 120° horizontal detection angle is about 1.46 times that of the equilibrium model. Using real coating emissivity, the 3–5 μm band radiation intensity is about 5% higher, and the 8–14 μm band is about 8.51% lower than the uniform emissivity model. When high-emissivity coating emissivity fails by 40%, the 3–5 μm band intensity decreases by about 15.38%, and the 8–14 μm band intensity decreases by about 12.67%.

Numerical investigation on the influence of swirl number to high-temperature zone evolution and outlet temperature distribution in multi-stage combustor

The outlet temperature distribution plays a crucial role in determining the lifespan and stability of downstream thermal components. It poses a significant challenge in the design and regulation of high-temperature combustors. This study aims to investigate the influence of swirl number on the outlet temperature distribution in multi-stage combustor. The impact of varying fourth-stage swirl numbers on flow structure, hot streaks, and temperature distribution is compared and analyzed. The results indicate that increasing the swirl number significantly improves the air–fuel mixing within the combustion zone, leading to the expansion of high-temperature regions. Additionally, the high-temperature region migrates upstream along the primary zone. With an increase in the swirl number, the relative area of hot spots within the combustor also increases. Remarkably, the swirler configuration with a swirl number of 1.5 exhibits the most minimal exit temperature gradient and the highest level of temperature distribution uniformity, quantified by an outlet temperature distribution factor (OTDF) of 0.26. Furthermore, changes in the swirl number affect the distribution of combustion modes within the combustor. Increasing the swirl number promotes the premixed combustion mode at the pilot stage exit and enhances temperature uniformity at the swirler exit. Thorough comprehension of the temperature distribution in the primary zone enables control of the outlet temperature distribution from its origin. Uniformity amplitude (UA) of HCO and primary temperature distribution factor (PTDF) are defined to describe the temperature distribution in primary zone. At a swirl number of 1.5, the UA and PTDF also shows the lowest values of 0.2754 and 0.3019, respectively. Overall, when improving the distribution of HCO under both premixed and non-premixed modes, a more uniform temperature distribution across the primary zone can be achieved, effectively optimize the exit temperature distribution.

Liquid state properties and amorphous solidification kinetics of multicomponent Fe_{50-x}Co_{x}Cr_{14}Mo_{14}C_{9}B_{8}Tm_{5} alloys investigated under containerless processing conditions.

The liquid state thermophysical properties and amorphous solidification kinetics of Fe_{50-x}Co_{x}Cr_{14}Mo_{14}C_{9}B_{8}Tm_{5} (x=10, 15, 20, and 25) alloys were explored by electromagnetic and electrostatic levitation techniques. It was found that the surface tension of liquid alloys with Fe contents below 30 at.% had a strong temperature dependence. The high surface tension led to a sharp increase in the interfacial free-energy penalty. A high nucleation barrier was formed inside the melt, which greatly inhibited the nucleation rate. The liquid viscosity revealed the strong liquid feature of this alloy series, which became the dominant kinetic factor determining their solidification mechanisms. The high viscosity hindered atomic diffusion and delayed nucleation. The long crystallizing incubation time of Fe_{30}Co_{20}Cr_{14}Mo_{14}C_{9}B_{8}Tm_{5} ensured a strong glass-forming ability and a low critical cooling rate of 2.11×10^{3}Ks^{-1}. As a result, a bulk metallic glass rod with an 8-mm diameter was successfully prepared by a convenient casting procedure. This rod could remain in a glassy state at a higher temperature and over a wider temperature range due to its high glass-transition temperature and large undercooled liquid region. The apparent activation energy for nucleation was derived as 463.8 kJmol^{-1} according to the nonisothermal crystallization kinetics, indicating that the bulk metallic glass had to absorb a large amount of energy to overcome the potential barriers before nucleation and thus exhibited excellent thermal stability.

The effect of cavity shape on critical high-temperature regions formation pattern and structural evolution difference in explosive particle bed

The occurrence of critical high-temperature regions within explosives is a crucial factor in initiating combustion and explosion. The primary causes of hot regions in cavity-containing explosives are dissipative processes occurring inside the explosive or at its interface, such as shear, friction, viscosity, plasticity, and heat transfer. Currently, there is limited research on the generation and dissipation mechanisms of critical high-temperature regions in cavity-containing explosives during impacts, based on the mesoscale analysis that describes particle motion. In this study, a model based on the discrete element method is employed allowing for an accurate analysis of the dynamic and thermodynamic processes among particles, which systematically considers the shear history between particles, elastoplastic collisions, and heat transfer between particles. The temporal evolution of particle temperature and dissipated work contours during the impact process of explosives with different cavity shapes is analyzed. Explosives with cavities undergo three stages during the impact process: the formation of the trapezoidal high-temperature region, cavity collapse, and dispersion of the particle bed, ultimately leading to the formation of two main critical high-temperature regions: near the cavity and near the wall. The temperature rise of hot particles near the wall can be roughly divided into two stages, and the temperature rise mechanisms are different. Initially, the temperature rise is primarily attributed to plastic dissipation of hot particles, and the same initial impact velocity results in the same temperature rise. During the cavity collapsing, continuous momentum transfer occurs between the hot particles near the wall and the particles in the high-temperature shear bands, deeply influencing the development of the high-temperature shear bands through the movement of the hot particles near the cavity. Therefore, the different collapse modes of hot particles near cavities of different shapes lead to differences in the later temperature evolution of the hot particles near the wall. The unique aggregation processes of hot particles near cavities of various shapes result in distinct high-temperature region morphologies: circular cavity exhibits approximate “——” shape, triangular cavity displays arch shape, and inverted triangular cavity presents “M” shape. Significantly, a new mechanism resulting from the interaction of two opposing “jet “flows for the formation of higher-temperature hot regions near the circular cavity has been discovered.