A hybrid level set method for modelling detonation and combustion problems in complex geometries

We present an accurate and fast wave tracking method that uses parametric representations of tracked fronts, combined with modifications of level set methods that use narrow bands. Our strategy generates accurate computations of the front curvature and other geometric properties of the front. We introduce data structures that can store discrete representations of the location of the moving fronts and boundaries, as well as the corresponding level set fields, that are designed to reduce computational overhead and memory storage. We present an algorithm we call stack sweeping to efficiently sort and store data that is used to represent orientable fronts. Our implementation features two reciprocal procedures, a forward ‘front parameterization’ that constructs a parameterization of a front given a level set field and a backward ‘field construction’ that constructs an approximation of the signed normal distance to the front, given a parameterized representation of the front. These reciprocal procedures are used to achieve and maintain high spatial accuracy. Close to the front, precise computation of the normal distance is carried out by requiring that displacement vectors from grid points to the front be along a normal direction. For front curves in two dimensions, a cubic interpolation scheme is used, and G ¹ surface parameterization based on triangular patches is used for the three-dimensional implementation to compute the distances from grid points near the front. To demonstrate this new, high accuracy method we present validations and show examples of combustion-like applications that include detonation shock dynamics, material interface motions in a compressible multi-material simulation and the Stephan problem associated with dendrite solidification.     

Deprotection‐free preparation of propargyl ether‐containing phosphinated benzoxazine and the structure–property relationship of the resulting thermosets

Generally, protection and deprotection procedures of amino groups are required in preparing propargyl ether‐containing benzoxazines. In this study, we report a facile, deprotection‐free preparation of a propargyl ether‐containing phosphinated benzoxazine (2) from the nucleophilic substitution of a phenolic OH‐containing phosphinated benzoxazine (1) and propargyl bromide in the catalysis of potassium carbonate. The structure of (2) was characterized and confirmed by a high‐resolution mass spectrum, ¹H, ¹³C, ¹H‐¹H, ¹H‐¹³C nuclear magnetic resonance (NMR) spectra, and X‐ray single crystal diffractogram. infrared (IR) and differential scanning calorimetry were used to monitor the ring‐opening of benzoxazine and crosslinking of propargyl ether. The microstructure and the structure–property relationship of the resulting homopolymers and copolymers are discussed. The T_g of homopolymer of (2) is 208 °C by dynamic mechanical analysis, the coefficient of thermal expansion is 43 ppm/°C, and T_d 5% (N₂) is 393 °C, respectively, which are higher than those of the homopolymer of (1) . Similar trends were observed in the copolymerization system. The results demonstrate the beneficial effect of crosslinking afforded by the propargyl ether group is higher than that by the phenolic OH group. © 2011 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem, 2012     

Modelling effects of subgrid-scale mixture fraction variance in LES of a piloted diffusion flame

A posteriori analysis of the statistics of two large-eddy simulation (LES) solutions describing a piloted methane–air (Sandia D) flame is performed on a series of grids with progressively increased resolution reaching about 10.5 million cells. Chemical compositions, density and temperature fields are modelled with a steady flamelet approach and parametrised by the mixture fraction. The difference between the LES solutions arises from a different numerical treatment of the subgrid scale (SGS) mixture fraction variance – an important quantity of interest in non-premixed combustion modelling. In the first case (model I), the variance transport equation is solved directly, while in the second (model II), an equation for the square of the mixture fraction is solved, and the variance is computed from its definition. The comparison of the LES solutions is based on the convergence properties of their statistics with respect to the turbulence resolution length scale. The dependence of the LES statistics is analysed for velocity and the mixture fraction fields, and tested for convergence. For the most part, the statistics converge for the finest grids, but the variance of the mixture fraction shows some residual grid dependence in the high-gradient regions of the jet near field. The SGS variance given by model I exhibits realisability everywhere, whereas in regions of the flame model II is non-realisable, predicting negative variances. Furthermore, the LES statistics of model I exhibit superior convergence behaviour.     

Solid support flame synthesis of 1-D and 3-D tungsten-oxide nanostructures

