非马尔科夫耗散系统长时演化下的极限环振荡现象

本文研究结构化环境中非马尔科夫耗散系统在长时演化下可能出现的极限环振荡现象. 对于欧姆型谱密度环境中的二能级系统, 由于体系只允许一个束缚态模, 给定初态系统在Bloch空间的长时演化将收敛于一个极限环. 研究揭示了极限环半径与环心位置同环境谱密度函数间的关系. 对于多带光子晶体环境中的二能级系统, 由于其可以存在多个束缚态, 研究展现了系统在长时演化下可能出现的收敛于环面或周期或准周期的振荡行为. 有关环面的特征量与环境谱密度间的量化关系同样得以刻画. 论文随后讨论了两比特系统关联量在局域非马尔科夫耗散环境中长时演化可能出现的特征行为.     

LES/PDF for premixed combustion in the DNS limit

We investigated the behaviour of the composition probability density function (PDF) model equations used in a large-eddy simulation (LES) of turbulent combustion in the direct numerical simulation (DNS) limit; that is, in the limit of the LES resolution length scale Δ (and the numerical mesh spacing h) being small compared to the smallest flow length scale, so that the resolution is sufficient to perform a DNS. The correct behaviour of a PDF model in the DNS limit is that the resolved composition fields satisfy the DNS equations, and there are no residual fluctuations (i.e. the PDF is everywhere a delta function). In the DNS limit, the treatment of molecular diffusion in the PDF equations is crucial, and both the ‘random-walk’ and ‘mean-drift’ models for molecular diffusion are investigated. Two test cases are considered, both of premixed laminar flames (of thickness δ_L). We examine the solutions of the model PDF equations for these test cases as functions of Δ/δ_L and h/δ_L. Each of the two PDF models has advantages and disadvantages. The mean-drift model behaves correctly in the DNS limit, but it is more difficult to implement and computationally more expensive. The random-walk model does not have the correct behaviour in the DNS limit in that it produces non-zero residual fluctuations. However, if the specified mixing rate Ω normalised by the reaction timescale τ_c is sufficiently large (Ωτ_c ? 1), then the residual fluctuations are less than 10% and the observed flame speed and thickness are close to their laminar values. Away from the DNS limit (i.e. h/δ_L ? 1), the observed flame thickness scales with the mesh spacing h, and the flame speed scales with Ωh. For this case it is possible to construct a non-general specification of the mixing rate Ω such that the flame speed matches the laminar flame speed.     

Differences between PREMIER combustion in a natural gas spark-ignition engine and knocking with pressure oscillations

PREMIER (PREmixed Mixture Ignition in the End-gas Region) combustion occurs with auto-ignition in the end-gas region when the main combustion flame propagation is nearly finished. Auto-ignition is triggered by the increases in pressure and temperature induced by the main combustion flame. Similarly to engine knocking, heat is released in two stages when engines undergo this type of combustion. This pattern of heat release does not occur during normal combustion. However, engine knocking induces pressure oscillations that cause fatal damage to engines, whereas PREMIER combustion does not. The purpose of this study was to elucidate PREMIER combustion in natural gas spark-ignition engines, and differentiate the causes of knocking and PREMIER combustion. We applied combustion visualization and in-cylinder pressure analysis using a compression–expansion machine (CEM) to investigate the auto-ignition characteristics in the end-gas region of a natural gas spark-ignition engine. We occasionally observed knocking accompanied by pressure oscillations under the spark timings and initial gas conditions used to generate PREMIER combustion. No pressure oscillations were observed during normal and PREMIER combustion. Auto-ignition in the end-gas region was found to induce a secondary increase in pressure before the combustion flame reached the cylinder wall, during both knocking and PREMIER combustion. The auto-ignited flame area spread faster during knocking than during PREMIER combustion. This caused a sudden pressure difference and imbalance between the flame propagation region and the end-gas region, followed by a pressure oscillation.     

