基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

New Falling Film Reactor for Melt Polycondensation Process

Summary: A grid falling film tower (GFFT) has been invented as an ideal polycondensation reactor. In this reactor, polymer melt flows through multi-layers grids from top to bottom to form falling film owing to gravity without agitation and shear; large gas-liquid interfacial area is generated; the grids are perpendicular between adjacent layers to ensure film renewal and to achieve uniformly flowing. The fluid flow in this reactor has little back-mixing and dead zone, which is near to plug flow. All melts are under the state of thin film which avoids the negative effect of hydrostatic head on the mass transfer impetus. Furthermore, the GFFT has wide operation flexibility as well as adjustable configuration parameters to meet different demands. A pilot scale GFFT with the height of 4.0 meters has been applied to polyester polycondensation process. The intrinsic viscosity of polyethylene terephthalate increased from 0.45 dl/g to 0.8–0.9 dl/g successfully. GFFT is supposed to be an universal apparatus for many devolatilization processes.     

基于混合网格Navier-Stokes方程的并行隐式计算方法研究

针对结构网格很难处理复杂外形和非结构网格无法计算具有边界层的粘性流动的缺点,发展了基于混合网格格点的隐式算法,成功地解决了在工程应用中难于处理的复杂外形粘性流场计算和效率问题。同时针对大规模的工程问题,发展了基于MPI通信技术的染色分层通讯并行计算方法。其中空间离散采用基于Roe格式发展的三阶迎风HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)或AUSM格式,湍流模型采用k??两方程湍流模型,时间推进考虑到LU-SGS并行等效较困难则采用基于DP-LUR(Data-Parallel Lower-Upper Relaxation)格式的隐式算法,计算CFL数可取到105量级,从2个到128个CPU的并行加速效率都保持在90%以上,大大提高了计算效率。算例对标模M6机翼模型流场进行计算,验证了方法的可靠性;然后对标模DLR-F6翼身组合体进行混合网格粘性与无粘计算结果进行比较,进一步验证混合网格方法;最后计算了DLR-WBNP外挂发动机翼身组合体模型,准确模拟了外挂和超临界机翼的相互干扰流动问题,采用4 CPU 16 CORE到24 CPU 96 CORE,2000步计算时间都不超...     

Two‐grid method for two‐dimensional nonlinear Schrödinger equation by finite element method

A conservative two‐grid finite element scheme is presented for the two‐dimensional nonlinear Schrödinger equation. One Newton iteration is applied on the fine grid to linearize the fully discrete problem using the coarse‐grid solution as the initial guess. Moreover, error estimates are conducted for the two‐grid method. It is shown that the coarse space can be extremely coarse, with no loss in the order of accuracy, and still achieve the asymptotically optimal approximation as long as the mesh sizes satisfy in the two‐grid method. The numerical results show that this method is very effective.     

多孔介质中两相不可压缩流体的时间局部网格加密有限差分格式

引入Charent压力变量,对于多孔介质中两相不可压缩流体的非混溶驱动问题,其模型表现为耦合的非线性偏微分方程组,一个是压力方程,另一个为饱和度方程．文中考虑一维问题且假定达西速度“已知,建立了在时间上进行局部加密的有限差分格式,给出了饱和度的最大模误差估计．最后给出了数值算例．     

一种基于并行混沌和SFLA的风电电网无功优化模型设计

针对已有方法在解决电网无功优化时,由于系统的无功不足和电网电压的不稳定,容易过早收敛到局部最优解的缺点,设计了一种基于并行混沌和混合蛙跳算法 (Shuffle Frog Leaping Algorithm, SFLA)的电网无功优化模型。首先,建立了最小化有功网损、最大化静态电压稳定裕度和最大化无功补偿单位投资收益的多目标数学优化模型,然后,对经典的SFLA进行改进,通过引入精英协同进化机制和划分种群的方式实现并行寻优,从而增加个体的多样性和加快最优解的求取速度,在不同种群中设计不同的适应度函数和个体更新进化方法。为了使得算法的初始解分布更为均匀,引入用混沌机制来对种群进行初始化,最后,对基于并行混沌和SFLA的总体算法进行了设计和分析。在Matlab环境下进行实验,实验结果表明文中方法得到的优化结果具有电网有功损耗小、单位投资收益高和静态电压稳定裕度大的优点,具有较强的可行性和适应性。     

平栅极结构碳纳米管场发射性能实验

研究了碳纳米管(CNT)场发射显示器(FED)三电极结构的平栅极结构,得到了进一步降低场致发射的开启电压和缩小动态调制电压范围的方法,同时也为相关的场发射安全操作提供了借鉴。实验表明:二极结构场发射调制电压范围较大,调制电压达上千伏,而在三电极的平栅极结构中通过调节阳极电压不仅可控制显示亮度,还对栅极调制电路有保护作用。适当升高阳极电压、适当缩短阴极和阳极之间的距离以及阴栅极经老化后可减小栅极调制电压, 同时还能有效的降低场致发射的动态调制电压的范围。这对新一代的显示器研制提供帮助。     

Schwarz domain decomposition for the incompressible Navier–Stokes equations in general co‐ordinates

This paper describes a domain decomposition method for the incompressible Navier–Stokes equations in general co‐ordinates. Domain decomposition techniques are needed for solving flow problems in complicated geometries while retaining structured grids on each of the subdomains. This is the so‐called block‐structured approach. It enables the use of fast vectorized iterative methods on the subdomains. The Navier–Stokes equations are discretized on a staggered grid using finite volumes. The pressure‐correction technique is used to solve the momentum equations together with incompressibility conditions. Schwarz domain decomposition is used to solve the momentum and pressure equations on the composite domain. Convergence of domain decomposition is accelerated by a GMRES Krylov subspace method. Computations are presented for a variety of flows. Copyright © 2000 John Wiley & Sons, Ltd.     

An unstructured grid,three‐dimensional model based on the shallow water equations

A semi‐implicit finite difference model based on the three‐dimensional shallow water equations is modified to use unstructured grids. There are obvious advantages in using unstructured grids in problems with a complicated geometry. In this development, the concept of unstructured orthogonal grids is introduced and applied to this model. The governing differential equations are discretized by means of a semi‐implicit algorithm that is robust, stable and very efficient. The resulting model is relatively simple, conserves mass, can fit complicated boundaries and yet is sufficiently flexible to permit local mesh refinements in areas of interest. Moreover, the simulation of the flooding and drying is included in a natural and straightforward manner. These features are illustrated by a test case for studies of convergence rates and by examples of flooding on a river plain and flow in a shallow estuary. Copyright © 2000 John Wiley & Sons, Ltd.