In this paper we report the growth of 1-D and 3-D tungsten-oxide nanostructures on tungsten wire probes inserted in an opposed-flow oxy-fuel flame. The probe diameter and oxygen content in the oxidizer were varied to study their influence on the growth of tungsten-oxide nanostructures. The introduction of a 1-mm diameter W probe into the flame environment with an oxidizer composition of 50%O₂ + 50%N₂, resulted in the formation of 1-D nanorods on the upper surface of the probe. The formation of triangular, rectangular, square, and cylindrical 3-D channels with completely hollow or semi-hollow morphology was achieved by reducing the probe diameter to 0.5 mm. Whereas, the increase of the O₂ content to 100% and the employment of a 1-mm probe resulted in the growth of ribbon-like micron-sized structures. The lattice spacing of ∼0.38 nm measured for the 1-D W-oxides closely matches a monoclinic WO₃ structure. X-ray photoelectron spectroscopy analysis revealed that the larger 3-D structures also consist of WO₃ confirming that the chemical composition of the structures remains the same while varying the probe and flame parameters. The proposed growth mechanism states that the 3-D WO₃ structures are formed through the lateral coalescence of 1-D W-oxide nanorods.     

Improving the Flame Retardancy of Polybenzoxazines with a Reactive Phosphorus-Containing Compound

A novel flame retardant resin was prepared by modifying benzoxazines with a reactive phosphorus-containing compound, 2-(6-Oxido-6H-dibenzo[c,e][1,2]oxa-phosphorin-6-yl) 1,4-benzenediol (ODOPB). The curing reactions of the benzoxazines containing ODOPB (ODOPB-BOZ) and a type of bifunctional benzoxazine (B-BOZ) were investigated via differential scanning calorimetry. The ODOPB-BOZ precursors exhibited lower curing temperature in comparison with the B-BOZ. Through dynamic mechanical analysis tests, it was found that the glass transition temperature and storage modulus of the cured ODOPB-BOZ were lower than those of the cured B-BOZ. Moreover, the tensile and flexural properties of the cured ODOPB-BOZ were also lower than those of the cured B-BOZ. The TGA tests proved that, owing to the presence of ODOPB, the char yield of the polybenzoxazines were remarkably raised while the decomposition rates at higher temperature were reduced. It was found that the flame retardancy of cured ODOPB-BOZ attained V-0 grade upon the UL-94 tests while that of cured B-BOZ was V-1 grade. SEM images showed that more compact char layers were formed in the cured ODOPB-BOZ after combustion, resulting in the improvement of flame retardance.     

“蜡烛火焰划分”实验的局限成因分析及改进研究

人教版初中化学首个探究性实验"蜡烛及其燃烧"明确指出火焰的划分方法，但在实际教学中肉眼观察法和火柴梗燃烧法对蜡烛火焰的划分存在一定的局限性。采用文献研究法和实验验证法分析其局限性的成因——燃料的析炭能力强弱、材料结构的不同等，并提出相应的改进建议。     

CH4/空气简化动力学机理数值验证

选择绕圆柱预混燃烧算例，验证CH₄/空气三种简化动力学机理（16s41r、15s19r和53s325r）.考虑均匀来流，忽略湍流和湍流与燃烧相互作用以及燃料扩散效应，假设层流有限反应速率，采用保自由流5阶WENO格式求解多组分Euler方程组，得到CH₄/空气预混燃烧流场温度等值线、沿驻点线压力和温度及其CH₄、CO和CO₂质量百分数分布.结果表明：三种简化动力学机理给出的流场均出现弓形激波和火焰面，弓形激波和火焰驻点距离及其形状、诱导区宽度和简化动力学机理相关.当圆柱直径增大，弓形激波和火焰向圆柱上游移动，对应的驻点距离均增大，诱导区宽度变短，点火延时变小，但火焰和弓形激波位置次序未变化.53s325r模型要比16s41r模型和15s19r模型精度要高，点火延时覆盖的压力和温度范围也较宽，所有简化机理均未完全反应，在较大圆柱直径下游达到化学平衡.