Langford系统Hopf分岔极限环幅值控制

系统地研究了Langford 系统Hopf分岔极限环幅值非线性反馈控制问题.根据中心流形定理和规范型降维理论, 推导含控制增益项的非线性曲率系数控制公式和幅值近似解,通过数值仿真验证并绘制极限环幅控关系曲线. 所推导的控制公式为Langford 系统极限环幅值控制提供了方便有效的控制方法.     

Amplitude control of limit cycle in a van der Pol Duffing system

This paper applies washout filter technology to amplitude control of limit cycles emerging from Hopf bifurcation of the van der Pol--Duffing system. The controlling parameters for the appearance of Hopf bifurcation are given by the Routh--Hurwitz criteria. Noticeably, numerical simulation indicates that the controllers control the amplitude of limit cycles not only of the weakly nonlinear van der Pol--Duffing system but also of the strongly nonlinear van der Pol--Duffing system. In particular, the emergence of Hopf bifurcation can be controlled by a suitable choice of controlling parameters. Gain-amplitude curves of controlled systems are also drawn.     

基于光纤气体检测技术的煤矿自然发火预测预报系统

 分析了近红外光谱吸收式气体检测的基本原理以及煤矿采空区自然发火特征气体的选择,通过使用近红外激光光谱吸收式光纤气体传感器对采空区火灾特征气体浓度检测,报道实现了一套基于多种气体参数的火灾预测预报系统。系统使用分布反馈式半导体激光器（DFB LD）作为光源,光源驱动为低频锯齿波信号叠加高频正弦波信号,通过锁相放大器实现信号提取。根据采空区三带静态分布理论,给出了判断采空区自然发火危险程度的判定方法。该系统全部使用光纤进行井下信号传输和气体检测,使得检测精度更高,不带电工作,实时性更强,系统更加稳定可靠。     

Asymptotic properties of coupled nonlinear langevin equations in the limit of weak noise. II: Transition to a limit cycle

We apply the singular perturbation technique, developed in the companion paper, to the study of the fluctuations at the onset of a limit cycle, both for the cases of a soft and a hard transition. The technique and results are illustrated on the Poincaré model (soft transition) and on the Van der Pol oscillator (hard transition).     

Validation of flow simulation and gas combustion sub-models for the CFD-based prediction of NOx formation in biomass grate furnaces

While reasonably accurate in simulating gas phase combustion in biomass grate furnaces, CFD tools based on simple turbulence–chemistry interaction models and global reaction mechanisms have been shown to lack in reliability regarding the prediction of NO_x formation. Coupling detailed NO_x reaction kinetics with advanced turbulence–chemistry interaction models is a promising alternative, yet computationally inefficient for engineering purposes. In the present work, a model is proposed to overcome these difficulties. The model is based on the Realizable k–? model for turbulence, Eddy Dissipation Concept for turbulence–chemistry interaction and the HK97 reaction mechanism. The assessment of the sub-models in terms of accuracy and computational effort was carried out on three laboratory-scale turbulent jet flames in comparison with the experimental data. Without taking NO_x formation into account, the accuracy of turbulence modelling and turbulence–chemistry interaction modelling was systematically examined on Sandia Flame D and Sandia CO/H₂/N₂ Flame B to support the choice of the associated models. As revealed by the Large Eddy Simulations of the former flame, the shortcomings of turbulence modelling by the Reynolds averaged Navier–Stokes (RANS) approach considerably influence the prediction of the mixing-dominated combustion process. This reduced the sensitivity of the RANS results to the variations of turbulence–chemistry interaction models and combustion kinetics. Issues related to the NO_x formation with a focus on fuel bound nitrogen sources were investigated on a NH₃-doped syngas flame. The experimentally observed trend in NO_x yield from NH₃ was correctly reproduced by HK97, whereas the replacement of its combustion subset by that of a detailed reaction scheme led to a more accurate agreement, but at increased computational costs. Moreover, based on results of simulations with HK97, the main features of the local course of the NO_x formation processes were identified by a detailed analysis of the interactions between the nitrogen chemistry and the underlying flow field.     

