Large eddy simulation of turbulent premixed flames using level-set G-equation

Level-set G-equation and stationary flamelet chemistry are used in large eddy simulation of a propane/air premixed turbulent flame stabilized by a bluff body. The aim was to study the interaction between the flame front and turbulent eddies, and in particular to examine the effect of sub-grid scale (SGS) eddies on the wrinkling of the flame surface. The results indicated that the two types of turbulence eddies—the resolved large scale eddies and the unresolved SGS eddies—have different effects on the flame. The fluctuation of the flame surface, which is responsible for the broadening of the time averaged mean flame brush by turbulence, depends on the large resolved turbulence eddies. Time averaged mean flow velocity, temperature, and major species concentrations mainly depend on the large scale resolved eddies. The unresolved SGS eddies contribute to the wrinkling at the SGS level and play an important role in the enhancement of the propagation speed of the resolved flame front. In addition, the spatially filtered intermediate species, such as radicals, and the spatially filtered reaction rates strongly depend on the small SGS eddies. The asymptotic behavior of flame wrinkling by the SGS eddies, with respect to the decrease in filter size and grid size, is investigated further using a simplified level-set equation in a model shear flow. It is shown that to minimize the influence of the SGS eddies, fine grid and filter size may have to be used.     

Investigation of lengthscales, scalar dissipation, and flame orientation in a piloted diffusion flame by LES

This work investigates the structure of a diffusion flame in terms of lengthscales, scalar dissipation, and flame orientation by using large eddy simulation. This has been performed for a turbulent, non-premixed, piloted methane/air jet flame (Flame D) at a Reynolds-number of 22,400. A steady flamelet model, which was represented by artificial neural networks, yields species mass fractions, density, and viscosity as a function of the mixture fraction. This will be shown to suffice to simulate such flames. To allow to examine scalar dissipation, a grid of 1.97 × 10⁶ nodes was applied that resolves more than 75% of the turbulent kinetic energy. The accuracy of the results is assessed by varying the grid-resolution and by comparison to experimental data by Barlow, Frank, Karpetis, Schneider (Sandia, Darmstadt), and others. The numerical procedure solves the filtered, incompressible transport equations for mass, momentum, and mixture fraction. For subgrid closure, an eddy viscosity/diffusivity approach is applied, relying on the dynamic Germano model. Artificial turbulent inflow velocities were generated to feature proper one- and two-point statistics. The results obtained for both the one- and two-point statistics were found in good agreement to the experimental data. The PDF of the flame orientation shows the tilting of the flame fronts towards the centerline. Finally, the steady flamelet approach was found to be sufficient for this type of flame unless slowly reacting species are of interest.     

Numerical simulation of the flow field in a reactor for titanium dioxide production

A physicochemical and fluid dynamic model is formulated for the numerical simulation of the flow field in a reactor for titanium dioxide production, the turbulence motion is described by theKε equation, the governing equations are solved by the SIMPLER algorithm devised by Partankar and Spalding. The velocity, tmperature and concentration fields are obtained for three cases: A) with chemical reaction and thermal insulation on the walls; B) with chemical reaction and wall temperature is 450K; C) without chemical reaction and thermal insulation on the walls, and the physicochemical numerical simulation for the titanium dioxide production has been done. The results of the paper can be used as a theoretical guide for the engineers in the design of such reactors.     

A fractional order hyperchaotic system derived from Liu system and its circuit realization

In this paper we propose a novel four-dimensional fractional order hyperchaotic system derived from Liu system. Electronics workbench (EWB) and Matlab simulations show the dynamical behavior of the proposed four-dimensional fractional order hyperchaotic system. Finally, after separately using EWB and Matlab, an electronic circuit is designed to realize the novel four-dimensional fractional order hyperchaotic system and the experimental circuit results are obtained which are identical to software simulations.     

高质量等曲率物参光像面数字全息显微系统

