热采耦合数学模型定量化研究及数值模拟

针对稠油油藏开发过程的薄弱环节,提出了油藏热质传递过程的黑箱分析法和白箱分析法,建立了包含蒸汽相的稠油热采耦合黑箱模型和白箱模型,使得储层温度场模型得以完整化和系统化。应用CMG软件进行数值模拟,验证了所建数学模型的合理性,并得到了孔隙度对油藏温度的影响关系。最后,对油藏热流体的指进现象进行了机理分析,并总结了热流体指进现象的产生范围。     

超临界流体Rayleigh-Bénard对流的数值模拟研究

物质气液临界点附近热物理性质发生剧烈变化,会出现一种对热力学平衡有显著加速作用的热声活塞效应。而在长时间尺度上,因重力作用而产生的Rayleigh-Bénard对流在活塞效应的影响下,其表现出来的物理特性与普通流体相比存在较大的差异。我们通过SIMPLE方法对超临界氮在不同临界距离下的自然对流发生过程进行数值模拟,结果显示当流体热力学状态接近临界点时,对流作用的发生取决于边界层内热羽流的形成,并且具有明显的湍流特征,而随着离开临界点的距离加大,流场的形成逐渐过渡到一般可压缩流体的情形。     

氢-空气绕流无限长斜劈斜爆震波数值模拟

采用二维多组分有化学反应的Euler方程,考虑单步化学反应,采用有限反应速率模型,对斜爆震波进行计算,通过将数值模拟结果与实验结果对比,验证了数值模拟方法的有效性。在此基础上,对氢空气绕流无限长斜劈的斜爆震过程进行了数值模拟,分析了斜劈角度对斜爆震波特性的影响。结果表明:当斜劈角度增大时,爆震波中的化学反应速率随之增大,相应的斜爆震波的角度也增加,波后的温度更高,同时波后法向马赫数减小。     

电弧放电等离子体诱导激波的计算

基于电弧放电物理过程,分析气动激励机理,建立用于电弧放电等离子体诱导激波数值模拟的爆炸丝传热模型.主要结论有:电弧放电等离子体气动激励的主要机理是热等离子体的热阻塞效应,热电弧放电对于超声速来流而言就像-个具有-定斜坡角度的虚拟突起;理论分析只适用于纵坐标较小的阶段;当传热的功率设为放电功率的10%时,本文所建立的模型能够用于电弧放电等离子体诱导激波的仿真研究;等离子体虚拟斜坡角度及其诱导激波角都随来流总压和速度的增大而减小,随着放电功率的增大而增大,在总压、速度和放电功率较小的阶段这种变化较明显,在总压、速度和放电功率较大的阶段这种变化较缓慢.     

Evolution of the ring-like vortices and spike structure in transitional boundary layers

At the late stage of transitional boundary layers, the nonlinear evolution of the ring-like vortices and spike structures and their effects on the surrounding flow were studied by means of direct numerical simulation with high order accuracy. A spatial transition of the flat-plate boundary layers in the compressible flow was conducted. Detailed numerical results with high resolution clearly represented the typical vortex structures, such as ring-like vortices and so on, and induced ejection and sweep events...     

量子阱层和垒层具有不同Al组分的270/290/330nm AlGaN基深紫外LED光电性能

为了研究AlGaN量子阱层和垒层中Al组分不同对AlGaN基深紫外发光二极管(LED)光电性能的影响,本文利用MOCVD生长、光刻和干法刻蚀工艺制备了AlGaN量子阱层和垒层具有不同Al组分的270/290/330nm深紫外LED,通过实验和数值模拟计算方法发现,量子阱层和垒层中具有低Al组分紫外LED的AlGaN材料具有较低的位错密度、较高的光输出功率和外量子效率。通过电流-电压(I-V)曲线拟合出的较大的理想因子(3.5)和能带结构图表明,AlGaN深紫外LED的电流产生是隧穿机制占据主导作用,这是因为高Al组分AlGaN量子阱中强极化场造成了有源层区域较大的能带弯曲和电势降。     

A high-order discontinuous Galerkin method for wave propagation through coupled elastic–acoustic media

We introduce a high-order discontinuous Galerkin (dG) scheme for the numerical solution of three-dimensional (3D) wave propagation problems in coupled elastic–acoustic media. A velocity–strain formulation is used, which allows for the solution of the acoustic and elastic wave equations within the same unified framework. Careful attention is directed at the derivation of a numerical flux that preserves high-order accuracy in the presence of material discontinuities, including elastic–acoustic interfaces. Explicit expressions for the 3D upwind numerical flux, derived as an exact solution for the relevant Riemann problem, are provided. The method supports h-non-conforming meshes, which are particularly effective at allowing local adaptation of the mesh size to resolve strong contrasts in the local wavelength, as well as dynamic adaptivity to track solution features. The use of high-order elements controls numerical dispersion, enabling propagation over many wave periods. We prove consistency and stability of the proposed dG scheme. To study the numerical accuracy and convergence of the proposed method, we compare against analytical solutions for wave propagation problems with interfaces, including Rayleigh, Lamb, Scholte, and Stoneley waves as well as plane waves impinging on an elastic–acoustic interface. Spectral rates of convergence are demonstrated for these problems, which include a non-conforming mesh case. Finally, we present scalability results for a parallel implementation of the proposed high-order dG scheme for large-scale seismic wave propagation in a simplified earth model, demonstrating high parallel efficiency for strong scaling to the full size of the Jaguar Cray XT5 supercomputer.     

Numerical solution for curved crack problem in elastic half-plane using hypersingular integral equation

A hypersingular integral equation for the curved crack problems of an elastic half-plane is introduced. Formulation of the equation is based on the usage of a modified complex potential. The potential is generally expressed in the form of a Cauchy-type integral. The modified complex potential is composed of the principal part and the complementary part. The principal part of the complex potential is actually equivalent to the original complex potential for the curved crack in an infinite plate. The role of the complementary part is to eliminate the boundary traction along the boundary of the half-plane caused by the principal part. From the assumed boundary traction condition, a hypersingular integral equation is obtained for the curved crack problems of an elastic half-plane. The curve length coordinate method is used to obtain a final solution. Several numerical examples are presented that prove the efficiency of the suggested method.     

Conditional dissipation of scalars in homogeneous turbulence: Closure for MMC modelling

While the mean and unconditional variance are to be predicted well by any reasonable turbulent combustion model, these are generally not sufficient for the accurate modelling of complex phenomena such as extinction/reignition. An additional criterion has been recently introduced: accurate modelling of the dissipation timescales associated with fluctuations of scalars about their conditional mean (conditional dissipation timescales). Analysis of Direct Numerical Simulation (DNS) results for a passive scalar shows that the conditional dissipation timescale is of the order of the integral timescale and smaller than the unconditional dissipation timescale. A model is proposed: the conditional dissipation timescale is proportional to the integral timescale. This model is used in Multiple Mapping Conditioning (MMC) modelling for a passive scalar case and a reactive scalar case, comparing to DNS results for both. The results show that this model improves the accuracy of MMC predictions so as to match the DNS results more closely using a relatively-coarse spatial resolution compared to other turbulent combustion models.