Evaluation of a flamelet/progress variable approach for pulverized coal combustion in a turbulent mixing layer

A steady flamelet/progress variable (FPV) approach for pulverized coal flames is employed to simulate coal particle burning in a turbulent shear and mixing layer. The configuration consists of a carrier-gas stream of air laden with coal particles that mixes with an oxidizer stream of hot products from lean combustion. Carrier-phase DNS (CP-DNS) are performed, where the turbulent flow field is fully resolved, whereas the coal is represented by Lagrangian point particles. CP-DNS with direct chemistry integration is performed first and provides state-of-the-art validation data for FPV modeling. In a second step the control variables for FPV are extracted from the CP-DNS and used to test if the tabulated manifold can correctly describe the reacting flow (a priorianalysis). Finally a fully coupled a posteriori FPV simulation is performed, where only the FPV control variables are transported, and the chemical state is retrieved from the table and fed back to the flow solver. The a priori results show that the FPV approach is suitable for modeling the complex reacting multiphase flow considered here. The a posteriori data is similarly in good agreement with the reference CP-DNS, although stronger deviations than a priori can be observed. These discrepancies mainly appear in the upper flame (of the present DNS), where premixing and highly unsteady extinction and re-ignition effects play a role, which are difficult to capture by steady non-premixed FPV modeling. However, the present FPV model accurately captures the lower, more stable flame that burns in non-premixed mode.     

Application and theoretical analysis of the flamelet model for supersonic turbulent combustion flows in the scramjet engine

In the framework of Reynolds-averaged Navier–Stokes simulation, supersonic turbulent combustion flows at the German Aerospace Centre (DLR) combustor and Japan Aerospace Exploration Agency (JAXA) integrated scramjet engine are numerically simulated using the flamelet model. Based on the DLR combustor case, theoretical analysis and numerical experiments conclude that: the finite rate model only implicitly considers the large-scale turbulent effect and, due to the lack of the small-scale non-equilibrium effect, it would overshoot the peak temperature compared to the flamelet model in general. Furthermore, high-Mach-number compressibility affects the flamelet model mainly through two ways: the spatial pressure variation and the static enthalpy variation due to the kinetic energy. In the flamelet library, the mass fractions of the intermediate species, e.g. OH, are more sensible to the above two effects than the main species such as H₂O. Additionally, in the combustion flowfield where the pressure is larger than the value adopted in the generation of the flamelet library or the conversion from the static enthalpy to the kinetic energy occurs, the temperature obtained by the flamelet model without taking compressibility effects into account would be undershot, and vice versa. The static enthalpy variation effect has only little influence on the temperature simulation of the flamelet model, while the effect of the spatial pressure variation may cause relatively large errors. From the JAXA case, it is found that the flamelet model cannot in general be used for an integrated scramjet engine. The existence of the inlet together with the transverse injection scheme could cause large spatial variations of pressure, so the pressure value adopted for the generation of a flamelet library should be fine-tuned according to a pre-simulation of pure mixing.     

Asymptotic behavior of fluctuations for the 1D Ising model in zero-temperature limit

Fluctuation of the average spin for one-dimensional Ising spins with nearest neighbor interactions are studied. The distribution function for the average spin is calculated for a finite volume, finite temperature, and finite magnetic field. As the volume increases and the temperature diminishes at zero magnetic field, there are two limits in which the probability distribution shows quite different behaviors: in the thermodynamic limit as the volume goes to infinity for finite temperature, small deviations of the fluctuations are described by a Gaussian distribution, and in the limit as the temperature vanishes for a finite volume, the ground states are realized with probability one. The crossover between these limits is analyzed via a ratio of the correlation length to the volume. The helix-coil transition in a polypeptide is discussed as an application.