对等曲率物参光像面数字全息显微成像系统进行研究,分析了光路配置方法,推导了系统的点扩散函数,并由此指出了决定系统成像分辨率的因素及系统的成像特点,最后讨论了再现像光场一次位相畸变校正的方法.结果表明,等曲率物参光像面数字全息具有最大的信息容量;该系统的成像分辨率取决于显微物镜的数值孔径和CCD的像元大小,与CCD的光敏面尺寸无关;物体各点中通过显微物镜的所有频率成分均能被系统完全记录与再现,样品被照亮区域的大小对记录条件和再现像质没有影响;等曲率物参光像面数字全息系统是一种优化的全息记录与再现系统,利用该系统可以实现高质量成像.实验结果验证了理论分析的正确性.     

宽光谱棱镜型太阳光谱仪设计

为实现大气层外太阳光谱辐照度(SSI)变化的长期例行监测,设计了一种星载宽光谱太阳光谱仪结构。全系统仅使用单片折反式曲面棱镜实现太阳光谱250～2500 nm的分光,并通过棱镜转动实现谱平面上多探测器的同步扫描探测;同时基于Huygens子波点扩展函数（PSF）仿真了光谱仪的光谱响应函数（SRF）和光谱分辨率。分光棱镜在±2.5°扫描转角内的全谱段子午像差小于8μm;光谱分辨率在紫外谱段(250～400 nm)为0.7～3.5 nm,可见/近红外谱段(400～1000 nm)为3.5～35.0 nm,短波红外谱段(1000～2500 nm)内为28.5～41.2 nm。整个系统结构简单紧凑,性能稳定可靠,分光和像差校正能力满足大气层外太阳光谱辐照度长期监测需求。     

大气湍流对天文望远镜光电导行精度的影响

讨论了由大气湍流造成的望远镜导行误差,这一误差将给望远镜光电导行的精度带来较大影响。其中使用单点源目标的导行精度受限于大气相干长度和大气相干时间,而多点源目标及面源的导行精度还与湍流随高度的变化有关。通过分析和仿真(数值模拟)结果表明,现代天文观测必须考虑大气湍流对光电导行精度的影响。通过延长曝光时间可以有效降低湍流大气的影响,同时也降低了光电导行的反馈控制频率。当导行信标为多星或面源则可在一定程度上降低高层大气湍流带来的光电导行误差,从而可以适度减少光电导行曝光时间,提高其反馈控制频率。     

二氧化铬阴离子光电子能谱的理论研究

凭借密度泛函理论,采用不同基组对中性分子CrO2的基态（ 3B1）以及阴离子CrO2‾的基态（ 4B 1）进行几何优化和振动频率分析;应用量化计算得到的力常数及结构和光谱参数,基于推得的两维四模Franck-Condon重叠积分的代数表示,对CrO2 ( 3B1)－CrO2‾ ( 4B1) 的光脱附过程进行Franck-Condon分析和光谱模拟,理论上得到光电子能谱的谱线相对强度及振动结构分布,理论谱与实验测得的二氧化铬阴离子光电子能谱达到一致,并对光电子能谱的振动结构进行归属及热带分析;另外,在光谱模拟过程中通过迭代Franck-Condon分析过程,推得CrO2‾（ 4B1）与CrO2（ 3B1）平衡几何结构之差：ΔR(Cr-O)= 0.05Å,Δ(O-Cr-O)=12o.     

D-T中子源n-γ密度测井规律的蒙特卡罗模拟

根据D-T反应中子的能谱和角分布数据,建立了中子源模型;根据石灰岩地层标准刻度井群数据,建立了井模型。采用MCNP程序模拟了井中中子和射线的输运,得到了不同地层密度、不同源距处NaI探测器中的混合能谱和非弹能谱。在混合能谱2.5~4.5 MeV能区开窗,混合射线相对计数随源距的变化曲线显示,源距应选择在20~80 cm,密度与混合射线计数之间呈现非线性关系。在非弹能谱1.0~8.0 MeV能区开窗,非弹射线相对计数随源距的变化数据显示,源距应选择在20~40 cm或80 cm附近,密度与非弹射线计数之间成近似线性关系。     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构